Differential mRNAs Research Articles

Phloridzin Reveals New Treatment Strategies for Liver Fibrosis.

Liver fibrosis is an urgent public health problem which is difficult to resolve. However, various drugs for the treatment of liver fibrosis in clinical practice have their own problems during use. In this study, we used phloridzin to treat hepatic fibrosis in the CCl4-induced C57/BL6N mouse model, which was extracted from lychee core, a traditional Chinese medicine. The therapeutic effect was evaluated by biochemical index detections and ultrasound detection. Furthermore, in order to determine the mechanism of phloridzin in the treatment of liver fibrosis, we performed high-throughput sequencing of mRNA and lncRNA in different groups of liver tissues. The results showed that compared with the model group, the phloridzin-treated groups revealed a significant decrease in collagen deposition and decreased levels of serum alanine aminotransferase, aspartate aminotransferase, laminin, and hyaluronic acid. GO and KEGG pathway enrichment analysis of the differential mRNAs was performed and revealed that phloridzin mainly affects cell ferroptosis. Gene co-expression analysis showed that the target genes of lncRNA were obvious in cell components such as focal adhesions, intercellular adhesion, and cell–substrate junctions and in metabolic pathways such as carbon metabolism. These results showed that phloridizin can effectively treat liver fibrosis, and the mechanism may involve ferroptosis, carbon metabolism, and related changes in biomechanics.

The Construction and Analysis of Infiltrating Immune Cell and ceRNA Networks in Diabetic Foot Ulcer.

BackgroundDiabetic foot ulcer (DFU) is a severe complication characterized by low-grade infectious inflammation and probably associated with specific competitive endogenous RNAs (ceRNAs) and infiltrating immune cells. Nonetheless, no reliable biomarkers are used for detecting infectious inflammation in DFU. Therefore, it is essential to explore potential biomarkers for the accurate diagnosis and treatment of DFU.MethodsThe gene expression profile was retrieved from Gene Expression Omnibus (GEO) database and divided into two groups, namely, standard samples and DFU samples. To establish the ceRNA networks, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were utilized to analyze differential expression genes (DEGs). The cell type identification was achieved by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm to screen-specific immune-infiltrating cells associated with DFU.ResultsA ceRNA network was constructed with 20 differential expression circRNA (DEcircRNAs), 11 differential expression microRNAs (DEmiRNAs), and 9 differential expression mRNAs (DEmRNAs). Functional enrichment analysis demonstrated that DFU was mainly enriched in vascular endothelial growth factor (VEGF) and T-cell receptor signaling. In addition, CIBERSORT estimation indicated that CD8+ T cells and Monocytes were significantly related to the expression of IL-6, a DFU-specific infectious inflammation factor.ConclusionThis study identified that some significant ceRNAs (JUNB, GATA3, hsa-circ-0049271 and hsa-circ-0074559) and infiltrating immune cells (CD8+ T cells and monocytes) might be related to DFU infectious inflammation.

Screening of Prognostic Markers for Hepatocellular Carcinoma Patients Based on Multichip Combined Analysis.

Methods GSE (14520, 36376, 57957, 76427) datasets were accessed from GEO database. 55 differential mRNAs (DEGs) were obtained by differential analysis based on the datasets. GO and KEGG analysis results indicated that the DEGs were enriched in xenobiotic metabolic process and other pathways. Expression profiles and clinical data of TCGA-LIHC mRNAs were from TCGA database. We established a prognostic model of HCC through univariate and multivariate Cox risk regression analyses. ROC curve analysis was used to examine the prognostic model performance. GSEA analysis was performed between the high- and low-risk score sample groups. Results A 4-gene HCC prognostic model was constructed, in which the gene expressions correlated to HCC patients' survival. The AUC value presented 0.734 in the ROC analysis for the prognostic model. Conclusion The four-gene model could be introduced as an independent prognostic factors to assess HCC patients' survival status.

Transcriptome Profiling of N7-Methylguanosine Modification of Messenger RNA in Drug-Resistant Acute Myeloid Leukemia.

Acute myeloid leukemia (AML) is an aggressive hematological tumor caused by the malignant transformation of myeloid progenitor cells. Although intensive chemotherapy leads to an initial therapeutic response, relapse due to drug resistance remains a significant challenge. In recent years, accumulating evidence has suggested that post-transcriptional methylation modifications are strongly associated with tumorigenesis. However, the mRNA profile of m7G modification in AML and its role in drug-resistant AML are unknown. In this study, we used MeRIP-seq technology to establish the first transcriptome-wide m7G methylome profile for AML and drug-resistant AML cells, and differences in m7G between the two groups were analyzed. In addition, bioinformatics analysis was conducted to explore the function of m7G-specific methylated transcripts. We found significant differences in m7G mRNA modification between AML and drug-resistant AML cells. Furthermore, bioinformatics analysis revealed that differential m7G-modified mRNAs were associated with a wide range of cellular functions. Importantly, down-methylated m7G modification was significantly enriched in ABC transporter-related mRNAs, which are widely recognized to play a key role in multidrug resistance. Our results provide new insights into a novel function of m7G methylation in drug resistance progression of AML.

Identification of key miRNAs regulating fat metabolism based on RNA-seq from fat-tailed sheep and F2 of wild Argali

The evolution mechanism of sheep tail fat has not yet been clear, many researches focus on this issue, yet there still many gaps to be filled in the targets and non-coding RNA regulation. In our study, the differential expression mRNAs and miRNAs were detected by RNA-seq and constructed a mRNA-miRNA network related to the lipid deposition in tail of fat-tailed sheep and F2 of Argali with domestic sheep (thin-tailed). Then 6 kinds of tissues from thin-tailed and control group were extracted for function validation of candidate genes and its regulator miRNAs. 125 differentially expressed miRNAs were identified by RNA-seq, and enrichment analysis of their target genes revealed 10 significantly enriched pathways related to lipid metabolism. In these pathways, 126 DE-miRNA target genes were also differentially expression in the same tissues in our previous transcriptomic data. In PPI network, 6 hubgenes (SCD, ACACA, GPD2, ELOVL6, ELOVL5, GPAM) were extracted using the cytoHubba application, and they may be target genes for 3 candidate DE-miRNAs (miR-320d, miR-151b, miR-6715). The validation results of RT-qPCR show: the expression trend of miR-320d is opposite to the target gene SCD, and that of miR-151b and the target gene ACACA are also opposite in 6 tissues, implying that they may have direct targeting relationships. Moreover, the expression of miR-320d in F2 tail fat was significantly higher than that in fat-tailed sheep (P < 0.05), and the expression of SCD in F2 tail fat was extremely significantly lower than that in fat-tailed sheep (P < 0.01). The expression of miR-151b in F2 tail fat and subcutaneous fat was significantly higher than that in fat-tailed sheep (P < 0.05), and the expression of ACACA in F2 subcutaneous fat was significantly lower than that in fat-tailed sheep. miR-320d may directly and negatively regulate tail fat deposition by targeting SCD, while miR-151b may indirectly and negatively regulate tail fat deposition by targeting ACACA.

MicroRNA-27a-3p inhibits lung and skin fibrosis of systemic sclerosis by negatively regulating SPP1

ObjectiveTo investigate the role and mechanism of microRNAs (miRNAs) in fibrotic processes involved in the pathology of systemic sclerosis (SSc). MethodsR language and bioinformatics methods were used to identify differential miRNAs and mRNAs and analyze their biological functions. Transfection experiments were performed to evaluate the function and regulatory mechanism of miR-27a-3p in vitro. Levels of fibrosis-related genes, SPP1 and cell proliferation were assessed. ResultsMiR-27a-3p is reduced both in SSc lung and skin tissues. Overexpression of miR-27a-3p significantly inhibited fibrosis-related genes expression and protein abundance and cell proliferation, whereas inhibition of miR-27a-3p significantly enhanced these phenomena. Moreover, miR-27a-3p exerts its anti-fibrosis effect by negatively regulating SPP1 and ERK signal, more prominent in fibroblasts. ConclusionsOur findings show that miR-27a-3p regulates a common mechanism in the process of SSc skin and lung fibrosis. MiR-27a-3p/SPP1/ERK1/2 axis may be an important target for delaying the progression of SSc fibrosis.

Potential Value of Circular RNA circTBC1D4 in Gastrointestinal Stromal Tumors.

Aims To explore the expression of circular RNA (circRNA) in gastrointestinal stromal tumors. Background Gastrointestinal stromal tumors (GIST) are mainly distributed in the stomach and small intestine. Recently, it has been verified that circular RNA (circRNA) has an important function in the regulation of GIST. Nevertheless, detailed investigations of circRNA-miRNA-mRNA regulatory networks in GIST are lacking. Objective To analyze the gastrointestinal stromal tumor circRNA-miRNA-mRNA network, assessing the effect of circle RNA in gastrointestinal stromal tumors. Method All the differential circRNAs and mRNAs were obtained from Gene Expression Omnibus (GEO) microarray data (GSE131481 and GSE147303, GSE131481, and GSE13861). Furthermore, a circRNA-miRNA-mRNA network was established. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment were used to reveal the correlation between the functions of signaling pathways and target genes. The hub genes of protein-protein interaction (PPI) network and cytoHubba were also defined. Quantitative real-time PCR (qRT-PCR) was used to measure the expression levels of hsa-circ-0002917 (circTBC1D4), hsa-miR-590-5p (miR-590-5p), and PLN. Results PPI network and Cytoscape showed that ATP1A2, PLN, KCNMA1, and SCNN1B were four central DEGs. GO analysis results revealed that DEGs were involved in negative management of myocardial contraction, regulation of myocardial cell contraction, ethanol oxidation, cellular potassium ion homeostasis, and relaxation of cardiac muscle, and KEGG analysis showed that major DEGs were with cGMP-PKG signaling pathway. Moreover, we obtained two pairs of axes, namely, hsa-circ-0039216/hsa-miR-338-3p/ATP1A2 and hsa-circ-0002917/hsa-miR-590-5p/PLN. The target of TBC1D4 is miR-590-5p, and miR-590-5p increased after knocking down TBC1D4. Moreover, PLN was the target of miR-590-5p, and miR-590-5p exerts antitumor effects by reducing PLN. Conclusions In this study, we constructed a circRNA-miRNA-mRNA management network interrelated with GIST and researched the potential roles of circRNA. Moreover, we discovered a new molecular landmarker for the prediction, diagnosis, and therapy of patients.

Thyroid Transcriptomic Profiling Reveals the Follicular Phase Differential Regulation of lncRNA and mRNA Related to Prolificacy in Small Tail Han Sheep with Two FecB Genotypes.

Long non-coding RNA (lncRNA) accounts for a large proportion of RNA in animals. The thyroid gland has been established as an important gland involved in animal reproduction, however, little is known of its gene expression patterns and potential roles in the sheep. Herein, RNA-Seq was used to detect reproduction-related differentially expressed lncRNAs (DELs) and mRNAs (DEGs) in the follicular phase (FT) FecBBB (MM) and FecB++ (ww) genotypes of Small Tail Han (STH) sheep thyroids. Overall, 29 DELs and 448 DEGs in thyroid between MM and ww sheep were screened. Moreover, GO and KEGG enrichment analysis showed that targets of DELs and DEGs were annotated in biological transitions, such as cell cycle, oocyte meiosis and methylation, which in turn affect reproductive performance in sheep. In addition, we constructed co-expression and networks of lncRNAs-mRNAs. Specifically, XLOC_075176 targeted MYB, XLOC_014695 targeted VCAN, 106991527 targeted CASR, XLOC_075176 targeted KIFC1, XLOC_360232 targeted BRCA2. All these differential lncRNAs and mRNAs expression profiles in the thyroid provide a new resource for elucidating the regulatory mechanism underlying STH sheep prolificacy.

Predict the role of lncRNA in kidney aging based on RNA sequencing

BackgroundLong noncoding RNAs (lncRNAs) are involved in physiological and pathological processes. However, no studies have been conducted on the relationship between lncRNAs and renal aging.ResultsFirst, we evaluated the histopathology of young (3-month-old) and old (24-month-old) C57BL/6J mouse kidneys. Masson trichrome staining and PAS staining showed interstitial collagen deposition and fibrosis, mesangial matrix expansion, a thicker basement membrane and renal interstitial fibrosis in old mouse kidneys. Senescence-associated β-galactosidase (SA-β-gal)-positive areas in the kidneys of old mice were significantly elevated compared to those of young mice. Then, we analyzed the differential expression of lncRNAs and mRNAs in the kidneys of young and old mouse kidneys by RNA-seq analysis. 42 known and 179 novel differentially expressed lncRNAs and 702 differential mRNAs were detected in the mouse kidney. Next, we focused on the differentially expressed mRNAs and lncRNAs by RNA-seq. GO and KEGG analyses were performed based on differentially expressed mRNAs between young and old mouse kidneys. Transregulation based on RIsearch and the correlation coefficient of mRNA-lncRNA were also calculated. The mRNA-lncRNA network was constructed by choosing a Spearman correlation coefficient > 0.9 or <-0.9. GO and KEGG pathway enrichment analyses revealed that differentially expressed mRNAs participated in aging-related pathways. A total of 10 lncRNAs and trans-regulated mRNAs were constructed. Finally, we validated the role of lncRNA Gm43360 by CCK-8, flow cytometry, western blot and SA-β-gal staining. The expression level of Adra1a was positively correlated and Csnk1a1 was negatively correlated with lncRNA Gm43360. The cell counting kit-8 (CCK-8) results showed that lncRNA Gm43360 promoted cell viability. LncRNA Gm43360 increased the percentage of S phase cells and decreased the percentage of G1 phase cells compared with the negative control. LncRNA Gm43360 decreased the expression of p53, p21 and SA-β-gal.ConclusionsLncRNA Gm43360 may play a protective role in kidney aging.

Identification of lncRNAs and Their Regulatory Relationships with mRNAs in Response to Cryptococcus neoformans Infection of THP-1 Cells.

Aims Cryptococcosis is an invasive fungal disease that is associated with an increasing prevalence along with a very high fatality and is primarily caused by Cryptococcus. However, its mechanism to cause pathogenicity is not yet completely understood. In this study, we aim to screen the lncRNA markers in human monocytic (THP-1) cells infected by Cryptococcus neoformans (C. neoformans) through high-throughput sequencing technology and to explore its effects on biological functions. Methods We initially conducted an lncRNA microarray analysis of the THP-1 cells infected by C. neoformans and normal THP-1 cells. Based upon these data, RT-qPCR was used to verify the expressions of the selected lncRNAs and mRNAs. We then performed functional and pathway enrichment analyses. Lastly, target prediction was performed by using the lncRNA target tool which was based on the differentially expressed lncRNAs. Results We determined 81 upregulated and 96 downregulated lncRNAs using microarray. In addition, the profiling data showed 42 upregulated and 57 downregulated genes and discovered that neuroactive ligand-receptor interaction, tyrosine metabolism, and phenylalanine metabolism are extremely impaired in the regulation of C. neoformans infection. GO enrichment analysis of the 99 differentially expressed mRNAs exhibited that these modules showed different signaling pathways and biological mechanisms like protein binding and metal ion binding. Moreover, lncRNAs and mRNAs were analyzed for their coexpression relations. A qRT-PCR analysis confirmed that the expression of the top 10 differently expressed mRNA and lincRNA. The expressions of the lncRNAs after C. neoformans infection in THP-1 cells were detected by RNA-sequence, suggesting that microarray analysis could reveal lncRNAs having functional significance that might be linked with the progression of patients. Conclusion The current study analyzed the differential lncRNAs and mRNAs in C. neoformans infection and predicted the corresponding pathways and their correlations that can offer new potential insights into the mechanistic basis of this condition.

Global analysis of miRNA-mRNA regulation pair in bladder cancer

PurposeMicroRNA (miRNA) is a class of short non-coding RNA molecules that functions in RNA silencing and post-transcriptional regulation of gene expression. This study aims to identify critical miRNA-mRNA regulation pairs contributing to bladder cancer (BLCA) pathogenesis.Patients and methodsMiRNA and mRNA microarray and RNA-sequencing datasets were downloaded from gene expression omnibus (GEO) and the cancer genome atlas (TCGA) databases. The tool of GEO2R and R packages were used to screen differential miRNAs (DE-miRNAs) and mRNAs (DE-mRNAs) and DAVID, DIANA, and Hiplot tools were used to perform gene enrichment analysis. The miRNA-mRNA regulation pair were screened from the experimentally validated miRNA-target interactions databases (miRTarbase and TarBase). Twenty-eight pairs of BLCA tissues were used to further verify the screened DE-miRNAs and DE-mRNAs by quantitative reverse transcription and polymerase chain reaction (qRT-PCR). The diagnostic value of the miRNA-mRNA regulation pairs was evaluated by receiver operating characteristic curve (ROC) and decision curve analysis (DCA). The correlation analysis between the selected miRNA-mRNAs regulation pair and clinical, survival and tumor-related phenotypes was performed in this study.ResultsAfter miRTarBase, the analysis of 2 miRNA datasets, 6 mRNA datasets, and TCGA-BLCA dataset, a total of 13 miRNAs (5 downregulated and 8 upregulated in BLCA tissues) and 181 mRNAs (72 upregulated and 109 downregulated in BLCA tissues) were screened out. The pairs of miR-17-5p (upregulated in BLCA tissues) and TGFBR2 (downregulated in BLCA tissues) were verified in the external validation cohort (28 BLCA vs. 28 NC) using qRT-PCR. Areas under the ROC curve of the miRNA-mRNA regulation pair panel were 0.929 (95% CI: 0.885–0.972, p < 0.0001) in TCGA-BLCA and 0.767 (95% CI: 0.643–0.891, p = 0.001) in the external validation. The DCA also showed that the miRNA-mRNA regulation pairs had an excellent diagnostic performance distinguishing BLCA from normal controls. Correlation analysis showed that miR-17-5p and TGFBR2 correlated with tumor immunity.ConclusionsThe research identified potential miRNA-mRNA regulation pairs, providing a new idea for exploring the genesis and development of BLCA.

The role of exosomal CDC6 in the hirudin-mediated suppression of the malignant phenotype of bladder cancer cells

BackgroundBladder cancer is a malignant tumor characterized by high recurrence and persistence due to the limited therapies that are currently available. Hirudin exerts a strong anticancer effect on several tumors. Thus, it is urgent to explore the biological function of hirudin in bladder cancer and the role of bladder cancer-derived exosomes in tumor inhibition. MethodsFirst, a network pharmacology analysis was performed to explore the relationships among hirudin, bladder cancer, and exosomes. Then, the effects of hirudin were examined by CCK-8 assay, flow cytometry, Transwell assay, and tumorigenic ability experiments in vitro. Exosomes derived from cells were identified with transmission electron microscopy, fluorescence labeling, and Western blotting and collected for further microarray analysis. Only CDC6 expression and mRNA abundance in hirudin-treated cells and exosomes was subjected to further validation using quantitative PCR and Western blotting. ResultsThrough network analysis, we found that hirudin affected bladder cancer, and this effect was related to exosomes. Our studies verified the effects of hirudin by revealing that hirudin inhibits malignant processes of bladder cancer cells in vitro, such as invasion, metastasis, and apoptosis. Similarly, the oncogenic effects of bladder cancer-derived exosomes were successfully isolated and identified. Via microarray assessment of the exosomes, we identified 600 differential mRNAs, of which the expression of the core target CDC6 was found to be significantly different in both The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. We further confirmed that hirudin suppresses CDC6 expression mRNA abundance in both cells and exosomes. ConclusionHirudin was able to decrease the expression of CDC6 in bladder cancer cells and exosomes, which effectively repressed the malignant processes of bladder cancer cells.

Plasma tRNA-derived small RNAs signature as a predictive and prognostic biomarker in lung adenocarcinoma

BackgroundThe prevalence of lung adenocarcinoma (LUAD) has increased, thus novel biomarkers for its early diagnosis is becoming more important than ever. tRNA-derived small RNA (tsRNA) is a new class of non-coding RNA which has important regulatory roles in cancer biology. This study was designed to identify novel predictive and prognostic tsRNA biomarkers.MethodstsRNAs were identified and performed differential expression analysis from 10 plasma samples (6 LUAD and 4 normal, SRP266333) and 96 tissue samples (48 LUAD and 48 normal, SRP133217). Then a tsRNA-mRNA regulatory network was constructed to find hub tsRNAs. Functional enrichment analysis was performed to infer the potential pathways associated with tsRNAs. Afterwards, a Support Vector Machine (SVM) algorithm was used to explore the potential biomarkers for diagnosing LUAD. Lastly, the function of tRF-21-RK9P4P9L0 was explored in A549 and H1299 cell lines.ResultsA significant difference of read distribution was observed between normal people and LUAD patients whether in plasma or tissue. A tsRNA-mRNA regulatory network consisting of 155 DEtsRNAs (differential expression tsRNAs) and 406 DEmRNAs (differential expression mRNAs) was established. Three tsRNAs (tRF-16-L85J3KE, tRF-21-RK9P4P9L0 and tRF-16-PSQP4PE) were identified as hub genes with degree > 100. We found Co-DEmRNAs (intersection of DEtsRNAs target mRNAs and differentially expressed mRNAs in LUAD) were engaged in a number of cancer pathways. The AUC of the three hub tsRNAs’ expression for diagnosing LUAD reached 0.92. Furthermore, the qPCR validation of the three hub tsRNAs in 37 paired normal and LUAD tissues was consistent with the RNA-Seq results. In addition, tRF-21-RK9P4P9L0 was negatively associated with LUAD prognosis. Inhibition of tRF-21-RK9P4P9L0 expression reduced the proliferation, migration and invasion ability of A549 and H1299 cell lines.ConclusionThese findings will help us further understand the molecular mechanisms of LUAD and contribute to novel diagnostic biomarkers and therapeutic target discovery.

Screening and identification of differentially expressed long non-coding RNAs in multidrug-resistant tuberculosis.

BackgroundEfforts to eradicate tuberculosis are largely threatened by drug-resistant tuberculosis, particularly, multidrug-resistant tuberculosis (MDR-TB). Screening and identification potential biomarkers for MDR-TB is crucial to diagnose early and reduce the incidence of MDR-TB.MethodsTo screen the differentially expressed long non-coding RNAs in MDR-TB, the lncRNA and mRNA expression profiles in serum derived from healthy controls (HCs), individuals with MDR-TB and drug-sensitive tuberculosis (DS-TB) were analyzed by microarray assay and 10 lncRNAs were randomly selected for further validation by reverse transcription-quantitative real-time PCR(RT-qPCR). The biological functions of differentially expressed mRNAs as well as relationships between genes and signaling pathways were investigated using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG), respectively.ResultsA total of 353 differentially expressed lncRNAs (312 upregulated) and 202 mRNAs (99 upregulated) were found in the MDR-TB group compared to HCs. And compared with the DS-TB group, 442 differentially expressed lncRNAs (115 upregulated) and 190 mRNAs (87 upregulated) were found in the MDR-TB group. The expression levels of lncRNA n335659 were found to differ significantly between each group by RT-qPCR. Compared with DS-TB group, the GO analysis showed that the differential mRNAs were mainly enriched in the processes associated with the detection of the chemical stimulus, the regulation of mRNA metabolic process and neutrophil activation in the MDR-TB group; the KEGG analysis indicated that the differential mRNAs between DS-TB and MDR-TB were mainly enriched in proteasome and Notch signaling pathway, which might reveal a fraction of the mechanism of MDR-TB. The discovery of the serum lncRNA n335659 might serve as a potential biomarker for MDR-TB and Notch signaling pathway provided a new clue for the investigation of the pathological mechanism of MDR-TB.

Integrated Dissection of lncRNA-miRNA-mRNA Pairs and Potential Regulatory Role of lncRNA PCAT19 in Lung Adenocarcinoma.

Lung adenocarcinoma (LUAD) is the main cause of cancer-related deaths worldwide. Long noncoding RNAs have been reported to play an important role in various cancers due to their special functions. Therefore, identifying the lncRNAs involved in LUAD tumorigenesis and development can help improve therapeutic strategies. The TCGA-LUAD RNA expression profile was downloaded from The Cancer Genome Atlas, and a total of 49 differential lncRNAs, 112 differential miRNAs, and 2,953 differential mRNAs were screened. Through Kaplan–Meier curves, interaction networks, hub RNAs (lncRNAs, miRNAs, and mRNAs) were obtained. These hub genes are mainly involved in cell proliferation, cell cycle, lung development, and tumor-related signaling pathways. Two lncRNAs (SMIM25 and PCAT19) more significantly related to the prognosis of LUAD were screened by univariate Cox analysis, multivariate Cox analysis, and risk model analysis. The qPCR results showed that the expression levels of SMIM25 and PCAT19 were downregulated in clinical tissues, A549 and SPC-A1 cells, which were consistent with the bioinformatics analysis results. Subsequently, the PCAT19/miR-143-3p pairs were screened through the weighted gene co-expression network analysis and miRNA-lncRNA regulatory network. Dual luciferase detection confirmed that miR-143-3p directly targets PCAT19, and qPCR results indicated that the expression of the two is positively correlated. Cell function tests showed that overexpression of PCAT19 could significantly inhibit the proliferation, migration, and invasion of A549 and SPC-A1 cells. In contrast, knockout of PCAT19 can better promote the proliferation and migration of A549 and SPC-A1 cells. The expression of PCAT19 was negatively correlated with tumor grade, histological grade, and tumor mutation load in LUAD. In addition, co-transfection experiments confirmed that the miR-143-3p mimic could partially reverse the effect of PCAT19 knockout on the proliferation of A549 and SPC-A1 cells. In summary, PCAT19 is an independent prognostic factor in patients with LUAD that can regulate the proliferation, migration, and invasion of LUAD cells and may be a potential biomarker for the diagnosis of LUAD. PCAT19/miR-143-3p plays a very important regulatory role in the occurrence and development of LUAD.

Integrative Analysis of Long Non-coding RNAs, Messenger RNAs, and MicroRNAs Indicates the Neurodevelopmental Dysfunction in the Hippocampus of Gut Microbiota-Dysbiosis Mice

Major depressive disorder is caused by gene–environment interactions and the gut microbiota plays a pivotal role in the development of depression. However, the underlying mechanisms remain elusive. Herein, the differentially expressed hippocampal long non-coding RNAs (lncRNAs), messenger RNAs (mRNAs), and microRNAs (miRNAs) between mice inoculated with gut microbiota from major depressive disorder patients or healthy controls were detected, to identify the effects of gut microbiota-dysbiosis on gene regulation patterns at the transcriptome level, and in further to explore the microbial-regulated pathological mechanisms of depression. As a result, 200 mRNAs, 358 lncRNAs, and 4 miRNAs were differentially expressed between the two groups. Functional analysis of these differential mRNAs indicated dysregulated inflammatory response to be the primary pathological change. Intersecting these differential mRNAs with targets of differentially expressed miRNAs identified 47 intersected mRNAs, which were mainly related to neurodevelopment. Additionally, a microbial-regulated lncRNA–miRNA–mRNA network based on RNA–RNA interactions was constructed. Subsequently, according to the competitive endogenous RNAs (ceRNA) hypothesis and the biological functions of these intersected genes, two neurodevelopmental ceRNA sub-networks implicating in depression were identified, one including two lncRNAs (4930417H01Rik and AI480526), one miRNA (mmu-miR-883b-3p) and two mRNAs (Adcy1 and Nr4a2), and the other including six lncRNAs (5930412G12Rik, 6430628N08Rik, A530013C23Rik, A930007I19Rik, Gm15489, and Gm16251), one miRNA (mmu-miR-377-3p) and three mRNAs (Six4, Stx16, and Ube3a), and these molecules could be recognized as potential genetic and epigenetic biomarkers in microbial-associated depression. This study provides new understanding of the pathogenesis of depression induced by gut microbiota-dysbiosis and may act as a theoretical basis for the development of gut microbiota-based antidepressants.

MicroRNA‐mRNA genetic networks in refractive eye development and myopia

AbstractPurposeMyopia is a complex disease controlled by both environmental and genetic factors. Recent genetic studies suggested that the development of myopia is controlled by numerous coding genes. However, the contribution of non‐coding genes, such as microRNAs, is poorly understood. The purpose of this study was to investigate the role of microRNAs in refractive eye development and the development of myopia.MethodsDifferential expression of both microRNAs and mRNAs was analyzed in the retina, RPE, choroid, and sclera of C57BL mice upon induction of experimental myopia with −20 D lenses using RNA‐seq and microRNA‐mRNA genetic networks were reconstructed. We also explored the impact of microRNAs on signaling pathways and evaluated the role of microRNAs in refractive eye development in vivo using conditional Dicer1 knockout mice.ResultsGene expression profiling revealed large number of differential mRNAs and microRNAs in the retina and RPE/choroid/sclera complexes. MicroRNAs were involved in the regulation of several key canonical pathways and biological processes previously implicated in refractive eye development. We then compared visual acuity, contrast sensitivity, and refractive errors in Dicer1 knockouts (Dicer1−/−) and wild‐type littermates (Dicer1+/+) at baseline. We found that loss of miRNAs in the postnatal eyes resulted in a large hyperopic shift in refractive error (+16.3 ± 3.6 D in Dicer1 knockouts versus +4.7 ± 3.0 D in wild‐type animals, p &lt; 1.44 × 10−4), which was accompanied by a small reduction in visual acuity (0.360 ± 0.006 cpd in Dicer1−/− versus 0.403 ± 0.002 cpd in Dicer1+/+, p &lt; 1.16 × 10−7) and contrast sensitivity (12.01 ± 0.27 in Dicer1−/− versus 13.65 ± 0.40 in Dicer1+/+, p &lt; 1.91 × 10−5). Loss of miRNAs also caused a statistically significant increase in susceptibility to form‐deprivation myopia in Dicer1 knockouts compared to wild‐type animals (−19.1 ± 6.8 D in Dicer1−/− versus −9.7 ± 2.7 D in Dicer1+/+, p &lt; 1.51 × 10−4).ConclusionsOur data suggest that microRNAs play significant role in refractive eye development and the development of myopia via modulation of several canonical signaling pathways. MicroRNAs appear to be integrated in the signaling cascade involved in the optical‐defocus‐regulated refractive eye development at multiple levels.

BRAFV600E mutant colorectal cancer cells mediate local immunosuppressive microenvironment through exosomal long noncoding RNAs.

BACKGROUNDBRAFV600E mutated colorectal cancer (CRC) is prone to peritoneal and distant lymph node metastasis and this correlates with a poor prognosis. The BRAFV600E mutation is closely related to the formation of an immunosuppressive microenvironment. However, the correlation between BRAFV600E mutation and changes in local immune microenvironment of CRC is not clear.AIMTo explore the effect and mechanism of BRAFV600E mutant on the immune microenvironment of CRC.METHODSThirty patients with CRC were included in this study: 20 in a control group and 10 in a treatment group. The density of microvessels and microlymphatic vessels, and M2 subtype macrophages in tumor tissues were detected by immunohistochemistry. Screening and functional analysis of exosomal long noncoding RNAs (lncRNAs) were performed by transcriptomics. The proliferation and migration of human umbilical vein endothelial cells (HUVECs) and human lymphatic endothelial cells (HLECs) were detected by CCK-8 assay and scratch test, respectively. The tube-forming ability of endothelial cells was detected by tube formation assay. The macrophage subtypes were obtained by flow cytometry. The expression of vascular endothelial growth factor (VEGF)-A, basic fibroblast growth factor (bFGF), transforming growth factor (TGF)-β1, VEGF-C, claudin-5, occludin, zonula occludens (ZO)-1, fibroblast activation protein, and α-smooth muscle actin was assessed by western blot analysis. The levels of cytokines interleukin (IL)-6, TGF-β1, and VEGF were assessed by enzyme-linked immunosorbent assay.RESULTSBRAFV600E mutation was positively correlated with the increase of preoperative serum carbohydrate antigen 19-9 (P < 0.05), and with poor tumor tissue differentiation in CRC (P < 0.01). Microvascular density and microlymphatic vessel density in BRAFV600E mutant CRC tissues were higher than those in BRAF wild-type CRC (P < 0.05). The number of CD163+ M2 macrophages in BRAFV600E mutant CRC tumor tissue was markedly increased (P < 0.05). Compared with exosomes from CRC cells with BRAF gene silencing, the expression of 13 lncRNAs and 192 mRNAs in the exosomes from BRAFV600E mutant CRC cells was upregulated, and the expression of 22 lncRNAs and 236 mRNAs was downregulated (P < 0.05). The biological functions and signaling pathways predicted by differential lncRNA target genes and differential mRNAs were closely related to angiogenesis, tumor cell proliferation, differentiation, metabolism, and changes in the microenvironment. The proliferation, migration, and tube formation ability of HUVECs and HLECs induced by exosomes in the 1627 cell group (HT29 cells with BRAF gene silencing) was greatly reduced compared with the HT29 cell group (P < 0.05). Compared with the HT29 cell group, the expression levels of VEGF-A, bFGF, TGF-β1, and VEGF-C in the exosomes derived from 1627 cells were reduced. The expression of ZO-1 in HUVECs, and claudin-5, occludin, and ZO-1 in HLECs of the 1627 cell group was higher. Compared with the 1627 cell group, the exosomes of the HT29 cell group promoted the expression of CD163 in macrophages (P < 0.05). IL-6 secretion by macrophages in the HT29 cell group was markedly elevated (P < 0.05), whereas TGF-β1 was decreased (P < 0.05). The levels of IL-6, TGF-β1, and VEGF secreted by fibroblasts in the 1627 cell group decreased, compared with the HT29 cell group (P < 0.05).CONCLUSIONBRAFV600E mutant CRC cells can reach the tumor microenvironment by releasing exosomal lncRNAs, and induce the formation of an immunosuppressive microenvironment.

Study on differential gene expression profile of serum exosomes in patients with acute cerebral infarction

ObjectiveTo analyze the differential gene expression profile of serum exosomes in patients with acute cerebral infarction (ACI) and clarify the changes in gene expression related to cerebral infarction injury and the potential serum markers. MethodsFour patients with ACI and five healthy people were enrolled in the Phase Ⅰ study. After serum isolation from peripheral blood, exosomes were extracted with exosomes kits, high-throughput detection of mRNA was performed with gene chips, and differentially expressed mRNAs were screened. Gene Ontology (GO) functional analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed simultaneously. Furthermore, real-time polymerase chain reaction (qRT-PCR) was used to verify the expression levels of the screened differential mRNAs in the serum exosomes collected in Phase Ⅱ from 32 patients each in the ACI case and normal control groups. ResultsIn the Phase Ⅰ study, there were 248 differentially expressed mRNAs (fold change ≥ 2.0, P < 0.05) among five patients in the normal control group and four patients in the case group, of which the expression of 242 was upregulated and that of six was downregulated. The results of GO functional enrichment analysis mainly included behavior regulation, cell connection, and antioxidant activity. The results of KEGG pathway enrichment analysis mainly included ribosomes, proteasomes, oxytocin signaling pathways, and oxidative phosphorylation. After researching and screening based on relevant literature, it was found that among the genes with significant differential expression, H3F3B mRNA may be associated with and might play an important role in ACI. The qRT-PCR method was used to detect the H3F3B mRNA expression in serum exosomes of 32 patients each in the normal control and case groups in Phase Ⅱ; the expression was significantly higher in serum exosomes of the case group than in those of the normal control group (P < 0.001). H3F3B mRNA expression in serum exosomes of the case group positively correlated with age, the National Institutes of Health Stroke Scale (NIHSS) score, and the maximum infarct size (P < 0.05). ConclusionACI can lead to changes in the serum exosomes mRNA expression profile, which may be closely related to the occurrence, development, and prognosis of this condition. These findings will provide direction for research on the molecular mechanism, diagnostic markers, and therapeutic targets of ACI.

Screening and Analysis of Key Genes in miRNA-mRNA Regulatory Network of Membranous Nephropathy.

Background MicroRNAs (miRNAs) are confirmed to participate in occurrence, development, and prevention of membranous nephropathy (MN), but their mechanism of action is unclear. Objective With the GEO database and the use of bioinformatics, miRNA-mRNA regulatory network genes relevant to MN were explored and their potential mechanism of action was explained. Methods The MN-related miRNA chip data set (GSE51674) and mRNA chip data set (GSE108109) were downloaded from the GEO database. Differential analysis was performed using the GEO2R online tool. TargetScan, miRTarBase, and StarBase databases were used to predict potential downstream target genes regulated by differentially expressed miRNAs, and the intersection with differential genes were taken to obtain candidate target genes. According to the regulatory relationship between miRNA and mRNA, the miRNA-mRNA relationship pair was clarified and Cytoscape was used to construct a miRNA-mRNA regulatory network. WebGestalt was used to conduct enrichment analysis of the biological process of differential mRNAs in the regulatory network; FunRich analyzes the differential mRNA pathways in the miRNA-mRNA regulatory network. And the STRING database was used to construct a PPI network for candidate target genes, and Cytoscape visually analyzes the PPI network. Results Experiments were conducted to screen differentially expressed miRNAs and mRNAs. There were 30 differentially expressed miRNAs, including 22 upregulated and 8 downregulated; and 1267 differentially expressed mRNAs, including 536 upregulated and 731 downregulated. Using TargetScan, miRTarBase, and StarBase databases to predict the downstream targets of differentially expressed miRNAs, 2957 downstream target genes coexisting in the 3 databases were predicted to intersect with differentially expressed mRNAs to obtain 175 candidate target genes. Finally, 36 miRNA-mRNA relationship pairs comprising 10 differentially expressed miRNAs and 27 differentially expressed mRNAs were screened out, and the regulatory network was constructed. Further analysis revealed that the miRNA regulatory network genes may be involved in the development of membranous nephropathy by mTOR, PDGFR-β, LKB1, and VEGF/VEGFR signaling pathways. Conclusion The miRNA regulatory network genes may participate in the regulation of podocyte autophagy, lipid metabolism, and renal fibrosis through mTOR, PDGFR-β, LKB1, and VEGF/VEGFR signaling pathways, thereby affecting the occurrence and development of membranous nephropathy.