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.

Uterine leiomyomata and keloids fibrosis origins: a mini-review of fibroproliferative diseases.

Diseases such as uterine leiomyomata (fibroids and benign tumors of the uterus) and keloids (raised scars) may share common etiology. Fibroids and keloids can co-occur in individuals, and both are highly heritable, suggesting they may share common genetic risk factors. Fibroproliferative diseases are common and characterized by scarring and overgrowth of connective tissue, impacting multiple organ systems. These conditions both have racial disparities in prevalence, with the highest prevalence observed among individuals of African ancestry. Several fibroproliferative diseases are more severe and common in populations of sub-Saharan Africa. This mini-review aims to provide a broad overview of the current knowledge of the evolutionary origins and causes of fibroproliferative diseases. We also discuss current hypotheses proposing that the increased prevalence of these diseases in African-derived populations is due to the selection for profibrotic alleles that are protective against helminth infections and provide examples from knowledge of uterine fibroid and keloid research.

Predictive models for abdominal aortic aneurysms using polygenic scores and PheWAS-derived risk factors.

Abdominal aortic aneurysms (AAA) are common enlargements of the abdominal aorta which can grow larger until rupture, often leading to death. Detection of AAA is often by ultrasonography and screening recommendations are mostly directed at men over 65 with a smoking history. Recent large-scale genome-wide association studies have identified genetic loci associated with AAA risk. We combined known risk factors, polygenic risk scores (PRS) and precedent clinical diagnoses from electronic health records (EHR) to develop predictive models for AAA, and compared performance against screening recommendations. The PRS included genome-wide summary statistics from the Million Veteran Program and FinnGen (10,467 cases, 378,713 controls of European ancestry), with optimization in Vanderbilt's BioVU and validated in the eMERGE Network, separately across both White and Black participants. Candidate diagnoses were identified through a temporally-oriented Phenome-wide association study in independent EHR data from Vanderbilt, and features were selected via elastic net. We calculated C-statistics in eMERGE for models including PRS, phecodes, and covariates using regression weights from BioVU. The AUC for the full model in the test set was 0.883 (95% CI 0.873-0.892), 0.844 (0.836-0.851) for covariates only, 0.613 (95% CI 0.604-0.622) when using primary USPSTF screening criteria, and 0.632 (95% CI 0.623-0.642) using primary and secondary criteria. Brier scores were between 0.003 and 0.023 for our models indicating good calibration, and net reclassification improvement over combined primary and secondary USPSTF criteria was 0.36-0.60. We provide PRS for AAA which are strongly associated with AAA risk and add to predictive model performance. These models substantially improve identification of people at risk of a AAA diagnosis compared with existing guidelines, with evidence of potential applicability in minority populations.

Uterine fibroid polygenic risk score (PRS) associates and predicts risk for uterine fibroid.

Uterine fibroids (UF) are common pelvic tumors in women, heritable, and genome-wide association studies (GWAS) have identified ~ 30 loci associated with increased risk in UF. Using summary statistics from a previously published UF GWAS performed in a non-Hispanic European Ancestry (NHW) female subset from the Electronic Medical Records and Genomics (eMERGE) Network, we constructed a polygenic risk score (PRS) for UF. UF-PRS was developed using PRSice and optimized in the separate clinical population of BioVU. PRS was validated using parallel methods of 10-fold cross-validation logistic regression and phenome-wide association study (PheWAS) in a seperate subset of eMERGE NHW females (validation set), excluding samples used in GWAS. PRSice determined pt < 0.001 and after linkage disequilibrium pruning (r2 < 0.2), 4458 variants were in the PRS which was significant (pseudo-R2 = 0.0018, p = 0.041). 10-fold cross-validation logistic regression modeling of validation set revealed the model had an area under the curve (AUC) value of 0.60 (95% confidence interval [CI] 0.58-0.62) when plotted in a receiver operator curve (ROC). PheWAS identified six phecodes associated with the PRS with the most significant phenotypes being 218 'benign neoplasm of uterus' and 218.1 'uterine leiomyoma' (p = 1.94 × 10-23, OR 1.31 [95% CI 1.26-1.37] and p = 3.50 × 10-23, OR 1.32 [95% CI 1.26-1.37]). We have developed and validated the first PRS for UF. We find our PRS has predictive ability for UF and captures genetic architecture of increased risk for UF that can be used in further studies.

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.

Blunted PTH response to vitamin D insufficiency/deficiency and colorectal neoplasia risk.

BACKGROUND & AIMS: In contrast to many observational studies, large-scale randomized trials do not support the protective role of vitamin D for the prevention of colorectal neoplasia. However, in previous studies, individuals with blunted parathyroid hormone (PTH) response to vitamin D insufficiency/deficiency (BPRVID), were not differentiated from those with high PTH response to vitamin D insufficiency/deficiency (HPRVID). Individuals with BPRVID are responsive to magnesium treatment, particularly treatment of magnesium plus vitamin D while those with HPRVID are responsive to vitamin D treatment. We prospectively compared these two distinct groups (i.e. BPRVID and HPRVID) for risk of incident adenoma, metachronous adenoma, and incident colorectal cancer (CRC) METHODS: Three nested case-control studies in the Prostate, Lung, Colorectal and Ovarian Cancer (PLCO) Screening Trial. RESULTS: We found optimal 25(OH)D levels were associated with a significantly reduced risk of CRC, primarily among women. The associations between 25(OH)D and CRC risk significantly differed by PTH levels, particularly among women. Compared to individuals with optimal levels for both 25(OH)D and PTH, all others were at an elevated risk of incident CRC, primarily in women. We found those with BPRVID had 2.56-fold significantly increased risk of CRC compared to 1.65-fold non-significantly increased risk for those with HPRVID. Among women, we observed those with BPRVID had 4.79-6.25-fold significantly increased risks of incident CRC and adenoma whereas those with HPRVID had 3.65-fold significantly increased risk of CRC. CONCLUSIONS: Individuals with BPRVID are at higher risks of incident adenoma and CRC compared to those with HPRVID, particularly among women.

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.

Heritability and genome-wide association study of benign prostatic hyperplasia (BPH) in the eMERGE network.

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.

Genetic Risk Score Is Associated With Prevalence of Advanced Neoplasms in a Colorectal Cancer Screening Population.

BACKGROUND & AIMS: The presence of specific single nucleotide polymorphisms (SNPs) can be used to calculate an individual's risk for colorectal cancer (CRC), called a genetic risk score (GRS). We investigated whether GRS can identify individuals with clinically relevant neoplasms in a screening colonoscopy population. METHODS: We derived a GRS based on 48 SNPs associated with CRC, identified in a comprehensive literature search. We obtained genetic data from 1043 participants (50-79 years old) in a screening colonoscopy study in Germany, recruited from 2005 through 2013 (294 with advanced neoplasms, 249 with non-advanced adenoma (NAAs), and 500 without neoplasms). Each participant was assigned a GRS by aggregating their risk alleles (0, 1, or 2). Risk of advanced neoplasms and NAA according to GRS was calculated by multiple logistic regression. Risk advancement periods were calculated. We replicated our findings using data from a subset of the Tennessee Colorectal Polyp Study. RESULTS: An increased GRS was associated with higher prevalence of advanced neoplasms, but not NAAs. Participants in the middle and upper tertiles of GRS had a 2.2-fold and 2.7-fold increase in risk, respectively, of advanced neoplasms compared to those in the lower tertile. Adjusted odds ratios (ORs) were 1.09 (95% confidence interval [CI], 0.76-1.57) for NAA in the middle tertile and 1.05 (95% CI, 0.70-1.55) for NAA in the upper tertile. The ORs were largest for proximal advanced neoplasms for participants in the middle tertile (OR, 3.55; 95% CI 1.85-6.82) and the upper tertile (OR, 3.61; 95% CI 1.84-7.10). The risk advancement period for medium vs low GRS was 13.4 years (95% CI 4.8-22.0) and for high vs low GRS was 17.5 years (95% CI, 7.8-27.3). CONCLUSIONS: In a genetic analysis of participants in a CRC screening study in Germany, an increased GRS (based on CRC-associated SNPs) was associated with increased prevalence of advanced neoplasms. These findings might be used in defining risk-adapted screening ages.

Gene-based evaluation of low-frequency variation and genetically-predicted gene expression impacting risk of keloid formation.

Keloids are benign dermal tumors occurring approximately 20 times more often in individuals of African descent as compared to individuals of European descent. While most keloids occur sporadically, a genetic predisposition is supported by both familial aggregation of some keloids and large differences in risk among populations. Despite Africans and African Americans being at increased risk over lighter-skinned individuals, little genetic research exists into this phenotype. Using a combination of admixture mapping and exome analysis, we reported multiple common variants within chr15q21.2-22.3 associated with risk of keloid formation in African Americans. Here we describe a gene-based association analysis using 478 African American samples with exome genotyping data to identify genes containing low-frequency variants associated with keloids, with evaluation of genetically-predicted gene expression in skin tissues using association summary statistics. The strongest signal from gene-based association was located in C15orf63 (P-value = 6.6 × 10-6 ) located at 15q15.3. The top result from gene expression was increased predicted DCAF4 expression (P-value = 5.5 × 10-4 ) in non-sun-exposed skin, followed by increased predicted OR10A3 expression in sun-exposed skin (P-value = 6.9 × 10-4 ). Our findings identify variation with putative roles in keloid formation, enhanced by the use of predicted gene expression to support the biological roles of variation identified only though genetic association studies.

Erratum to: A multi-stage genome-wide association study of uterine fibroids in African Americans.

The article "A multi-stage genome-wide association study of uterine fibroids in African Americans", written by Jacklyn N. Hellwege, was originally published Online First without open access. After publication in volume 136, issue 10, page 1363-1373 the author decided to opt for Open Choice and to make the article an open access publication. Therefore, the copyright of the article has been changed to

Race/Ethnicity and the Pharmacogenetics of Reported Suicidality With Efavirenz Among Clinical Trials Participants.

Background: We examined associations between suicidality and genotypes that predict plasma efavirenz exposure among AIDS Clinical Trials Group study participants in the United States. Methods: Four clinical trials randomly assigned treatment-naive participants to efavirenz-containing regimens; suicidality was defined as reported suicidal ideation or attempted or completed suicide. Genotypes that predict plasma efavirenz exposure were defined by CYP2B6 and CYP2A6 polymorphisms. Associations were evaluated with weighted Cox proportional hazards models stratified by race/ethnicity. Additional analyses adjusted for genetic ancestry and selected covariates. Results: Among 1833 participants, suicidality was documented in 41 in exposed analyses, and 34 in on-treatment analyses. In unadjusted analyses based on 12 genotype levels, suicidality increased per level in exposed (hazard ratio, 1.11; 95% confidence interval, .96-1.27) and on-treatment 1.16; 1.01-1.34) analyses. In the on-treatment analysis, the association was strongest among white but nearly null among black participants. Considering 3 metabolizer levels (extensive, intermediate and slow), slow metabolizers were at increased risk. Results were similar after baseline covariate-adjustment for genetic ancestry, sex, age, weight, injection drug use history, and psychiatric history or recent psychoactive medication. Conclusions: Genotypes that predict higher plasma efavirenz exposure were associated with increased risk of suicidality. Strength of association varied by race/ethnicity.

Evidence of selection as a cause for racial disparities in fibroproliferative disease.

Fibroproliferative diseases are common complex traits featuring scarring and overgrowth of connective tissue which vary widely in presentation because they affect many organ systems. Most fibroproliferative diseases are more prevalent in African-derived populations than in European populations, leading to pronounced health disparities. It is hypothesized that the increased prevalence of these diseases in African-derived populations is due to selection for pro-fibrotic alleles that are protective against helminth infections. We constructed a genetic risk score (GRS) of fibroproliferative disease risk-increasing alleles using 147 linkage disequilibrium-pruned variants identified through genome-wide association studies of seven fibroproliferative diseases with large African-European prevalence disparities. A comparison of the fibroproliferative disease GRS between 1000 Genomes Phase 3 populations detected a higher mean GRS in AFR (mean = 148 risk alleles) than EUR (mean = 136 risk alleles; T-test p-value = 1.75x10-123). To test whether differences in GRS burden are systematic and may be due to selection, we employed the quantitative trait loci (QTL) sign test. The QTL sign test result indicates that population differences in risk-increasing allele burdens at these fibroproliferative disease variants are systematic and support a model featuring selective pressure (p-value = 0.011). These observations were replicated in an independent sample and were more statistically significant (T-test p-value = 7.26x10-237, sign test p-value = 0.015). This evidence supports the role of selective pressure acting to increase frequency of fibroproliferative alleles in populations of African relative to European ancestry populations.

A multi-stage genome-wide association study of uterine fibroids in African Americans.

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.

Genome-wide linkage and association analysis of cardiometabolic phenotypes in Hispanic Americans.

Linkage studies of complex genetic diseases have been largely replaced by genome-wide association studies, due in part to limited success in complex trait discovery. However, recent interest in rare and low-frequency variants motivates re-examination of family-based methods. In this study, we investigated the performance of two-point linkage analysis for over 1.6 million single-nucleotide polymorphisms (SNPs) combined with single variant association analysis to identify high impact variants, which are both strongly linked and associated with cardiometabolic traits in up to 1414 Hispanics from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Evaluation of all 50 phenotypes yielded 83 557 000 LOD (logarithm of the odds) scores, with 9214 LOD scores ⩾3.0, 845 ⩾4.0 and 89 ⩾5.0, with a maximal LOD score of 6.49 (rs12956744 in the LAMA1 gene for tumor necrosis factor-α (TNFα) receptor 2). Twenty-seven variants were associated with P<0.005 as well as having an LOD score >4, including variants in the NFIB gene under a linkage peak with TNFα receptor 2 levels on chromosome 9. Linkage regions of interest included a broad peak (31 Mb) on chromosome 1q with acute insulin response (max LOD=5.37). This region was previously documented with type 2 diabetes in family-based studies, providing support for the validity of these results. Overall, we have demonstrated the utility of two-point linkage and association in comprehensive genome-wide array-based SNP genotypes.

Genetic variants in selenoprotein P plasma 1 gene (SEPP1) are associated with fasting insulin and first phase insulin response in Hispanics.

CONTEXT: Insulin resistance is not fully explained on a molecular level, though several genes and proteins have been tied to this defect. Knockdowns of the SEPP1 gene, which encodes the selenoprotein P (SeP) protein, have been shown to increase insulin sensitivity in mice. SeP is a liver-derived plasma protein and a major supplier of selenium, which is a proposed insulin mimetic and antidiabetic agent. OBJECTIVE: SEPP1 single nucleotide polymorphisms (SNPs) were selected for analysis with glucometabolic measures. PARTICIPANTS AND MEASURES: The study included1424 Hispanics from families in the Insulin Resistance Atherosclerosis Family Study (IRASFS). Additionally, the multi-ethnic Insulin Resistance Atherosclerosis Study was used. A frequently sampled intravenous glucose tolerance test was used to obtain precise measures of acute insulin response (AIR) and the insulin sensitivity index (SI). DESIGN: 21 SEPP1 SNPs (tagging SNPs (n=12) from HapMap, 4 coding variants and 6 SNPs in the promoter region) were genotyped and analyzed for association. RESULTS: Two highly correlated (r(2)=1) SNPs showed association with AIR (rs28919926; Cys368Arg; p=0.0028 and rs146125471; Ile293Met; p=0.0026) while rs16872779 (intronic) was associated with fasting insulin levels (p=0.0097). In the smaller IRAS Hispanic cohort, few of the associations seen in the IRASFS were replicated, but meta-analysis of IRASFS and all 3 IRAS cohorts (N=2446) supported association of rs28919926 and rs146125471 with AIR (p=0.013 and 0.0047, respectively) as well as rs7579 with SI (p=0.047). CONCLUSIONS: Overall, these results in a human sample are consistent with the literature suggesting a role for SEPP1 in insulin resistance.