Acute Kidney Injury, Systemic Inflammation and Long-term Cognitive Function: ASSESS-AKI

BACKGROUND

Cognitive dysfunction is a well-known complication of chronic kidney disease, but it is less known whether cognitive decline occurs in survivors after acute kidney injury (AKI). We hypothesized that an episode of AKI is associated with poorer cognitive function, mediated, at least in part, by persistent systemic inflammation.

METHODS

ASSESS-AKI enrolled patients surviving three months after hospitalization with and without AKI matched based on demographics, comorbidities, and baseline kidney function. A subset underwent cognitive testing using the modified mini-mental status examination (3MS) at 3, 12, and 36 months. We examined the association of AKI with 3MS scores using mixed linear models and assessed the proportion of risk mediated by systemic inflammatory biomarkers.

RESULTS

Among 1538 participants in ASSESS-AKI, 1420 (92%) completed the 3MS assessment at 3 months and had a corresponding matched participant. Participants with AKI had lower 3MS scores at three years (difference -1.1 (95% CI: -2.0, -0.3) P=0.009) compared to participants without AKI. A higher proportion of AKI participants had a clinically meaningful (≥ 5 point) reduction in 3MS scores at three years compared to participants without AKI (14% vs. 10%, P=0.04). In mediation analyses, plasma soluble tumor necrosis factor receptor-1 (sTNFR-1) at three months after AKI mediated 35% (P=0.02) of the AKI related risk for 3MS scores at three years.

CONCLUSIONS

AKI was associated with lower 3MS scores and sTNFR-1 concentrations appeared to mediate a significant proportion of the risk of long-term cognitive impairment. Further work is needed to determine if AKI is causal or a marker for cognitive impairment.

Identification of prognostic biomarkers for antibiotic associated nephrotoxicity in cystic fibrosis

BACKGROUND

Our objective was to discover novel urinary biomarkers of antibiotic-associated nephrotoxicity using an ex-vivo human microphysiological system (MPS) and to translate these findings to a prospectively enrolled cystic fibrosis (CF) population receiving aminoglycosides and/or polymyxin E (colistin) for a pulmonary exacerbation.

METHODS

We populated the MPS with primary human kidney proximal tubule epithelial cells (PTECs) from three donors and modeled nephrotoxin injury through exposure to 50 µg/mL polymyxin E for 72 h. We analyzed gene transcriptional responses by RNAseq and tested MPS effluents. We translated candidate biomarkers to a CF cohort via analysis of urine collected prior to, during and two weeks after antibiotics and patients were followed for a median of 3 years after antibiotic use.

RESULTS

Polymyxin E treatment resulted in a statistically significant increase in the pro-apoptotic Fas gene relative to control in RNAseq of MPS: fold-change = 1.63, FDR q-value = 7.29 × 10-5. Effluent analysis demonstrated an acute rise of soluble Fas (sFas) concentrations that correlated with cellular injury. In 16 patients with CF, urinary sFas concentrations were significantly elevated during antibiotic treatment, regardless of development of AKI. Over a median of three years of follow up, we identified seven cases of incident chronic kidney disease (CKD). Urinary sFas concentrations during antibiotic treatment were significantly associated with subsequent development of incident CKD (unadjusted relative risk = 2.02 per doubling of urinary sFas, 95 % CI = 1.40, 2.90, p < 0.001).

CONCLUSIONS

Using an ex-vivo MPS, we identified a novel biomarker of proximal tubule epithelial cell injury, sFas, and translated these findings to a clinical cohort of patients with CF.

Mitochondrial N-formyl methionine peptides contribute to exaggerated neutrophil activation in patients with COVID-19

Neutrophil dysregulation is well established in COVID-19. However, factors contributing to neutrophil activation in COVID-19 are not clear. We assessed if N-formyl methionine (fMet) contributes to neutrophil activation in COVID-19. Elevated levels of calprotectin, neutrophil extracellular traps (NETs) and fMet were observed in COVID-19 patients (n = 68), particularly in critically ill patients, as compared to HC (n = 19, p < 0.0001). Of note, the levels of NETs were higher in ICU patients with COVID-19 than in ICU patients without COVID-19 (p < 0.05), suggesting a prominent contribution of NETs in COVID-19. Additionally, plasma from COVID-19 patients with mild and moderate/severe symptoms induced in vitro neutrophil activation through fMet/FPR1 (formyl peptide receptor-1) dependent mechanisms (p < 0.0001). fMet levels correlated with calprotectin levels validating fMet-mediated neutrophil activation in COVID-19 patients (r = 0.60, p = 0.0007). Our data indicate that fMet is an important factor contributing to neutrophil activation in COVID-19 disease and may represent a potential target for therapeutic intervention.

Genome-wide Association Study for AKI

Key Points

Two genetic variants in the DISP1-TLR5 gene locus were associated with risk of AKI. DISP1 and TLR5 were differentially regulated in kidney biopsy tissue from patients with AKI compared with no AKI.

Background

Although common genetic risks for CKD are well established, genetic factors influencing risk for AKI in hospitalized patients are poorly understood.

Methods

We conducted a genome-wide association study in 1369 participants in the Assessment, Serial Evaluation, and Subsequent Sequelae of AKI Study; a multiethnic population of hospitalized participants with and without AKI matched on demographics, comorbidities, and kidney function before hospitalization. We then completed functional annotation of top-performing variants for AKI using single-cell RNA sequencing data from kidney biopsies in 12 patients with AKI and 18 healthy living donors from the Kidney Precision Medicine Project.

Results

No genome-wide significant associations with AKI risk were found in Assessment, Serial Evaluation, and Subsequent Sequelae of AKI (P < 5×10 −8 ). The top two variants with the strongest association with AKI mapped to the dispatched resistance-nodulation-division (RND) transporter family member 1 (DISP1) gene and toll-like receptor 5 (TLR5) gene locus, rs17538288 (odds ratio, 1.55; 95% confidence interval, 1.32 to 182; P = 9.47×10 −8 ) and rs7546189 (odds ratio, 1.53; 95% confidence interval, 1.30 to 1.81; P = 4.60×10 −7 ). In comparison with kidney tissue from healthy living donors, kidney biopsies in patients with AKI showed differential DISP1 expression in proximal tubular epithelial cells (adjusted P = 3.9× 10−2) and thick ascending limb of the loop of Henle (adjusted P = 8.7× 10−3) and differential TLR5 gene expression in thick ascending limb of the loop of Henle (adjusted P = 4.9× 10−30).

Conclusions

AKI is a heterogeneous clinical syndrome with various underlying risk factors, etiologies, and pathophysiology that may limit the identification of genetic variants. Although no variants reached genome-wide significance, we report two variants in the intergenic region between DISP1 and TLR5 , suggesting this region as a novel risk for AKI susceptibility.

Association of Trauma Molecular Endotypes With Differential Response to Transfusion Resuscitation Strategies

Importance

It is not clear which severely injured patients with hemorrhagic shock may benefit most from a 1:1:1 vs 1:1:2 (plasma:platelets:red blood cells) resuscitation strategy. Identification of trauma molecular endotypes may reveal subgroups of patients with differential treatment response to various resuscitation strategies.

Objective

To derive trauma endotypes (TEs) from molecular data and determine whether these endotypes are associated with mortality and differential treatment response to 1:1:1 vs 1:1:2 resuscitation strategies.

Design, Setting, and Participants

This was a secondary analysis of the Pragmatic, Randomized Optimal Platelet and Plasma Ratios (PROPPR) randomized clinical trial. The study cohort included individuals with severe injury from 12 North American trauma centers. The cohort was taken from the participants in the PROPPR trial who had complete plasma biomarker data available. Study data were analyzed on August 2, 2021, to October 25, 2022.

Exposures

TEs identified by K-means clustering of plasma biomarkers collected at hospital arrival.

Main Outcomes and Measures

An association between TEs and 30-day mortality was tested using multivariable relative risk (RR) regression adjusting for age, sex, trauma center, mechanism of injury, and injury severity score (ISS). Differential treatment response to transfusion strategy was assessed using an RR regression model for 30-day mortality by incorporating an interaction term for the product of endotype and treatment group adjusting for age, sex, trauma center, mechanism of injury, and ISS.

Results

A total of 478 participants (median [IQR] age, 34.5 [25-51] years; 384 male [80%]) of the 680 participants in the PROPPR trial were included in this study analysis. A 2-class model that had optimal performance in K-means clustering was found. TE-1 (n = 270) was characterized by higher plasma concentrations of inflammatory biomarkers (eg, interleukin 8 and tumor necrosis factor α) and significantly higher 30-day mortality compared with TE-2 (n = 208). There was a significant interaction between treatment arm and TE for 30-day mortality. Mortality in TE-1 was 28.6% with 1:1:2 treatment vs 32.6% with 1:1:1 treatment, whereas mortality in TE-2 was 24.5% with 1:1:2 treatment vs 7.3% with 1:1:1 treatment (P for interaction = .001).

Conclusions and Relevance

Results of this secondary analysis suggest that endotypes derived from plasma biomarkers in trauma patients at hospital arrival were associated with a differential response to 1:1:1 vs 1:1:2 resuscitation strategies in trauma patients with severe injury. These findings support the concept of molecular heterogeneity in critically ill trauma populations and have implications for tailoring therapy for patients at high risk for adverse outcomes.

Acarbose suppresses symptoms of mitochondrial disease in a mouse model of Leigh syndrome

Mitochondrial diseases represent a spectrum of disorders caused by impaired mitochondrial function, ranging in severity from mortality during infancy to progressive adult-onset disease. Mitochondrial dysfunction is also recognized as a molecular hallmark of the biological ageing process. Rapamycin, a drug that increases lifespan and health during normative ageing, also increases survival and reduces neurological symptoms in a mouse model of the severe mitochondrial disease Leigh syndrome. The Ndufs4 knockout (Ndufs4^-/-) mouse lacks the complex I subunit NDUFS4 and shows rapid onset and progression of neurodegeneration mimicking patients with Leigh syndrome. Here we show that another drug that extends lifespan and delays normative ageing in mice, acarbose, also suppresses symptoms of disease and improves survival of Ndufs4^-/- mice. Unlike rapamycin, acarbose rescues disease phenotypes independently of inhibition of the mechanistic target of rapamycin. Furthermore, rapamycin and acarbose have additive effects in delaying neurological symptoms and increasing maximum lifespan in Ndufs4^-/- mice. We find that acarbose remodels the intestinal microbiome and alters the production of short-chain fatty acids. Supplementation with tributyrin, a source of butyric acid, recapitulates some effects of acarbose on lifespan and disease progression, while depletion of the endogenous microbiome in Ndufs4^-/- mice appears to fully recapitulate the effects of acarbose on healthspan and lifespan in these animals. To our knowledge, this study provides the first evidence that alteration of the gut microbiome plays a significant role in severe mitochondrial disease and provides further support for the model that biological ageing and severe mitochondrial disorders share underlying common mechanisms.

Alteration of oral microbiome composition in children living with pesticide-exposed farm workers

Our prior work shows that azinphos-methyl pesticide exposure is associated with altered oral microbiomes in exposed farmworkers. Here we extend this analysis to show the same association pattern is also evident in their children. Oral buccal swab samples were analyzed at two time points, the apple thinning season in spring-summer 2005 for 78 children and 101 adults and the non-spray season in winter 2006 for 62 children and 82 adults. The pesticide exposure for the children were defined by the farmworker occupation of the cohabitating household adult and the blood azinphos-methyl detection of the cohabitating adult. Oral buccal swab 16S rRNA sequencing determined taxonomic microbiota proportional composition from concurrent samples from both adults and children. Analysis of the identified bacteria showed significant proportional changes for 12 of 23 common oral microbiome genera in association with azinphos-methyl detection and farmworker occupation. The most common significantly altered genera had reductions in the abundance of Streptococcus, suggesting an anti-microbial effect of the pesticide. Principal component analysis of the microbiome identified two primary clusters, with association of principal component 1 to azinphos-methyl blood detection and farmworker occupational status of the household. The children's buccal microbiota composition clustered with their household adult in ∼95% of the households. Household adult farmworker occupation and household pesticide exposure is associated with significant alterations in their children's oral microbiome composition. This suggests that parental occupational exposure and pesticide take-home exposure pathways elicit alteration of their children's microbiomes.

Genetic Predictors of Blood Pressure Traits are Associated with Preeclampsia

Background

Preeclampsia, a pregnancy complication characterized by hypertension after 20 gestational weeks, is a major cause of maternal and neonatal morbidity and mortality. The mechanisms leading to preeclampsia are unclear; however, there is evidence that preeclampsia is highly heritable. We evaluated the association of polygenic risk scores (PRS) for blood pressure traits and preeclampsia to assess whether there is shared genetic architecture.

Methods

Participants were obtained from Vanderbilt University's BioVU, the Electronic Medical Records and Genomics network, and the Penn Medicine Biobank. Non-Hispanic Black and White females of reproductive age with indications of pregnancy and genotype information were included. Preeclampsia was defined by ICD codes. Summary statistics for diastolic blood pressure (DBP), systolic blood pressure (SBP), and pulse pressure (PP) PRS were obtained from Giri et al 2019. Associations between preeclampsia and each PRS were evaluated separately by race and study population before evidence was meta-analyzed. Prediction models were developed and evaluated using 10-fold cross validation.

Results

In the 3,504 Black and 5,009 White individuals included, the rate of preeclampsia was 15.49%. The DBP and SBP PRSs were associated with preeclampsia in Whites but not Blacks. The PP PRS was significantly associated with preeclampsia in Blacks and Whites. In trans-ancestry meta-analysis, all PRSs were associated with preeclampsia (OR DBP =1.10, 95% CI=1.02-1.17, p =7.68×10 -3 ; OR SBP =1.16, 95% CI=1.09-1.23, p =2.23×10 -6 ; OR PP =1.14, 95% CI=1.07-1.27, p =9.86×10 -5 ). However, addition of PRSs to clinical prediction models did not improve predictive performance.

Conclusions

Genetic factors contributing to blood pressure regulation in the general population also predispose to preeclampsia.

Angiopoietin-Like4 Is a Novel Marker of COVID-19 Severity

Vascular dysfunction and capillary leak are common in critically ill COVID-19 patients, but identification of endothelial pathways involved in COVID-19 pathogenesis has been limited. Angiopoietin-like 4 (ANGPTL4) is a protein secreted in response to hypoxic and nutrient-poor conditions that has a variety of biological effects including vascular injury and capillary leak.

OBJECTIVES

To assess the role of ANGPTL4 in COVID-19-related outcomes.

DESIGN SETTING AND PARTICIPANTS

Two hundred twenty-five COVID-19 ICU patients were enrolled from April 2020 to May 2021 in a prospective, multicenter cohort study from three different medical centers, University of Washington, University of Southern California and New York University.

MAIN OUTCOMES AND MEASURES

Plasma ANGPTL4 was measured on days 1, 7, and 14 after ICU admission. We used previously published tissue proteomic data and lung single nucleus RNA (snRNA) sequencing data from specimens collected from COVID-19 patients to determine the tissues and cells that produce ANGPTL4.

RESULTS

Higher plasma ANGPTL4 concentrations were significantly associated with worse hospital mortality (adjusted odds ratio per log2 increase, 1.53; 95% CI, 1.17-2.00; p = 0.002). Higher ANGPTL4 concentrations were also associated with higher proportions of venous thromboembolism and acute respiratory distress syndrome. Longitudinal ANGPTL4 concentrations were significantly different during the first 2 weeks of hospitalization in patients who subsequently died compared with survivors (p for interaction = 8.1 × 10-5). Proteomics analysis demonstrated abundance of ANGPTL4 in lung tissue compared with other organs in COVID-19. ANGPTL4 single-nuclear RNA gene expression was significantly increased in pulmonary alveolar type 2 epithelial cells and fibroblasts in COVID-19 lung tissue compared with controls.

CONCLUSIONS AND RELEVANCE

ANGPTL4 is expressed in pulmonary epithelial cells and fibroblasts and is associated with clinical prognosis in critically ill COVID-19 patients.

Quantifying the phenome-wide disease burden of obesity using electronic health records and genomics

OBJECTIVE

High BMI is associated with many comorbidities and mortality. This study aimed to elucidate the overall clinical risk of obesity using a genome- and phenome-wide approach.

METHODS

This study performed a phenome-wide association study of BMI using a clinical cohort of 736,726 adults. This was followed by genetic association studies using two separate cohorts: one consisting of 65,174 adults in the Electronic Medical Records and Genomics (eMERGE) Network and another with 405,432 participants in the UK Biobank.

RESULTS

Class 3 obesity was associated with 433 phenotypes, representing 59.3% of all billing codes in individuals with severe obesity. A genome-wide polygenic risk score for BMI, accounting for 7.5% of variance in BMI, was associated with 296 clinical diseases, including strong associations with type 2 diabetes, sleep apnea, hypertension, and chronic liver disease. In all three cohorts, 199 phenotypes were associated with class 3 obesity and polygenic risk for obesity, including novel associations such as increased risk of renal failure, venous insufficiency, and gastroesophageal reflux.

CONCLUSIONS

This combined genomic and phenomic systematic approach demonstrated that obesity has a strong genetic predisposition and is associated with a considerable burden of disease across all disease classes.

Genetic regulation of serum IgA levels and susceptibility to common immune, infectious, kidney, and cardio-metabolic traits

Immunoglobulin A (IgA) mediates mucosal responses to food antigens and the intestinal microbiome and is involved in susceptibility to mucosal pathogens, celiac disease, inflammatory bowel disease, and IgA nephropathy. We performed a genome-wide association study of serum IgA levels in 41,263 individuals of diverse ancestries and identified 20 genome-wide significant loci, including 9 known and 11 novel loci. Co-localization analyses with expression QTLs prioritized candidate genes for 14 of 20 significant loci. Most loci encoded genes that produced immune defects and IgA abnormalities when genetically manipulated in mice. We also observed positive genetic correlations of serum IgA levels with IgA nephropathy, type 2 diabetes, and body mass index, and negative correlations with celiac disease, inflammatory bowel disease, and several infections. Mendelian randomization supported elevated serum IgA as a causal factor in IgA nephropathy. African ancestry was consistently associated with higher serum IgA levels and greater frequency of IgA-increasing alleles compared to other ancestries. Our findings provide novel insights into the genetic regulation of IgA levels and its potential role in human disease.

Large-scale genome-wide association study of coronary artery disease in genetically diverse populations

We report a genome-wide association study (GWAS) of coronary artery disease (CAD) incorporating nearly a quarter of a million cases, in which existing studies are integrated with data from cohorts of white, Black and Hispanic individuals from the Million Veteran Program. We document near equivalent heritability of CAD across multiple ancestral groups, identify 95 novel loci, including nine on the X chromosome, detect eight loci of genome-wide significance in Black and Hispanic individuals, and demonstrate that two common haplotypes at the 9p21 locus are responsible for risk stratification in all populations except those of African origin, in which these haplotypes are virtually absent. Moreover, in the largest GWAS for angiographically derived coronary atherosclerosis performed to date, we find 15 loci of genome-wide significance that robustly overlap with established loci for clinical CAD. Phenome-wide association analyses of novel loci and polygenic risk scores (PRSs) augment signals related to insulin resistance, extend pleiotropic associations of these loci to include smoking and family history, and precisely document the markedly reduced transferability of existing PRSs to Black individuals. Downstream integrative analyses reinforce the critical roles of vascular endothelial, fibroblast, and smooth muscle cells in CAD susceptibility, but also point to a shared biology between atherosclerosis and oncogenesis. This study highlights the value of diverse populations in further characterizing the genetic architecture of CAD.

Psychiatric manifestations of rare variation in medically actionable genes: a PheWAS approach

BACKGROUND

As genomic sequencing moves closer to clinical implementation, there has been an increasing acceptance of returning incidental findings to research participants and patients for mutations in highly penetrant, medically actionable genes. A curated list of genes has been recommended by the American College of Medical Genetics and Genomics (ACMG) for return of incidental findings. However, the pleiotropic effects of these genes are not fully known. Such effects could complicate genetic counseling when returning incidental findings. In particular, there has been no systematic evaluation of psychiatric manifestations associated with rare variation in these genes.

RESULTS

Here, we leveraged a targeted sequence panel and real-world electronic health records from the eMERGE network to assess the burden of rare variation in the ACMG-56 genes and two psychiatric-associated genes (CACNA1C  and TCF4) across common mental health conditions in 15,181 individuals of European descent. As a positive control, we showed that this approach replicated the established association between rare mutations in LDLR and hypercholesterolemia with no visible inflation from population stratification. However, we did not identify any genes significantly enriched with rare deleterious variants that confer risk for common psychiatric disorders after correction for multiple testing. Suggestive associations were observed between depression and rare coding variation in PTEN (P = 1.5 × 10^-4), LDLR (P = 3.6 × 10^-4), and CACNA1S (P = 5.8 × 10^-4). We also observed nominal associations between rare variants in KCNQ1 and substance use disorders (P = 2.4 × 10^-4), and APOB and tobacco use disorder (P = 1.1 × 10^-3).

CONCLUSIONS

Our results do not support an association between psychiatric disorders and incidental findings in medically actionable gene mutations, but power was limited with the available sample sizes. Given the phenotypic and genetic complexity of psychiatric phenotypes, future work will require a much larger sequencing dataset to determine whether incidental findings in these genes have implications for risk of psychopathology.

Angiopoietins as Prognostic Markers for Future Kidney Disease and Heart Failure Events after Acute Kidney Injury

BACKGROUND

The mechanisms underlying long-term sequelae after AKI remain unclear. Vessel instability, an early response to endothelial injury, may reflect a shared mechanism and early trigger for CKD and heart failure.

METHODS

To investigate whether plasma angiopoietins, markers of vessel homeostasis, are associated with CKD progression and heart failure admissions after hospitalization in patients with and without AKI, we conducted a prospective cohort study to analyze the balance between angiopoietin-1 (Angpt-1), which maintains vessel stability, and angiopoietin-2 (Angpt-2), which increases vessel destabilization. Three months after discharge, we evaluated the associations between angiopoietins and development of the primary outcomes of CKD progression and heart failure and the secondary outcome of all-cause mortality 3 months after discharge or later.

RESULTS

Median age for the 1503 participants was 65.8 years; 746 (50%) had AKI. Compared with the lowest quartile, the highest quartile of the Angpt-1:Angpt-2 ratio was associated with 72% lower risk of CKD progression (adjusted hazard ratio [aHR], 0.28; 95% confidence interval [CI], 0.15 to 0.51), 94% lower risk of heart failure (aHR, 0.06; 95% CI, 0.02 to 0.15), and 82% lower risk of mortality (aHR, 0.18; 95% CI, 0.09 to 0.35) for those with AKI. Among those without AKI, the highest quartile of Angpt-1:Angpt-2 ratio was associated with 71% lower risk of heart failure (aHR, 0.29; 95% CI, 0.12 to 0.69) and 68% less mortality (aHR, 0.32; 95% CI, 0.15 to 0.68). There were no associations with CKD progression.

CONCLUSIONS

A higher Angpt-1:Angpt-2 ratio was strongly associated with less CKD progression, heart failure, and mortality in the setting of AKI.

Cross-validation of SARS-CoV-2 responses in kidney organoids and clinical populations

Kidneys are critical target organs of COVID-19, but susceptibility and responses to infection remain poorly understood. Here, we combine SARS-CoV-2 variants with genome-edited kidney organoids and clinical data to investigate tropism, mechanism, and therapeutics. SARS-CoV-2 specifically infects organoid proximal tubules among diverse cell types. Infections produce replicating virus, apoptosis, and disrupted cell morphology, features of which are revealed in the context of polycystic kidney disease. Cross-validation of gene expression patterns in organoids reflects proteomic signatures of COVID-19 in the urine of critically ill patients indicating interferon pathway upregulation. SARS-CoV-2 viral variants alpha, beta, gamma, kappa, and delta exhibit comparable levels of infection in organoids. Infection is ameliorated in ACE2-/- organoids and blocked via treatment with de novo-designed spike binder peptides. Collectively, these studies clarify the impact of kidney infection in COVID-19 as reflected in organoids and clinical populations, enabling assessment of viral fitness and emerging therapies.

A genome-wide association study suggests correlations of common genetic variants with peritoneal solute transfer rates in patients with kidney failure receiving peritoneal dialysis

Movement of solutes across the peritoneum allows for the use of peritoneal dialysis to treat kidney failure. However, there is a large inter-individual variability in the peritoneal solute transfer rate (PSTR). Here, we tested the hypothesis that common genetic variants are associated with variability in PSTR. Of the 3561 participants from 69 centers in six countries, 2850 with complete data were included in a genome-wide association study. PSTR was defined as the four-hour dialysate/plasma creatinine ratio from the first peritoneal equilibration test after starting PD. Heritability of PSTR was estimated using genomic-restricted maximum-likelihood analysis, and the association of PSTR with a genome-wide polygenic risk score was also tested. The mean four-hour dialysate/plasma creatinine ratio in participants was 0.70. In 2212 participants of European ancestry, no signal reached genome-wide significance but 23 single nucleotide variants at four loci demonstrated suggestive associations with PSTR. Meta-analysis of ancestry-stratified regressions in 2850 participants revealed five single-nucleotide variants at four loci with suggestive correlations with PSTR. Association across ancestry strata was consistent for rs28644184 at the KDM2B locus. The estimated heritability of PSTR was 19%, and a permuted model polygenic risk score was significantly associated with PSTR. Thus, this genome-wide association study of patients receiving peritoneal dialysis bolsters evidence for a genetic contribution to inter-individual variability in PSTR.

Genome-wide association study of susceptibility to hospitalised respiratory infections

Background: Globally, respiratory infections contribute to significant morbidity and mortality. However, genetic determinants of respiratory infections are understudied and remain poorly understood. Methods: We conducted a genome-wide association study in 19,459 hospitalised respiratory infection cases and 101,438 controls from UK Biobank. We followed-up well-imputed top signals from the UK Biobank discovery analysis in 50,912 respiratory infection cases and 150,442 controls from 11 cohorts. We aggregated effect estimates across studies using inverse variance-weighted meta-analyses. Additionally, we investigated the function of the top signals in order to gain understanding of the underlying biological mechanisms. Results: In the discovery analysis, we report 56 signals at P<5×10-6, one of which was genome-wide significant (P<5×10-8). The genome-wide significant signal was in an intron of PBX3, a gene that encodes pre-B-cell leukaemia transcription factor 3, a homeodomain-containing transcription factor. Further, the genome-wide significant signal was found to colocalise with gene-specific expression quantitative trait loci (eQTLs) affecting expression of PBX3 in lung tissue, where the respiratory infection risk alleles were associated with decreased PBX3 expression in lung tissue, highlighting a possible biological mechanism. Of the 56 signals, 40 were well-imputed in UK Biobank and were investigated in the 11 follow-up cohorts. None of the 40 signals replicated, with effect estimates attenuated. Conclusions: Our discovery analysis implicated PBX3 as a candidate causal gene and suggests a possible role of transcription factor binding activity in respiratory infection susceptibility. However, the PBX3 signal, and the other well-imputed signals, did not replicate when aggregating effect estimates across 11 independent cohorts. Significant phenotypic heterogeneity and differences in study ascertainment may have contributed to this lack of statistical replication. Overall, our study highlighted putative associations and possible biological mechanisms that may provide insight into respiratory infection susceptibility.

Genomic and functional characterization of a mucosal symbiont involved in early-stage colorectal cancer

Colorectal cancer is a major health concern worldwide. Growing evidence for the role of the gut microbiota in the initiation of CRC has sparked interest in approaches that target these microorganisms. However, little is known about the composition and role of the microbiota associated with precancerous polyps. Here, we found distinct microbial signatures between patients with and without polyps and between polyp subtypes using sequencing and culturing techniques. We found a correlation between Bacteroides fragilis recovered and the level of inflammatory cytokines in the mucosa adjacent to the polyp. Additional analysis revealed that B. fragilis from patients with polyps are bft-negative, activate NF-κB through Toll-like receptor 4, induce a pro-inflammatory response, and are enriched in genes associated with LPS biosynthesis. This study provides fundamental insight into the microbial microenvironment of the pre-neoplastic polyp by highlighting strain-specific genomic and proteomic differences, as well as more broad compositional differences in the microbiome.

Penetrance of Breast Cancer Susceptibility Genes From the eMERGE III Network

Background

Unbiased estimates of penetrance are challenging but critically important to make informed choices about strategies for risk management through increased surveillance and risk-reducing interventions.

Methods

We studied the penetrance and clinical outcomes of 7 breast cancer susceptibility genes (BRCA1, BRCA2, TP53, CHEK2, ATM, PALB2, and PTEN) in almost 13 458 participants unselected for personal or family history of breast cancer. We identified 242 female participants with pathogenic or likely pathogenic variants in 1 of the 7 genes for penetrance analyses, and 147 women did not previously know their genetic results.

Results

Out of the 147 women, 32 women were diagnosed with breast cancer at an average age of 52.8 years. Estimated penetrance by age 60 years ranged from 17.8% to 43.8%, depending on the gene. In clinical-impact analysis, 42.3% (95% confidence interval = 31.3% to 53.3%) of women had taken actions related to their genetic results, and 2 new breast cancer cases were identified within the first 12 months after genetic results disclosure.

Conclusions

Our study provides population-based penetrance estimates for the understudied genes CHEK2, ATM, and PALB2 and highlights the importance of using unselected populations for penetrance studies. It also demonstrates the potential clinical impact of genetic testing to improve health care through early diagnosis and preventative screening.

Genome-wide association studies of low back pain and lumbar spinal disorders using electronic health record data identify a locus associated with lumbar spinal stenosis

ABSTRACT

Identifying genetic risk factors for lumbar spine disorders may lead to knowledge regarding underlying mechanisms and the development of new treatments. We conducted a genome-wide association study involving 100,811 participants with genotypes and longitudinal electronic health record data from the Electronic Medical Records and Genomics Network and Geisinger Health. Cases and controls were defined using validated algorithms and clinical diagnostic codes. Electronic health record-defined phenotypes included low back pain requiring healthcare utilization (LBP-HC), lumbosacral radicular syndrome (LSRS), and lumbar spinal stenosis (LSS). Genome-wide association study used logistic regression with additive genetic effects adjusting for age, sex, site-specific factors, and ancestry (principal components). A fixed-effect inverse-variance weighted meta-analysis was conducted. Genetic variants of genome-wide significance (P < 5 × 10-8) were carried forward for replication in an independent sample from UK Biobank. Phenotype prevalence was 48.8% for LBP-HC, 19.8% for LSRS, and 7.9% for LSS. No variants were significantly associated with LBP-HC. One locus was associated with LSRS (lead variant rs146153280:C>G, odds ratio [OR] = 1.17 for G, P = 2.1 × 10-9), but was not replicated. Another locus on chromosome 2 spanning GFPT1, NFU1, and AAK1 was associated with LSS (lead variant rs13427243:G>A, OR = 1.10 for A, P = 4.3 × 10-8) and replicated in UK Biobank (OR = 1.11, P = 5.4 × 10-5). This was the first genome-wide association study meta-analysis of lumbar spinal disorders using electronic health record data. We identified 2 novel associations with LSRS and LSS; the latter was replicated in an independent sample.

Medical Records-Based Genetic Studies of the Complement System

BACKGROUND

Genetic variants in complement genes have been associated with a wide range of human disease states, but well-powered genetic association studies of complement activation have not been performed in large multiethnic cohorts.

METHODS

We performed medical records-based genome-wide and phenome-wide association studies for plasma C3 and C4 levels among participants of the Electronic Medical Records and Genomics (eMERGE) network.

RESULTS

In a GWAS for C3 levels in 3949 individuals, we detected two genome-wide significant loci: chr.1q31.3 (CFH locus; rs3753396-A; β=0.20; 95% CI, 0.14 to 0.25; P=1.52x10-11) and chr.19p13.3 (C3 locus; rs11569470-G; β=0.19; 95% CI, 0.13 to 0.24; P=1.29x10-8). These two loci explained approximately 2% of variance in C3 levels. GWAS for C4 levels involved 3998 individuals and revealed a genome-wide significant locus at chr.6p21.32 (C4 locus; rs3135353-C; β=0.40; 95% CI, 0.34 to 0.45; P=4.58x10-35). This locus explained approximately 13% of variance in C4 levels. The multiallelic copy number variant analysis defined two structural genomic C4 variants with large effect on blood C4 levels: C4-BS (β=-0.36; 95% CI, -0.42 to -0.30; P=2.98x10-22) and C4-AL-BS (β=0.25; 95% CI, 0.21 to 0.29; P=8.11x10-23). Overall, C4 levels were strongly correlated with copy numbers of C4A and C4B genes. In comprehensive phenome-wide association studies involving 102,138 eMERGE participants, we cataloged a full spectrum of autoimmune, cardiometabolic, and kidney diseases genetically related to systemic complement activation.

CONCLUSIONS

We discovered genetic determinants of plasma C3 and C4 levels using eMERGE genomic data linked to electronic medical records. Genetic variants regulating C3 and C4 levels have large effects and multiple clinical correlations across the spectrum of complement-related diseases in humans.

Medical records-based chronic kidney disease phenotype for clinical care and "big data" observational and genetic studies

Chronic Kidney Disease (CKD) represents a slowly progressive disorder that is typically silent until late stages, but early intervention can significantly delay its progression. We designed a portable and scalable electronic CKD phenotype to facilitate early disease recognition and empower large-scale observational and genetic studies of kidney traits. The algorithm uses a combination of rule-based and machine-learning methods to automatically place patients on the staging grid of albuminuria by glomerular filtration rate ("A-by-G" grid). We manually validated the algorithm by 451 chart reviews across three medical systems, demonstrating overall positive predictive value of 95% for CKD cases and 97% for healthy controls. Independent case-control validation using 2350 patient records demonstrated diagnostic specificity of 97% and sensitivity of 87%. Application of the phenotype to 1.3 million patients demonstrated that over 80% of CKD cases are undetected using ICD codes alone. We also demonstrated several large-scale applications of the phenotype, including identifying stage-specific kidney disease comorbidities, in silico estimation of kidney trait heritability in thousands of pedigrees reconstructed from medical records, and biobank-based multicenter genome-wide and phenome-wide association studies.

Copy Number Variant Analysis and Genome-wide Association Study Identify Loci with Large Effect for Vesicoureteral Reflux

BACKGROUND

Vesicoureteral reflux (VUR) is a common, familial genitourinary disorder, and a major cause of pediatric urinary tract infection (UTI) and kidney failure. The genetic basis of VUR is not well understood.

METHODS

A diagnostic analysis sought rare, pathogenic copy number variant (CNV) disorders among 1737 patients with VUR. A GWAS was performed in 1395 patients and 5366 controls, of European ancestry.

RESULTS

Altogether, 3% of VUR patients harbored an undiagnosed rare CNV disorder, such as the 1q21.1, 16p11.2, 22q11.21, and triple X syndromes ((OR, 3.12; 95% CI, 2.10 to 4.54; P=6.35×10^-8) The GWAS identified three study-wide significant and five suggestive loci with large effects (ORs, 1.41-6.9), containing canonical developmental genes expressed in the developing urinary tract (WDPCP, OTX1, BMP5, VANGL1, and WNT5A). In particular, 3.3% of VUR patients were homozygous for an intronic variant in WDPCP (rs13013890; OR, 3.65; 95% CI, 2.39 to 5.56; P=1.86×10^-9). This locus was associated with multiple genitourinary phenotypes in the UK Biobank and eMERGE studies. Analysis of Wnt5a mutant mice confirmed the role of Wnt5a signaling in bladder and ureteric morphogenesis.

CONCLUSIONS

These data demonstrate the genetic heterogeneity of VUR. Altogether, 6% of patients with VUR harbored a rare CNV or a common variant genotype conferring an OR >3. Identification of these genetic risk factors has multiple implications for clinical care and for analysis of outcomes in VUR.

Association between triglycerides, known risk SNVs and conserved rare variation in SLC25A40 in a multi-ancestry cohort

Background

Elevated triglycerides (TG) are associated with, and may be causal for, cardiovascular disease (CVD), and co-morbidities such as type II diabetes and metabolic syndrome. Pathogenic variants in APOA5 and APOC3 as well as risk SNVs in other genes [APOE (rs429358, rs7412), APOA1/C3/A4/A5 gene cluster (rs964184), INSR (rs7248104), CETP (rs7205804), GCKR (rs1260326)] have been shown to affect TG levels. Knowledge of genetic causes for elevated TG may lead to early intervention and targeted treatment for CVD. We previously identified linkage and association of a rare, highly conserved missense variant in SLC25A40, rs762174003, with hypertriglyceridemia (HTG) in a single large family, and replicated this association with rare, highly conserved missense variants in a European American and African American sample.

Methods

Here, we analyzed a longitudinal mixed-ancestry cohort (European, African and Asian ancestry, N = 8966) from the Electronic Medical Record and Genomics (eMERGE) Network. We tested associations between median TG and the genes of interest, using linear regression, adjusting for sex, median age, median BMI, and the first two principal components of ancestry.

Results

We replicated the association between TG and APOC3, APOA5, and risk variation at APOE, APOA1/C3/A4/A5 gene cluster, and GCKR. We failed to replicate the association between rare, highly conserved variation at SLC25A40 and TG, as well as for risk variation at INSR and CETP.

Conclusions

Analysis using data from electronic health records presents challenges that need to be overcome. Although large amounts of genotype data is becoming increasingly accessible, usable phenotype data can be challenging to obtain. We were able to replicate known, strong associations, but were unable to replicate moderate associations due to the limited sample size and missing drug information.

Hemochromatosis risk genotype is not associated with colorectal cancer or age at its diagnosis

Homozygotes for the higher penetrance hemochromatosis risk allele, HFE c.845G>A (p.Cys282Tyr, or C282Y), have been reported to be at a 2- to 3-fold increased risk for colorectal cancer (CRC). These results have been reported for small sample size studies with no information about age at diagnosis for CRC. An association with age at diagnosis might alter CRC screening recommendations. We analyzed two large European ancestry datasets to assess the association of HFE genotype with CRC risk and age at CRC diagnosis. The first dataset included 59,733 CRC or advanced adenoma cases and 72,351 controls from a CRC epidemiological study consortium. The second dataset included 13,564 self-reported CRC cases and 2,880,218 controls from the personal genetics company, 23andMe. No association of the common hereditary hemochromatosis (HH) risk genotype and CRC was found in either dataset. The odds ratios (ORs) for the association of CRC and HFE C282Y homozygosity were 1.08 (95% confidence interval [CI], 0.91-1.29; p = 0.4) and 1.01 (95% CI, 0.78-1.31, p = 0.9) in the two cohorts, respectively. Age at CRC diagnosis also did not differ by HFE C282Y/C282Y genotype in either dataset. These results indicate no increased CRC risk in individuals with HH genotypes and suggest that persons with HH risk genotypes can follow population screening recommendations for CRC.

Loci identified by a genome-wide association study of carotid artery stenosis in the eMERGE network

Carotid artery atherosclerotic disease (CAAD) is a risk factor for stroke. We used a genome-wide association (GWAS) approach to discover genetic variants associated with CAAD in participants in the electronic Medical Records and Genomics (eMERGE) Network. We identified adult CAAD cases with unilateral or bilateral carotid artery stenosis and controls without evidence of stenosis from electronic health records at eight eMERGE sites. We performed GWAS with a model adjusting for age, sex, study site, and genetic principal components of ancestry. In eMERGE we found 1793 CAAD cases and 17,958 controls. Two loci reached genome-wide significance, on chr6 in LPA (rs10455872, odds ratio [OR] (95% confidence interval [CI]) = 1.50 (1.30-1.73), p = 2.1 × 10^-8 ) and on chr7, an intergenic single nucleotide variant (SNV; rs6952610, OR (95% CI) = 1.25 (1.16-1.36), p = 4.3 × 10^-8 ). The chr7 association remained significant in the presence of the LPA SNV as a covariate. The LPA SNV was also associated with coronary heart disease (CHD; 4199 cases and 11,679 controls) in this study (OR (95% CI) = 1.27 (1.13-1.43), p = 5 × 10^-5 ) but the chr7 SNV was not (OR (95% CI) = 1.03 (0.97-1.09), p = .37). Both variants replicated in UK Biobank. Elevated lipoprotein(a) concentrations ([Lp(a)]) and LPA variants associated with elevated [Lp(a)] have previously been associated with CAAD and CHD, including rs10455872. With electronic health record phenotypes in eMERGE and UKB, we replicated a previously known association and identified a novel locus associated with CAAD.

Genome-wide Modeling of Polygenic Risk Score in Colorectal Cancer Risk

Accurate colorectal cancer (CRC) risk prediction models are critical for identifying individuals at low and high risk of developing CRC, as they can then be offered targeted screening and interventions to address their risks of developing disease (if they are in a high-risk group) and avoid unnecessary screening and interventions (if they are in a low-risk group). As it is likely that thousands of genetic variants contribute to CRC risk, it is clinically important to investigate whether these genetic variants can be used jointly for CRC risk prediction. In this paper, we derived and compared different approaches to generating predictive polygenic risk scores (PRS) from genome-wide association studies (GWASs) including 55,105 CRC-affected case subjects and 65,079 control subjects of European ancestry. We built the PRS in three ways, using (1) 140 previously identified and validated CRC loci; (2) SNP selection based on linkage disequilibrium (LD) clumping followed by machine-learning approaches; and (3) LDpred, a Bayesian approach for genome-wide risk prediction. We tested the PRS in an independent cohort of 101,987 individuals with 1,699 CRC-affected case subjects. The discriminatory accuracy, calculated by the age- and sex-adjusted area under the receiver operating characteristics curve (AUC), was highest for the LDpred-derived PRS (AUC = 0.654) including nearly 1.2 M genetic variants (the proportion of causal genetic variants for CRC assumed to be 0.003), whereas the PRS of the 140 known variants identified from GWASs had the lowest AUC (AUC = 0.629). Based on the LDpred-derived PRS, we are able to identify 30% of individuals without a family history as having risk for CRC similar to those with a family history of CRC, whereas the PRS based on known GWAS variants identified only top 10% as having a similar relative risk. About 90% of these individuals have no family history and would have been considered average risk under current screening guidelines, but might benefit from earlier screening. The developed PRS offers a way for risk-stratified CRC screening and other targeted interventions.

A Polygenic and Phenotypic Risk Prediction for Polycystic Ovary Syndrome Evaluated by Phenome-Wide Association Studies

CONTEXT

As many as 75% of patients with polycystic ovary syndrome (PCOS) are estimated to be unidentified in clinical practice.

OBJECTIVE

Utilizing polygenic risk prediction, we aim to identify the phenome-wide comorbidity patterns characteristic of PCOS to improve accurate diagnosis and preventive treatment.

DESIGN, PATIENTS, AND METHODS

Leveraging the electronic health records (EHRs) of 124 852 individuals, we developed a PCOS risk prediction algorithm by combining polygenic risk scores (PRS) with PCOS component phenotypes into a polygenic and phenotypic risk score (PPRS). We evaluated its predictive capability across different ancestries and perform a PRS-based phenome-wide association study (PheWAS) to assess the phenomic expression of the heightened risk of PCOS.

RESULTS

The integrated polygenic prediction improved the average performance (pseudo-R2) for PCOS detection by 0.228 (61.5-fold), 0.224 (58.8-fold), 0.211 (57.0-fold) over the null model across European, African, and multi-ancestry participants respectively. The subsequent PRS-powered PheWAS identified a high level of shared biology between PCOS and a range of metabolic and endocrine outcomes, especially with obesity and diabetes: "morbid obesity", "type 2 diabetes", "hypercholesterolemia", "disorders of lipid metabolism", "hypertension", and "sleep apnea" reaching phenome-wide significance.

CONCLUSIONS

Our study has expanded the methodological utility of PRS in patient stratification and risk prediction, especially in a multifactorial condition like PCOS, across different genetic origins. By utilizing the individual genome-phenome data available from the EHR, our approach also demonstrates that polygenic prediction by PRS can provide valuable opportunities to discover the pleiotropic phenomic network associated with PCOS pathogenesis.

Association of Genetic Risk of Obesity with Postoperative Complications Using Mendelian Randomization

BACKGROUND

The extent to which obesity and genetics determine postoperative complications is incompletely understood.

METHODS

We performed a retrospective study using two population cohorts with electronic health record (EHR) data. The first included 736,726 adults with body mass index (BMI) recorded between 1990 and 2017 at Vanderbilt University Medical Center. The second cohort consisted of 65,174 individuals from 12 institutions contributing EHR and genome-wide genotyping data to the Electronic Medical Records and Genomics (eMERGE) Network. Pairwise logistic regression analyses were used to measure the association of BMI categories with postoperative complications derived from International Classification of Disease-9 codes, including postoperative infection, incisional hernia, and intestinal obstruction. A genetic risk score was constructed from 97 obesity-risk single-nucleotide polymorphisms for a Mendelian randomization study to determine the association of genetic risk of obesity on postoperative complications. Logistic regression analyses were adjusted for sex, age, site, and race/principal components.

RESULTS

Individuals with overweight or obese BMI (≥25 kg/m2) had increased risk of incisional hernia (odds ratio [OR] 1.7-5.5, p < 3.1 × 10-20), and people with obesity (BMI ≥ 30 kg/m2) had increased risk of postoperative infection (OR 1.2-2.3, p < 2.5 × 10-5). In the eMERGE cohort, genetically predicted BMI was associated with incisional hernia (OR 2.1 [95% CI 1.8-2.5], p = 1.4 × 10-6) and postoperative infection (OR 1.6 [95% CI 1.4-1.9], p = 3.1 × 10-6). Association findings were similar after limitation of the cohorts to those who underwent abdominal procedures.

CONCLUSIONS

Clinical and Mendelian randomization studies suggest that obesity, as measured by BMI, is associated with the development of postoperative incisional hernia and infection.

Anchoring a dynamic in vitro model of human neuronal differentiation to key processes of early brain development in vivo

We characterize temporal pathway dynamics of differentiation in an in vitro neurotoxicity model with the aim of informing design and interpretation of toxicological assays. Human neural progenitor cells (hNPCs) were cultured in differentiation conditions up to 21 days. Genes significantly changed through time were identified and grouped according to temporal dynamics. Quantitative pathway analysis identified gene ontology (GO) terms enriched among significantly changed genes and provided a temporal roadmap of pathway trends in vitro. Gene expression in hNPCs was compared with publicly available gene expression data from developing human brain tissue in vivo. Quantitative pathway analysis of significantly changed genes and targeted analysis of specific pathways of interest identified concordance between in vivo and in vitro expression associated with proliferation, migration, differentiation, synapse formation, and neurotransmission. Our analysis anchors gene expression patterns in vitro to sensitive windows of in vivo development, helping to define appropriate applications of the model.

Harmonizing Clinical Sequencing and Interpretation for the eMERGE III Network

The advancement of precision medicine requires new methods to coordinate and deliver genetic data from heterogeneous sources to physicians and patients. The eMERGE III Network enrolled >25,000 participants from biobank and prospective cohorts of predominantly healthy individuals for clinical genetic testing to determine clinically actionable findings. The network developed protocols linking together the 11 participant collection sites and 2 clinical genetic testing laboratories. DNA capture panels targeting 109 genes were used for testing of DNA and sample collection, data generation, interpretation, reporting, delivery, and storage were each harmonized. A compliant and secure network enabled ongoing review and reconciliation of clinical interpretations, while maintaining communication and data sharing between clinicians and investigators. A total of 202 individuals had positive diagnostic findings relevant to the indication for testing and 1,294 had additional/secondary findings of medical significance deemed to be returnable, establishing data return rates for other testing endeavors. This study accomplished integration of structured genomic results into multiple electronic health record (EHR) systems, setting the stage for clinical decision support to enable genomic medicine. Further, the established processes enable different sequencing sites to harmonize technical and interpretive aspects of sequencing tests, a critical achievement toward global standardization of genomic testing. The eMERGE protocols and tools are available for widespread dissemination.

GWAS and enrichment analyses of non-alcoholic fatty liver disease identify new trait-associated genes and pathways across eMERGE Network

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) is a common chronic liver illness with a genetically heterogeneous background that can be accompanied by considerable morbidity and attendant health care costs. The pathogenesis and progression of NAFLD is complex with many unanswered questions. We conducted genome-wide association studies (GWASs) using both adult and pediatric participants from the Electronic Medical Records and Genomics (eMERGE) Network to identify novel genetic contributors to this condition.

METHODS

First, a natural language processing (NLP) algorithm was developed, tested, and deployed at each site to identify 1106 NAFLD cases and 8571 controls and histological data from liver tissue in 235 available participants. These include 1242 pediatric participants (396 cases, 846 controls). The algorithm included billing codes, text queries, laboratory values, and medication records. Next, GWASs were performed on NAFLD cases and controls and case-only analyses using histologic scores and liver function tests adjusting for age, sex, site, ancestry, PC, and body mass index (BMI).

RESULTS

Consistent with previous results, a robust association was detected for the PNPLA3 gene cluster in participants with European ancestry. At the PNPLA3-SAMM50 region, three SNPs, rs738409, rs738408, and rs3747207, showed strongest association (best SNP rs738409 p = 1.70 × 10- 20). This effect was consistent in both pediatric (p = 9.92 × 10- 6) and adult (p = 9.73 × 10- 15) cohorts. Additionally, this variant was also associated with disease severity and NAFLD Activity Score (NAS) (p = 3.94 × 10- 8, beta = 0.85). PheWAS analysis link this locus to a spectrum of liver diseases beyond NAFLD with a novel negative correlation with gout (p = 1.09 × 10- 4). We also identified novel loci for NAFLD disease severity, including one novel locus for NAS score near IL17RA (rs5748926, p = 3.80 × 10- 8), and another near ZFP90-CDH1 for fibrosis (rs698718, p = 2.74 × 10- 11). Post-GWAS and gene-based analyses identified more than 300 genes that were used for functional and pathway enrichment analyses.

CONCLUSIONS

In summary, this study demonstrates clear confirmation of a previously described NAFLD risk locus and several novel associations. Further collaborative studies including an ethnically diverse population with well-characterized liver histologic features of NAFLD are needed to further validate the novel findings.

Detecting potential pleiotropy across cardiovascular and neurological diseases using univariate, bivariate, and multivariate methods on 43,870 individuals from the eMERGE network

The link between cardiovascular diseases and neurological disorders has been widely observed in the aging population. Disease prevention and treatment rely on understanding the potential genetic nexus of multiple diseases in these categories. In this study, we were interested in detecting pleiotropy, or the phenomenon in which a genetic variant influences more than one phenotype. Marker-phenotype association approaches can be grouped into univariate, bivariate, and multivariate categories based on the number of phenotypes considered at one time. Here we applied one statistical method per category followed by an eQTL colocalization analysis to identify potential pleiotropic variants that contribute to the link between cardiovascular and neurological diseases. We performed our analyses on ~530,000 common SNPs coupled with 65 electronic health record (EHR)-based phenotypes in 43,870 unrelated European adults from the Electronic Medical Records and Genomics (eMERGE) network. There were 31 variants identified by all three methods that showed significant associations across late onset cardiac- and neurologic- diseases. We further investigated functional implications of gene expression on the detected "lead SNPs" via colocalization analysis, providing a deeper understanding of the discovered associations. In summary, we present the framework and landscape for detecting potential pleiotropy using univariate, bivariate, multivariate, and colocalization methods. Further exploration of these potentially pleiotropic genetic variants will work toward understanding disease causing mechanisms across cardiovascular and neurological diseases and may assist in considering disease prevention as well as drug repositioning in future research.

The eMERGE genotype set of 83,717 subjects imputed to ~40 million variants genome wide and association with the herpes zoster medical record phenotype

The Electronic Medical Records and Genomics (eMERGE) network is a network of medical centers with electronic medical records linked to existing biorepository samples for genomic discovery and genomic medicine research. The network sought to unify the genetic results from 78 Illumina and Affymetrix genotype array batches from 12 contributing medical centers for joint association analysis of 83,717 human participants. In this report, we describe the imputation of eMERGE results and methods to create the unified imputed merged set of genome‐wide variant genotype data. We imputed the data using the Michigan Imputation Server, which provides a missing single‐nucleotide variant genotype imputation service using the minimac3 imputation algorithm with the Haplotype Reference Consortium genotype reference set. We describe the quality control and filtering steps used in the generation of this data set and suggest generalizable quality thresholds for imputation and phenotype association studies. To test the merged imputed genotype set, we replicated a previously reported chromosome 6 HLA‐B herpes zoster (shingles) association and discovered a novel zoster‐associated loci in an epigenetic binding site near the terminus of chromosome 3 (3p29).

Rare loss of function variants in candidate genes and risk of colorectal cancer

Although ~ 25% of colorectal cancer or polyp (CRC/P) cases show familial aggregation, current germline genetic testing identifies a causal genotype in the 16 major genes associated with high penetrance CRC/P in only 20% of these cases. As there are likely other genes underlying heritable CRC/P, we evaluated the association of variation at novel loci with CRC/P. We evaluated 158 a priori selected candidate genes by comparing the number of rare potentially disruptive variants (PDVs) found in 84 CRC/P cases without an identified CRC/P risk-associated variant and 2440 controls. We repeated this analysis using an additional 73 CRC/P cases. We also compared the frequency of PDVs in select genes among CRC/P cases with two publicly available data sets. We found a significant enrichment of PDVs in cases vs. controls: 20% of cases vs. 11.5% of controls with ≥ 1 PDV (OR = 1.9, p = 0.01) in the original set of cases. Among the second cohort of CRC/P cases, 18% had a PDV, significantly different from 11.5% (p = 0.02). Logistic regression, adjusting for ancestry and multiple testing, indicated association between CRC/P and PDVs in NTHL1 (p = 0.0001), BRCA2 (p = 0.01) and BRIP1 (p = 0.04). However, there was no significant difference in the frequency of PDVs at each of these genes between all 157 CRC/P cases and two publicly available data sets. These results suggest an increased presence of PDVs in CRC/P cases and support further investigation of the association of NTHL1, BRCA2 and BRIP1 variation with CRC/P.

An Atlas of Genetic Variation Linking Pathogen-Induced Cellular Traits to Human Disease

Pathogens have been a strong driving force for natural selection. Therefore, understanding how human genetic differences impact infection-related cellular traits can mechanistically link genetic variation to disease susceptibility. Here we report the Hi-HOST Phenome Project (H2P2): a catalog of cellular genome-wide association studies (GWAS) comprising 79 infection-related phenotypes in response to 8 pathogens in 528 lymphoblastoid cell lines. Seventeen loci surpass genome-wide significance for infection-associated phenotypes ranging from pathogen replication to cytokine production. We combined H2P2 with clinical association data from patients to identify a SNP near CXCL10 as a risk factor for inflammatory bowel disease. A SNP in the transcriptional repressor ZBTB20 demonstrated pleiotropy, likely through suppression of multiple target genes, and was associated with viral hepatitis. These data are available on a web portal to facilitate interpreting human genome variation through the lens of cell biology and should serve as a rich resource for the research community.

Association of rare missense variants in the second intracellular loop of NaV1.7 sodium channels with familial autism

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder often accompanied by intellectual disability, language impairment and medical co-morbidities. The heritability of autism is high and multiple genes have been implicated as causal. However, most of these genes have been identified in de novo cases. To further the understanding of familial autism, we performed whole-exome sequencing on five families in which second- and third-degree relatives were affected. By focusing on novel and protein-altering variants, we identified a small set of candidate genes. Among these, a novel private missense C1143F variant in the second intracellular loop of the voltage-gated sodium channel Na_V1.7, encoded by the SCN9A gene, was identified in one family. Through electrophysiological analysis, we show that Na_V1.7^C1143F exhibits partial loss-of-function effects, resulting in slower recovery from inactivation and decreased excitability in cultured cortical neurons. Furthermore, for the same intracellular loop of Na_V1.7, we found an excess of rare variants in a case-control variant-burden study. Functional analysis of one of these variants, M932L/V991L, also demonstrated reduced firing in cortical neurons. However, although this variant is rare in Caucasians, it is frequent in Latino population, suggesting that genetic background can alter its effects on phenotype. Although the involvement of the SCN1A and SCN2A genes encoding Na_V1.1 and Na_V1.2 channels in de novo ASD has previously been demonstrated, our study indicates the involvement of inherited SCN9A variants and partial loss-of-function of Na_V1.7 channels in the etiology of rare familial ASD.

Burden of potentially pathologic copy number variants is higher in children with isolated congenital heart disease and significantly impairs covariate-adjusted transplant-free survival

OBJECTIVES

Copy number variants (CNVs) are duplications or deletions of genomic regions. Large CNVs are potentially pathogenic and are overrepresented in children with congenital heart disease (CHD). We sought to determine the frequency of large CNVs in children with isolated CHD, and to evaluate the relationship of these potentially pathogenic CNVs with transplant-free survival.

METHODS

These cases are derived from a prospective cohort of patients with nonsyndromic CHD (n = 422) identified before first surgery. Healthy pediatric controls (n = 500) were obtained from the electronic Medical Records and Genetic Epidemiology Network, and CNV frequency was contrasted for CHD cases and controls. CNVs were determined algorithmically; subsequently screened for >95% overlap between 2 methods, size (>300 kb), quality score, overlap with a gene, and novelty (absent from databases of known, benign CNVs); and separately validated by quantitative polymerase chain reaction. Survival likelihoods for cases were calculated using Cox proportional hazards modeling to evaluate the joint effect of CNV burden and known confounders on transplant-free survival.

RESULTS

Children with nonsyndromic CHD had a higher burden of potentially pathogenic CNVs compared with pediatric controls (12.1% vs 5.0%; P = .00016). Presence of a CNV was associated with significantly decreased transplant-free survival after surgery (hazard ratio, 3.42; 95% confidence interval, 1.66-7.09; P = .00090) with confounder adjustment.

CONCLUSIONS

We confirm that children with isolated CHD have a greater burden of rare/large CNVs. We report a novel finding that these CNVs are associated with an adjusted 2.55-fold increased risk of death or transplant. These data suggest that CNV burden is an important modifier of survival after surgery for CHD.

Patient genotypes impact survival after surgery for isolated congenital heart disease

BACKGROUND

Survival after cardiac surgery in infancy requires adaptive responses from oxidative stress management and vascular regulation pathways. We tested the hypothesis that genetic variation in these pathways influences postoperative survival in nonsyndromic congenital heart disease children.

METHODS

This is an analysis of a cohort of nonsyndromic congenital heart disease patients who underwent cardiac surgery with cardiopulmonary bypass before 6 months of age (n=422). Six single nucleotide polymorphisms (SNPs) in six genes involved in oxidative stress and vascular response pathways, identified through a priori literature search, were tested for effects on transplant-free survival. Survival curves, adjusting for confounding covariates, were calculated using the Cox proportional hazard models.

RESULTS

Long-term survival was strongly associated with vascular endothelial growth factor A gene SNP rs833069 (p=7.03×10(-4)) and superoxide dismutase 2 gene SNP rs2758331 (p=0.019). To test for joint effects of the two SNPs on transplant-free survival, the genotypes were grouped to form a risk score reflecting the cumulative number of risk alleles (0 to 4 alleles per patient). A higher risk score based on the VEGFA and SOD2 SNP genotypes was associated with worse transplant-free survival (p=3.02×10(-4)) after confounder adjustment. The total burden of risk alleles was additive; subjects with the highest risk score of 4 (n=59 subjects, 14.2% of the cohort) had a total covariate-adjusted hazard ratio of 15.64 for worse transplant-free survival.

CONCLUSIONS

After cardiac surgery, infants who are homozygous for the high-risk alleles for both the VEGFA and SOD2 SNPs have an approximately 16-fold increased risk of death or heart transplant, suggesting that genetic variants are important modifiers of survival after surgery for congenital heart disease.

Identification of rare variants from exome sequence in a large pedigree with autism

We carried out analyses with the goal of identifying rare variants in exome sequence data that contribute to disease risk for a complex trait. We analyzed a large, 47-member, multigenerational pedigree with 11 cases of autism spectrum disorder, using genotypes from 3 technologies representing increasing resolution: a multiallelic linkage marker panel, a dense diallelic marker panel, and variants from exome sequencing. Genome-scan marker genotypes were available on most subjects, and exome sequence data was available on 5 subjects. We used genome-scan linkage analysis to identify and prioritize the chromosome 22 region of interest, and to select subjects for exome sequencing. Inheritance vectors (IVs) generated by Markov chain Monte Carlo analysis of multilocus marker data were the foundation of most analyses. Genotype imputation used IVs to determine which sequence variants reside on the haplotype that co-segregates with the autism diagnosis. Together with a rare-allele frequency filter, we identified only one rare variant on the risk haplotype, illustrating the potential of this approach to prioritize variants. The associated gene, MYH9, is biologically unlikely, and we speculate that for this complex trait, the key variants may lie outside the exome.

Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders

Exome sequencing studies of autism spectrum disorders (ASDs) have identified many de novo mutations but few recurrently disrupted genes. We therefore developed a modified molecular inversion probe method enabling ultra-low-cost candidate gene resequencing in very large cohorts. To demonstrate the power of this approach, we captured and sequenced 44 candidate genes in 2446 ASD probands. We discovered 27 de novo events in 16 genes, 59% of which are predicted to truncate proteins or disrupt splicing. We estimate that recurrent disruptive mutations in six genes-CHD8, DYRK1A, GRIN2B, TBR1, PTEN, and TBL1XR1-may contribute to 1% of sporadic ASDs. Our data support associations between specific genes and reciprocal subphenotypes (CHD8-macrocephaly and DYRK1A-microcephaly) and replicate the importance of a β-catenin-chromatin-remodeling network to ASD etiology.

Results of genome-wide analyses on neurodevelopmental phenotypes at four-year follow-up following cardiac surgery in infancy

BACKGROUND

Adverse neurodevelopmental sequelae are reported among children who undergo early cardiac surgery to repair congenital heart defects (CHD). APOE genotype has previously been determined to contribute to the prediction of these outcomes. Understanding further genetic causes for the development of poor neurobehavioral outcomes should enhance patient risk stratification and improve both prevention and treatment strategies.

METHODS

We performed a prospective observational study of children who underwent cardiac surgery before six months of age; this included a neurodevelopmental evaluation between their fourth and fifth birthdays. Attention and behavioral skills were assessed through parental report utilizing the Attention Deficit-Hyperactivity Disorder-IV scale preschool edition (ADHD-IV), and Child Behavior Checklist (CBCL/1.5-5), respectively. Of the seven investigated, three neurodevelopmental phenotypes met genomic quality control criteria. Linear regression was performed to determine the effect of genome-wide genetic variation on these three neurodevelopmental measures in 316 subjects.

RESULTS

This genome-wide association study identified single nucleotide polymorphisms (SNPs) associated with three neurobehavioral phenotypes in the postoperative children ADHD-IV Impulsivity/Hyperactivity, CBCL/1.5-5 PDPs, and CBCL/1.5-5 Total Problems. The most predictive SNPs for each phenotype were: a LGALS8 intronic SNP, rs4659682, associated with ADHD-IV Impulsivity (P=1.03 × 10(-6)); a PCSK5 intronic SNP, rs2261722, associated with CBCL/1.5-5 PDPs (P=1.11 × 10(-6)); and an intergenic SNP, rs11617488, 50 kb from FGF9, associated with CBCL/1.5-5 Total Problems (P=3.47 × 10(-7)). 10 SNPs (3 for ADHD-IV Impulsivity, 5 for CBCL/1.5-5 PDPs, and 2 for CBCL/1.5-5 Total Problems) had p<10(-5).

CONCLUSIONS

No SNPs met genome-wide significance for our three neurobehavioral phenotypes; however, 10 SNPs reached a threshold for suggestive significance (p<10(-5)). Given the unique nature of this cohort, larger studies and/or replication are not possible. Studies to further investigate the mechanisms through which these newly identified genes may influence neurodevelopment dysfunction are warranted.

Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations

It is well established that autism spectrum disorders (ASD) have a strong genetic component; however, for at least 70% of cases, the underlying genetic cause is unknown. Under the hypothesis that de novo mutations underlie a substantial fraction of the risk for developing ASD in families with no previous history of ASD or related phenotypes--so-called sporadic or simplex families--we sequenced all coding regions of the genome (the exome) for parent-child trios exhibiting sporadic ASD, including 189 new trios and 20 that were previously reported. Additionally, we also sequenced the exomes of 50 unaffected siblings corresponding to these new (n = 31) and previously reported trios (n = 19), for a total of 677 individual exomes from 209 families. Here we show that de novo point mutations are overwhelmingly paternal in origin (4:1 bias) and positively correlated with paternal age, consistent with the modest increased risk for children of older fathers to develop ASD. Moreover, 39% (49 of 126) of the most severe or disruptive de novo mutations map to a highly interconnected β-catenin/chromatin remodelling protein network ranked significantly for autism candidate genes. In proband exomes, recurrent protein-altering mutations were observed in two genes: CHD8 and NTNG1. Mutation screening of six candidate genes in 1,703 ASD probands identified additional de novo, protein-altering mutations in GRIN2B, LAMC3 and SCN1A. Combined with copy number variant (CNV) data, these results indicate extreme locus heterogeneity but also provide a target for future discovery, diagnostics and therapeutics.

IGF1R variants associated with isolated single suture craniosynostosis

The genetic contribution to the pathogenesis of isolated single suture craniosynostosis is poorly understood. The role of mutations in genes known to be associated with syndromic synostosis appears to be limited. We present our findings of a candidate gene resequencing approach to identify rare variants associated with the most common forms of isolated craniosynostosis. Resequencing of the coding regions, splice junction sites, and 5' and 3' untranslated regions of 27 candidate genes in 186 cases of isolated non-syndromic single suture synostosis revealed three novel and two rare sequence variants (R406H, R595H, N857S, P190S, M446V) in insulin-like growth factor I receptor (IGF1R) that are enriched relative to control samples. Mapping the resultant amino acid changes to the modeled homodimer protein structure suggests a structural basis for segregation between these and other disease-associated mutations found in IGF1R. These data suggest that IGF1R mutations may contribute to the risk and in some cases cause single suture craniosynostosis.

Identification, replication, and functional fine-mapping of expression quantitative trait loci in primary human liver tissue

The discovery of expression quantitative trait loci ("eQTLs") can help to unravel genetic contributions to complex traits. We identified genetic determinants of human liver gene expression variation using two independent collections of primary tissue profiled with Agilent (n = 206) and Illumina (n = 60) expression arrays and Illumina SNP genotyping (550K), and we also incorporated data from a published study (n = 266). We found that ∼30% of SNP-expression correlations in one study failed to replicate in either of the others, even at thresholds yielding high reproducibility in simulations, and we quantified numerous factors affecting reproducibility. Our data suggest that drug exposure, clinical descriptors, and unknown factors associated with tissue ascertainment and analysis have substantial effects on gene expression and that controlling for hidden confounding variables significantly increases replication rate. Furthermore, we found that reproducible eQTL SNPs were heavily enriched near gene starts and ends, and subsequently resequenced the promoters and 3'UTRs for 14 genes and tested the identified haplotypes using luciferase assays. For three genes, significant haplotype-specific in vitro functional differences correlated directly with expression levels, suggesting that many bona fide eQTLs result from functional variants that can be mechanistically isolated in a high-throughput fashion. Finally, given our study design, we were able to discover and validate hundreds of liver eQTLs. Many of these relate directly to complex traits for which liver-specific analyses are likely to be relevant, and we identified dozens of potential connections with disease-associated loci. These included previously characterized eQTL contributors to diabetes, drug response, and lipid levels, and they suggest novel candidates such as a role for NOD2 expression in leprosy risk and C2orf43 in prostate cancer. In general, the work presented here will be valuable for future efforts to precisely identify and functionally characterize genetic contributions to a variety of complex traits.

Single nucleotide polymorphism discovery in TBX1 in individuals with and without 22q11.2 deletion syndrome

BACKGROUND

Children with 22q11.2 deletion syndrome (22q11.2DS) have a wide range of clinical features. TBX1 has been proposed as a candidate gene for some of the features in this condition. Polymorphisms in the nondeleted TBX1, which may affect the function of the sole TBX1 gene in individuals with the 22q11.2DS, may be a key to understanding the phenotypic variability among individuals with a shared deletion. Comprehensive single nucleotide polymorphism (SNP) discovery by resequencing candidate genes can identify genetic variants that influence a given phenotype. The purpose of this study was to further characterize the sequence variability in TBX1 by identifying all common SNPs in this gene.

METHODS

We resequenced TBX1 in 29 children with a documented 22q11.2 deletion and 95 nondeleted, healthy individuals. We estimated allele frequencies, performed tagSNP selection, and inferred haplotypes. We also compared SNP frequencies between 22q11.2DS and control samples.

RESULTS

We identified 355 biallelic markers among the 190 chromosomes resequenced in the control panel. The vast majority of the markers identified were SNPs (n = 331), and the remainder indels (n = 24). We did not identify SNPs or indels in the cis- regulatory element (FOX-binding site) upstream of TBX1. In children with 22q11.2DS we detected 187 biallelic markers, six of which were indels. Four of the seven coding SNPs identified in the controls were identified in children with 22q11.2DS.

CONCLUSIONS

This comprehensive SNP discovery data can be used to select SNPs to genotype for future association studies assessing the role of TBX1 and phenotypic variability in individuals with 22q11.2DS.

The environmental genome project: reference polymorphisms for drug metabolism genes and genome-wide association studies

The Environmental Genome Project (EGP) has generated a comprehensive resource of commonly occurring single nucleotide polymorphisms (SNPs) in more than six hundred environmental response genes, including a number of drug metabolism genes. The gene-oriented sequence variation discovery carried out by the EGP is complementary to the HapMap and enables genome-wide association studies (GWAS) that survey a large portion of the known common variation. For GWAS focused on drug metabolism genes and phenotypes, it is important to know the proportion of common SNPs covered by the commercially available high-throughput genotyping chips. Herein, we review a subset of Phase I cytochrome P450 genes studied by the EGP, approaches to GWAS, and the sensitivity of available genotyping platforms to capture common sequence variation in this subset of drug metabolism genes.

Direct detection of null alleles in SNP genotyping data

Pinpointing genetic associations in the human genome relies heavily on the accuracy of the underlying genotype data. Null alleles can generate significant inaccuracies in genotype data and can negatively affect the statistical power of a study. Existing quality control (QC) tests, including tests of Hardy-Weinberg equilibrium, are not sensitive enough to detect the presence of even moderately frequent null alleles in the data. We show that direct analysis of raw data from a quantitative genotyping platform can detect up to 75% of null alleles, even at frequencies below the sensitivity of more traditional methods. Detecting unexpected null alleles not only has benefits in QC of genotype data but may also be valuable in detecting rare, functional null alleles that would otherwise be missed.