Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

Do research participants share genomic screening results with family members?

The public health impact of genomic screening can be enhanced by cascade testing. However, cascade testing depends on communication of results to family members. While the barriers and facilitators of family communication have been researched following clinical genetic testing, the factors impacting the dissemination of genomic screening results are unknown. Using the pragmatic Electronic Medical Records and Genomics Network-3 (eMERGE-3) study, we explored the reported sharing practices of participants who underwent genomic screening across the United States. Six eMERGE-3 sites returned genomic screening results for mostly dominant medically actionable disorders and surveyed adult participants regarding communication of results with first-degree relatives. Across the sites, 279 participants completed a 1-month and/or 6-month post-results survey. By 6 months, only 34% of the 156 respondents shared their results with all first-degree relatives and 4% did not share with any. Over a third (39%) first-degree relatives were not notified of the results. Half (53%) of participants who received their results from a genetics provider shared them with all first-degree relatives compared with 11% of participants who received their results from a non-genetics provider. The most frequent reasons for sharing were a feeling of obligation (72%) and that the information could help family members make medical decisions (72%). The most common reasons indicated for not sharing were that the family members were too young (38%), or they were not in contact (25%) or not close to them (25%). These data indicate that the professional returning the results may impact sharing patterns, suggesting that there is a need to continue to educate healthcare providers regarding approaches to facilitate sharing of genetic results within families. Finally, these data suggest that interventions to increase sharing may be universally effective regardless of the origin of the genetic result.

Generalizability of Polygenic Risk Scores for Breast Cancer Among Women With European, African, and Latinx Ancestry

IMPORTANCE

Multiple polygenic risk scores (PRSs) for breast cancer have been developed from large research consortia; however, their generalizability to diverse clinical settings is unknown.

OBJECTIVE

To examine the performance of previously developed breast cancer PRSs in a clinical setting for women of European, African, and Latinx ancestry.

DESIGN, SETTING, AND PARTICIPANTS

This cohort study using the Electronic Medical Records and Genomics (eMERGE) network data set included 39 591 women from 9 contributing medical centers in the US that had electronic medical records (EMR) linked to genotype data. Breast cancer cases and controls were identified through a validated EMR phenotyping algorithm.

MAIN OUTCOMES AND MEASURES

Multivariable logistic regression was used to assess the association between breast cancer risk and 7 previously developed PRSs, adjusting for age, study site, breast cancer family history, and first 3 ancestry informative principal components.

RESULTS

This study included 39 591 women: 33 594 with European, 3801 with African, and 2196 with Latinx ancestry. The mean (SD) age at breast cancer diagnosis was 60.7 (13.0), 58.8 (12.5), and 60.1 (13.0) years for women with European, African, and Latinx ancestry, respectively. PRSs derived from women with European ancestry were associated with breast cancer risk in women with European ancestry (highest odds ratio [OR] per 1-SD increase, 1.46; 95% CI, 1.41-1.51), women with Latinx ancestry (highest OR, 1.31; 95% CI, 1.09-1.58), and women with African ancestry (OR, 1.19; 95% CI, 1.05-1.35). For women with European ancestry, this association with breast cancer risk was largest in the extremes of the PRS distribution, with ORs ranging from 2.19 (95% CI, 1.84-2.53) to 2.48 (95% CI, 1.89-3.25) for the 3 different PRSs examined for those in the highest 1% of the PRS compared with those in the middle quantile. Among women with Latinx and African ancestries at the extremes of the PRS distribution, there were no statistically significant associations.

CONCLUSIONS AND RELEVANCE

This cohort study found that PRS models derived from women with European ancestry for breast cancer risk generalized well for women with European, Latinx, and African ancestries across different clinical settings, although the effect sizes for women with African ancestry were smaller, likely because of differences in risk allele frequencies and linkage disequilibrium patterns. These results highlight the need to improve representation of diverse population groups, particularly women with African ancestry, in genomic research cohorts.

Genomic considerations for FHIR®; eMERGE implementation lessons

Structured representation of clinical genetic results is necessary for advancing precision medicine. The Electronic Medical Records and Genomics (eMERGE) Network's Phase III program initially used a commercially developed XML message format for standardized and structured representation of genetic results for electronic health record (EHR) integration. In a desire to move towards a standard representation, the network created a new standardized format based upon Health Level Seven Fast Healthcare Interoperability Resources (HL7® FHIR®), to represent clinical genomics results. These new standards improve the utility of HL7® FHIR® as an international healthcare interoperability standard for management of genetic data from patients. This work advances the establishment of standards that are being designed for broad adoption in the current health information technology landscape.

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.

An electronic health record (EHR) log analysis shows limited clinician engagement with unsolicited genetic test results

How clinicians utilize medically actionable genomic information, displayed in the electronic health record (EHR), in medical decision-making remains unknown. Participating sites of the Electronic Medical Records and Genomics (eMERGE) Network have invested resources into EHR integration efforts to enable the display of genetic testing data across heterogeneous EHR systems. To assess clinicians’ engagement with unsolicited EHR-integrated genetic test results of eMERGE participants within a large tertiary care academic medical center, we analyzed automatically generated EHR access log data. We found that clinicians viewed only 1% of all the eMERGE genetic test results integrated in the EHR. Using a cluster analysis, we also identified different user traits associated with varying degrees of engagement with the EHR-integrated genomic data. These data contribute important empirical knowledge about clinicians limited and brief engagements with unsolicited EHR-integrated genetic test results of eMERGE participants. Appreciation for user-specific roles provide additional context for why certain users were more or less engaged with the unsolicited results. This study highlights opportunities to use EHR log data as a performance metric to more precisely inform ongoing EHR-integration efforts and decisions about the allocation of informatics resources in genomic research.

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.

Association of Thyroid Function Genetic Predictors With Atrial Fibrillation: A Phenome-Wide Association Study and Inverse-Variance Weighted Average Meta-analysis

Importance

Thyroid hormone levels are tightly regulated through feedback inhibition by thyrotropin, produced by the pituitary gland. Hyperthyroidism is overwhelmingly due to thyroid disorders and is well recognized to contribute to a wide spectrum of cardiovascular morbidity, particularly the increasingly common arrhythmia atrial fibrillation (AF).

Objective

To determine the association between genetically determined thyrotropin levels and AF.

Design, Setting, and Participants

This phenome-wide association study scanned 1318 phenotypes associated with a polygenic predictor of thyrotropin levels identified by a previously published genome-wide association study that included participants of European ancestry. North American individuals of European ancestry with longitudinal electronic health records were analyzed from May 2008 to November 2016. Analysis began March 2018.

Main Outcomes and Measures

Clinical diagnoses associated with a polygenic predictor of thyrotropin levels.

Exposures

Genetically determined thyrotropin levels.

Results

Of 37 154 individuals, 19 330 (52%) were men. The thyrotropin polygenic predictor was positively associated with hypothyroidism (odds ratio [OR], 1.10; 95% CI, 1.07-1.14; P = 5 × 10-11) and inversely associated with diagnoses related to hyperthyroidism (OR, 0.64; 95% CI, 0.54-0.74; P = 2 × 10-8 for toxic multinodular goiter). Among nonthyroid associations, the top association was AF/flutter (OR, 0.93; 95% CI, 0.9-0.95; P = 9 × 10-7). When the analyses were repeated excluding 9801 individuals with any diagnoses of a thyroid-related disease, the AF association persisted (OR, 0.91; 95% CI, 0.88-0.95; P = 2.9 × 10-6). To replicate this association, we conducted an inverse-variance weighted average meta-analysis using AF single-nucleotide variant weights from a genome-wide association study of 17 931 AF cases and 115 142 controls. As in the discovery analyses, each SD increase in predicted thyrotropin was associated with a decreased risk of AF (OR, 0.86; 95% CI, 0.79-0.93; P = 4.7 × 10-4). In a set of AF cases (n = 745) and controls (n = 1680) older than 55 years, directly measured thyrotropin levels that fell within the normal range were inversely associated with AF risk (OR, 0.91; 95% CI, 0.83-0.99; P = .04).

Conclusions and Relevance

This study suggests a role for genetically determined variation in thyroid function within a physiologically accepted normal range as a risk factor for AF. The clinical decision to treat subclinical thyroid disease should incorporate the risk for AF as antithyroid medications to treat hyperthyroidism may reduce AF risk and thyroid hormone replacement for hypothyroidism may increase AF risk.

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.

Correction: Evaluation of the cost and effectiveness of diverse recruitment methods for a genetic screening study

The original version of this Article contained an error in the undergraduate degree awarded to the author Ian Halim, which was incorrectly given as BS. This has now been corrected to BA in both the PDF and HTML versions of the Article.

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.

Torsade de pointes and systemic azole antifungal agents: Analysis of global spontaneous safety reports

Background: Literature about torsade de pointes induced by azole antifungal agents is scarce, despite the well-known association. Furthermore, little is known about the latency time between commencing an azole antifungal agent and developing torsade de pointes. The objectives of the present study were therefore to identify all cases of torsade de pointes associated with systemic azole antifungal use reported to the WHO monitoring centre (Uppsala, Sweden) and to determine the latency times between commencing the azole and developing torsade de pointes. Methods: Investigator-driven, retrospective, descriptive analysis of post-marketing pharmacovigilance data regarding systemic azole antifungal agents and the development of torsade de pointes reported to the WHO monitoring centre 1995-2015. Results: 191 cases were reported as follows: fluconazole 130, itraconazole 22, ketoconazole 5, posaconazole 1, voriconazole 33. More than half of all cases involved concomitant suspected or interacting drugs. The median latency times between starting the azole and developing torsade de pointes ranged from 1 (posaconazole) - 9.5 days (itraconazole), range <1-250). Conclusions: Clinicians should be aware of these features of azole-associated torsade de pointes, avoid interacting drugs if at all possible and monitor at-risk patients.

[Drug-drug interactions with HIV treatments]

We develop a differential diagnosis and diagnostic approach to a possible drug-drug interaction between a cART based on boosted atazanavir and newly given drugs. History taking should not only encompass prescribed but also over-the-counter medication. Exchange of information between attending physicians as to therapeutic changes, documented side effects and adherence, therapeutic drug monitoring of cART and evaluation by pharmacologist in complex situations are the diagnostic tools at hand.

Steady-state fully noninductive current driven by electron cyclotron waves in a magnetically confined plasma

A steady-state, fully noninductive plasma current has been sustained for the first time in a tokamak using electron cyclotron current drive only. In this discharge, 123 kA of current have been sustained for the entire gyrotron pulse duration of 2 s. Careful distribution across the plasma minor radius of the power deposited from three 0. 5-MW gyrotrons was essential for reaching steady-state conditions. With central current drive, up to 153 kA of current have been fully replaced transiently for 100 ms. The noninductive scenario is confirmed by the ability to recharge the Ohmic transformer. The dependence of the current drive efficiency on the minor radius is also demonstrated.