Identification and analysis of differentially expressed long non-coding RNAs of Chinese Holstein cattle responses to heat stress

Long non-coding RNA (LncRNA) and epigenetic factors: their role in regulating the adipocytes in bovine

Investigating the involvement of long non-coding RNAs (lncRNAs) and epigenetic processes in bovine adipocytes can provide valuable new insights into controlling adipogenesis in livestock. Long non-coding RNAs have been associated with forming chromatin loops that facilitate enhancer-promoter interactions during adipogenesis, as well as regulating important adipogenic transcription factors like C/EBPα and PPARγ. They significantly influence gene expression regulation at the post-transcriptional level and are extensively researched for their diverse roles in cellular functions. Epigenetic modifications such as chromatin reorganization, histone alterations, and DNA methylation subsequently affect the activation of genes related to adipogenesis and the progression of adipocyte differentiation. By investigating how fat deposition is epigenetically regulated in beef cattle, scientists aim to unravel molecular mechanisms, identify key regulatory genes and pathways, and develop targeted strategies for modifying fat deposition to enhance desirable traits such as marbling and meat tenderness. This review paper delves into lncRNAs and epigenetic factors and their role in regulating bovine adipocytes while focusing on their potential as targets for genetic improvement to increase production efficiency. Recent genomics advancements, including molecular markers and genetic variations, can boost animal productivity, meeting global demands for high-quality meat products. This review establishes a foundation for future research on understanding regulatory networks linked to lncRNAs and epigenetic changes, contributing to both scholarly knowledge advancement and practical applications within animal agriculture.

Long non-coding RNA variability in porcine skeletal muscle

Recently, numerous studies including various tissues have been carried out on long non-coding RNAs (lncRNAs), but still, its variability has not yet been fully understood. In this study, we characterised the inter-individual variability of lncRNAs in pigs, in the context of number, length and expression. Transcriptomes collected from muscle tissue belonging to six Polish Landrace boars (PL1-PL6), including half-brothers (PL1-PL3), were investigated using bioinformatics (lncRNA identification and functional analysis) and statistical (lncRNA variability) methods. The number of lncRNA ranged from 1289 to 3500 per animal, and the total number of common lncRNAs among all boars was 232. The number, length and expression of lncRNAs significantly varied between individuals, and no consistent pattern has been found between pairs of half-brothers. In detail, PL5 exhibits lower expression than the others, while PL4 has significantly higher expression than PL2-PL3 and PL5-PL6. Noteworthy, comparing the inter-individual variability of lncRNA and mRNA expression, they exhibited concordant patterns. The enrichment analysis for common lncRNA target genes determined a variety of biological processes that play fundamental roles in cell biology, and they were mostly related to whole-body homeostasis maintenance, energy and protein synthesis as well as dynamics of multiple nucleoprotein complexes. The high variability of lncRNA landscape in the porcine genome has been revealed in this study. The inter-individual differences have been found in the context of three aspects: the number, length and expression of lncRNAs, which contribute to a better understanding of its complex nature.

Effect of Curcumin on Hepatic mRNA and lncRNA Co-Expression in Heat-Stressed Laying Hens

Heat stress is an important factor affecting poultry production; birds have a range of inflammatory reactions under high-temperature environments. Curcumin has anti-inflammatory and antioxidant effects. The purpose of this experiment was to investigate the effect of dietary curcumin supplementation on the liver transcriptome of laying hens under heat stress conditions. In the animal experiment, a total of 240 Hy-Line brown hens aged 280 days were divided randomly into four different experimental diets with four replicates, and each replicate consisted of 15 hens during a 42-D experiment. The ambient temperature was adjusted to 34 ± 2 °C for 8 h per day, transiting to a range of 22 °C to 28 °C for the remaining 16 h. In the previous study of our lab, it was found that supplemental 150 mg/kg curcumin can improve production performance, antioxidant enzyme activity, and immune function in laying hens under heat stress. To further investigate the regulatory mechanism of curcumin on heat stress-related genes, in total, six samples of three liver tissues from each of 0 mg/kg and 150 mg/kg curcumin test groups were collected for RNA-seq analysis. In the transcriptome analysis, we reported for the first time that the genes related to heat stress of mRNA, such as HSPA8, HSPH1, HSPA2, and DNAJA4, were co-expressed with lncRNA such as XLOC010450, XLOC037987, XLOC053511, XLOC061207, and XLOC100318, and all of these genes are shown to be down-regulated. These findings provide a scientific basis for the possible benefits of dietary curcumin addition in heat-stressed laying hens.

Heat Stress: A Serious Disruptor of the Reproductive Physiology of Dairy Cows

Global warming is a significant threat to the sustainability and profitability of the dairy sector, not only in tropical or subtropical regions but also in temperate zones where extreme summer temperatures have become a new and challenging reality. Prolonged exposure of dairy cows to high temperatures compromises animal welfare, increases morbidity, and suppresses fertility, resulting in devastating economic losses for farmers. To counteract the deleterious effects of heat stress, cattl e employ various adaptive thermoregulatory mechanisms including molecular, endocrine, physiological, and behavioral responses. These adaptations involve the immediate secretion of heat shock proteins and cortisol, followed by a complex network of disrupted secretion of metabolic and reproductive hormones such as prolactin, ghrelin, ovarian steroid, and pituitary gonadotrophins. While the strategic heat stress mitigation measures can restore milk production through modifications of the microclimate and nutritional interventions, the summer fertility records remain at low levels compared to those of the thermoneutral periods of the year. This is because sustainment of high fertility is a multifaceted process that requires appropriate energy balance, undisrupted mode of various hormones secretion to sustain the maturation and fertilizing competence of the oocyte, the normal development of the early embryo and unhampered maternal-embryo crosstalk. In this review, we summarize the major molecular and endocrine responses to elevated temperatures in dairy cows, as well as the impacts on maturing oocytes and early embryos, and discuss the consequences that heat stress brings about in dairy cattle fertility.

Preliminary Transcriptome Analysis of Long Noncoding RNA in Hypothalamic-Pituitary-Mammary Gland Axis of Dairy Cows under Heat Stress

Heat stress (HS) is directly correlated to mammary gland dysfunction in dairy cows, especially in summer. The hypothalamic-pituitary-mammary gland axis (HPM axis) plays an important role in the regulation of stress response and lactation physiology in heat-stressed dairy cows. The aim of this study was to explore the lncRNA profile, and the competitive endogenous RNA (ceRNA) regulatory network in hypothalamus, pituitary, and mammary gland tissues of heat-stressed and normal dairy cows. We performed RNA sequencing (RNA-seq) to identify differentially expressed (DE) lncRNAs, and the ceRNA regulatory network was established in HPM-axis-related tissues. Our results showed that a total of 13, 702 and 202 DE lncRNAs were identified in hypothalamus, pituitary, and mammary glands, respectively. Of lncRNAs, 8, 209 and 45 were up-regulated, and 5, 493 and 157 lncRNAs were down-regulated. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that DE lncRNAs target genes that might play a role in hormone synthesis, secretion and action, apoptosis, mitogen-activated protein kinase (MAPK), AMP-activated protein kinase (AMPK), and mechanistic target of rapamycin (mTOR) signaling pathway. Moreover, the ceRNA regulatory network associated with the MAPK signaling pathway in HPM-axis-related tissues contains 3286 lncRNA-mRNA pairs. Furthermore, the ceRNA regulatory network associated with apoptosis, prolactin, AMPK, and mTOR signaling pathway in the mammary gland contains 772 lncRNA-mRNA pairs. Thus, some lncRNAs may be involved in the regulation of stress response and the physiological process of lactation. The changes in lncRNA expression profiles and ceRNAs (lncRNA-miRNA-mRNA) in HPM-axis-related tissues are the key to affect the stress response and lactation physiology of dairy cows under HS, which provide a theoretical basis for the molecular mechanism in the stress response of HPM-axis-related tissues in dairy cows under HS.

The Role of Heat Shock Proteins in the Pathogenesis of Polycystic Ovarian Syndrome: A Review of the Literature

Polycystic ovarian syndrome (PCOS) is the most common endocrine disorder in women of reproductive age and post-menopausal women. PCOS is a multifactorial heterogeneous disorder associated with a variety of etiologies, outcomes, and clinical manifestations. However, the pathophysiology of PCOS is still unclear. Heat shock proteins (HSPs) have recently been investigated for their role in the pathogenesis of PCOS. HSPs are a class of proteins that act as molecular chaperones and maintain cellular proteostasis. More recently, their actions beyond that of molecular chaperones have highlighted their pathogenic role in several diseases. In PCOS, different HSP family members show abnormal expression that affects the proliferation and apoptotic rates of ovarian cells as well as immunological processes. HSP dysregulation in the ovaries of PCOS subjects leads to a proliferation/apoptosis imbalance that mechanistically impacts follicle stage development, resulting in polycystic ovaries. Moreover, HSPs may play a role in the pathogenesis of PCOS-associated conditions. Recent studies on HSP activity during therapeutic interventions for PCOS suggest that modulating HSP activity may lead to novel treatment strategies. In this review, we summarize what is currently known regarding the role of HSPs in the pathogenesis of PCOS and their potential role in the treatment of PCOS, and we outline areas for future research.

The alteration of N6-methyladenosine (m6A) modification at the transcriptome-wide level in response of heat stress in bovine mammary epithelial cells

Background

Heat stress has a substantial negative economic impact on the dairy industry. N6-methyladenosine (m6A) is the most common internal RNA modification in eukaryotes and plays a key role in regulating heat stress response in animals. In dairy cows, however, this modification remains largely unexplored. Therefore, we examined the effects of heat stress on the m6A modification and gene expression in bovine mammary epithelial cells to elucidate the mechanism of heat stress response. In this study, Mammary alveolar cells-large T antigen (MAC-T) cells were incubated at 37 °C (non-heat stress group, NH) and 40 °C (heat stress group, H) for 2 hours, respectively. HSP70, HSF1, BAX and CASP3 were up regulated in H group compared with those in the NH group.

Results

Methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were conducted to identify m6A peaks and to produce gene expression data of MAC-T cells in the H and NH groups. In total, we identified 17,927 m6A peaks within 9355 genes in the H group, and 18,974 peaks within 9660 genes in the NH groups using MeRIP-seq. Compared with the NH group, 3005 significantly differentially enriched m6A peaks were identified, among which 1131 were up-regulated and 1874 were down-regulated. In addition, 1502 significantly differentially expressed genes were identified using RNA-seq, among which 796 were up-regulated and 706 were down-regulated in the H group compared to the NH group. Furthermore, 199 differentially expressed and synchronously differentially methylated genes were identified by conjoint analysis of the MeRIP-seq and RNA-seq data, which were subsequently divided into four groups: 47 hyper-up, 53 hyper-down, 59 hypo-up and 40 hypo-down genes. In addition, GO enrichment and KEGG analyses were used to analyzed the potential functions of the genes in each section.

Conclusion

The comparisons of m6A modification patterns and conjoint analyses of m6A modification and gene expression profiles suggest that m6A modification plays a critical role in the heat stress response by regulating gene expression.

Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases

The pituitary gland directly regulates the reproduction of domestic animals. Research has increasingly focused on the potential regulatory mechanism of non-coding RNA in pituitary development. Little is known about the differential expression pattern of lncRNAs in Hu sheep, a famous sheep breed with high fecundity, and its role in the pituitary gland between the follicular phase and luteal phase. Herein, to identify the transcriptomic differences of the sheep pituitary gland during the estrus cycle, RNA sequencing (RNA-Seq) was performed. The results showed that 3529 lncRNAs and 16,651 mRNAs were identified in the pituitary gland. Among of them, 144 differentially expressed (DE) lncRNA transcripts and 557 DE mRNA transcripts were screened in the follicular and luteal phases. Moreover, GO and KEGG analyses demonstrated that 39 downregulated and 22 upregulated genes interacted with pituitary functions and reproduction. Lastly, the interaction of the candidate lncRNA XR_001039544.4 and its targeted gene LHB were validated in sheep pituitary cells in vitro. LncRNA XR_001039544.4 and LHB showed high expression levels in the luteal phase in Hu sheep. LncRNA XR_001039544.4 is mainly located in the cytoplasm, as determined by FISH analysis, indicating that XR_001039544.4 might act as competing endogenous RNAs for miRNAs to regulate LHB. LncRNA XR_001039544.4 knockdown significantly inhibited LH secretion and cell proliferation. LncRNA XR_001039544.4 may regulate the secretion of LH in the luteal-phase pituitary gland via affecting cell proliferation. Taken together, these findings provided genome-wide lncRNA- and mRNA-expression profiles for the sheep pituitary gland between the follicular and luteal phases, thereby contributing to the elucidation of the molecular mechanisms of pituitary function.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

Climate-Resilient Dairy Cattle Production: Applications of Genomic Tools and Statistical Models

The current changing climate trend poses a threat to the productive efficacy and welfare of livestock across the globe. This review is an attempt to synthesize information pertaining to the applications of various genomic tools and statistical models that are available to identify climate-resilient dairy cows. The different functional and economical traits which govern milk production play a significant role in determining the cost of milk production. Thus, identification of these traits may revolutionize the breeding programs to develop climate-resilient dairy cattle. Moreover, the genotype–environment interaction also influences the performance of dairy cattle especially during a challenging situation. The recent advancement in molecular biology has led to the development of a few biotechnological tools and statistical models like next-generation sequencing (NGS), microarray technology, whole transcriptome analysis, and genome-wide association studies (GWAS) which can be used to quantify the molecular mechanisms which govern the climate resilience capacity of dairy cows. Among these, the most preferred option for researchers around the globe was GWAS as this approach jointly takes into account all the genotype, phenotype, and pedigree information of farm animals. Furthermore, selection signatures can also help to demarcate functionally important regions in the genome which can be used to detect potential loci and candidate genes that have undergone positive selection in complex milk production traits of dairy cattle. These identified biomarkers can be incorporated in the existing breeding policies using genomic selection to develop climate-resilient dairy cattle.

Genome-wide identification and functional prediction of long non-coding RNAs in Sprague-Dawley rats during heat stress

Background

Heat stress (HS) is a major stress event in the life of an animal, with detrimental upshots in production and health. Long-non-coding RNAs (lncRNAs) play an important role in many biological processes by transcriptional regulation. However, no research has been reported on the characterization and functionality of lncRNAs in heat-stressed rats.

Results

We studied expression levels of lncRNAs in rats during HS, using strand-specific RNA sequencing. Six rats, three in each of Control (22 ± 1 °C) and H120 (42 °C for 120 min) experimental groups, were used to screen for lncRNAs in their liver and adrenal glands. Totally, 4498 and 7627 putative lncRNAs were identified in liver and adrenal glands of the Control and H120 groups, respectively. The majority of lncRNAs were relatively shorter and contained fewer exons than protein-coding transcripts. In total, 482 (174 up-regulated and 308 down-regulated) and 271 (126 up-regulated and 145 down-regulated) differentially-expressed lncRNAs (DElncRNAs, P < 0.05) were identified in the liver and adrenal glands of the Control and H120 groups, respectively. Furthermore, 1274, 121, and 73 target differentially-expressed genes (DEGs) in the liver were predicted to interact with DElncRNAs based on trans−/cis- and sequence similarity regulatory modes. Functional annotation analyses indicated that these DEGs were mostly significantly enriched in insulin signalling, myeloid leukaemia, and glucagon signalling pathways. Similarly, 437, 73 and 41 target DEGs in the adrenal glands were mostly significantly enriched in the cell cycle (trans-prediction) and lysosome pathways (cis-prediction). The DElncRNAs interacting with DEGs that encode heat shock proteins (HSPs) may play an important role in HS response, which include Hsf4, Dnaja1, Dnajb4, Hsph1 and Hspb1 in the liver, and Dnajb13 and Hspb8 in the adrenal glands. The strand-specific RNA sequencing findings were also further verified through RT-qPCR.

Conclusions

This study is the first to provide a detailed characterization and functional analysis of expression levels of lncRNAs in liver and adrenal glands of heat-stressed rats, which provides basis for further studies on the biological functions of lncRNAs under heat stress in rats and other mammalian species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07421-8.

Genome-wide identification and functional prediction of circular RNAs in response to heat stress in Chinese Holstein cows

Heat stress (HS) leads to substantial economic loss of dairy industry each year. The negative effect of HS in dairy cows is becoming one of the more urgent issue due to accelerating side-effects of global warming. Various genes are involved in HS response but the information about the role of noncoding RNAs, especially circular RNAs (circRNAs) is largely unknown. In our study, we aimed to investigate the different expression profile of circRNAs between HS and Non-heat-stressed condition (NC) of Chinese Holstein cow's mammary gland. CircRNAs were identified using RNA sequencing and bioinformatics analysis. In total, 37405 circRNAs were detected and 95 were differentially expressed (DE), including 15 downregulated and 80 upregulated circRNAs in HS group compared to NC. Eight circRNAs were randomly selected to verify the RNA sequencing result. Further, Sanger sequencing validated the backsplicing site of the eight circRNAs. Moreover, results obtained from the Quantitative real time PCR (qRT-PCR) showed consistent expression trend with that of RNA sequencing. GO annotation and KEGG analysis suggested that these DE circRNAs probably involved in the energy metabolic regulation. Furthermore, we constructed ceRNA network and the result indicated that these DE circRNAs could regulate lactation through IGF1 and PRL signaling pathway.

Comprehensive MicroRNA Expression Profile of the Mammary Gland in Lactating Dairy Cows With Extremely Different Milk Protein and Fat Percentages

A total of 31 differentially expressed genes in the mammary glands were identified in our previous study using RNA sequencing (RNA-Seq), for lactating cows with extremely high and low milk protein and fat percentages. To determine the regulation of milk composition traits, we herein investigated the expression profiles of microRNA (miRNA) using small RNA sequencing based on the same samples as in the previous RNA-Seq experiment. A total of 497 known miRNAs (miRBase, release 22.1) and 49 novel miRNAs among the reads were identified. Among these miRNAs, 71 were found differentially expressed between the high and low groups (p < 0.05, q < 0.05). Furthermore, 21 of the differentially expressed genes reported in our previous RNA-Seq study were predicted as target genes for some of the 71 miRNAs. Gene ontology and KEGG pathway analyses showed that these targets were enriched for functions such as metabolism of protein and fat, and development of mammary gland, which indicating the critical role of these miRNAs in regulating the formation of milk protein and fat. With dual luciferase report assay, we further validated the regulatory role of 7 differentially expressed miRNAs through interaction with the specific sequences in 3'UTR of the targets. In conclusion, the current study investigated the complexity of the mammary gland transcriptome in dairy cattle using small RNA-seq. Comprehensive analysis of differential miRNAs expression and the data from previous study RNA-seq provided the opportunity to identify the key candidate genes for milk composition traits.

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

Simple Summary

Long non-coding RNAs (lncRNAs) play a key role in regulating the expression level of mRNAs. The expression profiles of ovine mammary gland were investigated in two sheep breeds with divergent milk phenotype using RNA-Seq. A total of 1894 lncRNAs were found to be expressed and 68 of these were differentially expressed between the two breeds. Some important Gene Ontogeny (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways that were related to lactation and mammary gland morphogenesis were found for the target genes of differentially expressed lncRNAs. This study can improve our understanding of the functions of lncRNAs in the regulation of lactation, milk yield, and milk components in sheep.

Abstract

Long non-coding RNAs (lncRNAs) are a kind of non-coding RNA with >200 nucleotides in length. Some lncRNAs have been proven to have clear regulatory functions in many biological processes of mammals. However, there have been no reports on the roles of lncRNAs in ovine mammary gland tissues. In the study, the expression profiles of lncRNAs were studied using RNA-Seq in mammary gland tissues from lactating Small-Tailed Han (STH) ewes and Gansu Alpine Merino (GAM) ewes with different milk yield and ingredients. A total of 1894 lncRNAs were found to be expressed. Compared with the GAM ewes, the expression levels of 31 lncRNAs were significantly up-regulated in the mammary gland tissues of STH ewes, while 37 lncRNAs were remarkably down-regulated. Gene Ontogeny (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the target genes of differentially expressed lncRNAs were enriched in the development and proliferation of mammary epithelial cells, morphogenesis of mammary gland, ErbB signaling pathway, and Wnt signaling pathway. Some miRNA sponges of differentially expressed lncRNAs, reported to be associated with lactation and mammary gland morphogenesis, were found in a lncRNA-miRNA network. This study reveals comprehensive lncRNAs expression profiles in ovine mammary gland tissues, thereby providing a further understanding of the functions of lncRNAs in the lactation and mammary gland development of sheep.

Exploring the Regulatory Potential of Long Non-Coding RNA in Feed Efficiency of Indicine Cattle

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from the liver, muscle, hypothalamus, pituitary gland and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management.