Examining interactions between polygenic scores and interpersonal trauma exposure on alcohol consumption and use disorder in an ancestrally diverse college cohort

Introduction: Genetic factors impact alcohol consumption and use disorder (AUD), with large-scale genome-wide association studies (GWAS) identifying numerous associated variants. Aggregate genetic methods in combination with important environmental factors (e.g., interpersonal trauma [IPT]) can be applied to expand our understanding of the ways by which genetic and environmental variables work together to influence alcohol consumption and disordered use. The present study aimed to detail the relationships between genome-wide polygenic scores (PGS) for alcohol phenotypes (i.e., alcohol consumption and AUD status) and IPT exposure as well as the interaction between them across ancestry. Methods: Data were drawn from the Spit for Science (S4S) study, a US college student population, where participants reported on IPT exposure prior to college and alcohol consumption and problems during college (N = 9,006; ancestry: 21.3% African [AFR], 12.5% Admixed Americas [AMR], 9.6% East Asian [EAS], 48.1% European [EUR], 8.6% South Asian [SAS]). Two trans-ancestry PGS were constructed, one for alcohol consumption and another for AUD, using large-scale GWAS summary statistics from multiple ancestries weighted using PRS-CSx. Regression models were applied to test for the presence of associations between alcohol-PGS and IPT main and interaction effects. Results: In the meta-analysis across ancestry groups, IPT exposure and PGS were significantly associated with alcohol consumption (βIPT = 0.31, P IPT = 0.0002; βPGS = 0.09, P PGS = 0.004) and AUD (ORIPT = 1.12, P IPT = 3.5 × 10-8; ORPGS = 1.02, P PGS = 0.002). No statistically significant interactions were detected between IPT and sex nor between IPT and PGS. When inspecting ancestry specific results, the alcohol consumption-PGS and AUD-PGS were only statistically significant in the EUR ancestry group (βPGS = 0.09, P PGS = 0.04; ORPGS = 1.02, P PGS = 0.022, respectively). Discussion: IPT exposure prior to college was strongly associated with alcohol outcomes in this college-age sample, which could be used as a preventative measure to identify students at high risk for problematic alcohol use. Additionally, results add to developing evidence of polygenic score association in meta-analyzed samples, highlighting the importance of continued efforts to increase ancestral representation in genetic studies and inclusive analytic approaches to increase the generalizability of results from genetic association studies.

Bridging the diversity gap: Analytical and study design considerations for improving the accuracy of trans-ancestry genetic prediction

Genetic prediction of common complex disease risk is an essential component of precision medicine. Currently, genome-wide association studies (GWASs) are mostly composed of European-ancestry samples and resulting polygenic scores (PGSs) have been shown to poorly transfer to other ancestries partly due to heterogeneity of allelic effects between populations. Fixed-effects (FETA) and random-effects (RETA) trans-ancestry meta-analyses do not model such ancestry-related heterogeneity, while ancestry-specific (AS) scores may suffer from low power due to low sample sizes. In contrast, trans-ancestry meta-regression (TAMR) builds ancestry-aware PGS that account for more complex trans-ancestry architectures. Here, we examine the predictive performance of these four PGSs under multiple genetic architectures and ancestry configurations. We show that the predictive performance of FETA and RETA is strongly affected by cross-ancestry genetic heterogeneity, while AS PGS performance decreases in under-represented target populations. TAMR PGS is also impacted by heterogeneity but maintains good prediction performance in most situations, especially in ancestry-diverse scenarios. In simulations of human complex traits, TAMR scores currently explain 25% more phenotypic variance than AS in triglyceride levels and 33% more phenotypic variance than FETA in type 2 diabetes in most non-European populations. Importantly, a high proportion of non-European-ancestry individuals is needed to reach prediction levels that are comparable in those populations to the one observed in European-ancestry studies. Our results highlight the need to rebalance the ancestral composition of GWAS to enable accurate prediction in non-European-ancestry groups, and demonstrate the relevance of meta-regression approaches for compensating some of the current population biases in GWAS.

Genome-wide association studies: utility and limitations for research in physiology

Physiological systems are subject to interindividual variation encoded by genetics. Genome-wide association studies (GWAS) operate by surveying thousands of genetic variants from a substantial number of individuals and assessing their association to a trait of interest, be it a physiological variable, a molecular phenotype (e.g. gene expression), or even a disease or condition. Through a myriad of methods, GWAS downstream analyses then explore the functional consequences of each variant and attempt to ascertain a causal relationship to the phenotype of interest, as well as to delve into its links to other traits. This type of investigation allows mechanistic insights into physiological functions, pathological disturbances and shared biological processes between traits (i.e. pleiotropy). An exciting example is the discovery of a new thyroid hormone transporter (SLC17A4) and hormone metabolising enzyme (AADAT) from a GWAS on free thyroxine levels. Therefore, GWAS have substantially contributed with insights into physiology and have been shown to be useful in unveiling the genetic control underlying complex traits and pathological conditions; they will continue to do so with global collaborations and advances in genotyping technology. Finally, the increasing number of trans-ancestry GWAS and initiatives to include ancestry diversity in genomics will boost the power for discoveries, making them also applicable to non-European populations.

Association between antihypertensive drugs and hepatocellular carcinoma: a trans-ancestry and drug-target Mendelian randomization study

BACKGROUND & AIMS

Antihypertensive drugs were recently reported to have an oncogenic role in common cancer, however, whether these drugs would affect the risk of hepatocellular carcinoma (HCC) remains unclear.

METHODS

A drug-target Mendelian randomization method was adopted to examine the long-term effect of 12 antihypertensive drugs classes on the risk of HCC in Europeans and East Asians. To proxy antihypertensive drugs, we leveraged genetic variants located near or within drug target genes that were associated with systolic blood pressure (SBP). Genetically proxied drugs associated with reduced risk of coronary artery disease were included in primary analysis. Genetic summary statistics of SBP and HCC were derived from publicly available large-scale genome-wide association studies in Europeans and East Asians, respectively. Expression quantitative trait loci (eQTLs) of drugs target genes were used to proxy drugs in a sensitivity analysis.

RESULTS

Genetically proxied thiazides and related diuretics were associated with decreased risk of HCC in both Europeans (OR [95% CI]:0.79 [0.73,0.86] per 1mmHg reduction in SBP; P<0.001) and East Asians (0.60 [0.45,0.82]; P=0.001). Genetically proxied beta-adrenoceptor blockers (BBs) were strongly associated with increased risk of HCC in Europeans (1.46 [1.12,1.91]; P=0.004). These findings were replicated in deCODE genetics study and remained consistent when using eQTLs to proxy antihypertensive drugs.

CONCLUSIONS

Our findings suggested that thiazides diuretics may lower the risk of HCC in both Europeans and East Asians, while BBs may increase the risk of HCC specifically in Europeans. Further studies are warranted to explore the potential of repurposing or retargeting antihypertensive drugs for HCC prevention.

Leveraging the local genetic structure for trans-ancestry association mapping

Over the past two decades, genome-wide association studies (GWASs) have successfully advanced our understanding of the genetic basis of complex traits. Despite the fruitful discovery of GWASs, most GWAS samples are collected from European populations, and these GWASs are often criticized for their lack of ancestry diversity. Trans-ancestry association mapping (TRAM) offers an exciting opportunity to fill the gap of disparities in genetic studies between non-Europeans and Europeans. Here, we propose a statistical method, LOG-TRAM, to leverage the local genetic architecture for TRAM. By using biobank-scale datasets, we showed that LOG-TRAM can greatly improve the statistical power of identifying risk variants in under-represented populations while producing well-calibrated p values. We applied LOG-TRAM to the GWAS summary statistics of various complex traits/diseases from BioBank Japan, UK Biobank, and African populations. We obtained substantial gains in power and achieved effective correction of confounding biases in TRAM. Finally, we showed that LOG-TRAM can be successfully applied to identify ancestry-specific loci and the LOG-TRAM output can be further used for construction of more accurate polygenic risk scores in under-represented populations.

Genome-wide Association Studies in Ancestrally Diverse Populations: Opportunities, Methods, Pitfalls, and Recommendations

Genome-wide association studies (GWASs) have focused primarily on populations of European descent, but it is essential that diverse populations become better represented. Increasing diversity among study participants will advance our understanding of genetic architecture in all populations and ensure that genetic research is broadly applicable. To facilitate and promote research in multi-ancestry and admixed cohorts, we outline key methodological considerations and highlight opportunities, challenges, solutions, and areas in need of development. Despite the perception that analyzing genetic data from diverse populations is difficult, it is scientifically and ethically imperative, and there is an expanding analytical toolbox to do it well.

Adipose Tissue Gene Expression Associations Reveal Hundreds of Candidate Genes for Cardiometabolic Traits

Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with cardiometabolic traits including type 2 diabetes (T2D), lipid levels, body fat distribution, and adiposity, although most causal genes remain unknown. We used subcutaneous adipose tissue RNA-seq data from 434 Finnish men from the METSIM study to identify 9,687 primary and 2,785 secondary cis-expression quantitative trait loci (eQTL; <1 Mb from TSS, FDR < 1%). Compared to primary eQTL signals, secondary eQTL signals were located further from transcription start sites, had smaller effect sizes, and were less enriched in adipose tissue regulatory elements compared to primary signals. Among 2,843 cardiometabolic GWAS signals, 262 colocalized by LD and conditional analysis with 318 transcripts as primary and conditionally distinct secondary cis-eQTLs, including some across ancestries. Of cardiometabolic traits examined for adipose tissue eQTL colocalizations, waist-hip ratio (WHR) and circulating lipid traits had the highest percentage of colocalized eQTLs (15% and 14%, respectively). Among alleles associated with increased cardiometabolic GWAS risk, approximately half (53%) were associated with decreased gene expression level. Mediation analyses of colocalized genes and cardiometabolic traits within the 434 individuals provided further evidence that gene expression influences variant-trait associations. These results identify hundreds of candidate genes that may act in adipose tissue to influence cardiometabolic traits.

HLA haplotypes in primary sclerosing cholangitis patients of admixed and non-European ancestry

Primary sclerosing cholangitis (PSC) is strongly associated with several human leukocyte antigen (HLA) haplotypes. Due to extensive linkage disequilibrium and multiple polymorphic candidate genes in the HLA complex, identifying the alleles responsible for these associations has proven difficult. We aimed to evaluate whether studying populations of admixed or non-European descent could help in defining the causative HLA alleles. When assessing haplotypes carrying HLA-DRB1*13:01 (hypothesized to specifically increase the susceptibility to chronic cholangitis), we observed that every haplotype in the Scandinavian PSC population carried HLA-DQB1*06:03. In contrast, only 65% of HLA-DRB1*13:01 haplotypes in an admixed/non-European PSC population carried this allele, suggesting that further assessments of the PSC-associated haplotype HLA-DRB1*13:01-DQA1*01:03-DQB1*06:03 in admixed or multi-ethnic populations could aid in identifying the causative allele.