Genome-wide Association Studies Findings Research Articles

Integrating variant functional annotation scores have varied abilities to improve power of genome-wide association studies

Functional annotations have the potential to increase power of genome-wide association studies (GWAS) by prioritizing variants according to their biological function, but this potential has not been well studied. We comprehensively evaluated all 1132 traits in the UK Biobank whose SNP-heritability estimates were given “medium” or “high” labels by Neale’s lab. For each trait, we integrated GWAS summary statistics of close to 8 million common variants (minor allele frequency >1%) with either their 75 individual functional scores or their meta-scores, using three different data-integration methods. Overall, the number of new genome-wide significant findings after data-integration increases as a trait SNP-heritability estimate increases. However, there is a trade-off between new findings and loss of baseline GWAS findings, resulting in similar total numbers of significant findings between using GWAS alone and integrating GWAS with functional scores, across all 1132 traits analyzed and all three data-integration methods considered. Our findings suggest that, even with the current biobank-level sample size, more informative functional scores and/or new data-integration methods are needed to further improve the power of GWAS of common variants. For example, studying variants in coding sequence and obtaining cell-type-specific scores are potential future directions.

PB1969: FUNCTIONAL CHARACTERIZATION OF SUSCEPTIBILITY LOCI TO MULTIPLE MYELOMA

Background: Multiple myeloma (MM) is the second most common blood malignancy, caused by an uncontrolled growth of plasma cells in the bone marrow, accounting for 20% of all newly diagnosed hematological cancers. Although the current 5-year survival rate is ranging between 40-60%, MM is still considered an incurable disease since most of the patients eventually relapse. While the causes of MM are incompletely understood, several genome-wide association studies (GWAS) have been conducted to identify germline variants that predispose to MM. Up to date, a total of 24 loci were found to be associated with MM risk, but very little information is available about their functional role. Aims: The principal goal is to explore in silico the function of the germline variants associated with MM risk. As we do not know if the GWAS-identified SNPs are the causal variants or just markers of risk, the functional characterization of the causal risk variants would lead to a better understanding of disease development. Methods: GWAS design takes advantage of the linkage disequilibrium (LD) structure of the human genome, thus the main GWAS findings are single-nucleotide polymorphisms (SNPs) that show the strongest association with MM risk (measured as the lowest p-values), but they are not necessarily the functionally causal variants. In this project, we used bioinformatics tools (GTEx, HaploReg v4.1, Roadmap, LDlinke, SNPnexus, RegulomeDB 2.0.3, SNP2TFBS, miRNASNP v3, GeneMANIA) to perform fine mapping of all GWAS-identified loci and to prioritize in each locus the polymorphism with the highest chance of being functionally relevant. In particular, we focused on the loci with the smallest number of SNPs in high LD (r2>0.8) in order to maximize the probability to capture the casual variant. Results: Four of the 24 MM risk loci had a relatively small number of SNPs in high LD and within them we found that the locus located at chromosome 16 contained the greatest number of functionally annotated SNPs. Of particular functional interest was rs3747481 (chr16:30666367 C/T) due to the following reasons: it is a missense variant (protein change: P359L), it is located in the PRR14 gene, that contributes to chromatin hierarchical organization and has a role in gene regulation, has a high CADD PHRED score (22.1). Additionally, according to GTEx portal, rs3747481 is associated with the expression level of the RNF40 gene (which plays a central role in histone code and gene regulation) in whole blood cells (p=3.02-12). It has a score of “1d” in RegulomeDB (eQTL+ TF binding + any motif + DNase peak), meaning that this variant has a high likelihood to affect binding of transcription factors. Some other SNPs (like rs35629860 and rs67128646) in the same LD block show the co-occurrence of H3K4me3 and H3K27me3 histone marks (associated with gene activation and repression, respectively) in promoters and enhancers in B lymphocytes. rs6565197 is predicted to affect the binding of the KLF4 and KLF5 transcription factors which play key roles in cell cycle regulation. Summary/Conclusion: Through a fine mapping of MM risk loci by bioinformatics tools, we found a variant in the locus 16p11.2 that shows in silico a very high probability to have biological role in the risk disease.

Alcohol use disorder is associated with DNA methylation-based shortening of telomere length and regulated by TESPA1: implications for aging

Chronic heavy alcohol consumption is associated with increased mortality and morbidity and often leads to premature aging; however, the mechanisms of alcohol-associated cellular aging are not well understood. In this study, we used DNA methylation derived telomere length (DNAmTL) as a novel approach to investigate the role of alcohol use on the aging process. DNAmTL was estimated by 140 cytosine phosphate guanines (CpG) sites in 372 individuals with alcohol use disorder (AUD) and 243 healthy controls (HC) and assessed using various endophenotypes and clinical biomarkers. Validation in an independent sample of DNAmTL on alcohol consumption was performed (N = 4219). Exploratory genome-wide association studies (GWAS) on DNAmTL were also performed to identify genetic variants contributing to DNAmTL shortening. Top GWAS findings were analyzed using in-silico expression quantitative trait loci analyses and related to structural MRI hippocampus volumes of individuals with AUD. DNAmTL was 0.11-kilobases shorter per year in AUD compared to HC after adjustment for age, sex, race, and blood cell composition (p = 4.0 × 10−12). This association was partially attenuated but remained significant after additionally adjusting for BMI, and smoking status (0.06 kilobases shorter per year, p = 0.002). DNAmTL shortening was strongly associated with chronic heavy alcohol use (ps < 0.001), elevated gamma-glutamyl transferase (GGT), and aspartate aminotransferase (AST) (ps < 0.004). Comparison of DNAmTL with PCR-based methods of assessing TL revealed positive correlations (R = 0.3, p = 2.2 × 10−5), highlighting the accuracy of DNAmTL as a biomarker. The GWAS meta-analysis identified a single nucleotide polymorphism (SNP), rs4374022 and 18 imputed ones in Thymocyte Expressed, Positive Selection Associated 1(TESPA1), at the genome-wide level (p = 3.75 × 10−8). The allele C of rs4374022 was associated with DNAmTL shortening, lower hippocampus volume (p < 0.01), and decreased mRNA expression in hippocampus tissue (p = 0.04). Our study demonstrates DNAmTL-related aging acceleration in AUD and suggests a functional role for TESPA1 in regulating DNAmTL length, possibly via the immune system with subsequent biological effects on brain regions negatively affected by alcohol and implicated in aging.

Abstract 1556: Long non-coding RNA WNT5A-AS1 shows sex-dimorphic differences in survival for males with glioblastoma

Abstract Glioblastoma (GBM) is the most common and aggressive primary brain tumor with a median survival time of 12-15 months and 5-year survival rate of 5%. Despite the application of advanced genetics and biological breakthroughs in GBM, the outcome for GBM remains dismal. Previous genome-wide association studies (GWAS) for GBM shows that the incidence rate for GBM in males is 1.6 times higher in men than women, with men displaying overall worse outcomes. In-vivo studies in gender-matched isogenic murine astrocytic cell lines show that the rate of astrocytic transformation is higher in males than in females, corroborating the GWAS findings in human GBM. Thus, understanding how sex effects GBM biology may provide insight into the mechanism of gliomagenesis and response to therapy. To explore the molecular mechanisms by which transformation of GBM increases in males, we used DESeq2 to assess differentially expressed genes between male and female patients. We found that a long non-coding RNA WNT5A-AS1 is overexpressed in males [log2(fold change) of 1.40; FDR = 0.036], regardless of cofounding factors such as IDH1 status, age, and race. Next, we examined if expression of WNT5A-AS1 is correlated with survival within male patients versus female patients in the TCGA-GBM dataset. In males, WNT5A-AS1’s expression is significantly correlated with patient overall survival (P = 0.042, N=109), but not in female (P = 0.585, N=60). Male patients with high expression of WNT5A-AS1 display poor survival outcomes compared to males with low expression. WNT5A-AS1 is located in the upstream, promoter region of WNT5A. Q-PCR quantification of WNT5A expression increases in male GBM specimens as compared to GBM specimens from females. The role that WNT5A-AS1 plays in regulating WNT5A expression is unclear to date. Interestingly, WNT5A expression shows a positive correlation with patient survival in males, with high expression displaying better survival outcomes (p=0.0272). We have data obtained from a subset of patient samples showing an inverse relationship between WNT5A and WNT5A-AS1 expression. This suggests that there may be an inverse correlation between WNT5A-AS1 and WNT5A expression, which elucidates that WNT5A-AS1 may be regulating WNT5A expression. Current ongoing studies will evaluate the role of WNT5A-AS1/WNT5A in sex-differences of GBM progression. Citation Format: Aunay Miller, Junwen Wang, Joseph Loftus, Nhan L. Tran. Long non-coding RNA WNT5A-AS1 shows sex-dimorphic differences in survival for males with glioblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1556.

Improving polygenic prediction in ancestrally diverse populations.

Polygenic risk scores (PRS) have attenuated cross-population predictive performance. As existing genome-wide association studies (GWAS) have been conducted predominantly in individuals of European descent, the limited transferability of PRS reduces their clinical value in non-European populations, and may exacerbate healthcare disparities. Recent efforts to level ancestry imbalance in genomic research have expanded the scale of non-European GWAS, although most remain underpowered. Here, we present a new PRS construction method, PRS-CSx, which improves cross-population polygenic prediction by integrating GWAS summary statistics from multiple populations. PRS-CSx couples genetic effects across populations via a shared continuous shrinkage (CS) prior, enabling more accurate effect size estimation by sharing information between summary statistics and leveraging linkage disequilibrium diversity across discovery samples, while inheriting computational efficiency and robustness from PRS-CS. We show that PRS-CSx outperforms alternative methods across traits with a wide range of genetic architectures, cross-population genetic overlaps and discovery GWAS sample sizes in simulations, and improves the prediction of quantitative traits and schizophrenia risk in non-European populations.

Mapping anorexia nervosa genes to clinical phenotypes.

Anorexia nervosa (AN) is a psychiatric disorder with complex etiology, with a significant portion of disease risk imparted by genetics. Traditional genome-wide association studies (GWAS) produce principal evidence for the association of genetic variants with disease. Transcriptomic imputation (TI) allows for the translation of those variants into regulatory mechanisms, which can then be used to assess the functional outcome of genetically regulated gene expression (GReX) in a broader setting through the use of phenome-wide association studies (pheWASs) in large and diverse clinical biobank populations with electronic health record phenotypes. Here, we applied TI using S-PrediXcan to translate the most recent PGC-ED AN GWAS findings into AN-GReX. For significant genes, we imputed AN-GReX in the Mount Sinai BioMe™ Biobank and performed pheWASs on over 2000 outcomes to test the clinical consequences of aberrant expression of these genes. We performed a secondary analysis to assess the impact of body mass index (BMI) and sex on AN-GReX clinical associations. Our S-PrediXcan analysis identified 53 genes associated with AN, including what is, to our knowledge, the first-genetic association of AN with the major histocompatibility complex. AN-GReX was associated with autoimmune, metabolic, and gastrointestinal diagnoses in our biobank cohort, as well as measures of cholesterol, medications, substance use, and pain. Additionally, our analyses showed moderation of AN-GReX associations with measures of cholesterol and substance use by BMI, and moderation of AN-GReX associations with celiac disease by sex. Our BMI-stratified results provide potential avenues of functional mechanism for AN-genes to investigate further.

Phenotypes, mechanisms and therapeutics: insights from bipolar disorder GWAS findings

Genome-wide association studies (GWAS) have reported substantial genomic loci significantly associated with clinical risk of bipolar disorder (BD), and studies combining techniques of genetics, neuroscience, neuroimaging, and pharmacology are believed to help tackle clinical problems (e.g., identifying novel therapeutic targets). However, translating findings of psychiatric genetics into biological mechanisms underlying BD pathogenesis remains less successful. Biological impacts of majority of BD GWAS risk loci are obscure, and the involvement of many GWAS risk genes in this illness is yet to be investigated. It is thus necessary to review the progress of applying BD GWAS risk genes in the research and intervention of the disorder. A comprehensive literature search found that a number of such risk genes had been investigated in cellular or animal models, even before they were highlighted in BD GWAS. Intriguingly, manipulation of many BD risk genes (e.g., ANK3, CACNA1C, CACNA1B, HOMER1, KCNB1, MCHR1, NCAN, SHANK2 etc.) resulted in altered murine behaviors largely restoring BD clinical manifestations, including mania-like symptoms such as hyperactivity, anxiolytic-like behavior, as well as antidepressant-like behavior, and these abnormalities could be attenuated by mood stabilizers. In addition to recapitulating phenotypic characteristics of BD, some GWAS risk genes further provided clues for the neurobiology of this illness, such as aberrant activation and functional connectivity of brain areas in the limbic system, and modulated dendritic spine morphogenesis as well as synaptic plasticity and transmission. Therefore, BD GWAS risk genes are undoubtedly pivotal resources for modeling this illness, and might be translational therapeutic targets in the future clinical management of BD. We discuss both promising prospects and cautions in utilizing the bulk of useful resources generated by GWAS studies. Systematic integrations of findings from genetic and neuroscience studies are called for to promote our understanding and intervention of BD.

A genome-wide association study in autoimmune neurological syndromes with anti-GAD65 autoantibodies.

Autoimmune neurological syndromes (AINS) with autoantibodies against the 65 kDa isoform of the glutamic acid decarboxylase (GAD65) present with limbic encephalitis, including temporal lobe seizures or epilepsy, cerebellitis with ataxia, and stiff-person-syndrome or overlap forms. Anti-GAD65 autoantibodies are also detected in autoimmune diabetes mellitus, which has a strong genetic susceptibility conferred by human leukocyte antigen (HLA) and non-HLA genomic regions. We investigated the genetic predisposition in patients with anti-GAD65 AINS. We performed a genome-wide association study (GWAS) and an association analysis of the HLA region in a large German cohort of 1214 individuals. These included 167 patients with anti-GAD65 AINS, recruited by the German Network for Research on Autoimmune Encephalitis (GENERATE), and 1047 individuals without neurological or endocrine disease as population-based controls. Predictions of protein expression changes based on GWAS findings were further explored and validated in the CSF proteome of a virtually independent cohort of 10 patients with GAD65-AINS and 10 controls. Our GWAS identified 16 genome-wide significant (P < 5 × 10-8) loci for the susceptibility to anti-GAD65 AINS. The top variant, rs2535288 [P = 4.42 × 10-16, odds ratio (OR) = 0.26, 95% confidence interval (CI) = 0.187-0.358], localized to an intergenic segment in the middle of the HLA class I region. The great majority of variants in these loci (>90%) mapped to non-coding regions of the genome. Over 40% of the variants have known regulatory functions on the expression of 48 genes in disease relevant cells and tissues, mainly CD4+ T cells and the cerebral cortex. The annotation of epigenomic marks suggested specificity for neural and immune cells. A network analysis of the implicated protein-coding genes highlighted the role of protein kinase C beta (PRKCB) and identified an enrichment of numerous biological pathways participating in immunity and neural function. Analysis of the classical HLA alleles and haplotypes showed no genome-wide significant associations. The strongest associations were found for the DQA1*03:01-DQB1*03:02-DRB1*04:01HLA haplotype (P = 4.39 × 10-4, OR = 2.5, 95%CI = 1.499-4.157) and DRB1*04:01 allele (P = 8.3 × 10-5, OR = 2.4, 95%CI = 1.548-3.682) identified in our cohort. As predicted, the CSF proteome showed differential levels of five proteins (HLA-A/B, C4A, ATG4D and NEO1) of expression quantitative trait loci genes from our GWAS in the CSF proteome of anti-GAD65 AINS. These findings suggest a strong genetic predisposition with direct functional implications for immunity and neural function in anti-GAD65 AINS, mainly conferred by genomic regions outside the classical HLA alleles.

Health influenced by genetics: A first comprehensive analysis of breast cancer high and moderate penetrance susceptibility genes in the Tunisian population.

Significant advances have been made to understand the genetic basis of breast cancer. High, moderate and low penetrance variants have been identified with inter-ethnic variability in mutation frequency and spectrum. Genome wide association studies (GWAS) are widely used to identify disease-associated SNPs. Understanding the functional impact of these risk-SNPs will help the translation of GWAS findings into clinical interventions. Here we aim to characterize the genetic patterns of high and moderate penetrance breast cancer susceptibility genes and to assess the functional impact of non-coding SNPs. We analyzed BRCA1/2, PTEN, STK11, TP53, ATM, BRIP1, CHEK2 and PALB2 genotype data obtained from 135 healthy participants genotyped using Affymetrix Genome-Wide Human SNP-Array 6.0. Haplotype analysis was performed using Haploview.V4.2 and PHASE.V2.1. Population structure and genetic differentiation were assessed using principal component analysis (PCA) and fixation index (FST). Functional annotation was performed using In Silico web-based tools including RegulomeDB and VARAdb. Haplotype analysis showed distinct LD patterns with high levels of recombination and haplotype blocks of moderate to small size. Our findings revealed also that the Tunisian population tends to have a mixed origin with European, South Asian and Mexican footprints. Functional annotation allowed the selection of 28 putative regulatory variants. Of special interest were BRCA1_ rs8176318 predicted to alter the binding sites of a tumor suppressor miRNA hsa-miR-149 and PALB2_ rs120963 located in tumorigenesis-associated enhancer and predicted to strongly affect the binding of P53. Significant differences in allele frequencies were observed with populations of African and European ancestries for rs8176318 and rs120963 respectively. Our findings will help to better understand the genetic basis of breast cancer by guiding upcoming genome wide studies in the Tunisian population. Putative functional SNPs may be used to develop an efficient polygenic risk score to predict breast cancer risk leading to better disease prevention and management.

Genome-Wide Association Study of Alzheimer's Disease Brain Imaging Biomarkers and Neuropsychological Phenotypes in the European Medical Information Framework for Alzheimer's Disease Multimodal Biomarker Discovery Dataset.

Alzheimer’s disease (AD) is the most frequent neurodegenerative disease with an increasing prevalence in industrialized, aging populations. AD susceptibility has an established genetic basis which has been the focus of a large number of genome-wide association studies (GWAS) published over the last decade. Most of these GWAS used dichotomized clinical diagnostic status, i.e., case vs. control classification, as outcome phenotypes, without the use of biomarkers. An alternative and potentially more powerful study design is afforded by using quantitative AD-related phenotypes as GWAS outcome traits, an analysis paradigm that we followed in this work. Specifically, we utilized genotype and phenotype data from n = 931 individuals collected under the auspices of the European Medical Information Framework for Alzheimer’s Disease Multimodal Biomarker Discovery (EMIF-AD MBD) study to perform a total of 19 separate GWAS analyses. As outcomes we used five magnetic resonance imaging (MRI) traits and seven cognitive performance traits. For the latter, longitudinal data from at least two timepoints were available in addition to cross-sectional assessments at baseline. Our GWAS analyses revealed several genome-wide significant associations for the neuropsychological performance measures, in particular those assayed longitudinally. Among the most noteworthy signals were associations in or near EHBP1 (EH domain binding protein 1; on chromosome 2p15) and CEP112 (centrosomal protein 112; 17q24.1) with delayed recall as well as SMOC2 (SPARC related modular calcium binding 2; 6p27) with immediate recall in a memory performance test. On the X chromosome, which is often excluded in other GWAS, we identified a genome-wide significant signal near IL1RAPL1 (interleukin 1 receptor accessory protein like 1; Xp21.3). While polygenic score (PGS) analyses showed the expected strong associations with SNPs highlighted in relevant previous GWAS on hippocampal volume and cognitive function, they did not show noteworthy associations with recent AD risk GWAS findings. In summary, our study highlights the power of using quantitative endophenotypes as outcome traits in AD-related GWAS analyses and nominates several new loci not previously implicated in cognitive decline.

Transcriptome‐wide association and prediction for carotenoids and tocochromanols in fresh sweet corn kernels

Sweet corn (Zea mays L.) is consistently one of the most highly consumed vegetables in the United States, providing a valuable opportunity to increase nutrient intake through biofortification. Significant variation for carotenoid (provitamin A, lutein, zeaxanthin) and tocochromanol (vitamin E, antioxidants) levels is present in temperate sweet corn germplasm, yet previous genome-wide association studies (GWAS) of these traits have been limited by low statistical power and mapping resolution. Here, we employed a high-quality transcriptomic dataset collected from fresh sweet corn kernels to conduct transcriptome-wide association studies (TWAS) and transcriptome prediction studies for 39 carotenoid and tocochromanol traits. In agreement with previous GWAS findings, TWAS detected significant associations for four causal genes, β-carotene hydroxylase (crtRB1), lycopene epsilon cyclase (lcyE), γ-tocopherol methyltransferase (vte4), and homogentisate geranylgeranyltransferase (hggt1) on a transcriptome-wide level. Pathway-level analysis revealed additional associations for deoxy-xylulose synthase2 (dxs2), diphosphocytidyl methyl erythritol synthase2 (dmes2), cytidine methyl kinase1 (cmk1), and geranylgeranyl hydrogenase1 (ggh1), of which, dmes2, cmk1, and ggh1 have not previously been identified through maize association studies. Evaluation of prediction models incorporating genome-wide markers and transcriptome-wide abundances revealed a trait-dependent benefit to the inclusion of both genomic and transcriptomic data over solely genomic data, but both transcriptome- and genome-wide datasets outperformed a priori candidate gene-targeted prediction models for most traits. Altogether, this study represents an important step toward understanding the role of regulatory variation in the accumulation of vitamins in fresh sweet corn kernels.

10 Years of GWAS discovery in endometrial cancer: Aetiology, function and translation.

SummaryEndometrial cancer is a common gynaecological cancer with increasing incidence and mortality. In the last decade, endometrial cancer genome-wide association studies (GWAS) have provided a resource to explore aetiology and for functional interpretation of heritable risk variation, informing endometrial cancer biology. Indeed, GWAS data have been used to assess relationships with other traits through correlation and Mendelian randomisation analyses, establishing genetic relationships and potential risk factors. Cross-trait GWAS analyses have increased statistical power and identified novel endometrial cancer risk variation related to other traits. Functional analysis of risk loci has helped prioritise candidate susceptibility genes, revealing molecular mechanisms and networks. Lastly, risk scores generated using endometrial cancer GWAS data may allow for clinical translation through identification of patients at high risk of disease. In the next decade, this knowledge base should enable substantial progress in our understanding of endometrial cancer and, potentially, new approaches for its screening and treatment.

Genetic Loci Associated With COVID-19 Positivity and Hospitalization in White, Black, and Hispanic Veterans of the VA Million Veteran Program.

SARS-CoV-2 has caused symptomatic COVID-19 and widespread death across the globe. We sought to determine genetic variants contributing to COVID-19 susceptibility and hospitalization in a large biobank linked to a national United States health system. We identified 19,168 (3.7%) lab-confirmed COVID-19 cases among Million Veteran Program participants between March 1, 2020, and February 2, 2021, including 11,778 Whites, 4,893 Blacks, and 2,497 Hispanics. A multi-population genome-wide association study (GWAS) for COVID-19 outcomes identified four independent genetic variants (rs8176719, rs73062389, rs60870724, and rs73910904) contributing to COVID-19 positivity, including one novel locus found exclusively among Hispanics. We replicated eight of nine previously reported genetic associations at an alpha of 0.05 in at least one population-specific or the multi-population meta-analysis for one of the four MVP COVID-19 outcomes. We used rs8176719 and three additional variants to accurately infer ABO blood types. We found that A, AB, and B blood types were associated with testing positive for COVID-19 compared with O blood type with the highest risk for the A blood group. We did not observe any genome-wide significant associations for COVID-19 severity outcomes among those testing positive. Our study replicates prior GWAS findings associated with testing positive for COVID-19 among mostly White samples and extends findings at three loci to Black and Hispanic individuals. We also report a new locus among Hispanics requiring further investigation. These findings may aid in the identification of novel therapeutic agents to decrease the morbidity and mortality of COVID-19 across all major ancestral populations.

Evaluating the interactive effects of dietary habits and human gut microbiome on the risks of depression and anxiety.

Gut microbiome and dietary patterns have been suggested to be associated with depression/anxiety. However, limited effort has been made to explore the effects of possible interactions between diet and microbiome on the risks of depression and anxiety. Using the latest genome-wide association studies findings in gut microbiome and dietary habits, polygenic risk scores (PRSs) analysis of gut microbiome and dietary habits was conducted in the UK Biobank cohort. Logistic/linear regression models were applied for evaluating the associations for gut microbiome-PRS, dietary habits-PRS, and their interactions with depression/anxiety status and Patient Health Questionnaire (PHQ-9)/Generalized Anxiety Disorder-7 (GAD-7) score by R software. We observed 51 common diet-gut microbiome interactions shared by both PHQ score and depression status, such as overall beef intake × genus Sporobacter [hurdle binary (HB)] (PPHQ = 7.88 × 10-4, Pdepression status = 5.86 × 10-4); carbohydrate × genus Lactococcus (HB) (PPHQ = 0.0295, Pdepression status = 0.0150). We detected 41 common diet-gut microbiome interactions shared by GAD score and anxiety status, such as sugar × genus Parasutterella (rank normal transformed) (PGAD = 5.15 × 10-3, Panxiety status = 0.0347); tablespoons of raw vegetables per day × family Coriobacteriaceae (HB) (PGAD = 6.02 × 10-4, Panxiety status = 0.0345). Some common significant interactions shared by depression and anxiety were identified, such as overall beef intake × genus Sporobacter (HB). Our study results expanded our understanding of how to comprehensively consider the relationships for dietary habits-gut microbiome interactions with depression and anxiety.

Stability of polygenic scores across discovery genome-wide association studies.

SummaryPolygenic scores (PGS) are commonly evaluated in terms of their predictive accuracy at the population level by the proportion of phenotypic variance they explain. To be useful for precision medicine applications, they also need to be evaluated at the individual level when phenotypes are not necessarily already known. We investigated the stability of PGS in European American (EUR) and African American (AFR)-ancestry individuals from the Philadelphia Neurodevelopmental Cohort and the Adolescent Brain Cognitive Development study using different discovery genome-wide association study (GWAS) results for post-traumatic stress disorder (PTSD), type 2 diabetes (T2D), and height. We found that pairs of EUR-ancestry GWAS for the same trait had genetic correlations >0.92. However, PGS calculated from pairs of same-ancestry and different-ancestry GWAS had correlations that ranged from <0.01 to 0.74. PGS stability was greater for height than for PTSD or T2D. A series of height GWAS in the UK Biobank suggested that correlation between PGS is strongly dependent on the extent of sample overlap between the discovery GWAS. Focusing on the upper end of the PGS distribution, different discovery GWAS do not consistently identify the same individuals in the upper quantiles, with the best case being 60% of individuals above the 80th percentile of PGS overlapping from one height GWAS to another. The degree of overlap decreases sharply as higher quantiles, less heritable traits, and different-ancestry GWAS are considered. PGS computed from different discovery GWAS have only modest correlation at the individual level, underscoring the need to proceed cautiously with integrating PGS into precision medicine applications.

Functional characterisation of the amyotrophic lateral sclerosis risk locus GPX3/TNIP1

BackgroundAmyotrophic lateral sclerosis (ALS) is a complex, late-onset, neurodegenerative disease with a genetic contribution to disease liability. Genome-wide association studies (GWAS) have identified ten risk loci to date, including the TNIP1/GPX3 locus on chromosome five. Given association analysis data alone cannot determine the most plausible risk gene for this locus, we undertook a comprehensive suite of in silico, in vivo and in vitro studies to address this.MethodsThe Functional Mapping and Annotation (FUMA) pipeline and five tools (conditional and joint analysis (GCTA-COJO), Stratified Linkage Disequilibrium Score Regression (S-LDSC), Polygenic Priority Scoring (PoPS), Summary-based Mendelian Randomisation (SMR-HEIDI) and transcriptome-wide association study (TWAS) analyses) were used to perform bioinformatic integration of GWAS data (Ncases = 20,806, Ncontrols = 59,804) with ‘omics reference datasets including the blood (eQTLgen consortium N = 31,684) and brain (N = 2581). This was followed up by specific expression studies in ALS case-control cohorts (microarray Ntotal = 942, protein Ntotal = 300) and gene knockdown (KD) studies of human neuronal iPSC cells and zebrafish-morpholinos (MO).ResultsSMR analyses implicated both TNIP1 and GPX3 (p < 1.15 × 10−6), but there was no simple SNP/expression relationship. Integrating multiple datasets using PoPS supported GPX3 but not TNIP1. In vivo expression analyses from blood in ALS cases identified that lower GPX3 expression correlated with a more progressed disease (ALS functional rating score, p = 5.5 × 10−3, adjusted R2 = 0.042, Beffect = 27.4 ± 13.3 ng/ml/ALSFRS unit) with microarray and protein data suggesting lower expression with risk allele (recessive model p = 0.06, p = 0.02 respectively). Validation in vivo indicated gpx3 KD caused significant motor deficits in zebrafish-MO (mean difference vs. control ± 95% CI, vs. control, swim distance = 112 ± 28 mm, time = 1.29 ± 0.59 s, speed = 32.0 ± 2.53 mm/s, respectively, p for all < 0.0001), which were rescued with gpx3 expression, with no phenotype identified with tnip1 KD or gpx3 overexpression.ConclusionsThese results support GPX3 as a lead ALS risk gene in this locus, with more data needed to confirm/reject a role for TNIP1. This has implications for understanding disease mechanisms (GPX3 acts in the same pathway as SOD1, a well-established ALS-associated gene) and identifying new therapeutic approaches. Few previous examples of in-depth investigations of risk loci in ALS exist and a similar approach could be applied to investigate future expected GWAS findings.

SNP characteristics and validation success in genome wide association studies.

Genome wide association studies (GWASs) have identified tens of thousands of single nucleotide polymorphisms (SNPs) associated with human diseases and characteristics. A significant fraction of GWAS findings can be false positives. The gold standard for true positives is an independent validation. The goal of this study was to identify SNP features associated with validation success. Summary statistics from the Catalog of Published GWASs were used in the analysis. Since our goal was an analysis of reproducibility, we focused on the diseases/phenotypes targeted by at least 10 GWASs. GWASs were arranged in discovery-validation pairs based on the time of publication, with the discovery GWAS published before validation. We used four definitions of the validation success that differ by stringency. Associations of SNP features with validation success were consistent across the definitions. The strongest predictor of SNP validation was the level of statistical significance in the discovery GWAS. The magnitude of the effect size was associated with validation success in a non-linear manner. SNPs with risk allele frequencies in the range 30-70% showed a higher validation success rate compared to rarer or more common SNPs. Missense, 5'UTR, stop gained, and SNPs located in transcription factor binding sites had a higher validation success rate compared to intergenic, intronic and synonymous SNPs. There was a positive association between validation success and the level of evolutionary conservation of the sites. In addition, validation success was higher when discovery and validation GWASs targeted the same ethnicity. All predictors of validation success remained significant in a multivariate logistic regression model indicating their independent contribution. To conclude, we identified SNP features predicting validation success of GWAS hits. These features can be used to select SNPs for validation and downstream functional studies.

Genetic Variation and Mendelian Randomization Approaches.

While genome-wide association studies (GWAS) on levels of nuclear receptors are sparse, the genetics of ligands of these receptors (steroid hormones, thyroid hormones, and liposoluble vitamins) have been extensively studied in GWAS of predominantly European populations. Hundreds of genetic variants across the genome have been associated with serum levels of nuclear receptor ligands, shedding light on the physiology of hormone metabolism. These GWAS findings have been used to explore causal associations of these hormones with complex human traits and diseases in Mendelian randomization (MR) studies, and in studies using polygenic risk scores to quantify the genetic predisposition to higher/lower hormone levels. As such, besides providing insights into hormonal pathophysiology and its causal relationship with clinical complications, GWAS-identified genetic markers could ultimately play an important role in the daily clinical management of patients. As large trans-ethnic GWAS on levels of nuclear receptor ligands emerge, and with the fast advances in genotyping techniques and constant decrease of the genotyping costs, studying an individual's genetically predicted hormonal profile could be the next step in personalizing the management of patients with pathologies related to nuclear receptors and their ligands.

Transcriptome-wide association study: Opportunity and challenges for cancer studies

Genome-wide association studies (GWAS) have uncovered thousands of single nucleotide polymorphism (SNP) loci that are associated with complex traits. However, the majority of GWAS discoveries are located in non-coding regions and the biological mechanisms behind these associations are not well understood. Transcriptome-wide association studies (TWAS) have gained popularity in recent years by generating biological interpretable discoveries and facilitating the identification of novel associations that have been missed by GWAS. TWAS has identified more than hundreds of susceptibility genes for many complex diseases and traits, including cancers. Here, in this review, we first summarize TWAS methods, then discuss the opportunities for cancer studies and finally review current challenges and future directions for this method.

Transcriptomics insights into interpreting AMD-GWAS discoveries for biological and clinical applications

Genome-wide association studies (GWAS) have been successful in identifying genetic risk factors for a large number of complex diseases, including age-related macular degeneration (AMD), which is a highly heritable complex disease affecting millions of elderly individuals. However, the progress of elucidating the functional relevance of genetic findings in AMD has been slow, as most risk factors are non-coding, and we have little insight into the causal genes and disease mechanisms. In the last few years, gene expression regulation is emerging as a dominant mechanism through which GWAS risk variants lead to the disease. The purpose of this review is to provide an overview of how transcriptome studies can help in identifying the genes, pathways and therapeutic targets underlying GWAS discoveries in AMD. These approaches help pave the road for mechanistic understanding of GWAS findings and drive translational advances that will lead to improved AMD management and treatment.