Transcriptome Analysis Reveals the Potential Role of Long Non-coding RNAs in Mammary Gland of Yak During Lactation and Dry Period

The landscape of long noncoding RNA expression in the goat brain

The brain regulates multiple metabolic processes, such as food intake, energy expenditure, insulin secretion, hepatic glucose production, and glucose and fatty acid metabolism in adipose tissue, which are fundamental for the maintenance of energy and glucose homeostasis during lactation and pregnancy. In addition, brain expression has a fundamental impact on the development of maternal behavior. Although brain functions are partly regulated by long noncoding RNAs (lncRNAs), their expression profiles have not been characterized in depth in any ruminant species. We have sequenced the transcriptome of 12 brain tissues from 3 goats that were 1 mo pregnant and 4 nonpregnant goats to investigate their lncRNA expression patterns. Between 4,363 (adenohypophysis) and 4,604 (olfactory bulb) lncRNAs were expressed in brain tissues, leading us to establish a set of 794 already annotated lncRNAs and 5,098 novel lncRNA candidates. The detected lncRNAs shared features with those of other mammals, and tissue-specific lncRNAs were enriched in brain development-related terms. Differential expression analyses between goats that were 1 mo pregnant and nonpregnant goats showed that the lncRNA expression profiles of certain brain regions experience substantial changes associated with early pregnancy (238 lncRNAs are differentially expressed in the olfactory bulb), but others do not. Enrichment analysis showed that differentially expressed lncRNAs from the olfactory bulb are co-expressed with genes previously linked to behavioral changes related to pregnancy. These findings provide a first characterization of the landscape of lncRNA expression in the goat brain and provides valuable clues to understand the molecular events triggered by early pregnancy in the central nervous system.

Association between Single Nucleotide Polymorphisms of PRKD1 and KCNQ3 Gene and Milk Quality Traits in Gannan Yak (Bos grunniens)

Protein kinase D1 (PRKD1) functions primarily in normal mammary cells, and the potassium voltage-gated channel subfamily Q member 3 (KCNQ3) gene plays an important role in controlling membrane potential and neuronal excitability, it has been found that this particular gene is linked to the percentage of milk fat in dairy cows. The purpose of this study was to investigate the relationship between nucleotide polymorphisms (SNPs) of PRKD1 and KCNQ3 genes and the milk quality of Gannan yak and to find molecular marker sites that may be used for milk quality breeding of Gannan yak. Three new SNPs were detected in the PRKD1 (g.283,619T>C, g.283,659C>A) and KCNQ3 gene (g.133,741T>C) of 172 Gannan lactating female yaks by Illumina yak cGPS 7K liquid-phase microarray technology. Milk composition was analyzed using a MilkoScanTM milk composition analyzer. We found that the mutations of these three loci significantly improved the lactose, milk fat, casein, protein, non-fat milk solid (SNF) content and acidity of Gannan yaks. The lactose content of the TC heterozygous genotype population at g.283,619T>C locus was significantly higher than that of the TT wild-type population (p < 0.05); the milk fat content of the CA heterozygous genotype population at g.283,659C>A locus was significantly higher than that of the CC wild-type and AA mutant populations (p < 0.05); the casein, protein and acidity of the CC mutant and TC heterozygous groups at the g.133,741T>C locus were significantly higher than those of the wild type (p < 0.05), and the SNF of the TC heterozygous group was significantly higher than that of the mutant group (p < 0.05). The results showed that PRKD1 and KCNQ3 genes could be used as candidate genes affecting the milk traits of Gannan yak.

Polymorphisms of CCSER1 Gene and Their Correlation with Milk Quality Traits in Gannan Yak (Bos grunniens)

Coiled-coil serine-rich protein 1 (CCSER 1) gene is a regulatory protein gene. This gene has been reported to be associated with various economic traits in large mammals in recent years. The aim of this study was to investigate the association between CCSER1 gene single nucleotide polymorphisms (SNPs) and Gannan yaks and to identify potential molecular marker loci for breeding milk quality in Gannan yaks. We genotyped 172 Gannan yaks using Illumina Yak cGPS 7K liquid microarrays and analyzed the correlation between the three SNPs loci of the CCSER1 gene and the milk qualities of Gannan yaks, including milk fat, protein and casein. It was found that mutations at the g.183,843A>G, g.222,717C>G and g.388,723G>T loci all affected the fat, protein, casein and lactose traits of Gannan yak milk to varying extents, and that the milk quality of individuals with mutant phenotypes was significantly improved. Among them, the milk fat content of AG heterozygous genotype population at g.183,843A>G locus was significantly higher than that of AA and GG genotype populations (p < 0.05); the casein and protein content of mutant GG and CG genotype populations at g.222,717C>G locus was significantly higher than that of wild-type CC genotype population (p < 0.05); and the g.388,723G>T locus of the casein and protein contents of the mutant TT genotype population were significantly higher (p < 0.05) than those of the wild-type GG genotype population. These results provide potential molecular marker sites for Gannan yak breeding.

Dynamic Profile of the Yak Mammary Transcriptome during the Lactation Cycle

The objective of this study was to assess the transcriptome of the mammary tissue of four yaks during the whole lactation cycle. For this purpose, biopsies of the mammary gland were performed at -30, -15, 1, 15, 30, 60, 120, 180, and 240 days relative to parturition (d). The transcriptome analysis was performed using a commercial bovine microarray platform and the results were analyzed using several bioinformatic tools. The statistical analysis using an overall false discovery rate ≤ 0.05 for the effect of whole lactation and p < 0.05 for each comparison identified >6000 differentially expressed genes (DEGs) throughout lactation, with a large number of DEGs observed at the onset (1 d vs. -15 d) and at the end of lactation (240 d vs. 180 d). Bioinformatics analysis revealed a major role of genes associated with BTA3, BTA4, BTA6, BTA9, BTA14, and BTA28 in lactation. Functional analysis of DEG underlined an overall induction of lipid metabolism, suggesting an increase in triglycerides synthesis, likely regulated by PPAR signaling. The same analysis revealed an induction of amino acid metabolism and secretion of protein, with a concomitant decrease in proteasome, indicating a major role of amino acid handling and reduced protein degradation in the synthesis and secretion of milk proteins. Glycan biosynthesis was induced for both N-glycan and O-glycan, suggesting increased glycan content in the milk. The cell cycle and immune response, especially antigen processing and presentation, were strongly inhibited during lactation, suggesting that morphological changes are minimized during lactation, while the mammary gland prevents immune hyper-response. Transcripts associated with response to radiation and low oxygen were enriched in the down-regulated DEG affected by the stage of lactation. Except for this last finding, the functions affected by the transcriptomic adaptation to lactation in mammary tissue of yak are very similar to those observed in dairy cows.

Characterization of RNA Editome in the Mammary Gland of Yaks during the Lactation and Dry Periods

Simple Summary

In order to study the influence of RNA editing sites on lactation and mammary gland development process in yaks, we comprehensively characterized the RNA editome of the yak mammary gland during the lactation period and dry period by using the transcriptome and genome sequencing data. The results revealed 82,872 nonredundant RNA editing sites, 14,159 of which were differentially edited between the lactation period and dry period. Enrichment analysis showed that the genes harboring differential editing sites were mainly associated with mammary gland development-related pathways, such as MAPK pathway, PI3K-Akt pathway, FoxO signaling pathway, GnRH signaling pathway, and focal adhesion pathway. Our findings offer some novel insights into the RNA editing function in the mammary gland of yaks.

Abstract

The mammary gland is a complicated organ comprising several types of cells, and it undergoes extensive morphogenetic and metabolic changes during the female reproductive cycle. RNA editing is a posttranscriptional modification event occurring at the RNA nucleotide level, and it drives transcriptomic and proteomic diversities, with potential functional consequences. RNA editing in the mammary gland of yaks, however, remains poorly understood. Here, we used REDItools to identify RNA editing sites in mammary gland tissues in yaks during the lactation period (LP, n = 2) and dry period (DP, n = 3). Totally, 82,872 unique RNA editing sites were identified, most of which were detected in the noncoding regions with a low editing degree. In the coding regions (CDS), we detected 5235 editing sites, among which 1884 caused nonsynonymous amino acid changes. Of these RNA editing sites, 486 were found to generate novel possible miRNA target sites or interfere with the initial miRNA binding sites, indicating that RNA editing was related to gene regulation mediated by miRNA. A total of 14,159 RNA editing sites (involving 3238 common genes) showed a significant differential editing level in the LP when compared with that in the DP through Tukey’s Honest Significant Difference method (p < 0.05). According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, genes that showed different RNA editing levels mainly participated in pathways highly related to mammary gland development, including MAPK, PI3K-Akt, FoxO, and GnRH signaling pathways. Collectively, this work demonstrated for the first time the dynamic RNA editome profiles in the mammary gland of yaks and shed more light on the mechanism that regulates lactation together with mammary gland development.

Transcriptome analysis reveals potential function of long non-coding RNAs in 20-hydroxyecdysone regulated autophagy in Bombyx mori

Background

20-hydroxyecdysone (20E) plays important roles in insect molting and metamorphosis. 20E-induced autophagy has been detected during the larval–pupal transition in different insects. In Bombyx mori, autophagy is induced by 20E in the larval fat body. Long non-coding RNAs (lncRNAs) function in various biological processes in many organisms, including insects. Many lncRNAs have been reported to be potential for autophagy occurrence in mammals, but it has not been investigated in insects.

Results

RNA libraries from the fat body of B. mori dissected at 2 and 6 h post-injection with 20E were constructed and sequenced, and comprehensive analysis of lncRNAs and mRNAs was performed. A total of 1035 lncRNAs were identified, including 905 lincRNAs and 130 antisense lncRNAs. Compared with mRNAs, lncRNAs had longer transcript length and fewer exons. 132 lncRNAs were found differentially expressed at 2 h post injection, compared with 64 lncRNAs at 6 h post injection. Thirty differentially expressed lncRNAs were common at 2 and 6 h post-injection, and were hypothesized to be associated with the 20E response. Target gene analysis predicted 6493 lncRNA-mRNA cis pairs and 42,797 lncRNA-mRNA trans pairs. The expression profiles of LNC_000560 were highly consistent with its potential target genes, Atg4B, and RNAi of LNC_000560 significantly decreased the expression of LNC_000560 and Atg4B. These results indicated that LNC_000560 was potentially involved in the 20E-induced autophagy of the fat body by regulating Atg4B.

Conclusions

This study provides the genome-wide identification and functional characterization of lncRNAs associated with 20E-induced autophagy in the fat body of B. mori. LNC_000560 and its potential target gene were identified to be related to 20-regulated autophagy in B. mori. These results will be helpful for further studying the regulatory mechanisms of lncRNAs in autophagy and other biological processes in this insect model.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07692-1.