Search

Whole-genome sequencing (WGS) and whole-exome sequencing studies have become increasingly available and are being used to identify rare genetic variants associated with health and disease outcomes. Investigators routinely use mixed models to account for genetic relatedness or other clustering variables (e.g., family or household) when testing genetic associations. However, no existing tests of the association of a rare variant with a binary outcome in the presence of correlated data control the type 1 error where there are (1) few individuals harboring the rare allele, (2) a small proportion of cases relative to controls, and (3) covariates to adjust for. Here, we address all three issues in developing a framework for testing rare variant association with a binary trait in individuals harboring at least one risk allele. In this framework, we estimate outcome probabilities under the null hypothesis and then use them, within the individuals with at least one risk allele, to test variant associations. We extend the BinomiRare test, which was previously proposed for independent observations, and develop the Conway-Maxwell-Poisson (CMP) test and study their properties in simulations. We show that the BinomiRare test always controls the type 1 error, while the CMP test sometimes does not. We then use the BinomiRare test to test the association of rare genetic variants in target genes with small-vessel disease (SVD) stroke, short sleep, and venous thromboembolism (VTE), in whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program.

BackgroundWe examined how expanding electrocardiographic trait genome-wide association studies to include ancestrally diverse populations, prioritize more precise phenotypic measures, and evaluate evidence for shared genetic effects enabled the detection and characterization of loci.MethodsWe decomposed 10 seconds, 12-lead electrocardiograms from 34 668 multi-ethnic participants (15% Black; 30% Hispanic/Latino) into 6 contiguous, physiologically distinct (P wave, PR segment, QRS interval, ST segment, T wave, and TP segment) and 2 composite, conventional (PR interval and QT interval) interval scale traits and conducted multivariable-adjusted, trait-specific univariate genome-wide association studies using 1000-G imputed single-nucleotide polymorphisms. Evidence of shared genetic effects was evaluated by aggregating meta-analyzed univariate results across the 6 continuous electrocardiographic traits using the combined phenotype adaptive sum of powered scores test.ResultsWe identified 6 novels (CD36, PITX2, EMB, ZNF592, YPEL2, and BC043580) and 87 known loci (adaptive sum of powered score test P

BackgroundThe electrocardiographically quantified QRS duration measures ventricular depolarization and conduction. QRS prolongation has been associated with poor heart failure prognosis and cardiovascular mortality, including sudden death. While previous genome-wide association studies (GWAS) have identified 32 QRS SNPs across 26 loci among European, African, and Asian-descent populations, the genetics of QRS among Hispanics/Latinos has not been previously explored.MethodsWe performed a GWAS of QRS duration among Hispanic/Latino ancestry populations (n = 15,124) from four studies using 1000 Genomes imputed genotype data (adjusted for age, sex, global ancestry, clinical and study-specific covariates). Study-specific results were combined using fixed-effects, inverse variance-weighted meta-analysis.ResultsWe identified six loci associated with QRS (P

ObjectivePR interval (PR) is a heritable electrocardiographic measure of atrial and atrioventricular nodal conduction. Changes in PR duration may be associated with atrial fibrillation, heart failure and all-cause mortality. Hispanic/Latino populations have high burdens of cardiovascular morbidity and mortality, are highly admixed and represent exceptional opportunities for novel locus identification. However, they remain chronically understudied. We present the first genome-wide association study (GWAS) of PR in 14 756 participants of Hispanic/Latino ancestry from three studies.MethodsStudy-specific summary results of the association between 1000 Genomes Phase 1 imputed single-nucleotide polymorphisms (SNPs) and PR assumed an additive genetic model and were adjusted for global ancestry, study centre/region and clinical covariates. Results were combined using fixed-effects, inverse variance weighted meta-analysis. Sequential conditional analyses were used to identify independent signals. Replication of novel loci was performed in populations of Asian, African and European descent. ENCODE and RoadMap data were used to annotate results.ResultsWe identified a novel genome-wide association (P

The genetic basis of supraventricular and ventricular ectopy (SVE, VE) remains largely uncharacterized, despite established genetic mechanisms of arrhythmogenesis. To identify novel genetic variants associated with SVE/VE in ancestrally diverse human populations, we conducted a genome-wide association study of electrocardiographically identified SVE and VE in five cohorts including approximately 43,000 participants of African, European and Hispanic/Latino ancestry. In thirteen ancestry-stratified subgroups, we tested multivariable-adjusted associations of SVE and VE with single nucleotide polymorphism (SNP) dosage. We combined subgroup-specific association estimates in inverse variance-weighted, fixed-effects and Bayesian meta-analyses. We also combined fixed-effects meta-analytic t-test statistics for SVE and VE in multi-trait SNP association analyses. No loci reached genome-wide significance in trans-ethnic meta-analyses. However, we found genome-wide significant SNPs intronic to an apoptosis-enhancing gene previously associated with QRS interval duration (FAF1; lead SNP rs7545860; effect allele frequency = 0.02; P = 2.0 × 10-8) in multi-trait analysis among European ancestry participants and near a locus encoding calcium-dependent glycoproteins (DSC3; lead SNP rs8086068; effect allele frequency = 0.17) in meta-analysis of SVE (P = 4.0 × 10-8) and multi-trait analysis (P = 2.9 × 10-9) among African ancestry participants. The novel findings suggest several mechanisms by which genetic variation may predispose to ectopy in humans and highlight the potential value of leveraging pleiotropy in future studies of ectopy-related phenotypes.

BackgroundIncreased heart rate and a prolonged QT interval are important risk factors for cardiovascular morbidity and mortality, and can be influenced by the use of various medications, including tricyclic/tetracyclic antidepressants (TCAs). We aim to identify genetic loci that modify the association between TCA use and RR and QT intervals.Methods and resultsWe conducted race/ethnic-specific genome-wide interaction analyses (with HapMap phase II imputed reference panel imputation) of TCAs and resting RR and QT intervals in cohorts of European (n=45 706; n=1417 TCA users), African (n=10 235; n=296 TCA users) and Hispanic/Latino (n=13 808; n=147 TCA users) ancestry, adjusted for clinical covariates. Among the populations of European ancestry, two genome-wide significant loci were identified for RR interval: rs6737205 in BRE (β=56.3, pinteraction=3.9e-9) and rs9830388 in UBE2E2 (β=25.2, pinteraction=1.7e-8). In Hispanic/Latino cohorts, rs2291477 in TGFBR3 significantly modified the association between TCAs and QT intervals (β=9.3, pinteraction=2.55e-8). In the meta-analyses of the other ethnicities, these loci either were excluded from the meta-analyses (as part of quality control), or their effects did not reach the level of nominal statistical significance (pinteraction>0.05). No new variants were identified in these ethnicities. No additional loci were identified after inverse-variance-weighted meta-analysis of the three ancestries.ConclusionsAmong Europeans, TCA interactions with variants in BRE and UBE2E2 were identified in relation to RR intervals. Among Hispanic/Latinos, variants in TGFBR3 modified the relation between TCAs and QT intervals. Future studies are required to confirm our results.

BackgroundThe electrocardiographically measured QT interval (QT) is heritable and its prolongation is an established risk factor for several cardiovascular diseases. Yet, most QT genetic studies have been performed in European ancestral populations, possibly reducing their global relevance.ObjectiveTo leverage diversity and improve biological insight, we fine mapped 16 of the 35 previously identified QT loci (46%) in populations of African American (n = 12,410) and Hispanic/Latino (n = 14,837) ancestry.MethodsRacial/ethnic-specific multiple linear regression analyses adjusted for heart rate and clinical covariates were examined separately and in combination after inverse-variance weighted trans-ethnic meta-analysis.ResultsThe 16 fine-mapped QT loci included on the Illumina Metabochip represented 21 independent signals, of which 16 (76%) were significantly (P-value≤9.1×10-5) associated with QT. Through sequential conditional analysis we also identified three trans-ethnic novel SNPs at ATP1B1, SCN5A-SCN10A, and KCNQ1 and three Hispanic/Latino-specific novel SNPs at NOS1AP and SCN5A-SCN10A (two novel SNPs) with evidence of associations with QT independent of previous identified GWAS lead SNPs. Linkage disequilibrium patterns helped to narrow the region likely to contain the functional variants at several loci, including NOS1AP, USP50-TRPM7, and PRKCA, although intervals surrounding SLC35F1-PLN and CNOT1 remained broad in size (>100 kb). Finally, bioinformatics-based functional characterization suggested a regulatory function in cardiac tissues for the majority of independent signals that generalized and the novel SNPs.ConclusionOur findings suggest that a majority of identified SNPs implicate gene regulatory dysfunction in QT prolongation, that the same loci influence variation in QT across global populations, and that additional, novel, population-specific QT signals exist.

QT interval prolongation is a heritable risk factor for ventricular arrhythmias and can predispose to sudden death. Most genome-wide association studies (GWAS) of QT were performed in European ancestral populations, leaving other groups uncharacterized. Herein we present the first QT GWAS of Hispanic/Latinos using data on 15,997 participants from four studies. Study-specific summary results of the association between 1000 Genomes Project (1000G) imputed SNPs and electrocardiographically measured QT were combined using fixed-effects meta-analysis. We identified 41 genome-wide significant SNPs that mapped to 13 previously identified QT loci. Conditional analyses distinguished six secondary signals at NOS1AP (n = 2), ATP1B1 (n = 2), SCN5A (n = 1), and KCNQ1 (n = 1). Comparison of linkage disequilibrium patterns between the 13 lead SNPs and six secondary signals with previously reported index SNPs in 1000G super populations suggested that the SCN5A and KCNE1 lead SNPs were potentially novel and population-specific. Finally, of the 42 suggestively associated loci, AJAP1 was suggestively associated with QT in a prior East Asian GWAS; in contrast BVES and CAP2 murine knockouts caused cardiac conduction defects. Our results indicate that whereas the same loci influence QT across populations, population-specific variation exists, motivating future trans-ethnic and ancestrally diverse QT GWAS.

US Hispanic/Latino individuals are diverse in genetic ancestry, culture, and environmental exposures. Here, we characterized and controlled for this diversity in genome-wide association studies (GWASs) for the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). We simultaneously estimated population-structure principal components (PCs) robust to familial relatedness and pairwise kinship coefficients (KCs) robust to population structure, admixture, and Hardy-Weinberg departures. The PCs revealed substantial genetic differentiation within and among six self-identified background groups (Cuban, Dominican, Puerto Rican, Mexican, and Central and South American). To control for variation among groups, we developed a multi-dimensional clustering method to define a "genetic-analysis group" variable that retains many properties of self-identified background while achieving substantially greater genetic homogeneity within groups and including participants with non-specific self-identification. In GWASs of 22 biomedical traits, we used a linear mixed model (LMM) including pairwise empirical KCs to account for familial relatedness, PCs for ancestry, and genetic-analysis groups for additional group-associated effects. Including the genetic-analysis group as a covariate accounted for significant trait variation in 8 of 22 traits, even after we fit 20 PCs. Additionally, genetic-analysis groups had significant heterogeneity of residual variance for 20 of 22 traits, and modeling this heteroscedasticity within the LMM reduced genomic inflation for 19 traits. Furthermore, fitting an LMM that utilized a genetic-analysis group rather than a self-identified background group achieved higher power to detect previously reported associations. We expect that the methods applied here will be useful in other studies with multiple ethnic groups, admixture, and relatedness.

Thiazide diuretics, commonly used antihypertensives, may cause QT interval (QT) prolongation, a risk factor for highly fatal and difficult to predict ventricular arrhythmias. We examined whether common SNPs modified the association between thiazide use and QT or its component parts (QRS interval, JT interval) by performing ancestry-specific, trans-ethnic, and cross-phenotype genome-wide analyses of European (66%), African American (15%), and Hispanic (19%) populations (N=78,199), leveraging longitudinal data, incorporating corrected standard errors to account for underestimation of interaction estimate variances and evaluating evidence for pathway enrichment. Although no loci achieved genome-wide significance (P<5×10−8), we found suggestive evidence (P<5×10−6) for SNPs modifying the thiazide-QT association at 22 loci, including ion transport loci (e.g. NELL1, KCNQ3). The biologic plausibility of our suggestive results and simulations demonstrating modest power to detect interaction effects at genome-wide significant levels indicate that larger studies and innovative statistical methods are warranted in future efforts evaluating thiazide-SNP interactions.

Sulfonylureas, a commonly-used class of medication used to treat type 2 diabetes, have been associated with an increased risk of cardiovascular disease. Their effects on QT interval duration and related electrocardiographic phenotypes are potential mechanisms for this adverse effect. In eleven ethnically diverse cohorts that included 71 857 European, African American, and Hispanic/Latino ancestry individuals with repeated measures of medication use and electrocardiogram (ECG) measurements, we conducted a pharmacogenomic genome-wide association study of sulfonylurea use and three ECG phenotypes: QT, JT, and QRS intervals. In ancestry-specific meta-analyses, 8 novel pharmacogenomic loci met the threshold for genome-wide significance (P < 5 x 10−8), and a pharmacokinetic variant in CYP2C9 (rs1057910) that has been associated with sulfonylurea-related treatment effects and other adverse drug reactions in previous studies was replicated. Additional research is needed to replicate the novel findings and to understand their biological basis.

Background: Activated partial thromboplastin time (aPTT) is a clinical test used to measure the clotting time between the activation of factor XII and the formation of a fibrin clot. Prolonged aPTT indicates a deficiency in the coagulation pathway while shortened aPTT is associated with prothrombotic risk factors and venus thromboembolism. Despite the high heritability of aPTT (40-80%), previous genome-wide association studies (GWAS) of aPTT have been limited to European descent populations, with only a smaller candidate gene study conducted in African Americans. Methods: We included 13,803 participants of European ancestry (EU) and 2,724 participants of African American ancestry (AA) from the Cohorts for Heart and Aging Research in Genetic Epidemiology (CHARGE) consortium. aPTT, in seconds (s), was measured using standard protocols in plasma with the use of automated coagulometers. Genotype data were imputed to the 1000 Genomes Phase 1 reference panel, and associations were examined using linear regression assuming an additive genetic model and adjusting for global ancestry, age, sex, and study design. Study-specific results were combined using a fixed-effects, inverse variance weighted meta-analysis. Results: We identified seven loci associated with aPTT in EU populations (F5, FRMD5, KNG1, F11, F12, HLA, ABO) and three loci associated with aPTT in AA populations (KNG1, F12, ABO) at genome-wide significant levels (P < 5x10-8). These results are consistent with previously reported genetic studies in EU and AA populations. Effect sizes were larger in AA populations (1.08 to 1.32 s) than in EU populations (0.40 to 1.00 s). Conclusions: We successfully replicated associations with aPTT at seven loci in EU populations and three loci in AA populations. Our results suggest that genetic determinants of aPTT are consistent across race/ethnicity but that studies in AA populations are currently underpowered relative to EU populations. Future work in aPTT genetics should consider more diverse populations. [ABSTRACT FROM AUTHOR]

PR interval (PR) prolongation is a predictor of atrial fibrillation, a common cardiac arrhythmia, and is an established risk factor for pacemaker implantation, heart failure, stroke, and all-cause mortality. Previous genome-wide association studies in European (EU) and African (AA) descent populations have identified multiple loci associated with PR. However, it is unclear whether these loci are relevant in other racial groups, including Hispanic/Latinos (HL). Extending studies to include admixed populations such as HL also provides the potential to narrow regions associated with PR and identify population specific variants. Therefore, we evaluated 1,135 SNPs from four previously identified PR loci (SCN5A, SCN10A, CAV1-CAV2, and TBX5-TBX3), fine-mapped on the Illumina Metabochip for association with PR. Associations were examined using linear regression assuming an additive genetic model and adjusting for global ancestry and clinical covariates (age, sex, RR interval, body mass index, height, systolic blood pressure) in 11,800 HL participants from participating studies of the Population Architecture using Genomics and Epidemiology (PAGE) consortia. Three of the four loci (SCN5A, SCN10A, and CAV1-CAV2) generalized to HL populations (P<7.1E-4 based on the number of SNPs correlated at r²>0.2 with EU-identified index SNPs). At each of these loci, we identified SNPs with stronger evidence of association with PR than the EU-identified index SNP when examined in the HL population (SCN5A: rs7374540; SCN10A: rs6801957; CAV-CAV2: rs717957). Conditional analyses also identified two additional independent SNPs at each of the three generalized loci (SCN5A: rs3796387, rs9861242; SCN10A: rs10428132, rs7433723; CAV1-CAV2: rs3801995, rs3807994), representing potential secondary signals. Furthermore, linkage disequilibrium (LD) patterns in Hispanics enabled the narrowing of the regions flanking the EU-identified index SNPs. For example, at the CAV1-CAV2 locus, HL LD patterns indicated that seven SNPs were in LD with rs3807989, the EU-identified index SNP, compared to 17 SNPs that were in LD with rs3807989 when applying EU LD patterns, reducing the size of the region by 71 kilobases. Finally, we examined evidence for population specific signals outside of the specific regions identified in EU populations which were associated with PR in Hispanics by focusing on SNPs uncorrelated with the EU index signals. At all four loci, multiple SNPs had evidence of association independent of the signals examined based on EU results. In conclusion, the same genetic loci that influence PR in European descent populations are also important in Hispanic populations, although we detected evidence for population specific SNPs. These results emphasize the importance of examining genetic associations in diverse populations in order to narrow the genomic interval for future functional interrogation. [ABSTRACT FROM AUTHOR]

Background: Risk for venous thromboembolism has a strong genetic component. Whole genome sequencing from the TOPMed program (Trans-Omics for Precision Medicine) allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies., Methods: The 3793 cases and 7834 controls (11.6% of cases were individuals of African, Hispanic/Latino, or Asian ancestry) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants)., Results: Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only PROC (odds ratio, 6.2 for carriers of rare variants; P =7.4×10 -14 ) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at PROC (odds ratio, 3.8; P =1.6×10 -14 ), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: PROS1 became significant (minimum P =1.8×10 -6 with the secondary filter), while SERPINC1 did not (minimum P =4.4×10 -5 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, MS4A1 , became significant ( P =4.4×10 -7 using all missense variants with minor allele frequency <0.0005)., Conclusions: Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel MS4A1 locus and to identify additional rare variation associated with venous thromboembolism.