Multi-Omics Analysis of Key microRNA-mRNA Metabolic Regulatory Networks in Skeletal Muscle of Obese Rabbits

Transcriptome analysis identifies potential molecular mechanisms for growth of fast muscle in Chinese perch juvenile

Chinese perch (Siniperca chuatsi) is an important commercial fish species in China. Understanding the molecular mechanisms of growth and development of skeletal muscle is helpful for selection breeding and improving the growth rate of Chinese perch. We analyzed histological and transcriptomic differences in fast muscle of Chinese perch between 30 days post hatching (dph) and 60 dph using histological sections and high-throughput RNA-Seq. The results showed that the diameter of muscle fibers in 30 dph Chinese perch was mainly distributed in range of 30 - 40 μm, and that of 60 dph was primarily in the range of 40 - 50 μm. 34 differentially expressed genes (DEGs) were identified in the fast muscle of Chinese perch between 30 and 60 dph, of which 9 were up-regulated and 25 were down-regulated (60 dph vs 30 dph). The DEGs, including MYH4, ENO3, Bag3, krt13 and krt18, are associated with muscle cell differentiation and fusion in Chinese perch. The analysis of the protein-protein interaction network of DEGs revealed that FOS, junb and EGR1 may involve in the development of fast muscle. KEGG enrichment results showed that the up-regulated genes in the 60 dph were associated with several pathways related to metabolism and protein synthesis, such as glycosphingolipid biosynthesis and aminoacyl-tRNA biosynthesis. The results suggest that the development of fast muscle in Chinese perch from 30 to 60 dph is accompanied by an increase in muscle fiber diameter and changes in gene expression related to muscle cell differentiation and protein synthesis. Abbreviations: DEGs, differentially expressed genes; dph, days post hatching; FC, fold change; FDR, False Discovery Rate; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; MyHCs, Myosin heavy chains; qRT-PCR, quantitative real-time PCR.

miRNAs in pancreatic cancer progression and metastasis

Small non-coding RNA or microRNA (miRNA) are critical regulators of eukaryotic cells. Dysregulation of miRNA expression and function has been linked to a variety of diseases including cancer. They play a complex role in cancers, having both tumour suppressor and promoter properties. In addition, a single miRNA can be involved in regulating several mRNAs or many miRNAs can regulate a single mRNA, therefore assessing these roles is essential to a better understanding in cancer initiation and development. Pancreatic cancer is a leading cause of cancer death worldwide, in part due to the lack of diagnostic tools and limited treatment options. The most common form of pancreatic cancer, pancreatic ductal adenocarcinoma (PDAC), is characterised by major genetic mutations that drive cancer initiation and progression. The regulation or interaction of miRNAs with these cancer driving mutations suggests a strong link between the two. Understanding this link between miRNA and PDAC progression may give rise to novel treatments or diagnostic tools. This review summarises the role of miRNAs in PDAC, the downstream signalling pathways that they play a role in, how these are being used and studied as therapeutic targets as well as prognostic/diagnostic tools to improve the clinical outcome of PDAC.

Mechanic study based on untargeted metabolomics of Pi-pa-run-fei-tang on pepper combined with ammonia induced chronic cough model mice

ETHNOPHARMACOLOGICAL RELEVANCE

Pi-pa-run-fei-tang (PPRFT), a traditional Chinese medicine formula with long-standing history, demonstrated beneficial effect on chronic cough. However, the mechanism underlying efficacy unclear. In current research, we explored the impact and molecular mechanism of chronic cough mouse stimulating with capsaicin combined with ammonia.

AIM OF THE STUDY

To investigate the metabolic modulating effects, and potential mechanisms underlying the therapeutic effect of PPRFT in chronic cough.

MATERIALS AND METHODS

Chronic cough mouse models were created by stimulating mice by capsaicin combined with ammonia. Number of coughs and cough latency within 2 min were recorded. With lung tissue and serum samples collected for histopathology, metabolomics, RT-qPCR, immunohistochemistry, and WB analysis. Lymphocytes were isolated and flow cytometric assays were conducted to evaluate the differentiation between Th17 and Treg cell among CD4+ cells.

RESULTS

Results indicated that PPRFT obviously reduced the number of coughs, prolonged cough latency, reduced inflammatory cell infiltration and lung tissues damage, and decreased the serum level of IL-6, IL-1β, TNF-α, and IL-17 while increasing IL-10 levels. Notably, PPRFT suppressed Th17 cell divergence and promoted Treg cell divergence. Furthermore, serum metabolomic assays showed that 46 metabolites differed significantly between group, with 35 pathways involved. Moreover, mRNA levels of IL-6, NF-κB, IL-17, RORγT, JAK2, STAT3, PI3K and AKT in lung tissues remarkably reduced and mRNA levels of IL-10 and FOXP3 were elevated after PPRFT pretreatment. Additionally, PPRFT treatments decreased the protein levels of IL-6, NF-κB, IL-17, RORγT, p-JAK2, p-STAT3, p-PI3K, and p-AKT and increased the protein levels of IL-10 and FOXP3, but no significantly effects to the levels on JAK2, STAT3, PI3K, and AKT in the lungs.

CONCLUSION

Conclusively, our result suggested the effect with PPRFT on chronic cough may be mediated through IL-6/JAK2/STAT3 and PI3K/AKT/NF-κB pathway, which regulate the differentiation between Th17 and Treg cell. This beneficial effect of PPRFT in capsaicin and ammonia-stimulated chronic cough mice indicates its potential application in treating chronic cough.

Integrated Analysis of lncRNA-miRNA-mRNA Regulatory Network in Rapamycin-Induced Cardioprotection against Ischemia/Reperfusion Injury in Diabetic Rabbits

The inhibition of mammalian target of rapamycin (mTOR) with rapamycin (RAPA) provides protection against myocardial ischemia/reperfusion (I/R) injury in diabetes. Since interactions between transcripts, including long non-coding RNA (lncRNA), microRNA(miRNA) and mRNA, regulate the pathophysiology of disease, we performed unbiased miRarray profiling in the heart of diabetic rabbits following I/R injury with/without RAPA treatment to identify differentially expressed (DE) miRNAs and their predicted targets of lncRNAs/mRNAs. Results showed that among the total of 806 unique miRNAs targets, 194 miRNAs were DE after I/R in diabetic rabbits. Specifically, eight miRNAs, including miR-199a-5p, miR-154-5p, miR-543-3p, miR-379-3p, miR-379-5p, miR-299-5p, miR-140-3p, and miR-497-5p, were upregulated and 10 miRNAs, including miR-1-3p, miR-1b, miR-29b-3p, miR-29c-3p, miR-30e-3p, miR-133c, miR-196c-3p, miR-322-5p, miR-499-5p, and miR-672-5p, were significantly downregulated after I/R injury. Interestingly, RAPA treatment significantly reversed these changes in miRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated the participation of miRNAs in the regulation of several signaling pathways related to I/R injury, including MAPK signaling and apoptosis. Furthermore, in diabetic hearts, the expression of lncRNAs, HOTAIR, and GAS5 were induced after I/R injury, but RAPA suppressed these lncRNAs. In contrast, MALAT1 was significantly reduced following I/R injury, with the increased expression of miR-199a-5p and suppression of its target, the anti-apoptotic protein Bcl-2. RAPA recovered MALAT1 expression with its sponging effect on miR-199-5p and restoration of Bcl-2 expression. The identification of novel targets from the transcriptome analysis in RAPA-treated diabetic hearts could potentially lead to the development of new therapeutic strategies for diabetic patients with myocardial infarction.

Genome-wide unraveling SNP pairwise epistatic effects associated with sheep body weight

Epistatic effects are an important part of the genetic effect of complex traits in livestock. In this study, we used 218 synthetic ewes from the Xinjiang Academy of Agricultural Reclamation in China to identify interacting paired with genome-wide single nucleotide polymorphisms (SNPs) associated with birth weight, weaning weight, and one-yearling weight. We detected 2 and 66 SNP-SNP interactions of sheep birth weight and weaning weight, respectively. No significant epistatic interaction of one-year-old body weight was detected. The genetic interaction of sheep body weight is dynamic and time-dependent. Most significant interactions of weaning body weight contributed 1% or higher. In the weaning weight trait, 66 significant SNP pairs consisted of 98 single SNPs covering 23 chromosomes, 5 of which were nonsynonymous SNPs (nsSNPs), resulting in single amino acid substitution. We found that genes that interact with transcription factors (TFs) are target genes for the corresponding TFs. Four epitron networks affecting weaning weight, including subnetworks of HIVEP3 and BACH2 transcription factors, constructed using significant SNP pairs, were also analyzed and annotated. These results suggest that transcription factors may play an important role in explaining epistatic effects. It provides a new idea to study the genetic mechanism of weight developing.

High expression of miR-30c-5p in satellite cells of high-fat diet-induced obese rabbits inhibited satellite cell proliferation and promoted differentiation

It was revealed in our previous study that the expression of miR-30c-5p in the skeletal muscle of rabbits fed high-fat diet was highly expressed. In the present study, we further investigated the function of miR-30c-5p in proliferation and differentiation of skeletal muscle satellite cell (SMSC). The results obtained in the present study showed that the skeletal muscle fibers of the rabbits fed the standard normal diet (SND) were orderly, regular, and uniform after HE staining, however, the muscle fibers of the rabbits fed the high-fat diet (HFD) were generally atrophied, some were arranged disorderly, the intercellular space was enlarged, the nucleus was increased, and the morphology and position were abnormal. Compared with the SND group, it was observed that the weekly weight gain and fat percentage were relatively higher, and also the levels of the serum biochemical indexes such as glucose, cholesterol, and triglyceride increased significantly in the rabbits fed with HFD. In addition, the results after the transfection of miR-30c-5p mimic, mimic NC (negative control), miR-30c-5p inhibitor, and inhibitor NC into the SMSCs showed that the cell counting kit-8 (CCK-8) proliferation experiment suggested that the number of cells in the over expression group was significantly lower than that in the mimic NC group at 48, 72, 96 h of cell proliferation. At 48, 72, 120 h, the number of cells in the inhibitor group was significantly higher than that in the mimic NC group. The number of EdU positive cells decreased significantly in the over expression group compared with the mimic NC group, however, it increased significantly in the inhibitor group compared with the inhibitor NC group. Moreover, compared with the mimic NC group, the myotube area increased significantly in the miR-30c-5p mimic group, whereas it decreased significantly in the miR-30c-5p inhibitor group as compared with the inhibitor NC group. In addition, we found that trinucleotide repeat containing adaptor 6A (TNRC6A) was successfully validated as a target gene for miR-30c-5p. The expression of TNRC6A in the miR-30c-5p mimic group was significantly lower than that in the mimic NC group. It was further observed that the expression of TNRC6A increased significantly in the miR-30c-5p inhibitor group as compared to that in the inhibitor NC group. Taken together, the results obtained in this study showed that miR-30c-5p promotes the differentiation as well as inhibits the proliferation of rabbit skeletal muscle satellite cells, and TNRC6A is a target gene of miR-30c-5p.

Intramuscular mitochondrial and lipid metabolic changes of rats after regular high-intensity interval training (HIIT) of different training periods

BACKGROUND

High-intensity Interval Training (HIIT) is a time-efficient form of exercise and has gained popularity in recent years. However, at molecular level, the understanding about the effects of HIIT is not comprehensive, and even less is elucidated about HIIT of different training duration cycles, although different durations always lead to different post-training consequences.

METHOD

In this study, by training SD rats using HIIT protocols lasting for different training duration cycles, we investigated the adaptive response of intramuscular triglyceride abundance as well as mitochondrial and lipid metabolic changes after HIIT training (2, 4, 6, 8, and 10 weeks). We selected 72 h after the last session of training as the time point of sacrifice.

RESULTS

The suppressed activation of the cAMP-PKA pathway indicates that skeletal muscle was in the recovery phase at this time point. Intramuscular triglyceride abundance was significantly elevated after 2, 4, and 10 weeks of HIIT. However, the lipid metabolism-related proteins inconsistently changed in a chaotic trend (see Table 1). The expression levels of PGC1-α and COX IV decreased after 2 and 4 weeks of training and raised after 6 and 8 weeks of training. The expression level of citrate synthase (CS) decreased after 2, 4, 8, and 10 weeks of training, and showed an upward trend after 6 weeks of training. While the activity of CS decreased after 2 and 8 weeks of training and showed an upward trend after 6 weeks of HIIT.

CONCLUSION

Given the abovementioned changing trends, we propose two speculations: (A) the damaged mitochondria oxidation capacity might be one of the causes of IMTG accumulation observed after 2 and 4 weeks of HIIT. This phase might be similar to the condition of type 2 diabetes. (B) after 6-week HIIT, mitochondria function and biogenesis might be improved and the IMTG contents declined to baseline. This might be explained as: mitochondrial enhancement increased the capacity of lipid oxidation and then offset the increase in IMTG achieved during the first 4 weeks. For HIIT Rat Modelling, if the aim is to observe HIIT-induced positive effects, caution should be exercised when considering 2 and 4 weeks of training under our HIIT frame. Also, implementing six-week training is at least effective for mitochondrial enhancement when using similar HIIT frame of this study.

Genome-wide association study of abdominal MRI-measured visceral fat: The multiethnic cohort adiposity phenotype study

Few studies have explored the genetic underpinnings of intra-abdominal visceral fat deposition, which varies substantially by sex and race/ethnicity. Among 1,787 participants in the Multiethnic Cohort (MEC)-Adiposity Phenotype Study (MEC-APS), we conducted a genome-wide association study (GWAS) of the percent visceral adiposity tissue (VAT) area out of the overall abdominal area, averaged across L1-L5 (%VAT), measured by abdominal magnetic resonance imaging (MRI). A genome-wide significant signal was found on chromosome 2q14.3 in the sex-combined GWAS (lead variant rs79837492: Beta per effect allele = -4.76; P = 2.62 × 10-8) and in the male-only GWAS (lead variant rs2968545: (Beta = -6.50; P = 1.09 × 10-9), and one suggestive variant was found at 13q12.11 in the female-only GWAS (rs79926925: Beta = 6.95; P = 8.15 × 10-8). The negatively associated variants were most common in European Americans (T allele of rs79837492; 5%) and African Americans (C allele of rs2968545; 5%) and not observed in Japanese Americans, whereas the positively associated variant was most common in Japanese Americans (C allele of rs79926925, 5%), which was all consistent with the racial/ethnic %VAT differences. In a validation step among UK Biobank participants (N = 23,699 of mainly British and Irish ancestry) with MRI-based VAT volume, both rs79837492 (Beta = -0.026, P = 0.019) and rs2968545 (Beta = -0.028, P = 0.010) were significantly associated in men only (n = 11,524). In the MEC-APS, the association between rs79926925 and plasma sex hormone binding globulin levels reached statistical significance in females, but not in males, with adjustment for total adiposity (Beta = -0.24; P = 0.028), on the log scale. Rs79837492 and rs2968545 are located in intron 5 of CNTNAP5, and rs79926925, in an intergenic region between GJB6 and CRYL1. These novel findings differing by sex and racial/ethnic group warrant replication in additional diverse studies with direct visceral fat measurements.

Astragalus Polysaccharide Alleviates Constipation in the Elderly Via Modification of Gut Microbiota and Fecal Metabolism

Constipation is one of the most common gastrointestinal disorders, whose incidence increasing with age. As one of the main components, Astragalus polysaccharide (APS) has been used to treat a variety of diseases. This study aimed to explore the effects of APS on the improvement of gastrointestinal functions and learning memory in elderly rats with constipation. In this study, both 16S rRNA sequencing-based microbiome and ¹H NMR-based metabolomics were applied to demonstrate the effects of APS on host metabolism and gut microbiota of the elderly rats with constipation. On top of this, we constructed both inter- and inner-layer networks, intuitively showing the correlations among behavioral indicators, intestinal bacteria, and differential metabolites. Our results showed that APS significantly ameliorated the constipation and the cognitive dysfunctions of rats. Microbiome analysis revealed that APS raised the relative abundance of Blautia, whereas decreased the relative abundance of Lactobacillus in the elderly rats with constipation. In addition, APS decreased the levels of acetate, butyrate, and propionate in the fecal samples, correspondingly regulating glycolysis/gluconeogenesis metabolism and pyruvate metabolism. These findings lay solid foundations for understanding the pathogenesis of constipation in the elderly, and also offer a promising new treatment strategy for constipation in the elderly.

Estimating the Prognosis of Low-Grade Glioma with Gene Attention Using Multi-Omics and Multi-Modal Schemes

The prognosis estimation of low-grade glioma (LGG) patients with deep learning models using gene expression data has been extensively studied in recent years. However, the deep learning models used in these studies do not utilize the latest deep learning techniques, such as residual learning and ensemble learning. To address this limitation, in this study, a deep learning model using multi-omics and multi-modal schemes, namely the Multi-Prognosis Estimation Network (Multi-PEN), is proposed. When using Multi-PEN, gene attention layers are employed for each datatype, including mRNA and miRNA, thereby allowing us to identify prognostic genes. Additionally, recent developments in deep learning, such as residual learning and layer normalization, are utilized. As a result, Multi-PEN demonstrates competitive performance compared to conventional models for prognosis estimation. Furthermore, the most significant prognostic mRNA and miRNA were identified using the attention layers in Multi-PEN. For instance, MYBL1 was identified as the most significant prognostic mRNA. Such a result accords with the findings in existing studies that have demonstrated that MYBL1 regulates cell survival, proliferation, and differentiation. Additionally, hsa-mir-421 was identified as the most significant prognostic miRNA, and it has been extensively reported that hsa-mir-421 is highly associated with various cancers. These results indicate that the estimations of Multi-PEN are valid and reliable and showcase Multi-PEN's capacity to present hypotheses regarding prognostic mRNAs and miRNAs.