A high resolution atlas of gene expression in the domestic sheep (Ovis aries)

Genetic signatures of high-altitude adaptation and geographic distribution in Tibetan sheep

Most sheep breeding programs designed for the tropics and sub-tropics have to take into account the impacts of environmental adaptive traits. However, the genetic mechanism regulating the multiple biological processes driving adaptive responses remains unclear. In this study, we applied a selective sweep analysis by combing 1% top values of Fst and ZHp on both altitude and geographic subpopulations (APS) in 636 indigenous Tibetan sheep breeds. Results show that 37 genes were identified within overlapped genomic regions regarding Fst significantly associated with APS. Out of the 37 genes, we found that 8, 3 and 6 genes at chromosomes (chr.) 13, 23 and 27, respectively, were identified in the genomic regions with 1% top values of ZHp. We further analyzed the INDEL variation of 6 genes at chr.27 (X chromosome) in APS together with corresponding orthologs of 6 genes in Capra, Pantholops, and Bos Taurus. We found that an INDEL was located within 5'UTR region of HAG1 gene. This INDEL of HAG1 was strongly associated with the variation of APS, which was further confirmed by qPCR. Sheep breeds carrying "C-INDEL" of HAG1 have significantly greater body weight, shear amount, corpuscular hemoglobin and globulin levels, but lower body height, than those carrying "CA-INDEL" of HAG1. We concluded that "C-INDEL" variation of HAG1 gene confers better hypoxia tolerance in the highlands of Tibetan and explains well geographic distributions in this population. These results contribute to our understanding of adaptive responses to altitude and geographic adaptation in Tibetan sheep populations and will help to guide future conservation programs for Tibetan sheep native to Qinghai-Tibetan Plateau.

Global Analysis of Transcription Start Sites in the New Ovine Reference Genome (Oar rambouillet v1.0)

The overall aim of the Ovine FAANG project is to provide a comprehensive annotation of the new highly contiguous sheep reference genome sequence (Oar rambouillet v1.0). Mapping of transcription start sites (TSS) is a key first step in understanding transcript regulation and diversity. Using 56 tissue samples collected from the reference ewe Benz2616, we have performed a global analysis of TSS and TSS-Enhancer clusters using Cap Analysis Gene Expression (CAGE) sequencing. CAGE measures RNA expression by 5' cap-trapping and has been specifically designed to allow the characterization of TSS within promoters to single-nucleotide resolution. We have adapted an analysis pipeline that uses TagDust2 for clean-up and trimming, Bowtie2 for mapping, CAGEfightR for clustering, and the Integrative Genomics Viewer (IGV) for visualization. Mapping of CAGE tags indicated that the expression levels of CAGE tag clusters varied across tissues. Expression profiles across tissues were validated using corresponding polyA+ mRNA-Seq data from the same samples. After removal of CAGE tags with <10 read counts, 39.3% of TSS overlapped with 5' ends of 31,113 transcripts that had been previously annotated by NCBI (out of a total of 56,308 from the NCBI annotation). For 25,195 of the transcripts, previously annotated by NCBI, no TSS meeting stringent criteria were identified. A further 14.7% of TSS mapped to within 50 bp of annotated promoter regions. Intersecting these predicted TSS regions with annotated promoter regions (±50 bp) revealed 46% of the predicted TSS were "novel" and previously un-annotated. Using whole-genome bisulfite sequencing data from the same tissues, we were able to determine that a proportion of these "novel" TSS were hypo-methylated (32.2%) indicating that they are likely to be reproducible rather than "noise". This global analysis of TSS in sheep will significantly enhance the annotation of gene models in the new ovine reference assembly. Our analyses provide one of the highest resolution annotations of transcript regulation and diversity in a livestock species to date.

Expression of mitochondrial protein genes encoded by nuclear and mitochondrial genomes correlate with energy metabolism in dairy cattle

Background

Mutations in the mitochondrial genome have been implicated in mitochondrial disease, often characterized by impaired cellular energy metabolism. Cellular energy metabolism in mitochondria involves mitochondrial proteins (MP) from both the nuclear (NuMP) and mitochondrial (MtMP) genomes. The expression of MP genes in tissues may be tissue specific to meet varying specific energy demands across the tissues. Currently, the characteristics of MP gene expression in tissues of dairy cattle are not well understood. In this study, we profile the expression of MP genes in 29 adult and six foetal tissues in dairy cattle using RNA sequencing and gene expression analyses: particularly differential gene expression and co-expression network analyses.

Results

MP genes were differentially expressed (DE; over-expressed or under-expressed) across tissues in cattle. All 29 tissues showed DE NuMP genes in varying proportions of over-expression and under-expression. On the other hand, DE of MtMP genes was observed in < 50% of tissues and notably MtMP genes within a tissue was either all over-expressed or all under-expressed. A high proportion of NuMP (up to 60%) and MtMP (up to 100%) genes were over-expressed in tissues with expected high metabolic demand; heart, skeletal muscles and tongue, and under-expressed (up to 45% of NuMP, 77% of MtMP genes) in tissues with expected low metabolic rates; leukocytes, thymus, and lymph nodes. These tissues also invariably had the expression of all MtMP genes in the direction of dominant NuMP genes expression. The NuMP and MtMP genes were highly co-expressed across tissues and co-expression of genes in a cluster were non-random and functionally enriched for energy generation pathway. The differential gene expression and co-expression patterns were validated in independent cow and sheep datasets.

Conclusions

The results of this study support the concept that there are biological interaction of MP genes from the mitochondrial and nuclear genomes given their over-expression in tissues with high energy demand and co-expression in tissues. This highlights the importance of considering MP genes from both genomes in future studies related to mitochondrial functions and traits related to energy metabolism.

Network analysis of transcriptomic diversity amongst resident tissue macrophages and dendritic cells in the mouse mononuclear phagocyte system

The mononuclear phagocyte system (MPS) is a family of cells including progenitors, circulating blood monocytes, resident tissue macrophages, and dendritic cells (DCs) present in every tissue in the body. To test the relationships between markers and transcriptomic diversity in the MPS, we collected from National Center for Biotechnology Information Gene Expression Omnibus (NCBI-GEO) a total of 466 quality RNA sequencing (RNA-seq) data sets generated from mouse MPS cells isolated from bone marrow, blood, and multiple tissues. The primary data were randomly downsized to a depth of 10 million reads and requantified. The resulting data set was clustered using the network analysis tool BioLayout. A sample-to-sample matrix revealed that MPS populations could be separated based upon tissue of origin. Cells identified as classical DC subsets, cDC1s and cDC2s, and lacking Fcgr1 (encoding the protein CD64) were contained within the MPS cluster, no more distinct than other MPS cells. A gene-to-gene correlation matrix identified large generic coexpression clusters associated with MPS maturation and innate immune function. Smaller coexpression gene clusters, including the transcription factors that drive them, showed higher expression within defined isolated cells, including monocytes, macrophages, and DCs isolated from specific tissues. They include a cluster containing Lyve1 that implies a function in endothelial cell (EC) homeostasis, a cluster of transcripts enriched in intestinal macrophages, and a generic lymphoid tissue cDC cluster associated with Ccr7. However, transcripts encoding Adgre1, Itgax, Itgam, Clec9a, Cd163, Mertk, Mrc1, Retnla, and H2-a/e (encoding class II major histocompatibility complex [MHC] proteins) and many other proposed macrophage subset and DC lineage markers each had idiosyncratic expression profiles. Coexpression of immediate early genes (for example, Egr1, Fos, Dusp1) and inflammatory cytokines and chemokines (tumour necrosis factor [Tnf], Il1b, Ccl3/4) indicated that all tissue disaggregation and separation protocols activate MPS cells. Tissue-specific expression clusters indicated that all cell isolation procedures also co-purify other unrelated cell types that may interact with MPS cells in vivo. Comparative analysis of RNA-seq and single-cell RNA-seq (scRNA-seq) data from the same lung cell populations indicated that MPS heterogeneity implied by global cluster analysis may be even greater at a single-cell level. This analysis highlights the power of large data sets to identify the diversity of MPS cellular phenotypes and the limited predictive value of surface markers to define lineages, functions, or subpopulations.

The equine mononuclear phagocyte system: The relevance of the horse as a model for understanding human innate immunity

The mononuclear phagocyte system (MPS) is a family of cells of related function that includes bone marrow progenitors, blood monocytes and resident tissue macrophages. Macrophages are effector cells in both innate and acquired immunity. They are a major resident cell population in every organ and their numbers increase in response to proinflammatory stimuli. Their function is highly regulated by a wide range of agonists, including lymphokines, cytokines and products of microorganisms. Macrophage biology has been studied most extensively in mice, yet direct comparisons of rodent and human macrophages have revealed many functional differences. In this review, we provide an overview of the equine MPS, describing the variation in the function and phenotype of macrophages depending on their location and the similarities and differences between the rodent, human and equine immune response. We discuss the use of the horse as a large animal model in which to study macrophage biology and pathological processes shared with humans. Finally, following the recent update to the horse genome, facilitating further comparative analysis of regulated gene expression between the species, we highlight the importance of future transcriptomic macrophage studies in the horse, the findings of which may also be applicable to human as well as veterinary research.

Functional contexts of adipose and gluteal muscle tissue gene co-expression networks in the domestic horse

A gene's response to an environment is tightly bound to the underlying genetic variation present in an individual's genome and varies greatly depending on the tissue it is being expressed in. Gene co-expression networks provide a mechanism to understand and interpret the collective transcriptional responses of genes. Here, we use the Camoco co-expression network framework to characterize the transcriptional landscape of adipose and gluteal muscle tissue in 83 domestic horses (Equus caballus) representing 5 different breeds. In each tissue, gene expression profiles, capturing transcriptional response due to variation across individuals, were used to build two separate, tissue-focused, genotypically-diverse gene co-expression networks. The aim of our study was to identify significantly co-expressed clusters of genes in each tissue, then compare the clusters across networks to quantify the extent that clusters were found in both networks as well as to identify clusters found in a single network. The known and unknown functions for each network were quantified using complementary, supervised and unsupervised approaches. First, supervised ontological enrichment was utilized to quantify biological functions represented by each network. Curated ontologies (GO and KEGG) were used to measure the known biological functions present in each tissue. Overall, a large percentage of terms (40.3% of GO and 41% of KEGG) were co-expressed in at least one tissue. Many terms were co-expressed in both tissues, however a small proportion of terms exhibited single tissue co-expression suggesting functional differentiation based on curated, functional annotation. To complement this, an unsupervised approach not relying on ontologies was employed. Strongly co-expressed sets of genes defined by Markov clustering identified sets of unannotated genes showing similar patterns of co-expression within a tissue. We compared gene sets across tissues and identified clusters of genes the either segregate in co-expression by tissue or exhibit high levels of co-expression in both tissues. Clusters were also integrated with GO and KEGG ontologies to identify gene sets containing previously curated annotations versus unannotated gene sets indicating potentially novel biological function. Coupling together these transcriptional datasets, we mapped the transcriptional landscape of muscle and adipose setting up a generalizable framework for interpreting gene function for additional tissues in the horse and other species.

Expression of Calcification and Extracellular Matrix Genes in the Cardiovascular System of the Healthy Domestic Sheep (Ovis aries)

The maintenance of a healthy cardiovascular system requires expression of genes that contribute to essential biological activities and repression of those that are associated with functions likely to be detrimental to cardiovascular homeostasis. Vascular calcification is a major disruption to cardiovascular homeostasis, where tissues of the cardiovascular system undergo ectopic calcification and consequent dysfunction, but little is known about the expression of calcification genes in the healthy cardiovascular system. Large animal models are of increasing importance in cardiovascular disease research as they demonstrate more similar cardiovascular features (in terms of anatomy, physiology and size) to humans than do rodent species. We used RNA sequencing results from the sheep, which has been utilized extensively to examine calcification of prosthetic cardiac valves, to explore the transcriptome of the heart and cardiac valves in this large animal, in particular looking at expression of calcification and extracellular matrix genes. We then examined genes implicated in the process of vascular calcification in a wide array of cardiovascular tissues and across multiple developmental stages, using RT-qPCR. Our results demonstrate that there is a balance between genes that promote and those that suppress mineralization during development and across cardiovascular tissues. We show extensive expression of genes encoding proteins involved in formation and maintenance of the extracellular matrix in cardiovascular tissues, and high expression of hematopoietic genes in the cardiac valves. Our analysis will support future research into the functions of implicated genes in the development of valve calcification, and increase the utility of the sheep as a large animal model for understanding ectopic calcification in cardiovascular disease. This study provides a foundation to explore the transcriptome of the developing cardiovascular system and is a valuable resource for the fields of mammalian genomics and cardiovascular research.

Species-Specificity of Transcriptional Regulation and the Response to Lipopolysaccharide in Mammalian Macrophages

Mammalian macrophages differ in their basal gene expression profiles and response to the toll-like receptor 4 (TLR4) agonist, lipopolysaccharide (LPS). In human macrophages, LPS elicits a temporal cascade of transient gene expression including feed forward activators and feedback regulators that limit the response. Here we present a transcriptional network analysis of the response of sheep bone marrow-derived macrophages (BMDM) to LPS based upon RNA-seq at 0, 2, 4, 7, and 24 h post-stimulation. The analysis reveals a conserved transcription factor network with humans, and rapid induction of feedback regulators that constrain the response at every level. The gene expression profiles of sheep BMDM at 0 and 7 h post LPS addition were compared to similar data obtained from goat, cow, water buffalo, horse, pig, mouse and rat BMDM. This comparison was based upon identification of 8,200 genes annotated in all species and detected at >10TPM in at least one sample. Analysis of expression of transcription factors revealed a conserved transcriptional millieu associated with macrophage differentiation and LPS response. The largest co-expression clusters, including genes encoding cell surface receptors, endosome–lysosome components and secretory activity, were also expressed in all species and the combined dataset defines a macrophage functional transcriptome. All of the large animals differed from rodents in lacking inducible expression of genes involved in arginine metabolism and nitric oxide production. Instead, they expressed inducible transporters and enzymes of tryptophan and kynurenine metabolism. BMDM from all species expressed high levels of transcripts encoding transporters and enzymes involved in glutamine metabolism suggesting that glutamine is a major metabolic fuel. We identify and discuss transcripts that were uniquely expressed or regulated in rodents compared to large animals including ACOD1, CXC and CC chemokines, CD163, CLEC4E, CPM, CSF1, CSF2, CTSK, MARCO, MMP9, SLC2A3, SLC7A7, and SUCNR1. Conversely, the data confirm the conserved regulation of multiple transcripts for which there is limited functional data from mouse models and knockouts. The data provide a resource for functional annotation and interpretation of loci involved in susceptibility to infectious and inflammatory disease in humans and large animal species.

Microevolutionary Dynamics of Chicken Genomes under Divergent Selection for Adiposity

Decades of artificial selection have significantly improved performance and efficiency of animal production systems. However, little is known about the microevolution of genomes due to intensive breeding. Using whole-genome sequencing, we document dynamic changes of chicken genomes under divergent selection on adiposity over 19 generations. Directional selection reduced within-line but increased between-line genomic differences. We observed that artificial selection tended to result in recruitment of preexisting variations of genes related to adipose tissue growth. In addition, novel mutations contributed to divergence of phenotypes under selection but contributed significantly less than preexisting genomic variants. Integration of 15 generations genome sequencing, genome-wide association study, and multi-omics data further identified that genes involved in signaling pathways important to adipogenesis, such as autophagy and lysosome (URI1, MBL2), neural system (CHAT), and endocrine (PCSK1) pathways, were under strong selection. Our study provides insights into the microevolutionary dynamics of domestic animal genomes under artificial selection.

The origin of domestication genes in goats

Goat domestication was critical for agriculture and civilization, but its underlying genetic changes and selection regimes remain unclear. Here, we analyze the genomes of worldwide domestic goats, wild caprid species, and historical remains, providing evidence of an ancient introgression event from a West Caucasian tur-like species to the ancestor of domestic goats. One introgressed locus with a strong signature of selection harbors the MUC6 gene, which encodes a gastrointestinally secreted mucin. Experiments revealed that the nearly fixed introgressed haplotype confers enhanced immune resistance to gastrointestinal pathogens. Another locus with a strong signal of selection may be related to behavior. The selected alleles at these two loci emerged in domestic goats at least 7200 and 8100 years ago, respectively, and increased to high frequencies concurrent with the expansion of the ubiquitous modern mitochondrial haplogroup A. Tracking these archaeologically cryptic evolutionary transformations provides new insights into the mechanisms of animal domestication.

Novel Variants in GDF9 Gene Affect Promoter Activity and Litter Size in Mongolia Sheep

Litter size is an economically important trait in sheep breeding. The objectives of this study were as follows: (1) to ascertain if any of the 19 known variants in the BMPRIB, BMP15, and GDF9 genes are present and associated with the litter size of Mongolia sheep; (2) to identify novel variants in GDF9 and perform association analysis; and (3) to validate the effects of these GDF9 promoter variants on the activity of the gene. The results of the 19 known variants showed that the FecB^B affected the litter size of Mongolia sheep (p < 0.001). The association analysis results of novel variants showed that the g.46544883A>G (GenBank accession: NC_040256, the same below) in the 3' untranslated region (3' UTR), the c.1040T>C (Phe347Ser) in the exon 2, and the g.46547859C>T SNP in the promotor of GDF9 were significantly associated with litter size of Mongolia ewes (p < 0.01, p < 0.05, and p < 0.001, respectively). In addition, the GDF9 promoter activity analysis showed that the C allele at the -332 position (g.46547859C>T) could decrease luciferase activity compared with the T allele (p < 0.01). Our findings may facilitate effective marker-assisted selection to increase litter size in Mongolia sheep populations, as well as bring new insights into GDF9 expression.

Novel insights into cardiac regeneration based on differential fetal and adult ovine heart transcriptomic analysis

The adult mammalian cardiomyocyte has a very limited capacity to reenter the cell cycle and advance into mitosis. Therefore, diseases characterized by lost contractile tissue usually evolve into myocardial remodeling and heart failure. Analyzing the cardiac transcriptome at different developmental stages in a large mammal closer to the human than laboratory rodents may serve to disclose positive and negative cardiomyocyte cell cycle regulators potentially targetable to induce cardiac regeneration in the clinical setting. Thus we aimed at characterizing the transcriptomic profiles of the early fetal, late fetal, and adult sheep heart by employing RNA-seq technique and bioinformatic analysis to detect protein-encoding genes that in some of the stages were turned off, turned on, or differentially expressed. Genes earlier proposed as positive cell cycle regulators such as cyclin A, cdk2, meis2, meis3, and PCNA showed higher expression in fetal hearts and lower in AH, as expected. In contrast, genes previously proposed as cell cycle inhibitors, such as meis1, p16, and sav1, tended to be higher in fetal than in adult hearts, suggesting that these genes are involved in cell processes other than cell cycle regulation. Additionally, we described Gene Ontology (GO) enrichment of different sets of genes. GO analysis revealed that differentially expressed gene sets were mainly associated with metabolic and cellular processes. The cell cycle-related genes fam64a, cdc20, and cdk1, and the metabolism-related genes pitx and adipoq showed strong differential expression between fetal and adult hearts, thus being potent candidates to be targeted in human cardiac regeneration strategies.NEW & NOTEWORTHY We characterized the transcriptomic profiles of the fetal and adult sheep hearts employing RNAseq technique and bioinformatic analyses to provide sets of transcripts whose variation in expression level may link them to a specific role in cell cycle regulation. It is important to remark that this study was performed in a large mammal closer to humans than laboratory rodents. In consequence, the results can be used for further translational studies in cardiac regeneration.

Functional Annotation of the Transcriptome of the Pig, Sus scrofa, Based Upon Network Analysis of an RNAseq Transcriptional Atlas

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.

Enrichment of ovine gonadotropes via adenovirus gene targeting enhances assessment of transcriptional changes in response to estradiol-17 beta†

Gonadotropes represent approximately 5-15% of the total endocrine cell population in the mammalian anterior pituitary. Therefore, assessing the effects of experimental manipulation on virtually any parameter of gonadotrope biology is difficult to detect and parse from background noise. In non-rodent species, applying techniques such as high-throughput ribonucleic acid (RNA) sequencing is problematic due to difficulty in isolating and analyzing individual endocrine cell populations. Herein, we exploited cell-specific properties inherent to the proximal promoter of the human glycoprotein hormone alpha subunit gene (CGA) to genetically target the expression of a fluorescent reporter (green fluorescent protein [GFP]) selectively to ovine gonadotropes. Dissociated ovine pituitary cells were cultured and infected with an adenoviral reporter vector (Ad-hαCGA-eGFP). We established efficient gene targeting by successfully enriching dispersed GFP-positive cells with flow cytometry. Confirming enrichment of gonadotropes specifically, we detected elevated levels of luteinizing hormone (LH) but not thyrotropin-stimulating hormone (TSH) in GFP-positive cell populations compared to GFP-negative populations. Subsequently, we used next-generation sequencing to obtain the transcriptional profile of GFP-positive ovine gonadotropes in the presence or absence of estradiol 17-beta (E2), a key modulator of gonadotrope function. Compared to non-sorted cells, enriched GFP-positive cells revealed a distinct transcriptional profile consistent with established patterns of gonadotrope gene expression. Importantly, we also detected nearly 200 E2-responsive genes in enriched gonadotropes, which were not apparent in parallel experiments on non-enriched cell populations. From these data, we conclude that CGA-targeted adenoviral gene transfer is an effective means for selectively labeling and enriching ovine gonadotropes suitable for investigation by numerous experimental approaches.

Genetic and genomic analyses underpin the feasibility of concomitant genetic improvement of milk yield and mastitis resistance in dairy sheep

Milk yield is the most important dairy sheep trait and constitutes the key genetic improvement goal via selective breeding. Mastitis is one of the most prevalent diseases, significantly impacting on animal welfare, milk yield and quality, while incurring substantial costs. Our objectives were to determine the feasibility of a concomitant genetic improvement programme for enhanced milk production and resistance to mastitis. Individual records for milk yield, and four mastitis-related traits (milk somatic cell count, California Mastitis Test score, total viable bacterial count in milk and clinical mastitis presence) were collected monthly throughout lactation for 609 ewes of the Chios breed. All ewes were genotyped with a mastitis specific custom-made 960 single nucleotide polymorphism (SNP) array. We performed targeted genomic association studies, (co)variance component estimation and pathway enrichment analysis, and characterised gene expression levels and the extent of allelic expression imbalance. Presence of heritable variation for milk yield was confirmed. There was no significant genetic correlation between milk yield and mastitis traits. Environmental factors appeared to favour both milk production and udder health. There were no overlapping of SNPs associated with mastitis resistance and milk yield in Chios sheep. Furthermore, four distinct Quantitative Trait Loci (QTLs) affecting milk yield were detected on chromosomes 2, 12, 16 and 19, in locations other than those previously identified to affect mastitis resistance. Five genes (DNAJA1, GHR, LYPLA1, NUP35 and OXCT1) located within the QTL regions were highly expressed in both the mammary gland and milk transcriptome, suggesting involvement in milk synthesis and production. Furthermore, the expression of two of these genes (NUP35 and OXCT1) was enriched in immune tissues implying a potentially pleiotropic effect or likely role in milk production during udder infection, which needs to be further elucidated in future studies. In conclusion, the absence of genetic antagonism between milk yield and mastitis resistance suggests that simultaneous genetic improvement of both traits be achievable.

A Mini-Atlas of Gene Expression for the Domestic Goat (Capra hircus)

Goats (Capra hircus) are an economically important livestock species providing meat and milk across the globe. They are of particular importance in tropical agri-systems contributing to sustainable agriculture, alleviation of poverty, social cohesion, and utilisation of marginal grazing. There are excellent genetic and genomic resources available for goats, including a highly contiguous reference genome (ARS1). However, gene expression information is limited in comparison to other ruminants. To support functional annotation of the genome and comparative transcriptomics, we created a mini-atlas of gene expression for the domestic goat. RNA-Seq analysis of 17 transcriptionally rich tissues and 3 cell-types detected the majority (90%) of predicted protein-coding transcripts and assigned informative gene names to more than 1000 previously unannotated protein-coding genes in the current reference genome for goat (ARS1). Using network-based cluster analysis, we grouped genes according to their expression patterns and assigned those groups of coexpressed genes to specific cell populations or pathways. We describe clusters of genes expressed in the gastro-intestinal tract and provide the expression profiles across tissues of a subset of genes associated with functional traits. Comparative analysis of the goat atlas with the larger sheep gene expression atlas dataset revealed transcriptional similarities between macrophage associated signatures in the sheep and goats sampled in this study. The goat transcriptomic resource complements the large gene expression dataset we have generated for sheep and contributes to the available genomic resources for interpretation of the relationship between genotype and phenotype in small ruminants.

Transcriptional Response of Ovine Lung to Infection with Jaagsiekte Sheep Retrovirus

Jaagsiekte sheep retrovirus (JSRV) is the etiologic agent of ovine pulmonary adenocarcinoma (OPA), a neoplastic lung disease of sheep. OPA is an important economic and welfare issue for sheep farmers and a valuable naturally occurring animal model for human lung adenocarcinoma. Here, we used RNA sequencing to study the transcriptional response of ovine lung tissue to infection by JSRV. We identified 1,971 ovine genes differentially expressed in JSRV-infected lung compared to noninfected lung, including many genes with roles in carcinogenesis and immunomodulation. The differential expression of selected genes was confirmed using immunohistochemistry and reverse transcription-quantitative PCR. A key finding was the activation of anterior gradient 2, yes-associated protein 1, and amphiregulin in OPA tumor cells, indicating a role for this oncogenic pathway in OPA. In addition, there was differential expression of genes related to innate immunity, including genes encoding cytokines, chemokines, and complement system proteins. In contrast, there was little evidence for the upregulation of genes involved in T-cell immunity. Many genes related to macrophage function were also differentially expressed, reflecting the increased abundance of these cells in OPA-affected lung tissue. Comparison of the genes differentially regulated in OPA with the transcriptional changes occurring in human lung cancer revealed important similarities and differences between OPA and human lung adenocarcinoma. This study provides valuable new information on the pathogenesis of OPA and strengthens the use of this naturally occurring animal model for human lung adenocarcinoma.IMPORTANCE Ovine pulmonary adenocarcinoma is a chronic respiratory disease of sheep caused by jaagsiekte sheep retrovirus (JSRV). OPA is a significant economic problem for sheep farmers in many countries and is a valuable animal model for some forms of human lung cancer. Here, we examined the changes in host gene expression that occur in the lung in response to JSRV infection. We identified a large number of genes with altered expression in infected lung, including factors with roles in cancer and immune system function. We also compared the data from OPA to previously published data from human lung adenocarcinoma and found a large degree of overlap in the genes that were dysregulated. The results of this study provide exciting new avenues for future studies of OPA and may have comparative relevance for understanding human lung cancer.

Quantifying the contribution of sequence variants with regulatory and evolutionary significance to 34 bovine complex traits

Many genome variants shaping mammalian phenotype are hypothesized to regulate gene transcription and/or to be under selection. However, most of the evidence to support this hypothesis comes from human studies. Systematic evidence for regulatory and evolutionary signals contributing to complex traits in a different mammalian model is needed. Sequence variants associated with gene expression (expression quantitative trait loci [eQTLs]) and concentration of metabolites (metabolic quantitative trait loci [mQTLs]) and under histone-modification marks in several tissues were discovered from multiomics data of over 400 cattle. Variants under selection and evolutionary constraint were identified using genome databases of multiple species. These analyses defined 30 sets of variants, and for each set, we estimated the genetic variance the set explained across 34 complex traits in 11,923 bulls and 32,347 cows with 17,669,372 imputed variants. The per-variant trait heritability of these sets across traits was highly consistent (r > 0.94) between bulls and cows. Based on the per-variant heritability, conserved sites across 100 vertebrate species and mQTLs ranked the highest, followed by eQTLs, young variants, those under histone-modification marks, and selection signatures. From these results, we defined a Functional-And-Evolutionary Trait Heritability (FAETH) score indicating the functionality and predicted heritability of each variant. In additional 7,551 cattle, the high FAETH-ranking variants had significantly increased genetic variances and genomic prediction accuracies in 3 production traits compared to the low FAETH-ranking variants. The FAETH framework combines the information of gene regulation, evolution, and trait heritability to rank variants, and the publicly available FAETH data provide a set of biological priors for cattle genomic selection worldwide.

Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep

Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.

Functional annotation of the cattle genome through systematic discovery and characterization of chromatin states and butyrate-induced variations

BACKGROUND

The functional annotation of genomes, including chromatin accessibility and modifications, is important for understanding and effectively utilizing the increased amount of genome sequences reported. However, while such annotation has been well explored in a diverse set of tissues and cell types in human and model organisms, relatively little data are available for livestock genomes, hindering our understanding of complex trait variation, domestication, and adaptive evolution. Here, we present the first complete global landscape of regulatory elements in cattle and explore the dynamics of chromatin states in rumen epithelial cells induced by the rumen developmental regulator-butyrate.

RESULTS

We established the first global map of regulatory elements (15 chromatin states) and defined their coordinated activities in cattle, through genome-wide profiling for six histone modifications, RNA polymerase II, CTCF-binding sites, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial primary cells (REPC), rumen tissues, and Madin-Darby bovine kidney epithelial cells (MDBK). We demonstrated that each chromatin state exhibited specific enrichment for sequence ontology, transcription, methylation, trait-associated variants, gene expression-associated variants, selection signatures, and evolutionarily conserved elements, implying distinct biological functions. After butyrate treatments, we observed that the weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states, occurred concomitantly with significant alterations in gene expression and DNA methylation, which was significantly associated with heifer conception rate and stature economic traits.

CONCLUSION

Our results demonstrate the crucial role of functional genome annotation for understanding genome regulation, complex trait variation, and adaptive evolution in livestock. Using butyrate to induce the dynamics of the epigenomic landscape, we were able to establish the correlation among nutritional elements, chromatin states, gene activities, and phenotypic outcomes.

Selection signature analysis reveals genes underlying sheep milking performance

Sheep milk is the most important feed resource for newborn lambs and an important food resource for humans. Sheep milk production and ingredients are influenced by genetic and environmental factors. In this study, we implemented selection signature analysis using Illumina Ovine SNP50 BeadChip data of 78 meat Lacaune and 103 milk Lacaune sheep, which have similar genetic backgrounds, from the Sheep HapMap project to identify candidate genes related to ovine milk traits. Since different methods can detect different variation types and complement each other, we used a haplotype-based method (hapFLK) to implement selection signature analysis. The results revealed six selection signature regions showing signs of being selected ( P < 0.001 ): chromosomes 1, 2, 3, 6, 13 and 18. In addition, 38 quantitative trait loci (QTLs) related to sheep milk performance were identified in selection signature regions, which contain 334 candidate genes. Of those, SUCNR1 (succinate receptor 1) and PPARGC1A (PPARG coactivator 1 alpha) may be the most significant genes that affect sheep milking performance, which supply a significant indication for future studies to investigate candidate genes that play an important role in milk production and quality.

A Gene Expression Atlas of the Domestic Water Buffalo (Bubalus bubalis)

The domestic water buffalo (Bubalus bubalis) makes a major contribution to the global agricultural economy in the form of milk, meat, hides, and draught power. The global water buffalo population is predominantly found in Asia, and per head of population more people depend upon the buffalo than on any other livestock species. Despite its agricultural importance, there are comparatively fewer genomic and transcriptomic resources available for buffalo than for other livestock species. We have generated a large-scale gene expression atlas covering multiple tissue and cell types from all major organ systems collected from three breeds of riverine water buffalo (Mediterranean, Pandharpuri and Bhadawari) and used the network analysis tool Graphia Professional to identify clusters of genes with similar expression profiles. Alongside similar data, we and others have generated for ruminants as part of the Functional Annotation of Animal Genomes Consortium; this comprehensive transcriptome supports functional annotation and comparative analysis of the water buffalo genome.

Phenotypic impacts of CSF1R deficiencies in humans and model organisms

Mϕ proliferation, differentiation, and survival are controlled by signals from the Mϕ CSF receptor (CSF1R). Mono-allelic gain-of-function mutations in CSF1R in humans are associated with an autosomal-dominant leukodystrophy and bi-allelic loss-of-function mutations with recessive skeletal dysplasia, brain disorders, and developmental anomalies. Most of the phenotypes observed in these human disease states are also observed in mice and rats with loss-of-function mutations in Csf1r or in Csf1 encoding one of its two ligands. Studies in rodent models also highlight the importance of genetic background and likely epistatic interactions between Csf1r and other loci. The impacts of Csf1r mutations on the brain are usually attributed solely to direct impacts on microglial number and function. However, analysis of hypomorphic Csf1r mutants in mice and several other lines of evidence suggest that primary hydrocephalus and loss of the physiological functions of Mϕs in the periphery contribute to the development of brain pathology. In this review, we outline the evidence that CSF1R is expressed exclusively in mononuclear phagocytes and explore the mechanisms linking CSF1R mutations to pleiotropic impacts on postnatal growth and development.

Genetic and epigenetic architecture of paternal origin contribute to gestation length in cattle

The length of gestation can affect offspring health and performance. Both maternal and fetal effects contribute to gestation length; however, paternal contributions to gestation length remain elusive. Using genome-wide association study (GWAS) in 27,214 Holstein bulls with millions of gestation records, here we identify nine paternal genomic loci associated with cattle gestation length. We demonstrate that these GWAS signals are enriched in pathways relevant to embryonic development, and in differentially methylated regions between sperm samples with long and short gestation length. We reveal that gestation length shares genetic and epigenetic architecture in sperm with calving ability, body depth, and conception rate. While several candidate genes are detected in our fine-mapping analysis, we provide evidence indicating ZNF613 as a promising candidate for cattle gestation length. Collectively, our findings support that the paternal genome and epigenome can impact gestation length potentially through regulation of the embryonic development.

Lingzhao Fang et al. studied the paternal genetic variants that affect gestational length in cattle. They found that paternal genes from pathways involved in embryonic development were associated with gestation length, and that these were often found in differentially methylated regions of the genome.

Evolutionary Loss of Genomic Proximity to Conserved Noncoding Elements Impacted the Gene Expression Dynamics During Mammalian Brain Development

Conserved noncoding elements (CNEs) have a significant regulatory influence on their neighboring genes. Loss of proximity to CNEs through genomic rearrangements can, therefore, impact the transcriptional states of the cognate genes. Yet, the evolutionary implications of such chromosomal alterations have not been studied. Through genome-wide analysis of CNEs and the cognate genes of representative species from five different mammalian orders, we observed a significant loss of genes' linear proximity to CNEs in the rat lineage. The CNEs and the genes losing proximity had a significant association with fetal, but not postnatal, brain development as assessed through ontology terms, developmental gene expression, chromatin marks, and genetic mutations. The loss of proximity to CNEs correlated with the independent evolutionary loss of fetus-specific upregulation of nearby genes in the rat brain. DNA breakpoints implicated in brain abnormalities of germline origin had significant representation between a CNE and the gene that exhibited loss of proximity, signifying the underlying developmental tolerance of genomic rearrangements that allowed the evolutionary splits of CNEs and the cognate genes in the rodent lineage. Our observations highlighted a nontrivial impact of chromosomal rearrangements in shaping the evolutionary dynamics of mammalian brain development and might explain the loss of brain traits, like cerebral folding of the cortex, in the rodent lineage.

Comprehensive Transcriptional Profiling of the Gastrointestinal Tract of Ruminants from Birth to Adulthood Reveals Strong Developmental Stage Specific Gene Expression

One of the most significant physiological challenges to neonatal and juvenile ruminants is the development and establishment of the rumen. Using a subset of RNA-Seq data from our high-resolution atlas of gene expression in sheep (Ovis aries) we have provided the first comprehensive characterization of transcription of the entire gastrointestinal (GI) tract during the transition from pre-ruminant to ruminant. The dataset comprises 164 tissue samples from sheep at four different time points (birth, one week, 8 weeks and adult). Using network cluster analysis we illustrate how the complexity of the GI tract is reflected in tissue- and developmental stage-specific differences in gene expression. The most significant transcriptional differences between neonatal and adult sheep were observed in the rumen complex. Comparative analysis of gene expression in three GI tract tissues from age-matched sheep and goats revealed species-specific differences in genes involved in immunity and metabolism. This study improves our understanding of the transcriptomic mechanisms involved in the transition from pre-ruminant to ruminant by identifying key genes involved in immunity, microbe recognition and metabolism. The results form a basis for future studies linking gene expression with microbial colonization of the developing GI tract and provide a foundation to improve ruminant efficiency and productivity through identifying potential targets for novel therapeutics and gene editing.

Dosage Compensation and Gene Expression of the X Chromosome in Sheep

Ohno's hypothesis predicts that the expression of the single X chromosome in males needs compensatory upregulation to balance its dosage with that of the diploid autosomes. Additionally, X chromosome inactivation ensures that quadruple expression of the two X chromosomes is avoided in females. These mechanisms have been actively studied in mice and humans but lag behind in domestic species. Using RNA sequencing data, we analyzed the X chromosome upregulation in sheep fetal tissues from day 135 of gestation under control, over or restricted maternal diets (100%, 140% and 60% of National Research Council Total Digestible Nutrients), and in conceptuses, juvenile, and adult somatic tissues. By computing the mean expression ratio of all X-linked genes to all autosomal genes (X:A), we found that all samples displayed some levels of X chromosome upregulation. The degrees of X upregulation were not significant (P-value = 0.74) between ovine females and males in the same somatic tissues. Brain, however, displayed complete X upregulation. Interestingly, the male and female reproduction-related tissues exhibited divergent X dosage upregulation. Moreover, expression upregulation of the X chromosome in fetal tissues was not affected by maternal diets. Maternal nutrition, however, did change expression levels of several X-linked genes, such as sex determination genes SOX3 and NR0B1 In summary, our results showed that X chromosome upregulation occurred in nearly all sheep somatic tissues analyzed, thus support Ohno's hypothesis in a new species. However, the levels of upregulation differed by different subgroups of genes such as those that are house-keeping and "dosage-sensitive".

The Structure, Evolution, and Gene Expression Within the Caprine Leukocyte Receptor Complex

The leukocyte receptor complex (LRC) encodes a large number of immunoglobulin (Ig)-like receptors involved in the immune response, particularly in modulating natural killer (NK) cell function. The killer cell Ig-like receptors (KIR), the leukocyte Ig-like receptors (LILR), and a recently described novel Ig-like receptor family are highly variable between species, which is consistent with rapid evolution driven by selection pressure from pathogens. Among the species studied to date, only simians (such as humans) and bovids (such as cattle and goats) have an expanded complement of KIR genes and represent an interesting model to study KIR evolution. Using recently improved genome assemblies and an assembly of bacterial artificial chromosomes, we describe the structure of the LRC, and the KIR region in particular, in goats and compare this to sheep as the assemblies allow. These species diverged from a common ancestor ~10 million years ago and from cattle ~25 million years ago. We identified conserved KIR genes common to both goats and sheep and confirm a partial sheep haplotype shared between the Rambouillet and Texel breeds. Goats and sheep have independently expanded two novel KIR subgroups, and unlike cattle or any other mammal, they do not appear to possess a functional 3DL-lineage KIR gene. Investigation of LRC gene expression using available transcriptomic data for various sheep and goat tissues largely confirmed putative gene annotation and revealed that a relatively conserved caprinae-specific KIR subgroup is expressed in macrophages. The LILR and novel Ig-like receptors were also highly expressed across a diverse range of tissues. This further step toward our understanding of the LRC receptor repertoire will help inform future studies investigating immune response variation in these species.

Integrating Selection Mapping With Genetic Mapping and Functional Genomics

Genomic scans for signatures of selection allow us to, in principle, detect variants and genes that underlie recent adaptations. By combining selection mapping with genetic mapping of traits known to be relevant to adaptation, we can simultaneously investigate whether genes and variants show signals of recent selection and whether they impact traits that have likely been selected. There are three ways to integrate selection mapping with genetic mapping or functional genomics: (1) To use genetic mapping data from other populations as a form of genome annotation. (2) To perform experimental evolution or artificial selection to be able to study selected variants when they segregate, either by performing genetic mapping before selection or by crossing the selected individuals to some reference population. (3) To perform a comparative study of related populations facing different selection regimes. This short review discusses these different ways of integrating selection mapping with genetic mapping and functional genomics, with examples of how each has been done.

Livestock 2.0 - genome editing for fitter, healthier, and more productive farmed animals

The human population is growing, and as a result we need to produce more food whilst reducing the impact of farming on the environment. Selective breeding and genomic selection have had a transformational impact on livestock productivity, and now transgenic and genome-editing technologies offer exciting opportunities for the production of fitter, healthier and more-productive livestock. Here, we review recent progress in the application of genome editing to farmed animal species and discuss the potential impact on our ability to produce food.

Functional Annotation of Animal Genomes (FAANG): Current Achievements and Roadmap

Functional annotation of genomes is a prerequisite for contemporary basic and applied genomic research, yet farmed animal genomics is deficient in such annotation. To address this, the FAANG (Functional Annotation of Animal Genomes) Consortium is producing genome-wide data sets on RNA expression, DNA methylation, and chromatin modification, as well as chromatin accessibility and interactions. In addition to informing our understanding of genome function, including comparative approaches to elucidate constrained sequence or epigenetic elements, these annotation maps will improve the precision and sensitivity of genomic selection for animal improvement. A scientific community-driven effort has already created a coordinated data collection and analysis enterprise crucial for the success of this global effort. Although it is early in this continuing process, functional data have already been produced and application to genetic improvement reported. The functional annotation delivered by the FAANG initiative will add value and utility to the greatly improved genome sequences being established for domesticated animal species.

A 5-Methylcytosine Site of Growth Differentiation Factor 9 (GDF9) Gene Affects Its Tissue-Specific Expression in Sheep

Simple Summary

Growth differentiation factor 9 (GDF9) is an important gene for ovine fertility. GDF9 is highly expressed in the ovary as opposed to other tissues, but the reason for this is unknown. Our study found this can be caused by the methylation level of the promoter CpG island mC-4 site. This finding contributes to the understanding of the regulatory mechanism of GDF9 gene in reproduction.

Abstract

Growth differentiation factor 9 (GDF9) plays an important role in the early folliculogenesis of sheep. This study investigated the mRNA expression of ovine GDF9 in different tissues by real-time PCR. GDF9 exhibits significantly higher levels of expression (p < 0.01) in the ovary, relative to other tissues, indicating that its expression is tissue specific. To explore the regulatory mechanism of this tissue-specific expression, the methylation level of one CpG island (−1453 to −1854) of GDF9 promoter in ovary and heart was determined. In this region (−1987 to −1750), only the mC-4 site was present in the Sp4 binding site showed differential methylation between the heart and ovary; with increased (p < 0.01) methylation being observed in the heart. Additionally, the methylation level was negatively correlated with GDF9 mRNA expression (R = −0.75, p = 0.012), indicating that the methylation of this site plays an important role in transcriptional regulation of GDF9. The methylation effect of the mC-4 site was confirmed by using dual-luciferase. Site-directed mutation (methylation) of mC-4 site significantly reduced (p < 0.05) basal transcriptional activity of GDF9 promoter in oocytes. These results imply that methylation of GDF9 promoter CpG island mC-4 site may affect the binding of the Sp4 transcription factor to the GDF9 promoter region in sheep, thereby regulating GDF9 expression and resulting in a tissue-specific expression.

ADGRE1 (EMR1, F4/80) Is a Rapidly-Evolving Gene Expressed in Mammalian Monocyte-Macrophages

The F4/80 antigen, encoded by the Adgre1 locus, has been widely-used as a monocyte-macrophage marker in mice, but its value as a macrophage marker in other species is unclear, and has even been questioned. ADGRE1 is a seven transmembrane G protein-coupled receptor with an extracellular domain containing repeated Epidermal Growth Factor (EGF)-like calcium binding domains. Using a new monoclonal antibody, we demonstrated that ADGRE1 is a myeloid differentiation marker in pigs, absent from progenitors in bone marrow, highly-expressed in mature granulocytes, monocytes, and tissue macrophages and induced by macrophage colony-stimulating factor (CSF1) treatment in vivo. Based upon these observations, we utilized RNA-Seq to assess the expression of ADGRE1 mRNA in bone marrow or monocyte-derived macrophages (MDM) and alveolar macrophages from 8 mammalian species including pig, human, rat, sheep, goat, cow, water buffalo, and horse. ADGRE1 mRNA was expressed by macrophages in each species, with inter-species variation both in expression level and response to lipopolysaccharide (LPS) stimulation. Analysis of the RNA-Seq data also revealed additional exons in several species compared to current Ensembl annotations. The ruminant species and horses appear to encode a complete duplication of the 7 EGF-like domains. In every species, Sashimi plots revealed evidence of exon skipping of the EGF-like domains, which are highly-variable between species and polymorphic in humans. Consistent with these expression patterns, key elements of the promoter and a putative enhancer are also conserved across all species. The rapid evolution of this molecule and related ADGRE family members suggests immune selection and a role in pathogen recognition.

Across-Experiment Transcriptomics of Sheep Rumen Identifies Expression of Lipid/Oxo-Acid Metabolism and Muscle Cell Junction Genes Associated With Variation in Methane-Related Phenotypes

Ruminants are significant contributors to the livestock generated component of the greenhouse gas, methane (CH₄). The CH₄ is primarily produced by the rumen microbes. Although the composition of the diet and animal intake amount have the largest effect on CH₄ production and yield (CH₄ production/dry matter intake, DMI), the host also influences CH₄ yield. Shorter rumen feed mean retention time (MRT) is associated with higher dry matter intake and lower CH₄ yield, but the molecular mechanism(s) by which the host affects CH₄ production remain unclear. We integrated rumen wall transcriptome data and CH₄ phenotypes from two independent experiments conducted with sheep in Australia (AUS, n = 62) and New Zealand (NZ, n = 24). The inclusion of the AUS data validated the previously identified clusters and gene sets representing rumen epithelial, metabolic and muscular functions. In addition, the expression of the cell cycle genes as a group was consistently positively correlated with acetate and butyrate concentrations (p < 0.05, based on AUS and NZ data together). The expression of a group of metabolic genes showed positive correlations in both AUS and NZ datasets with CH₄ production (p < 0.05) and yield (p < 0.01). These genes encode key enzymes in the ketone body synthesis pathway and included members of the poorly characterized aldo-keto reductase 1C (AKR1C) family. Several AKR1C family genes appear to have ruminant specific evolution patterns, supporting their specialized roles in the ruminants. Combining differential gene expression in the rumen wall muscle of the shortest and longest MRT AUS animals (no data available for the NZ animals) with correlation and network analysis, we identified a set of rumen muscle genes involved in cell junctions as potential regulators of MRT, presumably by influencing contraction rates of the smooth muscle component of the rumen wall. Higher rumen expression of these genes, including SYNPO (synaptopodin, p < 0.01) and NEXN (nexilin, p < 0.05), was associated with lower CH₄ yield in both AUS and NZ datasets. Unlike the metabolic genes, the variations in the expression of which may reflect the availability of rumen metabolites, the muscle genes are currently our best candidates for causal genes that influence CH₄ yield.

Combination of novel and public RNA-seq datasets to generate an mRNA expression atlas for the domestic chicken

Background

The domestic chicken (Gallus gallus) is widely used as a model in developmental biology and is also an important livestock species. We describe a novel approach to data integration to generate an mRNA expression atlas for the chicken spanning major tissue types and developmental stages, using a diverse range of publicly-archived RNA-seq datasets and new data derived from immune cells and tissues.

Results

Randomly down-sampling RNA-seq datasets to a common depth and quantifying expression against a reference transcriptome using the mRNA quantitation tool Kallisto ensured that disparate datasets explored comparable transcriptomic space. The network analysis tool Graphia was used to extract clusters of co-expressed genes from the resulting expression atlas, many of which were tissue or cell-type restricted, contained transcription factors that have previously been implicated in their regulation, or were otherwise associated with biological processes, such as the cell cycle. The atlas provides a resource for the functional annotation of genes that currently have only a locus ID. We cross-referenced the RNA-seq atlas to a publicly available embryonic Cap Analysis of Gene Expression (CAGE) dataset to infer the developmental time course of organ systems, and to identify a signature of the expansion of tissue macrophage populations during development.

Conclusion

Expression profiles obtained from public RNA-seq datasets – despite being generated by different laboratories using different methodologies – can be made comparable to each other. This meta-analytic approach to RNA-seq can be extended with new datasets from novel tissues, and is applicable to any species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4972-7) contains supplementary material, which is available to authorized users.

Expression patterns and genetic variation of the ovine skeletal muscle transcriptome of sheep from five Spanish meat breeds

The goal of the current study is to analyse the gene expression profile of the ovine skeletal muscle as well as to characterize the genetic variation of transcripts expressed in such tissue. This aim has been achieved by sequencing the longissimus dorsi transcriptomes of 50 sheep distributed in five pools representing the Canaria de Pelo, Roja Mallorquina, Gallega, Xisqueta and Ripollesa Spanish autochthonous breeds. Approximately, 363 million reads per pool have been produced and 71.9–82.9% have been successfully mapped to the ovine genome in a paired-end mode (2 × 75 bp). The 200 most expressed muscle transcripts (≈1% of the total transcript count) account for 51% (Canaria de Pelo) to 67% (Gallega) of the total ovine skeletal muscle mRNA expression. These highly expressed genes play key roles in pathways related with striated muscle contraction, gluconeogenesis, glycolysis, citric acid cycle and respiratory electron transport. RNA-Sequencing of muscle transcripts has also revealed that ~72% of the SNPs detected with this approach are shared by at least two pools, and 10% of them segregate in the five pools under analysis. Most of the substitutions detected by RNA-Seq are synonymous or missense and only a minority are predicted to have consequences on protein function.

Development of in vitro enteroids derived from bovine small intestinal crypts

Cattle are an economically important domestic animal species. In vitro 2D cultures of intestinal epithelial cells or epithelial cell lines have been widely used to study cell function and host-pathogen interactions in the bovine intestine. However, these cultures lack the cellular diversity encountered in the intestinal epithelium, and the physiological relevance of monocultures of transformed cell lines is uncertain. Little is also known of the factors that influence cell differentiation and homeostasis in the bovine intestinal epithelium, and few cell-specific markers that can distinguish the different intestinal epithelial cell lineages have been reported. Here we describe a simple and reliable procedure to establish in vitro 3D enteroid, or "mini gut", cultures from bovine small intestinal (ileal) crypts. These enteroids contained a continuous central lumen lined with a single layer of polarized enterocytes, bound by tight junctions with abundant microvilli on their apical surfaces. Histological and transcriptional analyses suggested that the enteroids comprised a mixed population of intestinal epithelial cell lineages including intestinal stem cells, enterocytes, Paneth cells, goblet cells and enteroendocrine cells. We show that bovine enteroids can be successfully maintained long-term through multiple serial passages without observable changes to their growth characteristics, morphology or transcriptome. Furthermore, the bovine enteroids can be cryopreserved and viable cultures recovered from frozen stocks. Our data suggest that these 3D bovine enteroid cultures represent a novel, physiologically-relevant and tractable in vitro system in which epithelial cell differentiation and function, and host-pathogen interactions in the bovine small intestine can be studied.

Comparative mRNA and miRNA expression in European mouflon (Ovis musimon) and sheep (Ovis aries) provides novel insights into the genetic mechanisms for female reproductive success

Prolific breeds of domestic sheep (Ovis aries) are important genetic resources due to their reproductive performance, which is characterized by multiple lambs per birth and out-of-season breeding. However, the lack of a comprehensive understanding of the genetic mechanisms underlying the important reproductive traits, particularly from the evolutionary genomics perspective, has impeded the efficient advancement of sheep breeding. Here, for the first time, by performing RNA-sequencing we built a de novo transcriptome assembly of ovarian and endometrial tissues in European mouflon (Ovis musimon) and performed an mRNA–miRNA integrated expression profiling analysis of the wild species and a highly prolific domestic sheep breed, the Finnsheep. We identified several novel genes with differentially expressed mRNAs (e.g., EREG, INHBA, SPP1, AMH, TDRD5, and ZP2) between the wild and domestic sheep, which are functionally involved in oocyte and follicle development and fertilization, and are significantly (adjusted P-value < 0.05) enriched in the Gene Ontology (GO) terms of various reproductive process, including the regulation of fertilization, oogenesis, ovarian follicle development, and sperm–egg recognition. Additionally, we characterized 58 differentially expressed miRNAs and 210 associated target genes that are essential for the regulation of female reproduction cycles through specific regulatory networks [e.g., (miR-136, miR-374a, miR-9-5p)-(EREG, INHBA)]. Furthermore, our integrated mRNA and miRNA expression profiling analysis elucidated novel direct and indirect miRNA/mRNA causal regulatory relationships related to the reproductive traits of the Ovis species. This study provides in-depth insights into the genomic evolution underlying the reproductive traits of the Ovis species and valuable resources for ovine genomics.

Cross-species inference of long non-coding RNAs greatly expands the ruminant transcriptome

Background

mRNA-like long non-coding RNAs (lncRNAs) are a significant component of mammalian transcriptomes, although most are expressed only at low levels, with high tissue-specificity and/or at specific developmental stages. Thus, in many cases lncRNA detection by RNA-sequencing (RNA-seq) is compromised by stochastic sampling. To account for this and create a catalogue of ruminant lncRNAs, we compared de novo assembled lncRNAs derived from large RNA-seq datasets in transcriptional atlas projects for sheep and goats with previous lncRNAs assembled in cattle and human. We then combined the novel lncRNAs with the sheep transcriptional atlas to identify co-regulated sets of protein-coding and non-coding loci.

Results

Few lncRNAs could be reproducibly assembled from a single dataset, even with deep sequencing of the same tissues from multiple animals. Furthermore, there was little sequence overlap between lncRNAs that were assembled from pooled RNA-seq data. We combined positional conservation (synteny) with cross-species mapping of candidate lncRNAs to identify a consensus set of ruminant lncRNAs and then used the RNA-seq data to demonstrate detectable and reproducible expression in each species. In sheep, 20 to 30% of lncRNAs were located close to protein-coding genes with which they are strongly co-expressed, which is consistent with the evolutionary origin of some ncRNAs in enhancer sequences. Nevertheless, most of the lncRNAs are not co-expressed with neighbouring protein-coding genes.

Conclusions

Alongside substantially expanding the ruminant lncRNA repertoire, the outcomes of our analysis demonstrate that stochastic sampling can be partly overcome by combining RNA-seq datasets from related species. This has practical implications for the future discovery of lncRNAs in other species.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0391-0) contains supplementary material, which is available to authorized users.

Long noncoding RNAs are associated with metabolic and cellular processes in the jejunum mucosa of pre-weaning calves in response to different diets

Long noncoding RNAs (lncRNAs) emerged as important regulatory component of mechanisms involved in gene expression, chromatin modification and epigenetic processes, but they are rarely annotated in the bovine genome.

Our study monitored the jejunum transcriptome of German Holstein calves fed two different milk diets using transcriptome sequencing (RNA-seq). To identify potential lncRNAs within the pool of unknown transcripts, four bioinformatic lncRNA prediction tools were applied. The intersection of the alignment-free lncRNA prediction tools (CNCI, PLEK and FEELnc) predicted 1,812 lncRNA transcripts concordantly comprising a catalogue of 1,042 putative lncRNA loci expressed in the calves’ intestinal mucosa.

Nine lncRNA loci were differentially expressed (DE lncRNAs) between both calf groups. To elucidate their biological function, we applied a systems biology approach that combines weighted gene co-expression network analysis with functional enrichment and biological pathway analysis. Four DE lncRNAs were found to be strongly correlated with a gene network module (GNM) enriched for genes from canonical pathways of remodeling of epithelial adherens junction, tight junction and integrin signaling. Another DE lncRNA was strongly correlated with a GNM enriched for genes associated with energy metabolism and maintaining of cellular homeostasis with a focus on mitochondrial processes.

Our data suggest that these DE lncRNAs may play potential regulatory roles in modulating biological processes associated with energy metabolism pathways and cellular signaling processes affecting the barrier function of intestinal epithelial cells of calves in response to different feeding regimens in the pre-weaning period.

Beta-hydroxybutyrate infusion identifies acutely differentially expressed genes related to metabolism and reproduction in the hypothalamus and pituitary of castrated male sheep

To identify molecular pathways that couple metabolic imbalances and reproduction, we randomly assigned 10 castrated male sheep to be centrally injected into the lateral ventricle through intracerebroventricular cannulas with 1 ml of β-hydroxybutyric acid sodium salt solution (BHB; 12,800 µmol/l) or saline solution (CON; 0.9% NaCl). Approximately 2 h postinjection, sheep were humanely euthanized, and hypothalamus and pituitary tissues were harvested for transcriptome characterization by RNA sequencing. RNA was extracted from the hypothalamus and pituitary and sequenced at a high depth (hypothalamus: 468,912,732 reads; pituitary: 515,106,092 reads) with the Illumina Hi-Seq 2500 platform and aligned to Bos taurus and Ovis aries genomes. Of the total raw reads, 87% (hypothalamus) and 90.5% (pituitary) mapped to the reference O. aries genome. Within these read sets, ~56% in hypothalamus and 69% in pituitary mapped to either known or putative protein coding genes. Fragments per kilobase of transcripts per million normalized counts were averaged and ranked to identify the transcript expression level. Gene Ontology analysis (DAVID Bioinformatics Resources) was utilized to identify biological process functions related to genes shared between tissues, as well as functional categories with tissue-specific enrichment. Between CON- and BHB-treated sheep, 11 and 44 genes were differentially expressed (adj. P < 0.05) within the pituitary and hypothalamus, respectively. Functional enrichment analyses revealed BHB altered expression of genes in pathways related to stimulus perception, inflammation, and cell cycle control. The set of genes altered by BHB creates a foundation from which to identify the signaling pathways that impact reproduction during metabolic imbalances.

Species-Specific Transcriptional Regulation of Genes Involved in Nitric Oxide Production and Arginine Metabolism in Macrophages

Activated mouse macrophages metabolize arginine via NO synthase (NOS2) to produce NO as an antimicrobial effector. Published gene expression datasets provide little support for the activation of this pathway in human macrophages. Generation of NO requires the coordinated regulation of multiple genes. We have generated RNA-sequencing data from bone marrow-derived macrophages from representative rodent (rat), monogastric (pig and horse), and ruminant (sheep, goat, cattle, and water buffalo) species, and analyzed the expression of genes involved in arginine metabolism in response to stimulation with LPS. In rats, as in mice, LPS strongly induced Nos2, the arginine transporter Slc7a2, arginase 1 (Arg1), GTP cyclohydrolase (Gch1), and argininosuccinate synthase (Ass1). None of these responses was conserved across species. Only cattle and water buffalo showed substantial NOS2 induction. The species studied also differed in expression and regulation of arginase (ARG2, rather than ARG1), and amino acid transporters. Variation between species was associated with rapid promoter evolution. Differential induction of NOS2 and ARG2 between the ruminant species was associated with insertions of the Bov-A2 retrotransposon in the promoter region. Bov-A2 was shown to possess LPS-inducible enhancer activity in transfected RAW264.7 macrophages. Consistent with a function in innate immunity, NO production and arginine metabolism vary greatly between species and differences may contribute to pathogen host restriction.

Classic and atypical scrapie - a genetic perspective

Scrapie was the first prion disease to be recognised and the study of this disease in sheep and goats has provided a wealth of information not only for scrapie but also for the other prion diseases. All prion diseases are under strong genetic control of the prion gene PRNP, independent of whether they are typical or atypical scrapie and which of the different prion strains is causing infection. Decades of studies using experimental disease challenges and field surveys have established disease association models, in which species-specific amino acid variations in the prion or PrP protein, encoded by the PRNP gene, can predict disease susceptibility or resistance. PRNP genetics represents an important and successful basis for implementing scrapie eradication strategies in sheep and goats. In general terms these studies have revealed that there appear to be many more amino acid changes in PrP leading to increased resistance than to higher susceptibility. Most changes are in the globular part of PrP protein and three regions appear to have major influence. This knowledge can be transferred into prion diseases of other species to facilitate genetic control strategies. However, an obstacle remains with the lack of fully understanding the underlying molecular mechanism, impeding our ability to deal with the difference in the genetic control between typical and atypical forms of scrapie or to predict association in newly infected species. This chapter will discuss the advances in both typical and atypical scrapie from a genetic perspective.

A tissue-based approach to selection of reference genes for quantitative real-time PCR in a sheep osteoporosis model

BACKGROUND

In order to better understand the multifactorial nature of osteoporosis, animal models are utilized and compared to healthy controls. Female sheep are well established as a model for osteoporosis induced by ovariectomy, calcium and vitamin D low diet, application of steroids, or a combination of these treatments. Transcriptional studies can be performed by applying quantitative real time PCR (RT-qPCR). RT-qPCR estimates mRNA-levels of target genes in relation to reference genes. A chosen set of reference genes should not show variation under experimental conditions. Currently, no standard reference genes are accepted for all tissue types and experimental conditions. Studies examining reference genes for sheep are rare and only one study described stable reference in mandibular bone. However, this type of bone differs from trabecular bone where most osteoporotic fractures occur. The present study aimed at identifying a set of reference genes for relative quantification of transcriptional activity of ovine spine bone and ovine in vitro differentiated mesenchymal stromal cells (MSC) for reliable comparability.

METHODS

Twelve candidate reference genes belonging to different functional classes were selected and their expression was measured from cultured ovMSCs (n = 18) and ovine bone samples (n = 16), respectively. RefFinder was used to rank the candidate genes.

RESULTS

We identified B2M, GAPDH, RPL19 and YWHAZ as the best combination of reference genes for normalization of RT-qPCR results for transcriptional analyses of these ovine samples.

CONCLUSION

This study demonstrates the importance of applying a set of reference genes for RT-qPCR analysis in sheep. Based on our data we recommend using four identified reference genes for relative quantification of gene expression studies in ovine bone or for in vitro experiments with osteogenically differentiated ovine MSCs.

The genomic architecture of mastitis resistance in dairy sheep

Background

Mastitis is the most prevalent disease in dairy sheep with major economic, hygienic and welfare implications. The disease persists in all dairy sheep production systems despite the implementation of improved management practises. Selective breeding for enhanced mastitis resistance may provide the means to further control the disease. In the present study, we investigated the genetic architecture of four mastitis traits in dairy sheep. Individual animal records for clinical mastitis occurrence and three mastitis indicator traits (milk somatic cell count, total viable bacterial count in milk and the California mastitis test) were collected monthly throughout lactation for 609 ewes of the Greek Chios breed. All animals were genotyped with a custom-made 960-single nucleotide polymorphism (SNP) DNA array based on markers located in quantitative trait loci (QTL) regions for mastitis resistance previously detected in three other distinct dairy sheep populations.

Results

Heritable variation and strong positive genetic correlations were estimated for clinical mastitis occurrence and the three mastitis indicator traits. SNP markers significantly associated with these mastitis traits were confirmed on chromosomes 2, 3, 5, 16 and 19. We identified pathways, molecular interaction networks and functional gene clusters for mastitis resistance. Candidate genes within the detected regions were identified based upon analysis of an ovine transcriptional atlas and transcriptome data derived from milk somatic cells. Relevant candidate genes implicated in innate immunity included SOCS2, CTLA4, C6, C7, C9, PTGER4, DAB2, CARD6, OSMR, PLXNC1, IDH1, ICOS, FYB, and LYFR.

Conclusions

The results confirmed the presence of animal genetic variability in mastitis resistance and identified genomic regions associated with specific mastitis traits in the Chios sheep. The conserved genetic architecture of mastitis resistance between distinct dairy sheep breeds suggests that across-breed selection programmes would be feasible.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3982-1) contains supplementary material, which is available to authorized users.

Integration of quantitated expression estimates from polyA-selected and rRNA-depleted RNA-seq libraries

BACKGROUND

The availability of fast alignment-free algorithms has greatly reduced the computational burden of RNA-seq processing, especially for relatively poorly assembled genomes. Using these approaches, previous RNA-seq datasets could potentially be processed and integrated with newly sequenced libraries. Confounding factors in such integration include sequencing depth and methods of RNA extraction and selection. Different selection methods (typically, either polyA-selection or rRNA-depletion) omit different RNAs, resulting in different fractions of the transcriptome being sequenced. In particular, rRNA-depleted libraries sample a broader fraction of the transcriptome than polyA-selected libraries. This study aimed to develop a systematic means of accounting for library type that allows data from these two methods to be compared.

RESULTS

The method was developed by comparing two RNA-seq datasets from ovine macrophages, identical except for RNA selection method. Gene-level expression estimates were obtained using a two-part process centred on the high-speed transcript quantification tool Kallisto. Firstly, a set of reference transcripts was defined that constitute a standardised RNA space, with expression from both datasets quantified against it. Secondly, a simple ratio-based correction was applied to the rRNA-depleted estimates. The outcome is an almost perfect correlation between gene expression estimates, independent of library type and across the full range of levels of expression.

CONCLUSION

A combination of reference transcriptome filtering and a ratio-based correction can create equivalent expression profiles from both polyA-selected and rRNA-depleted libraries. This approach will allow meta-analysis and integration of existing RNA-seq data into transcriptional atlas projects.