Search

Abstract Background Circular RNAs (circRNAs) function as essential regulatory elements with pivotal roles in various biological processes. However, their expression profiles and functional regulation during the differentiation of goat myoblasts have not been thoroughly explored. This study conducts an analysis of circRNA expression profiles during the proliferation phase (cultured in growth medium, GM) and differentiation phase (cultured in differentiation medium, DM1/DM5) of skeletal muscle satellite cells (MuSCs) in goats. Results A total of 2,094 circRNAs were identified, among which 84 were differentially expressed as determined by pairwise comparisons across three distinct groups. Validation of the expression levels of six randomly selected circRNAs was performed using reverse transcription PCR (RT-PCR) and quantitative RT-PCR (qRT-PCR), with confirmation of their back-splicing junction sites. Enrichment analysis of the host genes associated with differentially expressed circRNAs (DEcircRNAs) indicated significant involvement in biological processes such as muscle contraction, muscle hypertrophy, and muscle tissue development. Additionally, these host genes were implicated in key signaling pathways, including Hippo, TGF-beta, and MAPK pathways. Subsequently, employing Cytoscape, we developed a circRNA-miRNA interaction network to elucidate the complex regulatory mechanisms underlying goat muscle development, encompassing 21 circRNAs and 47 miRNAs. Functional assays demonstrated that circTGFβ2 enhances myogenic differentiation in goats, potentially through a miRNA sponge mechanism. Conclusion In conclusion, we identified the genome-wide expression profiles of circRNAs in goat MuSCs during both proliferation and differentiation phases, and established that circTGFβ2 plays a role in the regulation of myogenesis. This study offers a significant resource for the advanced exploration of the biological functions and mechanisms of circRNAs in the myogenesis of goats.

Abstract Lacking effective therapeutic targets heavily restricts the improvement of clinical prognosis for patients diagnosed with esophageal squamous cell carcinoma (ESCC). Ubiquitin Specific Peptidase 21 (USP21) is dysregulated in plenty of human cancers, however, its potential function and relevant molecular mechanisms in ESCC malignant progression as well as its value in clinical translation remain largely unknown. Here, in vitro and in vivo experiments revealed that aberrant upregulation of USP21 accelerated the proliferation and metastasis of ESCC in a deubiquitinase-dependent manner. Mechanistically, we found that USP21 binds to, deubiquitinates, and stabilizes the G3BP Stress Granule Assembly Factor 1 (G3BP1) protein, which is required for USP21-mediated ESCC progression. Further molecular studies demonstrated that the USP21/G3BP1 axis played a tumor-promoting role in ESCC progression by activating the Wnt/β-Catenin signaling pathway. Additionally, disulfiram (DSF), an inhibitor against USP21 deubiquitylation activity, markedly abolished the USP21-mediated stability of G3BP1 protein and significantly displayed an anti-tumor effect on USP21-driving ESCC progression. Finally, the regulatory axis of USP21/G3BP1 was demonstrated to be aberrantly activated in ESCC tumor tissues and closely associated with advanced clinical stages and unfavorable prognoses, which provides a promising therapeutic strategy targeting USP21/G3BP1 axis for ESCC patients.

Abstract 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

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

Understanding the colonization and change patterns of gut microbiota is pivotal for comprehending host health. As a newly cultured breed, the studies on the gut microbiota of Tianfu goats remain limited. This study aimed to address this gap by analyzing the microbial composition and colonization patterns of fecal samples collected from goat kids from birth to weaning. Fecal samples were collected on days 0, 7, 14, 21, 28, 35, 42, 49, 53, 55, 57, and 64, and the changes and colonization patterns of microorganisms were analyzed through high-throughput 16S rRNA sequencing. The results showed that the abundance of fecal microbiota in goat kids gradually increased over time, followed by a decrease after weaning and stabilization, with reduced individual differences. The colonization of fecal microorganisms mainly presented three different stages: days 0–14, days 21–49, and days 53–64. During the suckling period, the relative abundance of Proteobacteria (72.34%) was the highest, followed by Firmicutes (21.66%). From 21 days old, the microbiota in goat kids gradually to be diverse, with Lachnospiraceae and Ruminococcaceae being dominant. During post-weaning, Ruminococcaceae (30.98–33.34%) was becoming prominence which helpful for cellulose decomposition. LEfSe analyzed three important time points (d0 vs. d7, d7 vs. d14, d49 vs. d53, LDA score > 4 and p

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

The process of skeletal muscle development is intricate and involves the regulation of a diverse array of genes. Accurate gene expression profiles are crucial for studying muscle development, making it essential to choose the right reference genes for real-time quantitative PCR (RT-qPCR). In the present study, eight candidate reference genes were identified from our previous transcriptome sequencing analysis of caprine skeletal muscle satellite cells (MuSCs), and two traditional reference genes (ACTB and GAPDH) were assessed. The quantitative levels of the candidate reference genes were determined through the RT-qPCR technique, while the stability of their expression was evaluated utilizing the GeNorm, NormFinder, BestKeeper, and RefFinder programs. Furthermore, the chosen reference genes were utilized for the normalization of the gene expression levels of PCNA and Myf5. It was determined that conventional reference genes, including ACTB and GAPDH, were not appropriate for normalizing target gene expression. Conversely, RPL14 and RPS15A, identified through RNA sequencing analysis, exhibited minimal variability and were identified as the optimal reference genes for normalizing gene expression during the proliferation and differentiation of goat MuSCs. Our research offers a validated panel of optimal reference genes for the detection of differentially expressed genes in goat muscle satellite cells using RT-qPCR.

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.

Measurement(s) genome assemble Technology Type(s) Illumina paired-end sequencing, PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies Sample Characteristic - Organism Micropterus floridanus

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

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

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.

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

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

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.

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

Early weaning and milk substitutes increase the incidence of diarrhea in young ruminants, which may modify their gut microbiota, metabolism, immunity, and health. The aim of the study was to determine if early weaning and milk substitutes affect the gut microbiota, metabolism, and immunological status of goat kids suffering from diarrhea. The 16S rRNA gene and metagenomic sequencing in feces and serum metabolomics of early-weaned and artificially reared goat kids suffering from diarrhea (DK group) and healthy goat kids reared by their mothers (HK group) were analyzed. The serum biochemistry and immunoglobulin concentration were also determined. Several probiotics, such as Streptococcus and Lactobacillus, were higher in the feces of the DK group than in feces of the HK group. Ruminococcus sp. was elevated in the feces of HKs, likely being a biomarker for goat health. Taking all the carbohydrate-active enzyme (CAZyme) families into consideration, 20 CAZyme families were different between the groups. Compared with the DK group, the relative quantity of glycoside hydrolases (GH) and glycosyltransferase (GT) families in the HK group decreased. GT70 was only identified in HK kids participating in the activity of β-glucuronosyltransferase during the carbohydrate metabolism. Overall, 24 metabolites were different between the groups, which were mainly involved in protein digestion and absorption, cyanoamino acid metabolism, and cholesterol metabolism. The concentrations of immunoglobulins G and M were significantly lower in the DK than in the HK group. In conclusion, our study characterized the fecal microbiota, metabolism, and immunological status of early-weaned and artificially reared goat kids suffering from diarrhea.

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

Circular RNAs (circRNAs) are key regulatory factors with vital functions in various biological activities. However, little has been reported concerning the genetic regulation of circRNAs during rumen development in goats. The aim of this study was to identify the genome-wide expression profiles of circRNAs in the rumen of goats during fetal development and before and after weaning. Histological morphology showed that from the fetal period (days 60 and 135 of gestation) to the prepuberal period (days 60 and 150 of age) the rumen papilla developed gradually, and the thickness of the rumen muscular layer increased. A total of 11,149 circRNAs were identified in the four development stages by RNA-sequencing. From this, 1,518 were differentially expressed circRNAs (DECs). Fifty-eight DECs were up-regulated from 60 to 135 days of gestation, and 93 from day 135 of pregnancy to 30 days after birth. A large proportion (598) of DECs were down-regulated from day 135 of gestation to 30 days after birth. The expression levels of six randomly selected circRNAs were validated by qPCR, and their back-splicing junction (BSJ) sites were also confirmed. Ontology and pathway analyses revealed that the parental genes of DECs were mainly involved in the signaling pathways related to cell proliferation and apoptosis. The interaction network of circRNAs with their target miRNAs showed its involvement in cell proliferation and apoptosis signaling pathways. In conclusion, we identified the genome-wide expression profiles of circRNAs in the rumen of goats during fetal development and before and after weaning. These results provide a basis for further study on the regulatory effect of circRNAs on the development of rumen tissues.

Abstract Background Copy number variations (CNVs) are a major form of genetic variations and are involved in animal domestication and genetic adaptation to local environments. We investigated CNVs in the domestic goat (Capra hircus) using Illumina short-read sequencing data, by comparing our lab data for 38 goats from three Chinese breeds (Chengdu Brown, Jintang Black, and Tibetan Cashmere) to public data for 26 individuals from three other breeds (two Moroccan and one Chinese) and 21samples from Bezoar ibexes. Results We obtained a total of 2394 CNV regions (CNVRs) by merging 208,649 high-confidence CNVs, which spanned ~ 267 Mb of total length and accounted for 10.80% of the goat autosomal genome. Functional analyses showed that 2322 genes overlapping with the CNVRs were significantly enriched in 57 functional GO terms and KEGG pathways, most related to the nervous system, metabolic process, and reproduction system. Clustering patterns of all 85 samples generated separately from duplications and deletions were generally consistent with the results from SNPs, agreeing with the geographical origins of these goats. Based on genome-wide F ST at each CNV locus, some genes overlapping with the highly divergent CNVs between domestic and wild goats were mainly enriched for several immunity-related pathways, whereas the genes overlapping with the highly differentiated CNVs between highland and lowland goats were mainly related to vitamin and lipid metabolism. Remarkably, a 507-bp deletion at ~ 14 kb downstream of FGF5 on chromosome 6 showed highly divergent (F ST = 0.973) between the highland and lowland goats. Together with an enhancer activity of this sequence shown previously, the function of this duplication in regulating fiber growth deserved to be further investigated in detail. Conclusion We generated a comprehensive map of CNVs in goats. Many genetically differentiated CNVs among various goat populations might be associated with the population characteristics of domestic goat breeds.

The Jianchang duck is mainly distributed in Southwest China, and has the characteristics of fast growth rate and strong abilities in lipid deposition in the liver. In order to investigate the effects of domestication process on formation of the unique characteristics of Jianchang duck, the whole genome of sixteen individuals and three pooling of Jianchang duck were re-sequenced, and genome data of 70 mallards and 83 domestic ducks from thirteen different places in China were obtained from NCBI. The population stratification and evolution analysis showed gene exchanges existed between the Jianchang and other domestic duck populations, as well as Jianchang ducks and mallards. Genomic comparison between mallards and Jianchang ducks showed genes, including CNTN1, CHRNA9, and SHANK2, which is involved in brain and nerve development, experienced strong positive selection in the process of Jianchang duck domestication. The genomic comparison between Jianchang and domestic duck populations showed that HSD17B12 and ESM1, which affect lipid metabolism, experienced strong positive selection during the domestication process. FST analysis among populations of Jianchang duck with different plumage colors indicated that MITF was related to the phenotype of a white feather, while MC1R was related to the phenotype of hemp feather. Our results provided a base for the domestication process of Jianchang duck and the genomic genes for unique traits.

Understanding the genetic composition of indigenous goats is essential to promote the scientific conservation and sustainable utilization of these breeds. The Jianchang Black (JC) goat, a Chinese native breed, is solid black and exhibits crude feed tolerance, but is characterized by a low growth rate and small body size. Based on the whole-genome sequencing data for 30 JC, 41 Jintang Black (JT), and 40 Yunshang Black (YS) goats, and 21 Bezoar ibexes, here, we investigated the genetic composition of JC goats by conducting analyses of the population structure, runs of homozygosity (ROH), genomic inbreeding, and selection signature. Our results revealed that JT and YS showed a close genetic relationship with a non-negligible amount of gene flows but were genetically distant from JC, apart from Bezoars. An average of 2039 ROHs were present in the autosomal genome per individual. The ROH-based inbreeding estimates in JC goats generally showed moderate values ranging from 0.134 to 0.264, mainly due to rapid declines in the effective population size during recent generations. The annotated genes (e.g., IL2, IL7, and KIT) overlapping with ROH islands were significantly enriched in immune-related biological processes. Further, we found 61 genes (e.g., STIM1, MYO9A, and KHDRBS2) under positive selection in JC goats via three complementary approaches, which may underly genetic adaptations to local environmental conditions. Our findings provided references for the conservation and sustainable utilization of JC goats.

Abstract Background As one of the important livestock species around the world, goats provide abundant meat, milk, and fiber to fulfill basic human needs. However, the genetic loci that underlie phenotypic variations in domestic goats are largely unknown, particularly for economically important traits. In this study, we sequenced the whole genome of 38 goats from three Chinese breeds (Chengdu Brown, Jintang Black, and Tibetan Cashmere) and downloaded the genome sequence data of 30 goats from five other breeds (four non-Chinese and one Chinese breed) and 21 Bezoar ibexes to investigate the genetic composition and selection signatures of the Chinese goat breeds after domestication. Results Based on population structure analysis and F ST values (average F ST = 0.22), the genetic composition of Chengdu Brown goats differs considerably from that of Bezoar ibexes as a result of geographic isolation. Strikingly, the genes under selection that we identified in Tibetan Cashmere goats were significantly enriched in the categories hair growth and bone and nervous system development, possibly because they are involved in adaptation to high-altitude. In particular, we found a large difference in allele frequency of one novel SNP (c.-253G>A) in the 5′-UTR of FGF5 between Cashmere goats and goat breeds with short hair. The mutation at this site introduces a start codon that results in the occurrence of a premature FGF5 protein and is likely a natural causal variant that is involved in the long hair phenotype of cashmere goats. The haplotype tagged with the AGG-allele in exon 12 of DSG3, which encodes a cell adhesion molecule that is expressed mainly in the skin, was almost fixed in Tibetan Cashmere goats, whereas this locus still segregates in the lowland goat breeds. The pigmentation gene KITLG showed a strong signature of selection in Tibetan Cashmere goats. The genes ASIP and LCORL were identified as being under positive selection in Jintang Black goats. Conclusions After domestication, geographic isolation of some goat breeds has resulted in distinct genetic structures. Furthermore, our work highlights several positively selected genes that likely contributed to breed-related traits in domestic goats.

Skeletal myogenesis is a complicated biological event that involves a succession of tightly controlled gene expressions. In order to identify novel regulators of this process, we performed mRNA-Seq studies of goat skeletal muscle satellite cells (MuSCs) cultured under proliferation (GM) and differentiation (DM1/DM5) conditions. A total of 19,871 goat genes were expressed during these stages, 198 of which represented novel transcripts. Notably, in pairwise comparisons at the different stages, 2551 differentially expressed genes (DEGs) were identified (p < 0.05), including 1560 in GM vs. DM1, 1597 in GM vs. DM5, and 959 in DM1 vs. DM5 DEGs. The time-series expression profile analysis clustered the DEGs into eight gene groups, three of which had significantly upregulated and downregulated patterns (p < 0.05). Functional enrichment analysis showed that DEGs were enriched for essential biological processes such as muscle structure development, muscle contraction, muscle cell development, striated muscle cell differentiation, and myofibril assembly, and were involved in pathways such as the MAPK, Wnt and PPAR signaling pathways. Moreover, the expression of eight DEGs (MYL2, DES, MYOG, FAP, PLK2, ADAM, WWC1, and PRDX1) was validated. These findings offer novel insights into the transcriptional regulation of skeletal myogenesis in goats.

In livestock, a comprehensive understanding of the early-life establishment and acquisition of commensal gut microbiota allow us to develop better husbandry management operations and manipulate the gut microbiota for young animals, improving the efficiency of animal production. Here, we collected 123 microbial samples of 11 healthy goat kids and their mothers to investigate the colonization and acquisition of the gut microbiota and their correlations with volatile fatty acids (VFAs) in goat kids from birth to day 56. An age-dependent increasing and more homogeneous diversity were observed for the feces of goat kids. Overall, Firmicutes, Bacteroidetes, and Proteobacteria were the predominant phyla in the fecal microbiota of goat kids, but their relative abundance varied considerably with age. Accordingly, the colonization of the fecal microbiota in goat kids was divided into three distinct stages: newborn (day 0), non-rumination stage (days 7–21), and transition stages (days 28–56). LEfSe analysis revealed a total of 49 bacterial biomarkers that are stage-specific (LDA score > 3, P < 0.05). Significant Spearman correlations (P < 0.05) were observed between the abundances of several bacterial biomarkers and the VFA concentrations. Furthermore, a substantial difference in the fecal microbiota composition was present between 56-day-old goat kids and mothers, whereas there was a moderate difference in the rumen microbiota between them. Among four body sites (i.e., feces, oral cavity, vagina, and breast milk) of mothers, the maternal vaginal and breast milk microbiota were the major source of the fecal microbiota of goat kids in the first 56 days after birth, although their contributions decreased with age and unknown sources increased after day 28. In summary, we concluded that the gut bacterial community in goat kids after birth was mainly acquired from the maternal vagina and breast milk. Its colonization showed three distinct phases with dramatic shifts of composition mainly driven by age and diet changes. Our results provide a framework for a better understanding of the roles of the gut microbiota in young ruminants.

An amendment to this paper has been published and can be accessed via the original article.

Objective It is commonly accepted that adiponectin binds to its two receptors to regulate fatty acid metabolism in adipocytes. To better understand their functions in the regulation of intramuscular adipogenesis in goats, we cloned the three genes (adiponectin [AdipoQ], adiponectin receptor 1 [AdipoR1], and AdipoR2) encoding these proteins and detected their mRNA distribution in different tissues. We also determined the role of AdipoQ in the adipogenic differentiation of goat skeletal muscle satellite cells (SMSCs). Methods SMSCs were isolated using 1 mg/mL Pronase E from the longissimus dorsi muscles of 3-day-old female Nanjiang brown goats. Adipogenic differentiation was induced in satellite cells by transferring the cells to Dulbecco’s modified Eagle’s medium supplemented with an isobutylmethylxanthine, dexamethasone and insulin cocktail. The pEGFP-N1-AD plasmid was transfected into SMSCs using Lipofectamine 2000. Expression of adiponectin in tissues and SMSCs was detected by quantitative polymerase chain reaction and immunocytochemical staining. Results The three genes were predominantly expressed in adipose and skeletal muscle tissues. According to fluorescence and immunocytochemical analyses, adiponectin protein expression was only observed in the cytoplasm, suggesting that adiponectin is localized to the cytoplasm of goat SMSCs. In SMSCs overexpressing the AdipoQ gene, adiponectin promoted SMSC differentiation into adipocytes and significantly (p

As the largest chamber of the ruminant stomach, the rumen not only serves as the principal absorptive surface and nutrient transport pathway from the lumen into the animal, but also plays an important short-chain fatty acid (SCFA) metabolic role in addition to protective functions. Accurate characterization of the gene expression profiles of genes of interest is essential to the exploration of the intrinsic regulatory mechanisms of rumen development in goats. Thus, the selection of suitable reference genes (RGs) is an important prerequisite for real-time quantitative PCR (RT-qPCR). In the present study, 16 candidate RGs were identified from our previous transcriptome sequencing of caprine rumen tissues. The quantitative expressions of the candidate RGs were measured using the RT-qPCR method, and the expression stability of the RGs was assessed using the geNorm, NormFinder, and BestKeeper programs. GeNorm analysis showed that the M values were less than 0.5 for all the RGs except GAPT4, indicating that they were stably expressed in the rumen tissues throughout development. RPS4X and RPS6 were the two most stable RGs. Furthermore, the expressions of two randomly selected target genes (IGF1 and TOP2A), normalized by the selected most stable RGs (RPS4X and RPS6), were consistent with the results of RNA sequencing, while the use of GAPDH and ACTB as RGs resulted in altered profiles. Overall, RPS4X and RPS6 showed the highest expression stability and the lowest coefficients of variation, and could be used as the optimal reference combination for quantifying gene expression in rumen tissues via RT-qPCR analysis.

Long noncoding RNAs (lncRNAs) have emerged as essential regulators of skeletal myogenesis, but few myogenesis-associated lncRNAs have been identified and our understanding of their regulatory mechanisms remains limited, particularly in goat. Here, we identified a novel lncRNA, TCONS_00006810 (named lncR-125b), from our previous lncRNA sequencing data on fetal (45, 60, and 105 days of gestation, three biological replicates for each point) and postnatal (3 days after birth, n = 3) goat skeletal muscle, and found that it is highly expressed in skeletal muscle and gradually upregulated during skeletal muscle satellite cell (SMSC) differentiation in goat. Notably, overexpression of lncR-125b accelerated the expression of myogenic differentiation 1 (MyoD 1) and myogenin (MyoG), and the formation of myotubes, and knockdown of lncR-125b showed opposite effects in SMSCs. Results of dual-luciferase assay and quantitative real-time polymerase chain reaction revealed that lncR-125b acts as a molecular sponge for miR-125b and that insulin-like growth factor 2 (IGF2), a critical regulator of skeletal myogenesis, is a direct target gene of miR-125b. Further analyses showed that lncR-125b negatively regulates miR-125b expression and positively regulates IGF2 expression in SMSCs. Mechanistically, lncR-125b promotes SMSC differentiation by functioning as a competing endogenous RNA (ceRNA) for miR-125b to control IGF2 expression. These findings identify lncR-125b as a novel noncoding regulator of muscle cell differentiation and skeletal muscle development in goat.

Adipose tissues are mainly divided into brown adipose tissue (BAT) and white adipose tissue (WAT). WAT mainly functions to buffer excess calories, whereas BAT plays a role in the non-shivering thermogenesis to maintain body temperature and energy balance. Moreover, circRNAs play important roles in various biological processes. However, knowledge of the expression profile and function of circRNAs from BAT to WAT remains largely unknown. In this study, a total of 6610 unique circRNAs were identified in the perirenal adipose tissues of 1-day, 30-days, and 1-year goats. Functional annotation revealed that host genes of circRNAs were involved in some BAT-related pathways, such as the thyroid hormone signaling pathway, MAPK signaling pathway, and VEGF signaling pathway. Furthermore, a total of 61 DEcircRNAs were detected across three stages. Additionally, five selected circRNAs were validated by RNase R assay, qPCR, and Sanger sequencing. Finally, the circRNA–miRNA network was constructed between the DEcircRNAs and their miRNA binding sites.

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

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

In this study, we conducted comparative analyses to characterize the rumen microbiota and volatile fatty acid (VFA) profiles of weaned Nanjiang Yellow goat kids under shrub-grassland grazing (GR), shrub-grassland grazing and supplementary feeding (SF), and indoor feeding (IF) systems. We observed significant differences (p < 0.05) in the concentrations of total VFA and the proportions of acetate and butyrate in the rumen fluid among the three groups, whereas the proportions of propionate and the acetate/propionate ratio did not differ substantially. Alpha diversity of the rumen bacterial and archaeal populations in the GR and SF kids was significantly higher (p < 0.05) than that in the IF goat kids, and significant differences (p < 0.05) in similarity were observed in the comparisons of GR vs. IF and SF vs. IF. The most predominant bacterial phyla were Bacteroidetes and Firmicutes across the three groups, and the archaeal community was mainly composed of Euryarchaeota. At the genus and species levels, the cellulose-degrading bacteria, including Lachnospiraceae, Ruminococcaceae and Butyrivibrio fibrisolvens, were abundant in the GR and SF groups. Furthermore, 27 bacterial and 11 unique archaeal taxa, such as Lachnospiraceae, Butyrivibrio fibrisolvens, and Methanobrevibacter ruminantium, were identified as biomarkers, and showed significantly different (p < 0.05) abundances among the three groups. Significant Spearman correlations (p < 0.05), between the abundances of several microbial biomarkers and the concentrations of VFAs, were further observed. In summary, our results demonstrated that the adaptation to grazing required more rumen bacterial populations due to complex forage types in shrub-grassland, although the rumen fermentation pattern did not change substantially among the three feeding systems. Some microbial taxa could be used as biomarkers for different feeding systems, particularly cellulose-degrading bacteria associated with grazing.

The rumen is an important digestive organ in ruminants. Numerous regulatory factors including microRNAs (miRNAs) are involved in embryonic organ development. In the present study, miRNAs expressed in the rumens of goats (Capra hircus) and their potential roles in the pathways involved in rumen development were identified using high-throughput sequencing. Histological morphology revealed a distinct difference in each layer of rumen during the period from embryonic day 60 (E60) to embryonic day 135 (E135). We determined the expression profiles of miRNAs in the goat rumen, and identified 423 known miRNAs and 559 potentially novel miRNAs in the E60 and E135 embryonic rumen, respectively. Bioinformatics analysis annotated the 42 differentially expressed miRNAs and the top 10 most highly expressed miRNAs of the two libraries to 48 and 38 gene ontology categories, as well as to 168 and 71 Kyoto Encyclopedia of Genes and Genomes pathways, respectively. The expression patterns of eight randomly selected miRNAs were validated by stem-loop quantitative reverse transcription PCR, suggesting that the sequencing data were reliable. We profiled the genome-wide expression of rumen-expressed miRNAs at different prenatal stages of rumen tissues, revealing that a subset of miRNAs might play important roles in the formation of the rumen layers. Taken together, these findings will aid the investigation of dominant rumen-related miRNA sets and help understand the genetic control of rumen development in goats.

Ectopic lipid accumulation in muscle is important not only for obesity and myopathy treatment, but also for meat quality improvement in farm animals. However, the molecular mechanisms involved in lipid metabolism in muscle satellite cells are still elusive. In this study, SB216763 reduced GSK3β activation by increasing the level of pGSK3β (Ser9) and decreasing the level of total GSK3β protein. GSK3β inhibition decreased lipid accumulation and downregulated the expression level of lipogenesis-related genes in the adipogenic differentiation of goat muscle satellite cells. Furthermore, SB216763 treatment increased the levels of pAMPKα (T172) and pACC (Ser79). Further, we found that GSK3β inhibition promoted levels of LC3B-II and reduced the protein levels of p62 to induce the autophagy in muscle satellite cells. Taken together, our results provide new insight into a critical function for GSK3β: modulating lipid accumulation in goat muscle satellite cells through activating the AMPK pathway.

Long noncoding RNAs (lncRNAs) play an important role in the thermogenesis and energy storage of brown adipose tissue (BAT). However, knowledge of the cellular transition from BAT to white adipose tissue (WAT) and the potential role of lncRNAs in goat adipose tissue remains largely unknown. In this study, we analyzed the transformation from BAT to WAT using histological and uncoupling protein 1 (UCP1) gene analyses. Brown adipose tissue mainly existed within the goat perirenal fat at 1 day and there was obviously a transition from BAT to WAT from 1 day to 1 year. The RNA libraries constructed from the perirenal adipose tissues of 1 day, 30 days, and 1 year goats were sequenced. A total number of 21,232 lncRNAs from perirenal fat were identified, including 5393 intronic-lncRNAs and 3546 antisense-lncRNAs. Furthermore, a total of 548 differentially expressed lncRNAs were detected across three stages (fold change ≥ 2.0, false discovery rate (FDR) < 0.05), and six lncRNAs were validated by qPCR. Furthermore, trans analysis found lncRNAs that were transcribed close to 890 protein-coding genes. Additionally, a coexpression network suggested that 4519 lncRNAs and 5212 mRNAs were potentially in trans-regulatory relationships (r > 0.95 or r < −0.95). In addition, Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that the targeted genes were involved in the biosynthesis of unsaturated fatty acids, fatty acid elongation and metabolism, the citrate cycle, oxidative phosphorylation, the mitochondrial respiratory chain complex, and AMP-activated protein kinase (AMPK) signaling pathways. The present study provides a comprehensive catalog of lncRNAs involved in the transformation from BAT to WAT and provides insight into understanding the role of lncRNAs in goat brown adipogenesis.

Skeletal muscle development is a complex biological process regulated by numerous genes and non-coding RNAs, such as microRNAs (miRNAs). In the current study, we made use of the deep sequencing data from Jianzhou Da'er goat longissimus dorsi sampled on days 45, 60, and 105 of gestation, as well as day three after birth to identify miRNAs that regulate goat skeletal myogenesis, and examine their temporal expression profiles. A total of 410 known goat miRNAs, 752 miRNA homologs and 88 novel miRNAs were identified across four stages. Besides three myomiRs, the abundance of 17 miRNAs, including chi-miR-424, chi-miR-542-3p and chi-miR-136-5p was more than 10,000 reads per million mapped reads (RPM), on average. Furthermore, 50 miRNAs with more than 100 RPM clustered at the imprinted DLK1-DIO3 locus on chromosome 21 and showed similar expression patterns, indicating that these miRNAs played important roles in skeletal myogenesis of goats. Based on pairwise comparisons, 221 differentially expressed (DE), known miRNAs were identified across four stages. GO and KEGG analyses of the genes targeted by the DE miRNAs revealed the significantly enriched processes and pathways to be consistent with temporal changes of skeletal muscle development across all sampled stages. However, follow-up experimental studies were required to explore functions of these miRNAs and targets underlying skeletal myogenesis.