Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds.

Hu sheep (HS), a breed of sheep carrying the FecB mutation gene, is known for its "year-round estrus and multiple births" and is an ideal model for studying the high fecundity mechanisms of livestock. Through analyzing and comparing the genomic selection features of Hu sheep and other sheep breeds, we identified a series of candidate genes that may play a role in Hu sheep's high fecundity mechanisms. In this study, we conducted whole-genome resequencing on six breeds and screened key mutations significantly correlated with high reproductive traits in sheep. Notably, the CC2D1B gene was selected by the fixation index (FST) and the cross-population composite likelihood ratio (XP-CLR) methods in HS and other five breeds. It was worth noting that the CC2D1B gene in HS was different from that in other sheep breeds, and seven missense mutations have been identified. Furthermore, the linkage disequilibrium (LD) analysis revealed a strong linkage disequilibrium in this specific gene region. Subsequently, by performing different grouping based on FecB genotypes in Hu sheep, genome-wide selective signal analysis screened several genes related to reproduction, such as BMPR1B and PPM1K. Besides, FST analysis identified functional genes related to reproductive traits, including RHEB, HSPA2, PPP1CC, HVCN1, and CCDC63. Additionally, a missense mutation was found in the CCDC63 gene and the haplotype was different between the high reproduction (HR) group and low reproduction (LR) group in HS. In summary, we discovered genetic differentiation among six distinct breeding sheep breeds at the whole genome level. Additionally, we identified a set of genes which were associated with reproductive performance in Hu sheep and visualized how these genes differed in different breeds. These findings laid a theoretical foundation for understanding genetic mechanisms behind high prolific traits 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.

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.

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.

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.

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.

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.

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.

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.

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.

Systematic analyses reveal RNA editing events involved in skeletal muscle development of goat (Capra hircus).

RNA editing is a posttranscriptional molecular process involved with specific nucleic modification, which can enhance the diversity of gene products. Adenosine-to-inosine (A-to-I, I is read as guanosine by both splicing and translation machinery) is the main type of RNA editing in mammals, which manifested as AG (adenosine-to-guanosine) in sequence data. Here, we aimed to explore patterns of RNA editing using RNA sequencing data from skeletal muscle at four developmental stages (three fetal periods and one postnatal period) in goat. We found the occurrences of RNA editing events raised at fetal periods and declined at the postnatal period. Also, we observed distinct editing levels of AG editing across stages, and significant difference was found between postnatal period and fetal periods. AG editing patterns in newborn goats are similar to those of 45-day embryo compared with embryo at 105 days and embryo at 60 days. In this study, we found a total of 1415 significantly differential edited AG sites among four groups. Moreover, 420 sites were obviously clustered into six time-series profiles, and one profile had significant association between editing level and gene expression. Our findings provided some novel insights into understanding the molecular mechanism of muscle development in mammals.