Search

A comprehensive characterization of regulatory elements in the chicken genome across tissues will have substantial impacts on both fundamental and applied research. Here, we systematically identified and characterized regulatory elements in the chicken genome by integrating 377 genome-wide sequencing datasets from 23 adult tissues. In total, we annotated 1.57 million regulatory elements, representing 15 distinct chromatin states, and predicted about 1.2 million enhancer-gene pairs and 7662 super-enhancers. This functional annotation of the chicken genome should have wide utility on identifying regulatory elements accounting for gene regulation underlying domestication, selection, and complex trait regulation, which we explored. In short, this comprehensive atlas of regulatory elements provides the scientific community with a valuable resource for chicken genetics and genomics.

The assessment of genomic conservation between human and pig at the functional level can help understand and improve the potential of pig as a human biomedical model. To address this, we developed aDeeplearning-based approach to learn theGenomicConservation at theFunctional level (DeepGCF) between species by integrating 386 and 374 epigenome and transcriptome profiles from human and pig, respectively. DeepGCF demonstrated a better prediction performance compared to the previous functional conservation prediction method. In addition, we showed that the resulting DeepGCF score captures the functional conservation by examining DeepGCF on chromatin states, sequence ontologies, and regulatory variants. Regions with higher DeepGCF score play a more important role in regulatory activities and show heritability enrichment in human complex traits and diseases. Our DeepGCF approach shows a promising application on the comparison of cross-species functional conservation, and the model framework can be easily adapted to other species. By expanding the model to integrate the functional profiles of multiple species, including human, mouse, pig, cattle, and other livestock animals in the future, the functional conservation information will provide additional insight into the genetic and evolutionary mechanisms behind complex traits and diseases.

The Farm animal Genotype-Tissue Expression (FarmGTEx, https://www.farmgtex.org/) project has been established to develop a comprehensive public resource of genetic regulatory variants in domestic animal species, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biology discovery and exploitation in animal breeding and human biomedicine. Here we present results from the pilot phase of PigGTEx (http://piggtex.farmgtex.org/), where we processed 9,530 RNA-sequencing and 1,602 whole-genome sequencing samples from pigs. We build a pig genotype imputation panel, characterize the transcriptional landscape across over 100 tissues, and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We study interactions between genotype and breed/cell type, evaluate tissue specificity of regulatory effects, and elucidate the molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying about 80% of the genetic associations for 207 pig complex phenotypes, and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, corroborating the importance of pigs as a human biomedical model.

The objective of this study was to assess the genetic diversity and phylogenetic relationship of nine sheep populations, including two famous high prolific populations and seven popular mutton populations raised in China. Overall, these sheep populations in this study exhibited a rich genetic diversity. Both the expected heterozygosity and Nei's unbiased gene diversity ranged from 0.64 to 0.75, with the lowest value found in Dorset sheep (DST) and the highest in Hu sheep (HUS) and Ba Han sheep (BAS). The polymorphic information content (PIC) varied between 0.59 in DST and 0.71 in HUS and BAS. Specifically, for individual breeds, the small-tail Han sheep (STH) and the four introduced populations did not display the expected diversity; therefore more attention should be paid to the maintenance of diversity during management of these populations. The results of un-weighted pair-group method (UPGMA) phylogenetic tree and structure analysis indicated that the nine investigated populations can be divided into two groups. Suffolk (SUF) and DST were clustered in one group, and the other group can be further divided into three clusters: German Mutton Merino (GMM)–BAS–Bamei Mutton sheep (BAM), HUS–STH and Du Han (DOS)–Dorper (DOP). This clustering result is consistent with sheep breeding history. TreeMix analysis also hinted at the possible gene flow from GMM to SUF. Together, an in-depth view of genetic diversity and genetic relationship will have important implications for breed-specific management.

Transposable elements (TEs) are an extensive source of genetic polymorphisms and play an indispensable role in chromatin architecture, transcriptional regulatory networks, and genomic evolution. The pig is an important source of animal protein and serves as a biomedical model for humans, yet the functional role of TEs in pigs and their contributions to complex traits are largely unknown. Here, we built a comprehensive catalog of TEs (n = 3,087,929) in pigs by a newly developed pipeline. Through integrating multi-omics data from 21 tissues, we found that SINEs with different ages were significantly associated with genomic regions with distinct functions across tissues. The majority of young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, the expression of transcripts that were derived from the remaining active young SINEs exhibited strong tissue specificity through cross-examining 3,570 RNA-seq from 79 tissues and cell types. Furthermore, we detected 211,067 polymorphic SINEs (polySINEs) in 374 individuals genome-wide and found that they clearly recapitulated known patterns of population admixture in pigs. Out of them, 340 population-specific polySINEs were associated with local adaptation. Mapping these polySINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs in functional genomics and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

Gene regulatory elements are central drivers of phenotypic variation and thus of critical importance towards understanding the genetics of complex traits. The Functional Annotation of Animal Genomes consortium was formed to collaboratively annotate the functional elements in animal genomes, starting with domesticated animals. Here we present an expansive collection of datasets from eight diverse tissues in three important agricultural species: chicken (Gallus gallus), pig (Sus scrofa), and cattle (Bos taurus). Comparative analysis of these datasets and those from the human and mouse Encyclopedia of DNA Elements projects reveal that a core set of regulatory elements are functionally conserved independent of divergence between species, and that tissue-specific transcription factor occupancy at regulatory elements and their predicted target genes are also conserved. These datasets represent a unique opportunity for the emerging field of comparative epigenomics, as well as the agricultural research community, including species that are globally important food resources., In order to interpret non-coding variants, information about regulatory elements in the genome is essential. Here, the authors annotate regulatory elements in chicken, pig and cattle, and characterize conservation of these elements between species.

The functional annotation of livestock genomes is crucial for understanding the molecular mechanisms that underpin complex traits of economic importance, adaptive evolution and comparative genomics. Here, we provide the most comprehensive catalogue to date of regulatory elements in the pig (Sus scrofa) by integrating 223 epigenomic and transcriptomic data sets, representing 14 biologically important tissues. We systematically describe the dynamic epigenetic landscape across tissues by functionally annotating 15 different chromatin states and defining their tissue-specific regulatory activities. We demonstrate that genomic variants associated with complex traits and adaptive evolution in pig are significantly enriched in active promoters and enhancers. Furthermore, we reveal distinct tissue-specific regulatory selection between Asian and European pig domestication processes. Compared with human and mouse epigenomes, we show that porcine regulatory elements are more conserved in DNA sequence, under both rapid and slow evolution, than those under neutral evolution across pig, mouse, and human. Finally, we provide novel biological insights on tissue-specific regulatory conservation and demonstrate that, depending on the traits, mouse or pig might be more appropriate biomedical models for different complex traits and diseases in humans through integrating comparative epigenomes with 47 human genome-wide association studies.

Simple Summary Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. There is large variation in the reproductive abilities of Bamei mutton sheep. After recent artificial selection, the average lambing rate of the Bamei mutton nucleus group was over 150%. We used the FST (Fixation Index) and XP-EHH (The Cross-Population Extended Haplotype Homozygosity) statistical approach to detect the selective sweeps between high- and low-fecundity Bamei mutton sheep groups. JUN (JUN proto-oncogene, AP-1 transcription factor subunit), ITPR3 (inositol 1,4,5-trisphosphate receptor type 3, PLCB2 (phospholipase C beta 2), HERC5 (HECT and RLD domain containing E3 ubiquitin protein ligase 5), and KDM4B (lysine demethylase 4B) were detected that are potential responsible for litter size. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. Abstract Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. Here, we focused on detecting candidate genes associated with the increasing of the litter size in this breeds under recent artificial selection to improve the efficiency of mutton production. We selected five high- and five low-fecundity Bamei mutton sheep for whole-genome resequencing to identify candidate genes for sheep prolificacy. We used the FST and XP-EHH statistical approach to detect the selective sweeps between these two groups. Combining the two selective sweep methods, the reproduction-related genes JUN, ITPR3, PLCB2, HERC5, and KDM4B were detected. JUN, ITPR3, and PLCB2 play vital roles in GnRH (gonadotropin-releasing hormone), oxytocin, and estrogen signaling pathway. Moreover, KDM4B, which had the highest FST value, exhibits demethylase activity. It can affect reproduction by binding the promoters of estrogen-regulated genes, such as FOXA1 (forkhead box A1) and ESR1 (estrogen receptor 1). Notably, one nonsynonymous mutation (p.S936A) specific to the high-prolificacy group was identified at the TUDOR domain of KDM4B. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. The identified genomic regions that are responsible for litter size can in turn be used for further selection.

Background: Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. Here, we focused on increasing the litter size of this breed to improve the efficiency of mutton production.Results: We selected five high- and five low-fecundity Bamei mutton sheep for whole-genome resequencing to identify candidate genes for sheep prolificacy. We used the FST and XP-EHH statistical approach to detect the selective sweeps between these two groups. Combining the two selective sweep methods, the reproduction-related genes JUN, ITPR3, PLCB2, HERC5, and KDM4B were detected. JUN, ITPR3, and PLCB2 play vital roles in GnRH, oxytocin, and estrogen signaling pathway. Moreover, KDM4B, which had the highest FST value, exhibits demethylase activity. It can affect reproduction by binding the promoters of estrogen-regulated genes, such as FOXA1 and ESR1. Notably, one nonsynonymous mutation (p.S936A) specific to the high-prolificacy group was identified at the TUDOR domain of KDM4B.Conclusion: These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. The identified genomic regions of KDM4B that are responsible for litter size can in turn be used for further selection.

The objective of this study was to assess the genetic diversity and population structure of nine sheep populations (four famous mutton breeds worldwide, three popular mutton populations and two most famous high prolific breeds in China) and predict the heterosis among these populations based on microsatellite markers. Results showed the relatively rich genetic diversity in Chinese sheep populations. The results of UPGMA phylogenetic tree and Structure analysis indicated that the nine investigated populations can be divided into two groups. Suffolk and Dorset were clustered in one group, and the other group can be further divided into three clusters: GMM-BAS-BAM, HUS-STH and DOS-DOP. This is consistent with sheep breeding history. TreeMix analysis also hinted the possible gene flow from GMM to SUF. Additionally, based on the calculation of genetic distance, the optimal cross combinations are recommended for two popular mutton sheep populations (BAM and BAS) in the southern Mongolian Plateau and the most famous high prolific breeds (Small Tail Han sheep and Hu sheep) in China, respectively. Together, these results bring us an in-depth view of genetic diversity within sheep population and the genetic relationship among populations.

Figure S1. Evidence of selective sweep at chromosome 15 near PDGFD. Plots of selective sweep statistics at chromosome 15, from top to bottom: (1) LSBL; (2) pair-wise genetic distance dxy; (3) intra-lineage heterozygosity HP (standardized). Figure S2. Validation of LOC101117953 DNA sequence. (A) Electrophoresis of PCR product captured by LOC101117953-specific primers. (B) Comparison between the result of PCR product sequencing and the reference assembly of sheep (oviAri3). Figure S3. Dot-plots for alignment between eight gene sequences and the sheep genome at chromosome 13. Alignments generated by LASTZ were visualized as dot plots using R script. Exon positions on the normal paralog are showed in blue boxes. Figure S4. Features of different tail types in the hybrid population. Pictures taken for two individual sheep from the hybrid population with thin tail (left) and fat tail (right). Figure S5. Tissue-expression patterns of genes near IBH region. Primers were designed to capture cDNA of five genes from ovine transcriptome. M, marker; 1, pituitary; 2, hypothalamus; 3, cerebellum; 4, cerebrum; 5, ovary; 6, oviduct; 7, cornua uterus; 8, corpus uterus; 9, thyroid; 10, adrenal gland; 11, heart; 12, liver; 13, spleen; 14, lung; 15, kidney; 16, perirenal adipose; 17, tail adipose. Figure S6. Stop-gain mutation on the putative protein sequence of LOC101117953. The position and type of the mutation which truncates the putative protein encoded by LOC101117953 (or r-PPP1CC). The truncated codon is highlighted in red. Figure S7. Chromatin interactions detected in human ESC. Sheep IBH region is homologous to the displayed area on human chromosome 20, between BMP2 and HAO1. Interaction density is visualized by the 3D Genome Browser ( http://promoter.bx.psu.edu/hi-c/ ), using H1-ESC data and a resolution of 10â K. (DOC 1839 kb)

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 &lt, 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 (&minus, 1453 to &minus, 1854) of GDF9 promoter in ovary and heart was determined. In this region (&minus, 1987 to &minus, 1750), only the mC-4 site was present in the Sp4 binding site showed differential methylation between the heart and ovary, with increased (p &lt, 0.01) methylation being observed in the heart. Additionally, the methylation level was negatively correlated with GDF9 mRNA expression (R = &minus, 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 &lt, 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.

DIO2 was studied as a candidate gene for reproduction in sheep. Three pairs of primers were designed to detect single nucleotide polymorphisms (SNPs) of exon 2 of DIO2 in nonseasonal and seasonal reproduction breeds. The polymorphisms were detected using SeqMan II and analyzed by PCR-RFLP. The sheep DIO2 gene was cloned and its mutations were detected in five sheep breeds whose reproductive seasonality and litter sizes were different. Association analyses were performed to reveal the relationships between SNPs and reproductive seasonality as well as litter size in sheep. The results showed that a 1041-bp DNA fragment of the sheep DIO2 gene was obtained and a part of the coding sequence of 591 bp encodes 196 amino acids, having high homology with that of other mammals. Moreover, three mutations (G318T, G482C, and C885T) were detected. The test of genotype distribution in five sheep breeds showed that there were significant differences between nonseasonal and seasonal breeds for 318 and 482 loci. However, the associations between the three polymorphisms and sheep seasonal reproduction were not very distinct. Additionally, there was no significant associations between polymorphisms and litter sizes in Small Tail Han sheep. These results may provide the basic theory for further study.

Background Sheep have developed the ability to store fat in their tails, which is a unique way of reserving energy to survive a harsh environment. However, the mechanism underlying this adaptive trait remains largely unsolved. Results In the present study, we provide evidence for the genetic determinants of fat tails, based on whole genome sequences of 89 individual sheep. A genome-wide scan of selective sweep identified several candidate loci including a region at chromosome 13, a haplotype of which underwent rapid evolution and spread through fat-tailed populations in China and the Middle East. Sequence analysis revealed an inter-genic origin of this locus, which later became a hotspot of ruminant-specific retro-transposon named BovB. Additionally, the candidate locus was validated based on a fat- and thin-tailed cross population. The expression of an upstream gene BMP2 was differentially regulated between fat-tailed and thin-tailed individuals in tail adipose and several other tissue types. Conclusions Our findings suggest the fixation of fat tails in domestic sheep is caused by a selective sweep near a retro-transposable hotspot at chromosome 13, the diversity of which specifically affects the expression of BMP2. The present study has shed light onto the understanding of fat metabolism. Electronic supplementary material The online version of this article (10.1186/s12864-019-5620-6) contains supplementary material, which is available to authorized users.

Simple Summary Litter size is one of the most important reproductive traits in sheep. Two single nucleotide polymorphisms (SNPs), g.71874104G>A and g.71833755T>C, in the Histone Cell Cycle Regulator (HIRA) gene, were identified by whole-genome sequencing (WGS) and may be correlated with litter size in sheep. The two SNPs were genotyped and expression patterns of HIRA was determined in sheep breeds with different fecundity and in groups of Small Tail Han sheep producing large or small litters. Association analysis indicated that both SNPs were significantly correlated with litter size in Small Tail Han sheep. Furthermore, high levels of HIRA expression may have a negative effect on litter size in Small Tail Han sheep. Abstract Maintenance of appropriate levels of fecundity is critical for efficient sheep production. Opportunities to increase sheep litter size include identifying single gene mutations with major effects on ovulation rate and litter size. Whole-genome sequencing (WGS) data of 89 Chinese domestic sheep from nine different geographical locations and ten Australian sheep were analyzed to detect new polymorphisms affecting litter size. Comparative genomic analysis of sheep with contrasting litter size detected a novel set of candidate genes. Two SNPs, g.71874104G>A and g.71833755T>C, were genotyped in 760 Small Tail Han sheep and analyzed for association with litter size. The two SNPs were significantly associated with litter size, being in strong linkage disequilibrium in the region 71.80–71.87 Mb. This haplotype block contains one gene that may affect litter size, Histone Cell Cycle Regulator (HIRA). HIRA mRNA levels in sheep with different lambing ability were significantly higher in ovaries of Small Tail Han sheep (high fecundity) than in Sunite sheep (low fecundity). Moreover, the expression levels of HIRA in eight tissues of uniparous Small Tail Han sheep were significantly higher than in multiparous Small Tail Han sheep (p < 0.05). HIRA SNPs significantly affect litter size in sheep and are useful as genetic markers for litter size.

Background Animal domestication has been extensively studied, but the process of feralization remains poorly understood. Results Here, we performed whole-genome sequencing of 99 sheep and identified a primary genetic divergence between 2 heterogeneous populations in the Tibetan Plateau, including 1 semi-feral lineage. Selective sweep and candidate gene analysis revealed local adaptations of these sheep associated with sensory perception, muscle strength, eating habit, mating process, and aggressive behavior. In particular, a horn-related gene, RXFP2, showed signs of rapid evolution specifically in the semi-feral breeds. A unique haplotype and repressed horn-related tissue expression of RXFP2 were correlated with higher horn length, as well as spiral and horizontally extended horn shape. Conclusions Semi-feralization has an extensive impact on diverse phenotypic traits of sheep. By acquiring features like those of their wild ancestors, semi-feral sheep were able to regain fitness while in frequent contact with wild surroundings and rare human interventions. This study provides a new insight into the evolution of domestic animals when human interventions are no longer dominant.

Litter size is a favorable economic trait for the goat industry, but remains a complex trait controlled by multiple genes in multiple organs. Several genes have been identified that may affect embryo survival, follicular development, and the health of fetuses during pregnancy. Jining Grey goats demonstrate the largest litter size among goat breeds indigenous to China. In order to better understand the genetic basis of this trait, six suppression subtractive hybridization (SSH) cDNA libraries were constructed using pooled mRNAs from hypothalamuses, pituitaries, and ovaries of sexually mature and adult polytocous Jining Grey goats, as testers, versus the pooled corresponding mRNAs of monotocous Liaoning Cashmere goats, as drivers. A total of 1,458 true-positive clones--including 955 known genes and 481 known and 22 unknown expressed sequence tags--were obtained from the SSH libraries by sequencing and alignment. The known genes were categorized into cellular processes and signaling information storage and processing, and metabolism. Three genes (FTH1, GH, and SAA) were selected to validate the SSH results by quantitative real-time PCR; all three were up-regulated in the corresponding tissues in the tester group indicating that these are candidate genes associated with the large litter size of Jining Grey goats. Several other identified genes may affect embryo survival, follicular development, and health during pregnancy. This study provides insights into the mechanistic basis by which the caprine hypothalamic-pituitary-gonadal axis affects reproductive traits and provides a theoretical basis for goat production and breeding.

The existence of major genes affecting fecundity in sheep flocks throughout the world has been demonstrated. Three major genes whose mutations can increase ovulation rate have been discovered, and all related to the transforming growth factor β (TGF-β) superfamily. The mutant FecB of bone morphogenetic protein receptor 1B (BMPR1B) has an additive effect on ovulation rate. Six mutations (FecXI, FecXH, FecXG, FecXB, FecXL, FecXR) of bone morphogenetic protein 15 (BMP15) related with fertility have been identified that share the same mechanism. All the mutants can increase ovulation rate in heterozygotes and cause complete sterility in homozygotes. Homozygous ewes with two new mutations (FecXGr, FecXO) of BMP15 had increased ovulation rate without causing sterility. There are five mutations in growth differentiation factor 9 (GDF9) associated with sheep prolificacy where FecGE and FecGF have additive an effect on ovulation rate and litter size. The newly identified β-1,4-N-acetylgalactosaminyltransferase 2 (B4GALNT2) gene of FecL is proposed as a new mechanism of ovulation rate regulation in sheep. Woodlands is an X-linked maternally imprinted gene which increases ovulation rate. In addition, several putative major genes need to be verified. This review is focused on the identification of the mutations and mechanisms whereby the major genes affecting ovulation rate.

The aims of study were to reveal the changes in genetic diversity in reproduction selection of sheep. Genetic polymorphisms of microsatellite GC101 locus and its linkage with the FecB gene were detected in Chinese local sheep (Small Tail Han and Hu sheep) and western sheep (Texel and Dorset) breeds. There was the same FecB mutation (A746G) in both Small Tail Han and Hu sheep as that in Booroola Merino ewes, but not in Texel and Dorset sheep. Seven alleles, from 192 bp to 210 bp, and 20 genotypes were detected at the GC101 locus in 447 individuals totally. The preponderant allele was 200 bp for Chinese sheep and 202 bp for western sheep.? Linkage analysis showed that there was strong linkage disequilibrium between the 200-bp allele of the GC101 locus and the B allele of the FecB gene (r2 = 0.639). These results preliminarily implied that both the 200-bp allele of the GC101 locus and the B allele of FecB had been reinforced in reproduction selection of Chinese sheep. The differentiation of GC101 allele distribution between Chinese sheep and western sheep also indicated that the FecB of Chinese sheep might be derived from Chinese local areas.