Expression characteristics of pineal miRNAs at ovine different reproductive stages and the identification of miRNAs targeting the AANAT gene.

BACKGROUND: Many recent studies have shown that miRNAs play important roles in the regulation of animal reproduction, including seasonal reproduction. The pineal gland is a crucial hub in the regulation of seasonal reproduction. However, little is known about the expression characteristics of pineal miRNAs in different reproductive seasons (anestrus and breeding season). Therefore, the expression profiles and regulatory roles of ovine pineal miRNAs were investigated during different reproductive stages using Solexa sequencing technology and dual luciferase reporter assays. RESULTS: A total of 427 miRNAs were identified in the sheep pineal gland. Significant differences in miRNA expression were demonstrated between anestrus and the breeding season in terms of the frequency distributions of miRNA lengths, number of expressed miRNAs, and specifically and highly expressed miRNAs in each reproductive stage. KEGG analysis of the differentially expressed (DE) miRNAs between anestrus and the breeding season indicated that they are significantly enriched in pathways related to protein synthesis, secretion and uptake. Furthermore, transcriptome analysis revealed that many target genes of DE miRNAs in the ribosome pathway showed relatively low expression in the breeding season. On the other hand, analyses combining miRNA-gene expression data with target relationship validation in vitro implied that miR-89 may participate in the negative regulation of aralkylamine N-acetyltransferase (AANAT) mRNA expression by targeting its 3'UTR at a unique binding site. CONCLUSIONS: Our results provide new insights into the expression characteristics of sheep pineal miRNAs at different reproductive stages and into the negative regulatory effects of pineal miRNAs on AANAT mRNA expression.

Rapid evolution of a retro-transposable hotspot of ovine genome underlies the alteration of BMP2 expression and development of fat tails.

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.

Selection signature in domesticated animals.

Domesticated animals play an important role in the life of humanity. All these domesticated animals undergo same process, first domesticated from wild animals, then after long time natural and artificial selection, formed various breeds that adapted to the local environment and human needs. In this process, domestication, natural and artificial selection will leave the selection signal in the genome. The research on these selection signals can find functional genes directly, is one of the most important strategies in screening functional genes. The current studies of selection signal have been performed in pigs, chickens, cattle, sheep, goats, dogs and other domestic animals, and found a great deal of functional genes. This paper provided an overview of the types and the detected methods of selection signal, and outlined researches of selection signal in domestic animals, and discussed the key issues in selection signal analysis and its prospects.

Identification and verification of differentially expressed genes in the caprine hypothalamic-pituitary-gonadal axis that are associated with litter size.

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.

Expression, Polymorphism, and Potential Functional Sites of the BMPR1A Gene in the Sheep Horn.

Sheep horns are composed of bone and sheaths, and the BMPR1A gene is required for cartilage and osteogenic differentiation. Therefore, the BMPR1A gene may have a function related to the sheep horn, but its relationship with the sheep horn remains unclear. In this study, we first utilized RNA sequencing (RNA-seq) data to investigate the expression of the BMPR1A gene in different tissues and breeds of sheep. Second, whole-genome sequencing (WGS) data were used to explore the functional sites of the BMPR1A gene. Lastly, the allele-specific expression of the BMPR1A gene was explored. Our results indicate that BMPR1A gene expression is significantly higher in the normal horn groups than in the scurred groups. Importantly, this trend is consistent across several sheep breeds. Therefore, this finding suggests that the BMPR1A gene may be related to horn type. A total of 43 Single-Nucleotide Polymorphisms (SNPs) (F-statistics > 0.15) and 10 allele-specific expressions (ASEs) exhibited difference between the large and small horn populations. It is probable that these sites significantly impact the size of sheep horns. Compared to other polled species, we discovered ten amino acid sites that could influence horn presence. By combining RNA-seq and WGS functional loci results, we identified a functional site at position 40574836 on chromosome 25 that is both an SNP and exhibits allele-specific expression. In conclusion, we demonstrated that the BMPR1A gene is associated with horn type and identified some important functional sites which can be used as molecular markers in the breeding of sheep horns.

Mapping and functional characterization of structural variation in 1060 pig genomes.

BACKGROUND: Structural variations (SVs) have significant impacts on complex phenotypes by rearranging large amounts of DNA sequence. RESULTS: We present a comprehensive SV catalog based on the whole-genome sequence of 1060 pigs (Sus scrofa) representing 101 breeds, covering 9.6% of the pig genome. This catalog includes 42,487 deletions, 37,913 mobile element insertions, 3308 duplications, 1664 inversions, and 45,184 break ends. Estimates of breed ancestry and hybridization using genotyped SVs align well with those from single nucleotide polymorphisms. Geographically stratified deletions are observed, along with known duplications of the KIT gene, responsible for white coat color in European pigs. Additionally, we identify a recent SINE element insertion in MYO5A transcripts of European pigs, potentially influencing alternative splicing patterns and coat color alterations. Furthermore, a Yorkshire-specific copy number gain within ABCG2 is found, impacting chromatin interactions and gene expression across multiple tissues over a stretch of genomic region of ~200 kb. Preliminary investigations into SV's impact on gene expression and traits using the Pig Genotype-Tissue Expression (PigGTEx) data reveal SV associations with regulatory variants and gene-trait pairs. For instance, a 51-bp deletion is linked to the lead eQTL of the lipid metabolism regulating gene FADS3, whose expression in embryo may affect loin muscle area, as revealed by our transcriptome-wide association studies. CONCLUSIONS: This SV catalog serves as a valuable resource for studying diversity, evolutionary history, and functional shaping of the pig genome by processes like domestication, trait-based breeding, and adaptive evolution.

A compendium of genetic regulatory effects across pig tissues.

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

A body map of super-enhancers and their function in pig.

Introduction: Super-enhancers (SEs) are clusters of enhancers that act synergistically to drive the high-level expression of genes involved in cell identity and function. Although SEs have been extensively investigated in humans and mice, they have not been well characterized in pigs. Methods: Here, we identified 42,380 SEs in 14 pig tissues using chromatin immunoprecipitation sequencing, and statistics of its overall situation, studied the composition and characteristics of SE, and explored the influence of SEs characteristics on gene expression. Results: We observed that approximately 40% of normal enhancers (NEs) form SEs. Compared to NEs, we found that SEs were more likely to be enriched with an activated enhancer and show activated functions. Interestingly, SEs showed X chromosome depletion and short interspersed nuclear element enrichment, implying that SEs play an important role in sex traits and repeat evolution. Additionally, SE-associated genes exhibited higher expression levels and stronger conservation than NE-associated genes. However, genes with the largest SEs had higher expression levels than those with the smallest SEs, indicating that SE size may influence gene expression. Moreover, we observed a negative correlation between SE gene distance and gene expression, indicating that the proximity of SEs can affect gene activity. Gene ontology enrichment and motif analysis revealed that SEs have strong tissue-specific activity. For example, the CORO2B gene with a brain-specific SE shows strong brain-specific expression, and the phenylalanine hydroxylase gene with liver-specific SEs shows strong liver-specific expression. Discussion: In this study, we illustrated a body map of SEs and explored their functions in pigs, providing information on the composition and tissue-specific patterns of SEs. This study can serve as a valuable resource of gene regulatory and comparative analyses to the scientific community and provides a theoretical reference for genetic control mechanisms of important traits in pigs.

Genetic diversity, tissue-specific expression, and functional analysis of the ATP7A gene in sheep.

In humans, variation of the ATP7A gene may cause cranial exostosis, which is similar to "human horn," but the function of the ATP7A gene in sheep is still unknown. Tissue expression patterns and potential functional loci analysis of the ATP7A gene could help understand its function in sheep horn. In this study, we first identified tissue, sex, breed, and species-specific expression of the ATP7A gene in sheep based on the RNA-sequencing (RNA-seq) data. Second, the potential functional sites of the ATP7A gene were analyzed by using the whole genome sequencing (WGS) data of 99 sheep from 10 breeds. Last, the allele-specific expression of the ATP7A gene was explored. Our result showed the ATP7A gene has significantly higher expression in the big horn than in the small horn, and the ATP7A gene has high expression in the horn and skin, suggesting that this gene may be related to the horn. The PCA results show that the region around the ATP7A can distinguish horned and hornless groups to some extent, further indicating that the ATP7A may be related to horns. When compared with other species, we find seven ruminate specific amino acid sites of the ATP7A protein, which can be important to the ruminate horn. By analyzing WGS, we found 6 SNP sites with significant differences in frequency in horned and hornless populations, and most of these variants are present in the intron. But we still find some potential functional sites, including three missenses, three synonymous mutations, and four Indels. Finally, by combining the RNA-seq and WGS functional loci results, we find three mutations that showed allele-specific expression between big and small horns. This study shows that the ATP7A gene in sheep may be related to horn size, and several potential functional sites we identified here can be useful molecular markers for sheep horn breeding.

Association Analyses between Single Nucleotide Polymorphisms in ZFAT, FBN1, FAM184B Genes and Litter Size of Xinggao Mutton Sheep.

Previous studies have screened key candidate genes for litter size in sheep, including fibrillin-1 (FBN1), family with sequence similarity 184 member B (FAM184B) and zinc finger and AT-hook domain containing (ZFAT). Therefore, it is necessary to verify these genes in the Xinggao mutton sheep population and determine the associated loci for litter size. In this study, three loci (FBN1 g.160338382 T > C, FAM184B g.398531673 C > T and ZFAT g.20150315 C > T) were firstly screened based on the population differentiation coefficient between the polytocous and monotocous sheep groups. Then, population genetic analysis and association analysis were performed on these loci. The results revealed that the g.160338382 T > C in FBN1 was significantly associated with the litter size of sheep. Moreover, there was no significant interaction effect between the g.160338382 T > C locus and FecB on litter size. Notably, g.160338382 T > C is adjacent to the anterior border of exon 58 and belongs to a splice polypyrimidine tract variant, which may lead to alternative splicing and ultimately cause changes in the structure and function of the protein. In summary, our results provided a potentially effective genetic marker for improving the litter size of sheep.

Learning functional conservation between human and pig to decipher evolutionary mechanisms underlying gene expression and complex traits.

Assessment of genomic conservation between humans and pigs at the functional level can improve the potential of pigs as a human biomedical model. To address this, we developed a deep learning-based approach to learn the genomic conservation at the functional level (DeepGCF) between species by integrating 386 and 374 functional profiles from humans and pigs, respectively. DeepGCF demonstrated better prediction performance compared with the previous method. In addition, the resulting DeepGCF score captures the functional conservation between humans and pigs by examining chromatin states, sequence ontologies, and regulatory variants. We identified a core set of genomic regions as functionally conserved that plays key roles in gene regulation and is enriched for the heritability of complex traits and diseases in humans. Our results highlight the importance of cross-species functional comparison in illustrating the genetic and evolutionary basis of complex phenotypes.

Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits.

Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs 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 and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

An atlas of regulatory elements in chicken: A resource for chicken genetics and genomics.

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 effect of mutation at M307 in FUT1 gene on susceptibility of Escherichia coli F18 and gene expression in Sutai piglets.

Escherichia coli F18 (ECF18) is a common porcine enteric pathogen. The pathogenicity of ECF18 bacteria depends on the existence of ECF18 receptor in the brush border membranes of piglet's small intestinal mucosa. Alpha (1) fucosyltransferase gene (FUT1) has been identified as the candidate gene controlling the adhesion to ECF18 receptor. The genetic variations in the position of M307 nucleotide in open reading frame of FUT1 have been proposed as a marker for selecting resistant pigs. The piglets were divided into three groups, AA, AG and GG, according to the genotypes present at M307 of FUT1. Small intestinal epithelium cells of piglets with AA, AG and GG genotypes were selected to test the adhesion capability of the wild type E.coli expressing F18ab fimbriae, the recombinant E. coli expressing F18ac fimbriae or the recombinant E. coli secreting and surface-displaying the FedF subunit of F18ab fimbriae, respectively. Here, we examined the distribution and expression of porcine FUT1 mRNA in different tissues in Sutai pigs using real-time PCR. The results showed that piglets with AA genotype show resistance, whereas piglets with GG or AG genotypes are sensitive to the pathogenic E. coli F18 in Sutai piglets. FUT1 was expressed in all the tissues that were examined, and the gene's expression was highest in the lungs. There was no significant difference in expression level among the three genotypes in the liver, lung, stomach and duodenum, where the gene expression was relatively high. The present analysis suggested that mutation at M307 in FUT1 gene determines susceptibility of small intestinal epithelium to E. coli F18 adhesion in Sutai piglet and the expression of FUT1 gene may be regulated by other factors or the mutation was likely to be in linkage disequilibrium with some cis-regulatory variants.

Evaluation of M307 of FUT1 gene as a genetic marker for disease resistance breeding of Sutai pigs.

Alpha (1,2) fucosyltransferase (FUT1) gene has been identified as a candidate gene for controlling the expression of the receptor for ETEC F18. The genetic variations in the position of M307 nucleotide in open reading frame of FUT1 have been proposed as a marker for selecting ETEC F18 resistant pigs. The polymorphisms of M307 in FUT1 of breeding base group for ETEC F18 resistance of Sutai pigs (Duroc × Meishan) was detected and their correlations to some immune indexes, growth and development ability, carcass traits and meat quality were also analyzed, which aimed to investigate feasibility of further breeding for diseases resistance based on M307 of FUT1 for Sutai pigs. After digested by Hin6 I, M307 of FUT1 gene could be divided into three kinds of genotypes, AA, AG, and GG. The frequencies were 0.235, 0.609, and 0.156, respectively. The results indicated that Sutai pigs with the AA genotype in M307 of FUT1 gene not only have relatively strong general disease resistance ability in piglets, but also have higher growth and development ability and stable carcass traits and meat quality. It is entirely feasible to raise the new strains of Sutai pigs resistant to Escherichia coli F18 based on genetic marker of the M307 position in FUT1gene.

DNA methylation may affect beef tenderness through signal transduction in Bos indicus.

BACKGROUND: Beef tenderness is a complex trait of economic importance for the beef industry. Understanding the epigenetic mechanisms underlying this trait may help improve the accuracy of breeding programs. However, little is known about epigenetic effects on Bos taurus muscle and their implications in tenderness, and no studies have been conducted in Bos indicus. RESULTS: Comparing methylation profile of Bos indicus skeletal muscle with contrasting beef tenderness at 14 days after slaughter, we identified differentially methylated cytosines and regions associated with this trait. Interestingly, muscle that became tender beef had higher levels of hypermethylation compared to the tough group. Enrichment analysis of predicted target genes suggested that differences in methylation between tender and tough beef may affect signal transduction pathways, among which G protein signaling was a key pathway. In addition, different methylation levels were found associated with expression levels of GNAS, PDE4B, EPCAM and EBF3 genes. The differentially methylated elements correlated with EBF3 and GNAS genes overlapped CpG islands and regulatory elements. GNAS, a complex imprinted gene, has a key role on G protein signaling pathways. Moreover, both G protein signaling pathway and the EBF3 gene regulate muscle homeostasis, relaxation, and muscle cell-specificity. CONCLUSIONS: We present differentially methylated loci that may be of interest to decipher the epigenetic mechanisms affecting tenderness. Supported by the previous knowledge about regulatory elements and gene function, the methylation data suggests EBF3 and GNAS as potential candidate genes and G protein signaling as potential candidate pathway associated with beef tenderness via methylation.

Comparative transcriptome in large-scale human and cattle populations.

BACKGROUND: Cross-species comparison of transcriptomes is important for elucidating evolutionary molecular mechanisms underpinning phenotypic variation between and within species, yet to date it has been essentially limited to model organisms with relatively small sample sizes. RESULTS: Here, we systematically analyze and compare 10,830 and 4866 publicly available RNA-seq samples in humans and cattle, respectively, representing 20 common tissues. Focusing on 17,315 orthologous genes, we demonstrate that mean/median gene expression, inter-individual variation of expression, expression quantitative trait loci, and gene co-expression networks are generally conserved between humans and cattle. By examining large-scale genome-wide association studies for 46 human traits (average n = 327,973) and 45 cattle traits (average n = 24,635), we reveal that the heritability of complex traits in both species is significantly more enriched in transcriptionally conserved than diverged genes across tissues. CONCLUSIONS: In summary, our study provides a comprehensive comparison of transcriptomes between humans and cattle, which might help decipher the genetic and evolutionary basis of complex traits in both species.

[Isolation of new alleles of the swine TLR4 gene and analysis of its genetic variation].

Using the PCR-SSCP method, the genetic variation in exon 1 of the TLR4 gene was detected among 893 animals, including Asian wild boars, 3 imported commercial and 10 Chinese indigenous swine breeds. This was conducted to analyze the polymorphisms of exon 1 of TLR4 gene in native and foreign pig breeds and aimed at providing a theoretical foundation for further research on the role that TLR4 gene played in immune and defense system. New alleles were isolated for exon 1 of the swine TLR4 gene for the first time, There were 6 genotypes and 3 alleles, in which the Duroc appeared AA, BB, CC, AB, AC and BC genotypes; Sutai pig, which has Duroc pig origin, were detected to be BB, CC, and BC genotypes; Yorkshire and Landrace were detected to be CC and BC genotypes. Wild boar and all 10 Chinese native pig breeds appeared highly conserved in exon 1 of TLR4 gene, with only CC genotype. Among the 3 homozygous genotypes, the CC genotype matches the sequence in GenBank, while a G93C synonymous mutation and a G194A nonsense mutation were found in the BB and AA genotypes, respectively. The correlation between these two mutation points of TLR4 gene with resistance to stress and disease is worthy of further study.

[cDNA microarray on differently expressed genes in duodenum in porcine sensitive or resistant to Escherichia coli F18].

Based on the paired full-sib individuals selected from the established resource populations of Sutai pig that were characterized as resistant or sensitive to ETEC F18, Agilent double labeled cDNA microarray was used to identify the gene expression profiles in duodenum on purpose of investigating the genes related to Escherichia coli F18 receptor, which may cause edema disease and post-weaning diarrhea in piglets, as well as exploring the molecular mechanism about the differences involved in two different lineages. The results showed that thirteen differently expressed genes were found in one matched group including sensitive ones with GG genotype comparing with resistant ones with AA genotype at a two-fold filter, where there were 6 up-regulated genes and 7 down-regulated genes. In the other matched group composed of sensitive ones with AG genotype, 4 up-regulated genes and 2 down-regulated genes, 6 in total were screened out. GO analy-sis revealed that the differently expressed genes participated in many biological processes, such as immune response, ex-tracellular region, bacterial binding, response to external stimulus and so on. Meanwhile, these genes were mainly related to the Glycan Biosynthesis and Metabolism and Immune System pathways. Actually, the roles that they may play in edema disease and post-weaning diarrhea need further study and verification.

Genetic diversity and phylogenetic relationship of nine sheep populations based on microsatellite markers.

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.