Uganda Genome Resource Enables Insights into Population History and Genomic Discovery in Africa.

Genomic studies in African populations provide unique opportunities to understand disease etiology, human diversity, and population history. In the largest study of its kind, comprising genome-wide data from 6,400 individuals and whole-genome sequences from 1,978 individuals from rural Uganda, we find evidence of geographically correlated fine-scale population substructure. Historically, the ancestry of modern Ugandans was best represented by a mixture of ancient East African pastoralists. We demonstrate the value of the largest sequence panel from Africa to date as an imputation resource. Examining 34 cardiometabolic traits, we show systematic differences in trait heritability between European and African populations, probably reflecting the differential impact of genes and environment. In a multi-trait pan-African GWAS of up to 14,126 individuals, we identify novel loci associated with anthropometric, hematological, lipid, and glycemic traits. We find that several functionally important signals are driven by Africa-specific variants, highlighting the value of studying diverse populations across the region.

Whole-exome sequencing study identifies four novel gene loci associated with diabetic kidney disease.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease and diabetes. Our two-stage WES study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort and Atherosclerosis Risk in Communities studies (stage 1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine participants (stage 2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single-variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds [95% confidence interval (CI): 33.6, 1105] of DKD compared with noncarriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% CI: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

African ancestry-derived APOL1 risk genotypes show proximal epigenetic associations.

Apolipoprotein L1 (APOL1) coding variants, termed G1 and G2, are established genetic risk factors for a growing spectrum of diseases, including kidney disease, in individuals of African ancestry. Evidence suggests that the risk variants, which show a recessive mode of inheritance, lead to toxic gain-of-function changes of the APOL1 protein. Disease occurrence and presentation vary, likely due to modifiers or second hits. To understand the role of the epigenetic landscape in relation to APOL1 risk variants, we performed methylation quantitative trait locus (meQTL) analysis to identify differentially methylated CpGs influenced by APOL1 risk variants in 611 African American individuals. We identified five CpGs that were significantly associated with APOL1 risk alleles in discovery and replication studies, and one CpG-APOL1 association was independent of other genomic variants. Our study highlights proximal DNA methylation alterations that may help explain the variable disease risk and clinical manifestation of APOL1 variants.

Collectrin (Tmem27) deficiency in proximal tubules causes hypertension in mice and a TMEM27 variant associates with blood pressure in males in a Latino cohort.

Collectrin (Tmem27), an angiotensin-converting enzyme 2 homologue, is a chaperone of amino acid transporters in the kidney and endothelium. Global collectrin knockout (KO) mice have hypertension, endothelial dysfunction, exaggerated salt sensitivity, and diminished renal blood flow. This phenotype is associated with altered nitric oxide and superoxide balance and increased proximal tubule (PT) Na+/H+ exchanger isoform 3 (NHE3) expression. Collectrin is located on the X chromosome where genome-wide association population studies have largely been excluded. In the present study, we generated PT-specific collectrin KO (PT KO) mice to determine the precise contribution of PT collectrin in blood pressure homeostasis. We also examined the association of human TMEM27 single-nucleotide polymorphisms with blood pressure traits in 11,926 Hispanic Community Health Study/Study of Latinos (HCHS/SOL) Hispanic/Latino participants. PT KO mice exhibited hypertension, and this was associated with increased baseline NHE3 expression and diminished lithium excretion. However, PT KO mice did not display exaggerated salt sensitivity or a reduction in renal blood flow compared with control mice. Furthermore, PT KO mice exhibited enhanced endothelium-mediated dilation, suggesting a compensatory response to systemic hypertension induced by deficiency of collectrin in the PT. In HCHS/SOL participants, we observed sex-specific single-nucleotide polymorphism associations with diastolic blood pressure. In conclusion, loss of collectrin in the PT is sufficient to induce hypertension, at least in part, through activation of NHE3. Importantly, our model supports the notion that altered renal blood flow may be a determining factor for salt sensitivity. Further studies are needed to investigate the role of the TMEM27 locus on blood pressure and salt sensitivity in humans.NEW & NOTEWORTHY The findings of our study are significant in several ways: 1) loss of an amino acid chaperone in the proximal tubule is sufficient to cause hypertension, 2) the results in global and proximal tubule-specific collectrin knockout mice support the notion that vascular dysfunction is required for salt sensitivity or that impaired renal tubule function causes hypertension but is not sufficient to cause salt sensitivity, and 3) our study is the first to implicate a role of collectrin in human hypertension.

Discovery and fine-mapping of kidney function loci in first genome-wide association study in Africans.

Genome-wide association studies (GWAS) of kidney function have uncovered hundreds of loci, primarily in populations of European ancestry. We have undertaken the first continental African GWAS of estimated glomerular filtration rate (eGFR), a measure of kidney function used to define chronic kidney disease (CKD). We conducted GWAS of eGFR in 3288 East Africans from the Uganda General Population Cohort (GPC) and replicated in 8224 African Americans from the Women's Health Initiative. Loci attaining genome-wide significant evidence for association (P < 5 × 10-8) were followed up with Bayesian fine-mapping to localize potential causal variants. The predictive power of a genetic risk score (GRS) constructed from previously reported trans-ancestry eGFR lead single nucleotide polymorphism (SNPs) was evaluated in the Uganda GPC. We identified and validated two eGFR loci. At the glycine amidinotransferase (GATM) locus, the association signal (lead SNP rs2433603, P = 1.0 × 10-8) in the Uganda GPC GWAS was distinct from previously reported signals at this locus. At the haemoglobin beta (HBB) locus, the association signal (lead SNP rs141845179, P = 3.0 × 10-8) has been previously reported. The lead SNP at the HBB locus accounted for 88% of the posterior probability of causality after fine-mapping, but did not colocalise with kidney expression quantitative trait loci. The trans-ancestry GRS of eGFR was not significantly predictive into the Ugandan population. In the first GWAS of eGFR in continental Africa, we validated two previously reported loci at GATM and HBB. At the GATM locus, the association signal was distinct from that previously reported. These results demonstrate the value of performing GWAS in continental Africans, providing a rich genomic resource to larger consortia for further discovery and fine-mapping. The study emphasizes that additional large-scale efforts in Africa are warranted to gain further insight into the genetic architecture of CKD.

Clonal Hematopoiesis of Indeterminate Potential and Kidney Function Decline in the General Population.

RATIONALE & OBJECTIVE: Clonal hematopoiesis of indeterminate potential (CHIP), defined by the age-related ontogenesis of expanded leukemogenic variants indicative of a genetically distinct clonal leukocyte population, is associated with risk of hematologic malignancy and cardiovascular disease. In experimental models, recapitulation of CHIP promotes kidney interstitial fibrosis with direct tissue infiltration of donor macrophages. We tested the hypothesis that CHIP is associated with kidney function decline in the general population. STUDY DESIGN: Cohort study. SETTING & PARTICIPANTS: 12,004 individuals from 3 community-based cohorts in the TOPMed Consortium. EXPOSURE: CHIP status from whole-genome sequences obtained from DNA extracted from peripheral blood. OUTCOME: Risk of 30% decline in estimated glomerular filtration rate (eGFR) and percent eGFR decline per year during the follow-up period. ANALYTICAL APPROACH: Cox proportional hazards models for 30% eGFR decline end point and generalized estimating equations for annualized relative change in eGFR with meta-analysis. Study-specific estimates were combined using fixed-effect meta-analysis. RESULTS: The median baseline eGFR was 84mL/min/1.73m2. The prevalence of CHIP was 6.6%, 9.0%, and 12.2% in persons aged 50-60, 60-70, and>70 years, respectively. Over a median follow-up period of 8 years, for the 30% eGFR outcome 205 events occurred among 1,002 CHIP carriers (2.1 events per 100 person-years) and 2,041 events in persons without CHIP (1.7 events per 100 person-years). In meta-analysis, CHIP was associated with greater risk of a 30% eGFR decline (17% [95% CI, 1%-36%] higher; P=0.04). Differences were not observed between those with baseline eGFR above or below 60mL/min/1.73m2, of age above or below 60 years, or with or without diabetes. LIMITATIONS: Small number of participants with moderate-to-advanced kidney disease and restricted set of CHIP driver variants. CONCLUSIONS: We report an association between CHIP and eGFR decline in 3 general population cohorts without known kidney disease. Further studies are needed to investigate this novel condition and its potential impact among individuals with overt kidney disease.

Hypertension burden, treatment, and control among people with HIV at a clinical care center in the Southeastern US, 2014-2019.

Hypertension management outcomes in people with HIV (PWH) are not well characterized, despite high hypertension burden. We assessed hypertension prevalence, incidence, treatment, and outcomes among patients with HIV at a clinical center in the southeastern US, from 2014 to 2019. To identify characteristics associated with treatment and outcomes, we estimated adjusted risk ratios (aRR) and 95% confidence intervals (CI). Among 2274 patients, 72% were cisgender men, 56% non-Hispanic Black, median age 47 years, 48% MSM, 12% had CD4 cell count <200 cells/µl, 72% HIV RNA level <400 copies/mL and 39% prevalent hypertension. Hypertension incidence rate was 6.3/100 person-years (95% CI, 5.6-7.0). Among incident hypertension cases (n = 275), 16% (95% CI, 11-20) initiated an antihypertensive within one year. Compared to non-Hispanic white patients, Hispanic (aRR, 6.68; 95% CI, 1.50-29.74) and non-Hispanic Black patients (aRR, 2.18; 95% CI, 0.91-5.24) were more likely to initiate an antihypertensive. Among patients initiating an antihypertensive (n = 178), 63% (95% CI 56-70) experienced blood pressure control within one year. Patients with HIV experienced a high burden of hypertension with notable delays in antihypertensive initiation, as well as gaps in achieving blood pressure control, highlighting opportunities for interventions designed to minimize delays in controlling hypertension in this vulnerable population.

Genome-Wide Admixture Mapping of Estimated Glomerular Filtration Rate and Chronic Kidney Disease Identifies European and African Ancestry-of-Origin Loci in Hispanic and Latino Individuals in the United States.

BACKGROUND: Admixture mapping is a powerful approach for gene mapping of complex traits that leverages the diverse genetic ancestry in populations with recent admixture, such as Hispanic or Latino individuals in the United States. These individuals have an increased risk of CKD. METHODS: We performed genome-wide admixture mapping for both CKD and eGFR in a sample of 12,601 participants from the Hispanic Community Health Study/Study of Latinos, with validation in a sample of 8191 Black participants from the Women's Health Initiative (WHI). We also compared the findings with those from a conventional genome-wide association study. RESULTS: Three novel ancestry-of-origin loci were identified on chromosomes 2, 14, and 15 for CKD and eGFR. The chromosome 2 locus comprises two European ancestry regions encompassing the FSHR and NRXN1 genes, with European ancestry at this locus associated with increased CKD risk. The chromosome 14 locus, found within the DLK1-DIO3 imprinted domain, was associated with lower eGFR and driven by European ancestry. The eGFR-associated locus on chromosome 15 included intronic variants of RYR3 and was within an African-specific genomic region associated with higher eGFR. The genome-wide association study failed to identify significant associations in these regions. We validated the chromosome 14 and 15 loci associated with eGFR in the WHI Black participants. CONCLUSIONS: This study provides evidence of shared ancestry-specific genomic regions influencing eGFR in Hispanic or Latino individuals and Black individuals and illustrates the potential for leveraging genetic ancestry in recently admixed populations for the discovery of novel candidate loci for kidney phenotypes.

Rare coding variants in RCN3 are associated with blood pressure.

BACKGROUND: While large genome-wide association studies have identified nearly one thousand loci associated with variation in blood pressure, rare variant identification is still a challenge. In family-based cohorts, genome-wide linkage scans have been successful in identifying rare genetic variants for blood pressure. This study aims to identify low frequency and rare genetic variants within previously reported linkage regions on chromosomes 1 and 19 in African American families from the Trans-Omics for Precision Medicine (TOPMed) program. Genetic association analyses weighted by linkage evidence were completed with whole genome sequencing data within and across TOPMed ancestral groups consisting of 60,388 individuals of European, African, East Asian, Hispanic, and Samoan ancestries. RESULTS: Associations of low frequency and rare variants in RCN3 and multiple other genes were observed for blood pressure traits in TOPMed samples. The association of low frequency and rare coding variants in RCN3 was further replicated in UK Biobank samples (N = 403,522), and reached genome-wide significance for diastolic blood pressure (p = 2.01 × 10- 7). CONCLUSIONS: Low frequency and rare variants in RCN3 contributes blood pressure variation. This study demonstrates that focusing association analyses in linkage regions greatly reduces multiple-testing burden and improves power to identify novel rare variants associated with blood pressure traits.

Kidney omics in hypertension: from statistical associations to biological mechanisms and clinical applications.

Hypertension is a major cardiovascular disease risk factor and contributor to premature death globally. Family-based investigations confirmed a significant heritable component of blood pressure (BP), whereas genome-wide association studies revealed >1000 common and rare genetic variants associated with BP and/or hypertension. The kidney is not only an organ of key relevance to BP regulation and the development of hypertension, but it also acts as the tissue mediator of genetic predisposition to hypertension. The identity of kidney genes, pathways, and related mechanisms underlying the genetic associations with BP has started to emerge through integration of genomics with kidney transcriptomics, epigenomics, and other omics as well as through applications of causal inference, such as Mendelian randomization. Single-cell methods further enabled mapping of BP-associated kidney genes to cell types, and in conjunction with other omics, started to illuminate the biological mechanisms underpinning associations of BP-associated genetic variants and kidney genes. Polygenic risk scores derived from genome-wide association studies and refined on kidney omics hold the promise of enhanced diagnostic prediction, whereas kidney omics-informed drug discovery is likely to contribute new therapeutic opportunities for hypertension and hypertension-mediated kidney damage.

Metabolome-wide association study of estimated glomerular filtration rates in Hispanics.

Circulating metabolites are by-products of endogenous metabolism or exogenous sources and may inform disease states. Our study aimed to identify the source of variability in the association of metabolites with estimated glomerular filtration rate (eGFR) in Hispanics/Latinos with low chronic kidney disease prevalence by testing the association of 640 metabolites in 3,906 participants of the Hispanic Community Health Study/Study of Latinos. Metabolites were quantified in fasting serum through non-targeted mass spectrometry analysis. eGFR was regressed on inverse normally transformed metabolites in models accounting for study design and covariates. To identify the source of variation on eGFR associations, we tested the interaction of metabolites with lifestyle and clinical risk factors, and results were integrated with genotypes to identify metabolite genetic regulation. The mean age was 46 years, 43% were men, 22% were current smokers, 47% had a Caribbean Hispanic background, 19% had diabetes and the mean cohort eGFR was 96.4 ml/min/1.73 m2. We identified 404 eGFR-metabolite associations (False Discovery Rate under 0.05). Of these, 69 were previously reported, and 79 were novel associations with eGFR replicated in one or more published studies. There were significant interactions with lifestyle and clinical risk factors, with larger differences in eGFR-metabolite associations within strata of age, urine albumin to creatinine ratio, diabetes and Hispanic/Latino background. Several newly identified metabolites were genetically regulated, and variants were located at genomic regions previously associated with eGFR. Thus, our results suggest complex mechanisms contribute to the association of eGFR with metabolites and provide new insights into these associations.

Lipidomic profiling in the Strong Heart Study identified American Indians at risk of chronic kidney disease.

Dyslipidemia associates with and usually precedes the onset of chronic kidney disease (CKD), but a comprehensive assessment of molecular lipid species associated with risk of CKD is lacking. Here, we sought to identify fasting plasma lipids associated with risk of CKD among American Indians in the Strong Heart Family Study, a large-scale community-dwelling of individuals, followed by replication in Mexican Americans from the San Antonio Family Heart Study and Caucasians from the Australian Diabetes, Obesity and Lifestyle Study. We also performed repeated measurement analysis to examine the temporal relationship between the change in the lipidome and change in kidney function between baseline and follow-up of about five years apart. Network analysis was conducted to identify differential lipid classes associated with risk of CKD. In the discovery cohort, we found that higher baseline level of multiple lipid species, including glycerophospholipids, glycerolipids and sphingolipids, was significantly associated with increased risk of CKD, independent of age, sex, body mass index, diabetes and hypertension. Many lipid species were replicated in at least one external cohort at the individual lipid species and/or the class level. Longitudinal change in the plasma lipidome was significantly associated with change in the estimated glomerular filtration rate after adjusting for covariates, baseline lipids and the baseline rate. Network analysis identified distinct lipidomic signatures differentiating high from low-risk groups. Thus, our results demonstrated that disturbed lipid metabolism precedes the onset of CKD. These findings shed light on the mechanisms linking dyslipidemia to CKD and provide potential novel biomarkers for identifying individuals with early impaired kidney function at preclinical stages.

Analyses of biomarker traits in diverse UK biobank participants identify associations missed by European-centric analysis strategies.

Despite the dramatic underrepresentation of non-European populations in human genetics studies, researchers continue to exclude participants of non-European ancestry, as well as variants rare in European populations, even when these data are available. This practice perpetuates existing research disparities and can lead to important and large effect size associations being missed. Here, we conducted genome-wide association studies (GWAS) of 31 serum and urine biomarker quantitative traits in African (n = 9354), East Asian (n = 2559), and South Asian (n = 9823) ancestry UK Biobank (UKBB) participants. We adjusted for all known GWAS catalog variants for each trait, as well as novel signals identified in a recent European ancestry-focused analysis of UKBB participants. We identify 7 novel signals in African ancestry and 2 novel signals in South Asian ancestry participants (p < 1.61E-10). Many of these signals are highly plausible, including a cis pQTL for the gene encoding gamma-glutamyl transferase and PIEZO1 and G6PD variants with impacts on HbA1c through likely erythrocytic mechanisms. This work illustrates the importance of using the genetic data we already have in diverse populations, with novel discoveries possible in even modest sample sizes.

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.

Plasma Protein Profile of Carotid Artery Atherosclerosis and Atherosclerotic Outcomes: Meta-Analyses and Mendelian Randomization Analyses.

[Figure: see text].

Genetic and stress influences on the prevalence of hypertension among hispanics/latinos in the hispanic community health study/study of latinos (HCHS/SOL).

PURPOSE: The current study examined the effects of chronic stress and a genetic risk score on the presence of hypertension and elevated systolic blood pressure and diastolic blood pressure among Hispanics/Latinos in the target population of Hispanic Community Health Study/Study of Latinos. MATERIALS AND METHODS: Of the participants (N = 11,623) assessed during two clinic visits (Visit 1 2008-2013 & Visit 2 2014-2018), we analysed data from 7,429 adults (50.4% female), aged 18-74, who were genotyped and responded to chronic stress questionnaires. We calculated an unweighted genetic risk score using blood pressure increasing single nucleotide polymorphisms (SNPs) found to be generalisable to Hispanics/Latinos (10 SNPs). Linear and logistic regression models were used to estimate associations between chronic stress and genetic risk score and their interaction, with prevalent Visit 2 SBP or DBP, and hypertension, respectively. Models accounted for sampling weights, stratification, and cluster design. RESULTS: Chronic stress (adjusted OR = 1.18, 95%CI:1.15,1.22) and hypertension genetic risk score (adjusted OR = 1.04, 95%CI:1.01,1.07) were significantly associated with prevalent hypertension, but there was no significant interaction between the chronic stress and genetic risk score on hypertension (p = .49). genetic risk score (b = .32, 95%CI:.08, .55, R2 = .02) and chronic stress (b = .45, 95%CI:.19, .72, R2 = .11) were related to DBP, with no significant interaction (p = .62). Genetic risk score (b = .42, 95%CI:.08, .76, R2 = .01) and chronic stress (b = .80, 95%CI:.34,1.26, R2 = .11) were also related to SBP, with no significant interaction (p = .51). CONCLUSION: Results demonstrate the utility of a genetic risk score for blood pressure and are consistent with literature suggesting chronic stress has a strong, direct association with elevated blood pressure among U.S. Hispanics/Latinos.

A System for Phenotype Harmonization in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) Program.

Genotype-phenotype association studies often combine phenotype data from multiple studies to increase statistical power. Harmonization of the data usually requires substantial effort due to heterogeneity in phenotype definitions, study design, data collection procedures, and data-set organization. Here we describe a centralized system for phenotype harmonization that includes input from phenotype domain and study experts, quality control, documentation, reproducible results, and data-sharing mechanisms. This system was developed for the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program, which is generating genomic and other -omics data for more than 80 studies with extensive phenotype data. To date, 63 phenotypes have been harmonized across thousands of participants (recruited in 1948-2012) from up to 17 studies per phenotype. Here we discuss challenges in this undertaking and how they were addressed. The harmonized phenotype data and associated documentation have been submitted to National Institutes of Health data repositories for controlled access by the scientific community. We also provide materials to facilitate future harmonization efforts by the community, which include 1) the software code used to generate the 63 harmonized phenotypes, enabling others to reproduce, modify, or extend these harmonizations to additional studies, and 2) the results of labeling thousands of phenotype variables with controlled vocabulary terms.

A multi-ancestry genome-wide study incorporating gene-smoking interactions identifies multiple new loci for pulse pressure and mean arterial pressure.

Elevated blood pressure (BP), a leading cause of global morbidity and mortality, is influenced by both genetic and lifestyle factors. Cigarette smoking is one such lifestyle factor. Across five ancestries, we performed a genome-wide gene-smoking interaction study of mean arterial pressure (MAP) and pulse pressure (PP) in 129 913 individuals in stage 1 and follow-up analysis in 480 178 additional individuals in stage 2. We report here 136 loci significantly associated with MAP and/or PP. Of these, 61 were previously published through main-effect analysis of BP traits, 37 were recently reported by us for systolic BP and/or diastolic BP through gene-smoking interaction analysis and 38 were newly identified (P < 5 × 10-8, false discovery rate < 0.05). We also identified nine new signals near known loci. Of the 136 loci, 8 showed significant interaction with smoking status. They include CSMD1 previously reported for insulin resistance and BP in the spontaneously hypertensive rats. Many of the 38 new loci show biologic plausibility for a role in BP regulation. SLC26A7 encodes a chloride/bicarbonate exchanger expressed in the renal outer medullary collecting duct. AVPR1A is widely expressed, including in vascular smooth muscle cells, kidney, myocardium and brain. FHAD1 is a long non-coding RNA overexpressed in heart failure. TMEM51 was associated with contractile function in cardiomyocytes. CASP9 plays a central role in cardiomyocyte apoptosis. Identified only in African ancestry were 30 novel loci. Our findings highlight the value of multi-ancestry investigations, particularly in studies of interaction with lifestyle factors, where genomic and lifestyle differences may contribute to novel findings.

A Large-Scale Multi-ancestry Genome-wide Study Accounting for Smoking Behavior Identifies Multiple Significant Loci for Blood Pressure.

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).