High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

Red blood cell (RBC) traits are important heritable clinical biomarkers and modifiers of disease severity. To identify coding genetic variants associated with these traits, we conducted meta-analyses of seven RBC phenotypes in 130,273 multi-ethnic individuals from studies genotyped on an exome array. After conditional analyses and replication in 27,480 independent individuals, we identified 16 new RBC variants. We found low-frequency missense variants in MAP1A (rs55707100, minor allele frequency [MAF] = 3.3%, p = 2 × 10(-10) for hemoglobin [HGB]) and HNF4A (rs1800961, MAF = 2.4%, p < 3 × 10(-8) for hematocrit [HCT] and HGB). In African Americans, we identified a nonsense variant in CD36 associated with higher RBC distribution width (rs3211938, MAF = 8.7%, p = 7 × 10(-11)) and showed that it is associated with lower CD36 expression and strong allelic imbalance in ex vivo differentiated human erythroblasts. We also identified a rare missense variant in ALAS2 (rs201062903, MAF = 0.2%) associated with lower mean corpuscular volume and mean corpuscular hemoglobin (p < 8 × 10(-9)). Mendelian mutations in ALAS2 are a cause of sideroblastic anemia and erythropoietic protoporphyria. Gene-based testing highlighted three rare missense variants in PKLR, a gene mutated in Mendelian non-spherocytic hemolytic anemia, associated with HGB and HCT (SKAT p < 8 × 10(-7)). These rare, low-frequency, and common RBC variants showed pleiotropy, being also associated with platelet, white blood cell, and lipid traits. Our association results and functional annotation suggest the involvement of new genes in human erythropoiesis. We also confirm that rare and low-frequency variants play a role in the architecture of complex human traits, although their phenotypic effect is generally smaller than originally anticipated.

Evidence suggests European American (EA) women have two- to five-fold increased odds of having pelvic organ prolapse (POP) when compared with African American (AA) women. However, the role of genetic ancestry in relation to POP risk is not clear. Here we evaluate the association between genetic ancestry and POP in AA women from the Women's Health Initiative Hormone Therapy trial. Women with grade 1 or higher classification, and grade 2 or higher classification for uterine prolapse, cystocele or rectocele at baseline or during follow-up were considered to have any POP (N = 805) and moderate/severe POP (N = 156), respectively. Women with at least two pelvic exams with no indication for POP served as controls (N = 344). We performed case-only, and case-control admixture-mapping analyses using multiple logistic regression while adjusting for age, BMI, parity and global ancestry. We evaluated the association between global ancestry and POP using multiple logistic regression. European ancestry at the individual level was not associated with POP risk. Case-only and case-control local ancestry analyses identified two ancestry-specific loci that may be associated with POP. One locus (Chromosome 15q26.2) achieved empirically-estimated statistical significance and was associated with decreased POP odds (considering grade ≥2 POP) with each unit increase in European ancestry (OR: 0.35; 95% CI: 0.30, 0.57; p-value = 1.48x10-5). This region includes RGMA, a potent regulator of the BMP family of genes. The second locus (Chromosome 1q42.1-q42.3) was associated with increased POP odds with each unit increase in European ancestry (Odds ratio [OR]: 1.69; 95% confidence interval [CI]: 1.28, 2.22; p-value = 1.93x10-4). Although this region did not reach statistical significance after considering multiple comparisons, it includes potentially relevant genes including TBCE, and ACTA1. Unique non-overlapping European and African ancestry-specific susceptibility loci may be associated with increased POP risk.

Uterine fibroids are benign tumors of the uterus affecting up to 77% of women by menopause. They are the leading indication for hysterectomy, and account for $34 billion annually in the United States. Race/ethnicity and age are the strongest known risk factors. African American (AA) women have higher prevalence, earlier onset, and larger and more numerous fibroids than European American women. We conducted a multi-stage genome-wide association study (GWAS) of fibroid risk among AA women followed by in silico genetically predicted gene expression profiling of top hits. In Stage 1, cases and controls were confirmed by pelvic imaging, genotyped and imputed to 1000 Genomes. Stage 2 used self-reported fibroid and GWAS data from 23andMe, Inc. and the Black Women's Health Study. Associations with fibroid risk were modeled using logistic regression adjusted for principal components, followed by meta-analysis of results. We observed a significant association among 3399 AA cases and 4764 AA controls at rs739187 (risk-allele frequency = 0.27) in CYTH4 (OR (95% confidence interval) = 1.23 (1.16-1.30), p value = 7.82 × 10-9). Evaluation of the genetic association results with MetaXcan identified lower predicted gene expression of CYTH4 in thyroid tissue as significantly associated with fibroid risk (p value = 5.86 × 10-8). In this first multi-stage GWAS for fibroids among AA women, we identified a novel risk locus for fibroids within CYTH4 that impacts gene expression in thyroid and has potential biological relevance for fibroids.

Increasingly, clinical phenotypes with matched genetic data from bio-bank linked electronic health records (EHRs) have been used for pleiotropy analyses. Thus far, pleiotropy analysis using individual-level EHR data has been limited to data from one site. However, it is desirable to integrate EHR data from multiple sites to improve the detection power and generalizability of the results. Due to privacy concerns, individual-level patients' data are not easily shared across institutions. As a result, we introduce Sum-Share, a method designed to efficiently integrate EHR and genetic data from multiple sites to perform pleiotropy analysis. Sum-Share requires only summary-level data and one round of communication from each site, yet it produces identical test statistics compared with that of pooled individual-level data. Consequently, Sum-Share can achieve lossless integration of multiple datasets. Using real EHR data from eMERGE, Sum-Share is able to identify 1734 potential pleiotropic SNPs for five cardiovascular diseases.

Polycystic ovary syndrome (PCOS) is one of the leading causes of infertility, yet current diagnostic criteria are ineffective at identifying patients whose symptoms reside outside strict diagnostic criteria. As a result, PCOS is underdiagnosed and its etiology is poorly understood.

Plasma lipids are known heritable risk factors for cardiovascular disease, but increasing evidence also supports shared genetics with diseases of other organ systems. We devised a comprehensive three-phase framework to identify new lipid-associated genes and study the relationships among lipids, genotypes, gene expression and hundreds of complex human diseases from the Electronic Medical Records and Genomics (347 traits) and the UK Biobank (549 traits). Aside from 67 new lipid-associated genes with strong replication, we found evidence for pleiotropic SNPs/genes between lipids and diseases across the phenome. These include discordant pleiotropy in the HLA region between lipids and multiple sclerosis and putative causal paths between triglycerides and gout, among several others. Our findings give insights into the genetic basis of the relationship between plasma lipids and diseases on a phenome-wide scale and can provide context for future prevention and treatment strategies.

Multimorbidity, characterized by the simultaneous occurrence of multiple diseases in an individual, is an increasing global health concern, posing substantial challenges to healthcare systems. Comprehensive understanding of disease-disease interactions and intrinsic mechanisms behind multimorbidity can offer opportunities for innovative prevention strategies, targeted interventions, and personalized treatments. Yet, there exist limited tools and datasets that characterize multimorbidity patterns across different populations. To bridge this gap, we used large-scale electronic health record (EHR) systems to develop the Phenome-wide Multi-Institutional Multimorbidity Explorer (PheMIME), which facilitates research in exploring and comparing multimorbidity patterns among multiple institutions, potentially leading to the discovery of novel and robust disease associations and patterns that are interoperable across different systems and organizations.

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.

Benign prostatic hyperplasia (BPH) results in a significant public health burden due to the morbidity caused by the disease and many of the available remedies. As much as 70% of men over 70 will develop BPH. Few studies have been conducted to discover the genetic determinants of BPH risk. Understanding the biological basis for this condition may provide necessary insight for development of novel pharmaceutical therapies or risk prediction. We have evaluated SNP-based heritability of BPH in two cohorts and conducted a genome-wide association study (GWAS) of BPH risk using 2,656 cases and 7,763 controls identified from the Electronic Medical Records and Genomics (eMERGE) network. SNP-based heritability estimates suggest that roughly 60% of the phenotypic variation in BPH is accounted for by genetic factors. We used logistic regression to model BPH risk as a function of principal components of ancestry, age, and imputed genotype data, with meta-analysis performed using METAL. The top result was on chromosome 22 in SYN3 at rs2710383 (p-value = 4.6 × 10-7; Odds Ratio = 0.69, 95% confidence interval = 0.55-0.83). Other suggestive signals were near genes GLGC, UNCA13, SORCS1 and between BTBD3 and SPTLC3. We also evaluated genetically-predicted gene expression in prostate tissue. The most significant result was with increasing predicted expression of ETV4 (chr17; p-value = 0.0015). Overexpression of this gene has been associated with poor prognosis in prostate cancer. In conclusion, although there were no genome-wide significant variants identified for BPH susceptibility, we present evidence supporting the heritability of this phenotype, have identified suggestive signals, and evaluated the association between BPH and genetically-predicted gene expression in prostate.

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.

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.

We examined associations between suicidality and genotypes that predict plasma efavirenz exposure among AIDS Clinical Trials Group study participants in the United States.

Osteoporosis is a major public health problem associated with excess disability and mortality. It is estimated that 50-70% of the variation in osteoporotic fracture risk is attributable to genetic factors. The purpose of this hypothesis-generating study was to identify possible genetic determinants of fracture among African American (AA) women in a GWAS meta-analysis.

Endometriosis is a heritable hormone-dependent gynecological disorder, associated with severe pelvic pain and reduced fertility; however, its molecular mechanisms remain largely unknown. Here we perform a meta-analysis of 11 genome-wide association case-control data sets, totalling 17,045 endometriosis cases and 191,596 controls. In addition to replicating previously reported loci, we identify five novel loci significantly associated with endometriosis risk (P<5 × 10-8), implicating genes involved in sex steroid hormone pathways (FN1, CCDC170, ESR1, SYNE1 and FSHB). Conditional analysis identified five secondary association signals, including two at the ESR1 locus, resulting in 19 independent single nucleotide polymorphisms (SNPs) robustly associated with endometriosis, which together explain up to 5.19% of variance in endometriosis. These results highlight novel variants in or near specific genes with important roles in sex steroid hormone signalling and function, and offer unique opportunities for more targeted functional research efforts.

Resting metabolic rates (RMR) vary across individuals. Understanding the determinants of RMR could provide biological insight into obesity and its metabolic consequences such as type 2 diabetes and cardiovascular diseases.

Engagement of integrins by the extracellular matrix initiates signaling cascades that drive a variety of cellular functions, including neuronal migration and axonal pathfinding in the brain. Multiple lines of evidence link the ITGB3 gene encoding the integrin β3 subunit with the serotonin (5-HT) system, likely via its modulation of the 5-HT transporter (SERT). The ITGB3 coding polymorphism Leu33Pro (rs5918, PlA2) produces hyperactive αvβ3 receptors that influence whole-blood 5-HT levels and may influence the risk for autism spectrum disorder (ASD). Using a phenome-wide scan of psychiatric diagnoses, we found significant, male-specific associations between the Pro33 allele and attention-deficit hyperactivity disorder and ASDs. Here, we used knock-in (KI) mice expressing an Itgb3 variant that phenocopies the human Pro33 variant to elucidate the consequences of constitutively enhanced αvβ3 signaling to the 5-HT system in the brain. KI mice displayed deficits in multiple behaviors, including anxiety, repetitive, and social behaviors. Anatomical studies revealed a significant decrease in 5-HT synapses in the midbrain, accompanied by decreases in SERT activity and reduced localization of SERTs to integrin adhesion complexes in synapses of KI mice. Inhibition of focal adhesion kinase (FAK) rescued SERT function in synapses of KI mice, demonstrating that constitutive active FAK signaling downstream of the Pro32Pro33 integrin αvβ3 suppresses SERT activity. Our studies identify a complex regulation of 5-HT homeostasis and behaviors by integrin αvβ3, revealing an important role for integrins in modulating risk for neuropsychiatric disorders.SIGNIFICANCE STATEMENT The integrin β3 Leu33Pro coding polymorphism has been associated with autism spectrum disorders (ASDs) within a subgroup of patients with elevated blood 5-HT levels, linking integrin β3, 5-HT, and ASD risk. We capitalized on these interactions to demonstrate that the Pro33 coding variation in the murine integrin β3 recapitulates the sex-dependent neurochemical and behavioral attributes of ASD. Using state-of-the-art techniques, we show that presynaptic 5-HT function is altered in these mice, and that the localization of 5-HT transporters to specific compartments within the synapse, disrupted by the integrin β3 Pro33 mutation, is critical for appropriate reuptake of 5-HT. Our studies provide fundamental insight into the genetic network regulating 5-HT neurotransmission in the CNS that is also associated with ASD risk.

The monocyte chemotactic protein-1 (MCP-1) is a chemokine that plays an important role in the recruitment of monocytes to M. tuberculosis infection sites, and previous studies have reported that genetic variants in MCP1 are associated with differential susceptibility to pulmonary tuberculosis (PTB). We examined eight MCP1 single nucleotide polymorphisms (SNPs) in a multi-ethnic, case-control design that included: 321 cases and 346 controls from Guinea-Bissau, 258 cases and 271 controls from The Gambia, 295 cases and 179 controls from the U.S. (African-Americans), and an additional set of 237 cases and 144 controls of European ancestry from the U.S. and Argentina. Two locus interactions were also examined for polymorphisms in MCP1 and interleukin 12B (IL12B), another gene implicated in PTB risk. Examination of previously associated MCP1 SNPs rs1024611 (-2581A/G), rs2857656 (-362G/C) and rs4586 (+900C/T) did not show evidence for association. One interaction between rs2857656 and IL12B SNP rs2288831 was observed among Africans but the effect was in the opposite direction in Guineans (OR = 1.90, p = 0.001) and Gambians (OR = 0.64, p = 0.024). Our data indicate that the effect of genetic variation within MCP1 is not clear cut and additional studies will be needed to elucidate its role in TB susceptibility.

Meta-analyses of association results for blood pressure using exome-centric single-variant and gene-based tests identified 31 new loci in a discovery stage among 146,562 individuals, with follow-up and meta-analysis in 180,726 additional individuals (total n = 327,288). These blood pressure-associated loci are enriched for known variants for cardiometabolic traits. Associations were also observed for the aggregation of rare and low-frequency missense variants in three genes, NPR1, DBH, and PTPMT1. In addition, blood pressure associations at 39 previously reported loci were confirmed. The identified variants implicate biological pathways related to cardiometabolic traits, vascular function, and development. Several new variants are inferred to have roles in transcription or as hubs in protein-protein interaction networks. Genetic risk scores constructed from the identified variants were strongly associated with coronary disease and myocardial infarction. This large collection of blood pressure-associated loci suggests new therapeutic strategies for hypertension, emphasizing a link with cardiometabolic risk.

Increasingly, patients without clinical indications are undergoing genomic tests. The purpose of this study was to assess their appreciation and comprehension of their test results and their clinicians' reactions. We conducted 675 surveys with participants from the Vanderbilt Electronic Medical Records and Genomics (eMERGE) cohort. We interviewed 36 participants: 19 had received positive results, and 17 were self-identified racial minorities. Eleven clinicians who had patients who had participated in eMERGE were interviewed. A further 21 of these clinicians completed surveys. Participants spontaneously admitted to understanding little or none of the information returned to them from the eMERGE study. However, they simultaneously said that they generally found testing to be "helpful," even when it did not inform their health care. Primary care physicians expressed discomfort in being asked to interpret the results for their patients and described it as an undue burden. Providing genetic testing to otherwise healthy patients raises a number of ethical issues that warrant serious consideration. Although our participants were enthusiastic about enrolling and receiving their results, they express a limited understanding of what the results mean for their health care. This fact, coupled the clinicians' concern, urges greater caution when educating and enrolling participants in clinically non-indicated testing.