Evolutionary Selection and Constraint on Human Knee Chondrocyte Regulation Impacts Osteoarthritis Risk.

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.

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.

Maternal L-carnitine supplementation promotes brown adipose tissue thermogenesis of newborn goats after cold exposure.

Brown adipose tissue (BAT) is an important component of energy expenditure and necessary to maintain body temperature for newborn mammals. In the previous study, we found that L-carnitine was enriched in BAT and promoted BAT adipogenesis and thermogenesis in goat brown adipocytes. However, whether dietary L-carnitine regulates BAT heat production and energy expenditure in lambs remains unclear. In this study, maternal L-carnitine supplementation elevated the rectal temperature, as well as the expression of UCP1 and mitochondrial DNA content to promote BAT thermogenesis in newborn goats. Moreover, maternal L-carnitine supplementation increased the levels of triglycerides (TG), non-esterified fatty acids (NEFA), and lactate in plasma, as well as the content of lipid droplet and glycogen in BAT of newborn goats. Lipidomic analysis showed that maternal L-carnitine supplementation remodeled the lipid composition of BAT in newborn goats. L-carnitine significantly increased the levels of TG and diglyceride (DG) and decreased the levels of glycerophospholipids and sphingolipids in BAT. Further studies showed that L-carnitine promoted TG and glycogen deposition in brown adipocytes through AMPKα. Our results indicate that maternal L-carnitine supplementation promotes BAT development and thermogenesis in newborn goats and provides new evidence for newborn goats to maintain body temperature in response to cold exposure.

Expression patterns and DNA methylation profile of GTL2 gene in goats.

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.

The Profiles and Functions of RNA Editing Sites Associated with High-Altitude Adaptation in Goats.

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.

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.

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.

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.

Cold exposure induces lipid dynamics and thermogenesis in brown adipose tissue of goats.

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.

RNA-Seq reveals miRNA role in thermogenic regulation in brown adipose tissues of goats.

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.

A 13.42-kb tandem duplication at the ASIP locus is strongly associated with the depigmentation phenotype of non-classic Swiss markings in goats.

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.

Integrated application of metabolomics and RNA-seq reveals thermogenic regulation in goat brown adipose tissues.

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.

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.

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.

Genome-wide association study reveals 14 new SNPs and confirms two structural variants highly associated with the horned/polled phenotype in goats.

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.

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.

Dynamic transcriptomic analysis in hircine longissimus dorsi muscle from fetal to neonatal development stages.

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.