Fabry disease with atypical phenotype identified by massively parallel sequencing in early-onset kidney failure.

Background: The cause of chronic kidney disease (CKD) remains unknown in ∼20% of patients with kidney failure. Massively parallel sequencing (MPS) can be a valuable diagnostic tool in patients with unexplained CKD, with a diagnostic yield of 12%-56%. Here, we report the use of MPS to establish a genetic diagnosis in a 24-year-old index patient who presented with hypertension, nephrotic-range proteinuria and kidney failure of unknown origin. Additionally, we describe a second family with the same mutation presenting with early-onset CKD. Results: In Family 1, MPS identified a known pathogenic variant in GLA (p.Ile319Thr), and plasma globotriaosylsphingosine and α-galactosidase A activity were compatible with the diagnosis of Fabry disease (FD). Segregation analysis identified three other family members carrying the same pathogenic variant who had mild or absent kidney phenotypes. One family member was offered enzyme therapy. While FD could not be established with certainty as the cause of kidney failure in the index patient, no alternative explanation was found. In Family 2, the index patient had severe glomerulosclerosis and a kidney biopsy compatible with FD at the age of 30 years, along with cardiac involvement and a history of acroparesthesia since childhood, in keeping with a more classical Fabry phenotype. Conclusion: These findings highlight the large phenotypic heterogeneity associated with GLA mutations in FD and underline several important implications of MPS in the work-up of patients with unexplained kidney failure.

Diagnostic yield of massively parallel sequencing in patients with chronic kidney disease of unknown etiology: rationale and design of a national prospective cohort study.

INTRODUCTION: Chronic kidney disease (CKD) can be caused by a variety of systemic or primary renal diseases. The cause of CKD remains unexplained in approximately 20% of patients. Retrospective studies indicate that massively parallel sequencing (MPS)-based gene panel testing may lead to a genetic diagnosis in 12%-56% of patients with unexplained CKD, depending on patient profile. The diagnostic yield of MPS-based testing in a routine healthcare setting is unclear. Therefore, the primary aim of the VARIETY (Validation of algoRithms and IdEnTification of genes in Young patients with unexplained CKD) study is to prospectively address the diagnostic yield of MPS-based gene panel testing in patients with unexplained CKD and an estimated glomerular filtration rate (eGFR) <60 mL/min/1.73 m2 before the age of 50 years in clinical practice. METHODS AND ANALYSIS: The VARIETY study is an ongoing, prospective, nationwide observational cohort study to investigate the diagnostic yield of MPS-based testing in patients with unexplained CKD in a routine healthcare setting in the Netherlands. Patients are recruited from outpatient clinics in hospitals across the Netherlands. At least 282 patients will be included to meet the primary aim. Secondary analyses include subgroup analyses according to age and eGFR at first presentation, family history, and the presence of extrarenal symptoms. ETHICS AND DISSEMINATION: Ethical approval for the study has been obtained from the institutional review board of the University Medical Center Groningen. Study findings should inform physicians and policymakers towards optimal implementation of MPS-based diagnostic testing in patients with unexplained CKD.

Acute kidney injury in patients treated with immune checkpoint inhibitors.

BACKGROUND: Immune checkpoint inhibitor-associated acute kidney injury (ICPi-AKI) has emerged as an important toxicity among patients with cancer. METHODS: We collected data on 429 patients with ICPi-AKI and 429 control patients who received ICPis contemporaneously but who did not develop ICPi-AKI from 30 sites in 10 countries. Multivariable logistic regression was used to identify predictors of ICPi-AKI and its recovery. A multivariable Cox model was used to estimate the effect of ICPi rechallenge versus no rechallenge on survival following ICPi-AKI. RESULTS: ICPi-AKI occurred at a median of 16 weeks (IQR 8-32) following ICPi initiation. Lower baseline estimated glomerular filtration rate, proton pump inhibitor (PPI) use, and extrarenal immune-related adverse events (irAEs) were each associated with a higher risk of ICPi-AKI. Acute tubulointerstitial nephritis was the most common lesion on kidney biopsy (125/151 biopsied patients [82.7%]). Renal recovery occurred in 276 patients (64.3%) at a median of 7 weeks (IQR 3-10) following ICPi-AKI. Treatment with corticosteroids within 14 days following ICPi-AKI diagnosis was associated with higher odds of renal recovery (adjusted OR 2.64; 95% CI 1.58 to 4.41). Among patients treated with corticosteroids, early initiation of corticosteroids (within 3 days of ICPi-AKI) was associated with a higher odds of renal recovery compared with later initiation (more than 3 days following ICPi-AKI) (adjusted OR 2.09; 95% CI 1.16 to 3.79). Of 121 patients rechallenged, 20 (16.5%) developed recurrent ICPi-AKI. There was no difference in survival among patients rechallenged versus those not rechallenged following ICPi-AKI. CONCLUSIONS: Patients who developed ICPi-AKI were more likely to have impaired renal function at baseline, use a PPI, and have extrarenal irAEs. Two-thirds of patients had renal recovery following ICPi-AKI. Treatment with corticosteroids was associated with improved renal recovery.

Genetic Determinants of Serum Calcification Propensity and Cardiovascular Outcomes in the General Population.

BACKGROUND: Serum calciprotein particle maturation time (T50), a measure of vascular calcification propensity, is associated with cardiovascular morbidity and mortality. We aimed to identify genetic loci associated with serum T50 and study their association with cardiovascular disease and mortality. METHODS: We performed a genome-wide association study of serum T50 in 2,739 individuals of European descent participating in the Prevention of REnal and Vascular ENd-stage Disease (PREVEND) study, followed by a two-sample Mendelian randomization (MR) study to examine causal effects of T50 on cardiovascular outcomes. Finally, we examined associations between T50 loci and cardiovascular outcomes in 8,566 community-dwelling participants in the Rotterdam study. RESULTS: We identified three independent genome-wide significant single nucleotide polymorphism (SNPs) in the AHSG gene encoding fetuin-A: rs4917 (p = 1.72 × 10-101), rs2077119 (p = 3.34 × 10-18), and rs9870756 (p = 3.10 × 10-8), together explaining 18.3% of variation in serum T50. MR did not demonstrate a causal effect of T50 on cardiovascular outcomes in the general population. Patient-level analyses revealed that the minor allele of rs9870756, which explained 9.1% of variation in T50, was associated with a primary composite endpoint of all-cause mortality or cardiovascular disease [odds ratio (95% CI) 1.14 (1.01-1.28)] and all-cause mortality alone [1.14 (1.00-1.31)]. The other variants were not associated with clinical outcomes. In patients with type 2 diabetes or chronic kidney disease, the association between rs9870756 and the primary composite endpoint was stronger [OR 1.40 (1.06-1.84), relative excess risk due to interaction 0.54 (0.01-1.08)]. CONCLUSIONS: We identified three SNPs in the AHSG gene that explained 18.3% of variability in serum T50 levels. Only one SNP was associated with cardiovascular outcomes, particularly in individuals with type 2 diabetes or chronic kidney disease.

Corrigendum: Normal Values of Corrected Heart-Rate Variability in 10-Second Electrocardiograms for All Ages.

[This corrects the article DOI: 10.3389/fphys.2018.00424.].

Genome-wide association meta-analysis of 30,000 samples identifies seven novel loci for quantitative ECG traits.

Genome-wide association studies (GWAS) of quantitative electrocardiographic (ECG) traits in large consortia have identified more than 130 loci associated with QT interval, QRS duration, PR interval, and heart rate (RR interval). In the current study, we meta-analyzed genome-wide association results from 30,000 mostly Dutch samples on four ECG traits: PR interval, QRS duration, QT interval, and RR interval. SNP genotype data was imputed using the Genome of the Netherlands reference panel encompassing 19 million SNPs, including millions of rare SNPs (minor allele frequency < 5%). In addition to many known loci, we identified seven novel locus-trait associations: KCND3, NR3C1, and PLN for PR interval, KCNE1, SGIP1, and NFKB1 for QT interval, and ATP2A2 for QRS duration, of which six were successfully replicated. At these seven loci, we performed conditional analyses and annotated significant SNPs (in exons and regulatory regions), demonstrating involvement of cardiac-related pathways and regulation of nearby genes.

Diagnostic Yield of Next-Generation Sequencing in Patients With Chronic Kidney Disease of Unknown Etiology.

Advances in next-generation sequencing (NGS) techniques, including whole exome sequencing, have facilitated cost-effective sequencing of large regions of the genome, enabling the implementation of NGS in clinical practice. Chronic kidney disease (CKD) is a major contributor to global burden of disease and is associated with an increased risk of morbidity and mortality. CKD can be caused by a wide variety of primary renal disorders. In about one in five CKD patients, no primary renal disease diagnosis can be established. Moreover, recent studies indicate that the clinical diagnosis may be incorrect in a substantial number of patients. Both the absence of a diagnosis or an incorrect diagnosis can have therapeutic implications. Genetic testing might increase the diagnostic accuracy in patients with CKD, especially in patients with unknown etiology. The diagnostic utility of NGS has been shown mainly in pediatric CKD cohorts, while emerging data suggest that genetic testing can also be a valuable diagnostic tool in adults with CKD. In addition to its implications for unexplained CKD, NGS can contribute to the diagnostic process in kidney diseases with an atypical presentation, where it may lead to reclassification of the primary renal disease diagnosis. So far, only a few studies have reported on the diagnostic yield of NGS-based techniques in patients with unexplained CKD. Here, we will discuss the potential diagnostic role of gene panels and whole exome sequencing in pediatric and adult patients with unexplained and atypical CKD.

Exome-chip meta-analysis identifies novel loci associated with cardiac conduction, including ADAMTS6.

BACKGROUND: Genome-wide association studies conducted on QRS duration, an electrocardiographic measurement associated with heart failure and sudden cardiac death, have led to novel biological insights into cardiac function. However, the variants identified fall predominantly in non-coding regions and their underlying mechanisms remain unclear. RESULTS: Here, we identify putative functional coding variation associated with changes in the QRS interval duration by combining Illumina HumanExome BeadChip genotype data from 77,898 participants of European ancestry and 7695 of African descent in our discovery cohort, followed by replication in 111,874 individuals of European ancestry from the UK Biobank and deCODE cohorts. We identify ten novel loci, seven within coding regions, including ADAMTS6, significantly associated with QRS duration in gene-based analyses. ADAMTS6 encodes a secreted metalloprotease of currently unknown function. In vitro validation analysis shows that the QRS-associated variants lead to impaired ADAMTS6 secretion and loss-of function analysis in mice demonstrates a previously unappreciated role for ADAMTS6 in connexin 43 gap junction expression, which is essential for myocardial conduction. CONCLUSIONS: Our approach identifies novel coding and non-coding variants underlying ventricular depolarization and provides a possible mechanism for the ADAMTS6-associated conduction changes.

PR interval genome-wide association meta-analysis identifies 50 loci associated with atrial and atrioventricular electrical activity.

Electrocardiographic PR interval measures atrio-ventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. Our genome-wide association study of over 92,000 European-descent individuals identifies 44 PR interval loci (34 novel). Examination of these loci reveals known and previously not-yet-reported biological processes involved in cardiac atrial electrical activity. Genes in these loci are over-represented in cardiac disease processes including heart block and atrial fibrillation. Variants in over half of the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with missense variants. Six additional loci were identified either by meta-analysis of ~105,000 African and European-descent individuals and/or by pleiotropic analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation. These findings implicate developmental pathways, and identify transcription factors, ion-channel genes, and cell-junction/cell-signaling proteins in atrio-ventricular conduction, identifying potential targets for drug development.

Normal Values of Corrected Heart-Rate Variability in 10-Second Electrocardiograms for All Ages.

Purpose: Heart-rate variability (HRV) measured on standard 10-s electrocardiograms (ECGs) has been associated with increased risk of cardiac and all-cause mortality, but age- and sex-dependent normal values have not been established. Since heart rate strongly affects HRV, its effect should be taken into account. We determined a comprehensive set of normal values of heart-rate corrected HRV derived from 10-s ECGs for both children and adults, covering both sexes. Methods: Five population studies in the Netherlands (Pediatric Normal ECG Study, Leiden University Einthoven Science Project, Prevention of Renal and Vascular End-stage Disease Study, Utrecht Health Project, Rotterdam Study) provided 10-s, 12-lead ECGs. ECGs were stored digitally and analyzed by well-validated analysis software. We included cardiologically healthy participants, 42% being men. Their ages ranged from 11 days to 91 years. After quality control, 13,943 ECGs were available. Heart-rate correction formulas were derived using an exponential model. Two time-domain HRV markers were analyzed: the corrected standard deviation of the normal-to-normal RR intervals (SDNNc) and corrected root mean square of successive RR-interval differences (RMSSDc). Results: There was a considerable age effect. For both SDNNc and RMSSDc, the median and the lower limit of normal decreased steadily from birth until old age. The upper limit of normal decreased until the age of 60, but increased markedly after that age. Differences of the median were minimal between men and women. Conclusion: We report the first comprehensive set of normal values for heart-rate corrected 10-s HRV, which can be of value in clinical practice and in further research.

52 Genetic Loci Influencing Myocardial Mass.

BACKGROUND: Myocardial mass is a key determinant of cardiac muscle function and hypertrophy. Myocardial depolarization leading to cardiac muscle contraction is reflected by the amplitude and duration of the QRS complex on the electrocardiogram (ECG). Abnormal QRS amplitude or duration reflect changes in myocardial mass and conduction, and are associated with increased risk of heart failure and death. OBJECTIVES: This meta-analysis sought to gain insights into the genetic determinants of myocardial mass. METHODS: We carried out a genome-wide association meta-analysis of 4 QRS traits in up to 73,518 individuals of European ancestry, followed by extensive biological and functional assessment. RESULTS: We identified 52 genomic loci, of which 32 are novel, that are reliably associated with 1 or more QRS phenotypes at p < 1 × 10(-8). These loci are enriched in regions of open chromatin, histone modifications, and transcription factor binding, suggesting that they represent regions of the genome that are actively transcribed in the human heart. Pathway analyses provided evidence that these loci play a role in cardiac hypertrophy. We further highlighted 67 candidate genes at the identified loci that are preferentially expressed in cardiac tissue and associated with cardiac abnormalities in Drosophila melanogaster and Mus musculus. We validated the regulatory function of a novel variant in the SCN5A/SCN10A locus in vitro and in vivo. CONCLUSIONS: Taken together, our findings provide new insights into genes and biological pathways controlling myocardial mass and may help identify novel therapeutic targets.

Thyroid Function and Sudden Cardiac Death: A Prospective Population-Based Cohort Study.

BACKGROUND: The association between thyroid function and cardiovascular disease is well established, but no study to date has assessed whether it is a risk factor for sudden cardiac death (SCD). Therefore, we studied the association of thyroid function with SCD in a prospective population-based cohort. METHODS: Participants from the Rotterdam Study ≥45 years with thyroid-stimulating hormone or free thyroxine (FT4) measurements and clinical follow-up were eligible. We assessed the association of thyroid-stimulating hormone and FT4 with the risk of SCD by using an age- and sex-adjusted Cox proportional-hazards model, in all participants and also after restricting the analysis to euthyroid participants (defined by thyroid-stimulating hormone 0.4-4.0 mIU/L). Additional adjustment included cardiovascular risk factors, notably hypertension, serum cholesterol, and smoking. We stratified by age and sex and performed sensitivity analyses by excluding participants with abnormal FT4 values (reference range of 0.85-1.95 ng/dL) and including only witnessed SCDs as outcome. Absolute risks were calculated in a competing risk model by taking death by other causes into account. RESULTS: We included 10 318 participants with 261 incident SCDs (median follow-up, 9.1 years). Higher levels of FT4 were associated with an increased SCD risk, even in the normal range of thyroid function (hazard ratio, 2.28 per 1 ng/dL FT4; 95% confidence interval, 1.31-3.97). Stratification by age or sex and sensitivity analyses did not change the risk estimates substantially. The absolute 10-year risk of SCD increased in euthyroid participants from 1% to 4% with increasing FT4 levels. CONCLUSIONS: Higher FT4 levels are associated with an increased risk of SCD, even in euthyroid participants.

Twenty-eight genetic loci associated with ST-T-wave amplitudes of the electrocardiogram.

The ST-segment and adjacent T-wave (ST-T wave) amplitudes of the electrocardiogram are quantitative characteristics of cardiac repolarization. Repolarization abnormalities have been linked to ventricular arrhythmias and sudden cardiac death. We performed the first genome-wide association meta-analysis of ST-T-wave amplitudes in up to 37 977 individuals identifying 71 robust genotype-phenotype associations clustered within 28 independent loci. Fifty-four genes were prioritized as candidates underlying the phenotypes, including genes with established roles in the cardiac repolarization phase (SCN5A/SCN10A, KCND3, KCNB1, NOS1AP and HEY2) and others with as yet undefined cardiac function. These associations may provide insights in the spatiotemporal contribution of genetic variation influencing cardiac repolarization and provide novel leads for future functional follow-up.

Association of Lipid-Related Genetic Variants with the Incidence of Atrial Fibrillation: The AFGen Consortium.

BACKGROUND: Several studies have shown associations between blood lipid levels and the risk of atrial fibrillation (AF). To test the potential effect of blood lipids with AF risk, we assessed whether previously developed lipid gene scores, used as instrumental variables, are associated with the incidence of AF in 7 large cohorts. METHODS: We analyzed 64,901 individuals of European ancestry without previous AF at baseline and with lipid gene scores. Lipid-specific gene scores, based on loci significantly associated with lipid levels, were calculated. Additionally, non-pleiotropic gene scores for high-density lipoprotein cholesterol (HDLc) and low-density lipoprotein cholesterol (LDLc) were calculated using SNPs that were only associated with the specific lipid fraction. Cox models were used to estimate the hazard ratio (HR) and 95% confidence intervals (CI) of AF per 1-standard deviation (SD) increase of each lipid gene score. RESULTS: During a mean follow-up of 12.0 years, 5434 (8.4%) incident AF cases were identified. After meta-analysis, the HDLc, LDLc, total cholesterol, and triglyceride gene scores were not associated with incidence of AF. Multivariable-adjusted HR (95% CI) were 1.01 (0.98-1.03); 0.98 (0.96-1.01); 0.98 (0.95-1.02); 0.99 (0.97-1.02), respectively. Similarly, non-pleiotropic HDLc and LDLc gene scores showed no association with incident AF: HR (95% CI) = 1.00 (0.97-1.03); 1.01 (0.99-1.04). CONCLUSIONS: In this large cohort study of individuals of European ancestry, gene scores for lipid fractions were not associated with incident AF.

Subclinical Abnormalities in Echocardiographic Parameters and Risk of Sudden Cardiac Death in a General Population: The Rotterdam Study.

BACKGROUND: Subclinical cardiac dysfunction has been associated with increased mortality, and heart failure increases the risk of sudden cardiac death (SCD). Less well known is whether subclinical cardiac dysfunction is also a risk factor for SCD. Our objective was to assess the association between echocardiographic parameters and SCD in a community-dwelling population free of heart failure. METHODS AND RESULTS: We computed hazard ratios (HRs) for left atrium diameter, left ventricular (LV) end-diastolic dimension, LV end-systolic dimension, LV mass, qualitative LV systolic function, LV fractional shortening, and diastolic function. During a median follow-up of 6.3 years in 4,686 participants, 68 participants died because of SCD. Significant associations with SCD were observed for qualitative LV systolic function and LV fractional shortening. For moderate/poor qualitative LV systolic function, the HR for SCD was 2.54 (95% confidence interval [CI] 1.10-5.87). Each standard deviation decrease in LV fractional shortening was associated with an HR of 1.36 (95% CI 1.09-1.70). CONCLUSIONS: Subclinical abnormalities in LV systolic function were associated with SCD risk in this general population. Although prediction of SCD remains difficult and traditional cardiovascular risk factors are of greatest importance, this knowledge might guide future directions to prevent SCD in persons with subclinical cardiac dysfunction.

ABCB1 gene variants, digoxin and risk of sudden cardiac death in a general population.

OBJECTIVE: The ATP-binding cassette B1 (ABCB1) gene encodes P-glycoprotein, a transport protein, which plays an important role in the bioavailability of digoxin. We aimed to investigate the interaction between variants within the ABCB1 gene and digoxin on the risk of sudden cardiac death (SCD). METHODS: Within the Rotterdam Study, a population-based cohort study in persons 45 years of age and older, we used Cox regression to analyse the association between three polymorphisms that have been associated with digoxin bioavailability, extracted from 1000-Genomes imputed ABCB1 genotypes and the risk of SCD, stratified by digoxin use. RESULTS: In a total study population of 10,932 persons, 419 SCDs occurred during a median follow-up of 9.8 years. In non-users of digoxin, the risk of SCD was not different across genotypes. In digoxin users, homozygous T allele carriers of C1236T (HR 1.90; 95% CI 1.09 to 3.30; allele frequency 0.43), G2677T (HR 1.89; 95% CI 1.10 to 3.24; allele frequency 0.44) and C3435T (HR 1.72; 95% CI 1.03 to 2.87; allele frequency 0.53) had a significantly increased risk of SCD in a recessive model. Interaction between the ABCB1 polymorphisms and digoxin use was significant for C1236T and G2677T in the age-adjusted and sex-adjusted model. CONCLUSIONS: In this study, we showed that in digoxin users variant alleles at each of the three loci in the ABCB1 gene were associated with an increased risk of SCD compared with digoxin users with none or one T allele. If replicated, the findings imply that the ABCB1 genotype modifies the risk of cardiac digoxin toxicity.

Drugs and ventricular repolarization in a general population: the Rotterdam Study.

PURPOSE: Prolonged ventricular repolarization (measured as heart-rate corrected QT (QTc) prolongation or JT-interval prolongation) is a risk factor for ventricular arrhythmias and can be drug-induced. Drugs can be classified as having known or possible QTc-prolonging properties. Regulatory agencies recommend avoiding concomitant use of multiple QTc-prolonging drugs, but evidence is lacking to what degree ventricular repolarization is influenced by concomitant use of these drugs. METHODS: Within a population-based cohort of persons aged 45 years and older, with up to five electrocardiograms recorded per participant between 1991 and 2010, we used generalised estimating equations to study the association between concomitant use of multiple QTc-prolonging drugs and repolarization duration. RESULTS: The study population consisted of 13 009 participants with 26 908 electrocardiograms. With the addition of a second or third QTc-prolonging drug there was no substantial increase in QTc and JT interval and no increased risk of a prolonged QTc interval, compared to use of one QTc-prolonging drug. There was a large difference between the effect of one known or one possible QTc-prolonging drugs on QTc interval: 15 ms for known, and 3 ms for possible QTc-prolonging drugs. CONCLUSIONS: In this study, the added prolongation in users of two or three QTc-prolonging drugs on QTc was small. There was a large difference in QTc prolongation between known and possible QTc-prolonging drugs. Further research in larger or high-risk populations is needed to establish whether it is safe to use multiple QTc-prolonging drugs concomitantly to prevent that the current advice might unnecessarily withhold beneficial drugs from patients.

Proton pump inhibitors and hypomagnesemia in the general population: a population-based cohort study.

BACKGROUND: Proton pump inhibitor (PPI) use has been associated with hypomagnesemia in case reports and hospital-based cohort studies. Our objective was to determine whether PPI use is associated with hypomagnesemia in the general population and whether this is also found in histamine 2 receptor antagonist (H2RA) users. STUDY DESIGN: Prospective cohort study. SETTING & PARTICIPANTS: 9,818 individuals from the general population (Rotterdam Study). PREDICTOR: PPI use and H2RA use compared to no use. OUTCOMES & MEASUREMENTS: Serum magnesium and hypomagnesemia (serum magnesium ≤ 1.44 mEq/L). Analyses were adjusted for age, sex, body mass index, kidney function, comorbid conditions, and alcohol and diuretic use. RESULTS: Serum magnesium level was 0.022 mEq/L lower in PPI users (n=724; 95% CI, -0.032 to -0.014 mEq/L) versus those with no use. PPI use was associated with increased risk of hypomagnesemia (n=36; OR, 2.00; 95% CI, 1.36-2.93) compared to no use. Effect modification was found between the use of PPIs and loop diuretics; in participants using loop diuretics (n=270), PPI use was associated with a further increased risk of hypomagnesemia (n=5; OR, 7.22; 95% CI, 1.69-30.83) compared to no use. The increased risk with PPIs was only seen after prolonged use (range, 182-2,618 days; OR, 2.99; 95% CI, 1.73-5.15). Including dietary magnesium intake into the model did not alter results (available for 2,504 participants, including 231 PPI users). H2RA users (n=250) also had a lower serum magnesium level (-0.016 [95% CI, -0.032 to -0.002] mEq/L) and increased risk of hypomagnesemia (n=12; OR, 2.00; 95% CI, 1.08-3.72) compared to those with no use, but no interaction with loop diuretics. LIMITATIONS: Cross-sectional analysis with single serum magnesium measurement. CONCLUSIONS: PPI use is associated with hypomagnesemia in the general population. Prolonged PPI use and concomitant loop diuretic use are associated with a stronger risk increase. Similar but weaker associations were found in H2RA users, except for interaction with loop diuretics.