Search

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. [ABSTRACT FROM AUTHOR]

Given GDF5 involvement in hip development, and osteoarthritis (OA) and developmental hip dysplasia (DDH) risk, here we sought to assess the role(s) of GDF5 and its regulatory sequence on the development of hip morphology linked to injury risk. The brachypodism (bp) mouse, which harbors a Gdf5 inactivating mutation, was used to survey how Gdf5 loss of function impacts the development of hip morphology. Two transgenic Gdf5 reporter BAC lines were used to assess the spatiotemporal expression of Gdf5 regulatory sequences. Each BAC line was also used to assess the functional roles of upstream and downstream sequence on hip morphology. bp/bp mice had shorter femora with smaller femoral heads and necks as well as larger alpha angles, smaller anterior offsets, and smaller acetabula, compared to bp/+ mice (p<0.04). Regulatory sequences downstream of Gdf5 drove strong prenatal (E17) expression and low postnatal (6 months) expression across regions of femoral head and acetabulum. Conversely, upstream regulatory sequences drove very low expression at E17 and no detectable expression at 6 months. Importantly, downstream, but not upstream Gdf5 regulatory sequences fully restored all the key morphologic features disrupted in bp/bp mice. Hip morphology is profoundly affected by Gdf5 absence, and downstream regulatory sequences mediate its effects by controlling Gdf5 expression during development. This downstream region contains numerous enhancers harboring risk variants related to hip OA, DDH, and dislocation. We posit that subtle alterations to morphology driven by changes in downstream regulatory sequence underlie this locus’ role in hip injury risk. [ABSTRACT FROM AUTHOR]

Simple Summary: Investigating the proliferation and differentiation of skeletal muscle satellite cells (MuSCs) in vitro is crucial for elucidating the underlying mechanisms of skeletal muscle development in goats. Real-time quantitative PCR (RT-qPCR) is a widely utilized technique for quantifying the expression levels of target genes. A significant challenge associated with this method is the identification of optimal reference genes for accurate normalization. In this study, we evaluated ten candidate reference genes for the standardization of gene expression. Our results indicated that RPL14 and RPS15A constituted the most stable reference gene combination during the proliferation and differentiation of goat MuSCs. The process of skeletal muscle development is intricate and involves the regulation of a diverse array of genes. Accurate gene expression profiles are crucial for studying muscle development, making it essential to choose the right reference genes for real-time quantitative PCR (RT-qPCR). In the present study, eight candidate reference genes were identified from our previous transcriptome sequencing analysis of caprine skeletal muscle satellite cells (MuSCs), and two traditional reference genes (ACTB and GAPDH) were assessed. The quantitative levels of the candidate reference genes were determined through the RT-qPCR technique, while the stability of their expression was evaluated utilizing the GeNorm, NormFinder, BestKeeper, and RefFinder programs. Furthermore, the chosen reference genes were utilized for the normalization of the gene expression levels of PCNA and Myf5. It was determined that conventional reference genes, including ACTB and GAPDH, were not appropriate for normalizing target gene expression. Conversely, RPL14 and RPS15A, identified through RNA sequencing analysis, exhibited minimal variability and were identified as the optimal reference genes for normalizing gene expression during the proliferation and differentiation of goat MuSCs. Our research offers a validated panel of optimal reference genes for the detection of differentially expressed genes in goat muscle satellite cells using RT-qPCR. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

β-1,4-N-acetylgalactosamine transferase 2 (B4GALNT2) is a vital candidate gene that affects the growth traits in sheep. However, whether it has the same function in goats remains to be investigated further. This study selected 348 Nanjiang Yellow goats, screened all exons, and conserved non-coding regions of the B4GALNT2 gene for single-nucleotide polymorphisms (SNPs). Our results revealed the presence of a synonymous mutation, rs672215506, within the exon of the B4GALNT2 gene in the Nanjiang Yellow goat population. The mutation resulted in a decrease in the mRNA stability of the B4GALNT2 gene. The results of SNP detection of the conserved non-coding region of the B4GALNT2 gene showed five potential regulatory SNPs in the Nanjiang Yellow goat population. Except for rs66095343, the ~500 bp fragments of the other four SNPs (rs649127714, rs649573228, rs652899012, and rs639183528) significantly increased the luciferase activity both in goat skeletal muscle satellite cells (MuSCs) and 293T cells. The genetic diversity indexes indicated low or intermediate levels for all six SNPs analyzed, and the genotype frequencies were in Hardy–Weinberg equilibrium. Association analysis showed that rs660965343, rs649127714, and rs649573228 significantly correlate with growth traits in the later stage of growth and development of Nanjiang Yellow goats. The haplotype combinations of H2H3 and H2H2 had higher body weight and greater body size. Moreover, H2H2 haplotype combinations significantly correlated with the litter size of the Nanjiang Yellow goats. The results of our study demonstrate the potential role of the B4GALNT2 gene as a functional genetic marker in the breeding programs of Nanjiang Yellow goats. [ABSTRACT FROM AUTHOR]

Gene trap locus 2 (GTL2), a long non-coding paternal imprinting gene, participates in various biological processes, including cell proliferation, differentiation, and apoptosis, by regulating the transcription of target mRNA, which is tightly related to the growth of the organic and maintenance of function. In this study, DNA methylation patterns of CpG islands (CGI) of GTL2 were explored, and its expression level was quantified in six tissues, rumen epithelium cells, and skeletal muscle cells in goats. GTL2 expression levels were measured by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), and the methylation model was confirmed by bisulfite-sequencing PCR (BSP). CGI methylation of GTL2 indicated a moderate methylation (ranging from 81.42 to 86.83%) in the brain, heart, liver, kidney, lung, and longissimus dorsi. GTL2 is most highly expressed in brain tissues, but there is no significant difference in the other five tissues. In addition, in the rumen epithelium cell proliferation, GTL2 expression was highest at 60 h, followed by 72 h, and almost unchanged at 12–48 h. In the skeletal muscle cell differentiation, GTL2 expression was highest at 0 and 24 h, significantly decreasing at 72 and 128 h. Pearson correlation analysis did not indicate a clear relationship between methylation and GTL2 expression levels, suggesting that other regulatory factors may modulate GTL2 expression. This study will provide a better understanding of the expression regulation mechanism of genes in the delta-like homolog 1 gene (DLK1)-GTL2 domain. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Brown adipose tissue (BAT) plays a critical role in maintaining the body temperature in newborn lamb due to its unique non-shivering thermogenesis. Previous studies have found that BAT thermogenesis is regulated by several long non-coding RNAs (lncRNAs). Here, we identified a novel lncRNA, MSTRG.310246.1, which was enriched in BAT. MSTRG.310246.1 was localized in both the nuclear and cytoplasmic compartments. In addition, MSTRG.310246.1 expression was upregulated during brown adipocyte differentiation. Overexpression of MSTRG.310246.1 increased the differentiation and thermogenesis of goat brown adipocytes. On the contrary, the knockdown of MSTRG.310246.1 inhibited the differentiation and thermogenesis of goat brown adipocytes. However, MSTRG.310246.1 had no effect on goat white adipocyte differentiation and thermogenesis. Our results show that MSTRG.310246.1 is a BAT-enriched LncRNA that improves the differentiation and thermogenesis of goat brown adipocytes. [ABSTRACT FROM AUTHOR]

High-altitude environments dramatically influenced the genetic evolution of vertebrates. However, little is known about the role of RNA editing on high-altitude adaptation in non-model species. Here, we profiled the RNA editing sites (RESs) of heart, lung, kidney, and longissimus dorsi muscle from Tibetan cashmere goats (TBG, 4500 m) and Inner Mongolia cashmere goats (IMG, 1200 m) to reveal RNA editing-related functions of high-altitude adaptation in goats. We identified 84,132 high-quality RESs that were unevenly distributed across the autosomes in TBG and IMG, and more than half of the 10,842 non-redundant editing sites were clustered. The majority (62.61%) were adenosine-to-inosine (A-to-I) sites, followed by cytidine-to-uridine (C-to-U) sites (19.26%), and 32.5% of them had a significant correlation with the expression of catalytic genes. Moreover, A-to-I and C-to-U RNA editing sites had different flanking sequences, amino acid mutations, and alternative splicing activity. TBG had higher editing levels of A-to-I and C-to-U than IMG in the kidney, whereas a lower level was found in the longissimus dorsi muscle. Furthermore, we identified 29 IMG and 41 TBG population-specific editing sites (pSESs) and 53 population-differential editing sites (pDESs) that were functionally involved in altering RNA splicing or recoding protein products. It is worth noting that 73.3% population-differential, 73.2% TBG-specific, and 80% IMG-specific A-to-I sites were nonsynonymous sites. Moreover, the pSESs and pDESs editing-related genes play critical functions in energy metabolisms such as ATP binding molecular function, translation, and adaptive immune response, which may be linked to goat high-altitude adaptation. Our results provide valuable information for understanding the adaptive evolution of goats and studying plateau-related diseases. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

[Display omitted] • Cold chain storage darkens mutton color and increases C18:2n-6, C20:3n-6, and C23:0 contents. • 1,2-Dioleoyl- sn -glycero-3-phosphoethanolamine and alanylphenylalanine are the top two altered metabolites during cold storage. • Aeromonas , Brachymonas , Corynebacterium , and Steroidobacter are dominated in the frozen mutton. • The identified different metabolites and microorganisms in mutton are highly correlated during cold chain storage. During the cold chain storage process, changes in metabolites and microorganisms are highly likely to lead to changes in meat quality. To elucidate the changes in the composition of metabolites and microbiota during cold chain storage of mutton, this study utilized untargeted metabolome and 5R 16S rRNA sequencing analyses to investigate the changes in the longissimus dorsi under different cold chain temperatures (4 °C and −20 °C). With the extension of cold chain storage time, the meat color darkened and the content of C18:2n-6, C20:3n-6, and C23:0 were significantly increased in mutton. In this study, nine metabolites, including 1,2-Dioleoyl- sn -glycero-3-phosphoethanolamine, alanylphenylala-nine, indole-3-acrylic acid and the others, were significantly altered during cold chain storage. The abundance of the dominant microorganisms, including Brachymonas , Aeromonas , Corynebacterium and Steroidobacter , was significantly altered. Furthermore, a high correlation was observed between the different metabolites and microorganisms. These findings provide an in-depth understanding of the effects of different cold chain storage temperatures and times on the quality of mutton. [ABSTRACT FROM AUTHOR]

Background: Adaptive thermogenesis by brown adipose tissue (BAT) is important to the maintenance of temperature in newborn mammals. Cold exposure activates gene expression and lipid metabolism to provide energy for BAT thermogenesis. However, knowledge of BAT metabolism in large animals after cold exposure is still limited. Results: In this study, we found that cold exposure induced expression of BAT thermogenesis genes and increased the protein levels of UCP1 and PGC1α. Pathway analysis showed that cold exposure activated BAT metabolism, which involved in cGMP-PKG, TCA cycle, fatty acid elongation, and degradation pathways. These were accompanied by decreased triglyceride (TG) content and increased phosphatidylcholine (PC) and phosphatidylethanolamine (PE) content in BAT. Conclusion: These results demonstrate that cold exposure induces metabolites involved in glycerolipids and glycerophospholipids metabolism in BAT. The present study provides evidence for lipid composition associated with adaptive thermogenesis in goat BAT and metabolism pathways regulated by cold exposure. [ABSTRACT FROM AUTHOR]

Background: The pigmentation phenotype diversity is rich in domestic goats, and identification of the genetic loci affecting coat color in goats has long been of interest. Via the detections of selection signatures, a duplication upstream ASIP was previously reported to be a variant affecting the Swiss markings depigmentation phenotype in goats. Results: We conducted a genome-wide association study using whole-genome sequencing (WGS) data to identify the genetic loci and causal variants affecting the pigmentation phenotype in 65 Jintang black (JT) goats (i.e., 48 solid black vs. 17 non-classic Swiss markings). Although a single association peak harboring the ASIP gene at 52,619,845–72,176,538 bp on chromosome 13 was obtained using a linear mixed model approach, all the SNPs and indels in this region were excluded as causal variants for the pigmentation phenotype. We then found that all 17 individuals with non-classic Swiss markings carried a 13,420-bp duplication (CHI13:63,129,198–63,142,617 bp) nearly 101 kb upstream of ASIP, and this variant was strongly associated (P = 1.48 × 10− 12) with the coat color in the 65 JT goats. The copy numbers obtained from the WGS data also showed that the duplication was present in all 53 goats from three European breeds with Swiss markings and absent in 45 of 51 non-Swiss markings goats from four other breeds and 21 Bezoars, which was further validated in 314 samples from seven populations based on PCR amplification. The copy numbers of the duplication vary in different goat breeds with Swiss markings, indicating a threshold effect instead of a dose-response effect at the molecular level. Furthermore, breakpoint flanking repeat analysis revealed that the duplication was likely to be a result of the Bov-B-mediated nonallelic homologous recombination. Conclusion: We confirmed that a genomic region harboring the ASIP gene is a major locus affecting the coat color phenotype of Swiss markings in goats. Although the molecular genetic mechanisms remain unsolved, the 13,420-bp duplication upstream of ASIP is a necessary but not sufficient condition for this phenotype in goats. Moreover, the variations in the copy number of the duplication across different goat breeds do not lead to phenotypic heterogeneity. [ABSTRACT FROM AUTHOR]

• Cold exposure increased energy metabolism, resulting in decreased glycogen content in skeletal muscle of newborn goats. • Transcriptome analysis found that cold exposure up-regulated the expression levels of genes related to glucose uptake, glycolysis, and fatty acid and TG synthesis, promoting lipid deposition in skeletal muscle. • Cold exposure induced mitophagy in skeletal muscle. Cold exposure is a common stressor for newborn goats. Skeletal muscle plays an important role in maintaining whole-body homeostasis of glucose and lipid metabolism. However, the molecular mechanisms underlying regulation of skeletal muscle of newborn goats by cold exposure remains unclear. In this study, we found a significant increase (P < 0.01) in serum glucagon levels after 24 h of cold exposure (COLD, 6°C), while glucose and insulin concentrations were significantly decreased (P < 0.01) compared to room temperature (RT, 25°C). Additionally, we found that cold exposure reduced glycogen content (P < 0.01) in skeletal muscle. Pathway enrichment analysis revealed that cold exposure activated skeletal muscle glucose metabolism pathways (including insulin resistance and the insulin signaling pathway) and mitophagy-related pathways. Cold exposure up-regulated the expression of genes involved in fatty acid and triglyceride synthesis, promoting skeletal muscle lipid deposition. Notably, cold exposure induced mitophagy in skeletal muscle. [ABSTRACT FROM AUTHOR]

Cell Res. 2002; 12: 9-18. 10.1038/sj.cr.7290105 41 Xie S.J., Li J.H., Chen H.F., Tan Y.Y., Liu S.R., Zhang Y., Xu H., Yang J.H., Liu S., Zheng L.L. Inhibition of the JNK/MAPK signaling pathway by myogenesis-associated miRNAs is required for skeletal muscle development. Many gene expression profiling studies in agriculturally significant animals have been carried out since the introduction of high-throughput sequencing tools, and many genes or factors associated to skeletal muscle development have been found [[8], [10]]. Although progress has been made in understanding how genes work during myogenesis in pigs [[8], [10]], cattle [[27]], goats [[11], [28]], and sheep [[9]], our understanding of the genes involved in goat skeletal muscle development is still limited and incomplete. Among them, we indeed identified many terms related to muscle development, such as muscle structure development, muscle tissue development, muscle contraction, muscle cell development, and muscle cell differentiation. [Extracted from the article]

Myogenesis is a complex process controlled by several coding and non-coding RNAs (ncRNAs), such as circular RNAs (circRNAs) that are known to function as endogenous microRNAs (miRNAs) sponges. Cerebellar Degeneration-Related protein 1 antisense (CDR1as) is the most spotlighted circRNA that is known as an miR-7 sponge, which has bloomed circRNAs' research in animal disease and physiology. Here, we screened for miRNAs and mRNA associated with CDR1as and further characterized their regulatory function during muscle differentiation. We found that a total of 43 miRNAs (including miR-107-3p, miR-125b-5p, miR-140-5p, miR-29a-3p, and miR-27a-3p upregulated) and 789 mRNAs (including ANGPT1, E2F2, CCN1, FGFR1, and MEF2C downregulated) were differentially expressed in goat skeletal muscle satellite cells (SMSCs). Further, knockdown of CDR1as and ANGPT1 inhibited SMSCs differentiation. miR-27a-3p was differentially upregulated after the knockdown of CDR1as in SMSCs. Overexpressed miR-27a-3p decreased SMSCs differentiation. Via RNAhybrid and luciferase, miR-27a-3p was identified to regulate ANGPT1. We discovered that miR-27a-3p has an inverse relationship with CDR1as and decreases the expression level of ANGPT1 during SMSCs differentiation. In summary, our study demonstrates that siCDR1as inhibits myoblast differentiation by downregulating ANGPT1 mRNA via miR-27a-3p in SMSCs. [ABSTRACT FROM AUTHOR]

Background: MicroRNAs (miRNAs) are a family of short non-coding RNA molecules and play important roles in various biological processes. However, knowledge of the expression profiles and function of miRNAs on the regulation of brown adipose tissue (BAT) thermogenesis remains largely unknown. Results: In this study, we found that brown adipose tissue (BAT) existed within the perirenal fat at 1 day after birth (D1) and transferred into white adipose tissue (WAT) at 30 days after birth (D30) by UCP1 protein expression and immunohistochemistry analysis. After that, we performed RNA sequencing on six libraries of goat BAT and WAT. A total of 238 known miRNAs and 1834 goat novel miRNAs were identified. Moreover, 395 differentially expressed miRNAs including 167 up-regulated and 228 down-regulated miRNAs were obtained in BAT. For the known BAT enriched miRNA, 30 miRNAs were enriched in goat BAT but not in mouse BAT. In addition, miR-433 was enriched in goat BAT but not in mouse BAT. Gain- and loss-of-function experiments reveal that miR-433 reduced the lipid accumulation of brown adipocytes and decreased the expression of BAT marker and mitochondrial related genes. However, miR-433 had no effect on lipid accumulation and thermogenesis in white adipocytes. In addition, miR-433 inhibited the expression of MAPK8 by targeting to the 3'UTR of MAPK8 gene. These data demonstrate that miR-433 acts as a negative regulator in controlling brown adipocytes differentiation and thermogenesis. Conclusion: The present study provides a detailed miRNAs expression landscape in BAT and WAT. Furthermore, we found that miR-433, which was highly expressed on BAT had a negative regulatory function on the thermogenesis and adipogenesis in goat brown adipocytes. This study provides evidence for understanding the role of miRNAs in regulating BAT thermogenesis and energy expenditure in goats. [ABSTRACT FROM AUTHOR]

Background: There is a long-term interest in investigating the genetic basis of the horned/polled phenotype in domestic goats. Here, we report a genome-wide association study (GWAS) to detect the genetic loci affecting the polled phenotype in goats. Results: We obtained a total of 13,980,209 biallelic SNPs, using the genotyping-by-sequencing data from 45 Jintang Black (JT) goats, which included 32 female and nine male goats, and four individuals with the polled intersex syndrome (PIS). Using a mixed-model based GWAS, we identified two association signals, which were located at 150,334,857–150,817,260 bp (P = 5.15 × 10− 119) and 128,286,704–131,306,537 bp (P = 2.74 × 10− 15) on chromosome 1. The genotype distributions of the 14 most significantly associated SNPs were completely correlated with horn status in goats, based on the whole-genome sequencing (WGS) data from JT and two other Chinese horned breeds. However, variant annotation suggested that none of the detected SNPs within the associated regions were plausible causal mutations. Via additional read-depth analyses and visual inspections of WGS data, we found a 10.1-kb deletion (CHI1:g. 129424781_129434939del) and a 480-kb duplication (CHI1:150,334,286–150,818,098 bp) encompassing two genes KCNJ15 and ERG in the associated regions of polled and PIS-affected goats. Notably, the 10.1-kb deletion also served as the insertion site for the 480-kb duplication, as validated by PCR and Sanger sequencing. Our WGS genotyping showed that all horned goats were homozygous for the reference alleles without either the structural variants (SVs), whereas the PIS-affected goats were homozygous for both the SVs. We also demonstrated that horned, polled, and PIS-affected individuals among 333 goats from JT and three other Chinese horned breeds can be accurately classified via PCR amplification and agarose gel electrophoresis of two fragments in both SVs. Conclusion: Our results revealed that two genomic regions on chromosome 1 are major loci affecting the polled phenotypes in goats. We provided a diagnostic PCR to accurately classify horned, polled, and PIS-affected goats, which will enable a reliable genetic test for the early-in-life prediction of horn status in goats. [ABSTRACT FROM AUTHOR]

Brown adipose tissue (BAT) plays an important role on no shivering thermogenesis during cold exposure to maintain animal body temperature and energy homeostasis. However, knowledge of the cellular transition from white adipose tissue (WAT) to BAT is still limited. In this study, we provided a comprehensive metabolomics and transcriptional signatures of goat BAT and WAT. A total of 157 metabolites were significantly changed, including 81 upregulated and 76 downregulated metabolites. In addition, we identified the citric acid cycle, fatty acid elongation, and degradation pathways as coordinately activated in BAT. Interestingly, five unsaturated fatty acids (Eicosadienoic Acid, C20:2; γ‐Linolenic acid, C20:3; Arachidonic Acid, C20:4; Adrenic acid, C22:4; Docosahexaenoic acid, C22:6), Succinate, L‐carnitine, and L‐palmitoyl‐carnitine were found to be abundant in BAT. Furthermore, L‐carnitine, an intermediate of fatty acid degradation, is required for goat brown adipocyte differentiation and thermogenesis through activating AMPK pathway. However, L‐carnitine decreased lipid accumulation through inducing lipolysis and thermogenesis in white adipocytes. These results revealed that there are the significant alterations in transcriptomic and metabolomic profiles between goat WAT and BAT, which may contribute to better understanding the roles of metabolites in BAT thermogenesis process. [ABSTRACT FROM AUTHOR]

Insulin‐like growth factor 2 mRNA‐binding protein 1 (IGF2BP1) plays essential roles in the proliferation of skeletal muscle satellite cells (MuSCs). Increasing evidence has shown that IGF2BP1 regulates the expression of noncoding RNAs and mRNAs. However, the related molecular network remains to be fully understood. Therefore, we performed RNA sequencing and analyzed the microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and mRNAs differentially expressed in goat MuSCs treated with IGF2BP1 overexpressing and empty vectors. A total of 36 miRNAs, 59 lncRNAs, and 44 mRNAs were differentially expressed caused by IGF2BP1. Expectedly, they were enriched in muscle development‐related Rap1, PI3K‐AKT, and FoxO signaling pathways. Finally, we constructed a lncRNA‐miRNA‐mRNA interaction network containing 30 lncRNAs, 15 miRNAs, and 34 mRNAs, in which several miRNAs, including miR‐133a‐3p, miR‐204‐5p, miR‐125a‐3p, miR‐145‐3p, and miR‐423‐5p, relate with cell growth and participate in muscle development. Overall, we constructed an IGF2BP1‐related network, which provides new insight into the myogenic proliferation of goat. [ABSTRACT FROM AUTHOR]

During human evolution, the knee adapted to the biomechanical demands of bipedalism by altering chondrocyte developmental programs. This adaptive process was likely not without deleterious consequences to health. Today, osteoarthritis occurs in 250 million people, with risk variants enriched in non-coding sequences near chondrocyte genes, loci that likely became optimized during knee evolution. We explore this relationship by epigenetically profiling joint chondrocytes, revealing ancient selection and recent constraint and drift on knee regulatory elements, which also overlap osteoarthritis variants that contribute to disease heritability by tending to modify constrained functional sequence. We propose a model whereby genetic violations to regulatory constraint, tolerated during knee development, lead to adult pathology. In support, we discover a causal enhancer variant (rs6060369) present in billions of people at a risk locus (GDF5 - UQCC1), showing how it impacts mouse knee-shape and osteoarthritis. Overall, our methods link an evolutionarily novel aspect of human anatomy to its pathogenesis. • Human/mouse long bone chondrocyte ATAC-seq reveals site-specific regulatory elements • Ancient positive and purifying selection shape human knee regulation and morphology • Modern human drift and pleiotropy violate regulatory constraint conferring OA risk • Validation at GDF5 reveals a single regulatory variant impacting both shape and OA Epigenetic profiling of chondrocytes reveals selection and constraint on knee regulatory elements including those overlapping osteoarthritis risk variants, including a functionally verified common enhancer variant sufficient to affect knee shape in mouse. [ABSTRACT FROM AUTHOR]

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. • H19 promotes differentiation of goat skeletal muscle satellite cells (MuSCs). • H19 partially colocalizes with MyoD transcripts. • H19 significantly stabilizes MyoD mRNA. • MyoD activates H19 by binding on its promoter and exons. [ABSTRACT FROM AUTHOR]

As a well-conserved microRNA, miR-183 is ubiquitously expressed in many tissues and cells including backfat and the 3T3-L1 adipocytes; however, the mechanisms regulating miR-183 in adipogenesis remain poorly understood. Here, we explored the expression pattern and role of miR-183 in adipogenesis using hircine preadipocytes. The results showed that miR-183 was up-regulated during preadipocyte differentiation, and overexpression of miR-183 enhanced lipid accumulation and dramatically increased the mRNA expression of the adipogenic genes PPARγ , C/EBPα , SREBP-1c , FAS , and ACC . Using bioinformatics tools, we predicted Smad4 to be a target of miR-183. This was subsequently validated with a luciferase reporter assay. Overexpression of miR-183 suppressed the mRNA and protein levels of Smad4 significantly, whereas inhibiting miR-183 had the opposite effect. However, inhibition of Smad4 greatly accelerated lipid deposition and increased the expression of adipogenic genes which consists with the results of miR-183 overexpression. In conclusion, these results indicate that miR-183 promotes hircine preadipocyte differentiation by targeting Smad4 . [ABSTRACT FROM AUTHOR]

Muscle growth and development from fetal to neonatal stages consist of a series of delicately regulated and orchestrated changes in expression of genes. In this study, we performed whole transcriptome profiling based on RNA-Seq of caprine longissimus dorsi muscle tissue obtained from prenatal stages (days 45, 60, and 105 of gestation) and neonatal stage (the 3-day-old newborn) to identify genes that are differentially expressed and investigate their temporal expression profiles. A total of 3276 differentially expressed genes (DEGs) were identified ( Q value < 0.01). Time-series expression profile clustering analysis indicated that DEGs were significantly clustered into eight clusters which can be divided into two classes ( Q value < 0.05), class I profiles with downregulated patterns and class II profiles with upregulated patterns. Based on cluster analysis, GO enrichment analysis found that 75, 25, and 8 terms to be significantly enriched in biological process (BP), cellular component (CC), and molecular function (MF) categories in class I profiles, while 35, 21, and 8 terms to be significantly enriched in BP, CC, and MF in class II profiles. KEGG pathway analysis revealed that DEGs from class I profiles were significantly enriched in 22 pathways and the most enriched pathway was Rap1 signaling pathway. DEGs from class II profiles were significantly enriched in 17 pathways and the mainly enriched pathway was AMPK signaling pathway. Finally, six selected DEGs from our sequencing results were confirmed by qPCR. Our study provides a comprehensive understanding of the molecular mechanisms during goat skeletal muscle development from fetal to neonatal stages and valuable information for future studies of muscle development in goats. [ABSTRACT FROM AUTHOR]

Circular RNA (circRNA) is involved in various biological processes. However, the regulatory roles of circRNA on brown adipose tissue (BAT) thermogenesis still remain to be elucidated. In this study, we characterized circZEB1, which was enriched in BAT and derived from exon 2 and exon 3 of goat ZEB1 gene. Sanger sequencing and RNase R assay demonstrated that circZEB1 was a head-to-tail structure with a length of 298 bp. Besides, the expression levels of circZEB1 were upregulated during brown adipocytes differentiation and circZEB1 was mainly located in cytoplasm. Furthermore, circZEB1 promoted brown adipocytes differentiation and thermogenesis. Luciferase reporter assay and RNA-binding protein immunoprecipitation (RIP) assay indicated that circZEB1 could bind to miR-326–3p in brown adipocytes. In addition, miR-326–3p reduced the lipid accumulation of brown adipocytes, whereas overexpression of circZEB1 reversed the inhibition effects of miR-326–3p on lipid deposition. The results demonstrated that circZEB1 regulated brown adipocytes differentiation and thermogenesis through circZEB1/miR-326–3p pathway. • We characterized circZEB1, which was enriched in BAT and derived from exon 2 and exon 3 of goat ZEB1 gene. • CircZEB1 promoted brown adipocytes differentiation and thermogenesis. • CircZEB1 regulated brown adipocytes differentiation and thermogenesis through circZEB1/miR-326–3p pathway. [ABSTRACT FROM AUTHOR]

An amendment to this paper has been published and can be accessed via the original article. [ABSTRACT FROM AUTHOR]

Simple Summary: The rumen plays an essential role as a digestive organ and serves as the primary site of energy substrate absorption for the productive ruminants. Understanding gene expression profiles is necessary to explore the intrinsic regulatory mechanisms of rumen development in goats. The selection of suitable reference genes (RGs) was the primary assay before the real-time quantitative PCR (RT-qPCR). We identified sixteen genome-wide candidate RGs for normalization of gene expression assessments in goat rumen tissues. We demonstrate that the RGs selected (RPS4X and RPS6) were more stably expressed than the commonly used HKGs (ACTB and GAPDH) in goat rumen tissues, suggesting that the ribosomal protein gene family may be another source for the RG pool. As the largest chamber of the ruminant stomach, the rumen not only serves as the principal absorptive surface and nutrient transport pathway from the lumen into the animal, but also plays an important short-chain fatty acid (SCFA) metabolic role in addition to protective functions. Accurate characterization of the gene expression profiles of genes of interest is essential to the exploration of the intrinsic regulatory mechanisms of rumen development in goats. Thus, the selection of suitable reference genes (RGs) is an important prerequisite for real-time quantitative PCR (RT-qPCR). In the present study, 16 candidate RGs were identified from our previous transcriptome sequencing of caprine rumen tissues. The quantitative expressions of the candidate RGs were measured using the RT-qPCR method, and the expression stability of the RGs was assessed using the geNorm, NormFinder, and BestKeeper programs. GeNorm analysis showed that the M values were less than 0.5 for all the RGs except GAPT4, indicating that they were stably expressed in the rumen tissues throughout development. RPS4X and RPS6 were the two most stable RGs. Furthermore, the expressions of two randomly selected target genes (IGF1 and TOP2A), normalized by the selected most stable RGs (RPS4X and RPS6), were consistent with the results of RNA sequencing, while the use of GAPDH and ACTB as RGs resulted in altered profiles. Overall, RPS4X and RPS6 showed the highest expression stability and the lowest coefficients of variation, and could be used as the optimal reference combination for quantifying gene expression in rumen tissues via RT-qPCR analysis. [ABSTRACT FROM AUTHOR]

Simple Summary: In our study, we launched RNA-seq in order to investigate the potential functions of circRNA during brown adipose tissue (BAT) to white adipose tissue (WAT) transformation. As a result, 6610 circRNAs and 61 differentially expressed circRNAs (DEcircRNAs) were identified. Moreover, 65 miRNAs were detected that could potentially interact with DEcircRNAs. The present study provides a detailed circRNA expression landscape and evidence for circRNA functions in the transformation from BAT to WAT. Adipose tissues are mainly divided into brown adipose tissue (BAT) and white adipose tissue (WAT). WAT mainly functions to buffer excess calories, whereas BAT plays a role in the non-shivering thermogenesis to maintain body temperature and energy balance. Moreover, circRNAs play important roles in various biological processes. However, knowledge of the expression profile and function of circRNAs from BAT to WAT remains largely unknown. In this study, a total of 6610 unique circRNAs were identified in the perirenal adipose tissues of 1-day, 30-days, and 1-year goats. Functional annotation revealed that host genes of circRNAs were involved in some BAT-related pathways, such as the thyroid hormone signaling pathway, MAPK signaling pathway, and VEGF signaling pathway. Furthermore, a total of 61 DEcircRNAs were detected across three stages. Additionally, five selected circRNAs were validated by RNase R assay, qPCR, and Sanger sequencing. Finally, the circRNA–miRNA network was constructed between the DEcircRNAs and their miRNA binding sites. [ABSTRACT FROM AUTHOR]

Simple Summary: In ruminants, the rumen epithelium plays an important role in nutrient absorption, metabolism and transport. MicroRNAs (miRNAs) have been reported to regulate the proliferation of diverse epithelial cells. In this study, we profiled the miRNA transcriptomes of goat rumens at four development stages and screened for candidate miRNAs related to rumen development. MiR-148a-3p was found to be highly expressed in the rumen tissues and induced the proliferation of GES-1 cells by targeting QKI. Our findings provide some insights into the functional roles of miRNAs in rumen growth and functional development in ruminants. The rumen is an important digestive organ of ruminants. From the fetal to adult stage, the morphology, structure and function of the rumen change significantly. However, the knowledge of the intrinsic genetic regulation of these changes is still limited. We previously reported a genome-wide expression profile of miRNAs in pre-natal goat rumens. In this study, we combined and analyzed the transcriptomes of rumen miRNAs during pre-natal (E60 and E135) and post-natal (D30 and D150) stages. A total of 66 differentially expressed miRNAs (DEMs) were identified in the rumen tissues from D30 and D150 goats. Of these, 17 DEMs were consistently highly expressed in the rumens at the pre-weaning stages (E60, E135 and D30), while down-regulated at D150. Noteworthy, annotation analysis revealed that the target genes regulated by the DEMs were mainly enriched in MAPK signaling pathway, Jak-STAT signaling pathway and Ras signaling pathway. Interestingly, the expression of miR-148a-3p was significantly high in the embryonic stage and down-regulated at D150. The potential binding sites of miR-148a-3p in the 3′-UTR of QKI were predicted by the TargetScan and verified by the dual luciferase report assay. The co-localization of miR-148a-3p and QKI through in situ hybridization was observed in the rumen tissues but not in the intestinal tracts. Moreover, the expression of miR-148a-3p in the epithelium was significantly higher than that in the other layers of the rumen, suggesting that miR-148a-3p is involved in the development of the rumen epithelial cells by targeting QKI. Subsequently, miR-148a-3p inhibitor was found to induce the proliferation of GES-1 cells. Taken together, our study identified DEMs involved in the development of the rumen and provides insights into the regulation mechanism of rumen development in goats. [ABSTRACT FROM AUTHOR]

Simple Summary: Brown adipose tissue (BAT) plays important roles in unique non-shivering thermogenesis. It is necessary to select reference genes during the transition process from brown (BAT) to white adipose tissue (WAT) for Quantitative PCR (qPCR) analysis. In this study, CTNNB, PFDN5 and EIF3M, selected from RNA sequencing data, were the most suitable reference genes. The present study provides a detailed analysis of the expression stability of reference genes for the study of gene expression profiling during the transition process from BAT to WAT. Brown adipose tissues have unique non-shivering thermogenesis functions, can be found in newborn ruminate animals, and then are gradually replaced by white adipose tissues in adulthood. For the purpose of exploring the intrinsic mechanism underlying the conversion process from brown (BAT) to white adipose tissue (WAT), it is necessary to utilize Quantitative PCR (qPCR) to study gene expression profiling. In this study, we identified reference genes that were consistently expressed during the transformation from goat BAT to WAT using RNA-seq data. Then, twelve genes were evaluated as candidate reference genes for qPCR in goat perirenal adipose tissue using three tools (geNorm, Normfinder, and BestKeeper). In addition, the selected reference genes were used to normalize the gene expression of PGC-1α and GPAT4. It was found that traditional reference genes, such as GAPDH, RPLP0, HPRT1, and PPIA were not suitable for target gene normalization. In contrast, CTNNB, PFDN5, and EIF3M, selected from RNA sequencing data, showed the least variation and were recommended as the best reference genes during the transformation from BAT to WAT. [ABSTRACT FROM AUTHOR]