Expression correlation attenuates within and between key signaling pathways in chronic kidney disease

Nicotine exacerbates diabetic nephropathy through upregulation of Grem1 expression

BACKGROUND

Diabetic nephropathy (DN) is a major complication of diabetes mellitus. Clinical reports indicate that smoking is a significant risk factor for chronic kidney disease, and the tobacco epidemic exacerbates kidney damage in patients with DN. However, the underlying molecular mechanisms remain unclear.

METHOD

In the present study, we used a diabetic mouse model to investigate the molecular mechanisms for nicotine-exacerbated DN. Twelve-week-old female mice were injected with streptozotocin (STZ) to establish a hyperglycemic diabetic model. After four months, the control and hyperglycemic diabetic mice were further divided into four groups (control, nicotine, diabetic mellitus, nicotine + diabetic mellitus) by intraperitoneal injection of nicotine or PBS. After two months, urine and blood were collected for kidney injury assay, and renal tissues were harvested for further molecular assays using RNA-seq analysis, real-time PCR, Western blot, and immunohistochemistry. In vitro studies, we used siRNA to suppress Grem1 expression in human podocytes. Then we treated them with nicotine and high glucose to compare podocyte injury.

RESULT

Nicotine administration alone did not cause apparent kidney injury, but it significantly increased hyperglycemia-induced albuminuria, BUN, plasma creatinine, and the kidney tissue mRNA expression of KIM-1 and NGAL. Results from RNA-seq analysis, real-time PCR, Western blot, and immunohistochemistry analysis revealed that, compared to hyperglycemia or nicotine alone, the combination of nicotine treatment and hyperglycemia significantly increased the expression of Grem1 and worsened DN. In vitro experiments, suppression of Grem1 expression attenuated nicotine-exacerbated podocyte injury.

CONCLUSION

Grem1 plays a vital role in nicotine-exacerbated DN. Grem1 may be a potential therapeutic target for chronic smokers with DN.

Novel targets in renal fibrosis based on bioinformatic analysis

Background: Renal fibrosis is a widely used pathological indicator of progressive chronic kidney disease (CKD), and renal fibrosis mediates most progressive renal diseases as a final pathway. Nevertheless, the key genes related to the host response are still unclear. In this study, the potential gene network, signaling pathways, and key genes under unilateral ureteral obstruction (UUO) model in mouse kidneys were investigated by integrating two transcriptional data profiles. Methods: The mice were exposed to UUO surgery in two independent experiments. After 7 days, two datasets were sequenced from mice kidney tissues, respectively, and the transcriptome data were analyzed to identify the differentially expressed genes (DEGs). Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were executed. A Protein-Protein Interaction (PPI) network was constructed based on an online database STRING. Additionally, hub genes were identified and shown, and their expression levels were investigated in a public dataset and confirmed by quantitative real time-PCR (qRT-PCR) in vivo. Results: A total of 537 DEGs were shared by the two datasets. GO and the KEGG analysis showed that DEGs were typically enriched in seven pathways. Specifically, five hub genes (Bmp1, CD74, Fcer1g, Icam1, H2-Eb1) were identified by performing the 12 scoring methods in cytoHubba, and the receiver operating characteristic (ROC) curve indicated that the hub genes could be served as biomarkers. Conclusion: A gene network reflecting the transcriptome signature in CKD was established. The five hub genes identified in this study are potentially useful for the treatment and/or diagnosis CKD as biomarkers.

A Novel Strategy to Identify Prognosis-Relevant Gene Sets in Cancers

Molecular prognosis markers hold promise for improved prediction of patient survival, and a pathway or gene set may add mechanistic interpretation to their prognostic prediction power. In this study, we demonstrated a novel strategy to identify prognosis-relevant gene sets in cancers. Our study consists of a first round of gene-level analyses and a second round of gene-set-level analyses, in which the Composite Gene Expression Score critically summarizes a surrogate expression value at gene set level and a permutation procedure is exerted to assess prognostic significance of gene sets. An optional differential coexpression module is appended to the two phases of survival analyses to corroborate and refine prognostic gene sets. Our strategy was demonstrated in 33 cancer types across 32,234 gene sets. We found oncogenic gene sets accounted for an increased proportion among the final gene sets, and genes involved in DNA replication and DNA repair have ubiquitous prognositic value for multiple cancer types. In summary, we carried out the largest gene set based prognosis study to date. Compared to previous similar studies, our approach offered multiple improvements in design and methodology implementation. Functionally relevant gene sets of ubiquitous prognostic significance in multiple cancer types were identified.

Conditional transcriptional relationships may serve as cancer prognostic markers

BACKGROUND

While most differential coexpression (DC) methods are bound to quantify a single correlation value for a gene pair across multiple samples, a newly devised approach under the name Correlation by Individual Level Product (CILP) revolutionarily projects the summary correlation value to individual product correlation values for separate samples. CILP greatly widened DC analysis opportunities by allowing integration of non-compromised statistical methods.

METHODS

Here, we performed a study to verify our hypothesis that conditional relationships, i.e., gene pairs of remarkable differential coexpression, may be sought as quantitative prognostic markers for human cancers. Alongside the seeking of prognostic gene links in a pan-cancer setting, we also examined whether a trend of global expression correlation loss appeared in a wide panel of cancer types and revisited the controversial subject of mutual relationship between the DE approach and the DC approach.

RESULTS

By integrating CILP with classical univariate survival analysis, we identified up to 244 conditional gene links as potential prognostic markers in five cancer types. In particular, five prognostic gene links for kidney renal papillary cell carcinoma tended to condense around cancer gene ESPL1, and the transcriptional synchrony between ESPL1 and PTTG1 tended to be elevated in patients of adverse prognosis. In addition, we extended the observation of global trend of correlation loss in more than ten cancer types and empirically proved DC analysis results were independent of gene differential expression in five cancer types.

CONCLUSIONS

Combining the power of CILP and the classical survival analysis, we successfully fetched conditional transcriptional relationships that conferred prognosis power for five cancer types. Despite a general trend of global correlation loss in tumor transcriptomes, most of these prognosis conditional links demonstrated stronger expression correlation in tumors, and their stronger coexpression was associated with poor survival.

RNA editing affects cis-regulatory elements and predicts adverse cancer survival

BACKGROUND

RNA editing exerts critical impacts on numerous biological processes and thus are implicated in crucial human phenotypes, including tumorigenesis and prognosis. While previous studies have analyzed aggregate RNA editing activity at the sample level and associated it with overall cancer survival, there is not yet a large-scale disease-specific survival study to examine genome-wide RNA editing sites' prognostic value taking into account the host gene expression and clinical variables.

METHODS

In this study, we solved comprehensive Cox proportional models of disease-specific survival on individual RNA-editing sites plus host gene expression and critical demographic covariates. This allowed us to interrogate the prognostic value of a large number of RNA-editing sites at single-nucleotide resolution.

RESULTS

As a result, we identified 402 gene-proximal RNA-editing sites that generally predict adverse cancer survival. For example, an RNA-editing site residing in ZNF264 indicates poor survival of uterine corpus endometrial carcinoma, with a hazard ratio of 2.13 and an adjusted p-value of 4.07 × 10-7 . Some of these prognostic RNA-editing sites mediate the binding of RNA binding proteins and microRNAs, thus propagating their impacts to extensive regulatory targets.

CONCLUSIONS

In conclusion, RNA editing affects cis-regulatory elements and predicts adverse cancer survival.

MetaGSCA: A tool for meta-analysis of gene set differential coexpression

Analyses of gene set differential coexpression may shed light on molecular mechanisms underlying phenotypes and diseases. However, differential coexpression analyses of conceptually similar individual studies are often inconsistent and underpowered to provide definitive results. Researchers can greatly benefit from an open-source application facilitating the aggregation of evidence of differential coexpression across studies and the estimation of more robust common effects. We developed Meta Gene Set Coexpression Analysis (MetaGSCA), an analytical tool to systematically assess differential coexpression of an a priori defined gene set by aggregating evidence across studies to provide a definitive result. In the kernel, a nonparametric approach that accounts for the gene-gene correlation structure is used to test whether the gene set is differentially coexpressed between two comparative conditions, from which a permutation test p-statistic is computed for each individual study. A meta-analysis is then performed to combine individual study results with one of two options: a random-intercept logistic regression model or the inverse variance method. We demonstrated MetaGSCA in case studies investigating two human diseases and identified pathways highly relevant to each disease across studies. We further applied MetaGSCA in a pan-cancer analysis with hundreds of major cellular pathways in 11 cancer types. The results indicated that a majority of the pathways identified were dysregulated in the pan-cancer scenario, many of which have been previously reported in the cancer literature. Our analysis with randomly generated gene sets showed excellent specificity, indicating that the significant pathways/gene sets identified by MetaGSCA are unlikely false positives. MetaGSCA is a user-friendly tool implemented in both forms of a Web-based application and an R package "MetaGSCA". It enables comprehensive meta-analyses of gene set differential coexpression data, with an optional module of post hoc pathway crosstalk network analysis to identify and visualize pathways having similar coexpression profiles.

Rewired Pathways and Disrupted Pathway Crosstalk in Schizophrenia Transcriptomes by Multiple Differential Coexpression Methods

Transcriptomic studies of mental disorders using the human brain tissues have been limited, and gene expression signatures in schizophrenia (SCZ) remain elusive. In this study, we applied three differential co-expression methods to analyze five transcriptomic datasets (three RNA-Seq and two microarray datasets) derived from SCZ and matched normal postmortem brain samples. We aimed to uncover biological pathways where internal correlation structure was rewired or inter-coordination was disrupted in SCZ. In total, we identified 60 rewired pathways, many of which were related to neurotransmitter, synapse, immune, and cell adhesion. We found the hub genes, which were on the center of rewired pathways, were highly mutually consistent among the five datasets. The combinatory list of 92 hub genes was generally multi-functional, suggesting their complex and dynamic roles in SCZ pathophysiology. In our constructed pathway crosstalk network, we found "Clostridium neurotoxicity" and "signaling events mediated by focal adhesion kinase" had the highest interactions. We further identified disconnected gene links underlying the disrupted pathway crosstalk. Among them, four gene pairs (PAK1:SYT1, PAK1:RFC5, DCTN1:STX1A, and GRIA1:MAP2K4) were normally correlated in universal contexts. In summary, we systematically identified rewired pathways, disrupted pathway crosstalk circuits, and critical genes and gene links in schizophrenia transcriptomes.

The International Conference on Intelligent Biology and Medicine 2019 (ICIBM 2019): computational methods and applications in medical genomics

In this editorial, we briefly summarized the International Conference on Intelligent Biology and Medicine 2019 (ICIBM 2019) that was held on June 9-11, 2019 at Columbus, Ohio, USA. We further introduced the 19 research articles included in this supplement issue, covering four major areas, namely computational method development, genomics analysis, network-based analysis and biomarker prediction. The selected papers perform cutting edge computational research applied to a broad range of human diseases such as cancer, neural degenerative and chronic inflammatory disease. They also proposed solutions for fundamental medical genomics problems range from basic data processing and quality control to functional interpretation, biomarker and drug prediction, and database releasing.