Systematic Analysis of mRNAs and ncRNAs in BMSCs of Senile Osteoporosis Patients

Senile osteoporosis (SOP) is a worldwide age-related disease characterized by the loss of bone mass and decrease in bone strength. Bone mesenchymal stem cells (BMSCs) play an important role in the pathology of senile osteoporosis. Abnormal expression and regulation of non-coding RNA (ncRNA) are involved in a variety of human diseases. In the present study, we aimed to identify differentially expressed mRNAs and ncRNAs in senile osteoporosis patient-derived BMSCs via high-throughput transcriptome sequencing in combination with bioinformatics analysis. As a result, 415 mRNAs, 30 lncRNAs, 6 circRNAs and 27 miRNAs were found to be significantly changed in the senile osteoporosis group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were applied to analyze the function of differentially expressed mRNAs and ncRNAs. The circRNA–miRNA–mRNA regulatory network was constructed using the cytoHubba plugin based on the Cytoscape software. Interestingly, circRNA008876-miR-150-5p-mRNA was the sole predicted circRNA-miRNA-mRNA network. The differential expression profile of this ceRNA network was further verified by qRT-PCR. The biological function of this network was validated by overexpression and knockdown experiments. In conclusion, circRNA008876-miR-150-5p-mRNA could be an important ceRNA network involved in senile osteoporosis, which provides potential biomarkers and therapeutic targets for senile osteoporosis.

Glutathione, carbohydrate and other metabolites of Larix olgensis A. Henry reponse to polyethylene glycol-simulated drought stress

Drought stress in trees limits their growth, survival, and productivity and it negatively affects the afforestation survival rate. Our study focused on the molecular responses to drought stress in a coniferous species Larix olgensis A. Henry. Drought stress was simulated in one-year-old seedlings using 25% polyethylene glycol 6000. The drought stress response in these seedlings was assessed by analyzing select biochemical parameters, along with gene expression and metabolite profiles. The soluble protein content, peroxidase activity, and malondialdehyde content of L. olgensis were significantly changed during drought stress. Quantitative gene expression analysis identified a total of 8172 differentially expressed genes in seedlings processed after 24 h, 48 h, and 96 h of drought stress treatment. Compared with the gene expression profile of the untreated control, the number of up-regulated genes was higher than that of down-regulated genes, indicating that L. olgensis mainly responded to drought stress through positive regulation. Metabolite analysis of the control and stress-treated samples showed that under drought stress, the increased abundance of linoleic acid was the highest among up-regulated metabolites, which also included some saccharides. A combined analysis of the transcriptome and metabolome revealed that genes dominating the differential expression profile were involved in glutathione metabolism, galactose metabolism, and starch and sucrose metabolism. Moreover, the relative abundance of specific metabolites of these pathways was also altered. Thus, our results indicated that L. olgensis prevented free radical-induced damage through glutathione metabolism and responded to drought through sugar accumulation.

Functional analysis of co-expression networks of zebrafish ace2 reveals enrichment of pathways associated with development and disease.

Human Angiotensin I Converting Enzyme 2 (ACE2) plays essential roles in blood pressure regulation and SARS-CoV-2 entry. ACE2 has a highly conserved, one-to-one ortholog, called ace2, in zebrafish, which is an important model for human diseases. However, zebrafish ace2 expression profile has not yet been studied during early development, between genders, across different genotypes, or in disease. Moreover, a network-based meta-analysis for the extraction of functionally enriched pathways associated with differential ace2 expression is lacking in the literature. Herein, we first identified significant development-, tissue-, genotype- and gender-specific modulations in ace2 expression via meta-analysis of zebrafish Affymetrix transcriptomics datasets (ndatasets=107); and the correlation analysis of ace2 meta-differential expression profile revealed distinct positively and negatively correlated local functionally enriched gene networks. Moreover, we demonstrated that ace2 expression was significantly modulated under different physiological and pathological conditions related to development, tissue, gender, diet, infection, and inflammation using additional RNAseq datasets. Our findings implicate a novel translational role for zebrafish ace2 in organ differentiation and pathologies observed in the intestines and liver.

A novel insight into differential expression profiles of sporadic cerebral cavernous malformation patients with different symptoms

Cerebral cavernous malformation (CCM) is a vascular lesion of the central nervous system that may lead to distinct symptoms among patients including cerebral hemorrhages, epileptic seizures, focal neurologic deficits, and/or headaches. Disease-related mutations were identified previously in one of the three CCM genes: CCM1, CCM2, and CCM3. However, the rate of these mutations in sporadic cases is relatively low, and new studies report that mutations in CCM genes may not be sufficient to initiate the lesions. Despite the growing body of research on CCM, the underlying molecular mechanism has remained largely elusive. In order to provide a novel insight considering the specific manifested symptoms, CCM patients were classified into two groups (as Epilepsy and Hemorrhage). Since the studied patients experience various symptoms, we hypothesized that the underlying cause for the disease may also differ between those groups. To this end, the respective transcriptomes were compared to the transcriptomes of the control brain tissues and among each other. This resulted into the identification of the differentially expressed coding genes and the delineation of the corresponding differential expression profile for each comparison. Notably, some of those differentially expressed genes were previously implicated in epilepsy, cell structure formation, and cell metabolism. However, no CCM1-3 gene deregulation was detected. Interestingly, we observed that when compared to the normal controls, the expression of some identified genes was only significantly altered either in Epilepsy (EGLN1, ELAVL4, and NFE2l2) or Hemorrhage (USP22, EYA1, SIX1, OAS3, SRMS) groups. To the best of our knowledge, this is the first such effort focusing on CCM patients with epileptic and hemorrhagic symptoms with the purpose of uncovering the potential CCM-related genes. It is also the first report that presents a gene expression dataset on Turkish CCM patients. The results suggest that the new candidate genes should be explored to further elucidate the CCM pathology. Overall, this work constitutes a step towards the identification of novel potential genetic targets for the development of possible future therapies.

Differential Expression Profile of lncRNA in Glioma Cells and the Effect of lncRNA NKX3-1 on Glioma Cells Through Fem1b/SPDEF Pathway

ObjectiveTo investigate the differential expression of lncRNA in glioma cells, as well as the effect of lncRNA NKX3-1 on glioma cells.MethodsGlioma-related data were first downloaded from the TCGA database and analyzed using bioinformatics, after which the lncRNA NKX3-1 was chosen for further experiments. The expression of the lncRNA NKX3-1 in glioma tumor samples was detected using qRT-PCR. The subcellular localization of lncRNA NKX3-1 was determined using fluorescence in situ hybridization (FISH). CCK-8, flow cytometry, cell scratch, and transwell assays were used to detect cell proliferation, apoptosis, and invasion. The downstream pathway of lncRNA NKX3-1 was investigated using luciferase assays and detected using western blot, transwell, and cell scratch assays.ResultsThe differential expression profile of lncRNA in glioma was obtained. NKX3-1 lncRNA was found to be significantly increased in glioma tumor tissues. LncRNA NKX3-1 was found in the nucleus. Proliferation, invasion, and migration of glioma cells were significantly increased (P <0.05) in the lncRNA NKX3-1 overexpression group, while apoptosis ability was significantly decreased (P <0.05). Tumor volume and weight were significantly increased in the lncRNA NKX3-1 overexpression group in nude mice (P <0.05). LncRNA NKX3-1 significantly increased the luciferase activity of Fem1b 3’-UTR-WT reporter genes (P <0.05) as well as the levels of SPDEF protein (P <0.05). The protein level of FEM1B was significantly reduced. Cell invasion and migration were significantly increased (P <0.05) in the lncRNA NKX3-1 overexpression group plus SPDEF group.ConclusionWe investigated the differential expression profile of lncRNAs in glioma and discovered that the lncRNA NKX3-1 plays an important role in cancer promotion via the Fem1b/SPDEF pathway.

OMRT-2. Liquid biopsy for patient stratification and monitoring of dacomitinib clinical trial in patients with EGFR amplified recurrent glioblastoma

Introduction Liquid biopsy for the detection and monitoring of brain tumors is of significant clinical interest. The ability to non-invasively profile tumors can avoid a risky biopsy and opens avenues for testing novel therapies by accurately stratifying patients to receive the right therapy. Here, we provide evidence of EV RNA-based diagnosis, patient stratification, and assessment of response to therapy in the setting of a clinical trial evaluating the efficacy of dacomitinib, an EGFR tyrosine kinase inhibitor in patients with recurrent, EGFR amplified GBM(NCT01112527). Methods We performed RNASeq on long RNA extracted from the serum samples, pre-treatment and 1-month post-treatment. Results Firstly, longRNASeq allowed the detection of thousands of mRNA, lincRNAs and antisense RNAs enabling the study of a wider repertoire of potential RNA based biomarkers. Secondly, we observed a differential expression profile in serum EV RNA of GBM patients and healthy controls. Combining our findings with TCGA data and literature screening, we generated a 25 gene signature representative of critical pathways in several hallmarks of cancer. Thirdly, we observed a differential expression profile in serum EV RNA of responders to dacomitinib compared to non-responders in pre-treatment serum. Specifically, the EV mRNAs ZNF35 and LAMTOR2 distinguish responders from non-responders (p-adjusted = 2.6E-8 and 2.4E-6, respectively) allowing potential patient stratification. Finally, we observed a differential expression profile in serum EV RNA of responders to dacomitinib compared to non-responders in post-treatment serum. EV mRNA DNMT3A is significantly enriched (p-adjusted = 1.8E-4) in post-treatment serum of responders compared to non-responders to dacomitinib allowing potential monitoring of response to therapy. Conclusion This study represents the first longitudinal profiling of the EV transcriptome in a cohort of genomically selected GBM patients. These findings are a tantalizing step toward liquid biopsy-based biomarkers for the detection of GBM, as well as patient stratification and monitoring.

Identifying prognostic lncRNAs based on a ceRNA regulatory network in laryngeal squamous cell carcinoma

Purpose Growing evidence demonstrates that long non-coding RNAs (lncRNAs) play a crucial role as competing endogenous RNAs (ceRNAs) in tumor occurrence. The lncRNAs’ functions and clinical significance in laryngeal squamous cell carcinoma (LSCC) remain unclear. The study aims to reveal the lncRNA-associated ceRNA regulatory network of LSCC and clarify its clinical relevance. Methods Here, we obtained LSCC transcriptome data from The Cancer Genome Atlas (TCGA) database and identified the differential expression profile of lncRNAs, miRNAs, and mRNAs by the EdgeR R package. The function enrichment analysis of mRNAs was performed using clusterProfiler R package and GSEA3.0. Then, we constructed a ceRNA network and prognosis model based on lncRNAs through bioinformatic methods. Moreover, we explored the functions of prognosis-related lncRNA in LSCC by CCK-8 and transwell assay. Results 1961 lncRNAs, 69 miRNAs, and 2224 mRNAs were identified as differentially expressed genes in LSCC tissues. According to the transcriptome differential expression profile, a ceRNA network containing 61 lncRNAs, 21 miRNAs, and 77 mRNAs was established. Then, four lncRNAs (AC011933.2, FAM30A, LINC02086, LINC02575) were identified from the ceRNA network to build a prognosis model for LSCC patients. And we found that LINC02086 and LINC02575 promoted the proliferation, migration, and invasion of LSCC cells while AC011933.2 and FAM30A inhibited these biological functions in vitro. Furthermore, we validated that LINc02086/miR-770-5p/SLC26A2 axis promoted migration in LSCC. Conclusion Four lncRNAs of the ceRNA network were abnormally expressed and related to patient prognosis in LSCC. They played a significant role in the progress of LSCC via affecting the proliferation and metastasis of tumor cells.

Identifying Prognostic Lncrnas Based on A ceRNA Regulatory Network in Laryngeal Squamous Cell Carcinoma

Purpose Growing evidence demonstrates that long non-coding RNAs (lncRNAs) play a crucial role as competing endogenous RNAs (ceRNAs) in tumor occurrence. The lncRNAs' function and clinical significance in laryngeal squamous cell carcinoma (LSCC) remain unclear. The study aims to reveal the lncRNA-associated ceRNA regulatory network in LSCC and clarify its clinical relevance.Methods Here, we obtained LSCC transcriptome sequencing data from The Cancer Genome Atlas (TCGA) database and identified the differential expression profile of lncRNAs, miRNAs, and mRNAs by the EdgeR R package. The function enrichment analysis of mRNAs was performed using clusterProfiler R package and GSEA3.0. Then, we constructed a ceRNA network and a prognosis model based on lncRNAs through bioinformatic methods. Moreover, we explored the functions of prognosis-related lncRNA in LSCC by CCK-8 and transwell assay. Results 1961 lncRNAs, 69 miRNAs, and 2224 mRNAs were identified as differentially expressed genes in LSCC tissues. According to the transcriptome differential expression profile, a ceRNA network containing 61 lncRNAs, 21 miRNAs, and 77 mRNAs was established. Then, four lncRNAs (AC011933.2, FAM30A, LINC02086, LINC02575) were extracted from the ceRNA network to build a prognosis model for LSCC patients. Furthermore, we found that LINC02086 and LINC02575 promoted the proliferation, migration, and invasion of LSCC cells while AC011933.2 and FAM30A inhibited these biological functions in vitro. Conclusion Four lncRNAs from the ceRNA network were abnormally expressed and related to patient prognosis in LSCC. They played a significant role in the progress of LSCC via affecting the proliferation and metastasis of tumor cells.