Bedside Back to Bench: Building Bridges between Basic and Clinical Genomic Research.

Genome sequencing has revolutionized the diagnosis of genetic diseases. Close collaborations between basic scientists and clinical genomicists are now needed to link genetic variants with disease causation. To facilitate such collaborations, we recommend prioritizing clinically relevant genes for functional studies, developing reference variant-phenotype databases, adopting phenotype description standards, and promoting data sharing.

A nondegenerate code of deleterious variants in Mendelian loci contributes to complex disease risk.

Although countless highly penetrant variants have been associated with Mendelian disorders, the genetic etiologies underlying complex diseases remain largely unresolved. By mining the medical records of over 110 million patients, we examine the extent to which Mendelian variation contributes to complex disease risk. We detect thousands of associations between Mendelian and complex diseases, revealing a nondegenerate, phenotypic code that links each complex disorder to a unique collection of Mendelian loci. Using genome-wide association results, we demonstrate that common variants associated with complex diseases are enriched in the genes indicated by this "Mendelian code." Finally, we detect hundreds of comorbidity associations among Mendelian disorders, and we use probabilistic genetic modeling to demonstrate that Mendelian variants likely contribute nonadditively to the risk for a subset of complex diseases. Overall, this study illustrates a complementary approach for mapping complex disease loci and provides unique predictions concerning the etiologies of specific diseases.

A brief history of human disease genetics.

A primary goal of human genetics is to identify DNA sequence variants that influence biomedical traits, particularly those related to the onset and progression of human disease. Over the past 25 years, progress in realizing this objective has been transformed by advances in technology, foundational genomic resources and analytical tools, and by access to vast amounts of genotype and phenotype data. Genetic discoveries have substantially improved our understanding of the mechanisms responsible for many rare and common diseases and driven development of novel preventative and therapeutic strategies. Medical innovation will increasingly focus on delivering care tailored to individual patterns of genetic predisposition.

Characterizing genetic profiles for high triglyceride levels in U.S. patients of African ancestry.

Hypertriglyceridemia (HTG) is a common cardiovascular risk factor characterized by elevated triglyceride (TG) levels. Researchers have assessed the genetic factors that influence HTG in studies focused predominantly on individuals of European ancestry. However, relatively little is known about the contribution of genetic variation of HTG in people of African ancestry (AA), potentially constraining research and treatment opportunities. Our objective was to characterize genetic profiles among individuals of AA with mild-to-moderate HTG and severe HTG versus those with normal TGs by leveraging whole-genome sequencing data and longitudinal electronic health records available in the All of Us program. We compared the enrichment of functional variants within five canonical TG metabolism genes, an AA-specific polygenic risk score for TGs, and frequencies of 145 known potentially causal TG variants between HTG patients and normal TG among a cohort of AA patients (N = 15,373). Those with mild-to-moderate HTG (N = 342) and severe HTG (N ≤ 20) were more likely to carry APOA5 p.S19W (odds ratio = 1.94, 95% confidence interval = [1.48-2.54], P = 1.63 × 10-6 and OR = 3.65, 95% confidence interval: [1.22-10.93], P = 0.02, respectively) than those with normal TG. They were also more likely to have an elevated (top 10%) polygenic risk score, elevated carriage of potentially causal variant alleles, and carry any genetic risk factor. Alternative definitions of HTG yielded comparable results. In conclusion, individuals of AA with HTG were enriched for genetic risk factors compared to individuals with normal TGs.

Host genetic effects in pneumonia.

Given the coronavirus disease 2019 (COVID-19) pandemic, investigations into host susceptibility to infectious diseases and downstream sequelae have never been more relevant. Pneumonia is a lung disease that can cause respiratory failure and hypoxia and is a common complication of infectious diseases, including COVID-19. Few genome-wide association studies (GWASs) of host susceptibility and severity of pneumonia have been conducted. We performed GWASs of pneumonia susceptibility and severity in the Vanderbilt University biobank (BioVU) with linked electronic health records (EHRs), including Illumina Expanded Multi-Ethnic Global Array (MEGAEX)-genotyped European ancestry (EA, n= 69,819) and African ancestry (AA, n = 15,603) individuals. Two regions of large effect were identified: the CFTR locus in EA (rs113827944; OR = 1.84, p value = 1.2 × 10-36) and HBB in AA (rs334 [p.Glu7Val]; OR = 1.63, p value = 3.5 × 10-13). Mutations in these genes cause cystic fibrosis (CF) and sickle cell disease (SCD), respectively. After removing individuals diagnosed with CF and SCD, we assessed heterozygosity effects at our lead variants. Further GWASs after removing individuals with CF uncovered an additional association in R3HCC1L (rs10786398; OR = 1.22, p value = 3.5 × 10-8), which was replicated in two independent datasets: UK Biobank (n = 459,741) and 7,985 non-overlapping BioVU subjects, who are genotyped on arrays other than MEGAEX. This variant was also validated in GWASs of COVID-19 hospitalization and lung function. Our results highlight the importance of the host genome in infectious disease susceptibility and severity and offer crucial insight into genetic effects that could potentially influence severity of COVID-19 sequelae.

Interactive network-based clustering and investigation of multimorbidity association matrices with associationSubgraphs.

MOTIVATION: Making sense of networked multivariate association patterns is vitally important to many areas of high-dimensional analysis. Unfortunately, as the data-space dimensions grow, the number of association pairs increases in O(n2); this means that traditional visualizations such as heatmaps quickly become too complicated to parse effectively. RESULTS: Here, we present associationSubgraphs: a new interactive visualization method to quickly and intuitively explore high-dimensional association datasets using network percolation and clustering. The goal is to provide an efficient investigation of association subgraphs, each containing a subset of variables with stronger and more frequent associations among themselves than the remaining variables outside the subset, by showing the entire clustering dynamics and providing subgraphs under all possible cutoff values at once. Particularly, we apply associationSubgraphs to a phenome-wide multimorbidity association matrix generated from an electronic health record and provide an online, interactive demonstration for exploring multimorbidity subgraphs. AVAILABILITY AND IMPLEMENTATION: An R package implementing both the algorithm and visualization components of associationSubgraphs is available at https://github.com/tbilab/associationsubgraphs. Online documentation is available at https://prod.tbilab.org/associationsubgraphs_info/. A demo using a multimorbidity association matrix is available at https://prod.tbilab.org/associationsubgraphs-example/.

High-throughput framework for genetic analyses of adverse drug reactions using electronic health records.

Understanding the contribution of genetic variation to drug response can improve the delivery of precision medicine. However, genome-wide association studies (GWAS) for drug response are uncommon and are often hindered by small sample sizes. We present a high-throughput framework to efficiently identify eligible patients for genetic studies of adverse drug reactions (ADRs) using "drug allergy" labels from electronic health records (EHRs). As a proof-of-concept, we conducted GWAS for ADRs to 14 common drug/drug groups with 81,739 individuals from Vanderbilt University Medical Center's BioVU DNA Biobank. We identified 7 genetic loci associated with ADRs at P < 5 × 10-8, including known genetic associations such as CYP2D6 and OPRM1 for CYP2D6-metabolized opioid ADR. Additional expression quantitative trait loci and phenome-wide association analyses added evidence to the observed associations. Our high-throughput framework is both scalable and portable, enabling impactful pharmacogenomic research to improve precision medicine.

Analysis of genetically determined gene expression suggests role of inflammatory processes in exfoliation syndrome.

BACKGROUND: Exfoliation syndrome (XFS) is an age-related systemic disorder characterized by excessive production and progressive accumulation of abnormal extracellular material, with pathognomonic ocular manifestations. It is the most common cause of secondary glaucoma, resulting in widespread global blindness. The largest global meta-analysis of XFS in 123,457 multi-ethnic individuals from 24 countries identified seven loci with the strongest association signal in chr15q22-25 region near LOXL1. Expression analysis have so far correlated coding and a few non-coding variants in the region with LOXL1 expression levels, but functional effects of these variants is unclear. We hypothesize that analysis of the contribution of the genetically determined component of gene expression to XFS risk can provide a powerful method to elucidate potential roles of additional genes and clarify biology that underlie XFS. RESULTS: Transcriptomic Wide Association Studies (TWAS) using PrediXcan models trained in 48 GTEx tissues leveraging on results from the multi-ethnic and European ancestry GWAS were performed. To eliminate the possibility of false-positive results due to Linkage Disequilibrium (LD) contamination, we i) performed PrediXcan analysis in reduced models removing variants in LD with LOXL1 missense variants associated with XFS, and variants in LOXL1 models in both multiethnic and European ancestry individuals, ii) conducted conditional analysis of the significant signals in European ancestry individuals, and iii) filtered signals based on correlated gene expression, LD and shared eQTLs, iv) conducted expression validation analysis in human iris tissues. We observed twenty-eight genes in chr15q22-25 region that showed statistically significant associations, which were whittled down to ten genes after statistical validations. In experimental analysis, mRNA transcript levels for ARID3B, CD276, LOXL1, NEO1, SCAMP2, and UBL7 were significantly decreased in iris tissues from XFS patients compared to control samples. TWAS genes for XFS were significantly enriched for genes associated with inflammatory conditions. We also observed a higher incidence of XFS comorbidity with inflammatory and connective tissue diseases. CONCLUSION: Our results implicate a role for connective tissues and inflammation pathways in the etiology of XFS. Targeting the inflammatory pathway may be a potential therapeutic option to reduce progression in XFS.

A transcriptome-wide association study identifies novel blood-based gene biomarker candidates for Alzheimer's disease risk.

Alzheimer's disease (ad) adversely affects the health, quality of life and independence of patients. There is a critical need to identify novel blood gene biomarkers for ad risk assessment. We performed a transcriptome-wide association study to identify biomarker candidates for ad risk. We leveraged two sets of gene expression prediction models of blood developed using different reference panels and modeling strategies. By applying the prediction models to a meta-GWAS including 71 880 (proxy) cases and 383 378 (proxy) controls, we identified significant associations of genetically determined expression of 108 genes in blood with ad risk. Of these, 15 genes were differentially expressed between ad patients and controls with concordant directions in measured expression data. With evidence from the analyses based on both genetic instruments and directly measured expression levels, this study identifies 15 genes with strong support as biomarkers in blood for ad risk, which may enhance ad risk assessment and mechanism-focused studies.

Fate or coincidence: do COPD and major depression share genetic risk factors?

Major depressive disorder (MDD) is a common comorbidity in chronic obstructive pulmonary disease (COPD), affecting up to 57% of patients with COPD. Although the comorbidity of COPD and MDD is well established, the causal relationship between these two diseases is unclear. A large-scale electronic health record clinical biobank and genome-wide association study summary statistics for MDD and lung function traits were used to investigate potential shared underlying genetic susceptibility between COPD and MDD. Linkage disequilibrium score regression was used to estimate genetic correlation between phenotypes. Polygenic risk scores (PRS) for MDD and lung function traits were developed and used to perform a phenome-wide association study (PheWAS). Multi-trait-based conditional and joint analysis identified single-nucleotide polymorphisms (SNPs) influencing both lung function and MDD. We found genetic correlations between MDD and all lung function traits were small and not statistically significant. A PRS-MDD was significantly associated with an increased risk of COPD in a PheWAS [odds ratio (OR) = 1.12, 95% confidence interval (CI): 1.09-1.16] when adjusting for age, sex and genetic ancestry, but this relationship became attenuated when controlling for smoking history (OR = 1.08, 95% CI: 1.04-1.13). No significant associations were found between the lung function PRS and MDD. Multi-trait-based conditional and joint analysis identified three SNPs that may contribute to both traits, two of which were previously associated with mood disorders and COPD. Our findings suggest that the observed relationship between COPD and MDD may not be driven by a strong shared genetic architecture.

Cross-talks between gut microbiota and tobacco smoking: a two-sample Mendelian randomization study.

BACKGROUND: Considerable evidence has been reported that tobacco use could cause alterations in gut microbiota composition. The microbiota-gut-brain axis also in turn hinted at a possible contribution of the gut microbiota to smoking. However, population-level studies with a higher evidence level for causality are lacking. METHODS: This study utilized the summary-level data of respective genome-wide association study (GWAS) for 211 gut microbial taxa and five smoking phenotypes to reveal the causal association between the gut microbiota and tobacco smoking. Two-sample bidirectional Mendelian randomization (MR) design was deployed and comprehensively sensitive analyses were followed to validate the robustness of results. We further performed multivariable MR to evaluate the effect of neurotransmitter-associated metabolites on observed associations. RESULTS: Our univariable MR results confirmed the effects of smoking on three taxa (Intestinimonas, Catenibacterium, and Ruminococcaceae, observed from previous studies) with boosted evidence level and identified another 13 taxa which may be causally affected by tobacco smoking. As for the other direction, we revealed that smoking behaviors could be potential consequence of specific taxa abundance. Combining with existing observational evidence, we provided novel insights regarding a positive feedback loop of smoking through Actinobacteria and indicated a potential mechanism for the link between parental smoking and early smoking initiation of their children driven by Bifidobacterium. The multivariable MR results suggested that neurotransmitter-associated metabolites (tryptophan and tyrosine, also supported by previous studies) probably played a role in the action pathway from the gut microbiota to smoking, especially for Actinobacteria and Peptococcus. CONCLUSIONS: In summary, the current study suggested the role of the specific gut microbes on the risk for cigarette smoking (likely involving alterations in metabolites) and in turn smoking on specific gut microbes. Our findings highlighted the hazards of tobacco use for gut flora dysbiosis and shed light on the potential role of specific gut microbiota for smoking behaviors.

APOL1 and the risk of adverse renal outcomes in patients of African ancestry with systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) disproportionately affects individuals of African ancestry (AA) compared to European ancestry (EA). In the general population, high risk (HR) variants in the apolipoprotein L1 (APOL1) gene increase the risk of renal and hypertensive disorders in individuals of AA. Since SLE is characterized by an interferon signature and APOL1 expression is driven by interferon, we examined the hypothesis that APOL1 HR genotypes predominantly drive higher rates of renal and hypertensive-related comorbidities observed in SLE patients of AA versus those of EA. METHODS: We performed a retrospective cohort study in patients with SLE of EA and AA using a genetic biobank linked to de-identified electronic health records. APOL1 HR genotypes were defined as G1/G1, G2/G2, or G1/G2 and low risk (LR) genotypes as 1 or 0 copies of the G1 and G2 alleles. To identify renal and hypertensive-related disorders that differed in prevalence by ancestry, we used a phenome-wide association approach. We then used logistic regression to compare the prevalence of renal and hypertensive-related disorders in EA and AA patients, both including and excluding patients with the APOL1 HR genotype. In a sensitivity analysis, we examined the association of end stage renal disease secondary to lupus nephritis (LN-related ESRD) with ancestry and the APOL1 genotype. RESULTS: We studied 784 patients with SLE; 195 (24.9%) were of AA, of whom 27 (13.8%) had APOL1 HR genotypes. Eighteen renal and hypertensive-related phenotypes were more common in AA than EA patients (p-value ≤ 1.4E-4). All phenotypes remained significantly different after exclusion of patients with APOL1 HR genotypes, and most point odds ratios (ORs) decreased only slightly. Even among ORs with the greatest decrease, risk for AA patients without the APOL1 HR genotype remained significantly elevated compared to EA patients. In the sensitivity analysis, LN-related ESRD was more prevalent in SLE patients of AA versus EA and AA patients with the APOL1 HR genotype versus LR (p-value < .05 for both). CONCLUSION: The higher prevalence of renal and hypertensive disorders in SLE patients of AA compared to those of EA is not fully explained by the presence of APOL1 high risk variants.

Sex differences in the genetic architecture of cognitive resilience to Alzheimer's disease.

Approximately 30% of elderly adults are cognitively unimpaired at time of death despite the presence of Alzheimer's disease neuropathology at autopsy. Studying individuals who are resilient to the cognitive consequences of Alzheimer's disease neuropathology may uncover novel therapeutic targets to treat Alzheimer's disease. It is well established that there are sex differences in response to Alzheimer's disease pathology, and growing evidence suggests that genetic factors may contribute to these differences. Taken together, we sought to elucidate sex-specific genetic drivers of resilience. We extended our recent large scale genomic analysis of resilience in which we harmonized cognitive data across four cohorts of cognitive ageing, in vivo amyloid PET across two cohorts, and autopsy measures of amyloid neuritic plaque burden across two cohorts. These data were leveraged to build robust, continuous resilience phenotypes. With these phenotypes, we performed sex-stratified [n (males) = 2093, n (females) = 2931] and sex-interaction [n (both sexes) = 5024] genome-wide association studies (GWAS), gene and pathway-based tests, and genetic correlation analyses to clarify the variants, genes and molecular pathways that relate to resilience in a sex-specific manner. Estimated among cognitively normal individuals of both sexes, resilience was 20-25% heritable, and when estimated in either sex among cognitively normal individuals, resilience was 15-44% heritable. In our GWAS, we identified a female-specific locus on chromosome 10 [rs827389, ß (females) = 0.08, P (females) = 5.76 × 10-09, ß (males) = -0.01, P(males) = 0.70, ß (interaction) = 0.09, P (interaction) = 1.01 × 10-04] in which the minor allele was associated with higher resilience scores among females. This locus is located within chromatin loops that interact with promoters of genes involved in RNA processing, including GATA3. Finally, our genetic correlation analyses revealed shared genetic architecture between resilience phenotypes and other complex traits, including a female-specific association with frontotemporal dementia and male-specific associations with heart rate variability traits. We also observed opposing associations between sexes for multiple sclerosis, such that more resilient females had a lower genetic susceptibility to multiple sclerosis, and more resilient males had a higher genetic susceptibility to multiple sclerosis. Overall, we identified sex differences in the genetic architecture of resilience, identified a female-specific resilience locus and highlighted numerous sex-specific molecular pathways that may underly resilience to Alzheimer's disease pathology. This study illustrates the need to conduct sex-aware genomic analyses to identify novel targets that are unidentified in sex-agnostic models. Our findings support the theory that the most successful treatment for an individual with Alzheimer's disease may be personalized based on their biological sex and genetic context.

Race, Genotype, and Azathioprine Discontinuation : A Cohort Study.

BACKGROUND: Thiopurines are an important class of immunosuppressants despite their risk for hematopoietic toxicity and narrow therapeutic indices. Benign neutropenia related to an ACKR1 variant (rs2814778-CC) is common among persons of African ancestries. OBJECTIVE: To test whether rs2814778-CC was associated with azathioprine discontinuation attributed to hematopoietic toxicity and lower thiopurine dosing. DESIGN: Retrospective cohort study. SETTING: Two tertiary care centers. PATIENTS: Thiopurine users with White or Black race. MEASUREMENTS: Azathioprine discontinuation attributed to hematopoietic toxicity. Secondary outcomes included weight-adjusted final dose, leukocyte count, and change in leukocyte count. RESULTS: The rate of azathioprine discontinuation attributed to hematopoietic toxicity was 3.92 per 100 person-years among patients with the CC genotype (n = 101) and 1.34 per 100 person-years among those with the TT or TC genotype (n = 1365) (hazard ratio [HR] from competing-risk model, 2.92 [95% CI, 1.57 to 5.41]). The risk remained significant after adjustment for race (HR, 2.61 [CI, 1.01 to 6.71]). The risk associated with race alone (HR, 2.13 [CI, 1.21 to 3.75]) was abrogated by adjustment for genotype (HR, 1.13 [CI, 0.48 to 2.69]). Lower last leukocyte count and lower dosing were significant among patients with the CC genotype. Lower dosing was validated in an external cohort of 94 children of African ancestries prescribed the thiopurine 6-mercaptopurine (6-MP) for acute lymphoblastic leukemia. The CC genotype was independently associated with lower 6-MP dose intensity relative to the target daily dose of 75 mg/m2 (median, 0.83 [IQR, 0.70 to 0.94] for the CC genotype vs. 0.94 [IQR, 0.72 to 1.13] for the TT or TC genotype; P = 0.013). LIMITATIONS: Unmeasured confounding; data limited to tertiary centers. CONCLUSION: Patients with the CC genotype had higher risk for azathioprine discontinuation attributed to hematopoietic toxicity and lower thiopurine doses. Genotype was associated with those risks, even after adjustment for race. PRIMARY FUNDING SOURCE: National Institutes of Health.

Recommendations for Statistical Reporting in Cardiovascular Medicine: A Special Report From the American Heart Association.

Statistical analyses are a crucial component of the biomedical research process and are necessary to draw inferences from biomedical research data. The application of sound statistical methodology is a prerequisite for publication in the American Heart Association (AHA) journal portfolio. The objective of this document is to summarize key aspects of statistical reporting that might be most relevant to the authors, reviewers, and readership of AHA journals. The AHA Scientific Publication Committee convened a task force to inventory existing statistical standards for publication in biomedical journals and to identify approaches suitable for the AHA journal portfolio. The experts on the task force were selected by the AHA Scientific Publication Committee, who identified 12 key topics that serve as the section headers for this document. For each topic, the members of the writing group identified relevant references and evaluated them as a resource to make the standards summarized herein. Each section was independently reviewed by an expert reviewer who was not part of the task force. Expert reviewers were also permitted to comment on other sections if they chose. Differences of opinion were adjudicated by consensus. The standards presented in this report are intended to serve as a guide for high-quality reporting of statistical analyses methods and results.

Genome-wide association analyses of common infections in a large practice-based biobank.

INTRODUCTION: Infectious diseases are common causes of morbidity and mortality worldwide. Susceptibility to infection is highly heritable; however, little has been done to identify the genetic determinants underlying common infectious diseases. One GWAS was performed using 23andMe information about self-reported infections; we set out to confirm previous loci and identify new ones using medically diagnosed infections. METHODS: We used the electronic health record (EHR)-based biobank at Vanderbilt and diagnosis codes to identify cases of 12 infectious diseases in white patients: urinary tract infection, pneumonia, chronic sinus infections, otitis media, candidiasis, streptococcal pharyngitis, herpes zoster, herpes labialis, hepatitis B, infectious mononucleosis, tuberculosis (TB) or a positive TB test, and hepatitis C. We selected controls from patients with no diagnosis code for the candidate disease and matched by year of birth, sex, and calendar year at first and last EHR visits. We conducted GWAS using SAIGE and transcriptome-wide analysis (TWAS) using S-PrediXcan. We also conducted phenome-wide association study to understand associations between identified genetic variants and clinical phenotypes. RESULTS: We replicated three 23andMe loci (p ≤ 0.05): herpes zoster and rs7047299-A (p = 2.6 × 10-3) and rs2808290-C (p = 9.6 × 10-3;); otitis media and rs114947103-C (p = 0.04). We also identified 2 novel regions (p ≤ 5 × 10-8): rs113235453-G for otitis media (p = 3.04 × 10-8), and rs10422015-T for candidiasis (p = 3.11 × 10-8). In TWAS, four gene-disease associations were significant: SLC30A9 for otitis media (p = 8.06 × 10-7); LRP3 and WDR88 for candidiasis (p = 3.91 × 10-7 and p = 1.95 × 10-6); and AAMDC for hepatitis B (p = 1.51 × 10-6). CONCLUSION: We conducted GWAS and TWAS for 12 infectious diseases and identified novel genetic contributors to the susceptibility of infectious diseases.

A transcriptome-wide association study identifies novel candidate susceptibility genes for prostate cancer risk.

A large proportion of heritability for prostate cancer risk remains unknown. Transcriptome-wide association study combined with validation comparing overall levels will help to identify candidate genes potentially playing a role in prostate cancer development. Using data from the Genotype-Tissue Expression Project, we built genetic models to predict normal prostate tissue gene expression using the statistical framework PrediXcan, a modified version of the unified test for molecular signatures and Joint-Tissue Imputation. We applied these prediction models to the genetic data of 79 194 prostate cancer cases and 61 112 controls to investigate the associations of genetically determined gene expression with prostate cancer risk. Focusing on associated genes, we compared their expression in prostate tumor vs normal prostate tissue, compared methylation of CpG sites located at these loci in prostate tumor vs normal tissue, and assessed the correlations between the differentiated genes' expression and the methylation of corresponding CpG sites, by analyzing The Cancer Genome Atlas (TCGA) data. We identified 573 genes showing an association with prostate cancer risk at a false discovery rate (FDR) ≤ 0.05, including 451 novel genes and 122 previously reported genes. Of the 573 genes, 152 showed differential expression in prostate tumor vs normal tissue samples. At loci of 57 genes, 151 CpG sites showed differential methylation in prostate tumor vs normal tissue samples. Of these, 20 CpG sites were correlated with expression of 11 corresponding genes. In this TWAS, we identified novel candidate susceptibility genes for prostate cancer risk, providing new insights into prostate cancer genetics and biology.

Integration of DNA sequencing with population pharmacokinetics to improve the prediction of irinotecan exposure in cancer patients.

BACKGROUND: Irinotecan (CPT-11) is an anticancer agent widely used to treat adult solid tumours. Large interindividual variability in the clearance of irinotecan and SN-38, its active and toxic metabolite, results in highly unpredictable toxicity. METHODS: In 217 cancer patients treated with intravenous irinotecan single agent or in combination, germline DNA was used to interrogate the variation in 84 genes by next-generation sequencing. A stepwise analytical framework including a population pharmacokinetic model with SNP- and gene-based testing was used to identify demographic/clinical/genetic factors that influence the clearance of irinotecan and SN-38. RESULTS: Irinotecan clearance was influenced by rs4149057 in SLCO1B1, body surface area, and co-administration of 5-fluorouracil/leucovorin/bevacizumab. SN-38 clearance was influenced by rs887829 in UGT1A1, pre-treatment total bilirubin, and EGFR rare variant burden. Within each UGT1A1 genotype group, elevated pre-treatment total bilirubin and/or presence of at least one rare variant in EGFR resulted in significantly lower SN-38 clearance. The model reduced the interindividual variability in irinotecan clearance from 38 to 34% and SN-38 clearance from 49 to 32%. CONCLUSIONS: This new model significantly reduced the interindividual variability in the clearance of irinotecan and SN-38. New genetic factors of variability in clearance have been identified.

Two polymorphic gene loci associated with treprostinil dose in pulmonary arterial hypertension.

OBJECTIVE: Prostacyclin infusion for pulmonary arterial hypertension (PAH) is an effective therapy with varied dosing requirements and clinical response. The major aim of this study was to determine new biologically-based predictors of prostacyclin treatment response heterogeneity. METHODS: Ninety-eight patients with hemodynamically defined PAH at two academic medical centers volunteered for registry studies. A stable dose of treprostinil was the quantitative phenotype for the genome-wide association study (GWAS). Candidate genes with the largest effect sizes and strongest statistical associations were further characterized with in silico and in-vitro assays to confirm mechanistic hypotheses. The clinical significance of these candidate predictors was assessed for mechanistically consistent physiologic effects in an independent cohort of patients. RESULTS: GWAS identified three loci for association with P < 10-6. All three loci had clinically significant effect sizes. Specific single-nucleotide polymorphisms (SNPs) at two of the loci: rs11078738 in phosphoribosylformylglycinamidine synthase and rs10023113 in CAMK2D encoded sequence changes with clear predicted consequences. Production of the primary mediator of prostacyclin-induced vasodilation, cyclic AMP, was reduced in human cell lines by the missense variant rs11078738 (p.L621P). Located in the promoter of CAMK2D, the allele of rs10023113 associated with a higher treprostinil dose has higher ventricular transcription of CAMK2δ. At initial diagnostic catheterization in a separate cohort of patients, the same allele of rs10023113 was associated with elevated right mean atrial and ventricular diastolic pressures. CONCLUSIONS: The quantitative phenotype of stable treprostinil dose identified two gene loci associated with pharmacodynamic response and right ventricular function in PAH worth further investigation.

GRIK5 Genetically Regulated Expression Associated with Eye and Vascular Phenomes: Discovery through Iteration among Biobanks, Electronic Health Records, and Zebrafish.

Although the use of model systems for studying the mechanism of mutations that have a large effect is common, we highlight here the ways that zebrafish-model-system studies of a gene, GRIK5, that contributes to the polygenic liability to develop eye diseases have helped to illuminate a mechanism that implicates vascular biology in eye disease. A gene-expression prediction derived from a reference transcriptome panel applied to BioVU, a large electronic health record (EHR)-linked biobank at Vanderbilt University Medical Center, implicated reduced GRIK5 expression in diverse eye diseases. We tested the function of GRIK5 by depletion of its ortholog in zebrafish, and we observed reduced blood vessel numbers and integrity in the eye and increased vascular permeability. Analyses of EHRs in >2.6 million Vanderbilt subjects revealed significant comorbidity of eye and vascular diseases (relative risks 2-15); this comorbidity was confirmed in 150 million individuals from a large insurance claims dataset. Subsequent studies in >60,000 genotyped BioVU participants confirmed the association of reduced genetically predicted expression of GRIK5 with comorbid vascular and eye diseases. Our studies pioneer an approach that allows a rapid iteration of the discovery of gene-phenotype relationships to the primary genetic mechanism contributing to the pathophysiology of human disease. Our findings also add dimension to the understanding of the biology driven by glutamate receptors such as GRIK5 (also referred to as GLUK5 in protein form) and to mechanisms contributing to human eye diseases.