MicroRNA sequence variation can impact interactions with target mRNA in cattle.

BACKGROUND: Small non-coding microRNAs (miRNAs) are important modulators at post-transcriptional levels. Single-nucleotide polymorphisms (SNPs) located in miRNA genes can alter the secondary structure of pre-miRNA to either impair or promote the miRNA maturation processes. Furthermore, SNPs located in the miRNA seed regions can stabilize or disturb miRNA-target interactions, thereby, quantitatively influence the expression of target genes. Therefore, the main objective of this study was to detect SNPs in bovine miRNAs using the whole-genome re-sequencing datasets of 1632 cattle of five breeds from the 1000 bull genomes project. RESULTS: In total, our study identified 1109, 334, and 130 SNPs in the miRNA precursor, mature, and seed regions, respectively. The heterozygosity values were generally less than 0.3, and the minor allele frequencies (MAFs) were mainly less than 0.1. Most SNPs were in Hardy-Weinberg equilibrium (HWE) (HWE-P > 0.05). Furthermore, we found that the majority of SNPs (MAF > 0.1 and HWE-P > 0.05) in the miRNA seed regions altered the repertoire of target genes, which in turn were enriched in different KEGG pathways or GO terms. Thus target prediction for bta-miR-2888 revealed loss of 309 target genes and gain of 691 target genes. The 691 gained target genes were significantly enriched in 60 KEGG pathways and 21 GO terms. CONCLUSION: In summary, our study identified candidate SNPs in miRNA precursor, mature, and seed regions that are likely to affect target RNA interactions, thereby potentially influencing cattle phenotypic traits.

A frameshift deletion in the GDF9 gene in Icelandic Loa prolific sheep.

Several genetic variants have been shown to affect the mean number of offspring in different sheep breeds. Here, we analyzed samples from Icelandic sheep with the aim of identifying the genetic cause of the Icelandic Loa phenotype using three previously identified prolificacy genes as candidates. We demonstrate that a 4-bp frameshift deletion positioned in the mature region of the GDF9 protein in the Loa animals is a likely causal mutation for the observed increase in prolificacy; however, sequencing showed that not all ewes with a high number of offspring carried the deletion, suggesting the presence of a second mutation segregating within this group of animals.

Integrated Analysis of mRNA and MicroRNA Co-expressed Network for the Differentiation of Bovine Skeletal Muscle Cells After Polyphenol Resveratrol Treatment.

Resveratrol (RSV) has been confirmed to benefit human health. Resveratrol supplemented in the feeds of animals improved pork, chicken, and duck meat qualities. In this study, we identified differentially expressed (DE) messenger RNAs (mRNAs) (n = 3,856) and microRNAs (miRNAs) (n = 93) for the weighted gene co-expression network analysis (WGCNA) to investigate the co-expressed DE mRNAs and DE miRNAs in the primary bovine myoblasts after RSV treatment. The mRNA results indicated that RSV treatments had high correlations with turquoise module (0.91, P-value = 0.01) and blue module (0.93, P-value < 0.01), while only the turquoise module (0.96, P-value < 0.01) was highly correlated with the treatment status using miRNA data. After biological enrichment analysis, the 2,579 DE genes in the turquoise module were significantly enriched in the Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The top two GO terms were actin filament-based process (GO:0030029) and actin cytoskeleton organization (GO:0030036). The top two KEGG pathways were regulation of actin cytoskeleton (bta04810) and tight junction (bta04530). Then, we constructed the DE mRNA co-expression and DE miRNA co-expression networks in the turquoise module and the mRNA-miRNA targeting networks based on their co-expressions in the key module. In summary, the RSV-induced miRNAs participated in the co-expression networks that could affect mRNA expressions to regulate the primary myoblast differentiation. Our study provided a better understanding of the roles of RSV in inducing miRNA and of the characteristics of DE miRNAs in the key co-expressed module in regulation of mRNAs and revealed new candidate regulatory miRNAs and genes for the beef quality traits.

Meta-analysis of genome-wide association studies for cattle stature identifies common genes that regulate body size in mammals.

Stature is affected by many polymorphisms of small effect in humans 1 . In contrast, variation in dogs, even within breeds, has been suggested to be largely due to variants in a small number of genes2,3. Here we use data from cattle to compare the genetic architecture of stature to those in humans and dogs. We conducted a meta-analysis for stature using 58,265 cattle from 17 populations with 25.4 million imputed whole-genome sequence variants. Results showed that the genetic architecture of stature in cattle is similar to that in humans, as the lead variants in 163 significantly associated genomic regions (P < 5 × 10-8) explained at most 13.8% of the phenotypic variance. Most of these variants were noncoding, including variants that were also expression quantitative trait loci (eQTLs) and in ChIP-seq peaks. There was significant overlap in loci for stature with humans and dogs, suggesting that a set of common genes regulates body size in mammals.

Deep sequencing of Danish Holstein dairy cattle for variant detection and insight into potential loss-of-function variants in protein coding genes.

BACKGROUND: Over the last few years, continuous development of high-throughput sequencing platforms and sequence analysis tools has facilitated reliable identification and characterization of genetic variants in many cattle breeds. Deep sequencing of entire genomes within a cattle breed that has not been thoroughly investigated would be imagined to discover functional variants that are underlying phenotypic differences. Here, we sequenced to a high coverage the Danish Holstein cattle breed to detect and characterize single nucleotide polymorphisms (SNPs), insertion/deletions (Indels), and loss-of-function (LoF) variants in protein-coding genes in order to provide a comprehensive resource for subsequent detection of causal variants for recessive traits. RESULTS: We sequenced four genetically unrelated Danish Holstein cows with a mean coverage of 27X using an Illumina Hiseq 2000. Multi-sample SNP calling identified 10,796,794 SNPs and 1,295,036 indels whereof 482,835 (4.5 %) SNPs and 231,359 (17.9 %) indels were novel. A comparison between sequencing-derived SNPs and genotyping from the BovineHD BeadChip revealed a concordance rate of 99.6-99.8 % for homozygous SNPs and 93.3-96.5 % for heterozygous SNPs. Annotation of the SNPs discovered 74,886 SNPs and 1937 indels affecting coding sequences with 2145 being LoF mutations. The frequency of LoF variants differed greatly across the genome, a hot spot with a strikingly high density was observed in a 6 Mb region on BTA18. LoF affected genes were enriched for functional categories related to olfactory reception and underrepresented for genes related to key cellular constituents and cellular and biological process regulation. Filtering using sequence derived genotype data for 288 Holstein animals from the 1000 bull genomes project removing variants containing homozygous individuals retained 345 of the LoF variants as putatively deleterious. A substantial number of the putative deleterious LoF variants had a minor allele frequency >0.05 in the 1000 bull genomes data set. CONCLUSIONS: Deep sequencing of Danish Holstein genomes enabled us to identify 12.1 million variants. An investigation into LoF variants discovered a set of variants predicted to disrupt protein-coding genes. This catalog of variants will be a resource for future studies to understand variation underlying important phenotypes, particularly recessively inherited lethal phenotypes.

Copy number variation in the bovine genome.

BACKGROUND: Copy number variations (CNVs), which represent a significant source of genetic diversity in mammals, have been shown to be associated with phenotypes of clinical relevance and to be causative of disease. Notwithstanding, little is known about the extent to which CNV contributes to genetic variation in cattle. RESULTS: We designed and used a set of NimbleGen CGH arrays that tile across the assayable portion of the cattle genome with approximately 6.3 million probes, at a median probe spacing of 301 bp. This study reports the highest resolution map of copy number variation in the cattle genome, with 304 CNV regions (CNVRs) being identified among the genomes of 20 bovine samples from 4 dairy and beef breeds. The CNVRs identified covered 0.68% (22 Mb) of the genome, and ranged in size from 1.7 to 2,031 kb (median size 16.7 kb). About 20% of the CNVs co-localized with segmental duplications, while 30% encompass genes, of which the majority is involved in environmental response. About 10% of the human orthologous of these genes are associated with human disease susceptibility and, hence, may have important phenotypic consequences. CONCLUSIONS: Together, this analysis provides a useful resource for assessment of the impact of CNVs regarding variation in bovine health and production traits.

Homozygosity for a single base-pair mutation in the oocyte-specific GDF9 gene results in sterility in Thoka sheep.

The control of fecundity is critical in determining mammalian offspring survival. It is regulated principally by the ovulation rate, so that primates and large farm species commonly have a single offspring. Previously, several mutations have been identified in sheep which increase the naturally low ovulation rate; although in some cases homozygous ewes are infertile. In the present study we present a detailed characterization of a novel mutation in growth differentiation factor 9 (GDF9), found in Icelandic Thoka sheep. This mutation is a single base change (A1279C) resulting in a nonconservative amino acid change (S109R) in the C-terminus of the mature GDF9 protein, which is normally expressed in oocytes at all stages of development. Genotyping all animals for which reproductive records were available confirmed this mutation to be associated with increased fecundity in heterozygous ewes and infertility in homozygotes. Analysis of homozygote ovarian morphology and a number of genes normally activated in growing follicles showed that GDF9 was not involved in oocyte activation, but in subsequent development of the follicle. This study highlights the importance of oocyte factors in regulating fertility and provides new information for structural analysis and investigation of the potentially important sites of dimerization or translational modifications required to produce biologically active GDF9. It also provides the basis for the utilization of these animals to enhance sheep production.

Revealing the history of sheep domestication using retrovirus integrations.

The domestication of livestock represented a crucial step in human history. By using endogenous retroviruses as genetic markers, we found that sheep differentiated on the basis of their "retrotype" and morphological traits dispersed across Eurasia and Africa via separate migratory episodes. Relicts of the first migrations include the Mouflon, as well as breeds previously recognized as "primitive" on the basis of their morphology, such as the Orkney, Soay, and the Nordic short-tailed sheep now confined to the periphery of northwest Europe. A later migratory episode, involving sheep with improved production traits, shaped the great majority of present-day breeds. The ability to differentiate genetically primitive sheep from more modern breeds provides valuable insights into the history of sheep domestication.

A snapshot of CNVs in the pig genome.

Recent studies of mammalian genomes have uncovered the extent of copy number variation (CNV) that contributes to phenotypic diversity, including health and disease status. Here we report a first account of CNVs in the pig genome covering part of the chromosomes 4, 7, 14, and 17 already sequenced and assembled. A custom tiling oligonucleotide array was used with a median probe spacing of 409 bp for screening 12 unrelated Duroc boars that are founders of a large family material. After a strict CNV calling pipeline, 37 copy number variable regions (CNVRs) across all four chromosomes were identified, with five CNVRs overlapping segmental duplications, three overlapping pig unigenes and one overlapping a RefSeq pig mRNA. This CNV snapshot analysis is the first of its kind in the porcine genome and constitutes the basis for a better understanding of porcine phenotypes and genotypes with the prospect of identifying important economic traits.