Interactive molecular causal networks of hypertension using a fast machine learning algorithm MRdualPC.

BACKGROUND: Understanding the complex interactions between genes and their causal effects on diseases is crucial for developing targeted treatments and gaining insight into biological mechanisms. However, the analysis of molecular networks, especially in the context of high-dimensional data, presents significant challenges. METHODS: This study introduces MRdualPC, a computationally tractable algorithm based on the MRPC approach, to infer large-scale causal molecular networks. We apply MRdualPC to investigate the upstream causal transcriptomics influencing hypertension using a comprehensive dataset of kidney genome and transcriptome data. RESULTS: Our algorithm proves to be 100 times faster than MRPC on average in identifying transcriptomics drivers of hypertension. Through clustering, we identify 63 modules with causal driver genes, including 17 modules with extensive causal networks. Notably, we find that genes within one of the causal networks are associated with the electron transport chain and oxidative phosphorylation, previously linked to hypertension. Moreover, the identified causal ancestor genes show an over-representation of blood pressure-related genes. CONCLUSIONS: MRdualPC has the potential for broader applications beyond gene expression data, including multi-omics integration. While there are limitations, such as the need for clustering in large gene expression datasets, our study represents a significant advancement in building causal molecular networks, offering researchers a valuable tool for analyzing big data and investigating complex diseases.

Kidney omics in hypertension: from statistical associations to biological mechanisms and clinical applications.

Hypertension is a major cardiovascular disease risk factor and contributor to premature death globally. Family-based investigations confirmed a significant heritable component of blood pressure (BP), whereas genome-wide association studies revealed >1000 common and rare genetic variants associated with BP and/or hypertension. The kidney is not only an organ of key relevance to BP regulation and the development of hypertension, but it also acts as the tissue mediator of genetic predisposition to hypertension. The identity of kidney genes, pathways, and related mechanisms underlying the genetic associations with BP has started to emerge through integration of genomics with kidney transcriptomics, epigenomics, and other omics as well as through applications of causal inference, such as Mendelian randomization. Single-cell methods further enabled mapping of BP-associated kidney genes to cell types, and in conjunction with other omics, started to illuminate the biological mechanisms underpinning associations of BP-associated genetic variants and kidney genes. Polygenic risk scores derived from genome-wide association studies and refined on kidney omics hold the promise of enhanced diagnostic prediction, whereas kidney omics-informed drug discovery is likely to contribute new therapeutic opportunities for hypertension and hypertension-mediated kidney damage.

Plasma Proteomics of Renal Function: A Transethnic Meta-Analysis and Mendelian Randomization Study.

BACKGROUND: Studies on the relationship between renal function and the human plasma proteome have identified several potential biomarkers. However, investigations have been conducted largely in European populations, and causality of the associations between plasma proteins and kidney function has never been addressed. METHODS: A cross-sectional study of 993 plasma proteins among 2882 participants in four studies of European and admixed ancestries (KORA, INTERVAL, HUNT, QMDiab) identified transethnic associations between eGFR/CKD and proteomic biomarkers. For the replicated associations, two-sample bidirectional Mendelian randomization (MR) was used to investigate potential causal relationships. Publicly available datasets and transcriptomic data from independent studies were used to examine the association between gene expression in kidney tissue and eGFR. RESULTS: In total, 57 plasma proteins were associated with eGFR, including one novel protein. Of these, 23 were additionally associated with CKD. The strongest inferred causal effect was the positive effect of eGFR on testican-2, in line with the known biological role of this protein and the expression of its protein-coding gene (SPOCK2) in renal tissue. We also observed suggestive evidence of an effect of melanoma inhibitory activity (MIA), carbonic anhydrase III, and cystatin-M on eGFR. CONCLUSIONS: In a discovery-replication setting, we identified 57 proteins transethnically associated with eGFR. The revealed causal relationships are an important stepping stone in establishing testican-2 as a clinically relevant physiological marker of kidney disease progression, and point to additional proteins warranting further investigation.

Hypertension and renin-angiotensin system blockers are not associated with expression of angiotensin-converting enzyme 2 (ACE2) in the kidney.

AIMS: Angiotensin-converting enzyme 2 (ACE2) is the cellular entry point for severe acute respiratory syndrome coronavirus (SARS-CoV-2)-the cause of coronavirus disease 2019 (COVID-19). However, the effect of renin-angiotensin system (RAS)-inhibition on ACE2 expression in human tissues of key relevance to blood pressure regulation and COVID-19 infection has not previously been reported. METHODS AND RESULTS: We examined how hypertension, its major metabolic co-phenotypes, and antihypertensive medications relate to ACE2 renal expression using information from up to 436 patients whose kidney transcriptomes were characterized by RNA-sequencing. We further validated some of the key observations in other human tissues and/or a controlled experimental model. Our data reveal increasing expression of ACE2 with age in both human lungs and the kidney. We show no association between renal expression of ACE2 and either hypertension or common types of RAS inhibiting drugs. We demonstrate that renal abundance of ACE2 is positively associated with a biochemical index of kidney function and show a strong enrichment for genes responsible for kidney health and disease in ACE2 co-expression analysis. CONCLUSION: Our results indicate that neither hypertension nor antihypertensive treatment is likely to alter the expression of the key entry receptor for SARS-CoV-2 in the human kidney. Our data further suggest that in the absence of SARS-CoV-2 infection, kidney ACE2 is most likely nephro-protective but the age-related increase in its expression within lungs and kidneys may be relevant to the risk of SARS-CoV-2 infection.

Human Y Chromosome Exerts Pleiotropic Effects on Susceptibility to Atherosclerosis.

OBJECTIVE: The male-specific region of the Y chromosome (MSY) remains one of the most unexplored regions of the genome. We sought to examine how the genetic variants of the MSY influence male susceptibility to coronary artery disease (CAD) and atherosclerosis. Approach and Results: Analysis of 129 133 men from UK Biobank revealed that only one of 7 common MSY haplogroups (haplogroup I1) was associated with CAD-carriers of haplogroup I1 had ≈11% increase in risk of CAD when compared with all other haplogroups combined (odds ratio, 1.11; 95% CI, 1.04-1.18; P=6.8×10-4). Targeted MSY sequencing uncovered 235 variants exclusive to this haplogroup. The haplogroup I1-specific variants showed 2.45- and 1.56-fold respective enrichment for promoter and enhancer chromatin states, in cells/tissues relevant to atherosclerosis, when compared with other MSY variants. Gene set enrichment analysis in CAD-relevant tissues showed that haplogroup I1 was associated with changes in pathways responsible for early and late stages of atherosclerosis development including defence against pathogens, immunity, oxidative phosphorylation, mitochondrial respiration, lipids, coagulation, and extracellular matrix remodeling. UTY was the only Y chromosome gene whose blood expression was associated with haplogroup I1. Experimental reduction of UTY expression in macrophages led to changes in expression of 59 pathways (28 of which overlapped with those associated with haplogroup I1) and a significant reduction in the immune costimulatory signal. CONCLUSIONS: Haplogroup I1 is enriched for regulatory chromatin variants in numerous cells of relevance to CAD and increases cardiovascular risk through proatherosclerotic reprogramming of the transcriptome, partly through UTY.

Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets.

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

Contributions of obesity to kidney health and disease: insights from Mendelian randomization and the human kidney transcriptomics.

AIMS: Obesity and kidney diseases are common complex disorders with an increasing clinical and economic impact on healthcare around the globe. Our objective was to examine if modifiable anthropometric obesity indices show putatively causal association with kidney health and disease and highlight biological mechanisms of potential relevance to the association between obesity and the kidney. METHODS AND RESULTS: We performed observational, one-sample, two-sample Mendelian randomization (MR) and multivariable MR studies in ∼300 000 participants of white-British ancestry from UK Biobank and participants of predominantly European ancestry from genome-wide association studies. The MR analyses revealed that increasing values of genetically predicted body mass index and waist circumference were causally associated with biochemical indices of renal function, kidney health index (a composite renal outcome derived from blood biochemistry, urine analysis, and International Classification of Disease-based kidney disease diagnoses), and both acute and chronic kidney diseases of different aetiologies including hypertensive renal disease and diabetic nephropathy. Approximately 13-16% and 21-26% of the potentially causal effect of obesity indices on kidney health were mediated by blood pressure and type 2 diabetes, respectively. A total of 61 pathways mapping primarily onto transcriptional/translational regulation, innate and adaptive immunity, and extracellular matrix and metabolism were associated with obesity measures in gene set enrichment analysis in up to 467 kidney transcriptomes. CONCLUSIONS: Our data show that a putatively causal association of obesity with renal health is largely independent of blood pressure and type 2 diabetes and uncover the signatures of obesity on the transcriptome of human kidney.

Uncovering genetic mechanisms of hypertension through multi-omic analysis of the kidney.

The kidney is an organ of key relevance to blood pressure (BP) regulation, hypertension and antihypertensive treatment. However, genetically mediated renal mechanisms underlying susceptibility to hypertension remain poorly understood. We integrated genotype, gene expression, alternative splicing and DNA methylation profiles of up to 430 human kidneys to characterize the effects of BP index variants from genome-wide association studies (GWASs) on renal transcriptome and epigenome. We uncovered kidney targets for 479 (58.3%) BP-GWAS variants and paired 49 BP-GWAS kidney genes with 210 licensed drugs. Our colocalization and Mendelian randomization analyses identified 179 unique kidney genes with evidence of putatively causal effects on BP. Through Mendelian randomization, we also uncovered effects of BP on renal outcomes commonly affecting patients with hypertension. Collectively, our studies identified genetic variants, kidney genes, molecular mechanisms and biological pathways of key relevance to the genetic regulation of BP and inherited susceptibility to hypertension.

Epigenome-wide association study of kidney function identifies trans-ethnic and ethnic-specific loci.

BACKGROUND: DNA methylation (DNAm) is associated with gene regulation and estimated glomerular filtration rate (eGFR), a measure of kidney function. Decreased eGFR is more common among US Hispanics and African Americans. The causes for this are poorly understood. We aimed to identify trans-ethnic and ethnic-specific differentially methylated positions (DMPs) associated with eGFR using an agnostic, genome-wide approach. METHODS: The study included up to 5428 participants from multi-ethnic studies for discovery and 8109 participants for replication. We tested the associations between whole blood DNAm and eGFR using beta values from Illumina 450K or EPIC arrays. Ethnicity-stratified analyses were performed using linear mixed models adjusting for age, sex, smoking, and study-specific and technical variables. Summary results were meta-analyzed within and across ethnicities. Findings were assessed using integrative epigenomics methods and pathway analyses. RESULTS: We identified 93 DMPs associated with eGFR at an FDR of 0.05 and replicated 13 and 1 DMPs across independent samples in trans-ethnic and African American meta-analyses, respectively. The study also validated 6 previously published DMPs. Identified DMPs showed significant overlap enrichment with DNase I hypersensitive sites in kidney tissue, sites associated with the expression of proximal genes, and transcription factor motifs and pathways associated with kidney tissue and kidney development. CONCLUSIONS: We uncovered trans-ethnic and ethnic-specific DMPs associated with eGFR, including DMPs enriched in regulatory elements in kidney tissue and pathways related to kidney development. These findings shed light on epigenetic mechanisms associated with kidney function, bridging the gap between population-specific eGFR-associated DNAm and tissue-specific regulatory context.

Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals.

Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to ~1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequency ≤ 0.01) variant BP associations (P < 5 × 10-8), of which 32 were in new BP-associated loci and 55 were independent BP-associated single-nucleotide variants within known BP-associated regions. Average effects of rare variants (44% coding) were ~8 times larger than common variant effects and indicate potential candidate causal genes at new and known loci (for example, GATA5 and PLCB3). BP-associated variants (including rare and common) were enriched in regions of active chromatin in fetal tissues, potentially linking fetal development with BP regulation in later life. Multivariable Mendelian randomization suggested possible inverse effects of elevated systolic and diastolic BP on large artery stroke. Our study demonstrates the utility of rare-variant analyses for identifying candidate genes and the results highlight potential therapeutic targets.

Methodology capture: discriminating between the "best" and the rest of community practice.

BACKGROUND: The methodologies we use both enable and help define our research. However, as experimental complexity has increased the choice of appropriate methodologies has become an increasingly difficult task. This makes it difficult to keep track of available bioinformatics software, let alone the most suitable protocols in a specific research area. To remedy this we present an approach for capturing methodology from literature in order to identify and, thus, define best practice within a field. RESULTS: Our approach is to implement data extraction techniques on the full-text of scientific articles to obtain the set of experimental protocols used by an entire scientific discipline, molecular phylogenetics. Our methodology for identifying methodologies could in principle be applied to any scientific discipline, whether or not computer-based. We find a number of issues related to the nature of best practice, as opposed to community practice. We find that there is much heterogeneity in the use of molecular phylogenetic methods and software, some of which is related to poor specification of protocols. We also find that phylogenetic practice exhibits field-specific tendencies that have increased through time, despite the generic nature of the available software. We used the practice of highly published and widely collaborative researchers ("expert" researchers) to analyse the influence of authority on community practice. We find expert authors exhibit patterns of practice common to their field and therefore act as useful field-specific practice indicators. CONCLUSION: We have identified a structured community of phylogenetic researchers performing analyses that are customary in their own local community and significantly different from those in other areas. Best practice information can help to bridge such subtle differences by increasing communication of protocols to a wider audience. We propose that the practice of expert authors from the field of evolutionary biology is the closest to contemporary best practice in phylogenetic experimental design. Capturing best practice is, however, a complex task and should also acknowledge the differences between fields such as the specific context of the analysis.

Molecular insights into genome-wide association studies of chronic kidney disease-defining traits.

Genome-wide association studies (GWAS) have identified >100 loci of chronic kidney disease-defining traits (CKD-dt). Molecular mechanisms underlying these associations remain elusive. Using 280 kidney transcriptomes and 9958 gene expression profiles from 44 non-renal tissues we uncover gene expression partners (eGenes) for 88.9% of CKD-dt GWAS loci. Through epigenomic chromatin segmentation analysis and variant effect prediction we annotate functional consequences to 74% of these loci. Our colocalisation analysis and Mendelian randomisation in >130,000 subjects demonstrate causal effects of three eGenes (NAT8B, CASP9 and MUC1) on estimated glomerular filtration rate. We identify a common alternative splice variant in MUC1 (a gene responsible for rare Mendelian form of kidney disease) and observe increased renal expression of a specific MUC1 mRNA isoform as a plausible molecular mechanism of the GWAS association signal. These data highlight the variants and genes underpinning the associations uncovered in GWAS of CKD-dt.

A Novel Y-Specific Long Non-Coding RNA Associated with Cellular Lipid Accumulation in HepG2 cells and Atherosclerosis-related Genes.

There is an increasing appreciation for the role of the human Y chromosome in phenotypic differences between the sexes in health and disease. Previous studies have shown that genetic variation within the Y chromosome is associated with cholesterol levels, which is an established risk factor for atherosclerosis, the underlying cause of coronary artery disease (CAD), a major cause of morbidity and mortality worldwide. However, the exact mechanism and potential genes implicated are still unidentified. To date, Y chromosome-linked long non-coding RNAs (lncRNAs) are poorly characterized and the potential link between these new regulatory RNA molecules and hepatic function in men has not been investigated. Advanced technologies of lncRNA subcellular localization and silencing were used to identify a novel intergenic Y-linked lncRNA, named lnc-KDM5D-4, and investigate its role in fatty liver-associated atherosclerosis. We found that lnc-KDM5D-4 is retained within the nucleus in hepatocytes. Its knockdown leads to changes in genes leading to increased lipid droplets formation in hepatocytes resulting in a downstream effect contributing to the chronic inflammatory process that underpin CAD. Our findings provide the first evidence for the implication of lnc-KDM5D-4 in key processes related to fatty liver and cellular inflammation associated with atherosclerosis and CAD in men.

RIPOSTE: a framework for improving the design and analysis of laboratory-based research.

Lack of reproducibility is an ongoing problem in some areas of the biomedical sciences. Poor experimental design and a failure to engage with experienced statisticians at key stages in the design and analysis of experiments are two factors that contribute to this problem. The RIPOSTE (Reducing IrreProducibility in labOratory STudiEs) framework has been developed to support early and regular discussions between scientists and statisticians in order to improve the design, conduct and analysis of laboratory studies and, therefore, to reduce irreproducibility. This framework is intended for use during the early stages of a research project, when specific questions or hypotheses are proposed. The essential points within the framework are explained and illustrated using three examples (a medical equipment test, a macrophage study and a gene expression study). Sound study design minimises the possibility of bias being introduced into experiments and leads to higher quality research with more reproducible results.