mRNA Co-expression Network Research Articles

Gene expression-based modeling of overall survival in Black or African American patients with lung adenocarcinoma

IntroductionLung cancer is a leading cause of cancer-related deaths worldwide. Black/African American (B/AA) populations, in particular, exhibit the highest incidence and mortality rates of lung adenocarcinoma (LUAD) in the United States.MethodsThis study aims to explore gene expression patterns linked to LUAD in B/AA and case-matched white patients, with the goal of developing predictive models for prognosis. Leveraging RNA sequencing data from The Cancer Genome Atlas (TCGA) database, genes and pathways associated with overall survival (OS) were identified.ResultsThe OS-associated genes in B/AA patients were distinct from those in white patients, showing predominant enrichment in immune-related pathways. Furthermore, mRNA co-expression network analysis revealed that OS-associated genes in B/AA patients had higher levels of interaction with various pathways, including those related to immunity, cell-ECM interaction, and specific intracellular signaling pathways. Notably, a potential B/AA-specific biomarker, C9orf64, demonstrated significant correlations with genes involved in immune response. Unsupervised machine learning algorithms stratified B/AA patients into groups with distinct survival outcomes, while supervised algorithms demonstrated a higher accuracy in predicting survival for B/AA LUAD patients compared to white patients.DiscussionIn total, this study explored OS-associated genes and pathways specific for B/AA LUAD patients. Further validation and clinical application of these findings are warranted to address disparities and improve outcomes in diverse patient populations.

Comprehensive analysis of an mRNA co-expression network and a ceRNA network reveals potential prognostic biomarkers in oral squamous cell carcinoma

BackgroundOral squamous cell carcinoma (OSCC) is a prevalent and aggressive oral cancer with a poor prognosis. Its polygenic risk is likely influenced by complex transcriptional disorders involving networks of co-expressed and functionally related genes, though such investigations are limited in OSCC.MethodsWe analyzed the GSE37991 dataset, comprising 40 OSCC and 40 normal oral tissue samples from the Gene Expression Omnibus. Tumor-specific modules were identified using weighted correlation network analysis (WGCNA), leading to the selection of hub mRNAs and lncRNAs. These lncRNAs were used to construct lncRNA–mRNA and competing endogenous RNA networks. We further examined the expression profiles and survival data of these genes from the Cancer Genome Atlas. Prognostic markers were identified and validated through 5-year survival analysis and Cox proportional hazards modeling. RT-qPCR was used to validate the expression levels in clinical OSCC tissues.ResultsWe identified 1847 differentially expressed genes in OSCC tissues. WGCNA revealed four OSCC-specific modules, screening 120 hub mRNAs and five hub lncRNAs. Two prognostic markers (AQP5, IL-26) from hub mRNAs and three (FRMD5, INHBB, GUCY1A3) from the lncRNA–mRNA network were associated with survival. Validation showed lower expression of AQP5 and GUCY1A3, and higher expression of FRMD5 and INHBB in OSCC compared to normal tissues.ConclusionThis study enhances our understanding of transcriptional dysregulation in OSCC and may highlights AQP5, IL-26, FRMD5, INHBB, and GUCY1A3 as promising prognostic biomarkers.

Nasal microRNA signatures for disease severity in infants with respiratory syncytial virus bronchiolitis: a multicentre prospective study

BackgroundRespiratory syncytial virus (RSV) bronchiolitis contributes to a large morbidity and mortality burden globally. While emerging evidence suggests that airway microRNA (miRNA) is involved in the pathobiology of RSV infection,...

RNA sequencing reveals the implication of the circRNA-associated ceRNA network in oesophageal squamous cell carcinoma.

Circular RNAs (circRNAs) have attracted increasing attention in cancer research. However, there are few studies about the high-throughput sequencing for clinical cohorts focussing on the expression characteristics and regulatory networks of circRNAs in oesophageal squamous cell carcinoma (ESCC) until now. Present study aim to comprehensively recognize the functional and mechanistic patterns of circRNA through constructing a circRNA-related competing endogenous RNA (ceRNA) network in ESCC. Summarily, RNA high-throughput sequencing was adopted to assess the circRNA, miRNA and mRNA expression profiles in ESCC. Through bioinformatics methods, a circRNA-miRNA-mRNA coexpression network was constructed and hub genes was identified. Finally, cellular function experiments combined with bioinformatics analysis were conducted to verify the identified circRNA was involved in the progression of ESCC through ceRNA mechanism. In this study, we established a ceRNA regulatory network, including 5 circRNAs, 7 miRNAs and 197 target mRNAs, and 20 hub genes were screened and identified to exert important roles in the progression of ESCC. As a verification, hsa_circ_0002470 (circIFI6) was revealed to be highly expressed in ESCC and regulate the expression of hub genes by absorbing miR-497-5p and miR-195-5p through ceRNA mechanism. Our results further indicated that silencing of circIFI6 repressed proliferation and migration of ESCC cells, highlighting the tumour promotion effects of circIFI6 in ESCC. Collectively, our study contributes a new insight into the progression of ESCC from the perspective of the circRNA-miRNA-mRNA network, shedding light on the circRNA research in ESCC.

Characterization of the circulating transcriptome expression profile and identification of novel miRNA biomarkers in hypertrophic cardiomyopathy

BackgroundHypertrophic cardiomyopathy (HCM), one of the most common genetic cardiovascular diseases, but cannot be explained by single genetic factors. Circulating microRNAs (miRNAs) are stable and highly conserved. Inflammation and immune response participate in HCM pathophysiology, but whether the miRNA profile changes correspondingly in human peripheral blood mononuclear cells (PBMCs) with HCM is unclear. Herein, we aimed to investigate the circulating non-coding RNA (ncRNA) expression profile in PBMCs and identify potential miRNAs for HCM biomarkers.MethodsA Custom CeRNA Human Gene Expression Microarray was used to identify differentially expressed (DE) mRNAs, miRNAs, and ncRNAs (including circRNA and lncRNA) in HCM PBMCs. Weighted correlation network analysis (WGCNA) was used to identify HCM-related miRNA and mRNA modules. The mRNAs and miRNAs from the key modules were used to construct a co-expression network. Three separate machine learning algorithms (random forest, support vector machine, and logistic regression) were applied to identify potential biomarkers based on miRNAs from the HCM co-expression network. Gene Expression Omnibus (GEO) database (GSE188324) and experimental samples were used for further verification. Gene set enrichment analysis (GSEA) and competing endogenous RNA (ceRNA) network was used to determine the potential functions of the selected miRNAs in HCM.ResultsWe identified 1194 DE-mRNAs, 232 DE-miRNAs and 7696 DE-ncRNAs in HCM samples compared with normal controls from the microarray data sets. WGCNA identified key miRNA modules and mRNA modules evidently associated with HCM. We constructed a miRNA‒mRNA co-expression network based on these modules. A total of three hub miRNAs (miR-924, miR-98 and miR-1) were identified by random forest, and the areas under the receiver operator characteristic curves of miR-924, miR-98 and miR-1 were 0.829, 0.866, and 0.866, respectively.ConclusionsWe elucidated the transcriptome expression profile in PBMCs and identified three hub miRNAs (miR-924, miR-98 and miR-1) as potential biomarkers for HCM detection.

Long non-coding RNAs regulate heavy metal-induced apoptosis in embryo-derived cells

Heavy metal pollution has been a worldwide prevalent problem, and particularly a threat to ecosystem integrity and animals’ health. Previous studies on the mechanisms of heavy metal toxicity have focused on protein-coding genes, whereas most genomic transcripts are long non-coding RNAs (lncRNAs). Although lncRNAs are known to play important regulatory roles in biological processes, their role in heavy metal stress regulation is still not fully understood. We here developed an insect embryo cell model for studying metal toxicity and the underlying regulatory mechanisms. We performed genome-wide screening and functional characterization of lncRNAs induced by two essential and two non-essential heavy metals in Drosophila embryo-derived S2 cells. We identified 4894 lncRNAs, of which 1410 were novel. Forty-one lncRNAs, together with 328 mRNAs, were induced by all the four heavy metals. LncRNA–mRNA co-expression network and pathway enrichment analysis showed that detoxification metabolism, circadian rhythm, and apoptosis regulation pathways were activated in response to heavy metal stress. LncRNA CR44138 was remarkably upregulated in cells exposed to the four heavy metals and was associated with the apoptosis pathway. Expression interference confirmed that CR44138 aggravated cytotoxicity-induced apoptosis in cells under heavy metals stress. This study highlights the important role of lncRNAs in regulating the cellular response to heavy metals. This study also lays the foundation for discovering the novel regulatory mechanisms and developing diagnostic biomarkers of the toxic effects of heavy metal pollutants on organisms.

Potential Function of Testicular MicroRNAs in Heat-Stress-Induced Spermatogenesis Disorders.

Spermatogenesis is temperature-dependent, and the increase in testicular temperature seriously affects mammalian spermatogenesis and semen quality. In this study, the testicular heat stress model of mice was made with a 43 °C water bath for 25 min, and the effects of heat stress on semen quality and spermatogenesis-related regulators were analyzed. On the 7th day after heat stress, testis weight shrank to 68.45% and sperm density dropped to 33.20%. High-throughput sequencing analysis showed that 98 microRNAs (miRNAs) and 369 mRNAs were down-regulated, while 77 miRNAs and 1424 mRNAs were up-regulated after heat stress. Through gene ontology (GO) analysis of differentially expressed genes and miRNA-mRNA co-expression networks, it was found that heat stress may be involved in the regulation of testicular atrophy and spermatogenesis disorders by affecting cell meiosis process and cell cycle. In addition, through functional enrichment analysis, co-expression regulatory network, correlation analysis and in vitro experiment, it was found that miR-143-3p may be a representative potential key regulatory factor affecting spermatogenesis under heat stress. In summary, our results enrich the understanding of miRNAs in testicular heat stress and provide a reference for the prevention and treatment of heat-stress-induced spermatogenesis disorders.

MRNA and lncRNA co-expression network in mice of acute intracerebral hemorrhage.

Intracerebral hemorrhage (ICH) is a severe subtype of stroke lacking effective pharmacological targets. Long noncoding RNA (lncRNA) has been confirmed to participate in the pathophysiological progress of various neurological disorders. However, how lncRNA affects ICH outcomes in the acute phase is not completely clear. In this study, we aimed to reveal the relationship of lncRNA-miRNA-mRNA following ICH. We conducted the autologous blood injection ICH model and extracted total RNAs on day 7. Microarray scanning was used to obtain mRNA and lncRNA profiles, which were validated by RT-qPCR. GO/KEGG analysis of differentially expressed mRNAs was performed using the Metascape platform. We calculated the Pearson correlation coefficients (PCCs) of lncRNA-mRNA for co-expression network construction. A competitive endogenous (Ce-RNA) network was established based on DIANALncBase and miRDB database. Finally, the Ce-RNA network was visualized and analyzed by Cytoscape. In total, 570 differentially expressed mRNAs and 313 differentially expressed lncRNAs were identified (FC ≥ 2 and value of p <0.05). The function of differentially expressed mRNAs was mainly enriched in immune response, inflammation, apoptosis, ferroptosis, and other typical pathways. The lncRNA-mRNA co-expression network contained 57 nodes (21 lncRNAs and 36 mRNAs) and 38 lncRNA-mRNA pairs. The ce-RNA network was generated with 303 nodes (29 lncRNAs, 163 mRNAs, and 111 miRNAs) and 906 edges. Three hub clusters were selected to indicate the most significant lncRNA-miRNA-mRNA interactions. Our study suggests that the top differentially expressed RNA molecules may be the biomarker of acute ICH. Furthermore, the hub lncRNA-mRNA pairs and lncRNA-miRNA-mRNA correlations may provide new clues for ICH treatment.

Analysis of the regulatory role of miR-34a-5p/PLCD3 in the progression of osteoarthritis

Osteoarthritis is a heterogeneous disease with a complex etiology. However, there is no effective treatment strategy at present. The purpose of this study was to explore the miRNA‒mRNA regulatory network and molecular mechanism that regulate the progression of osteoarthritis. In this article, we downloaded datasets (GSE55457, GSE82107, GSE143514 and GSE55235) from Gene Expression Omnibus (GEO) to screen differentially expressed mRNAs in osteoarthritis. Then, through weighted gene coexpression network (WGCNA), functional enrichment, protein‒protein interaction (PPI) network, miRNA‒mRNA coexpression network, ROC curve, and immune infiltration analyses and qPCR, the mRNA PLCD3, which was highly expressed in osteoarthritis and had clinical predictive value, was screened. We found that PLCD3 directly targets miR-34a-5p through DIANA and dual-luciferase experiments. The expression levels of PLCD3 and miR-34a-5p were negatively correlated. In addition, CCK-8 and wound healing assays showed that the miR-34a-5p mimic inhibited hFLS-OA cell proliferation and promoted hFLS-OA cell migration. PLCD3 overexpression showed the opposite trend. Western blotting further found that overexpression of miR-34a-5p reduced the protein expression levels of p-PI3K and p-AKT, while overexpression of PLCD3 showed the opposite trend. In addition, combined with the effect of the PI3K/AKT pathway inhibitor BIO (IC50 = 5.95μM), the results showed that overexpression of miR-34a-5p increased the inhibitory effects of BIO on p-PI3K and p-AKT protein expression, while overexpression of PLCD3 significantly reversed these inhibitory effects. Overall, the miR-34a-5p/PLCD3 axis may mediate the PI3K/AKT pathway in regulating cartilage homeostasis in synovial osteoarthritis. These data indicate that miR-34a-5p/PLCD3 may be a new prognostic factor in the pathology of synovial osteoarthritis.

Identification of Recurrence-Related mRNAs and Noncoding RNAs in Hepatocellular Carcinoma Following Liver Transplantation.

This study aimed to determine whether altered mRNA, long non-coding RNA, and circular RNA expression is related to hepatocellular carcinoma recurrence after liver transplantation. Five recurrent and 5 non-recurrent primary hepatocellular carcinoma samples were used to perform RNA sequencing. Then, differentially expressed mRNAs, differentially expressed long non-coding RNAs, and differentially expressed circular RNAs between recurrent group and non-recurrent group were identified. In addition, differentially expressed long non-coding RNA/differentially expressed circular RNA-differentially expressed mRNA co-expression network, and competing endogenous RNA (differentially expressed circular RNA/differentially expressed long non-coding RNA-miRNA-differentially expressed mRNA) regulatory network were constructed. Finally, real-time quantitative polymerase chain reaction was performed for verification. Five hundred forty-one differentially expressed mRNAs, 239 differentially expressed long non-coding RNAs, and 16 differentially expressed circular RNAs in the recurrent group were obtained. Gene set enrichment analysis indicated that these differentially expressed mRNAs may affect hepatocellular carcinoma recurrence through multiple pathways, such as E2F, epithelial-mesenchymal transition, G2M, and oxidative phosphorylation. Then, 993 differentially expressed long non-coding RNA-differentially expressed mRNA co-expression pairs and 28 differentially expressed circular RNA/differentially expressed mRNA co-expression pairs were obtained. The competing endogenous RNA network contained 10 circular RNA-miRNA pairs, 12 long non-coding RNA-miRNA pairs, and 36 miRNA- mRNA pairs. Real-time quantitative polymerase chain reaction results showed the same expression trend as RNA-seq results. Our results reveal key mRNAs, long non-coding RNAs, and circular RNAs associated with recurrence in hepatocellular carcinoma patients after liver transplantation.

Genome-wide analysis of circular RNA-mediated ceRNA regulation in porcine skeletal muscle development

BackgroundAs a diverse and abundant class of endogenous RNAs, circular RNAs (circRNAs) participate in various biological processes including cell proliferation and apoptosis. Nevertheless, few researchers have investigated the role of circRNAs in muscle development in cultivated pigs.ResultsIn this study, we used RNA-seq to construct circRNA expression profiles in skeletal muscle of Jinfen White pigs at the age of 1, 90, and 180 days. Among the 16,990 identified circRNAs, 584 circRNAs were differentially expressed. Moreover, the enrichment analysis of DE circRNA host genes showed that they were mainly involved in muscle contraction, muscle organ development and muscle system processes, as well as AMPK and cAMP-related signal pathways. We also constructed a circRNA–miRNA–mRNA co-expression network to find key circRNAs which many involved in the regulation of porcine skeletal muscle development through the competitive endogenous RNA (ceRNA) mechanism. It is noteworthy that circ_0018595/miR-1343/PGM1 axis may play a regulatory role in the development of porcine skeletal muscle.ConclusionsThis study identified the circRNAs and present the circRNA expression profile in the development of pigs, revealed that DE circRNA host genes participate in different cell fates and enriched the porcine ceRNA network. Thus, this work will become a valuable resource for further in-depth study of the regulatory mechanism of circRNA in the development of porcine skeletal muscle.

Co-expression network of mRNA and DNA methylation in first-episode and drug-naive adolescents with major depressive disorder.

We explored the DNA methylation and messenger RNA (mRNA) co-expression network and hub genes in first-episode, drug-naive adolescents with major depressive disorder (MDD). To preliminarily explore whether adolescent MDD has unique mechanisms compared with adult MDD. We compared DNA methylation and mRNA profiles of peripheral blood mononuclear cells from four first-episode and drug-naive adolescents with MDD and five healthy adolescent controls (HCs). We performed differential expression analysis, constructed co-expression network, and screened the hub genes. And enrichment analysis was performed based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). We also downloaded DNA methylation and mRNA datasets of adults with MDD (GSE113725/GSE38206) from the GEO database, and performed differential expression and enrichment analysis. Our clinical data showed that 3034 methylation sites and 4190 mRNAs were differentially expressed in first-episode, drug-naive adolescents MDD patients compared with HCs. 19 hub genes were screened out according to the high degree value in the co-expression network. The results from the GEO database showed that compared with adult HCs, there were 290 methylation sites and 127 mRNAs were differentially expressed in adult MDD patients. Compared with adolescent HCs and adult MDD patients, the DNA methylation and mRNA expression patterns of first-episode, drug-naive adolescent MDD patients were different. The co-expression network of DNA methylation and mRNA and the screened hub genes may play an important role in the pathogenesis of MDD in first-episode, drug-naive adolescents. Compared with adult MDD, adolescent MDD is more enriched in metabolism in terms of function and pathways.

An exosome-related long non-coding RNAs risk model could predict survival outcomes in patients with breast cancer

Breast cancer (BC) is one of the most frequent malignancies among women worldwide. Accumulating evidence indicates that long non-coding RNA (lncRNA) may affect BC progression. Exosomes, a class of small membrane vesicles, have been reported to promote tumor progression through transporting proteins, mRNAs, lncRNAs and some other small molecules. However, the interaction between exosome-related lncRNAs and the microenvironment of malignancies is unclear. Hence, we proceeded to investigate the relationship between exosome-related lncRNAs and BC microenvironment. 121 exosome-associated genes were extracted from ExoBCD database. Then, the Pearson analysis was used to screened out the exosome-related lncRNAs. After that, 15 exosome-related differentially expressed lncRNAs were identified by the correlation with BC prognosis. According to the sum of the expression of these 15 lncRNAs, extracted from The Cancer Genome Atlas, and the regression coefficients, an exosome-related lncRNAs signature was developed by using Cox regression analysis. With the median risk score of the training set, the patients in training and validation sets were separated to low-risk group and high-risk group. Subsequently, the lncRNA–mRNA co-expression network was constructed. The distinct enrichment pathways were compared among the different risk groups by using the R package clusterProfiler. The ESTIMATE method and ssGESA database were adopted to study the ESTIMATE Score and immune cell infiltration. Eventually, the expression of immune checkpoint associated genes, microsatellite instable and the immunophenoscore were further analyzed between different risk groups. Different risk groups exhibited different prognosis, with lower survival rate in the high-risk group. The differentially expressed genes between the different risk groups were enriched in biological processes pathways as well as immune responses. BC patients in high-risk group were identified with lower scores of ESTIMATE scores. Subsequently, we noticed that the infiltrating levels of aDCs, B cells, CD8+ T cells, iDCs, DCs, Neutrophils, macrophages, NK cells, pDCs, Tfh, T helper cells, TIL and Tregs were obvious elevated with the decreased risk score in training and validation cohorts. And some immune signatures were significantly activated with the decreased risk score in both cohorts. Eventually, the exosome-associated lncRNAs risk model was demonstrated to accurately predict immunotherapy response in patients with BC. The results of our study suggest that exosome-related lncRNAs risk model has close relationship with prognosis and immune cells infiltration in BC patients. These findings could make a great contribution to improving BC immunotherapy.

Comprehensive analysis of m5C-Related lncRNAs in the prognosis and immune landscape of hepatocellular carcinoma.

5-Methyladenosine (m5C) is a type of epigenetic modification involved in the progression of various cancers. To investigate the role of m5C-related long non-coding RNAs (lncRNAs) in the prognosis and immune cell infiltration in hepatocellular carcinoma (HCC), we obtained patients' clinical information and transcriptome data of HCC from the Cancer Genome Atlas (TCGA) database. We applied Pearson correlation analysis to construct an m5C-related lncRNA-messenger RNA (mRNA) co-expression network. Univariate Cox analysis, least absolute shrinkage and selection operator (LASSO), and multivariate Cox analysis were employed to establish an m5C-related lncRNA prognostic risk model. We then verified the model using Kaplan-Meier analysis, principal component analysis, as well as univariate and multivariate Cox analyses. The expression of m5C-related lncRNAs was validated in HCC tissues and different cell lines. Combining the risk score and clinicopathological features, a nomogram was established for predicting the overall survival (OS) of HCC patients. Furthermore, gene set enrichment analysis (GSEA) revealed that some tumor-associated pathways were significantly enriched in the high-risk group. Immune cell infiltration analysis demonstrated that the levels of Treg cells, neutrophils, and M2 macrophages were higher in the high-risk group. In addition, patients with high tumor mutation burden (TMB) had worse OS than those with low TMB. We also assessed the immune checkpoint level and chemotherapeutic agent sensibility. Then in vitro experiments were performed to examine the biological function of MKLN1-AS in HCC cells and found that knockdown of MKLN1-AS suppressed the proliferation, migration, and invasion. In conclusion, m5C-related lncRNAs played a critical role in predicting the prognosis of patients with HCC and may serve as new therapeutic targets for HCC patients.

Immune-relevant genes of systemic lupus erythematosus by transcriptome profiling analysis

ObjectiveOur study aimed to reveal the roles of long noncoding RNA (lncRNA) and immune-relevant genes in systemic lupus erythematosus (SLE). MethodsExpression profiling dataset GSE65391 consisted of 72 healthy blood samples and 924 SLE blood samples was downloaded from the Gene Expression Omnibus. Differentially expressed RNAs (DERs) in SLE samples were identified using the Limma package. Immune-relevant DERs were uncovered after assessing the immune cell infiltration types of different samples. Modules significantly related to the disease were extracted via weighted gene co-expression network analysis (WGCNA), followed by module–trait analysis. LncRNA–mRNA co-expression network was constructed for DERs in immune-relevant modules. Genes related to SLE were further revealed via the Comparative Toxicogenomics Database (CTD). Validation assay was conducted by another independent expression profiling dataset, GSE46907 (consisted of 5 healthy blood samples and 5 SLE blood samples) and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) analysis of 15 healthy blood samples and 15 SLE blood samples. ResultsA total of 606 DERs, including 19 lncRNAs and 587 mRNAs, were identified. Obvious differences were observed in the proportions of 20 immune cell types, and 378 immune-relevant DERs were revealed. Genes in the lncRNA–mRNA co-expression network were significantly enriched in primary immunodeficiency, hematopoietic cell lineage, B cell receptor signaling pathway, SLE (up-regulated FCGR3A, down-regulated CD40LG, up-regulated HIST1H2BE, down-regulated human leukocyte antigen [HLA]-DOA, down-regulated HLA-DOB, up-regulated HIST1H3D, and up-regulated HIST1H2BC), and intestinal immune network for IgA production. SLE-related gene, CD40LG, retrieved from CTD has co-expression interactions with seven differentially expressed lncRNAs (HCG27, LINC02555, LINC02210, DHRS4-AS1, MIR600HG, DANCR, and LINC01278). CD40LG and HLA-DOA, were observed down-regulated, and FCGR3A, and HIST1H2BE were up-regulated in validation dataset GES46907. Moreover, CD40LG was validated to be down-regulated using qRT-PCR. ConclusionsOur results provide new insight in understanding the pathogenesis of SLE and might be helpful for developing new therapeutic approaches of SLE by modulating immune response.

RNA-Seq Comprehensive Analysis Reveals the Long Noncoding RNA Expression Profile and Coexpressed mRNA in Adult Degenerative Scoliosis.

Objective: Owing to the intensification of the aging process worldwide, the prevalence of adult degenerative scoliosis (ADS) is increasing at an alarming rate. However, genomic research related to the etiology of ADS is rarely reported worldwide. Since long noncoding RNAs (lncRNAs) play a pivotal role in the progression of human diseases, this study aimed to investigate ADS-associated messenger RNAs (mRNAs) and lncRNAs by RNA sequencing (RNA-seq), as well as performed comprehensive bioinformatics analysis based on the lncRNA–mRNA coexpression network and protein–protein interaction (PPI) network. Methods: Initially, six whole blood (WB) samples were obtained from three ADS and three nondegenerative lumbar trauma patients who underwent surgical operation for RNA-seq exploration to construct differential mRNA and lncRNA expression profiles. Subsequently, quantitative RT-PCR (qRT-PCR) was performed to validate three randomly selected differentially expressed mRNAs and lncRNAs derived from the nucleus pulposus (NP) tissue of 14 other subjects (seven ADS patients and seven nondegenerative lumbar trauma patients), respectively. Results: A total of 1,651 upregulated and 1,524 downregulated mRNAs and 147 upregulated and 83 downregulated lncRNAs were screened out from the RNA-Seq data, which constructed coexpression networks to investigate their regulatory interactions further. GO gene function prediction revealed that lncRNA-targeted genes might play a vital role in ADS via participation in multiple biological processes such as the AMPK signaling pathway, lysosomes, and ubiquitin-mediated proteolysis, as well as cellular metabolic processes. Moreover, the expression levels of three selected lncRNAs and mRNAs were validated by qRT-PCR, respectively, demonstrating that the relative expression levels were consistent with the RNA-seq data. Notably, the dysregulated RNAs, AKT1, UBA52, PTPN12, and CLEC16A, were significantly differentially expressed in ADS WB samples and might serve as potentially regulated genes for research in the future. Conclusions: This study provides the first insight into the altered transcriptome profile of long-stranded noncoding RNAs associated with ADS, which paves the way for further exploration of the clinical biomarkers and molecular regulatory mechanisms for this poorly understood degenerative disease. However, the detailed biological mechanisms underlying these candidate lncRNAs in ADS necessitate further elucidation in future studies.

Construction of lncRNA and mRNA co-expression network associated with nasopharyngeal carcinoma progression.

Nasopharyngeal carcinoma is a type of head and neck cancer with a high incidence in men. In the past decades, the survival rate of NPC has remained around 70%, but it often leads to treatment failure due to its distant metastasis or recurrence. The lncRNA-mRNA regulatory network has not been fully elucidated. We downloaded the NPC-related gene expression datasets GSE53819 and GSE12452 from the Gene Expression Omnibus database; GSE53819 included 18 NPC tissues and 18 normal tissues, and GSE12452 included 31 NPC tissues and 10 normal tissues. Weighted gene co-expression network analysis was performed on mRNA and lncRNA to screen out modules that were highly correlated with tumor progression. The two datasets were subjected to differential analysis after removing batch effects, and then Venn diagrams were used to screen for overlapping genes in the module genes and differential genes. The lncRNA-mRNA co-expression network was then constructed, and key mRNAs were identified by MCODE analysis and expression analysis. GSEA analysis and qRT-PCR were performed on key mRNAs. Through a series of analyses, we speculated that BTK, CD72, PTPN6, and VAV1 may be independent predictors of the prognosis of NPC patients.Taken together, our study provides potential candidate biomarkers for NPC diagnosis, prognosis, or precise treatment.

CircRNA, lncRNA, and mRNA profiles of umbilical cord blood exosomes from preterm newborns showing bronchopulmonary dysplasia

Bronchopulmonary dysplasia (BPD) represents a multifactorial chronic pulmonary pathology and a major factor causing premature illness and death. The therapeutic role of exosomes in BPD has been feverishly investigated. Meanwhile, the potential roles of exosomal circRNAs, lncRNAs, and mRNAs in umbilical cord blood (UCB) serum have not been studied. This study aimed to detect the expression profiles of circRNAs, lncRNAs, and mRNAs in UCB-derived exosomes of infants with BPD. Microarray analysis was performed to compare the RNA profiles of UCB-derived exosomes of a preterm newborn with (BPD group) and without (non-BPD, NBPD group) BPD. Then, circRNA/lncRNA–miRNA–mRNA co-expression networks were built to determine their association with BPD. In addition, cell counting kit-8 (CCK-8) assay was used to evaluate the proliferation of lipopolysaccharide (LPS)-induced human bronchial epithelial cells (BEAS-2B cells) and human umbilical vein endothelial cells (HUVECs). The levels of tumor necrosis factor (TNF)-α and interleukin (IL)-1β in LPS-induced BEAS-2B cells and HUVECs were assessed through Western blot analysis. Then, quantitative reverse transcription–polymerase chain reaction assay was used to evaluate the expression levels of four differentially expressed circRNAs (hsa_circ_0086913, hsa_circ_0049170, hsa_circ_0087059, and hsa_circ_0065188) and two lncRNAs (small nucleolar RNA host gene 20 (SNHG20) and LINC00582) detected in LPS-induced BEAS-2B cells or HUVECs. A total of 317 circRNAs, 104 lncRNAs, and 135 mRNAs showed significant differential expression in UCB-derived exosomes of preterm infants with BPD compared with those with NBPD. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted to examine differentially expressed exosomal circRNAs, lncRNAs, and mRNAs. The results showed that the GO terms and KEGG pathways mostly involving differentially expressed exosomal RNAs were closely associated with endothelial or epithelial cell development. In vitro, CCK-8 and Western blot assays revealed that LPS remarkably inhibited the viability and promoted inflammatory responses (TNF-α and IL-1β) of BEAS-2B cells or HUVECs. The expression levels of circRNAs hsa_circ_0049170 and hsa_circ_0087059 were upregulated in LPS-induced BEAS-2B cells; the expression level of hsa_circ_0086913 was upregulated and that of hsa_circ_0065188 was downregulated in LPS-induced HUVECs. Moreover, the expression level of lncRNA SNHG20 was upregulated and that of LINC00582 was downregulated in LPS-induced BEAS-2B cells. Further, 455 circRNA/lncRNA–miRNA–mRNA interaction networks were predicted, including hsa_circ_0086913/hsa-miR-103a-3p/transmembrane 4 L six family member 1 (TM4SF1) and lncRNA-SNHG20/hsa-miR-6720-5p/spermine synthase (SMS) networks, which may take part in BPD.Conclusion: This study provided a systematic perspective on UCB-derived exosomal circRNAs and lncRNAs and laid an important foundation for further investigating the potential biological functions of exosomal circRNAs and lncRNAs in BPD.What is Known:• BPD represents a multifactorial chronic pulmonary pathology and a major factor causing premature illness and death.• The therapeutic role of exosomes in BPD has been feverishly investigated, and exosomal RNAs were ignored.What is New:• The profiles of UCB-derived exosomal circRNAs, lncRNAs, and mRNAs were performed.• Several differentially expressed circRNAs and lncRNAs were identified in LPS-induced BEAS-2B cells and HUVECs.

Analysis of lncRNA and mRNA Repertoires in Lung From BAFF-R-Deficient Pneumocystis-Infected Mice.

BackgroundPneumocystis pneumonia (PCP) is a common medical issue in immunosuppressive patients. Increasing evidence supports that B cells may play an essential role in PCP individuals. The present study aims to integrate lncRNA and mRNA expression profiles and further investigate the molecular function of mature B cells in PCP.MethodsThe lung tissue of wild-type (WT) mice and B-cell-activating factor receptor–deficient (mature B-cell deficiency, BAFF-R–/–) mice were harvested at 3 weeks after being infected with pneumocystis. After total RNAs were extracted, transcriptome profiling was performed following the Illumina HiSeq 3000 protocol. lncRNA-targeted miRNA pairs were predicted using the online databases. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment pathways were analyzed to functionally annotate these differentially expressed genes. Additionally, the immune-related lncRNA–miRNA–mRNA–ceRNA network was subsequently performed. The quantitative real-time PCR (RT-PCR) analysis was conducted to evaluate the lncRNA and mRNA expression profiles in WT-PCP mice and BAFF-R–/– PCP mice.ResultsCompared with the control group, 166 mRNAs were observed to be aberrantly expressed (fold change value ≥2; P <0.05) in the BAFF-R–/– PCP group, including 39 upregulated and 127 downregulated genes, while there were 69 lncRNAs differently expressed in the BAFF-R–/– PCP group, including 15 upregulated and 54 downregulated genes. In addition, GO and KEGG pathway analyses showed that BAFF-R deficiency played an important role in the primary and adaptive immune responses in PCP. Furthermore, the lncRNA and mRNA co-expression network was established. We noted that the core network of lncRNA-TF (transcription factor) pairs could be classified into the categories including infection and immunity pathways.ConclusionIn summary, in this study, we further explored the role of mature B cells in the pathogenesis and progression of PCP and the data demonstrated that BAFF-R deficiency could play a significant role in immune regulation in the PCP population.

A novel endoplasmic reticulum stress-related lncRNA prognostic risk model for cutaneous melanoma.

Endoplasmic reticulum stress (ERS) and long non-coding RNAs (lncRNAs) are important in melanoma development and progression. This study aimed to explore the prognostic value of ERS-associated lncRNA profiles in cutaneous melanoma (CM). The Cancer Genome Atlas (TCGA) provides the raw data of CM. GSEA website was used to obtain ERS-related genes, and mRNA and LncRNA co-expression network were used to obtain ERS-related lncRNAs. A Lasso regression analysis was used to identify a prognostic risk model for the composition of ERS-related lncRNAs. Patients were divided into high- and low-risk groups based on the model's risk score. The researchers then compared the two groups' survival rates, immune infiltration, chemotherapeutic drug sensitivity, and immune checkpoint gene expression. Thirty-nine ERS-related lncRNAs were discovered to be prognostic. A prognostic risk model made up of ten ERS-related lncRNAs was discovered. Patients in the low-risk group had a better prognosis than those in the high-risk group. An examination of tumor microenvironment revealed that risk scores correlated with immune cell infiltration in eight cases. Dacarbazine, paclitaxel, and cisplatin, three chemotherapy drugs, were more sensitive in the low-risk group than in the high-risk group. This study identified a risk model of ten ERS-related lncRNAs that have significant prognostic value in CM and could help guide clinical treatment.