Co-expression network of long non-coding RNA and mRNA reveals molecular phenotype changes in kidney development of prenatal chlorpyrifos exposure in a mouse model

Gene network landscape of mouse splenocytes reveals integrin complex as the A151 ODN-responsive hub molecule in the immune transcriptome

Homeostatic restoration of an inflammatory response requires quenching of the immune system after pathogen threats vanish. A continued assault orchestrated by host defense results in tissue destruction or autoimmunity. A151 is the epitome of synthetic oligodeoxynucleotides (ODNs) that curb the immune response by a subset of white corpuscles through repetitive telomere-derived TTAGGG sequences. Currently, the genuine effect of A151 on the immune cell transcriptome remains unknown. Here, we leveraged an integrative approach where weighted gene co-expression network analysis (WGCNA), differential gene expression analysis, and gene set enrichment analysis (GSEA) of our in-house microarray datasets aided our understanding of how A151 ODN suppresses the immune response in mouse splenocytes. Our bioinformatics results, together with experimental validations, indicated that A151 ODN acts on components of integrin complexes, Itgam and Itga6, to interfere with immune cell adhesion and thereby suppresses the immune response in mice. Moreover, independent lines of evidence in this work converged on the observation that cell adhesion by integrin complexes serves as a focal point for cellular response to A151 ODN treatment in immune cells. Taken together, the outcome of this study sheds light on the molecular basis of immune suppression by a clinically useful DNA-based therapeutic agent.

Maternal NO2 exposure disturbs the long noncoding RNA expression profile in the lungs of offspring in time-series patterns

Gestation is a sensitive window to nitrogen dioxide (NO₂) exposure, which may disturb fetal lung development and lung function later in life. Animal and epidemiological studies indicated that long noncoding RNAs (lncRNAs) participate in abnormal lung development induced by environmental pollutant exposure. In the present study, pregnant C57BL/6J mice were exposed to 2.5 ppm NO₂ (mimicking indoor occupational exposure) or clean air, and lncRNAs expression profiles in the lungs of offspring mice were determined by lncRNA-seq on embryonic day 13.5 (E13.5), E18.5, postnatal day 1 (P1), and P14. The lung histopathology examination of offspring was performed, followed by weighted gene coexpression network analysis (WGCNA), prediction of lncRNAs-target genes, and the biological processes enrichment analysis of lncRNAs. Our results indicated that maternal NO₂ exposure induced hypoalveolarization on P14 and differentially expressed lncRNAs showed a time-series pattern. Following WGCNA and enrichment analysis, 2 modules participated in development-related pathways. Importantly, the expressions of related genes were altered, some of which were confirmed to be related to abnormal vascular development and even lung diseases. The research points out that the maternal NO₂ exposure leads to abnormal lung development in offspring that might be related to altered lncRNAs expression profiles with time-series-pattern.

Identification of crucial lncRNAs and mRNAs in liver regeneration after portal vein ligation through weighted gene correlation network analysis

BACKGROUND

Portal vein ligation (PVL)-induced liver hypertrophy increases future liver remnant (FLR) volume and improves resectability of large hepatic carcinoma. However, the molecular mechanism by which PVL facilitates liver hypertrophy remains poorly understood.

METHODS

To gain mechanistic insight, we established a rat PVL model and carried out a comprehensive transcriptome analyses of hepatic lobes preserving portal blood supply at 0, 1, 7, and 14-day after PVL. The differentially expressed (DE) long-non coding RNAs (lncRNAs) and mRNAs were applied to conduct weighted gene co-expression network analysis (WGCNA). LncRNA-mRNA co-expression network was constructed in the most significant module. The modules and genes associated with PVL-induced liver hypertrophy were assessed through quantitative real-time PCR.

RESULTS

A total of 4213 DElncRNAs and 6809 DEmRNAs probesets, identified by transcriptome analyses, were used to carry out WGCNA, by which 10 modules were generated. The largest and most significant module (marked in black_M6) was selected for further analysis. Gene Ontology (GO) analysis of the module exhibited several key biological processes associated with liver regeneration such as complement activation, IL-6 production, Wnt signaling pathway, autophagy, etc. Sixteen mRNAs (Notch1, Grb2, IL-4, Cops4, Stxbp1, Khdrbs2, Hdac2, Gnb3, Gng10, Tlr2, Sod1, Gosr2, Rbbp5, Map3k3, Golga2, and Rev3l) and ten lncRNAs (BC092620, AB190508, EF076772, BC088302, BC158675, BC100646, BC089934, L20987, BC091187, and M23890) were identified as hub genes in accordance with gene significance value, module membership value, protein-protein interaction (PPI) and lncRNA-mRNA co-expression network. Furthermore, the overexpression of 3 mRNAs (Notch1, Grb2 and IL-4) and 4 lncRNAs (BC089934, EF076772, BC092620, and BC088302) was validated in hypertrophic liver lobe tissues from PVL rats and patients undergoing hepatectomy after portal vein embolization (PVE).

CONCLUSIONS

Microarray and WGCNA analysis revealed that the 3 mRNAs (Notch1, Grb2 and IL-4) and the 4 lncRNAs (BC089934, EF076772, BC092620 and BC088302) may be promising targets for accelerating liver regeneration before extensive hepatectomy.

Identification of Novel Prognostic Markers Associated With Laryngeal Squamous Cell Carcinoma Using Comprehensive Analysis

Background

Laryngeal squamous cell carcinoma (LSCC) is a leading malignant cancer of the head and neck. Patients with LSCC, in which the cancer has infiltrated and metastasized, have a poor prognosis. Therefore, there is an urgent need to identify more potential targets for drugs and biomarkers for early diagnosis.

Methods

RNA sequence data from LSCC and patients’ clinical traits were obtained from the Gene Expression Omnibus (GEO) (GSE142083) and The Cancer Genome Atlas (TCGA) database. Differentially expressed gene (DEG) analysis and weighted gene co-expression network analysis (WGCNA) were performed to identify hub genes. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, prognostic value analysis, receiver operating characteristic (ROC) curve analysis, gene mutation analysis, tumor-infiltrating immune cell abundance profile estimation, gene set variation analysis (GSVA), and gene set enrichment analysis (GSEA) were performed. Single-gene RNA sequencing data were obtained from the GSE150321 dataset. Cell proliferation and viability were confirmed by the CCK-8 assay and real-time PCR.

Results

A total of 701 DEGs, including 329 upregulated and 372 downregulated genes, were screened in the GSE142083 dataset. Using WGCNA, three modules were identified to be closely related to LSCC. After intersecting the DEGs and performing univariate and multivariate Cox analyses, a novel prognostic model based on three genes (SLC35C1, HOXB7, and TEDC2) for LSCC was established. Interfering TEDC2 expression inhibited tumor cell proliferation and migration.

Conclusions

Our results show that SLC35C1, HOXB7, and TEDC2 have the potential to become new therapeutic targets and prognostic biomarkers for LSCC.

FNDC3B and BPGM Are Involved in Human Papillomavirus-Mediated Carcinogenesis of Cervical Cancer

Human papillomavirus (HPV)-mediated cervical carcinogenesis is a multistep progressing from persistent infection, precancerous lesion to cervical cancer (CCa). Although molecular alterations driven by viral oncoproteins are necessary in cervical carcinogenesis, the key regulators behind the multistep process remain not well understood. It is pivotal to identify the key genes involved in the process for early diagnosis and treatment of this disease. Here we analyzed the mRNA expression profiles in cervical samples including normal, cervical intraepithelial neoplasia (CIN), and CCa. A co-expression network was constructed using weighted gene co-expression network analysis (WGCNA) to reveal the crucial modules in the dynamic process from HPV infection to CCa development. Furthermore, the differentially expressed genes (DEGs) that could distinguish all stages of progression of CCa were screened. The key genes involved in HPV-CCa were identified. It was found that the genes involved in DNA replication/repair and cell cycle were upregulated in CIN compared with normal control, and sustained in CCa, accompanied by substantial metabolic shifts. We found that upregulated fibronectin type III domain-containing 3B (FNDC3B) and downregulated bisphosphoglycerate mutase (BPGM) could differentiate all stages of CCa progression. In patients with CCa, a higher expression of FNDC3B or lower expression of BPGM was closely correlated with a shorter overall survival (OS) and disease-free survival (DFS). A receiver operating characteristic (ROC) analysis of CIN and CCa showed that FNDC3B had the highest sensitivity and specificity for predicting CCa development. Taken together, the current data showed that FNDC3B and BPGM were key genes involved in HPV-mediated transformation from normal epithelium to precancerous lesions and CCa.

The cytoskeleton as a non-cholinergic target of organophosphate compounds

Current organophosphate (OP) toxicity research now considers potential non-cholinergic mechanisms for these compounds, since the inhibition of acetylcholinesterase (AChE) cannot completely explain all the adverse biological effects of OP. Thanks to the development of new strategies for OP detection, some potential molecular targets have been identified. Among these molecules are several cytoskeletal proteins, including actin, tubulin, intermediate filament proteins, and associated proteins, such as motor proteins, microtubule-associated proteins (MAPs), and cofilin. in vitro, ex vivo, and some in vivo reports have identified alterations in the cytoskeleton following OP exposure, including cell morphology defects, cells detachments, intracellular transport disruption, aberrant mitotic spindle formation, modification of cell motility, and reduced phagocytic capability, which implicate the cytoskeleton in OP toxicity. Here, we reviewed the evidence indicating the cytoskeletal targets of OP compounds, including their strategies, the potential effects of their alterations, and their possible participation in neurotoxicity, embryonic development, cell division, and immunotoxicity related to OP compounds exposure.