Fat Mass- and Obesity-Associated Protein (FTO) Promotes the Proliferation of Goat Skeletal Muscle Satellite Cells by Stabilizing DAG1 mRNA in an IGF2BP1-Related m6A Manner.

Skeletal muscle development is spotlighted in mammals since it closely relates to animal health and economic benefits to the breeding industry. Researchers have successfully unveiled many regulatory factors and mechanisms involving myogenesis. However, the effect of N6-methyladenosine (m6A) modification, especially demethylase and its regulated genes, on muscle development remains to be further explored. Here, we found that the typical demethylase FTO (fat mass- and obesity-associated protein) was highly enriched in goats' longissimus dorsi (LD) muscles. In addition, the level of m6A modification on transcripts was negatively regulated by FTO during the proliferation of goat skeletal muscle satellite cells (MuSCs). Moreover, a deficiency of FTO in MuSCs significantly retarded their proliferation and promoted the expression of dystrophin-associated protein 1 (DAG1). m6A modifications of DAG1 mRNA were efficiently altered by FTO. Intriguingly, the results of DAG1 levels and its m6A enrichment from FB23-2 (FTO demethylase inhibitor)-treated cells were consistent with those of the FTO knockdown, indicating that the regulation of FTO on DAG1 depended on m6A modification. Further experiments showed that interfering FTO improved m6A modification at site DAG1-122, recognized by Insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) and consequently stabilized DAG1 transcripts. Our study suggests that FTO promotes the proliferation of MuSCs by regulating the expression of DAG1 through m6A modification. This will extend our knowledge of the m6A-related mechanism of skeletal muscle development in animals.

Transcriptomic and Metabolomic Analyses Reveal Molecular Regulatory Networks for Pigmentation Deposition in Sheep.

Domestic animals have multiple phenotypes of skin and coat color, which arise from different genes and their products, such as proteins and metabolites responsible with melanin deposition. However, the complex regulatory network of melanin synthesis remains to be fully unraveled. Here, the skin and tongue tissues of Liangshan black sheep (black group) and Liangshan semi-fine-wool sheep (pink group) were collected, stained with hematoxylin-eosin (HE) and Masson-Fontana, and the transcriptomic and metabolomic data were further analyzed. We found a large deposit of melanin granules in the epidermis of the black skin and tongue. Transcriptome and metabolome analysis identified 744 differentially expressed genes (DEGs) and 443 differentially expressed metabolites (DEMs) between the pink and black groups. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses revealed the DEGs and DEMs were mainly enriched in the pathways of secondary metabolic processes, melanin biosynthesis processes, melanin metabolism processes, melanosome membranes, pigment granule membranes, melanosome, tyrosine metabolism, and melanogenesis. Notably, we revealed the gene ENSARG00020006042 may be a family member of YWHAs and involved in regulating melanin deposition. Furthermore, several essential genes (TYR, TYRP1, DCT, PMEL, MLANA, SLC45A2) were significantly associated with metabolite prostaglandins and compounds involved in sheep pigmentation. These findings provide new evidence of the strong correlation between prostaglandins and related compounds and key genes that regulate sheep melanin synthesis, furthering our understanding of the regulatory mechanisms and molecular breeding of pigmentation in sheep.

Transcriptome analysis reveals long non-coding natural antisense transcripts involved in muscle development in fetal goat (Capra hircus).

Non-coding RNAs have been shown to play vital roles in muscle development. However, the biological roles of long non-coding natural antisense transcripts, antisense lncRNAs (ASlncRNAs), are largely unknown in embryonic muscle development. Here, we identified a total of 466 ASlncRNAs in the longissimus dorsi muscle. And 48 differentially expressed ASlncRNAs were identified based on pairwise comparisons (P < 0.05), sixteen of which were validated by qPCR. Additionally, 466 ASlncRNAs were predicted to target 335 protein coding genes based on complementary base-pairing. Enrichment analysis suggests that ASlncRNAs may be involved in muscle development by negatively regulating the expression of target genes. Furthermore, 170 ASlncRNAs were identified as potential miRNA precursors, suggesting that these ASlncRNAs may be involved in the regulation of muscle development by producing miRNA precursors. Our results provide a catalog of goat muscle-related ASlncRNAs, and will contribute to a fuller understanding of the roles of ASlncRNAs in muscle development.

HuR Promotes the Differentiation of Goat Skeletal Muscle Satellite Cells by Regulating Myomaker mRNA Stability.

Human antigen R (HuR) is an RNA-binding protein that contributes to a wide variety of biological processes and diseases. HuR has been demonstrated to regulate muscle growth and development, but its regulatory mechanisms are not well understood, especially in goats. In this study, we found that HuR was highly expressed in the skeletal muscle of goats, and its expression levels changed during longissimus dorsi muscle development in goats. The effects of HuR on goat skeletal muscle development were explored using skeletal muscle satellite cells (MuSCs) as a model. The overexpression of HuR accelerated the expression of myogenic differentiation 1 (MyoD), Myogenin (MyoG), myosin heavy chain (MyHC), and the formation of myotubes, while the knockdown of HuR showed opposite effects in MuSCs. In addition, the inhibition of HuR expression significantly reduced the mRNA stability of MyoD and MyoG. To determine the downstream genes affected by HuR at the differentiation stage, we conducted RNA-Seq using MuSCs treated with small interfering RNA, targeting HuR. The RNA-Seq screened 31 upregulated and 113 downregulated differentially expressed genes (DEGs) in which 11 DEGs related to muscle differentiation were screened for quantitative real-time PCR (qRT-PCR) detection. Compared to the control group, the expression of three DEGs (Myomaker, CHRNA1, and CAPN6) was significantly reduced in the siRNA-HuR group (p < 0.01). In this mechanism, HuR bound to Myomaker and increased the mRNA stability of Myomaker. It then positively regulated the expression of Myomaker. Moreover, the rescue experiments indicated that the overexpression of HuR may reverse the inhibitory impact of Myomaker on myoblast differentiation. Together, our findings reveal a novel role for HuR in promoting muscle differentiation in goats by increasing the stability of Myomaker mRNA.

Transcriptome Analysis Reveals the Profile of Long Non-Coding RNAs during Myogenic Differentiation in Goats.

The long non-coding RNAs (lncRNAs) are emerging as essential regulators of the growth and development of skeletal muscles. However, little is known about the expression profiles of lncRNAs during the proliferation and differentiation of skeletal muscle satellite cells (MuSCs) in goats. In this study, we investigate potential regulatory lncRNAs that govern muscle development by performing lncRNA expression profiling analysis during the proliferation (cultured in the growth medium, GM) and differentiation (cultured in the differentiation medium, DM1/DM5) of MuSCs. In total, 1001 lncRNAs were identified in MuSC samples, and 314 differentially expressed (DE) (FDR < 0.05, |log2FC| > 1) lncRNAs were screened by pairwise comparisons from three comparison groups (GM-vs-DM1, GM-vs-DM5, DM1-vs-DM5). Moreover, we identified the cis-, trans-, and antisense-regulatory target genes of DE lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that these target genes were significantly enriched in muscle development-related GO terms and KEGG pathways. In addition, the network of interactions between DE lncRNAs and their target genes was identified, which included well-known myogenesis regulators such as Myogenic differentiation 1 (MyoD), Myogenin (MyoG), and Myosin heavy chain (MyHC). Meanwhile, competing endogenous RNA (ceRNA) network analysis showed that 237 DE lncRNAs could bind to 329 microRNAs (miRNAs), while miRNAs could target 564 mRNAs. Together, our results provide a genome-wide resource of lncRNAs that may contribute to myogenic differentiation in goats and lay the groundwork for future investigation into their functions during skeletal muscle development.

Combining 16S rRNA Sequencing and Metabolomics Data to Decipher the Interactions between Gut Microbiota, Host Immunity, and Metabolites in Diarrheic Young Small Ruminants.

Diarrhea is associated with gut microbiota, immunity, and metabolic alterations in goat kids and lambs. This study used 28 lambs (11 healthy and 17 diarrheic) and 20 goat kids (10 healthy and 10 diarrheic) to investigate the association between diarrhea occurrence and changes in gut microbiota, metabolism, and immunity in goat kids and lambs. The results revealed that Firmicutes, Proteobacteria, and Bacteroidetes were the dominant phyla in goat kids and lambs. In addition, Enterobacteriaceae and Lachnospiraceae families were identified in both diarrheic goat kids and lambs. Furthermore, functional prediction of microbiota showed that it was involved in cell motility and cancer pathways. The identified differential metabolites were implicated in the bile secretion pathway. Lambs had significant differences in immunoglobulin G (IgG), immunoglobulin M (IgM), interleukin-1ß (IL-1ß), and tumor necrosis factor-alpha (TNF-α) compared to goat kids. IgG and IL-1ß were positively correlated to Patescibacteria, Clostridiaceae, and unclassified_Muribaculaceae in both diarrheic goat kids and lambs. In addition, weighted gene co-expression network analysis (WGCNA) revealed that the MEgreen module was positively associated with IgG, IgM, IL-1ß, TNF-α, and triglyceride (TG). In conclusion, our results characterized the gut microbiota, metabolism, and immune status of lambs and goat kids suffering from diarrhea.

METTL3 Promotes the Differentiation of Goat Skeletal Muscle Satellite Cells by Regulating MEF2C mRNA Stability in a m6A-Dependent Manner.

The development of mammalian skeletal muscle is a highly complex process involving multiple molecular interactions. As a prevalent RNA modification, N6-methyladenosine (m6A) regulates the expression of target genes to affect mammalian development. Nevertheless, it remains unclear how m6A participates in the development of goat muscle. In this study, methyltransferase 3 (METTL3) was significantly enriched in goat longissimus dorsi (LD) tissue. In addition, the global m6A modification level and differentiation of skeletal muscle satellite cells (MuSCs) were regulated by METTL3. By performing mRNA-seq analysis, 8050 candidate genes exhibited significant changes in expression level after the knockdown of METTL3 in MuSCs. Additionally, methylated RNA immunoprecipitation sequencing (MeRIP-seq) illustrated that myocyte enhancer factor 2c (MEF2C) mRNA contained m6A modification. Further experiments demonstrated that METTL3 enhanced the differentiation of MuSCs by upregulating m6A levels and expression of MEF2C. Moreover, the m6A reader YTH N6-methyladenosine RNA binding protein C1 (YTHDC1) was bound and stabilized to MEF2C mRNA. The present study reveals that METTL3 enhances myogenic differentiation in MuSCs by regulating MEF2C and provides evidence of a post-transcriptional mechanism in the development of goat skeletal muscle.

miR-145-3p Inhibits MuSCs Proliferation and Mitochondria Mass via Targeting MYBL1 in Jianzhou Big-Eared Goats.

Muscle growth and injury-induced regeneration are controlled by skeletal muscle satellite cells (MuSCs) through myogenesis in postnatal animals. Meanwhile, myogenesis is accompanied by mitochondrial function and enzyme activity. Nevertheless, the underlying molecular mechanisms involving non-coding RNAs including circular RNAs (circRNAs) and microRNAs (miRNAs) remain largely unsolved. Here, we explored the myogenic roles of miR-145-3p and MYBL1 on muscle development and mitochondrial mass. We noticed that overexpression of miR-145-3p inhibited MuSCs proliferation and reduced the number of viable cells. Meanwhile, deficiency of miR-145-3p caused by LNAantimiR-145-3p or an inhibitor retarded the differentiation of MuSCs. miR-145-3p altered the mitochondrial mass in MuSCs. Moreover, miR-145-3p targeted and negatively regulated the expression of CDR1as and MYBL1. The knockdown of the MYBL1 using ASO-2'MOE modification simulated the inhibitory function of miR-145-3p on cell proliferation. Additionally, MYBL1 mediated the regulation of miR-145-3p on Vexin, VCPIP1, COX1, COX2, and Pax7. These imply that CDR1as/miR-145-3p/MYBL1/COX1, COX2, VCPIP1/Vexin expression at least partly results in a reduction in mitochondrial mass and MuSCs proliferation. These novel findings confirm the importance of mitochondrial mass during myogenesis and the boosting of muscle/meat development in mammals.

miR-193b-3p Promotes Proliferation of Goat Skeletal Muscle Satellite Cells through Activating IGF2BP1.

As a well-known cancer-related miRNA, miR-193b-3p is enriched in skeletal muscle and dysregulated in muscle disease. However, the mechanism underpinning this has not been addressed so far. Here, we probed the impact of miR-193b-3p on myogenesis by mainly using goat tissues and skeletal muscle satellite cells (MuSCs), compared with mouse C2C12 myoblasts. miR-193b-3p is highly expressed in goat skeletal muscles, and ectopic miR-193b-3p promotes MuSCs proliferation and differentiation. Moreover, insulin-like growth factor-2 mRNA-binding protein 1 (IGF2BP1) is the most activated insulin signaling gene when there is overexpression of miR-193b-3p; the miRNA recognition element (MRE) within the IGF1BP1 3' untranslated region (UTR) is indispensable for its activation. Consistently, expression patterns and functions of IGF2BP1 were similar to those of miR-193b-3p in tissues and MuSCs. In comparison, ectopic miR-193b-3p failed to induce PAX7 expression and myoblast proliferation when there was IGF2BP1 knockdown. Furthermore, miR-193b-3p destabilized IGF2BP1 mRNA, but unexpectedly promoted levels of IGF2BP1 heteronuclear RNA (hnRNA), dramatically. Moreover, miR-193b-3p could induce its neighboring genes. However, miR-193b-3p inversely regulated IGF2BP1 and myoblast proliferation in the mouse C2C12 myoblast. These data unveil that goat miR-193b-3p promotes myoblast proliferation via activating IGF2BP1 by binding to its 3' UTR. Our novel findings highlight the positive regulation between miRNA and its target genes in muscle development, which further extends the repertoire of miRNA functions.

CircTCF4 Suppresses Proliferation and Differentiation of Goat Skeletal Muscle Satellite Cells Independent from AGO2 Binding.

The proliferation and differentiation of mammalian skeletal muscle satellite cells (MuSCs) are highly complicated. Apart from the regulatory signaling cascade driven by the protein-coding genes, non-coding RNAs such as microRNAs (miRNA) and circular RNAs (circRNAs) play essential roles in this biological process. However, circRNA functions in MuSCs proliferation and differentiation remain largely to be elucidated. Here, we screened for an exonic circTCF4 based on our previous RNA-Seq data, specifically expressed during the development of the longest dorsal muscle in goats. Subsequently, the circular structure and whole sequence of circTCF4 were verified using Sanger sequencing. Besides, circTCF4 was spatiotemporally expressed in multiple tissues from goats but strikingly enriched in muscles. Furthermore, circTCF4 suppressed MuSCs proliferation and differentiation, independent of AGO2 binding. Finally, we conducted Poly(A) RNA-Seq using cells treated with small interfering RNA targeting circTCF4 and found that circTCF4 would affect multiple signaling pathways, including the insulin signaling pathway and AMPK signaling pathway related to muscle differentiation. Our results provide additional solid evidence for circRNA regulating skeletal muscle formation.

Enhancer RNAs: transcriptional regulators and workmates of NamiRNAs in myogenesis.

miRNAs are well known to be gene repressors. A newly identified class of miRNAs termed nuclear activating miRNAs (NamiRNAs), transcribed from miRNA loci that exhibit enhancer features, promote gene expression via binding to the promoter and enhancer marker regions of the target genes. Meanwhile, activated enhancers produce endogenous non-coding RNAs (named enhancer RNAs, eRNAs) to activate gene expression. During chromatin looping, transcribed eRNAs interact with NamiRNAs through enhancer-promoter interaction to perform similar functions. Here, we review the functional differences and similarities between eRNAs and NamiRNAs in myogenesis and disease. We also propose models demonstrating their mutual mechanism and function. We conclude that eRNAs are active molecules, transcriptional regulators, and partners of NamiRNAs, rather than mere RNAs produced during enhancer activation.

CircRNA-Protein Interactions in Muscle Development and Diseases.

Circular RNA (circRNA) is a kind of novel endogenous noncoding RNA formed through back-splicing of mRNA precursor. The biogenesis, degradation, nucleus-cytoplasm transport, location, and even translation of circRNA are controlled by RNA-binding proteins (RBPs). Therefore, circRNAs and the chaperoned RBPs play critical roles in biological functions that significantly contribute to normal animal development and disease. In this review, we systematically characterize the possible molecular mechanism of circRNA-protein interactions, summarize the latest research on circRNA-protein interactions in muscle development and myocardial disease, and discuss the future application of circRNA in treating muscle diseases. Finally, we provide several valid prediction methods and experimental verification approaches. Our review reveals the significance of circRNAs and their protein chaperones and provides a reference for further study in this field.

Exploration of miRNA and mRNA Profiles in Fresh and Frozen-Thawed Boar Sperm by Transcriptome and Small RNA Sequencing.

Due to lower farrowing rate and reduced litter size with frozen-thawed semen, over 90% of artificial insemination (AI) is conducted using liquid stored boar semen. Although substantial progress has been made towards optimizing the cryopreservation protocols for boar sperm, the influencing factors and underlying mechanisms related to cryoinjury and freeze tolerance of boar sperm remain largely unknown. In this study, we report the differential expression of mRNAs and miRNAs between fresh and frozen-thawed boar sperm using high-throughput RNA sequencing. Our results showed that 567 mRNAs and 135 miRNAs were differentially expressed (DE) in fresh and frozen-thawed boar sperm. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the majority of DE mRNAs were enriched in environmental information processing such as cytokine-cytokine receptor interactions, PI3K-Akt signaling, cell adhesion, MAPK, and calcium signaling pathways. Moreover, the targets of DE miRNAs were enriched in significant GO terms such as cell process, protein binding, and response to stimuli. In conclusion, we speculate that DE mRNAs and miRNAs are heavily involved in boar sperm response to environment stimuli, apoptosis, and metabolic activities. The differences in expression also reflect the various structural and functional changes in sperm during cryopreservation.

Cryopreservation of boar sperm induces differential microRNAs expression.

Lower conception rates and litter sizes limit the wide use of artificial insemination with frozen-thawed boar sperm, due to a lack of understanding of the mechanisms that cause cryodamage and cryoinjury to sperm during cryopreservation. CryoMiRs, a family of freeze-related microRNAs (miRNAs), are associated with freeze tolerance, and regulate metabolism in mammalian hibernators and insects. Thus, we speculate that miRNAs maybe involved in the regulation of the freeze-thaw process and may affect boar sperm function. In this study, we studied the differential expression of 46 miRNAs that have roles in spermatogenesis, sperm maturation, and sperm quality in response to cryopreservation (with or without 3% glycerol). The results indicated that, in response to cryopreservation with 3% glycerol, 14 miRNAs were significantly up-regulated, but only two miRNAs (miR-22 and miR-450b-5p) were significantly down-regulated, relative to fresh sperm. Preservation with 3% glycerol caused up-regulation of 17 miRNAs, but only caused down-regulation of one miRNA (miR-24), relative to sperm cryopreserved without glycerol. Functional annotations of these differentially expressed miRNAs indicated that these miRNAs and their targets are mainly associated with metabolic and cellular processes. Therefore, our findings show that cryopreservation results in changes in miRNA expression, and suggest that the anti-freeze mechanisms of boar sperm need to be studied further.

ADAR1 Promotes Myogenic Proliferation and Differentiation of Goat Skeletal Muscle Satellite Cells.

As one of the most important economic traits for domestic animal husbandry, skeletal muscle is regulated by an intricate molecular network. Adenosine deaminase acting on RNA (ADAR1) involves various physiological processes and diseases, such as innate immunity and the development of lung adenocarcinoma, breast cancer, gastric cancer, etc. However, its role in skeletal muscle growth requires further clarification. Here, we explored the functions of ADAR1 in the myogenic process of goat skeletal muscle satellite cells (MuSCs). The ADAR1 transcripts were noticeably enriched in goat visceral tissues compared to skeletal muscle. Additionally, its levels in slow oxidative muscles like the psoas major and minor muscles were higher than in the fast oxidative glycolytic and fast glycolytic muscles. Among the two common isoforms from ADAR1, p110 is more abundant than p150. Moreover, overexpressing ADAR1 enhanced the proliferation and myogenic differentiation of MuSCs. The mRNA-seq performed on MuSCs' knockdown of ADAR1 obtained 146 differentially expressed genes (DEGs), 87 upregulated and 59 downregulated. These DEGs were concentrated in muscle development and process pathways, such as the MAPK and cAMP signaling pathways. Furthermore, many DEGs as the key nodes defined by protein-protein interaction networks (PPI), including STAT3, MYH3/8, TGFß2, and ACTN4, were closely related to the myogenic process. Finally, RNA immunoprecipitation combined with qPCR (RIP-qPCR) showed that ADAR1 binds to PAX7 and MyoD mRNA. This study indicates that ADAR1 promotes the myogenic development of goat MuSCs, which provides a useful scientific reference for further exploring the ADAR1-related regulatory networks underlying mammal skeletal muscle growth.

Global A-to-I RNA editing during myogenic differentiation of goat MuSCs.

Background: RNA editing, especially A-to-I editing sites, is a common RNA modification critical for stem cell differentiation, muscle development, and disease occurrence. Unveiling comprehensive RNA A-to-I editing events associated with myogenesis of the skeletal muscle satellite cells (MuSCs) is essential for extending our knowledge of the mechanism underpinning muscle development. Results: A total of 9,632 RNA editing sites (RESs) were screened in the myoblasts (GM), myocytes (DM1), and myotubes (DM5) samples. Among these sites, 4,559 A-to-I edits were classified and further analyzed. There were 3,266 A-to-I sites in the protein-coding region, out of which 113 missense sites recoded protein. Notably, five A-to-I sites in the 3' UTR of four genes (TRAF6, NALF1, SLC38A1, ENSCHIG00000019092) altered their targeted miRNAs. Furthermore, a total of 370 A-to-I sites with different editing levels were detected, including FBN1, MYH10, GSK3B, CSNK1D, and PRKACB genes. These genes were predominantly enriched in the cytoskeleton in muscle cells, the hippo signaling pathway, and the tight junction. Furthermore, we identified 14 hub genes (TUFM, GSK3B, JAK2, RPSA, YARS1, CDH2, PRKACB, RUNX1, NOTCH2, CDC23, VCP, FBN1, RARS1, MEF2C) that potentially related to muscle development. Additionally, 123 stage-specific A-to-I editing sites were identified, with 43 sites in GM, 25 in DM1, and 55 in DM5 samples. These stage-specific edited genes significantly enriched essential biological pathways, including the cell cycle, oocyte meiosis, motor proteins, and hedgehog signaling pathway. Conclusion: We systematically identified the RNA editing events in proliferating and differentiating goat MuSCs, which was crucial for expanding our understanding of the regulatory mechanisms of muscle development.

The novel RNA-RNA activation of H19 on MyoD transcripts promoting myogenic differentiation of goat muscle satellite cells.

The elaborate interplay of coding and noncoding factors governs muscle growth and development. Here, we reported a mutual activation between long noncoding RNA (lncRNA) H19 and MyoD (myogenic determination gene number 1) in the muscle process. We successfully cloned the two isoforms of goat H19, which were significantly enriched and positively correlated with MyoD transcripts in skeletal muscles or differentiating muscle satellite cells (MuSCs). To systematically screen genes altered by H19, we performed RNA-seq using cDNA libraries of differentiating H19-deficiency MuSCs and consequently anchored MyoD as the critical genes in mediating H19 function. Intriguingly, some transcripts of MyoD and H19 overlapped in the cytoplasm, which was dramatically damaged when the core complementary nucleotides were mutated. Meanwhile, MyoD RNA successfully pulled down H19 in MS2-RIP experiments. Furthermore, HuR could bind both H19 and MyoD transcripts, while H19 or its truncated mutants successfully stabilized MyoD mRNA, with or without HuR deficiency. In turn, novel functional MyoD protein-binding sites were identified in the promoter and exons of the H19 gene. Our results suggest that MyoD activates H19 transcriptionally, and RNA-RNA hybridization is critical for H19-promoted MyoD expression, which extends our knowledge of the hierarchy of regulatory networks in muscle growth.

LncDGAT2 is a novel positive regulator of the goat adipocyte thermogenic gene program.

Brown and beige adipose thermogenesis are important for newborn mammals to maintain their body temperature. In addition, these thermogenic fats are regulated by multiple molecular interactions. How the long non-coding RNAs (lncRNAs) regulate adipose thermogenesis in newborn mammals upon cold exposure remains unexplored. Here, we identified lncRNAs induced by cold exposure in brown adipose tissue (BAT) of newborn goats and found that lncDGAT2 was enriched in BAT after cold exposure. Functional studies revealed that lncDGAT2 promoted brown and white adipocyte differentiation as well as thermogenic gene expression. Additionally, PRDM4 directly bound the lncDGAT2 promoter to activate the transcription of lncDGAT2 and the PRDM4-lncDGAT2 axis was essential for the brown adipocyte thermogenic gene program. These findings provide evidence for lncRNA and transcription factor regulatory functions in controlling adipose thermogenesis and energy metabolism of newborn goats.

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array.

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5-10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.

LncRNA-mRNA modules involved in goat rumen development: Insights from genome-wide transcriptome profiling.

The rumen is an essential digestive and absorption organ of ruminants. During fetal life, lactation, and post-weaning period, goat rumen undergoes drastic morphological and metabolic-functional changes triggered by potential regulated genes and non-coding RNA molecules. As the essential regulatory factors, long non-coding RNAs (lncRNAs) have vital functions in various biological activities. However, their roles during rumen development are still poorly explored in ruminants. To explore the genome-wide expression profiles of lncRNAs and mRNAs in the goat rumens, we generated 5,007 lncRNAs and 19,738 mRNAs identified during the fetal and prepubertal stages by the high-throughput RNA sequencing. Notably, 365 lncRNAs and 2,877 mRNAs were considered to be differentially expressed. The weighted gene co-expression network analysis and functional analysis were performed to explore the regulatory roles of those differentially expressed molecules. The cis-and trans-target genes of differently expressed lncRNAs were enriched for pathways related to focal adhesion, cGMP-PKG signaling pathway, alpha-linolenic acid metabolism, arachidonic acid metabolism, and fat digestion and absorption. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes analyses showed that the differently expressed genes mainly participated in mitotic cytokinesis, desmosome, fatty acid degradation, cell adhesion molecules, and fatty acid metabolism. The prediction of lncRNA-mRNA interaction networks further revealed transcripts potentially involved in rumen development. The present study profiles a global overview of lncRNAs and mRNAs during rumen development. Our findings provide valuable resources for genetic regulation and molecular mechanisms of rumen development in ruminants.