Haplotype-resolved 3D chromatin architecture of the hybrid pig.

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Comparative three-dimensional genome architectures of adipose tissues provide insight into human-specific regulation of metabolic homeostasis.

Elucidating the regulatory mechanisms of human adipose tissues (ATs) evolution is essential for understanding human-specific metabolic regulation, but the functional importance and evolutionary dynamics of three-dimensional (3D) genome organizations of ATs are not well defined. Here, we compared the 3D genome architectures of anatomically distinct ATs from humans and six representative mammalian models. We recognized evolutionarily conserved and human-specific chromatin conformation in ATs at multiple scales, including compartmentalization, topologically associating domain (TAD), and promoter-enhancer interactions (PEI), which have not been described previously. We found PEI are much more evolutionarily dynamic with respect to compartmentalization and topologically associating domain. Compared to conserved PEIs, human-specific PEIs are enriched for human-specific sequence, and the binding motifs of their potential mediators (transcription factors) are less conserved. Our data also demonstrated that genes involved in the evolutionary dynamics of chromatin organization have weaker transcriptional conservation than those associated with conserved chromatin organization. Furthermore, the genes involved in energy metabolism and the maintenance of metabolic homeostasis are enriched in human-specific chromatin organization, while housekeeping genes, health-related genes, and genetic variations are enriched in evolutionarily conserved compared to human-specific chromatin organization. Finally, we showed extensively divergent human-specific 3D genome organizations among one subcutaneous and three visceral ATs. Together, these findings provide a global overview of 3D genome architecture dynamics between ATs from human and mammalian models and new insights into understanding the regulatory evolution of human ATs.

Acetate Alleviates Gut Microbiota Depletion-Induced Retardation of Skeletal Muscle Growth and Development in Young Mice.

The normal growth and development of skeletal muscle is essential for the health of the body. The regulation of skeletal muscle by intestinal microorganisms and their metabolites has been continuously demonstrated. Acetate is the predominant short-chain fatty acids synthesized by gut microbiota through the fermentation of dietary fiber; however, the underlying molecular mechanisms governing the interaction between acetate and skeletal muscle during the rapid growth stage remains to be further elucidated. Herein, specific pathogen-free (SPF) mice, germ-free (GF) mice, and germ-free mice supplemented with sodium acetate (GS) were used to evaluate the effects of acetate on the skeletal muscle growth and development of young mice with gut microbiota deficiency. We found that the concentration of serum acetate, body mass gain, succinate dehydrogenase activity, and expression of the myogenesis maker gene of skeletal muscle in the GS group were higher than those in the GF group, following sodium acetate supplementation. Furthermore, the transcriptome analysis revealed that acetate activated the biological processes that regulate skeletal muscle growth and development in the GF group, which are otherwise inhibited due to a gut microbiota deficiency. The in vitro experiment showed that acetate up-regulated Gm16062 to promote skeletal muscle cell differentiation. Overall, our findings proved that acetate promotes skeletal muscle growth and development in young mice via increasing Gm16062 expression.

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation.

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene's promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

Comparative 3D genome architecture in vertebrates.

BACKGROUND: The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species. RESULTS: This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture. CONCLUSIONS: Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle.

BACKGROUND: Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined. RESULTS: In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers. CONCLUSIONS: This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Comprehensive variation discovery and recovery of missing sequence in the pig genome using multiple de novo assemblies.

Uncovering genetic variation through resequencing is limited by the fact that only sequences with similarity to the reference genome are examined. Reference genomes are often incomplete and cannot represent the full range of genetic diversity as a result of geographical divergence and independent demographic events. To more comprehensively characterize genetic variation of pigs (Sus scrofa), we generated de novo assemblies of nine geographically and phenotypically representative pigs from Eurasia. By comparing them to the reference pig assembly, we uncovered a substantial number of novel SNPs and structural variants, as well as 137.02-Mb sequences harboring 1737 protein-coding genes that were absent in the reference assembly, revealing variants left by selection. Our results illustrate the power of whole-genome de novo sequencing relative to resequencing and provide valuable genetic resources that enable effective use of pigs in both agricultural production and biomedical research.

Interactive effect of thermal and hypoxia on largemouth bass (Micropterus salmoides) gill and liver: Aggravation of oxidative stress, inhibition of immunity and promotion of cell apoptosis.

High temperatures and low oxygen in aquatic environments, such as intensive aquaculture or in natural watersheds, inevitably cause stress in fish. Fish are exposed to high temperatures during the summer, which exacerbates hypoxia. Hypoxia (1.2 ± 0.2 mg/L) under 20 °C (20 HG) and 26 °C (26 HG) was simulated to induce stress in largemouth bass (Micropterus salmoides). Related enzymes and genes involved in antioxidant, immune, and apoptotic responses were selected to explore the interactive effects of temperature and hypoxia on largemouth bass. The results showed that malondialdehyde (MDA) levels in plasma, gill, and liver increased in the 26 HG (p < 0.05). Liver superoxide dismutase (SOD) activity increased in the 26 HG. Peak SOD (SOD1, SOD2, SOD3a, and SOD3b), CAT, and GSH-Px mRNA levels in the gill and liver were observed at 12-24 h of stress. The levels of gill and liver total antioxidant capacity, catalase (CAT), glutathione peroxidase (GSH-Px) activities and other enzyme activities and genes in the 26 HG were higher than those in the 20 HG (p < 0.05). The gill and liver acid phosphatase and alkaline phosphatase activities increased with time in the 26 HG (p < 0.05), while gill and liver lysozyme activities in the 26 HG were lower than those in the 20 HG (p < 0.05). Tumor necrosis factor-α mRNA level was upregulated in the gill and downregulated in the liver at 24 h in the 26 HG. Interleukin (IL)-1ß and IL-8 mRNA levels were upregulated in the gill and liver in the 26 HG at 24 h, whereas IL-15 mRNA level was downregulated in the 26 HG at 12 h. Transforming growth factor-ß1 mRNA level was upregulated in the gill in the 20 HG at 24 h, but downregulated in gill and liver in the 26 HG at 24 h. Similarly, IL-10, Hepcidin-1, and Hepcidin-2 showed lower expression levels in the 26 HG. Gill and liver caspase-3 activities were higher in the 26 HG (p < 0.05), and gill caspase-3 activity was higher than that in the liver. The mRNA levels of proapoptotic genes (caspase-3, caspase-8, and caspase-9) were higher in the 26 HG. The present study demonstrates the interactive effects of temperature and hypoxia on stress in largemouth bass gill and liver. These results will be helpful to understand the mechanisms of stress induced by temperature and hypoxia in fish and provide a theoretical basis for aquaculture management.

SNPs associated with body weight and backfat thickness in two pig breeds identified by a genome-wide association study.

Growth and fat deposition are important economic traits due to the influence on production in pigs. In this study, a dataset of 1200 pigs with 345,570 SNPs genotyped by sequencing (GBS) was used to conduct a GWAS with single-marker regression method to identify SNPs associated with body weight and backfat thickness (BFT) and to search for candidate genes in Landrace and Yorkshire pigs. A total of 27 and 13 significant SNPs were associated with body weight and BFT, respectively. In the region of 149.85-149.89 Mb on SSC6, the SNP (SSC6: 149876737) for body weight and the SNP (SSC6: 149876507) for BFT were in the same locus region (a gap of 230 bp). Two SNPs were located in the DOCK7 gene, which is a protein-coding gene that plays an important role in pigmentation. Two SNPs located on SSC8: 54567459 and SSC11: 33043081 were found to overlap weight and BFT; however, no candidate gene was found in these regions. In addition, based on other significant SNPs, two positional candidate genes, NSRP1 and CADPS, were proposed to influence weight. In conclusion, this is the first study report using GBS data to identify the significant SNPs for weight and BFT. A total of four particularly interesting SNPs and one potential candidate genes (DOCK7) were found for these traits in domestic pigs. This study improves our knowledge to better understand the complex genetic architecture of weight and BFT, but further validation studies of these candidate loci and genes are recommended in pigs.

Integrated analysis of RNA-seq and microRNA-seq depicts miRNA-mRNA networks involved in stripe patterns of Botia superciliaris skin.

Botia superciliaris, an endemic cobitid fish in China, is widely accepted by Chinese consumers because its edibility. Recently, the black and yellow stripes of B. superciliaris skin have made this species increasingly popular as a novel ornamental fish. However, the genetic basis of the stripe patterns in B. superciliaris skin has not been extensively studied. In this study, Illumina sequencing was employed to identify the mRNAs and miRNAs involved in stripe pattern formation in B. superciliaris skin. A total of 147.25 and 155.15 million (M) high-quality transcriptome reads were generated from three black and yellow skin libraries respectively, which resulted in 159,327 unigenes that were used as reference sequences. A total of 3169 genes exhibited significantly differential expression patterns (fold-change ≥ 2 or ≤ 0.5 and q ≤ 0.05), including 1891 upregulated genes (59.67%) and 1278 downregulated genes (40.33%) in black vs yellow skin. These genes were enriched in 50 GO terms and 10 KEGG pathways (q ≤ 0.05), including melanogenesis, with 21 upregulated genes and 5 downregulated genes in black vs yellow skin. Based on the zebrafish genome, miRNA-seq identified a total of 355 miRNAs, which included 38 novel miRNAs. Furthermore, 87 differentially expressed miRNAs including 50 upregulated and 37 downregulated miRNAs were identified in different color skin (fold-change ≥ 2 or ≤ 0.5 and q ≤ 0.05). Then, target prediction revealed a variety of putative target genes; differentially expressed mRNAs and miRNAs patterns of high-throughput sequencing were validated in 5 mRNAs and miR-217-5p by qRT-PCR. In vivo tests and dual-luciferase reporter assay revealed that overexpression of miR-217-5p can inhibit pheomelanin formation by targeting Zgc. In this study, a comparative analysis was conducted to profile the transcriptome of black and yellow skin for B. superciliaris, and we detected key genes and important miRNAs involved in the B. superciliaris skin pigmentation process. These results will enhance understanding of molecular mechanisms underlying skin pigmentation and facilitate molecular-assisted selection of highly valued skin colors.

Single step genome-wide association studies based on genotyping by sequence data reveals novel loci for the litter traits of domestic pigs.

In this study, data genotyping by sequence (GBS) was used to perform single step GWAS (ssGWAS) to identify SNPs associated with the litter traits in domestic pigs and search for candidate genes in the region of significant SNPs. After quality control, 167,355 high-quality SNPs from 532 pigs were obtained. Phenotypic traits on 2112 gilt litters from 532 pigs were recorded including total number born (TNB), number born alive (NBA), and litter weight born alive (LWB). A single-step genomic BLUP approach (ssGBLUP) was used to implement the genome-wide association analysis at a 5% genome-wide significance level. A total of 8, 23 and 20 significant SNPs were associated with TNB, NBA, and LWB, respectively, and these significant SNPs accounted for 62.78%, 79.75%, and 58.79% of genetic variance. Furthermore, 1 (SSC14: 16314857), 4 (SSC1: 81986236, SSC1: 66599775, SSC1: 161999013, and SSC1: 267883107), and 5 (SSC9: 29030061, SSC2: 32368561, SSC5: 110375350, SSC13: 45619882 and SSC13: 45647829) significant SNPs for TNB, NBA, and LWB were inferred to be novel loci. At SSC1, the AIM1 and FOXO3 genes were found to be associated with NBA; these genes increase ovarian reproductive capacity and follicle number and decrease gonadotropin levels. The genes SLC36A4 and INTU are involved in cell growth, cytogenesis and development were found to be associated with LWB. These significant SNPs can be used as an indication for regions in the Sus scrofa genome for variability in litter traits, but further studies are expected to confirm causative mutations.

miR-27a-5p Attenuates Hypoxia-induced Rat Cardiomyocyte Injury by Inhibiting Atg7.

Acute myocardial infarction (AMI) is an ischemic heart disease with high mortality worldwide. AMI triggers a hypoxic microenvironment and induces extensive myocardial injury, including autophagy and apoptosis. MiRNAs, which are a class of posttranscriptional regulators, have been shown to be involved in the development of ischemic heart diseases. We have previously reported that hypoxia significantly alters the miRNA transcriptome in rat cardiomyoblast cells (H9c2), including miR-27a-5p. In the present study, we further investigated the potential function of miR-27a-5p in the cardiomyocyte response to hypoxia, and showed that miR-27a-5p expression was downregulated in the H9c2 cells at different hypoxia-exposed timepoints and the myocardium of a rat AMI model. Follow-up experiments revealed that miR-27a-5p attenuated hypoxia-induced cardiomyocyte injury by regulating autophagy and apoptosis via Atg7, which partly elucidated the anti-hypoxic injury effects of miR-27a-5p. Taken together, this study shows that miR-27a-5p has a cardioprotective effect on hypoxia-induced H9c2 cell injury, suggesting it may be a novel target for the treatment of hypoxia-related heart diseases.

Identifying SNPs and candidate genes for three litter traits using single-step GWAS across six parities in Landrace and Large White pigs.

Total number born (TNB), number born alive (NBA), and litter weight born alive (LWB) are critically important traits in pig production. The sow's parity is one of the major factors influencing litter traits. Because of monogenic or polygenic contributions and the presence of temporal gene effects in different sows' parities, it is difficult to clarify the biological and genetic background. To systematically explore the genetic mechanism of litter traits, we conducted 18 GWASs using single-step GWAS (ssGWAS) based on two breeds (908 Landrace and 1,130 Large White sow litter records) for each litter trait in different parities. A total of 300 Landrace and 300 Large White sows were genotyped by sequencing (GBS). ssGWAS was performed separately for each breed and each parity due to population stratification and temporal gene effect. In summary, we identified 80 (15 for Landrace and 65 for Large White), 227 (52 for Landrace, 175 for Large White), and 187 (34 for Landrace, 153 for Large White) single nucleotide polymorphisms (SNPs) affecting TNB, NBA, and LWB, respectively. Of them, we suggest that a total of 22 loci (SSC1: 125098202, SSC1: 117560058, SSC14: 147794697, SSC8: 84823302, SSC9: 143554876, and SSC9: 138766097 for Landrace; SSC1: 4023577, SSC1: 3859573, SSC1: 4891063, SSC16: 5197665, SSC10: 32050819, SSC13: 13552924, SSC13: 92819, SSC17: 3579607, SSC13: 196698221, SSC7: 30918403, SSC16: 46221484, SSC16: 46169204, SSC2: 41988642, SSC2: 44475457, SSC2: 42521875, and SSC7: 58411951 for Large White) are shared by TNB, NBA, and LWB. These results indicate the existence of gene temporal effect in each parity. Furthermore, our findings suggest four interesting candidate genes (FBXL7, ALDH1A2, LEPR, and DDX1) associated with litter traits in different parities that have a major effect on embryonic development progression. In conclusion, 22 crucial SNPs and four interesting candidate genes were identified for three litter traits across six parities. These findings advance our understanding of the genetic architecture of litter traits and confirm the presence of temporal gene effects in different parities. Importantly, functional validation studies for findings of particular interest are recommended in litter traits.

Exploration of exosomal microRNA expression profiles in pigeon 'Milk' during the lactation period.

BACKGROUND: Pigeon crop has the unique ability to produce a nutrient rich substance termed pigeon 'milk' (PM), which has functional resemblance with the mammalian milk. Previous researches have demonstrated that a large number of exosomes and exosomal miRNAs exist in mammalian milk, and many of them are associated with immunity, growth and development. However, to date, little is known about the exosomes and exosomal miRNAs in PM. RESULTS: In this study, we isolated the exosomes from PM and used small RNA sequencing to investigate the distribution and expression profiles of exosomal miRNAs. A total of 301 mature miRNAs including 248 conserved and 53 novel miRNAs were identified in five lactation stages i.e. 1d, 5d, 10d, 15d, and 20d. From these, four top 10 conserved miRNAs (cli-miR-21-5p, cli-miR-148a-3p, cli-miR-10a-5p and cli-miR-26a-5p) were co-expressed in all five stages. We speculate that these miRNAs may have important role in the biosynthesis and metabolism of PM. Moreover, similar to the mammalian milk, a significant proportion of immune and growth-related miRNAs were also present and enriched in PM exosomes. Furthermore, we also identified 41 orthologous miRNAs group (giving rise to 81 mature miRNA) commonly shared with PM, human, bovine and porcine breast milk. Additionally, functional enrichment analysis revealed the role of exosomal miRNAs in organ development and in growth-related pathways including the MAPK, Wnt and insulin pathways. CONCLUSIONS: To sum-up, this comprehensive analysis will contribute to a better understanding of the underlying functions and regulatory mechanisms of PM in squabs.

MicroRNA-125a-5p Affects Adipocytes Proliferation, Differentiation and Fatty Acid Composition of Porcine Intramuscular Fat.

Intramuscular fat (IMF) content and composition are considered crucial indicators of porcine meat quality. However, the molecular mechanism of porcine IMF development is still mostly unclear. Recently, new evidence suggested that microRNA (miRNAs) play important roles in porcine intramuscular adipogenesis. Previously, microRNA-125a-5p (miR-125a-5p) was identified as an important regulator of adipogenesis. In the present study, we found that the expression of miR-125a-5p is dynamically regulated during porcine intramuscular preadipocytes differentiation and that its expression levels in different porcine muscle tissues were negatively involved with IMF content. To investigate the potential function role of miR-125a-5p in IMF development, porcine intramuscular preadipocytes were collected and transfected with miR-125a-5p mimics, inhibitors, or a negative control (NC), respectively. The results showed that overexpression of miR-125a-5p promoted proliferation and inhibited differentiation of porcine intramuscular preadipocytes while inhibition of miR-125a-5p had the opposite effects. Furthermore, a luciferase reporter assay demonstrated that porcine kruppel like factor 3 (KLF13) is a target gene of miR-125a-5p during porcine intramuscular preadipocytes differentiation. Interestingly, porcine ELOVL fatty acid elongase 6 (ELOVL6), a regulator of fatty acid composition, was also identified as a target gene of miR-125a-5p during porcine intramuscular adipogenesis. Further studies show that miR-125a-5p overexpression reduced total saturated fatty acids (SFA) content and monounsaturated fatty acids (MUFA)/SFA ratios while having no significant impact on polyunsaturated fatty acids (PUFA)/SFA and n-6/n-3 ratios. Taken together, our results identified that miR-125a-5p may be a novel regulator of porcine intramuscular adipogenesis and the fatty acid composition of porcine IMF.

Guanidinoacetic Acid Regulates Myogenic Differentiation and Muscle Growth Through miR-133a-3p and miR-1a-3p Co-mediated Akt/mTOR/S6K Signaling Pathway.

Guanidinoacetic acid (GAA), an amino acid derivative that is endogenous to animal tissues including muscle and nerve, has been reported to enhance muscular performance. MicroRNA (miRNA) is a post-transcriptional regulator that plays a key role in nutrient-mediated myogenesis. However, the effects of GAA on myogenic differentiation and skeletal muscle growth, and the potential regulatory mechanisms of miRNA in these processes have not been elucidated. In this study, we investigated the effects of GAA on proliferation, differentiation, and growth in C2C12 cells and mice. The results showed that GAA markedly inhibited the proliferation of myoblasts, along with the down-regulation of cyclin D1 (CCND1) and cyclin dependent kinase 4 (CDK4) mRNA expression, and the upregulation of cyclin dependent kinase inhibitor 1A (P21) mRNA expression. We also demonstrated that GAA treatment stimulated myogenic differentiation 1 (MyoD) and myogenin (MyoG) mRNA expression, resulting in an increase in the myotube fusion rate. Meanwhile, GAA supplementation promoted myotube growth through increase in total myosin heavy chain (MyHC) protein level, myotubes thickness and gastrocnemius muscle cross-sectional area. Furthermore, small RNA sequencing revealed that a total of eight miRNAs, including miR-133a-3p and miR-1a-3p cluster, showed differential expression after GAA supplementation. To further study the function of miR-133a-3p and miR-1a-3p in GAA-induced skeletal muscle growth, we transfected miR-133a-3p and miR-1a-3p mimics into myotube, which also induced muscle growth. Through bioinformatics and a dual-luciferase reporter system, the target genes of miR-133a-3p and miR-1a-3p were determined. These two miRNAs were shown to modulate the Akt/mTOR/S6K signaling pathway by restraining target gene expression. Taken together, these findings suggest that GAA supplementation can promote myoblast differentiation and skeletal muscle growth through miR-133a-3p- and miR-1a-3p-induced activation of the AKT/mTOR/S6K signaling pathway.

Testosterone-Dependent miR-26a-5p and let-7g-5p Act as Signaling Mediators to Regulate Sperm Apoptosis via Targeting PTEN and PMAIP1.

Recent evidence suggests that testosterone deficiency can dramatically decrease the quality of sperm. MicroRNAs (miRNAs) are conserved mediators of post-transcriptional gene regulation in eukaryotes. However, the systemic regulation and function of miRNAs in sperm quality decline induced by testosterone deficiency has not been investigated. Here, we found that the sperm apoptosis was significantly enhanced and the sperm motility was dramatically decreased in hemicastrated pigs. We then used small RNA sequencing to detect miRNA profiles of sperm from pigs with prepubertal hemicastration (HC) and compared them with control libraries. We identified 16 differentially expressed (DE) miRNAs between the sperm of prepubertal HC and control (CT) pigs. Functional enrichment analysis indicated that the target genes of these DE miRNAs were mainly enriched in apoptosis-related pathways including the p53, mitogen-activated protein kinase (MAPK), and mammalian target of rapamycin (mTOR) pathways. Furthermore, gain- and loss-of-function analyses demonstrated potential anti-apoptotic effects of the DE miRNAs miR-26a-5p and let-7g-5p on sperm cells. The luciferase reporter assay confirmed that PTEN and PMAIP1 are targets of miR-26a-5p and let-7g-5p, respectively. Spearman’s correlation analysis revealed significantly positive correlations between the sperm and its corresponding seminal plasma exosomes regarding the miRNA expression levels. In conclusion, testosterone deficiency-induced changes in the miRNA components of seminal plasma exosomes secreted by the genital tract may partially elucidate sperm miRNAome alterations, which are further responsible for the decline of sperm motility.

Hemicastration induced spermatogenesis-related DNA methylation and gene expression changes in mice testis.

OBJECTIVE: Hemicastration is a unilateral orchiectomy to remove an injured testis, which can induce hormonal changes and compensatory hypertrophy of the remaining testis, and may influence spermatogenesis. However, the underlying molecular mechanisms are poorly understood. Here, we investigated the impact of hemicastration on remaining testicular function. METHODS: Prepubertal mice (age 24 days) were hemicastrated, and their growth was monitored until they reached physical maturity (age 72 days). Subsequently, we determined testis DNA methylation patterns using reduced representation bisulfite sequencing of normal and hemicastrated mice. Moreover, we profiled the testicular gene expression patterns by RNA sequencing (RNA-seq) to examine whether methylation changes affected gene expression in hemicastrated mice. RESULTS: Hemicastration did not significantly affect growth or testosterone (p>0.05) compared with control. The genome-wide DNA methylation pattern of remaining testis suggested that substantial genes harbored differentially methylated regions (1,139) in gene bodies, which were enriched in process of protein binding and cell adhesion. Moreover, RNA-seq results indicated that 46 differentially expressed genes (DEGs) involved in meiotic cell cycle, synaptonemal complex assembly and spermatogenesis were upregulated in the hemicastration group, while 197 DEGs were downregulated, which were related to arachidonic acid metabolism. Integrative analysis revealed that proteasome 26S subunit ATPase 3 interacting protein gene, which encodes a protein crucial for homologous recombination in spermatocytes, exhibited promoter hypomethylation and higher expression level in hemicastrated mice. CONCLUSION: Global profiling of DNA methylation and gene expression demonstrated that hemicastration-induced compensatory response maintained normal growth and testicular morphological structure in mice.

Transcriptome-wide N 6 -methyladenosine methylome profiling of porcine muscle and adipose tissues reveals a potential mechanism for transcriptional regulation and differential methylation pattern.

BACKGROUND: N 6 -methyladenosine (m6A) is the most prevalent internal form of modification in messenger RNA in higher eukaryotes and potential regulatory functions of reversible m6A methylation on mRNA have been revealed by mapping of m6A methylomes in several species. m6A modification in active gene regulation manifests itself as altered methylation profiles in a tissue-specific manner or in response to changing cellular or species living environment. However, up to date, there has no data on m6A porcine transcriptome-wide map and its potential biological roles in adipose deposition and muscle growth. METHODS: In this work, we used methylated RNA immunoprecipitation with next-generation sequencing (MeRIP-Seq) technique to acquire the first ever m6A porcine transcriptome-wide map. Transcriptomes of muscle and adipose tissues from three different pig breeds, the wild boar, Landrace, and Rongchang pig, were used to generate these maps. RESULTS: Our findings show that there were 5,872 and 2,826 m6A peaks respectively, in the porcine muscle and adipose tissue transcriptomes. Stop codons, 3'-untranslated regions, and coding regions were found to be mainly enriched for m6A peaks. Gene ontology analysis revealed that common m6A peaks in nuclear genes are associated with transcriptional factors, suggestive of a relationship between m6A mRNA methylation and nuclear genome transcription. Some genes showed tissue- and breed-differential methylation, and have novel biological functions. We also found a relationship between the m6A methylation extent and the transcript level, suggesting a regulatory role for m6A in gene expression. CONCLUSION: This comprehensive map provides a solid basis for the determination of potential functional roles for RNA m6A modification in adipose deposition and muscle growth.

miR-199a-3p affects adipocytes differentiation and fatty acid composition through targeting SCD.

Body fat mass is closely associated to diseases related to obesity. MicroRNAs (miRNAs, miR) are important regulatory molecules that function as post-transcriptional gene regulators of adipocyte development. In the current study, we revealed that reduced expression of miR-199a-3p in adipose tissue resulting from high fat diet (HFD)-induced obesity in mice. Overexpression of miR-199a-3p promoted adipocyte proliferation by regulating the expression of regulating factors of the cell cycle. Furthermore, miR-199a-3p blunted lipid accumulation in 3T3-L1 adipocytes. This was accompanied by a marked decrease in the expression of adipocyte-specific genes involved in lipogenic transcription, fatty acid synthesis, and fatty acid transportation. Furthermore, the fatty acid oxidation process was enhanced. Luciferase activity assays confirmed that miR-199a-3p regulates adipocyte differentiation by directly targeting the 3'-untranslated region (3'-UTR) of stearoyl-CoA desaturase (SCD). Moreover, miR-199a-3p regulates fatty acid composition by decreasing the ratio of unsaturated fatty acids (UFAs) in adipocytes transfected with miR-199a-3p mimics. These results suggest that miR-199a-3p may promote adipocyte proliferation, while also repressing adipocyte differentiation by down-regulating SCD and changing fatty acid composition during adipogenesis.