Integrated Small RNA Sequencing, Transcriptome and GWAS Data Reveal microRNA Regulation in Response to Milk Protein Traits in Chinese Holstein Cattle

Abstract

Milk protein is one of the most important economic traits in the dairy industry. Yet, the regulatory network of miRNAs for the synthesis of milk protein in mammary is poorly understood. Samples from 12 Chinese Holstein cows with three high ( ≥ 3.5%) and three low ( ≤ 3.0%) phenotypic values for milk protein percentage in lactation and non-lactation were examined through deep small RNA sequencing. We characterized 388 known and 212 novel miRNAs in the mammary gland. Differentially expressed analysis detected 28 miRNAs in lactation and 52 miRNAs in the non-lactating period with a highly significant correlation with milk protein concentration. Target prediction and correlation analysis identified some key miRNAs and their targets potentially involved in the synthesis of milk protein. We analyzed for enrichments of GWAS signals in miRNAs and their correlated targets. Our results demonstrated that genomic regions harboring DE miRNA genes in lactation were significantly enriched with GWAS signals for milk protein percentage traits and that enrichments within DE miRNA targets were significantly higher than in random gene sets for the majority of milk production traits. This integrated study on the transcriptome and posttranscriptional regulatory profiles between significantly differential phenotypes of milk protein concentration provides new insights into the mechanism of milk protein synthesis, which should reveal the regulatory mechanisms of milk secretion.

Transcriptome Analysis of Bovine Rumen Tissue in Three Developmental Stages

“…Chromosome locations of PAQR3 co-located with QTL associated with fat deposition and one missense variant likely affect intramuscular fat via the LNPEP gene ( Pena et al, 2013 ; Derks et al, 2021 ). We integrated differential genes obtained from the transcriptome with GWAS data using sum-based marker-set test method, which have been shown to be more potent or at least equivalent to most commonly used marker-set test methods for polygenic traits ( Sørensen et al, 2017 ; Cai et al, 2021 ). In birth vs. youth comparison, the GWAS signals of MS trait were significantly enriched in fold change 0∼1 upregulated genes, which implied that these slightly upregulated genes were related to beef quality.…”

“…To investigate whether DEGs and gene co-expression modules were enriched with GWAS signals of carcass and other traits, we applied a sum-based method for GWAS signals enrichment analyses (sumGSE, https://github. com/WentaoCai/GWAS_enrichment) (Cai et al, 2021), which used signals of all markers within a pre-defined list of DEGs and then calculated the following summary statistics for the DEGs:…”

Transcriptional states and chromatin accessibility during bovine myoblasts proliferation and myogenic differentiation

“…To check whether the SNP effects were more enriched in chromatin peaks than background regions, GWAS signals enrichment analysis were performed using Perl scripts (sumGSE, https://github.com/WentaoCai/GWAS_enrichment). 51 …”

Integrated analysis of transcriptome and metabolome revealed biological basis of sows from estrus to lactation

“…For small RNAs, by Cutadapt 2.8, we trimmed the quality and removed the adaptor of the reads, and identified the miRNAs according to the description method in the report research. 43 The reference genome was Scrofa11.1, and the novel miRNAs were detected by combining unannotated sequences with the known miRNA annotation from Ovis aries , Bos taurus , and Gallus gallus .…”