With the emergence of large-scale sequencing data, methods for improving power in rare variant association tests are needed. Here we show that adjusting for common variant polygenic scores improves yield in gene-based rare variant association tests across 65 quantitative traits in the UK Biobank (up to 20% increase at α = 2.6 × 10-6), without marked increases in false-positive rates or genomic inflation. Benefits were seen for various models, with the largest improvements seen for efficient sparse mixed-effects models. Our results illustrate how polygenic score adjustment can efficiently improve power in rare variant association discovery.

Myocardial interstitial fibrosis is associated with cardiovascular disease and adverse prognosis. Here, to investigate the biological pathways that underlie fibrosis in the human heart, we developed a machine learning model to measure native myocardial T1 time, a marker of myocardial fibrosis, in 41,505 UK Biobank participants who underwent cardiac magnetic resonance imaging. Greater T1 time was associated with diabetes mellitus, renal disease, aortic stenosis, cardiomyopathy, heart failure, atrial fibrillation, conduction disease and rheumatoid arthritis. Genome-wide association analysis identified 11 independent loci associated with T1 time. The identified loci implicated genes involved in glucose transport (SLC2A12), iron homeostasis (HFE, TMPRSS6), tissue repair (ADAMTSL1, VEGFC), oxidative stress (SOD2), cardiac hypertrophy (MYH7B) and calcium signaling (CAMK2D). Using a transforming growth factor β1-mediated cardiac fibroblast activation assay, we found that 9 of the 11 loci consisted of genes that exhibited temporal changes in expression or open chromatin conformation supporting their biological relevance to myofibroblast cell state acquisition. By harnessing machine learning to perform large-scale quantification of myocardial interstitial fibrosis using cardiac imaging, we validate associations between cardiac fibrosis and disease, and identify new biologically relevant pathways underlying fibrosis.

BACKGROUND: Absence of a dicrotic notch on finger photoplethysmography is an easily ascertainable and inexpensive trait that has been associated with age and prevalent cardiovascular disease. However, the trait exists along a continuum, and little is known about its genetic underpinnings or prognostic value for incident cardiovascular disease. METHODS: In 169 787 participants in the UK Biobank, we identified absent dicrotic notch on photoplethysmography and created a novel continuous trait reflecting notch smoothness using machine learning. Next, we determined the heritability, genetic basis, polygenic risk, and clinical relations for the binary absent notch trait and the newly derived continuous notch smoothness trait. RESULTS: Heritability of the continuous notch smoothness trait was 7.5%, compared with 5.6% for the binary absent notch trait. A genome-wide association study of notch smoothness identified 15 significant loci, implicating genes including NT5C2 (P=1.2×10-26), IGFBP3 (P=4.8×10-18), and PHACTR1 (P=1.4×10-13), compared with 6 loci for the binary absent notch trait. Notch smoothness stratified risk of incident myocardial infarction or coronary artery disease, stroke, heart failure, and aortic stenosis. A polygenic risk score for notch smoothness was associated with incident cardiovascular disease and all-cause death in UK Biobank participants without available photoplethysmography data. CONCLUSIONS: We found that a machine learning derived continuous trait reflecting dicrotic notch smoothness on photoplethysmography was heritable and associated with genes involved in vascular stiffness. Greater notch smoothness was associated with greater risk of incident cardiovascular disease. Raw digital phenotyping may identify individuals at risk for disease via specific genetic pathways.

BACKGROUND: Common genetic variation has been shown to account for a large proportion of ASD heritability. Polygenic scores generated for autism spectrum disorder (ASD-PGS) using the most recent discovery data, however, explain less variance than expected, despite reporting significant associations with ASD and other ASD-related traits. Here, we investigate the extent to which information loss on the target study genome-wide microarray weakens the predictive power of the ASD-PGS. METHODS: We studied genotype data from three cohorts of individuals with high familial liability for ASD: The Early Autism Risk Longitudinal Investigation (EARLI), Markers of Autism Risk in Babies-Learning Early Signs (MARBLES), and the Infant Brain Imaging Study (IBIS), and one population-based sample, Study to Explore Early Development Phase I (SEED I). Individuals were genotyped on different microarrays ranging from 1 to 5 million sites. Coverage of the top 88 genome-wide suggestive variants implicated in the discovery was evaluated in all four studies before quality control (QC), after QC, and after imputation. We then created a novel method to assess coverage on the resulting ASD-PGS by correlating a PGS informed by a comprehensive list of variants to a PGS informed with only the available variants. RESULTS: Prior to imputations, None of the four cohorts directly or indirectly covered all 88 variants among the measured genotype data. After imputation, the two cohorts genotyped on 5-million arrays reached full coverage. Analysis of our novel metric showed generally high genome-wide coverage across all four studies, but a greater number of SNPs informing the ASD-PGS did not result in improved coverage according to our metric. LIMITATIONS: The studies we analyzed contained modest sample sizes. Our analyses included microarrays with more than 1-million sites, so smaller arrays such as Global Diversity and the PsychArray were not included. Our PGS metric for ASD is only generalizable to samples of European ancestries, though the coverage metric can be computed for traits that have sufficiently large-sized discovery findings in other ancestries. CONCLUSIONS: We show that commonly used genotyping microarrays have incomplete coverage for common ASD variants, and imputation cannot always recover lost information. Our novel metric provides an intuitive approach to reporting information loss in PGS and an alternative to reporting the total number of SNPs included in the PGS. While applied only to ASD here, this metric can easily be used with other traits.

Background

Rare sequence variation in genes underlying cardiac repolarization and common polygenic variation influence QT interval duration. However, current clinical genetic testing of individuals with unexplained QT prolongation is restricted to examination of monogenic rare variants. The recent emergence of large-scale biorepositories with sequence data enables examination of the joint contribution of rare and common variations to the QT interval in the population.

Methods

We performed a genome-wide association study of the QTc in 84 630 UK Biobank participants and created a polygenic risk score (PRS). Among 26 976 participants with whole-genome sequencing and ECG data in the TOPMed (Trans-Omics for Precision Medicine) program, we identified 160 carriers of putative pathogenic rare variants in 10 genes known to be associated with the QT interval. We examined QTc associations with the PRS and with rare variants in TOPMed.

Results

Fifty-four independent loci were identified by genome-wide association study in the UK Biobank. Twenty-one loci were novel, of which 12 were replicated in TOPMed. The PRS composed of 1 110 494 common variants was significantly associated with the QTc in TOPMed (ΔQTc_{/decile of PRS}=1.4 ms [95% CI, 1.3 to 1.5]; P=1.1×10^-196). Carriers of putative pathogenic rare variants had longer QTc than noncarriers (ΔQTc=10.9 ms [95% CI, 7.4 to 14.4]). Of individuals with QTc>480 ms, 23.7% carried either a monogenic rare variant or had a PRS in the top decile (3.4% monogenic, 21% top decile of PRS).

Conclusions

QTc duration in the population is influenced by both rare variants in genes underlying cardiac repolarization and polygenic risk, with a sizeable contribution from polygenic risk. Comprehensive assessment of the genetic determinants of QTc prolongation includes incorporation of both polygenic and monogenic risk.

Background

Alterations in electrocardiographic (ECG) intervals are well-known markers for arrhythmia and sudden cardiac death (SCD) risk. While the genetics of arrhythmia syndromes have been studied, relations between electrocardiographic intervals and rare genetic variation at a population level are poorly understood.

Methods

Using a discovery sample of 29 000 individuals with whole-genome sequencing from Trans-Omics in Precision Medicine and replication in nearly 100 000 with whole-exome sequencing from the UK Biobank and MyCode, we examined associations between low-frequency and rare coding variants with 5 routinely measured electrocardiographic traits (RR, P-wave, PR, and QRS intervals and corrected QT interval).

Results

We found that rare variants associated with population-based electrocardiographic intervals identify established monogenic SCD genes (KCNQ1, KCNH2, and SCN5A), a controversial monogenic SCD gene (KCNE1), and novel genes (PAM and MFGE8) involved in cardiac conduction. Loss-of-function and pathogenic SCN5A variants, carried by 0.1% of individuals, were associated with a nearly 6-fold increased odds of the first-degree atrioventricular block (P=8.4×10^-5). Similar variants in KCNQ1 and KCNH2 (0.2% of individuals) were associated with a 23-fold increased odds of marked corrected QT interval prolongation (P=4×10^-25), a marker of SCD risk. Incomplete penetrance of such deleterious variation was common as over 70% of carriers had normal electrocardiographic intervals.

Conclusions

Our findings indicate that large-scale high-depth sequence data and electrocardiographic analysis identifies monogenic arrhythmia susceptibility genes and rare variants with large effects. Known pathogenic variation in conventional arrhythmia and SCD genes exhibited incomplete penetrance and accounted for only a small fraction of marked electrocardiographic interval prolongation.

Accurate and efficient classification of variant pathogenicity is critical for research and clinical care. Using data from three large studies, we demonstrate that population-based associations between rare variants and quantitative endophenotypes for three monogenic diseases (low-density-lipoprotein cholesterol for familial hypercholesterolemia, electrocardiographic QTc interval for long QT syndrome, and glycosylated hemoglobin for maturity-onset diabetes of the young) provide evidence for variant pathogenicity. Effect sizes are associated with pathogenic ClinVar assertions (P < 0.001 for each trait) and discriminate pathogenic from non-pathogenic variants (area under the curve 0.82-0.84 across endophenotypes). An effect size threshold of ≥ 0.5 times the endophenotype standard deviation nominates up to 35% of rare variants of uncertain significance or not in ClinVar in disease susceptibility genes with pathogenic potential. We propose that variant associations with quantitative endophenotypes for monogenic diseases can provide evidence supporting pathogenicity.