Search

Background: Infectious diseases of farmed and wild animals pose a recurrent threat to food security and human health. The macrophage, a key component of the innate immune system, is the first line of defence against many infectious agents and plays a major role in shaping the adaptive immune response. However, this phagocyte is a target and host for many pathogens. Understanding the molecular basis of interactions between macrophages and pathogens is therefore crucial for the development of effective strategies to combat important infectious diseases., Results: We explored how porcine pluripotent stem cells (PSCs) can provide a limitless in vitro supply of genetically and experimentally tractable macrophages. Porcine PSC-derived macrophages (PSCdMs) exhibited molecular and functional characteristics of ex vivo primary macrophages and were productively infected by pig pathogens, including porcine reproductive and respiratory syndrome virus (PRRSV) and African swine fever virus (ASFV), two of the most economically important and devastating viruses in pig farming. Moreover, porcine PSCdMs were readily amenable to genetic modification by CRISPR/Cas9 gene editing applied either in parental stem cells or directly in the macrophages by lentiviral vector transduction., Conclusions: We show that porcine PSCdMs exhibit key macrophage characteristics, including infection by a range of commercially relevant pig pathogens. In addition, genetic engineering of PSCs and PSCdMs affords new opportunities for functional analysis of macrophage biology in an important livestock species. PSCs and differentiated derivatives should therefore represent a useful and ethical experimental platform to investigate the genetic and molecular basis of host-pathogen interactions in pigs, and also have wider applications in livestock., (© 2021. The Author(s).)

Background: The domestic pig (Sus scrofa) is important both as a food source and as a biomedical model given its similarity in size, anatomy, physiology, metabolism, pathology, and pharmacology to humans. The draft reference genome (Sscrofa10.2) of a purebred Duroc female pig established using older clone-based sequencing methods was incomplete, and unresolved redundancies, short-range order and orientation errors, and associated misassembled genes limited its utility., Results: We present 2 annotated highly contiguous chromosome-level genome assemblies created with more recent long-read technologies and a whole-genome shotgun strategy, 1 for the same Duroc female (Sscrofa11.1) and 1 for an outbred, composite-breed male (USMARCv1.0). Both assemblies are of substantially higher (>90-fold) continuity and accuracy than Sscrofa10.2., Conclusions: These highly contiguous assemblies plus annotation of a further 11 short-read assemblies provide an unprecedented view of the genetic make-up of this important agricultural and biomedical model species. We propose that the improved Duroc assembly (Sscrofa11.1) become the reference genome for genomic research in pigs., (© The Author(s) 2020. Published by Oxford University Press.)

The domestic pig ( Sus scrofa ) is both an economically important livestock species and a model for biomedical research. Two highly contiguous pig reference genomes have recently been released. To support functional annotation of the pig genomes and comparative analysis with large human transcriptomic data sets, we aimed to create a pig gene expression atlas. To achieve this objective, we extended a previous approach developed for the chicken. We downloaded RNAseq data sets from public repositories, down-sampled to a common depth, and quantified expression against a reference transcriptome using the mRNA quantitation tool, Kallisto. We then used the network analysis tool Graphia to identify clusters of transcripts that were coexpressed across the merged data set. Consistent with the principle of guilt-by-association, we identified coexpression clusters that were highly tissue or cell-type restricted and contained transcription factors that have previously been implicated in lineage determination. Other clusters were enriched for transcripts associated with biological processes, such as the cell cycle and oxidative phosphorylation. The same approach was used to identify coexpression clusters within RNAseq data from multiple individual liver and brain samples, highlighting cell type, process, and region-specific gene expression. Evidence of conserved expression can add confidence to assignment of orthology between pig and human genes. Many transcripts currently identified as novel genes with ENSSSCG or LOC IDs were found to be coexpressed with annotated neighbouring transcripts in the same orientation, indicating they may be products of the same transcriptional unit. The meta-analytic approach to utilising public RNAseq data is extendable to include new data sets and new species and provides a framework to support the Functional Annotation of Animals Genomes (FAANG) initiative., (Copyright © 2020 Summers, Bush, Wu, Su, Muriuki, Clark, Finlayson, Eory, Waddell, Talbot, Archibald and Hume.)

Background: The Japanese quail (Coturnix japonica) is a popular domestic poultry species and an increasingly significant model species in avian developmental, behavioural and disease research., Results: We have produced a high-quality quail genome sequence, spanning 0.93 Gb assigned to 33 chromosomes. In terms of contiguity, assembly statistics, gene content and chromosomal organisation, the quail genome shows high similarity to the chicken genome. We demonstrate the utility of this genome through three diverse applications. First, we identify selection signatures and candidate genes associated with social behaviour in the quail genome, an important agricultural and domestication trait. Second, we investigate the effects and interaction of photoperiod and temperature on the transcriptome of the quail medial basal hypothalamus, revealing key mechanisms of photoperiodism. Finally, we investigate the response of quail to H5N1 influenza infection. In quail lung, many critical immune genes and pathways were downregulated after H5N1 infection, and this may be key to the susceptibility of quail to H5N1., Conclusions: We have produced a high-quality genome of the quail which will facilitate further studies into diverse research questions using the quail as a model avian species.

Background: Despite the significance of chicken as a model organism, our understanding of the chicken transcriptome is limited compared to human. This issue is common to all non-human vertebrate annotations due to the difficulty in transcript identification from short read RNAseq data. While previous studies have used single molecule long read sequencing for transcript discovery, they did not perform RNA normalization and 5'-cap selection which may have resulted in lower transcriptome coverage and truncated transcript sequences., Results: We sequenced normalised chicken brain and embryo RNA libraries with Pacific Bioscience Iso-Seq. 5' cap selection was performed on the embryo library to provide methodological comparison. From these Iso-Seq sequencing projects, we have identified 60 k transcripts and 29 k genes within the chicken transcriptome. Of these, more than 20 k are novel lncRNA transcripts with ~3 k classified as sense exonic overlapping lncRNA, which is a class that is underrepresented in many vertebrate annotations. The relative proportion of alternative transcription events revealed striking similarities between the chicken and human transcriptomes while also providing explanations for previously observed genomic differences., Conclusions: Our results indicate that the chicken transcriptome is similar in complexity compared to human, and provide insights into other vertebrate biology. Our methodology demonstrates the potential of Iso-Seq sequencing to rapidly expand our knowledge of transcriptomics.

Purpose: The aim of this study was to investigate the effect of VEGF-targeted therapy (sunitinib) on molecular intratumoral heterogeneity (ITH) in metastatic clear cell renal cancer (mccRCC)., Experimental Design: Multiple tumor samples (n = 187 samples) were taken from the primary renal tumors of patients with mccRCC who were sunitinib treated (n = 23, SuMR clinical trial) or untreated (n = 23, SCOTRRCC study). ITH of pathologic grade, DNA (aCGH), mRNA (Illumina Beadarray) and candidate proteins (reverse phase protein array) were evaluated using unsupervised and supervised analyses (driver mutations, hypoxia, and stromal-related genes). ITH was analyzed using intratumoral protein variance distributions and distribution of individual patient aCGH and gene-expression clustering., Results: Tumor grade heterogeneity was greater in treated compared with untreated tumors (P = 0.002). In unsupervised analysis, sunitinib therapy was not associated with increased ITH in DNA or mRNA. However, there was an increase in ITH for the driver mutation gene signature (DNA and mRNA) as well as increasing variability of protein expression with treatment (P < 0.05). Despite this variability, significant chromosomal and transcript changes to key targets of sunitinib, such as VHL, PBRM1, and CAIX, occurred in the treated samples., Conclusions: These findings suggest that sunitinib treatment has significant effects on the expression and ITH of key tumor and treatment specific genes/proteins in mccRCC. The results, based on primary tumor analysis, do not support the hypothesis that resistant clones are selected and predominate following targeted therapy., (©2015 American Association for Cancer Research.)

Background: Small insertions and deletions (InDels) constitute the second most abundant class of genetic variants and have been found to be associated with many traits and diseases. The present study reports on the detection and characterisation of about 883 K high quality InDels from the whole-genome analysis of several modern layer chicken lines from diverse breeds., Results: To reduce the error rates seen in InDel detection, this study used the consensus set from two InDel-calling packages: SAMtools and Dindel, as well as stringent post-filtering criteria. By analysing sequence data from 163 chickens from 11 commercial and 5 experimental layer lines, this study detected about 883 K high quality consensus InDels with 93% validation rate and an average density of 0.78 InDels/kb over the genome. Certain chromosomes, viz, GGAZ, 16, 22 and 25 showed very low densities of InDels whereas the highest rate was observed on GGA6. In spite of the higher recombination rates on microchromosomes, the InDel density on these chromosomes was generally lower relative to macrochromosomes possibly due to their higher gene density. About 43-87% of the InDels were found to be fixed within each line. The majority of detected InDels (86%) were 1-5 bases and about 63% were non-repetitive in nature while the rest were tandem repeats of various motif types. Functional annotation identified 613 frameshift, 465 non-frameshift and 10 stop-gain/loss InDels. Apart from the frameshift and stopgain/loss InDels that are expected to affect the translation of protein sequences and their biological activity, 33% of the non-frameshift were predicted as evolutionary intolerant with potential impact on protein functions. Moreover, about 2.5% of the InDels coincided with the most-conserved elements previously mapped on the chicken genome and are likely to define functional elements. InDels potentially affecting protein function were found to be enriched for certain gene-classes e.g. those associated with cell proliferation, chromosome and Golgi organization, spermatogenesis, and muscle contraction., Conclusions: The large catalogue of InDels presented in this study along with their associated information such as functional annotation, estimated allele frequency, etc. are expected to serve as a rich resource for application in future research and breeding in the chicken.

Background: There is a lack of biomarkers to predict outcome with targeted therapy in metastatic clear cell renal cancer (mccRCC). This may be because dynamic molecular changes occur with therapy., Objective: To explore if dynamic, targeted-therapy-driven molecular changes correlate with mccRCC outcome., Design, Setting, and Participants: Multiple frozen samples from primary tumours were taken from sunitinib-naïve (n=22) and sunitinib-treated mccRCC patients (n=23) for protein analysis. A cohort (n=86) of paired, untreated and sunitinib/pazopanib-treated mccRCC samples was used for validation. Array comparative genomic hybridisation (CGH) analysis and RNA interference (RNAi) was used to support the findings., Intervention: Three cycles of sunitinib 50mg (4 wk on, 2 wk off)., Outcome Measurements and Statistical Analysis: Reverse phase protein arrays (training set) and immunofluorescence automated quantitative analysis (validation set) assessed protein expression., Results and Limitations: Differential expression between sunitinib-naïve and treated samples was seen in 30 of 55 proteins (p<0.05 for each). The proteins B-cell CLL/lymphoma 2 (BCL2), mutL homolog 1 (MLH1), carbonic anhydrase 9 (CA9), and mechanistic target of rapamycin (mTOR) (serine/threonine kinase) had both increased intratumoural variance and significant differential expression with therapy. The validation cohort confirmed increased CA9 expression with therapy. Multivariate analysis showed high CA9 expression after treatment was associated with longer survival (hazard ratio: 0.48; 95% confidence interval, 0.26-0.87; p=0.02). Array CGH profiles revealed sunitinib was associated with significant CA9 region loss. RNAi CA9 silencing in two cell lines inhibited the antiproliferative effects of sunitinib. Shortcomings of the study include selection of a specific protein for analysis, and the specific time points at which the treated tissue was analysed., Conclusions: CA9 levels increase with targeted therapy in mccRCC. Lower CA9 levels are associated with a poor prognosis and possible resistance, as indicated by the validation cohort., Patient Summary: Drug treatment of advanced kidney cancer alters molecular markers of treatment resistance. Measuring carbonic anhydrase 9 levels may be helpful in determining which patients benefit from therapy., (Copyright © 2014 European Association of Urology. Published by Elsevier B.V. All rights reserved.)

Currently there is no curative treatment for metastatic clear cell renal cell cancer, the commonest variant of the disease. A key factor in this treatment resistance is thought to be the molecular complexity of the disease. Targeted therapy such as the tyrosine kinase inhibitor (TKI)-sunitinib have been utilized, but only 40% of patients will respond, with the overwhelming majority of these patients relapsing within 1 year. As such the question of intrinsic and acquired resistance in renal cell cancer patients is highly relevant. In order to study resistance to TKIs, with the ultimate goal of developing effective, personalized treatments, sequential tissue after a specific period of targeted therapy is required, an approach which had proved successful in chronic myeloid leukaemia. However the application of such a strategy in renal cell carcinoma is complicated by the high level of both inter- and intratumoral heterogeneity, which is a feature of renal cell carcinoma as well as other solid tumors. Intertumoral heterogeneity due to transcriptomic and genetic differences is well established even in patients with similar presentation, stage and grade of tumor. In addition it is clear that there is great morphological (intratumoral) heterogeneity in RCC, which is likely to represent even greater molecular heterogeneity. Detailed mapping and categorization of RCC tumors by combined morphological analysis and Fuhrman grading allows the selection of representative areas for proteomic analysis. Protein based analysis of RCC is attractive due to its widespread availability in pathology laboratories; however, its application can be problematic due to the limited availability of specific antibodies. Due to the dot blot nature of the Reverse Phase Protein Arrays (RPPA), antibody specificity must be pre-validated; as such strict quality control of antibodies used is of paramount importance. Despite this limitation the dot blot format does allow assay miniaturization, allowing for the printing of hundreds of samples onto a single nitrocellulose slide. Printed slides can then be analyzed in a similar fashion to Western analysis with the use of target specific primary antibodies and fluorescently labelled secondary antibodies, allowing for multiplexing. Differential protein expression across all the samples on a slide can then be analyzed simultaneously by comparing the relative level of fluorescence in a more cost-effective and high-throughput manner.

Protein-coding sequences make up only about 1% of the mammalian genome. Much of the remaining 99% has been long assumed to be junk DNA, with little or no functional significance. Here, we show that in hominids, a group with historically low effective population sizes, all classes of noncoding DNA evolve more slowly than ancestral transposable elements and so appear to be subject to significant evolutionary constraints. Under the nearly neutral theory, we expected to see lower levels of selective constraints on most sequence types in hominids than murids, a group that is thought to have a higher effective population size. We found that this is the case for many sequence types examined, the most extreme example being 5'UTRs, for which constraint in hominids is only about one-third that of murids. Surprisingly, however, we observed higher constraints for some sequence types in hominids, notably 4-fold sites, where constraint is more than twice as high as in murids. This implies that more than about one-fifth of mutations at 4-fold sites are effectively selected against in hominids. The higher constraint at 4-fold sites in hominids suggests a more complex protein-coding gene structure than murids and indicates that methods for detecting selection on protein-coding sequences (e.g., using the d(N)/d(S) ratio), with 4-fold sites as a neutral standard, may lead to biased estimates, particularly in hominids. Our constraint estimates imply that 5.4% of nucleotide sites in the human genome are subject to effective negative selection and that there are three times as many constrained sites within noncoding sequences as within protein-coding sequences. Including coding and noncoding sites, we estimate that the genomic deleterious mutation rate U = 4.2. The mutational load predicted under a multiplicative model is therefore about 99% in hominids.