Compendium of 5810 genomes of sheep and goat gut microbiomes provides new insights into the glycan and mucin utilization.

BACKGROUND: Ruminant gut microbiota are critical in ecological adaptation, evolution, and nutrition utilization because it regulates energy metabolism, promotes nutrient absorption, and improves immune function. To study the functional roles of key gut microbiota in sheep and goats, it is essential to construct reference microbial gene catalogs and high-quality microbial genomes database. RESULTS: A total of 320 fecal samples were collected from 21 different sheep and goat breeds, originating from 32 distinct farms. Metagenomic deep sequencing and binning assembly were utilized to construct a comprehensive microbial genome information database for the gut microbiota. We successfully generated the largest reference gene catalogs for gut microbiota in sheep and goats, containing over 162 million and 82 million nonredundant predicted genes, respectively, with 49 million shared nonredundant predicted genes and 1138 shared species. We found that the rearing environment has a greater impact on microbial composition and function than the host's species effect. Through subsequent assembly, we obtained 5810 medium- and high-quality metagenome-assembled genomes (MAGs), out of which 2661 were yet unidentified species. Among these MAGs, we identified 91 bacterial taxa that specifically colonize the sheep gut, which encode polysaccharide utilization loci for glycan and mucin degradation. CONCLUSIONS: By shedding light on the co-symbiotic microbial communities in the gut of small ruminants, our study significantly enhances the understanding of their nutrient degradation and disease susceptibility. Our findings emphasize the vast potential of untapped resources in functional bacterial species within ruminants, further expanding our knowledge of how the ruminant gut microbiota recognizes and processes glycan and mucins. Video Abstract.

Ame-miR-1-3p of bee venom reduced cell viability through the AZIN1/OAZ1-ODC1-polyamines pathway and enhanced the defense ability of honeybee (Apis mellifera L.).

Bee venom serves as an essential defensive weapon for bees and also finds application as a medicinal drug. MicroRNAs (miRNAs) serve as critical regulators and have been demonstrated to perform a variety of biological functions. However, the presence of miRNAs in bee venom needs to be confirmed. Therefore, we conducted small RNA sequencing and identified 158 known miRNAs, 15 conserved miRNAs and 4 novel miRNAs. It is noteworthy that ame-miR-1-3p, the most abundant among them, accounted for over a quarter of all miRNA reads. To validate the function of ame-miR-1-3p, we screened 28 candidate target genes using transcriptome sequencing and three target gene prediction software (miRanda, PITA and TargetScan) for ame-miR-1-3p. Subsequently, we employed real-time quantitative reverse transcription PCR (qRT-PCR), Western blot and other technologies to confirm that ame-miR-1-3p inhibits the relative expression of antizyme inhibitor 1 (AZIN1) by targeting the 3' untranslated region (UTR) of AZIN1. This, in turn, caused ODC antizyme 1 (OAZ1) to bind to ornithine decarboxylase 1 (ODC1) and mark ODC1 for proteolytic destruction. The reduction in functional ODC1 ultimately resulted in a decrease in polyamine biosynthesis. Furthermore, we determined that ame-miR-1-3p accelerates cell death through the AZIN1/OAZ1-ODC1-polyamines pathway. Our studies demonstrate that ame-miR-1-3p diminishes cell viability and it may collaborate with sPLA2 to enhance the defence capabilities of honeybees (Apis mellifera L.). Collectively, these data further elucidate the defence mechanism of bee venom and expand the potential applications of bee venom in medical treatment.

An Identification of Functional Genetic Variants in B4GALNT2 Gene and Their Association with Growth Traits in Goats.

ß-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.

Cold exposure affects glucose metabolism, lipid droplet deposition and mitophagy in skeletal muscle of newborn goats.

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.

Whole-Genome Resequencing Reveals Genetic Diversity and Wool Trait-Related Genes in Liangshan Semi-Fine-Wool Sheep.

Understanding the genetic makeup of local sheep breeds is essential for their scientific conservation and sustainable utilization. The Liangshan semi-fine-wool sheep (LSS), a Chinese semi-fine-wool breed renowned for its soft wool, was analyzed using whole-genome sequencing data including 35 LSS, 84 sheep from other domestic breeds, and 20 Asiatic mouflons. We investigated the genetic composition of LSS by conducting analyses of the population structure, runs of homozygosity, genomic inbreeding coefficients, and selection signature. Our findings indicated that LSS shares greater genetic similarity with Border Leicester and Romney sheep than with Tibetan (TIB), Yunnan (YNS), and Chinese Merino sheep. Genomic analysis indicated low to moderate inbreeding coefficients, ranging from 0.014 to 0.154. In identifying selection signals across the LSS genome, we pinpointed 195 candidate regions housing 74 annotated genes (e.g., IRF2BP2, BVES, and ALOX5). We also found the overlaps between the candidate regions and several known quantitative trait loci related to wool traits, such as the wool staple length and wool fiber diameter. A selective sweep region, marked by the highest value of cross-population extended haplotype homozygosity, encompassed IRF2BP2-an influential candidate gene affecting fleece fiber traits. Furthermore, notable differences in genotype frequency at a mutation site (c.1051 + 46T > C, Chr25: 6,784,190 bp) within IRF2BP2 were observed between LSS and TIB and YNS sheep (Fisher's exact test, p < 2.2 × 10-16). Taken together, these findings offer insights crucial for the conservation and breeding enhancement of LSS.

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

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

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.

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

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

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.

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.

A Novel LncRNA MSTRG.310246.1 Promotes Differentiation and Thermogenesis in Goat Brown Adipocytes.

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.

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.

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.

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

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

Diarrhea in suckling lambs is associated with changes in gut microbiota, serum immunological and biochemical parameters in an intensive production system.

The incidence of diarrhea in lambs is frequent in large-scale sheep farms, which greatly impacts the growth and health of lambs. The aim of this study was to assess the changes in serum biochemical and immunological parameters and gut microbiome composition in suckling lambs suffering from diarrhea or not, reared on an intensive commercial farm. We found a reduced diversity of intestinal bacteria in suckling lambs suffering from diarrhea. Firmicutes and Bacteroidetes were the dominant flora in both groups of lambs, while the Bacteroidetes decreased in diarrheic lambs, no changes were reported in Firmicutes. Compared with healthy lambs, the proportion of aerobic bacteria, facultative anaerobic bacteria, and stress tolerant bacteria increased in lambs suffering from diarrhea, while that of anaerobic bacteria and potentially pathogenic bacteria decreased slightly. In addition, the contents of total cholesterol, immunoglobulins (Ig) G, and IgM in the serum of lambs suffering from diarrhea were lower than those of healthy lambs. This study explored the association between diarrhea occurrence, intestinal microbial community structure, and metabolic and immunological status in Hu lambs.

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

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

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

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

Markhor-derived Introgression of a Genomic Region Encompassing PAPSS2 Confers High-altitude Adaptability in Tibetan Goats.

Understanding the genetic mechanism of how animals adapt to extreme conditions is fundamental to determine the relationship between molecular evolution and changing environments. Goat is one of the first domesticated species and has evolved rapidly to adapt to diverse environments, including harsh high-altitude conditions with low temperature and poor oxygen supply but strong ultraviolet radiation. Here, we analyzed 331 genomes of domestic goats and wild caprid species living at varying altitudes (high > 3000 m above sea level and low < 1200 m), along with a reference-guided chromosome-scale assembly (contig-N50: 90.4 Mb) of a female Tibetan goat genome based on PacBio HiFi long reads, to dissect the genetic determinants underlying their adaptation to harsh conditions on the Qinghai-Tibetan Plateau (QTP). Population genomic analyses combined with genome-wide association studies (GWAS) revealed a genomic region harboring the 3'-phosphoadenosine 5'-phosphosulfate synthase 2 (PAPSS2) gene showing strong association with high-altitude adaptability (PGWAS = 3.62 × 10-25) in Tibetan goats. Transcriptomic data from 13 tissues revealed that PAPSS2 was implicated in hypoxia-related pathways in Tibetan goats. We further verified potential functional role of PAPSS2 in response to hypoxia in PAPSS2-deficient cells. Introgression analyses suggested that the PAPSS2 haplotype conferring the high-altitude adaptability in Tibetan goats originated from a recent hybridization between goats and a wild caprid species, the markhor (Capra falconeri). In conclusion, our results uncover a hitherto unknown contribution of PAPSS2 to high-altitude adaptability in Tibetan goats on QTP, following interspecific introgression and natural selection.