A network-based variable selection approach for identification of modules and biomarker genes associated with end-stage kidney disease

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Urine peptidome in combination with transcriptomics analysis highlights MMP7, MMP14 and PCSK5 for further investigation in chronic kidney disease

Chronic kidney disease (CKD) is characterized by the loss of kidney function. The molecular mechanisms underlying the development and progression of CKD are still not fully understood. Among others, the urinary peptidome has been extensively studied, with several urinary peptides effectively detecting disease progression. However, their link to proteolytic events has not been made yet. This study aimed to predict the proteases involved in the generation of CKD-associated urinary excreted peptides in a well-matched (for age, sex, lack of heart disease) case-control study. The urinary peptide profiles from CKD (n = 241) and controls (n = 240) were compared and statistically analyzed. The in-silico analysis of the involved proteases was performed using Proteasix and proteases activity was predicted based on the abundance changes of the associated peptides. Predictions were cross-correlated to transcriptomics datasets by using the Nephroseq database. Information on the respective protease inhibitors was also retrieved from the MEROPS database. Totally, 303 urinary peptides were significantly associated with CKD. Among the most frequently observed were fragments of collagen types I, II and III, uromodulin, albumin and beta-2-microglobulin. Proteasix predicted 16 proteases involved in their generation. Through investigating CKD-associated transcriptomics datasets, several proteases are highlighted including members of matrix metalloproteinases (MMP7, MMP14) and serine proteases (PCSK5); laying the foundation for further studies towards elucidating their role in CKD pathophysiology.

Urine peptidome analysis in cardiorenal syndrome reflects molecular processes

The cardiorenal syndrome (CRS) is defined as the confluence of heart-kidney dysfunction. This study investigates the molecular differences at the level of the urinary peptidome between CRS patients and controls and their association to disease pathophysiology. The urinary peptidome of CRS patients (n = 353) was matched for age and sex with controls (n = 356) at a 1:1 ratio. Changes in the CRS peptidome versus controls were identified after applying the Mann-Whitney test, followed by correction for multiple testing. Proteasix tool was applied to investigate predicted proteases involved in CRS-associated peptide generation. Overall, 559 differentially excreted urinary peptides were associated with CRS patients. Of these, 193 peptides were specifically found in CRS when comparing with heart failure and chronic kidney disease urinary peptide profiles. Proteasix predicted 18 proteases involved in > 1% of proteolytic cleavage events including multiple forms of MMPs, proprotein convertases, cathepsins and kallikrein 4. Forty-four percent of the cleavage events were produced by 3 proteases including MMP13, MMP9 and MMP2. Pathway enrichment analysis supported that ECM-related pathways, fibrosis and inflammation were represented. Collectively, our study describes the changes in urinary peptides of CRS patients and potential proteases involved in their generation, laying the basis for further validation.

Identification of the hub genes in gastric cancer through weighted gene co-expression network analysis

Background

Gastric cancer is one of the most lethal tumors and is characterized by poor prognosis and lack of effective diagnostic or therapeutic biomarkers. The aim of this study was to find hub genes serving as biomarkers in gastric cancer diagnosis and therapy.

Methods

GSE66229 from Gene Expression Omnibus (GEO) was used as training set. Genes bearing the top 25% standard deviations among all the samples in training set were performed to systematic weighted gene co-expression network analysis (WGCNA) to find candidate genes. Then, hub genes were further screened by using the “least absolute shrinkage and selection operator” (LASSO) logistic regression. Finally, hub genes were validated in the GSE54129 dataset from GEO by supervised learning method artificial neural network (ANN) algorithm.

Results

Twelve modules with strong preservation were identified by using WGCNA methods in training set. Of which, five modules significantly related to gastric cancer were selected as clinically significant modules, and 713 candidate genes were identified from these five modules. Then, ADIPOQ, ARHGAP39, ATAD3A, C1orf95, CWH43, GRIK3, INHBA, RDH12, SCNN1G, SIGLEC11 and LYVE1 were screened as the hub genes. These hub genes successfully differentiated the tumor samples from the healthy tissues in an independent testing set through artificial neural network algorithm with the area under the receiver operating characteristic curve at 0.946.

Conclusions

These hub genes bearing diagnostic and therapeutic values, and our results may provide a novel prospect for the diagnosis and treatment of gastric cancer in the future.