Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes.

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

Genome-wide association analysis of composite sleep health scores in 413,904 individuals.

Recent genome-wide association studies (GWASs) of several individual sleep traits have identified hundreds of genetic loci, suggesting diverse mechanisms. Moreover, sleep traits are moderately correlated, and together may provide a more complete picture of sleep health, while also illuminating distinct domains. Here we construct novel sleep health scores (SHSs) incorporating five core self-report measures: sleep duration, insomnia symptoms, chronotype, snoring, and daytime sleepiness, using additive (SHS-ADD) and five principal components-based (SHS-PCs) approaches. GWASs of these six SHSs identify 28 significant novel loci adjusting for multiple testing on six traits (p<8.3e-9), along with 341 previously reported loci (p<5e-08). The heritability of the first three SHS-PCs equals or exceeds that of SHS-ADD (SNP-h2=0.094), while revealing sleep-domain-specific genetic discoveries. Significant loci enrich in multiple brain tissues and in metabolic and neuronal pathways. Post GWAS analyses uncover novel genetic mechanisms underlying sleep health and reveal connections to behavioral, psychological, and cardiometabolic traits.

Heritability and genetic correlations for sleep apnea, insomnia, and hypersomnia in a large clinical biobank.

RATIONALE: Comorbid insomnia and sleep apnea is reported to have worse outcomes than either condition alone. The local genetic correlations of these disorders are unknown. OBJECTIVES: To identify local genomic regions with heritability for clinically diagnosed sleep apnea and insomnia, and to identify local genetic correlations between these disorders and/or hypersomnia. METHODS: Fifty thousand two hundred seventeen patients of European ancestry were examined. Global and local heritability and genetic correlations for independent regions were calculated, adjusting for obesity and other covariates. RESULTS: Sleep apnea and insomnia were significantly globally heritable and had 118 and 168 genetic regions with local heritability p-values <.05, respectively. One region had a significant genetic correlation for sleep apnea and hypersomnia (p-value = 9.85 × 10-4). CONCLUSIONS: Clinically diagnosed sleep apnea and insomnia have minimal shared genetic architecture, supporting genetically distinct comorbid insomnia and sleep apnea components. However, additional correlated regions may be identified with additional sample size and methodological improvements.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Characterization of sleep apnea among a sample of adults from Samoa.

Sleep apnea is a public health concern around the world, but little research has been dedicated to examining this issue in low- and middle-income countries, including Samoa. Using data collected through the Soifua Manuia ("Good Health") study, which aimed to investigate the impact of the body mass index (BMI)-associated genetic variant rs373863828 in CREB3 Regulatory Factor ( CREBRF ) on metabolic traits in Samoan adults, we examined the sample prevalence and characteristics of sleep apnea using data collected with a validated home sleep apnea device (WatchPAT, Itamar). A total of 330 participants (sampled to overrepresent the obesity-risk allele of interest) had sleep data available. Participants (53.3% female) had a mean (SD) age of 52.0 (9.9) years and BMI of 35.5 (7.5) kg/m 2 and 36.3% of the sample had type 2 diabetes. Based on the 3% and 4% apnea hypopnea indices (AHI) and the 4% oxygen desaturation index (ODI), descriptive analyses revealed that many participants had potentially actionable sleep apnea defined as >5 events/hr (87.9%, 68.5%, and 71.2%, respectively) or clinically actionable sleep apnea defined as ≥15 events/hr (54.9%, 31.5%, and 34.5%, respectively). Sleep apnea was more severe in men; for example, clinically actionable sleep apnea (≥15) based on the AHI 3% definition was observed in 61.7% of men and 48.9% of women. Correction for non-representational sampling related to the CREBRF obesity-risk allele resulted in only slightly lower estimates. Across the AHI 3%, AHI 4%, and ODI 4%, multiple linear regression revealed associations between a greater number of events/hr and higher age, male sex, higher body mass index, higher abdominal-hip circumference ratio, and geographic region of residence. Our study identified a much higher frequency of sleep apnea in Samoa compared with published data from other studies, but similar predictors. Continued research addressing generalizability of these findings, as well as a specific focus on diagnosis and affordable and equitable access to treatment, is needed to alleviate the burden of sleep apnea in Samoa and around the world.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

The public health impact of poor sleep on severe COVID-19, influenza and upper respiratory infections.

BACKGROUND: Poor sleep is associated with an increased risk of infections and all-cause mortality but the causal direction between poor sleep and respiratory infections has remained unclear. We examined if poor sleep contributes as a causal risk factor to respiratory infections. METHODS: We used data on insomnia, influenza and upper respiratory infections (URIs) from primary care and hospital records in the UK Biobank (N ≈ 231,000) and FinnGen (N ≈ 392,000). We computed logistic regression to assess association between poor sleep and infections, disease free survival hazard ratios, and performed Mendelian randomization analyses to assess causality. FINDINGS: Utilizing 23 years of registry data and follow-up, we discovered that insomnia diagnosis associated with increased risk for infections (FinnGen influenza Cox's proportional hazard (CPH) HR = 4.34 [3.90, 4.83], P = 4.16 × 10-159, UK Biobank influenza CPH HR = 1.54 [1.37, 1.73], P = 2.49 × 10-13). Mendelian randomization indicated that insomnia causally predisposed to influenza (inverse-variance weighted (IVW) OR = 1.65, P = 5.86 × 10-7), URI (IVW OR = 1.94, P = 8.14 × 10-31), COVID-19 infection (IVW OR = 1.08, P = 0.037) and risk of hospitalization from COVID-19 (IVW OR = 1.47, P = 4.96 × 10-5). INTERPRETATION: Our findings indicate that chronic poor sleep is a causal risk factor for contracting respiratory infections, and in addition contributes to the severity of respiratory infections. These findings highlight the role of sleep in maintaining sufficient immune response against pathogens. FUNDING: Instrumentarium Science Foundation, Academy of Finland, Signe and Ane Gyllenberg Foundation, National Institutes of Health.

The genetic determinants of recurrent somatic mutations in 43,693 blood genomes.

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

Aberrant activation of TCL1A promotes stem cell expansion in clonal haematopoiesis.

Mutations in a diverse set of driver genes increase the fitness of haematopoietic stem cells (HSCs), leading to clonal haematopoiesis1. These lesions are precursors for blood cancers2-6, but the basis of their fitness advantage remains largely unknown, partly owing to a paucity of large cohorts in which the clonal expansion rate has been assessed by longitudinal sampling. Here, to circumvent this limitation, we developed a method to infer the expansion rate from data from a single time point. We applied this method to 5,071 people with clonal haematopoiesis. A genome-wide association study revealed that a common inherited polymorphism in the TCL1A promoter was associated with a slower expansion rate in clonal haematopoiesis overall, but the effect varied by driver gene. Those carrying this protective allele exhibited markedly reduced growth rates or prevalence of clones with driver mutations in TET2, ASXL1, SF3B1 and SRSF2, but this effect was not seen in clones with driver mutations in DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 led to the expression of TCL1A protein and the expansion of HSCs in vitro. The protective allele restricted TCL1A expression and expansion of mutant HSCs, as did experimental knockdown of TCL1A expression. Forced expression of TCL1A promoted the expansion of human HSCs in vitro and mouse HSCs in vivo. Our results indicate that the fitness advantage of several commonly mutated driver genes in clonal haematopoiesis may be mediated by TCL1A activation.

Impact of Pulmonary and Sleep Disorders on COVID-19 Infection Severity in a Large Clinical Biobank.

Rationale: Multiple pulmonary, sleep, and other disorders are associated with the severity of Covid-19 infections but may or may not directly affect the etiology of acute Covid-19 infection. Identifying the relative importance of concurrent risk factors may prioritize respiratory disease outbreaks research. Objectives: To identify associations of common preexisting pulmonary and sleep disease on acute Covid-19 infection severity, investigate the relative contributions of each disease and selected risk factors, identify sex-specific effects, and examine whether additional electronic health record (EHR) information would affect these associations. Methods: 45 pulmonary and 6 sleep diseases were examined in 37,020 patients with Covid-19. We analyzed three outcomes: death; a composite measure of mechanical ventilation and/or ICU admission; and inpatient admission. The relative contribution of pre-infection covariates including other diseases, laboratory tests, clinical procedures, and clinical note terms was calculated using LASSO. Each pulmonary/sleep disease model was then further adjusted for covariates. Measurements and main results: 37 pulmonary/sleep diseases were associated with at least one outcome at Bonferroni significance, 6 of which had increased relative risk in LASSO analyses. Multiple prospectively collected non-pulmonary/sleep diseases, EHR terms and laboratory results attenuated the associations between preexisting disease and Covid-19 infection severity. Adjustment for counts of prior "blood urea nitrogen" phrases in clinical notes attenuated the odds ratio point estimates of 12 pulmonary disease associations with death in women by ≥1. Conclusions: Pulmonary diseases are commonly associated with Covid-19 infection severity. Associations are partially attenuated by prospectively-collected EHR data, which may aid in risk stratification and physiological studies.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Obstructive Sleep Apnea and Longitudinal Changes in Interstitial Lung Imaging and Lung Function: The MESA Study.

Rationale: Obstructive sleep apnea (OSA) has been hypothesized to be a risk factor in interstitial lung disease (ILD) and is associated with radiological markers that may represent the earlier stages of ILD. Prior studies have been limited by their cross-sectional design and potential confounding by body habitus. Objectives: To test the hypothesis that OSA severity is associated with more high-attenuation areas (HAAs) on computed tomography and worse lung function over time among older community-dwelling adults. Methods: We used data from participants in the MESA (Multi-Ethnic Study of Atherosclerosis) who had apnea-hypopnea index (AHI) measured from polysomnography (2010-2013), high attenuation areas (HAAs, -600 to -250 Hounsfield units, n = 784), assessments from exams 5 (2010-2012) and 6 (2016-2018) full-lung computed tomography scans, and spirometry assessments (n = 677). Linear mixed-effects models with random intercept were used to examine associations of OSA severity (i.e., AHI and hypoxic burden) with changes in HAAs, total lung volumes, and forced vital capacity (FVC) between exams 5 and 6. Potential confounders were adjusted for in the model, including age, sex, smoking history, height, and weight. Results: Among those with a higher AHI there were more men and a higher body mass index. Participants with AHI ⩾ 15 events/h and in the highest hypoxic burden quartile each had increases in HAAs of 11.30% (95% confidence interval [CI], 3.74-19.35%) and 9.85% (95% CI, 1.40-19.01%) per 10 years, respectively. There was a more rapid decline in total lung volumes imaged and FVC among those with AHI ⩾ 15 events/h of 220.2 ml (95% CI, 47.8-392.5 ml) and 3.63% (95% CI, 0.43-6.83%) per 10 years, respectively. Conclusions: A greater burden of hypoxia related to obstructive events during sleep was associated with increased lung densities over time and a more rapid decline in lung volumes regardless of body habitus. Our findings suggest OSA may be a contributing factor in the early stages of ILD.

An integrated multi-omics analysis of sleep-disordered breathing traits implicates P2XR4 purinergic signaling.

Sleep Disordered Breathing (SDB) is a common disease associated with increased risk for cardiometabolic, cardiovascular, and cognitive diseases. How SDB affects the molecular environment is still poorly understood. We study the association of three SDB measures with gene expression measured using RNA-seq in multiple blood tissues from the Multi-Ethnic Study of Atherosclerosis. We develop genetic instrumental variables for the associated transcripts as polygenic risk scores (tPRS), then generalize and validate the tPRS in the Women's Health Initiative. We measure the associations of the validated tPRS with SDB and serum metabolites in Hispanic Community Health Study/Study of Latinos. Here we find differential gene expression by blood cell type in relation to SDB traits and link P2XR4 expression to average oxyhemoglobin saturation during sleep and butyrylcarnitine (C4) levels. These findings can be used to develop interventions to alleviate the effect of SDB on the human molecular environment.

The genetic etiology of periodic limb movement in sleep.

STUDY OBJECTIVES: Periodic limb movement in sleep is a common sleep phenotype characterized by repetitive leg movements that occur during or before sleep. We conducted a genome-wide association study (GWAS) of periodic limb movements in sleep (PLMS) using a joint analysis (i.e., discovery, replication, and joint meta-analysis) of four cohorts (MrOS, the Wisconsin Sleep Cohort Study, HypnoLaus, and MESA), comprised of 6843 total subjects. METHODS: The MrOS study and Wisconsin Sleep Cohort Study (N = 1745 cases) were used for discovery. Replication in the HypnoLaus and MESA cohorts (1002 cases) preceded joint meta-analysis. We also performed LD score regression, estimated heritability, and computed genetic correlations between potentially associated traits such as restless leg syndrome (RLS) and insomnia. The causality and direction of the relationships between PLMS and RLS was evaluated using Mendelian randomization. RESULTS: We found 2 independent loci were significantly associated with PLMS: rs113851554 (p = 3.51 × 10-12, ß = 0.486), an SNP located in a putative regulatory element of intron eight of MEIS1 (2p14); and rs9369062 (p = 3.06 × 10-22, ß = 0.2093), a SNP located in the intron region of BTBD9 (6p12); both of which were also lead signals in RLS GWAS. PLMS is genetically correlated with insomnia, risk of stroke, and RLS, but not with iron deficiency. Pleiotropy adjusted Mendelian randomization analysis identified a causal effect of RLS on PLMS. CONCLUSIONS: Because PLMS is more common than RLS, PLMS may have multiple causes and additional studies are needed to further validate these findings.

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies.

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.

Whole genome sequence analysis of apparent treatment resistant hypertension status in participants from the Trans-Omics for Precision Medicine program.

Introduction: Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data. Methods: Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90 mmHg, or four or more medications regardless of BP (n = 1,705). A normotensive control group was defined as individuals with BP < 140/90 mmHg (n = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90 mmHg (n = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American). Results: One variant in the known HTN locus, KCNK3, was a top finding in the multi-ethnic analysis (p = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74-0.88)]. This variant was replicated in the Vanderbilt University Medical Center's DNA repository data. Aggregate gene-based signals included the genes AGTPBP, MYL4, PDCD4, BBS9, ERG, and IER3. Discussion: Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.