RNA-Seq Reveals the Expression Profiles of Long Non-Coding RNAs in Lactating Mammary Gland from Two Sheep Breeds with Divergent Milk Phenotype

Transcriptome analysis of the Bactrian camel (Camelus bactrianus) reveals candidate genes affecting milk production traits

BACKGROUND

Milk production traits are complex traits with vital economic importance in the camel industry. However, the genetic mechanisms regulating milk production traits in camels remain poorly understood. Therefore, we aimed to identify candidate genes and metabolic pathways that affect milk production traits in Bactrian camels.

METHODS

We classified camels (fourth parity) as low- or high-yield, examined pregnant camels using B-mode ultrasonography, observed the microscopic changes in the mammary gland using hematoxylin and eosin (HE) staining, and used RNA sequencing to identify differentially expressed genes (DEGs) and pathways.

RESULTS

The average standard milk yield over the 300 days during parity was recorded as 470.18 ± 9.75 and 978.34 ± 3.80 kg in low- and high-performance camels, respectively. Nine female Junggar Bactrian camels were subjected to transcriptome sequencing, and 609 and 393 DEGs were identified in the low-yield vs. high-yield (WDL vs. WGH) and pregnancy versus colostrum period (RSQ vs. CRQ) comparison groups, respectively. The DEGs were compared with genes associated with milk production traits in the Animal Quantitative Trait Loci database and in Alashan Bactrian camels, and 65 and 46 overlapping candidate genes were obtained, respectively. Functional enrichment and protein-protein interaction network analyses of the DEGs and candidate genes were conducted. After comparing our results with those of other livestock studies, we identified 16 signaling pathways and 27 core candidate genes associated with maternal parturition, estrogen regulation, initiation of lactation, and milk production traits. The pathways suggest that emerged milk production involves the regulation of multiple complex metabolic and cellular developmental processes in camels. Finally, the RNA sequencing results were validated using quantitative real-time PCR; the 15 selected genes exhibited consistent expression changes.

CONCLUSIONS

This study identified DEGs and metabolic pathways affecting maternal parturition and milk production traits. The results provides a theoretical foundation for further research on the molecular mechanism of genes related to milk production traits in camels. Furthermore, these findings will help improve breeding strategies to achieve the desired milk yield in camels.

Editorial: Role of Non-Coding RNAs in Animals

The importance of non-coding RNAs (ncRNAs), such as microRNAs (miRNA), long non-coding RNAs (lncRNA), and circular RNAs (circRNA), in gene regulation is increasingly being appreciated in many species [...].

RNA-Seq Reveals the Roles of Long Non-Coding RNAs (lncRNAs) in Cashmere Fiber Production Performance of Cashmere Goats in China

Long non-coding RNAs (lncRNAs) are a kind of non-coding RNA being >200 nucleotides in length, and they are found to participate in hair follicle growth and development and wool fiber traits regulation. However, there are limited studies reporting the role of lncRNAs in cashmere fiber production in cashmere goats. In this study, Liaoning cashmere (LC) goats (n = 6) and Ziwuling black (ZB) goats (n = 6) with remarkable divergences in cashmere yield, cashmere fiber diameter, and cashmere color were selected for the construction of expression profiles of lncRNAs in skin tissue using RNA sequencing (RNA-seq). According to our previous report about the expression profiles of mRNAs originated from the same skin tissue as those used in the study, the cis and trans target genes of differentially expressed lncRNAs between the two caprine breeds were screened, resulting in a lncRNA-mRNA network. A total of 129 lncRNAs were differentially expressed in caprine skin tissue samples between LC goats and ZB goats. The presence of 2 cis target genes and 48 trans target genes for the differentially expressed lncRNAs resulted in 2 lncRNA-cis target gene pairs and 93 lncRNA-trans target gene pairs. The target genes concentrated on signaling pathways that were related to fiber follicle development, cashmere fiber diameter, and cashmere fiber color, including PPAR signaling pathway, metabolic pathways, fatty acid metabolism, fatty acid biosynthesis, tyrosine metabolism, and melanogenesis. A lncRNA-mRNA network revealed 22 lncRNA-trans target gene pairs for seven differentially expressed lncRNAs selected, of which 13 trans target genes contributed to regulation of cashmere fiber diameter, while nine trans target genes were responsible for cashmere fiber color. This study brings a clear explanation about the influences of lncRNAs over cashmere fiber traits in cashmere goats.

Potential Novel Biomarkers for Mastitis Diagnosis in Sheep

Simple Summary

Inflammation of the mammary gland (mastitis) is an important disease of dairy sheep. Mastitis management depends mainly on the diagnosis. Conventional diagnostic methods including somatic cell count, California Mastitis Test, and microbial culture have limitations. Therefore researchers are looking for new diagnostic biomarkers of mastitis including specific proteins produced by the liver in case of disease (acute phase proteins), unique genetic sequences (miRNAs), or antimicrobial peptides produced by immune cells during inflammation (cathelicidines).

Abstract

This review aims to characterize promising novel markers of ovine mastitis. Mastitis is considered as one of the primary factors for premature culling in dairy sheep and has noticeable financial, productional, and animal welfare-related implications. Furthermore, clinical, and subclinical mammary infections negatively affect milk yield and alter the milk composition, thereby leading to lowered quality of dairy products. It is, therefore, crucial to control and prevent mastitis through proper diagnosis, treatment or culling, and appropriate udder health management particularly at the end of the lactation period. The clinical form of mastitis is characterized by abnormalities in milk and mammary gland tissue alteration or systemic symptoms consequently causing minor diagnostic difficulties. However, to identify ewes with subclinical mastitis, laboratory diagnostics is crucial. Mastitis control is primarily dependent on determining somatic cell count (SCC) and the California Mastitis Test (CMT), which aim to detect the quantity of cells in the milk sample. The other useful diagnostic tool is microbial culture, which complements SCC and CMT. However, all mentioned diagnostic methods have their limitations and therefore novel biomarkers of ovine subclinical mastitis are highly desired. These sensitive indicators include acute-phase proteins, miRNA, and cathelicidins measurements, which could be determined in ovine serum and/or milk and in the future may become useful in early mastitis diagnostics as well as a preventive tool. This may contribute to increased detection of ovine mammary gland inflammation in sheep, especially in subclinical form, and consequently improves milk quality and quantity.

Comparative Transcriptome Analysis of Milk Somatic Cells During Lactation Between Two Intensively Reared Dairy Sheep Breeds

In dairy sheep industry, milk production dictates the value of a ewe. Milk production is directly related to the morphology and physiology of the mammary gland; both being designated targets of breeding strategies. Although within a flock breeding parameters are mutual, large differences in milk production among individual ewes are usually observed. In this work, we tested two of the most productive dairy sheep breeds reared intensively in Greece, one local the Chios breed and one foreign the Lacaune breed. We used transcriptome sequencing to reveal molecular mechanisms that render the mammary gland highly productive or not. While highly expressed genes (caseins and major whey protein genes) were common among breeds, differences were observed in differentially expressed genes. ENSOARG00000008077, as a member of ribosomal protein 14 family, together with LPCAT2, CCR3, GPSM2, ZNF131, and ASIP were among the genes significantly differentiating mammary gland's productivity in high yielding ewes. Gene ontology terms were mainly linked to the inherent transcriptional activity of the mammary gland (GO:0005524, GO:0030552, GO:0016740, GO:0004842), lipid transfer activity (GO:0005319) and innate immunity (GO:0002376, GO:0075528, GO:0002520). In addition, clusters of genes affecting zinc and iron trafficking into mitochondria were highlighted for high yielding ewes (GO:0071294, GO:0010043). Our analyses provide insights into the molecular pathways involved in lactation between ewes of different performances. Results revealed management issues that should be addressed by breeders in order to move toward increased milk yields through selection of the desired phenotypes. Our results will also contribute toward the selection of the most resilient and productive ewes, thus, will strengthen the existing breeding systems against a spectrum of environmental threats.