Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Polygenic Risk Score Associates With Atherosclerotic Plaque Characteristics at Autopsy.

BACKGROUND: Polygenic risk scores (PRSs) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically. METHODS: From 4327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner for sudden death between 1994 and 2015, 2455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas, and thrombotic CAD. RESULTS: After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age, 48.8±14.7 years; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared with subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% versus 50.4%±38.7%; adjusted P<0.001) and a higher frequency of calcification (69.6% versus 35.8%; adjusted P=0.004) and thin-cap fibroatheroma (26.7% versus 9.5%; adjusted P=0.007). Even after adjustment for traditional CAD risk factors, subjects within the highest PRS quintile had higher odds of severe atherosclerosis (ie, ≥75% stenosis; adjusted odds ratio, 3.77 [95% CI, 2.10-6.78]; P<0.001) and plaque rupture (adjusted odds ratio, 4.05 [95% CI, 2.26-7.24]; P<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged ≤50 years (adjusted odds ratio, 4.08 [95% CI, 2.01-8.30]; P<0.001). No statistically significant associations were observed with plaque erosion after adjusting for covariates. CONCLUSIONS: This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects.

Mitochondrial genetic variation and risk of chronic kidney disease and acute kidney injury in UK Biobank participants.

Experimental models suggest an important role for mitochondrial dysfunction in the pathogenesis of chronic kidney disease (CKD) and acute kidney injury (AKI), but little is known regarding the impact of common mitochondrial genetic variation on kidney health. We sought to evaluate associations of inherited mitochondrial DNA (mtDNA) variation with risk of CKD and AKI in a large population-based cohort. We categorized UK Biobank participants who self-identified as white into eight distinct mtDNA haplotypes, which were previously identified based on their associations with phenotypes associated with mitochondrial DNA copy number, a measure of mitochondrial function. We used linear and logistic regression models to evaluate associations of these mtDNA haplotypes with estimated glomerular filtration rate by serum creatinine and cystatin C (eGFRCr-CysC, N = 362,802), prevalent (N = 416 cases) and incident (N = 405 cases) end-stage kidney disease (ESKD), AKI defined by diagnostic codes (N = 14,170 cases), and urine albumin/creatinine ratio (ACR, N = 114,662). The mean age was 57 ± 8 years and the mean eGFR was 90 ± 14 ml/min/1.73 m2. MtDNA haplotype was significantly associated with eGFR (p = 2.8E-12), but not with prevalent ESKD (p = 5.9E-2), incident ESKD (p = 0.93), AKI (p = 0.26), or urine ACR (p = 0.54). The association of mtDNA haplotype with eGFR remained significant after adjustment for diabetes mellitus and hypertension (p = 1.2E-10). When compared to the reference haplotype, mtDNA haplotypes I (ß = 0.402, standard error (SE) = 0.111; p = 2.7E-4), IV (ß = 0.430, SE = 0.073; p = 4.2E-9), and V (ß = 0.233, SE = 0.050; p = 2.7E-6) were each associated with higher eGFR. Among self-identified white UK Biobank participants, mtDNA haplotype was associated with eGFR, but not with ESKD, AKI or albuminuria.

Blood-derived mitochondrial DNA copy number is associated with gene expression across multiple tissues and is predictive for incident neurodegenerative disease.

Mitochondrial DNA copy number (mtDNA-CN) is a proxy for mitochondrial function and is associated with aging-related diseases. However, it is unclear how mtDNA-CN measured in blood can reflect diseases that primarily manifest in other tissues. Using the Genotype-Tissue Expression Project, we interrogated relationships between mtDNA-CN measured in whole blood and gene expression from whole blood and 47 additional tissues in 419 individuals. mtDNA-CN was significantly associated with expression of 700 genes in whole blood, including nuclear genes required for mtDNA replication. Significant enrichment was observed for splicing and ubiquitin-mediated proteolysis pathways, as well as target genes for the mitochondrial transcription factor NRF1. In nonblood tissues, there were more significantly associated genes than expected in 30 tissues, suggesting that global gene expression in those tissues is correlated with blood-derived mtDNA-CN. Neurodegenerative disease pathways were significantly associated in multiple tissues, and in an independent data set, the UK Biobank, we observed that higher mtDNA-CN was significantly associated with lower rates of both prevalent (OR = 0.89, CI = 0.83; 0.96) and incident neurodegenerative disease (HR = 0.95, 95% CI = 0.91;0.98). The observation that mtDNA-CN measured in blood is associated with gene expression in other tissues suggests that blood-derived mtDNA-CN can reflect metabolic health across multiple tissues. Identification of key pathways including splicing, RNA binding, and catalysis reinforces the importance of mitochondria in maintaining cellular homeostasis. Finally, validation of the role of mtDNA CN in neurodegenerative disease in a large independent cohort study solidifies the link between blood-derived mtDNA-CN, altered gene expression in multiple tissues, and aging-related disease.

Priors, population sizes, and power in genome-wide hypothesis tests.

BACKGROUND: Genome-wide tests, including genome-wide association studies (GWAS) of germ-line genetic variants, driver tests of cancer somatic mutations, and transcriptome-wide association tests of RNAseq data, carry a high multiple testing burden. This burden can be overcome by enrolling larger cohorts or alleviated by using prior biological knowledge to favor some hypotheses over others. Here we compare these two methods in terms of their abilities to boost the power of hypothesis testing. RESULTS: We provide a quantitative estimate for progress in cohort sizes and present a theoretical analysis of the power of oracular hard priors: priors that select a subset of hypotheses for testing, with an oracular guarantee that all true positives are within the tested subset. This theory demonstrates that for GWAS, strong priors that limit testing to 100-1000 genes provide less power than typical annual 20-40% increases in cohort sizes. Furthermore, non-oracular priors that exclude even a small fraction of true positives from the tested set can perform worse than not using a prior at all. CONCLUSION: Our results provide a theoretical explanation for the continued dominance of simple, unbiased univariate hypothesis tests for GWAS: if a statistical question can be answered by larger cohort sizes, it should be answered by larger cohort sizes rather than by more complicated biased methods involving priors. We suggest that priors are better suited for non-statistical aspects of biology, such as pathway structure and causality, that are not yet easily captured by standard hypothesis tests.

Monogenic and Polygenic Contributions to QTc Prolongation in the Population.

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.

Epigenome-wide association study of mitochondrial genome copy number.

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD.

Genetic variant annotation scores in congenital long QT syndrome.

BACKGROUND: Congenital Long QT Syndrome (LQTS) is a hereditary arrhythmic disorder. We aimed to assess the performance of current genetic variant annotation scores among LQTS patients and their predictive impact. METHODS: We evaluated 2025 patients with unique mutations for LQT1-LQT3. A patient-specific score was calculated for each of four established genetic variant annotation algorithms: CADD, SIFT, REVEL, and PolyPhen-2. The scores were tested for the identification of LQTS and their predictive performance for cardiac events (CE) and life-threatening events (LTE) and then compared with the predictive performance of LQTS categorization based on mutation location/function. Score performance was tested using Harrell's C-index. RESULTS: A total of 917 subjects were classified as LQT1, 838 as LQT2, and 270 as LQT3. The identification of a pathogenic variant occurred in 99% with CADD, 92% with SIFT, 100% with REVEL, and 86% with PolyPhen-2. However, none of the genetic scores correlated with the risk of CE (Harrell's C-index: CADD = 0.50, SIFT = 0.51, REVEL = 0.50, and PolyPhen-2 = 0.52) or LTE (Harrell's C-index: CADD = 0.50, SIFT = 0.53, REVEL = 0.54, and PolyPhen-2 = 0.52). In contrast, high-risk mutation categorization based on location/function was a powerful independent predictor of CE (HR = 1.88; p < .001) and LTE (HR = 1.89, p < .001). CONCLUSION: In congenital LQTS patients, well-established algorithms (CADD, SIFT, REVEL, and PolyPhen-2) were able to identify the majority of the causal variants as pathogenic. However, the scores did not predict clinical outcomes. These results indicate that mutation location/functional assays are essential for accurate interpretation of the risk associated with LQTS mutations.

Epigenetic Age and the Risk of Incident Atrial Fibrillation.

BACKGROUND: The most prominent risk factor for atrial fibrillation (AF) is chronological age; however, underlying mechanisms are unexplained. Algorithms using epigenetic modifications to the human genome effectively predict chronological age. Chronological and epigenetic predicted ages may diverge in a phenomenon referred to as epigenetic age acceleration (EAA), which may reflect accelerated biological aging. We sought to evaluate for associations between epigenetic age measures and incident AF. METHODS: Measures for 4 epigenetic clocks (Horvath, Hannum, DNA methylation [DNAm] PhenoAge, and DNAm GrimAge) and an epigenetic predictor of PAI-1 (plasminogen activator inhibitor-1) levels (ie, DNAm PAI-1) were determined for study participants from 3 population-based cohort studies. Cox models evaluated for associations with incident AF and results were combined via random-effects meta-analyses. Two-sample summary-level Mendelian randomization analyses evaluated for associations between genetic instruments of the EAA measures and AF. RESULTS: Among 5600 participants (mean age, 65.5 years; female, 60.1%; Black, 50.7%), there were 905 incident AF cases during a mean follow-up of 12.9 years. Unadjusted analyses revealed all 4 epigenetic clocks and the DNAm PAI-1 predictor were associated with statistically significant higher hazards of incident AF, though the magnitudes of their point estimates were smaller relative to the associations observed for chronological age. The pooled EAA estimates for each epigenetic measure, with the exception of Horvath EAA, were associated with incident AF in models adjusted for chronological age, race, sex, and smoking variables. After multivariable adjustment for additional known AF risk factors that could also potentially function as mediators, pooled EAA measures for 2 clocks remained statistically significant. Five-year increases in EAA measures for DNAm GrimAge and DNAm PhenoAge were associated with 19% (adjusted hazard ratio [HR], 1.19 [95% CI, 1.09-1.31]; P<0.01) and 15% (adjusted HR, 1.15 [95% CI, 1.05-1.25]; P<0.01) higher hazards of incident AF, respectively. Mendelian randomization analyses for the 5 EAA measures did not reveal statistically significant associations with AF. CONCLUSIONS: Our study identified adjusted associations between EAA measures and incident AF, suggesting that biological aging plays an important role independent of chronological age, though a potential underlying causal relationship remains unclear. These aging processes may be modifiable and not constrained by the immutable factor of time.

Identification of 969 protein quantitative trait loci in an African American population with kidney disease attributed to hypertension.

Investigations into the causal underpinnings of disease processes can be aided by the incorporation of genetic information. Genetic studies require populations varied in both ancestry and prevalent disease in order to optimize discovery and ensure generalizability of findings to the global population. Here, we report the genetic determinants of the serum proteome in 466 African Americans with chronic kidney disease attributed to hypertension from the richly phenotyped African American Study of Kidney Disease and Hypertension (AASK) study. Using the largest aptamer-based protein profiling platform to date (6,790 proteins or protein complexes), we identified 969 genetic associations with 900 unique proteins; including 52 novel cis (local) associations and 379 novel trans (distant) associations. The genetic effects of previously published cis-protein quantitative trait loci (pQTLs) were found to be highly reproducible, and we found evidence that our novel genetic signals colocalize with gene expression and disease processes. Many trans- pQTLs were found to reflect associations mediated by the circulating cis protein, and the common trans-pQTLs are enriched for processes involving extracellular vesicles, highlighting a plausible mechanism for distal regulation of the levels of secreted proteins. Thus, our study generates a valuable resource of genetic associations linking variants to protein levels and disease in an understudied patient population to inform future studies of drug targets and physiology.

Trans-ethnic genome-wide association study of blood metabolites in the Chronic Renal Insufficiency Cohort (CRIC) study.

Metabolomics genome wide association study (GWAS) help outline the genetic contribution to human metabolism. However, studies to date have focused on relatively healthy, population-based samples of White individuals. Here, we conducted a GWAS of 537 blood metabolites measured in the Chronic Renal Insufficiency Cohort (CRIC) Study, with separate analyses in 822 White and 687 Black study participants. Trans-ethnic meta-analysis was then applied to improve fine-mapping of potential causal variants. Mean estimated glomerular filtration rate was 44.4 and 41.5 mL/min/1.73m2 in the White and Black participants, respectively. There were 45 significant metabolite associations at 19 loci, including novel associations at PYROXD2, PHYHD1, FADS1-3, ACOT2, MYRF, FAAH, and LIPC. The strength of associations was unchanged in models additionally adjusted for estimated glomerular filtration rate and proteinuria, consistent with a direct biochemical effect of gene products on associated metabolites. At several loci, trans-ethnic meta-analysis, which leverages differences in linkage disequilibrium across populations, reduced the number and/or genomic interval spanned by potentially causal single nucleotide polymorphisms compared to fine-mapping in the White participant cohort alone. Across all validated associations, we found strong concordance in effect sizes of the potentially causal single nucleotide polymorphisms between White and Black study participants. Thus, our study identifies novel genetic determinants of blood metabolites in chronic kidney disease, demonstrates the value of diverse cohorts to improve causal inference in metabolomics GWAS, and underscores the shared genetic basis of metabolism across race.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits.

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Pleiotropic Meta-Analysis of Cognition, Education, and Schizophrenia Differentiates Roles of Early Neurodevelopmental and Adult Synaptic Pathways.

Susceptibility to schizophrenia is inversely correlated with general cognitive ability at both the phenotypic and the genetic level. Paradoxically, a modest but consistent positive genetic correlation has been reported between schizophrenia and educational attainment, despite the strong positive genetic correlation between cognitive ability and educational attainment. Here we leverage published genome-wide association studies (GWASs) in cognitive ability, education, and schizophrenia to parse biological mechanisms underlying these results. Association analysis based on subsets (ASSET), a pleiotropic meta-analytic technique, allowed jointly associated loci to be identified and characterized. Specifically, we identified subsets of variants associated in the expected ("concordant") direction across all three phenotypes (i.e., greater risk for schizophrenia, lower cognitive ability, and lower educational attainment); these were contrasted with variants that demonstrated the counterintuitive ("discordant") relationship between education and schizophrenia (i.e., greater risk for schizophrenia and higher educational attainment). ASSET analysis revealed 235 independent loci associated with cognitive ability, education, and/or schizophrenia at p < 5 × 10-8. Pleiotropic analysis successfully identified more than 100 loci that were not significant in the input GWASs. Many of these have been validated by larger, more recent single-phenotype GWASs. Leveraging the joint genetic correlations of cognitive ability, education, and schizophrenia, we were able to dissociate two distinct biological mechanisms-early neurodevelopmental pathways that characterize concordant allelic variation and adulthood synaptic pruning pathways-that were linked to the paradoxical positive genetic association between education and schizophrenia. Furthermore, genetic correlation analyses revealed that these mechanisms contribute not only to the etiopathogenesis of schizophrenia but also to the broader biological dimensions implicated in both general health outcomes and psychiatric illness.

Multi-ancestry genome-wide gene-sleep interactions identify novel loci for blood pressure.

Long and short sleep duration are associated with elevated blood pressure (BP), possibly through effects on molecular pathways that influence neuroendocrine and vascular systems. To gain new insights into the genetic basis of sleep-related BP variation, we performed genome-wide gene by short or long sleep duration interaction analyses on four BP traits (systolic BP, diastolic BP, mean arterial pressure, and pulse pressure) across five ancestry groups in two stages using 2 degree of freedom (df) joint test followed by 1df test of interaction effects. Primary multi-ancestry analysis in 62,969 individuals in stage 1 identified three novel gene by sleep interactions that were replicated in an additional 59,296 individuals in stage 2 (stage 1 + 2 Pjoint < 5 × 10-8), including rs7955964 (FIGNL2/ANKRD33) that increases BP among long sleepers, and rs73493041 (SNORA26/C9orf170) and rs10406644 (KCTD15/LSM14A) that increase BP among short sleepers (Pint < 5 × 10-8). Secondary ancestry-specific analysis identified another novel gene by long sleep interaction at rs111887471 (TRPC3/KIAA1109) in individuals of African ancestry (Pint = 2 × 10-6). Combined stage 1 and 2 analyses additionally identified significant gene by long sleep interactions at 10 loci including MKLN1 and RGL3/ELAVL3 previously associated with BP, and significant gene by short sleep interactions at 10 loci including C2orf43 previously associated with BP (Pint < 10-3). 2df test also identified novel loci for BP after modeling sleep that has known functions in sleep-wake regulation, nervous and cardiometabolic systems. This study indicates that sleep and primary mechanisms regulating BP may interact to elevate BP level, suggesting novel insights into sleep-related BP regulation.

Multiple SCN5A variant enhancers modulate its cardiac gene expression and the QT interval.

The rationale for genome-wide association study (GWAS) results is sequence variation in cis-regulatory elements (CREs) modulating a target gene's expression as the major cause of trait variation. To understand the complete molecular landscape of one of these GWAS loci, we performed in vitro reporter screens in cardiomyocyte cell lines for CREs overlapping nearly all common variants associated with any of five independent QT interval (QTi)-associated GWAS hits at the SCN5A-SCN10A locus. We identified 13 causal CRE variants using allelic reporter activity, cardiomyocyte nuclear extract-based binding assays, overlap with human cardiac tissue DNaseI hypersensitive regions, and predicted impact of sequence variants on DNaseI sensitivity. Our analyses identified at least one high-confidence causal CRE variant for each of the five sentinel hits that could collectively predict SCN5A cardiac gene expression and QTi association. Although all 13 variants could explain SCN5A gene expression, the highest statistical significance was obtained with seven variants (inclusive of the five above). Thus, multiple, causal, mutually associated CRE variants can underlie GWAS signals.

Polygenic Risk Scores for Kidney Function and Their Associations with Circulating Proteome, and Incident Kidney Diseases.

BACKGROUND: Genome-wide association studies (GWAS) have revealed numerous loci for kidney function (eGFR). The relationship between polygenic predictors of eGFR, risk of incident adverse kidney outcomes, and the plasma proteome is not known. METHODS: We developed a genome-wide polygenic risk score (PRS) for eGFR by applying the LDpred algorithm to summary statistics generated from a multiethnic meta-analysis of CKDGen Consortium GWAS ( n =765,348) and UK Biobank GWAS (90% of the cohort; n =451,508), followed by best-parameter selection using the remaining 10% of UK Biobank data ( n =45,158). We then tested the association of the PRS in the Atherosclerosis Risk in Communities (ARIC) study ( n =8866) with incident CKD, ESKD, kidney failure, and AKI. We also examined associations between the PRS and 4877 plasma proteins measured at middle age and older adulthood and evaluated mediation of PRS associations by eGFR. RESULTS: The developed PRS showed a significant association with all outcomes. Hazard ratios per 1 SD lower PRS ranged from 1.06 (95% CI, 1.01 to 1.11) to 1.33 (95% CI, 1.28 to 1.37). The PRS was significantly associated with 132 proteins at both time points. The strongest associations were with cystatin C, collagen α -1(XV) chain, and desmocollin-2. Most proteins were higher at lower kidney function, except for five proteins, including testican-2. Most correlations of the genetic PRS with proteins were mediated by eGFR. CONCLUSIONS: A PRS for eGFR is now sufficiently strong to capture risk for a spectrum of incident kidney diseases and broadly influences the plasma proteome, primarily mediated by eGFR.

Proteins Associated with Risk of Kidney Function Decline in the General Population.

BACKGROUND: Proteomic profiling may allow identification of plasma proteins that associate with subsequent changesin kidney function, elucidating biologic processes underlying the development and progression of CKD. METHODS: We quantified the association between 4877 plasma proteins and a composite outcome of ESKD or decline in eGFR by ≥50% among 9406 participants in the Atherosclerosis Risk in Communities (ARIC) Study (visit 3; mean age, 60 years) who were followed for a median of 14.4 years. We performed separate analyses for these proteins in a subset of 4378 participants (visit 5), who were followed at a later time point, for a median of 4.4 years. For validation, we evaluated proteins with significant associations (false discovery rate <5%) in both time periods in 3249 participants in the Chronic Renal Insufficiency Cohort (CRIC) and 703 participants in the African American Study of Kidney Disease and Hypertension (AASK). We also compared the genetic determinants of protein levels with those from a meta-analysis genome-wide association study of eGFR. RESULTS: In models adjusted for multiple covariates, including baseline eGFR and albuminuria, we identified 13 distinct proteins that were significantly associated with the composite end point in both time periods, including TNF receptor superfamily members 1A and 1B, trefoil factor 3, and ß-trace protein. Of these proteins, 12 were also significantly associated in CRIC, and nine were significantly associated in AASK. Higher levels of each protein associated with higher risk of 50% eGFR decline or ESKD. We found genetic evidence for a causal role for one protein, lectin mannose-binding 2 protein (LMAN2). CONCLUSIONS: Large-scale proteomic analysis identified both known and novel proteomic risk factors for eGFR decline.

Genome-wide association study of serum metabolites in the African American Study of Kidney Disease and Hypertension.

The genome-wide association study (GWAS) is a powerful means to study genetic determinants of disease traits and generate insights into disease pathophysiology. To date, few GWAS of circulating metabolite levels have been performed in African Americans with chronic kidney disease. Hypothesizing that novel genetic-metabolite associations may be identified in a unique population of African Americans with a lower glomerular filtration rate (GFR), we conducted a GWAS of 652 serum metabolites in 619 participants (mean measured glomerular filtration rate 45 mL/min/1.73m2) in the African American Study of Kidney Disease and Hypertension, a clinical trial of blood pressure lowering and antihypertensive medication in African Americans with chronic kidney disease. We identified 42 significant variant metabolite associations. Twenty associations had been previously identified in published GWAS, and eleven novel associations were replicated in a separate cohort of 818 African Americans with genetic and metabolomic data from the Atherosclerosis Risk in Communities Study. The replicated novel variant-metabolite associations comprised eight metabolites and eleven distinct genomic loci. Nine of the replicated associations represented clear enzyme-metabolite interactions, with high expression in the kidneys as well as the liver. Three loci (ACY1, ACY3, and NAT8) were associated with a common pool of metabolites, acetylated amino acids, but with different individual affinities. Thus, extensive metabolite profiling in an African American population with chronic kidney disease aided identification of novel genome-wide metabolite associations, providing clues about substrate specificity and the key roles of enzymes in modulating systemic levels of metabolites.

Functionally oriented analysis of cardiometabolic traits in a trans-ethnic sample.

Interpretation of genetic association results is difficult because signals often lack biological context. To generate hypotheses of the functional genetic etiology of complex cardiometabolic traits, we estimated the genetically determined component of gene expression from common variants using PrediXcan (1) and determined genes with differential predicted expression by trait. PrediXcan imputes tissue-specific expression levels from genetic variation using variant-level effect on gene expression in transcriptome data. To explore the value of imputed genetically regulated gene expression (GReX) models across different ancestral populations, we evaluated imputed expression levels for predictive accuracy genome-wide in RNA sequence data in samples drawn from European-ancestry and African-ancestry populations and identified substantial predictive power using European-derived models in a non-European target population. We then tested the association of GReX on 15 cardiometabolic traits including blood lipid levels, body mass index, height, blood pressure, fasting glucose and insulin, RR interval, fibrinogen level, factor VII level and white blood cell and platelet counts in 15 755 individuals across three ancestry groups, resulting in 20 novel gene-phenotype associations reaching experiment-wide significance across ancestries. In addition, we identified 18 significant novel gene-phenotype associations in our ancestry-specific analyses. Top associations were assessed for additional support via query of S-PrediXcan (2) results derived from publicly available genome-wide association studies summary data. Collectively, these findings illustrate the utility of transcriptome-based imputation models for discovery of cardiometabolic effect genes in a diverse dataset.

Variation in a Left Ventricle-Specific Hand1 Enhancer Impairs GATA Transcription Factor Binding and Disrupts Conduction System Development and Function.

RATIONALE: The ventricular conduction system (VCS) rapidly propagates electrical impulses through the working myocardium of the ventricles to coordinate chamber contraction. GWAS (Genome-wide association studies) have associated nucleotide polymorphisms, most are located within regulatory intergenic or intronic sequences, with variation in VCS function. Two highly correlated polymorphisms (r2>0.99) associated with VCS functional variation (rs13165478 and rs13185595) occur 5' to the gene encoding the basic helix-loop-helix transcription factor HAND1 (heart- and neural crest derivatives-expressed protein 1). OBJECTIVE: Here, we test the hypothesis that these polymorphisms influence HAND1 transcription thereby influencing VCS development and function. METHODS AND RESULTS: We employed transgenic mouse models to identify an enhancer that is sufficient for left ventricle (LV) cis-regulatory activity. Two evolutionarily conserved GATA transcription factor cis-binding elements within this enhancer are bound by GATA4 and are necessary for cis-regulatory activity, as shown by in vitro DNA binding assays. CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9-mediated deletion of this enhancer dramatically reduces Hand1 expression solely within the LV but does not phenocopy previously published mouse models of cardiac Hand1 loss-of-function. Electrophysiological and morphological analyses reveals that mice homozygous for this deleted enhancer display a morphologically abnormal VCS and a conduction system phenotype consistent with right bundle branch block. Using 1000 Genomes Project data, we identify 3 additional single nucleotide polymorphisms (SNPs), located within the Hand1 LV enhancer, that compose a haplotype with rs13165478 and rs13185595. One of these SNPs, rs10054375, overlaps with a critical GATA cis-regulatory element within the Hand1 LV enhancer. This SNP, when tested in electrophoretic mobility shift assays, disrupts GATA4 DNA-binding. Modeling 2 of these SNPs in mice causes diminished Hand1 expression and mice present with abnormal VCS function. CONCLUSIONS: Together, these findings reveal that SNP rs10054375, which is located within a necessary and sufficient LV-specific Hand1 enhancer, exhibits reduces GATA DNA-binding in electrophoretic mobility shift assay, and this enhancer in total, is required for VCS development and function in mice and perhaps humans.