Publications

Genome-wide association analysis advanced understanding of blood pressure (BP), a major risk factor for vascular conditions such as coronary heart disease and stroke. Accounting for smoking behavior may help identify BP loci and extend our knowledge of its genetic architecture. We performed genome-wide association meta-analyses of systolic and diastolic BP incorporating gene-smoking interactions in 610,091 individuals. Stage 1 analysis examined ∼18.8 million SNPs and small insertion/deletion variants in 129,913 individuals from four ancestries (European, African, Asian, and Hispanic) with follow-up analysis of promising variants in 480,178 additional individuals from five ancestries. We identified 15 loci that were genome-wide significant (p < 5 × 10-8) in stage 1 and formally replicated in stage 2. A combined stage 1 and 2 meta-analysis identified 66 additional genome-wide significant loci (13, 35, and 18 loci in European, African, and trans-ancestry, respectively). A total of 56 known BP loci were also identified by our results (p < 5 × 10-8). Of the newly identified loci, ten showed significant interaction with smoking status, but none of them were replicated in stage 2. Several loci were identified in African ancestry, highlighting the importance of genetic studies in diverse populations. The identified loci show strong evidence for regulatory features and support shared pathophysiology with cardiometabolic and addiction traits. They also highlight a role in BP regulation for biological candidates such as modulators of vascular structure and function (CDKN1B, BCAR1-CFDP1, PXDN, EEA1), ciliopathies (SDCCAG8, RPGRIP1L), telomere maintenance (TNKS, PINX1, AKTIP), and central dopaminergic signaling (MSRA, EBF2).

Obstructive sleep apnea (OSA) is a common heritable disorder displaying marked sexual dimorphism in disease prevalence and progression. Previous genetic association studies have identified a few genetic loci associated with OSA and related quantitative traits, but they have only focused on single ethnic groups, and a large proportion of the heritability remains unexplained. The apnea-hypopnea index (AHI) is a commonly used quantitative measure characterizing OSA severity. Because OSA differs by sex, and the pathophysiology of obstructive events differ in rapid eye movement (REM) and non-REM (NREM) sleep, we hypothesized that additional genetic association signals would be identified by analyzing the NREM/REM-specific AHI and by conducting sex-specific analyses in multiethnic samples. We performed genome-wide association tests for up to 19,733 participants of African, Asian, European, and Hispanic/Latino American ancestry in 7 studies. We identified rs12936587 on chromosome 17 as a possible quantitative trait locus for NREM AHI in men (N = 6,737; P = 1.7 × 10-8) but not in women (P = 0.77). The association with NREM AHI was replicated in a physiological research study (N = 67; P = 0.047). This locus overlapping the RAI1 gene and encompassing genes PEMT1, SREBF1, and RASD1 was previously reported to be associated with coronary artery disease, lipid metabolism, and implicated in Potocki-Lupski syndrome and Smith-Magenis syndrome, which are characterized by abnormal sleep phenotypes. We also identified gene-by-sex interactions in suggestive association regions, suggesting that genetic variants for AHI appear to vary by sex, consistent with the clinical observations of strong sexual dimorphism.

INTRODUCTION: Genetic variants associated with nicotine dependence have previously been identified, primarily in European-ancestry populations. No genome-wide association studies (GWAS) have been reported for smoking behaviors in Hispanics/Latinos in the United States and Latin America, who are of mixed ancestry with European, African, and American Indigenous components.

METHODS: We examined genetic associations with smoking behaviors in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) (N = 12 741 with smoking data, 5119 ever-smokers), using  2.3 million genotyped variants imputed to the 1000 Genomes Project phase 3. Mixed logistic regression models accounted for population structure, sampling, relatedness, sex, and age.

RESULTS: The known region of CHRNA5, which encodes the α5 cholinergic nicotinic receptor subunit, was associated with heavy smoking at genome-wide significance (p ≤ 5 × 10-8) in a comparison of 1929 ever-smokers reporting cigarettes per day (CPD) > 10 versus 3156 reporting CPD ≤ 10. The functional variant rs16969968 in CHRNA5 had a p value of 2.20 × 10-7 and odds ratio (OR) of 1.32 for the minor allele (A); its minor allele frequency was 0.22 overall and similar across Hispanic/Latino background groups (Central American = 0.17; South American = 0.19; Mexican = 0.18; Puerto Rican = 0.22; Cuban = 0.29; Dominican = 0.19). CHRNA4 on chromosome 20 attained p < 10-4, supporting prior findings in non-Hispanics. For nondaily smoking, which is prevalent in Hispanic/Latino smokers, compared to daily smoking, loci on chromosomes 2 and 4 achieved genome-wide significance; replication attempts were limited by small Hispanic/Latino sample sizes.

CONCLUSIONS: Associations of nicotinic receptor gene variants with smoking, first reported in non-Hispanic European-ancestry populations, generalized to Hispanics/Latinos despite different patterns of smoking behavior.

IMPLICATIONS: We conducted the first large-scale genome-wide association study (GWAS) of smoking behavior in a US Hispanic/Latino cohort, and the first GWAS of daily/nondaily smoking in any population. Results show that the region of the nicotinic receptor subunit gene CHRNA5, which in non-Hispanic European-ancestry smokers has been associated with heavy smoking as well as cessation and treatment efficacy, is also significantly associated with heavy smoking in this Hispanic/Latino cohort. The results are an important addition to understanding the impact of genetic variants in understudied Hispanic/Latino smokers.

Although genome-wide association studies (GWAS) have identified several variants linked to depression, few GWAS of non-European populations have been performed. We conducted a genome-wide analysis of depression in a large, population-based sample of Hispanics/Latinos. Data came from 12,310 adults in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). Past-week depressive symptoms were assessed using the 10-item Center for Epidemiological Studies of Depression Scale. Three phenotypes were examined: a total depression score, a total score modified to account for psychiatric medication use, and a score excluding anti-depressant medication users. We estimated heritability due to common variants (h2SNP), and performed a GWAS of the three phenotypes. Replication was attempted in three independent Hispanic/Latino cohorts. We also performed sex-stratified analyses, analyzed a binary trait indicating probable depression, and conducted three trans-ethnic analyses. The three phenotypes exhibited significant heritability (h2SNP = 6.3-6.9%; p = .002) in the total sample. No SNPs were genome-wide significant in analyses of the three phenotypes or the binary indicator of probable depression. In sex-stratified analyses, seven genome-wide significant SNPs (one in females; six in males) were identified, though none were supported through replication. Four out of 24 loci identified in prior GWAS were nominally associated in HCHS/SOL. There was no evidence of overlap in genetic risk factors across ancestry groups, though this may have been due to low power. We conducted the largest GWAS of depression-related phenotypes in Hispanic/Latino adults. Results underscore the genetic complexity of depressive symptoms as a phenotype in this population and suggest the need for much larger samples.

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 single-nucleotide polymorphisms (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 (for example, 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.

[This corrects the article DOI: 10.1371/journal.pgen.1006728.].

Heavy alcohol consumption is an established risk factor for hypertension; the mechanism by which alcohol consumption impact blood pressure (BP) regulation remains unknown. We hypothesized that a genome-wide association study accounting for gene-alcohol consumption interaction for BP might identify additional BP loci and contribute to the understanding of alcohol-related BP regulation. We conducted a large two-stage investigation incorporating joint testing of main genetic effects and single nucleotide variant (SNV)-alcohol consumption interactions. In Stage 1, genome-wide discovery meta-analyses in ≈131K individuals across several ancestry groups yielded 3,514 SNVs (245 loci) with suggestive evidence of association (P < 1.0 x 10-5). In Stage 2, these SNVs were tested for independent external replication in ≈440K individuals across multiple ancestries. We identified and replicated (at Bonferroni correction threshold) five novel BP loci (380 SNVs in 21 genes) and 49 previously reported BP loci (2,159 SNVs in 109 genes) in European ancestry, and in multi-ancestry meta-analyses (P < 5.0 x 10-8). For African ancestry samples, we detected 18 potentially novel BP loci (P < 5.0 x 10-8) in Stage 1 that warrant further replication. Additionally, correlated meta-analysis identified eight novel BP loci (11 genes). Several genes in these loci (e.g., PINX1, GATA4, BLK, FTO and GABBR2) have been previously reported to be associated with alcohol consumption. These findings provide insights into the role of alcohol consumption in the genetic architecture of hypertension.

Studies using self-reported data suggest a gene-physical activity interaction on obesity, yet the influence of sedentary behavior, distinct from a lack of physical activity, on genetic associations with obesity remains unclear. We analyzed interactions of accelerometer-measured moderate to vigorous physical activity (MVPA) and time spent sedentary with genetic variants on obesity among 9,645 U.S. Hispanics/Latinos. An overall genetic risk score (GRS), a central nervous system (CNS)-related GRS, and a non-CNS-related GRS were calculated based on 97 BMI-associated single nucleotide polymorphisms (SNPs). Genetic association with BMI was stronger in individuals with lower MVPA (first tertile) versus higher MVPA (third tertile) (β = 0.78 kg/m2 [SE, 0.10 kg/m2] vs. 0.39 kg/m2 [0.09 kg/m2] per SD increment of GRS; Pinteraction = 0.005), and in those with more time spent sedentary (third tertile) versus less time spent sedentary (first tertile) (β = 0.73 kg/m2 [SE, 0.10 kg/m2] vs. 0.44 kg/m2 [0.09 kg/m2]; Pinteraction = 0.006). Similar significant interaction patterns were observed for obesity risk, body fat mass, fat percentage, fat mass index, and waist circumference, but not for fat-free mass. The CNS-related GRS, but not the non-CNS-related GRS, showed significant interactions with MVPA and sedentary behavior, with effects on BMI and other adiposity traits. Our data suggest that both increasing physical activity and reducing sedentary behavior may attenuate genetic associations with obesity, although the independence of these interaction effects needs to be investigated further.

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.

Most regression-based tests of the association between a low-count variant and a binary outcome do not protect type 1 error, especially when tests are rejected based on a very low significance threshold. Noted exception is the Firth test. However, it was recently shown that in meta-analyzing multiple studies all asymptotic, regression-based tests, including the Firth, may not control type 1 error in some settings, and the Firth test may suffer a substantial loss of power. The problem is exacerbated when the case-control proportions differ between studies. I propose the BinomiRare exact test that circumvents the calibration problems of regression-based estimators. It quantifies the strength of association between the variant and the disease outcome based on the departure of the number of diseased individuals carrying the variant from the expected distribution of disease probability, under the null hypothesis of no association between the disease outcome and the rare variant. I provide a meta-analytic strategy to combine tests across multiple cohorts, which requires that each cohort provides the disease probabilities of all carriers of the variant in question, and the number of diseased individuals among the carriers. I show that BinomiRare controls type 1 error in meta-analysis even when the case-control proportions differ between the studies, and does not lose power compared to pooled analysis. I demonstrate the test in studying the association of rare variants with asthma in the Hispanic Community Health Study/Study of Latinos.