Mendelian randomization refers to an analytic approach to assess the causality of an observed association between a modifiable exposure or risk factor and a clinically relevant outcome. It presents a ...valuable tool, especially when randomized controlled trials to examine causality are not feasible and observational studies provide biased associations because of confounding or reverse causality. These issues are addressed by using genetic variants as instrumental variables for the tested exposure: the alleles of this exposure-associated genetic variant are randomly allocated and not subject to reverse causation. This, together with the wide availability of published genetic associations to screen for suitable genetic instrumental variables make Mendelian randomization a time- and cost-efficient approach and contribute to its increasing popularity for assessing and screening for potentially causal associations. An observed association between the genetic instrumental variable and the outcome supports the hypothesis that the exposure in question is causally related to the outcome. This review provides an overview of the Mendelian randomization method, addresses assumptions and implications, and includes illustrative examples. We also discuss special issues in nephrology, such as inverse risk factor associations in advanced disease, and outline opportunities to design Mendelian randomization studies around kidney function and disease.
Blood pressure and kidney function have a bidirectional relation. Hypertension has long been considered as a risk factor for kidney function decline. However, whether intensive blood pressure control ...could promote kidney health has been uncertain. The kidney is known to have a major role in affecting blood pressure through sodium extraction and regulating electrolyte balance. This bidirectional relation makes causal inference between these two traits difficult. Therefore, to examine the causal relations between these two traits, we performed two-sample Mendelian randomization analyses using summary statistics of large-scale genome-wide association studies. We selected genetic instruments more likely to be specific for kidney function using meta-analyses of complementary kidney function biomarkers (glomerular filtration rate estimated from serum creatinine eGFRcr, and blood urea nitrogen from the CKDGen Consortium). Systolic and diastolic blood pressure summary statistics were from the International Consortium for Blood Pressure and UK Biobank. Significant evidence supported the causal effects of higher kidney function on lower blood pressure. Based on the mode-based Mendelian randomization method, the effect estimates for one standard deviation (SD) higher in log-transformed eGFRcr was -0.17 SD unit (95 % confidence interval: -0.09 to -0.24) in systolic blood pressure and -0.15 SD unit (95% confidence interval: -0.07 to -0.22) in diastolic blood pressure. In contrast, the causal effects of blood pressure on kidney function were not statistically significant. Thus, our results support causal effects of higher kidney function on lower blood pressure and suggest preventing kidney function decline can reduce the public health burden of hypertension.
Display omitted
The new sequencing technologies enable to scan very long and dense genetic sequences, obtaining datasets of genetic markers that are an order of magnitude larger than previously available. Such ...genetic sequences are characterized by common alleles interspersed with multiple rarer alleles. This situation has renewed the interest for the identification of haplotypes carrying the rare risk alleles. However, large scale explorations of the linkage-disequilibrium (LD) pattern to identify haplotype blocks are not easy to perform, because traditional algorithms have at least Θ(n2) time and memory complexity.
We derived three incremental optimizations of the widely used haplotype block recognition algorithm proposed by Gabriel et al. in 2002. Our most efficient solution, called MIG ++, has only Θ(n) memory complexity and, on a genome-wide scale, it omits >80% of the calculations, which makes it an order of magnitude faster than the original algorithm. Differently from the existing software, the MIG ++ analyzes the LD between SNPs at any distance, avoiding restrictions on the maximal block length. The haplotype block partition of the entire HapMap II CEPH dataset was obtained in 457 hours. By replacing the standard likelihood-based D' variance estimator with an approximated estimator, the runtime was further improved. While producing a coarser partition, the approximate method allowed to obtain the full-genome haplotype block partition of the entire 1000 Genomes Project CEPH dataset in 44 hours, with no restrictions on allele frequency or long-range correlations. These experiments showed that LD-based haplotype blocks can span more than one million base-pairs in both HapMap II and 1000 Genomes datasets. An application to the North American Rheumatoid Arthritis Consortium (NARAC) dataset shows how the MIG ++ can support genome-wide haplotype association studies.
The MIG ++ enables to perform LD-based haplotype block recognition on genetic sequences of any length and density. In the new generation sequencing era, this can help identify haplotypes that carry rare variants of interest. The low computational requirements open the possibility to include the haplotype block structure into genome-wide association scans, downstream analyses, and visual interfaces for online genome browsers.
Levodopa is the standard long-term dopamine replacement therapy to treat Parkinson's disease (PD) symptoms. With time, levodopa may induce debilitating dyskinesias (LID), the treatment of which ...represents a large clinically unmet need. However, time-to-LID onset varies between patients, reflecting a possible genetic component. We performed an hypothesis-free whole-exome sequencing (WES)-based screening of time-to-LID onset and attempted replication of previously published candidate gene studies. A WES association analysis was carried out in 134 PD patients in a meta-analytical framework. Replication was attempted in an independent study of 97 PD patients. Variants from previously reported candidate genes (OPRM1, COMT, BDNF) were also specifically examined. We significantly replicated, for the first time, an association of variant rs1799971 in the OPRM1 gene with time-to-LID onset. Furthermore, we identified two novel potentially functional variants, in the MAD2L2 (rs2233019) and MAP7 (rs35350783) genes, which were significantly associated at the discovery stage. In the replication study, the two variants showed direction-consistent effects but did not achieve the replication significance threshold. Our study provides the first WES results for time-to-LID onset, where we replicate association at OPRM1, and suggest new variants in MAD2L2 and MAP7 genes that are significant in discovery, but require larger datasets for replication. The results are being made publicly available to allow for independent external validation.
Lower kidney function is known to enhance cardiovascular disease (CVD) risk. It is unclear which estimated glomerular filtration rate (eGFR) equation best predict an increased CVD risk and if ...prediction can be improved by integration of multiple kidney function markers. We performed structural equation modeling (SEM) of kidney markers and compared the performance of the resulting pooled indexes with established eGFR equations to predict CVD risk in a 10-year longitudinal population-based design. We split the study sample into a set of participants with only baseline data (n = 647; model-building set) and a set with longitudinal data (n = 670; longitudinal set). In the model-building set, we fitted five SEM models based on serum creatinine or creatinine-based eGFR (eGFRcre), cystatin C or cystatin-based eGFR (eGFRcys), uric acid (UA), and blood urea nitrogen (BUN). In the longitudinal set, 10-year incident CVD risk was defined as a Framingham risk score (FRS)>5% and a pooled cohort equation (PCE)>5%. Predictive performances of the different kidney function indexes were compared using the C-statistic and the DeLong test. In the longitudinal set, a SEM-based estimate of latent kidney function based on eGFRcre, eGFRcys, UA, and BUN showed better prediction performance for both FRS>5% (C-statistic: 0.70; 95% CI: 0.65-0.74) and PCE>5% (C-statistic: 0.75; 95%CI: 0.71-0.79) than other SEM models and different eGFR formulas (DeLong test p-values<3.21×10-6 for FRS>5% and <1.49×10-9 for PCE>5%, respectively). However, the new derived marker could not outperform eGFRcys (DeLong test p-values = 0.88 for FRS>5% and 0.20 for PCE>5%, respectively). SEM is a promising approach to identify latent kidney function signatures. However, for incident CVD risk prediction, eGFRcys could still be preferrable given its simpler derivation.
Heart rate variability (HRV) reflects the autonomous nervous system modulation on heart rate and is associated with several pathologies, including cardiac mortality. While mechanistic studies show ...that smoking is associated with lower HRV, population-based studies present conflicting results.
We assessed the mutual effects of active smoking status, cumulative smoking history, and current smoking intensity, on HRV among 4751 adults from the Cooperative Health Research In South Tyrol (CHRIS) study. The HRV metrics standard deviation of normal-to-normal (NN) inter-beat intervals (SDNN), square root of the mean squared differences of consecutive NN intervals (RMSSD), total power (TP), low (LF) and high frequency (HF) power, and their ratio (LF/HF), were derived from 20-minute electrocardiograms. Smoking status, pack-years (PY), and tobacco grams/day from standardized questionnaires were the main exposures. We fitted linear mixed models to account for relatedness, non-linearity, and moderating effects, and including fractional polynomials.
Past smokers had higher HRV levels than never smokers, independently of PY. The association of HRV with current smoking became apparent when accounting for the interaction between smoking status and PY. In current smokers, but not in past smokers, we observed HRV reductions between 2.0% (SDNN) and 4.9% (TP) every 5 PY increase. Furthermore, current smokers were characterized by dose-response reductions of 9.8% (SDNN), 8.9% (RMSSD), 20.1% (TP), 17.7% (LF), and 19.1% (HF), respectively, every 10 grams/day of smoked tobacco, independently of common cardiometabolic conditions and HRV-modifying drugs. The LF/HF ratio was not associated with smoking status, history, or intensity.
Smoking cessation was associated with higher HRV levels. In current smokers, heavier smoking intensity appears gradually detrimental on HRV, corroborating previous evidence. By affecting both the sympathetic and parasympathetic nervous system indexes, but not the LF/HF balance, smoking intensity seems to exert a systemic dysautonomic effect.
Genome-wide association studies have discovered hundreds of associations between common genotypes and kidney function but cannot comprehensively investigate rare coding variants. Here, we apply a ...genotype imputation approach to whole exome sequencing data from the UK Biobank to increase sample size from 166,891 to 408,511. We detect 158 rare variants and 105 genes significantly associated with one or more of five kidney function traits, including genes not previously linked to kidney disease in humans. The imputation-powered findings derive support from clinical record-based kidney disease information, such as for a previously unreported splice allele in PKD2, and from functional studies of a previously unreported frameshift allele in CLDN10. This cost-efficient approach boosts statistical power to detect and characterize both known and novel disease susceptibility variants and genes, can be generalized to larger future studies, and generates a comprehensive resource ( https://ckdgen-ukbb.gm.eurac.edu/ ) to direct experimental and clinical studies of kidney disease.
The oral microbiota plays an important role in the exogenous nitrate reduction pathway and is associated with heart and periodontal disease and cigarette smoking. We describe smoking-related changes ...in oral microbiota composition and resulting potential metabolic pathway changes that may explain smoking-related changes in disease risk. We analyzed health information and salivary microbiota composition among 1601 Cooperative Health Research in South Tyrol participants collected 2017-2018. Salivary microbiota taxa were assigned from amplicon sequences of the 16S-V4 rRNA and used to describe microbiota composition and predict metabolic pathways. Aerobic taxa relative abundance decreased with daily smoking intensity and increased with years since cessation, as did inferred nitrate reduction. Former smokers tended to be more similar to Never smokers than to Current smokers, especially those who had quit for longer than 5 years. Cigarette smoking has a consistent, generalizable association on oral microbiota composition and predicted metabolic pathways, some of which associate in a dose-dependent fashion. Smokers who quit for longer than 5 years tend to have salivary microbiota profiles comparable to never smokers.
The identification of genetic factors associated with kidney disease has the potential to provide critical insights into disease mechanisms. Genome-wide association studies have uncovered genomic ...regions associated with renal function metrics and risk of CKD.
is among the most outstanding loci associated with CKD in the general population, because it has a large effect on eGFR and CKD risk that is consistent across different ethnic groups. The relevance of
for CKD is clear, because the encoded protein, uromodulin (Tamm-Horsfall protein), is exclusively produced by the kidney tubule and has specific biochemical properties that mediate important functions in the kidney and urine. Rare mutations in
are the major cause of autosomal dominant tubulointerstitial kidney disease, a condition that leads to CKD and ESRD. In this brief review, we use the
paradigm to describe how population genetic studies can yield insight into the pathogenesis and prognosis of kidney diseases.
Mitochondrial DNA copy number (mtDNA-CN) is a biomarker for mitochondrial dysfunction associated with several diseases. Previous genome-wide association studies (GWAS) have been performed to unravel ...underlying mechanisms of mtDNA-CN regulation. However, the identified gene regions explain only a small fraction of mtDNA-CN variability. Most of this data has been estimated from microarrays based on various pipelines. In the present study we aimed to (1) identify genetic loci for qPCR-measured mtDNA-CN from three studies (16,130 participants) using GWAS, (2) identify potential systematic differences between our qPCR derived mtDNA-CN measurements compared to the published microarray intensity-based estimates, and (3) disentangle the nuclear from mitochondrial regulation of the mtDNA-CN phenotype. We identified two genome-wide significant autosomal loci associated with qPCR-measured mtDNA-CN: at HBS1L (rs4895440, p = 3.39 × 10
) and GSDMA (rs56030650, p = 4.85 × 10
) genes. Moreover, 113/115 of the previously published SNPs identified by microarray-based analyses were significantly equivalent with our findings. In our study, the mitochondrial genome itself contributed only marginally to mtDNA-CN regulation as we only detected a single rare mitochondrial variant associated with mtDNA-CN. Furthermore, we incorporated mitochondrial haplogroups into our analyses to explore their potential impact on mtDNA-CN. However, our findings indicate that they do not exert any significant influence on our results.