A new test for trait mean and variance detects unreported loci for blood-pressure variation.

Variability in quantitative traits has clinical, ecological, and evolutionary significance. Most genetic variants identified for complex quantitative traits have only a detectable effect on the mean of trait. We have developed the mean-variance test (MVtest) to simultaneously model the mean and log-variance of a quantitative trait as functions of genotypes and covariates by using estimating equations. The advantages of MVtest include the facts that it can detect effect modification, that multiple testing can follow conventional thresholds, that it is robust to non-normal outcomes, and that association statistics can be meta-analyzed. In simulations, we show control of type I error of MVtest over several alternatives. We identified 51 and 37 previously unreported associations for effects on blood-pressure variance and mean, respectively, in the UK Biobank. Transcriptome-wide association studies revealed 633 significant unique gene associations with blood-pressure mean variance. MVtest is broadly applicable to studies of complex quantitative traits and provides an important opportunity to detect novel loci.

Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits.

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

EVALUATING THE RELATIONSHIPS BETWEEN GENETIC ANCESTRY AND THE CLINICAL PHENOME.

There is a desire in research to move away from the concept of race as a clinical factor because it is a societal construct used as an imprecise proxy for geographic ancestry. In this study, we leverage the biobank from Vanderbilt University Medical Center, BioVU, to investigate relationships between genetic ancestry proportion and the clinical phenome. For all samples in BioVU, we calculated six ancestry proportions based on 1000 Genomes references: eastern African (EAFR), western African (WAFR), northern European (NEUR), southern European (SEUR), eastern Asian (EAS), and southern Asian (SAS). From PheWAS, we found phecode categories significantly enriched neoplasms for EAFR, WAFR, and SEUR, and pregnancy complication in SEUR, NEUR, SAS, and EAS (p < 0.003). We then selected phenotypes hypertension (HTN) and atrial fibrillation (AFib) to further investigate the relationships between these phenotypes and EAFR, WAFR, SEUR, and NEUR using logistic regression modeling and non-linear restricted cubic spline modeling (RCS). For EAS and SAS, we chose renal failure (RF) for further modeling. The relationships between HTN and AFib and the ancestries EAFR, WAFR, and SEUR were best fit by the linear model (beta p < 1x10-4 for all) while the relationships with NEUR were best fit with RCS (HTN ANOVA p = 0.001, AFib ANOVA p < 1x10-4). For RF, the relationship with SAS was best fit with a linear model (beta p < 1x10-4) while RCS model was a better fit for EAS (ANOVA p < 1x10-4). In this study, we identify relationships between genetic ancestry and phenotypes that are best fit with non-linear modeling techniques. The assumption of linearity for regression modeling is integral for proper fitting of a model and there is no knowing a priori to modeling if the relationship is truly linear.

Systematic replication of smoking disease associations using survey responses and EHR data in the All of Us Research Program.

OBJECTIVE: The All of Us Research Program (All of Us) aims to recruit over a million participants to further precision medicine. Essential to the verification of biobanks is a replication of known associations to establish validity. Here, we evaluated how well All of Us data replicated known cigarette smoking associations. MATERIALS AND METHODS: We defined smoking exposure as follows: (1) an EHR Smoking exposure that used International Classification of Disease codes; (2) participant provided information (PPI) Ever Smoking; and, (3) PPI Current Smoking, both from the lifestyle survey. We performed a phenome-wide association study (PheWAS) for each smoking exposure measurement type. For each, we compared the effect sizes derived from the PheWAS to published meta-analyses that studied cigarette smoking from PubMed. We defined two levels of replication of meta-analyses: (1) nominally replicated: which required agreement of direction of effect size, and (2) fully replicated: which required overlap of confidence intervals. RESULTS: PheWASes with EHR Smoking, PPI Ever Smoking, and PPI Current Smoking revealed 736, 492, and 639 phenome-wide significant associations, respectively. We identified 165 meta-analyses representing 99 distinct phenotypes that could be matched to EHR phenotypes. At P < .05, 74 were nominally replicated and 55 were fully replicated. At P < 2.68 × 10-5 (Bonferroni threshold), 58 were nominally replicated and 40 were fully replicated. DISCUSSION: Most phenotypes found in published meta-analyses associated with smoking were nominally replicated in All of Us. Both survey and EHR definitions for smoking produced similar results. CONCLUSION: This study demonstrated the feasibility of studying common exposures using All of Us data.

Genetic architecture of heart failure with preserved versus reduced ejection fraction.

Pharmacologic clinical trials for heart failure with preserved ejection fraction have been largely unsuccessful as compared to those for heart failure with reduced ejection fraction. Whether differences in the genetic underpinnings of these major heart failure subtypes may provide insights into the disparate outcomes of clinical trials remains unknown. We utilize a large, uniformly phenotyped, single cohort of heart failure sub-classified into heart failure with reduced and with preserved ejection fractions based on current clinical definitions, to conduct detailed genetic analyses of the two heart failure sub-types. We find different genetic architectures and distinct genetic association profiles between heart failure with reduced and with preserved ejection fraction suggesting differences in underlying pathobiology. The modest genetic discovery for heart failure with preserved ejection fraction (one locus) compared to heart failure with reduced ejection fraction (13 loci) despite comparable sample sizes indicates that clinically defined heart failure with preserved ejection fraction likely represents the amalgamation of several, distinct pathobiological entities. Development of consensus sub-phenotyping of heart failure with preserved ejection fraction is paramount to better dissect the underlying genetic signals and contributors to this highly prevalent condition.

Large-Scale Multi-Omics Studies Provide New Insights into Blood Pressure Regulation.

Recent genome-wide association studies uncovered part of blood pressure's heritability. However, there is still a vast gap between genetics and biology that needs to be bridged. Here, we followed up blood pressure genome-wide summary statistics of over 750,000 individuals, leveraging comprehensive epigenomic and transcriptomic data from blood with a follow-up in cardiovascular tissues to prioritise likely causal genes and underlying blood pressure mechanisms. We first prioritised genes based on coding consequences, multilayer molecular associations, blood pressure-associated expression levels, and coregulation evidence. Next, we followed up the prioritised genes in multilayer studies of genomics, epigenomics, and transcriptomics, functional enrichment, and their potential suitability as drug targets. Our analyses yielded 1880 likely causal genes for blood pressure, tens of which are targets of the available licensed drugs. We identified 34 novel genes for blood pressure, supported by more than one source of biological evidence. Twenty-eight (82%) of these new genes were successfully replicated by transcriptome-wide association analyses in a large independent cohort (n = ~220,000). We also found a substantial mediating role for epigenetic regulation of the prioritised genes. Our results provide new insights into genetic regulation of blood pressure in terms of likely causal genes and involved biological pathways offering opportunities for future translation into clinical practice.

Using Mendelian randomisation to identify opportunities for type 2 diabetes prevention by repurposing medications used for lipid management.

BACKGROUND: Maintaining a healthy lifestyle to reduce type 2 diabetes (T2D) risk is challenging and additional strategies for T2D prevention are needed. We evaluated several lipid control medications as potential therapeutic options for T2D prevention using tissue-specific predicted gene expression summary statistics in a two-sample Mendelian randomisation (MR) design. METHODS: Large-scale European genome-wide summary statistics for lipids and T2D were leveraged in our multi-stage analysis to estimate changes in either lipid levels or T2D risk driven by tissue-specific predicted gene expression. We incorporated tissue-specific predicted gene expression summary statistics to proxy therapeutic effects of three lipid control medications [i.e., statins, icosapent ethyl (IPE), and proprotein convertase subtilisin/kexin type-9 inhibitors (PCSK-9i)] on T2D susceptibility using two-sample Mendelian randomisation (MR). FINDINGS: IPE, as proxied via increased FADS1 expression, was predicted to lower triglycerides and was associated with a 53% reduced risk of T2D. Statins and PCSK-9i, as proxied by reduced HMGCR and PCSK9 expression, respectively, were predicted to lower LDL-C levels but were not associated with T2D susceptibility. INTERPRETATION: Triglyceride lowering via IPE may reduce the risk of developing T2D in populations of European ancestry. However, experimental validation using animal models is needed to substantiate our results and to motivate randomized control trials (RCTs) for IPE as putative treatment for T2D prevention. FUNDING: Only summary statistics were used in this analysis. Funding information is detailed under Acknowledgments.

Diastolic Blood Pressure Alleles Improve Congenital Heart Defect Repair Outcomes.

BACKGROUND: Congenital heart defects (CHDs) affect 40 000 US births per year, half of which require surgical intervention. Individual differences in surgical outcomes including mortality and complications are not well understood but may be due to genetic variability. We hypothesized that polygenic risk scores (PRSs) for blood pressure in adults are associated with treatments and postsurgical outcomes in children with CHD, as CHD survivors are at higher risk of negative cardiometabolic disease. METHODS: We used imputed genotype data from pediatric participants requiring surgery for CHD (median age at surgery, 201 days; nmax=2498). Base data for the systolic and diastolic blood pressure PRSs (nmax=760 226) came from published genome-wide association study. The blood pressure PRSs were tested for association with postsurgical outcomes. All effects presented are per SD increase in PRS and adjusted for age, sex, body mass index, surgical complexity score, and first 10 principal components of ancestry. RESULTS: A higher diastolic blood pressure PRS was associated with decreased in-hospital mortality risk (odds ratio, 0.57 [0.39-0.82]; P=0.0022). Additional analyses suggest an interaction between diastolic blood pressure PRS and vasopressor dose. Those with a diastolic blood pressure PRS 1 SD above the mean, receiving a vasopressor dose in the top tertile, were estimated to have 52% (32%-66%) lower risk of in-hospital mortality compared with those with a vasopressor dose in the bottom tertile. CONCLUSIONS: These results suggest a genetically determined postsurgical survival advantage for CHD patients with blood pressure increasing alleles. Further study may reveal novel mechanisms contributing to postoperative morbidity and mortality, and this approach may assist in early identification of children at risk for adverse postoperative outcomes.

Population pharmacokinetic analysis of dexmedetomidine in children using real-world data from electronic health records and remnant specimens.

AIMS: Our objectives were to perform a population pharmacokinetic analysis of dexmedetomidine in children using remnant specimens and electronic health records (EHRs) and explore the impact of patient's characteristics and pharmacogenetics on dexmedetomidine clearance. METHODS: Dexmedetomidine dosing and patient data were gathered from EHRs and combined with opportunistically sampled remnant specimens. Population pharmacokinetic models were developed using nonlinear mixed-effects modelling. Stage 1 developed a model without genotype variables; Stage 2 added pharmacogenetic effects. RESULTS: Our final study population included 354 post-cardiac surgery patients aged 0-22 years (median 16 mo). The data were best described with a 2-compartment model with allometric scaling for weight and Hill maturation function for age. Population parameter estimates and 95% confidence intervals were 27.3 L/h (24.0-31.1 L/h) for total clearance, 161 L (139-187 L) for central compartment volume of distribution, 26.0 L/h (22.5-30.0 L/h) for intercompartmental clearance and 7903 L (5617-11 119 L) for peripheral compartment volume of distribution. The estimate for postmenstrual age when 50% of adult clearance is achieved was 42.0 weeks (41.5-42.5 weeks) and the Hill coefficient estimate was 7.04 (6.99-7.08). Genotype was not statistically or clinically significant. CONCLUSION: Our study demonstrates the use of real-world EHR data and remnant specimens to perform a population pharmacokinetic analysis and investigate covariate effects in a large paediatric population. Weight and age were important predictors of clearance. We did not find evidence for pharmacogenetic effects of UGT1A4 or UGT2B10 genotype or CYP2A6 risk score.

Evidence that geographic variation in genetic ancestry associates with uterine fibroids.

Uterine fibroids disproportionately impact Black women. Evidence suggests Black women have earlier onset and higher cumulative risk. This risk disparity may be due an imbalance of risk alleles in one parental geographic ancestry subgroup relative to others. We investigated ancestry proportions for the 1000 Genomes phase 3 populations clustered into six geographic groups for association with fibroid traits in Black women (n = 583 cases, 797 controls) and White women (n = 1195 cases, 1164 controls). Global ancestry proportions were estimated using ADMIXTURE. Dichotomous (fibroids status and multiple fibroid status) and continuous outcomes (volume and largest dimension) were modeled for association with ancestry proportions using logistic and linear regression adjusting for age. Effect estimates are reported per 10% increase in genetically inferred ancestry proportion. Among Black women, West African (WAFR) ancestry was associated with fibroid risk, East African ancestry was associated with risk of multiple fibroids, Northern European (NEUR) ancestry was protective for multiple fibroids, Southern European ancestry was protective for fibroids and multiple fibroids, and South Asian (SAS) ancestry was positively associated with volume and largest dimension. In White women, NEUR ancestry was protective for fibroids, SAS ancestry was associated with fibroid risk, and WAFR ancestry was positively associated with volume and largest dimension. These results suggest that a proportion of fibroid risk and fibroid trait racial disparities are due to genetic differences between geographic groups. Further investigation at the local ancestry and single variant levels may yield novel insights into disease architecture and genetic mechanisms underlying ethnic disparities in fibroid risk.

Association Between Genetic Variation in Blood Pressure and Increased Lifetime Risk of Peripheral Artery Disease.

[Figure: see text].

Associations of biogeographic ancestry with hypertension traits.

OBJECTIVES: Ethnic disparities in hypertension prevalence are well documented, though the influence of genetic ancestry is unclear. The aim of this study was to evaluate associations of geographic genetic ancestry with hypertension and underlying blood pressure traits. METHODS: We tested genetically inferred ancestry proportions from five 1000 Genomes reference populations (GBR, PEL, YRI, CHB, and LWK) for association with four continuous blood pressure (BP) traits (SBP, DBP, PP, MAP) and the dichotomous outcomes hypertension and apparent treatment-resistant hypertension in 220 495 European American, 59 927 African American, and 21 273 Hispanic American individuals from the Million Veteran Program. Ethnicity stratified results were meta-analyzed to report effect estimates per 10% difference for a given ancestry proportion in all samples. RESULTS: Percentage GBR was negatively associated with BP (P = 2.13 × 10-19, 7.92 × 10-8, 4.41 × 10-11, and 3.57 × 10-13 for SBP, DBP, PP, and MAP, respectively; coefficient range -0.10 to -0.21 mmHg per 10% increase in ancestry proportion) and was protective against hypertension [P = 2.59 × 10-5, odds ratio (OR) = 0.98] relative to other ancestries. YRI percentage was positively associated with BP (P = 1.63 × 10-23, 1.94 × 10-26, 0.012, and 3.26 × 10-29 for SBP, DBP, PP, and MAP, respectively; coefficient range 0.06-0.32 mmHg per 10% increase in ancestry proportion) and was positively associated with hypertension risk (P = 3.10 × 10-11, OR = 1.04) and apparent treatment-resistant hypertension risk (P = 1.86 × 10-4, OR = 1.04) compared with other ancestries. Percentage PEL was inversely associated with DBP (P = 2.84 × 10-5, beta = -0.11 mmHg per 10% increase in ancestry proportion). CONCLUSION: These results demonstrate that risk for BP traits varies significantly by genetic ancestry. Our findings provide insight into the geographic origin of genetic factors underlying hypertension risk and establish that a portion of BP trait ethnic disparities are because of genetic differences between ancestries.

Discovery of 318 new risk loci for type 2 diabetes and related vascular outcomes among 1.4 million participants in a multi-ancestry meta-analysis.

We investigated type 2 diabetes (T2D) genetic susceptibility via multi-ancestry meta-analysis of 228,499 cases and 1,178,783 controls in the Million Veteran Program (MVP), DIAMANTE, Biobank Japan and other studies. We report 568 associations, including 286 autosomal, 7 X-chromosomal and 25 identified in ancestry-specific analyses that were previously unreported. Transcriptome-wide association analysis detected 3,568 T2D associations with genetically predicted gene expression in 687 novel genes; of these, 54 are known to interact with FDA-approved drugs. A polygenic risk score (PRS) was strongly associated with increased risk of T2D-related retinopathy and modestly associated with chronic kidney disease (CKD), peripheral artery disease (PAD) and neuropathy. We investigated the genetic etiology of T2D-related vascular outcomes in the MVP and observed statistical SNP-T2D interactions at 13 variants, including coronary heart disease (CHD), CKD, PAD and neuropathy. These findings may help to identify potential therapeutic targets for T2D and genomic pathways that link T2D to vascular outcomes.

A genome-wide association study of polycystic ovary syndrome identified from electronic health records.

BACKGROUND: Polycystic ovary syndrome is the most common endocrine disorder affecting women of reproductive age. A number of criteria have been developed for clinical diagnosis of polycystic ovary syndrome, with the Rotterdam criteria being the most inclusive. Evidence suggests that polycystic ovary syndrome is significantly heritable, and previous studies have identified genetic variants associated with polycystic ovary syndrome diagnosed using different criteria. The widely adopted electronic health record system provides an opportunity to identify patients with polycystic ovary syndrome using the Rotterdam criteria for genetic studies. OBJECTIVE: To identify novel associated genetic variants under the same phenotype definition, we extracted polycystic ovary syndrome cases and unaffected controls based on the Rotterdam criteria from the electronic health records and performed a discovery-validation genome-wide association study. STUDY DESIGN: We developed a polycystic ovary syndrome phenotyping algorithm on the basis of the Rotterdam criteria and applied it to 3 electronic health record-linked biobanks to identify cases and controls for genetic study. In the discovery phase, we performed an individual genome-wide association study using the Geisinger MyCode and the Electronic Medical Records and Genomics cohorts, which were then meta-analyzed. We attempted validation of the significant association loci (P<1×10-6) in the BioVU cohort. All association analyses used logistic regression, assuming an additive genetic model, and adjusted for principal components to control for population stratification. An inverse-variance fixed-effect model was adopted for meta-analysis. In addition, we examined the top variants to evaluate their associations with each criterion in the phenotyping algorithm. We used the STRING database to characterize protein-protein interaction network. RESULTS: Using the same algorithm based on the Rotterdam criteria, we identified 2995 patients with polycystic ovary syndrome and 53,599 population controls in total (2742 cases and 51,438 controls from the discovery phase; 253 cases and 2161 controls in the validation phase). We identified 1 novel genome-wide significant variant rs17186366 (odds ratio [OR]=1.37 [1.23, 1.54], P=2.8×10-8) located near SOD2. In addition, 2 loci with suggestive association were also identified: rs113168128 (OR=1.72 [1.42, 2.10], P=5.2×10-8), an intronic variant of ERBB4 that is independent from the previously published variants, and rs144248326 (OR=2.13 [1.52, 2.86], P=8.45×10-7), a novel intronic variant in WWTR1. In the further association tests of the top 3 single-nucleotide polymorphisms with each criterion in the polycystic ovary syndrome algorithm, we found that rs17186366 (SOD2) was associated with polycystic ovaries and hyperandrogenism, whereas rs11316812 (ERBB4) and rs144248326 (WWTR1) were mainly associated with oligomenorrhea or infertility. We also validated the previously reported association with DENND1A1. Using the STRING database to characterize protein-protein interactions, we found both ERBB4 and WWTR1 can interact with YAP1, which has been previously associated with polycystic ovary syndrome. CONCLUSION: Through a discovery-validation genome-wide association study on polycystic ovary syndrome identified from electronic health records using an algorithm based on Rotterdam criteria, we identified and validated a novel genome-wide significant association with a variant near SOD2. We also identified a novel independent variant within ERBB4 and a suggestive association with WWTR1. With previously identified polycystic ovary syndrome gene YAP1, the ERBB4-YAP1-WWTR1 network suggests involvement of the epidermal growth factor receptor and the Hippo pathway in the multifactorial etiology of polycystic ovary syndrome.

Mapping eGFR loci to the renal transcriptome and phenome in the VA Million Veteran Program.

Chronic kidney disease (CKD), defined by low estimated glomerular filtration rate (eGFR), contributes to global morbidity and mortality. Here we conduct a transethnic Genome-Wide Association Study of eGFR in 280,722 participants of the Million Veteran Program (MVP), with replication in 765,289 participants from the Chronic Kidney Disease Genetics (CKDGen) Consortium. We identify 82 previously unreported variants, confirm 54 loci, and report interesting findings including association of the sickle cell allele of betaglobin among non-Hispanic blacks. Our transcriptome-wide association study of kidney function in healthy kidney tissue identifies 36 previously unreported and nine known genes, and maps gene expression to renal cell types. In a Phenome-Wide Association Study in 192,868 MVP participants using a weighted genetic score we detect associations with CKD stages and complications and kidney stones. This investigation reinterprets the genetic architecture of kidney function to identify the gene, tissue, and anatomical context of renal homeostasis and the clinical consequences of dysregulation.

A Trans-Ethnic Genome-Wide Association Study of Uterine Fibroids.

Uterine fibroids affect up to 77% of women by menopause and account for up to $34 billion in healthcare costs each year. Although fibroid risk is heritable, genetic risk for fibroids is not well understood. We conducted a two-stage case-control meta-analysis of genetic variants in European and African ancestry women with and without fibroids classified by a previously published algorithm requiring pelvic imaging or confirmed diagnosis. Women from seven electronic Medical Records and Genomics (eMERGE) network sites (3,704 imaging-confirmed cases and 5,591 imaging-confirmed controls) and women of African and European ancestry from UK Biobank (UKB, 5,772 cases and 61,457 controls) were included in the discovery genome-wide association study (GWAS) meta-analysis. Variants showing evidence of association in Stage I GWAS (P < 1 × 10-5) were targeted in an independent replication sample of African and European ancestry individuals from the UKB (Stage II) (12,358 cases and 138,477 controls). Logistic regression models were fit with genetic markers imputed to a 1000 Genomes reference and adjusted for principal components for each race- and site-specific dataset, followed by fixed-effects meta-analysis. Final analysis with 21,804 cases and 205,525 controls identified 326 genome-wide significant variants in 11 loci, with three novel loci at chromosome 1q24 (sentinel-SNP rs14361789; P = 4.7 × 10-8), chromosome 16q12.1 (sentinel-SNP rs4785384; P = 1.5 × 10-9) and chromosome 20q13.1 (sentinel-SNP rs6094982; P = 2.6 × 10-8). Our statistically significant findings further support previously reported loci including SNPs near WT1, TNRC6B, SYNE1, BET1L, and CDC42/WNT4. We report evidence of ancestry-specific findings for sentinel-SNP rs10917151 in the CDC42/WNT4 locus (P = 1.76 × 10-24). Ancestry-specific effect-estimates for rs10917151 were in opposite directions (P-Het-between-groups = 0.04) for predominantly African (OR = 0.84) and predominantly European women (OR = 1.16). Genetically-predicted gene expression of several genes including LUZP1 in vagina (P = 4.6 × 10-8), OBFC1 in esophageal mucosa (P = 8.7 × 10-8), NUDT13 in multiple tissues including subcutaneous adipose tissue (P = 3.3 × 10-6), and HEATR3 in skeletal muscle tissue (P = 5.8 × 10-6) were associated with fibroids. The finding for HEATR3 was supported by SNP-based summary Mendelian randomization analysis. Our study suggests that fibroid risk variants act through regulatory mechanisms affecting gene expression and are comprised of alleles that are both ancestry-specific and shared across continental ancestries.

Genome-wide association study identifies novel loci for type 2 diabetes-attributed end-stage kidney disease in African Americans.

BACKGROUND: End-stage kidney disease (ESKD) is a significant public health concern disproportionately affecting African Americans (AAs). Type 2 diabetes (T2D) is the leading cause of ESKD in the USA, and efforts to uncover genetic susceptibility to diabetic kidney disease (DKD) have had limited success. A prior genome-wide association study (GWAS) in AAs with T2D-ESKD was expanded with additional AA cases and controls and genotypes imputed to the higher density 1000 Genomes reference panel. The discovery analysis included 3432 T2D-ESKD cases and 6977 non-diabetic non-nephropathy controls (N = 10,409), followed by a discrimination analysis in 2756 T2D non-nephropathy controls to exclude T2D-associated variants. RESULTS: Six independent variants located in or near RND3/RBM43, SLITRK3, ENPP7, GNG7, and APOL1 achieved genome-wide significant association (P < 5 × 10-8) with T2D-ESKD. Following extension analyses in 1910 non-diabetic ESKD cases and 908 non-diabetic non-nephropathy controls, a meta-analysis of 5342 AA all-cause ESKD cases and 6977 AA non-diabetic non-nephropathy controls revealed an additional novel all-cause ESKD locus at EFNB2 (rs77113398; P = 9.84 × 10-9; OR = 1.94). Exclusion of APOL1 renal-risk genotype carriers identified two additional genome-wide significant T2D-ESKD-associated loci at GRAMD3 and MGAT4C. A second variant at GNG7 (rs373971520; P = 2.17 × 10-8, OR = 1.46) remained associated with all-cause ESKD in the APOL1-negative analysis. CONCLUSIONS: Findings provide further evidence for genetic factors associated with advanced kidney disease in AAs with T2D.

Trans-ethnic association study of blood pressure determinants in over 750,000 individuals.

In this trans-ethnic multi-omic study, we reinterpret the genetic architecture of blood pressure to identify genes, tissues, phenomes and medication contexts of blood pressure homeostasis. We discovered 208 novel common blood pressure SNPs and 53 rare variants in genome-wide association studies of systolic, diastolic and pulse pressure in up to 776,078 participants from the Million Veteran Program (MVP) and collaborating studies, with analysis of the blood pressure clinical phenome in MVP. Our transcriptome-wide association study detected 4,043 blood pressure associations with genetically predicted gene expression of 840 genes in 45 tissues, and mouse renal single-cell RNA sequencing identified upregulated blood pressure genes in kidney tubule cells.

An Exome-wide Association Study for Type 2 Diabetes-Attributed End-Stage Kidney Disease in African Americans.

INTRODUCTION: Compared with European Americans, African Americans (AAs) are at higher risk for developing end-stage kidney disease (ESKD). Genome-wide association studies (GWAS) have identified >70 genetic variants associated with kidney function and chronic kidney disease (CKD) in patients with and without diabetes. However, these variants explain a small proportion of disease liability. This study examined the contribution of coding genetic variants for risk of type 2 diabetes (T2D)-attributed ESKD and advanced CKD in AAs. METHODS: Exome sequencing was performed in 456 AA T2D-ESKD cases, and 936 AA nondiabetic, non-nephropathy control individuals at the discovery stage. A mixed logistic regression model was used for association analysis. Nominal associations (P < 0.05) were replicated in an additional 2020 T2D-ESKD cases and 1121 nondiabetic, non-nephropathy control individuals. A meta-analysis combining 4533 discovery and replication samples was performed. Putative T2D-ESKD associations were tested in additional 1910 nondiabetic ESKD and 219 T2D-ESKD cases, as well as 912 AA nondiabetic non-nephropathy control individuals. RESULTS: A total of 11 suggestive T2D-ESKD associations (P < 1 x 10-4) from 8 loci (PLEKHN1, NADK, RAD51AP2, RREB1, PEX6, GRM8, PRX, APOL1) were apparent in the meta-analysis. Exclusion of APOL1 renal-risk genotype carriers identified 3 additional suggestive loci (OTUD7B, IFITM3, DLGAP5). Rs41302867 in RREB1 displayed consistent association with T2D-ESKD and nondiabetic ESKD (odds ratio: 0.47; P = 1.2 x 10-6 in 4605 all-cause ESKD and 2969 nondiabetic non-nephropathy control individuals). CONCLUSION: Our findings suggest that coding genetic variants are implicated in predisposition to T2D-ESKD in AAs.

Genome-wide interaction with the insulin secretion locus MTNR1B reveals CMIP as a novel type 2 diabetes susceptibility gene in African Americans.

Although type 2 diabetes (T2D) results from metabolic defects in insulin secretion and insulin sensitivity, most of the genetic risk loci identified to date relates to insulin secretion. We reported that T2D loci influencing insulin sensitivity may be identified through interactions with insulin secretion loci, thereby leading to T2D. Here, we hypothesize that joint testing of variant main effects and interaction effects with an insulin secretion locus increases power to identify genetic interactions leading to T2D. We tested this hypothesis with an intronic MTNR1B SNP, rs10830963, which is associated with acute insulin response to glucose, a dynamic measure of insulin secretion. rs10830963 was tested for interaction and joint (main + interaction) effects with genome-wide data in African Americans (2,452 cases and 3,772 controls) from five cohorts. Genome-wide genotype data (Affymetrix Human Genome 6.0 array) was imputed to a 1000 Genomes Project reference panel. T2D risk was modeled using logistic regression with rs10830963 dosage, age, sex, and principal component as predictors. Joint effects were captured using the Kraft two degrees of freedom test. Genome-wide significant (P < 5 × 10-8 ) interaction with MTNR1B and joint effects were detected for CMIP intronic SNP rs17197883 (Pinteraction  = 1.43 × 10-8 ; Pjoint  = 4.70 × 10-8 ). CMIP variants have been nominally associated with T2D, fasting glucose, and adiponectin in individuals of East Asian ancestry, with high-density lipoprotein, and with waist-to-hip ratio adjusted for body mass index in Europeans. These data support the hypothesis that additional genetic factors contributing to T2D risk, including insulin sensitivity loci, can be identified through interactions with insulin secretion loci.