Ancestry- and sex-specific effects underlying inguinal hernia susceptibility identified in a multiethnic genome-wide association study meta-analysis.

Inguinal hernias are some of the most frequently diagnosed conditions in clinical practice and inguinal hernia repair is the most common procedure performed by general surgeons. Studies of inguinal hernias in non-European populations are lacking, though it is expected that such studies could identify novel loci. Further, the cumulative lifetime incidence of inguinal hernia is nine times greater in men than women, however, it is not clear why this difference exists. We conducted a genome-wide association meta-analysis of inguinal hernia risk across 513 120 individuals (35 774 cases and 477 346 controls) of Hispanic/Latino, African, Asian and European descent, with replication in 728 418 participants (33 491 cases and 694 927 controls) from the 23andMe, Inc dataset. We identified 63 genome-wide significant loci (P < 5 × 10-8), including 41 novel. Ancestry-specific analyses identified two loci (LYPLAL1-AS1/SLC30A10 and STXBP6-NOVA1) in African ancestry individuals. Sex-stratified analyses identified two loci (MYO1D and ZBTB7C) that are specific to women, and four (EBF2, EMX2/RAB11FIP2, VCL and FAM9A/FAM9B) that are specific to men. Functional experiments demonstrated that several of the associated regions (EFEMP1 and LYPLAL1-SLC30A10) function as enhancers and show differential activity between risk and reference alleles. Our study highlights the importance of large-scale genomic studies in ancestrally diverse populations for identifying ancestry-specific inguinal hernia susceptibility loci and provides novel biological insights into inguinal hernia etiology.

Multitrait genetic association analysis identifies 50 new risk loci for gastro-oesophageal reflux, seven new loci for Barrett's oesophagus and provides insights into clinical heterogeneity in reflux diagnosis.

OBJECTIVE: Gastro-oesophageal reflux disease (GERD) has heterogeneous aetiology primarily attributable to its symptom-based definitions. GERD genome-wide association studies (GWASs) have shown strong genetic overlaps with established risk factors such as obesity and depression. We hypothesised that the shared genetic architecture between GERD and these risk factors can be leveraged to (1) identify new GERD and Barrett's oesophagus (BE) risk loci and (2) explore potentially heterogeneous pathways leading to GERD and oesophageal complications. DESIGN: We applied multitrait GWAS models combining GERD (78 707 cases; 288 734 controls) and genetically correlated traits including education attainment, depression and body mass index. We also used multitrait analysis to identify BE risk loci. Top hits were replicated in 23andMe (462 753 GERD cases, 24 099 BE cases, 1 484 025 controls). We additionally dissected the GERD loci into obesity-driven and depression-driven subgroups. These subgroups were investigated to determine how they relate to tissue-specific gene expression and to risk of serious oesophageal disease (BE and/or oesophageal adenocarcinoma, EA). RESULTS: We identified 88 loci associated with GERD, with 59 replicating in 23andMe after multiple testing corrections. Our BE analysis identified seven novel loci. Additionally we showed that only the obesity-driven GERD loci (but not the depression-driven loci) were associated with genes enriched in oesophageal tissues and successfully predicted BE/EA. CONCLUSION: Our multitrait model identified many novel risk loci for GERD and BE. We present strong evidence for a genetic underpinning of disease heterogeneity in GERD and show that GERD loci associated with depressive symptoms are not strong predictors of BE/EA relative to obesity-driven GERD loci.

Insights into the genetic basis of retinal detachment.

Retinal detachment (RD) is a serious and common condition, but genetic studies to date have been hampered by the small size of the assembled cohorts. In the UK Biobank data set, where RD was ascertained by self-report or hospital records, genetic correlations between RD and high myopia or cataract operation were, respectively, 0.46 (SE = 0.08) and 0.44 (SE = 0.07). These correlations are consistent with known epidemiological associations. Through meta-analysis of genome-wide association studies using UK Biobank RD cases (N = 3 977) and two cohorts, each comprising ~1 000 clinically ascertained rhegmatogenous RD patients, we uncovered 11 genome-wide significant association signals. These are near or within ZC3H11B, BMP3, COL22A1, DLG5, PLCE1, EFEMP2, TYR, FAT3, TRIM29, COL2A1 and LOXL1. Replication in the 23andMe data set, where RD is self-reported by participants, firmly establishes six RD risk loci: FAT3, COL22A1, TYR, BMP3, ZC3H11B and PLCE1. Based on the genetic associations with eye traits described to date, the first two specifically impact risk of a RD, whereas the last four point to shared aetiologies with macular condition, myopia and glaucoma. Fine-mapping prioritized the lead common missense variant (TYR S192Y) as causal variant at the TYR locus and a small set of credible causal variants at the FAT3 locus. The larger study size presented here, enabled by resources linked to health records or self-report, provides novel insights into RD aetiology and underlying pathological pathways.

Analysis of putative cis-regulatory elements regulating blood pressure variation.

Hundreds of loci have been associated with blood pressure (BP) traits from many genome-wide association studies. We identified an enrichment of these loci in aorta and tibial artery expression quantitative trait loci in our previous work in ~100 000 Genetic Epidemiology Research on Aging study participants. In the present study, we sought to fine-map known loci and identify novel genes by determining putative regulatory regions for these and other tissues relevant to BP. We constructed maps of putative cis-regulatory elements (CREs) using publicly available open chromatin data for the heart, aorta and tibial arteries, and multiple kidney cell types. Variants within these regions may be evaluated quantitatively for their tissue- or cell-type-specific regulatory impact using deltaSVM functional scores, as described in our previous work. We aggregate variants within these putative CREs within 50 Kb of the start or end of 'expressed' genes in these tissues or cell types using public expression data and use deltaSVM scores as weights in the group-wise sequence kernel association test to identify candidates. We test for association with both BP traits and expression within these tissues or cell types of interest and identify the candidates MTHFR, C10orf32, CSK, NOV, ULK4, SDCCAG8, SCAMP5, RPP25, HDGFRP3, VPS37B and PPCDC. Additionally, we examined two known QT interval genes, SCN5A and NOS1AP, in the Atherosclerosis Risk in Communities Study, as a positive control, and observed the expected heart-specific effect. Thus, our method identifies variants and genes for further functional testing using tissue- or cell-type-specific putative regulatory information.

Gene- and tissue-level interactions in normal gastrointestinal development and Hirschsprung disease.

The development of the gut from endodermal tissue to an organ with multiple distinct structures and functions occurs over a prolonged time during embryonic days E10.5-E14.5 in the mouse. During this process, one major event is innervation of the gut by enteric neural crest cells (ENCCs) to establish the enteric nervous system (ENS). To understand the molecular processes underpinning gut and ENS development, we generated RNA-sequencing profiles from wild-type mouse guts at E10.5, E12.5, and E14.5 from both sexes. We also generated these profiles from homozygous Ret null embryos, a model for Hirschsprung disease (HSCR), in which the ENS is absent. These data reveal 4 major features: 1) between E10.5 and E14.5 the developmental genetic programs change from expression of major transcription factors and its modifiers to genes controlling tissue (epithelium, muscle, endothelium) specialization; 2) the major effect of Ret is not only on ENCC differentiation to enteric neurons but also on the enteric mesenchyme and epithelium; 3) a muscle genetic program exerts significant effects on ENS development; and 4) sex differences in gut development profiles are minor. The genetic programs identified, and their changes across development, suggest that both cell autonomous and nonautonomous factors, and interactions between the different developing gut tissues, are important for normal ENS development and its disorders.

Rare variants in fox-1 homolog A (RBFOX1) are associated with lower blood pressure.

Many large genome-wide association studies (GWAS) have identified common blood pressure (BP) variants. However, most of the identified BP variants do not overlap with the linkage evidence observed from family studies. We thus hypothesize that multiple rare variants contribute to the observed linkage evidence. We performed linkage analysis using 517 individuals in 130 European families from the Cleveland Family Study (CFS) who have been genotyped on the Illumina OmniExpress Exome array. The largest linkage peak was observed on chromosome 16p13 (MLOD = 2.81) for systolic blood pressure (SBP). Follow-up conditional linkage and association analyses in the linkage region identified multiple rare, coding variants in RBFOX1 associated with reduced SBP. In a 17-member CFS family, carriers of the missense variant rs149974858 are normotensive despite being obese (average BMI = 60 kg/m2). Gene-based association test of rare variants using SKAT-O showed significant association with SBP (p-value = 0.00403) and DBP (p-value = 0.0258) in the CFS participants and the association was replicated in large independent replication studies (N = 57,234, p-value = 0.013 for SBP, 0.0023 for PP). RBFOX1 is expressed in brain tissues, the atrial appendage and left ventricle in the heart, and in skeletal muscle tissues, organs/tissues which are potentially related to blood pressure. Our study showed that associations of rare variants could be efficiently detected using family information.

Trans-ethnic meta-analysis of genome-wide association studies for Hirschsprung disease.

Hirschsprung disease (HSCR) is the most common cause of neonatal intestinal obstruction. It is characterized by the absence of ganglia in the nerve plexuses of the lower gastrointestinal tract. So far, three common disease-susceptibility variants at the RET, SEMA3 and NRG1 loci have been detected through genome-wide association studies (GWAS) in Europeans and Asians to understand its genetic etiologies. Here we present a trans-ethnic meta-analysis of 507 HSCR cases and 1191 controls, combining all published GWAS results on HSCR to fine-map these loci and narrow down the putatively causal variants to 99% credible sets. We also demonstrate that the effects of RET and NRG1 are universal across European and Asian ancestries. In contrast, we detected a European-specific association of a low-frequency variant, rs80227144, in SEMA3 [odds ratio (OR) = 5.2, P = 4.7 × 10-10]. Conditional analyses on the lead SNPs revealed a secondary association signal, corresponding to an Asian-specific, low-frequency missense variant encoding RET p.Asp489Asn (rs9282834, conditional OR = 20.3, conditional P = 4.1 × 10-14). When in trans with the RET intron 1 enhancer risk allele, rs9282834 increases the risk of HSCR from 1.1 to 26.7. Overall, our study provides further insights into the genetic architecture of HSCR and has profound implications for future study designs.

The role of rare variants in systolic blood pressure: analysis of ExomeChip data in HyperGEN African Americans.

Cardiovascular diseases are among the most significant health problems in the United States today, with their major risk factor, hypertension, disproportionately affecting African Americans (AAs). Although GWAS have identified dozens of common variants associated with blood pressure (BP) and hypertension in European Americans, these variants collectively explain <2.5% of BP variance, and most of the genetic variants remain yet to be identified. Here, we report the results from rare-variant analysis of systolic BP using 94,595 rare and low-frequency variants (minor allele frequency, MAF, <5%) from the Illumina exome array genotyped in 2,045 HyperGEN AAs. In addition to single-variant analysis, 4 gene-level association tests were used for analysis: burden and family-based SKAT tests using MAF cutoffs of 1 and 5%. The gene-based methods often provided lower p values than the single-variant approach. Some consistency was observed across these 4 gene-based analysis options. While neither the gene-based analyses nor the single-variant analysis produced genome-wide significant results, the top signals, which had supporting evidence from multiple gene-based methods, were of borderline significance. Though additional molecular validations are required, 6 of the 16 most promising genes are biologically plausible with physiological connections to BP regulation.

Genome-wide meta-analyses of restless legs syndrome yield insights into genetic architecture, disease biology and risk prediction.

Restless legs syndrome (RLS) affects up to 10% of older adults. Their healthcare is impeded by delayed diagnosis and insufficient treatment. To advance disease prediction and find new entry points for therapy, we performed meta-analyses of genome-wide association studies in 116,647 individuals with RLS (cases) and 1,546,466 controls of European ancestry. The pooled analysis increased the number of risk loci eightfold to 164, including three on chromosome X. Sex-specific meta-analyses revealed largely overlapping genetic predispositions of the sexes (rg = 0.96). Locus annotation prioritized druggable genes such as glutamate receptors 1 and 4, and Mendelian randomization indicated RLS as a causal risk factor for diabetes. Machine learning approaches combining genetic and nongenetic information performed best in risk prediction (area under the curve (AUC) = 0.82-0.91). In summary, we identified targets for drug development and repurposing, prioritized potential causal relationships between RLS and relevant comorbidities and risk factors for follow-up and provided evidence that nonlinear interactions are likely relevant to RLS risk prediction.

Genome-wide meta-analysis of insomnia prioritizes genes associated with metabolic and psychiatric pathways.

Insomnia is a heritable, highly prevalent sleep disorder for which no sufficient treatment currently exists. Previous genome-wide association studies with up to 1.3 million subjects identified over 200 associated loci. This extreme polygenicity suggested that many more loci remain to be discovered. The current study almost doubled the sample size to 593,724 cases and 1,771,286 controls, thereby increasing statistical power, and identified 554 risk loci (including 364 novel loci). To capitalize on this large number of loci, we propose a novel strategy to prioritize genes using external biological resources and functional interactions between genes across risk loci. Of all 3,898 genes naively implicated from the risk loci, we prioritize 289 and find brain-tissue expression specificity and enrichment in specific gene sets of synaptic signaling functions and neuronal differentiation. We show that this novel gene prioritization strategy yields specific hypotheses on underlying mechanisms of insomnia that would have been missed by traditional approaches.

Biophysical and structural considerations for protein sequence evolution.

BACKGROUND: Protein sequence evolution is constrained by the biophysics of folding and function, causing interdependence between interacting sites in the sequence. However, current site-independent models of sequence evolutions do not take this into account. Recent attempts to integrate the influence of structure and biophysics into phylogenetic models via statistical/informational approaches have not resulted in expected improvements in model performance. This suggests that further innovations are needed for progress in this field. RESULTS: Here we develop a coarse-grained physics-based model of protein folding and binding function, and compare it to a popular informational model. We find that both models violate the assumption of the native sequence being close to a thermodynamic optimum, causing directional selection away from the native state. Sampling and simulation show that the physics-based model is more specific for fold-defining interactions that vary less among residue type. The informational model diffuses further in sequence space with fewer barriers and tends to provide less support for an invariant sites model, although amino acid substitutions are generally conservative. Both approaches produce sequences with natural features like dN/dS < 1 and gamma-distributed rates across sites. CONCLUSIONS: Simple coarse-grained models of protein folding can describe some natural features of evolving proteins but are currently not accurate enough to use in evolutionary inference. This is partly due to improper packing of the hydrophobic core. We suggest possible improvements on the representation of structure, folding energy, and binding function, as regards both native and non-native conformations, and describe a large number of possible applications for such a model.

Nuclear genome-wide associations with mitochondrial heteroplasmy.

The role of the nuclear genome in maintaining the stability of the mitochondrial genome (mtDNA) is incompletely known. mtDNA sequence variants can exist in a state of heteroplasmy, which denotes the coexistence of organellar genomes with different sequences. Heteroplasmic variants that impair mitochondrial capacity cause disease, and the state of heteroplasmy itself is deleterious. However, mitochondrial heteroplasmy may provide an intermediate state in the emergence of novel mitochondrial haplogroups. We used genome-wide genotyping data from 982,072 European ancestry individuals to evaluate variation in mitochondrial heteroplasmy and to identify the regions of the nuclear genome that affect it. Age, sex, and mitochondrial haplogroup were associated with the extent of heteroplasmy. GWAS identified 20 loci for heteroplasmy that exceeded genome-wide significance. This included a region overlapping mitochondrial transcription factor A (TFAM), which has multiple roles in mtDNA packaging, replication, and transcription. These results show that mitochondrial heteroplasmy has a heritable nuclear component.

Multiple, independent, common variants at RET, SEMA3 and NRG1 gut enhancers specify Hirschsprung disease risk in European ancestry subjects.

PURPOSE: Hirschsprung disease (HSCR) is a developmental disorder of the enteric nervous system (ENS) characterized by congenital aganglionosis arising from coding variants in ENS genes causing partial or total loss-of-function. Low-penetrance, common, noncoding variants at RET, SEMA3 and NRG1 loci are also associated with HSCR, with small-to-moderate loss of gene expression mediated through sequence variants in cis-regulatory elements (CRE) as another causal mechanism. Since these latter variants are common, many individuals carry multiple risk variants. However, the extent and combinatorial effects of all putative CRE variants within and across these loci on HSCR is unknown. METHODS: Using 583 HSCR subjects, one of the largest samples of European ancestry studied, and genotyping 56 tag variants, we evaluated association of all common variants overlapping putative gut CREs and fine-mapped causal variants at RET, SEMA3 and NRG1. RESULTS: We demonstrate that 28 and 8 tag variants, several of which are genetically independent, overlap putative-enhancers at the RET and SEMA3 loci, respectively, as well as two fine-mapped tag variants at the NRG1 locus, are significantly associated with HSCR. Importantly, disease risk increases with increasing numbers of risk alleles from multiple variants within and across these loci, varying >25-fold across individuals. CONCLUSION: This increasing allele number-dependent risk, we hypothesize, arises from HSCR-relevant ENS cells sensing the reduced gene expression at multiple ENS genes since their developmental effects are integrated through gene regulatory networks.

Genetic analyses identify widespread sex-differential participation bias.

Genetic association results are often interpreted with the assumption that study participation does not affect downstream analyses. Understanding the genetic basis of participation bias is challenging since it requires the genotypes of unseen individuals. Here we demonstrate that it is possible to estimate comparative biases by performing a genome-wide association study contrasting one subgroup versus another. For example, we showed that sex exhibits artifactual autosomal heritability in the presence of sex-differential participation bias. By performing a genome-wide association study of sex in approximately 3.3 million males and females, we identified over 158 autosomal loci spuriously associated with sex and highlighted complex traits underpinning differences in study participation between the sexes. For example, the body mass index-increasing allele at FTO was observed at higher frequency in males compared to females (odds ratio = 1.02, P = 4.4 × 10-36). Finally, we demonstrated how these biases can potentially lead to incorrect inferences in downstream analyses and propose a conceptual framework for addressing such biases. Our findings highlight a new challenge that genetic studies may face as sample sizes continue to grow.

Combined linkage and association analysis identifies rare and low frequency variants for blood pressure at 1q31.

High blood pressure (BP) is a major risk factor for cardiovascular disease (CVD) and is more prevalent in African Americans as compared to other US groups. Although large, population-based genome-wide association studies (GWAS) have identified over 300 common polymorphisms modulating inter-individual BP variation, largely in European ancestry subjects, most of them do not localize to regions previously identified through family-based linkage studies. This discrepancy has remained unexplained despite the statistical power differences between current GWAS and prior linkage studies. To address this issue, we performed genome-wide linkage analysis of BP traits in African-American families from the Family Blood Pressure Program (FBPP) and genotyped on the Illumina Human Exome BeadChip v1.1. We identified a genomic region on chromosome 1q31 with LOD score 3.8 for pulse pressure (PP), a region we previously implicated in DBP studies of European ancestry families. Although no reported GWAS variants map to this region, combined linkage and association analysis of PP identified 81 rare and low frequency exonic variants accounting for the linkage evidence. Replication analysis in eight independent African ancestry cohorts (N = 16,968) supports this specific association with PP (P = 0.0509). Additional association and network analyses identified multiple potential candidate genes in this region expressed in multiple tissues and with a strong biological support for a role in BP. In conclusion, multiple genes and rare variants on 1q31 contribute to PP variation. Beyond producing new insights into PP, we demonstrate how family-based linkage and association studies can implicate specific rare and low frequency variants for complex traits.

Contributions of rare coding variants in hypotension syndrome genes to population blood pressure variation.

Rare variants, in particular renal salt handling genes, contribute to monogenic forms of hypertension and hypotension syndromes with electrolyte abnormalities. A study by Ji et al (2008) demonstrated this effect for rare loss-of-function coding variants in SLC12A3 (NCCT), SLC12A1 (NKCC2), and KCNJ1 (ROMK) that led to reduction of ∼6 mm Hg for SBP and ∼3 mm Hg for DBP among carriers in 2492 European ancestry Framingham Heart Study (FHS) subjects. These findings support a potentially large role for these variants in interindividual variation in systolic and diastolic blood pressure (SBP, DBP) in the population. The present study focuses on replicating the analyses completed by Ji et al to identify effects of rare variants in the population-based Atherosclerosis Risk in Communities (ARIC) study.We attempted to replicate the findings by Ji et al by applying their criteria to identify putative loss-of-function variants with allele frequency <0.001 and complete conservation across a set of orthologs, to exome sequencing data from 7444 European ancestry participants of the ARIC study.Although we failed to replicate the previous findings when applying their methods to the ARIC study data, we observed a similar effect when we restricted analyses to the subset of variants they observed.These results simultaneously support the utility of exome sequencing data for studying extremely rare coding variants in hypertension and underscore the need for improved filtering methods for identifying functional variants in human sequences.

Genome-wide association study of Hirschsprung disease detects a novel low-frequency variant at the RET locus.

Hirschsprung disease (HSCR) is a congenital disorder with a population incidence of ~1/5000 live births, defined by an absence of enteric ganglia along variable lengths of the colon. HSCR genome-wide association studies (GWAS) have found common associated variants at RET, SEMA3, and NRG1, but they still fail to explain all of its heritability. To enhance gene discovery, we performed a GWAS of 170 cases identified from the Danish nationwide pathology registry with 4717 controls, based on 6.2 million variants imputed from the haplotype reference consortium panel. We found a novel low-frequency variant (rs144432435), which, when conditioning on the lead RET single-nucleotide polymorphism (SNP), was of genome-wide significance in the discovery analysis. This conditional association signal was replicated in a Swedish HSCR cohort with discovery plus replication meta-analysis conditional odds ratio of 6.6 (P = 7.7 × 10-10; 322 cases and 4893 controls). The conditional signal was, however, not replicated in two HSCR cohorts from USA and Finland, leading to the hypothesis that rs144432435 tags a rare haplotype present in Denmark and Sweden. Using the genome-wide complex trait analysis method, we estimated the SNP heritability of HSCR to be 88%, close to estimates based on classical family studies. Moreover, by using Lasso (least absolute shrinkage and selection operator) regression we were able to construct a genetic HSCR predictor with a area under the receiver operator characteristics curve of 76% in an independent validation set. In conclusion, we combined the largest collection of sporadic Hirschsprung cases to date (586 cases) to further elucidate HSCR's genetic architecture.

MicroRNAs in the miR-17 and miR-15 families are downregulated in chronic kidney disease with hypertension.

BACKGROUND: In older adults (aged 70-74 years), African-Americans have 4-fold higher risk of developing hypertension-attributed end-stage renal disease (ESRD) than European-Americans. A hypothesized mechanism linking hypertension and progressive chronic kidney disease (CKD) is the innate immune response and inflammation. Persons with CKD are also more susceptible to infection. Gene expression in peripheral blood can provide a view of the innate immune activation profile. We aimed to identify differentially expressed genes, microRNAs, and pathways in peripheral blood between cases with CKD and high blood pressure under hypertension treatment versus controls without CKD and with controlled blood pressure in African Americans. METHODS: Case and control pairs (N = 15x2) were selected from those without diabetes and were matched for age, sex, body mass index, APOL1 risk allele count, and hypertension medication use. High blood pressure under hypertension treatment was defined as hypertension medication use and systolic blood pressure (SBP) ≥ 145 mmHg. CKD was defined as estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73m2. Cases were selected from those with CKD and high blood pressure under hypertension treatment, and controls were selected from those without CKD (eGFR: 75-120 mL/min/1.73m2 and urine albumin-to-creatinine ratio < 30mg/g) and with blood pressure controlled by hypertension medication use (SBP < 135 mmHg and D(diastolic)BP < 90 mm Hg). We perform RNA sequencing of mRNA and microRNA and conducted differential expression and co-expression network analysis. RESULTS: Of 347 miRNAs included in the analysis, 14 were significantly associated with case status (Benjamini-Hochberg adjusted p-value [BH p] < 0.05). Of these, ten were downregulated in cases: three of each belong to the miR-17 and miR-15 families. In co-expression network analysis of miRNA, one module, which included 13 of the 14 significant miRNAs, had significant association with case status. Of the 14,488 genes and 41,739 transcripts included in the analysis, none had significant association with case status. Gene co-expression network analyses did not yield any significant associations for mRNA. CONCLUSION: We have identified 14 differentially expressed miRNAs in the peripheral blood of CKD cases with high blood pressure under hypertension treatment as compared to appropriate controls. Most of the significant miRNAs were downregulated and have critical roles in immune cell functions. Future studies are needed to replicate our findings and determine whether the downregulation of these miRNAs in immune cells may influence CKD progression or complications.

Genome-wide association analyses using electronic health records identify new loci influencing blood pressure variation.

Longitudinal electronic health records on 99,785 Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort individuals provided 1,342,814 systolic and diastolic blood pressure measurements for a genome-wide association study on long-term average systolic, diastolic, and pulse pressure. We identified 39 new loci among 75 genome-wide significant loci (P ≤ 5 × 10-8), with most replicating in the combined International Consortium for Blood Pressure (ICBP; n = 69,396) and UK Biobank (UKB; n = 152,081) studies. Combining GERA with ICBP yielded 36 additional new loci, with most replicating in UKB. Combining all three studies (n = 321,262) yielded 241 additional genome-wide significant loci, although no replication sample was available for these. All associated loci explained 2.9%, 2.5%, and 3.1% of variation in systolic, diastolic, and pulse pressure, respectively, in GERA non-Hispanic whites. Using multiple blood pressure measurements in GERA doubled the variance explained. A normalized risk score was associated with time to onset of hypertension (hazards ratio = 1.18, P = 8.2 × 10-45). Expression quantitative trait locus analysis of blood pressure loci showed enrichment in aorta and tibial artery.

Rare coding variants associated with blood pressure variation in 15 914 individuals of African ancestry.

OBJECTIVES: Hypertension is a major risk factor for all cardiovascular diseases, especially among African Americans. This study focuses on identifying specific blood pressure (BP) genes using 15 914 individuals of African ancestry from eight cohorts (Africa America Diabetes Mellitus, Atherosclerosis Risk in Communities Study, Coronary Artery Risk Development in young Adults, Genetics Network, Genetic Epidemiology Network of Arteriopathy, Howard University Family Study, Hypertension Genetic Epidemiology Network, and Loyola University Chicago Cohort) to further genetic findings in this population which has generally been underrepresented in BP studies. METHODS: We genotyped and performed various single variant and gene-based exome-wide analyses on 15 914 individuals on the Illumina HumanExome Beadchip v1.0 or v1.1 to test association with SBP and DBP long-term average residuals that were adjusted for age, age-squared, sex, and BMI. RESULTS: We identified rare variants affecting SBP and DBP in 10 genes: AFF1, GAPDHS, SLC28A3, COL6A1, CRYBA2, KRBA1, SEL1L3, YOD1, CCDC13, and QSOX1. Prior experimental evidence for six of these 10 candidate genes supports their involvement in cardiovascular mechanisms, corroborating their potential roles in BP regulation. CONCLUSION: Although our results require replication or validation due to their low numbers of carriers, and an ethnicity-specific genotyping array may be more informative, this study, which has identified several candidate genes in this population most susceptible to hypertension, presents one of the largest African-ancestry BP studies to date and the largest including analysis of rare variants.