Across-cohort QC analyses of GWAS summary statistics from complex traits

Decoration of Fe3O4-vitamin C nanoparticles on alginate-chitosan nanocomplex: Characterization, safety, bioacessibility boost and Iron Nanofortification in A2 goat milk gels

In this study, an alginate-chitosan (AL-CS) nanocomplex decorated with vitamin C coated iron oxide nanoparticles (Fe3O4-vit C NPs) was investigated as a novel nanoiron fortification agent. The Fe3O4-vit C NPs decorated on AL-CS nanocomplex underwent comprehensive characterization, including zeta potential, fourier transform infrared spectroscopy, X-ray diffraction, and UV-vis spectroscopy. The transmission electron microscopy (TEM) analysis confirmed the decoration of Fe3O4-vit C NPs on AL-CS nanocomplex. The dynamic light scattering and thermogravimetric analysis showed enhanced thermal properties of decorated nanocomplex than the undecorated control. Biocompatibility testing on HepG2 cell lines revealed improved compatibility, while intestinal Caco2 cell lines showed approximately 51 % greater bioacessibility than controls. Further, 8 mg of Fe3O4-vit C NPs decorated AL-CS nanocomplex nanofortified 80 g of A2 goat milk gels (GMGs) which provided 0.072 mg/g of nanoiron without showing significant changes in texture and color compared to the control A2 GMGs. The PCA analysis helped to identify the impact of various factors for the preparation of decorated nanocomplex.

Shared genetic architecture and bidirectional clinical risks within the psycho-metabolic nexus

BACKGROUND

Increasing evidence suggests a complex interplay between psychiatric disorders and metabolic dysregulations. However, most research has been limited to specific disorder pairs, leaving a significant gap in our understanding of the broader psycho-metabolic nexus.

METHODS

This study leveraged large-scale cohort data and genome-wide association study (GWAS) summary statistics, covering 8 common psychiatric disorders and 43 metabolic traits. We introduced a comprehensive analytical strategy to identify shared genetic bases sequentially, from key genetic correlation regions to local pleiotropy and pleiotropic genes. Finally, we developed polygenic risk score (PRS) models to translate these findings into clinical applications.

FINDINGS

We identified significant bidirectional clinical risks between psychiatric disorders and metabolic dysregulations among 310,848 participants from the UK Biobank. Genetic correlation analysis confirmed 104 robust trait pairs, revealing 1088 key genomic regions, including critical hotspots such as chr3: 47588462-50387742. Cross-trait meta-analysis uncovered 388 pleiotropic single nucleotide variants (SNVs) and 126 shared causal variants. Among variants, 45 novel SNVs were associated with psychiatric disorders and 75 novel SNVs were associated with metabolic traits, shedding light on new targets to unravel the mechanism of comorbidity. Notably, RBM6, a gene involved in alternative splicing and cellular stress response regulation, emerged as a key pleiotropic gene. When psychiatric and metabolic genetic information were integrated, PRS models demonstrated enhanced predictive power.

INTERPRETATION

The study highlights the intertwined genetic and clinical relationships between psychiatric disorders and metabolic dysregulations, emphasising the need for integrated approaches in diagnosis and treatment.

FUNDING

The National Key Research and Development Program of China (2023YFC2506200, SHH). The National Natural Science Foundation of China (82273741, SY).

Two-stage strategy using denoising autoencoders for robust reference-free genotype imputation with missing input genotypes

Widely used genotype imputation methods are based on the Li and Stephens model, which assumes that new haplotypes can be represented by modifying existing haplotypes in a reference panel through mutations and recombinations. These methods use genotypes from SNP arrays as inputs to estimate haplotypes that align with the input genotypes by analyzing recombination patterns within a reference panel, and then infer unobserved variants. While these methods require reference panels in an identifiable form, their public use is limited due to privacy and consent concerns. One strategy to overcome these limitations is to use de-identified haplotype information, such as summary statistics or model parameters. Advances in deep learning (DL) offer the potential to develop imputation methods that use haplotype information in a reference-free manner by handling it as model parameters, while maintaining comparable imputation accuracy to methods based on the Li and Stephens model. Here, we provide a brief introduction to DL-based reference-free genotype imputation methods, including RNN-IMP, developed by our research group. We then evaluate the performance of RNN-IMP against widely-used Li and Stephens model-based imputation methods in terms of accuracy (R2), using the 1000 Genomes Project Phase 3 dataset and corresponding simulated Omni2.5 SNP genotype data. Although RNN-IMP is sensitive to missing values in input genotypes, we propose a two-stage imputation strategy: missing genotypes are first imputed using denoising autoencoders; RNN-IMP then processes these imputed genotypes. This approach restores the imputation accuracy that is degraded by missing values, enhancing the practical use of RNN-IMP.

The goldmine of GWAS summary statistics: a systematic review of methods and tools

Genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of complex traits and diseases. GWAS summary statistics have become essential tools for various genetic analyses, including meta-analysis, fine-mapping, and risk prediction. However, the increasing number of GWAS summary statistics and the diversity of software tools available for their analysis can make it challenging for researchers to select the most appropriate tools for their specific needs. This systematic review aims to provide a comprehensive overview of the currently available software tools and databases for GWAS summary statistics analysis. We conducted a comprehensive literature search to identify relevant software tools and databases. We categorized the tools and databases by their functionality, including data management, quality control, single-trait analysis, and multiple-trait analysis. We also compared the tools and databases based on their features, limitations, and user-friendliness. Our review identified a total of 305 functioning software tools and databases dedicated to GWAS summary statistics, each with unique strengths and limitations. We provide descriptions of the key features of each tool and database, including their input/output formats, data types, and computational requirements. We also discuss the overall usability and applicability of each tool for different research scenarios. This comprehensive review will serve as a valuable resource for researchers who are interested in using GWAS summary statistics to investigate the genetic basis of complex traits and diseases. By providing a detailed overview of the available tools and databases, we aim to facilitate informed tool selection and maximize the effectiveness of GWAS summary statistics analysis.

Polygenic risk score for acute rejection based on donor-recipient non-HLA genotype mismatch

BACKGROUND

Acute rejection (AR) after kidney transplantation is an important allograft complication. To reduce the risk of post-transplant AR, determination of kidney transplant donor-recipient mismatching focuses on blood type and human leukocyte antigens (HLA), while it remains unclear whether non-HLA genetic mismatching is related to post-transplant complications.

METHODS

We carried out a genome-wide scan (HLA and non-HLA regions) on AR with a large kidney transplant cohort of 784 living donor-recipient pairs of European ancestry. An AR polygenic risk score (PRS) was constructed with the non-HLA single nucleotide polymorphisms (SNPs) filtered by independence (r2 < 0.2) and P-value (< 1×10-3) criteria. The PRS was validated in an independent cohort of 352 living donor-recipient pairs.

RESULTS

By the genome-wide scan, we identified one significant SNP rs6749137 with HR = 2.49 and P-value = 2.15×10-8. 1,307 non-HLA PRS SNPs passed the clumping plus thresholding and the PRS exhibited significant association with the AR in the validation cohort (HR = 1.54, 95% CI = (1.07, 2.22), p = 0.019). Further pathway analysis attributed the PRS genes into 13 categories, and the over-representation test identified 42 significant biological processes, the most significant of which is the cell morphogenesis (GO:0000902), with 4.08 fold of the percentage from homo species reference and FDR-adjusted P-value = 8.6×10-4.

CONCLUSIONS

Our results show the importance of donor-recipient mismatching in non-HLA regions. Additional work will be needed to understand the role of SNPs included in the PRS and to further improve donor-recipient genetic matching algorithms. Trial registry: Deterioration of Kidney Allograft Function Genomics (NCT00270712) and Genomics of Kidney Transplantation (NCT01714440) are registered on ClinicalTrials.gov.

Searching across-cohort relatives in 54,092 GWAS samples via encrypted genotype regression

Explicitly sharing individual level data in genomics studies has many merits comparing to sharing summary statistics, including more strict QCs, common statistical analyses, relative identification and improved statistical power in GWAS, but it is hampered by privacy or ethical constraints. In this study, we developed encG-reg, a regression approach that can detect relatives of various degrees based on encrypted genomic data, which is immune of ethical constraints. The encryption properties of encG-reg are based on the random matrix theory by masking the original genotypic matrix without sacrificing precision of individual-level genotype data. We established a connection between the dimension of a random matrix, which masked genotype matrices, and the required precision of a study for encrypted genotype data. encG-reg has false positive and false negative rates equivalent to sharing original individual level data, and is computationally efficient when searching relatives. We split the UK Biobank into their respective centers, and then encrypted the genotype data. We observed that the relatives estimated using encG-reg was equivalently accurate with the estimation by KING, which is a widely used software but requires original genotype data. In a more complex application, we launched a finely devised multi-center collaboration across 5 research institutes in China, covering 9 cohorts of 54,092 GWAS samples. encG-reg again identified true relatives existing across the cohorts with even different ethnic backgrounds and genotypic qualities. Our study clearly demonstrates that encrypted genomic data can be used for data sharing without loss of information or data sharing barrier.

Multi-ancestry and multi-trait genome-wide association meta-analyses inform clinical risk prediction for systemic lupus erythematosus

Systemic lupus erythematosus is a heritable autoimmune disease that predominantly affects young women. To improve our understanding of genetic etiology, we conduct multi-ancestry and multi-trait meta-analysis of genome-wide association studies, encompassing 12 systemic lupus erythematosus cohorts from 3 different ancestries and 10 genetically correlated autoimmune diseases, and identify 16 novel loci. We also perform transcriptome-wide association studies, computational drug repurposing analysis, and cell type enrichment analysis. We discover putative drug classes, including a histone deacetylase inhibitor that could be repurposed to treat lupus. We also identify multiple cell types enriched with putative target genes, such as non-classical monocytes and B cells, which may be targeted for future therapeutics. Using this newly assembled result, we further construct polygenic risk score models and demonstrate that integrating polygenic risk score with clinical lab biomarkers improves the diagnostic accuracy of systemic lupus erythematosus using the Vanderbilt BioVU and Michigan Genomics Initiative biobanks.

Principal Component Analyses (PCA)-based findings in population genetic studies are highly biased and must be reevaluated

Principal Component Analysis (PCA) is a multivariate analysis that reduces the complexity of datasets while preserving data covariance. The outcome can be visualized on colorful scatterplots, ideally with only a minimal loss of information. PCA applications, implemented in well-cited packages like EIGENSOFT and PLINK, are extensively used as the foremost analyses in population genetics and related fields (e.g., animal and plant or medical genetics). PCA outcomes are used to shape study design, identify, and characterize individuals and populations, and draw historical and ethnobiological conclusions on origins, evolution, dispersion, and relatedness. The replicability crisis in science has prompted us to evaluate whether PCA results are reliable, robust, and replicable. We analyzed twelve common test cases using an intuitive color-based model alongside human population data. We demonstrate that PCA results can be artifacts of the data and can be easily manipulated to generate desired outcomes. PCA adjustment also yielded unfavorable outcomes in association studies. PCA results may not be reliable, robust, or replicable as the field assumes. Our findings raise concerns about the validity of results reported in the population genetics literature and related fields that place a disproportionate reliance upon PCA outcomes and the insights derived from them. We conclude that PCA may have a biasing role in genetic investigations and that 32,000-216,000 genetic studies should be reevaluated. An alternative mixed-admixture population genetic model is discussed.

Shared Genetic Liability and Causal Associations Between Major Depressive Disorder and Cardiovascular Diseases

Major depressive disorder (MDD) is phenotypically associated with cardiovascular diseases (CVD). We aim to investigate mechanisms underlying relationships between MDD and CVD in the context of shared genetic variations. Polygenic overlap analysis was used to test genetic correlation and to analyze shared genetic variations between MDD and seven cardiovascular outcomes (coronary artery disease (CAD), heart failure, atrial fibrillation, stroke, systolic blood pressure, diastolic blood pressure, and pulse pressure measurement). Mendelian randomization analysis was used to uncover causal relationships between MDD and cardiovascular traits. By cross-trait meta-analysis, we identified a set of genomic loci shared between the traits of MDD and stroke. Putative causal genes for MDD and stroke were prioritized by fine-mapping of transcriptome-wide associations. Polygenic overlap analysis pointed toward substantial genetic variation overlap between MDD and CVD. Mendelian randomization analysis indicated that genetic liability to MDD has a causal effect on CAD and stroke. Comparison of genome-wide genes shared by MDD and CVD suggests 20q12 as a pleiotropic region conferring risk for both MDD and CVD. Cross-trait meta-analyses and fine-mapping of transcriptome-wide association signals identified novel risk genes for MDD and stroke, including RPL31P12, BORSC7, PNPT11, and PGF. Many genetic variations associated with MDD and CVD outcomes are shared, thus, pointing that genetic liability to MDD may also confer risk for stroke and CAD. Presented results shed light on mechanistic connections between MDD and CVD phenotypes.

Genetic evidence suggests posttraumatic stress disorder as a subtype of major depressive disorder

BACKGROUND

Major depressive disorder (MDD) and posttraumatic stress disorder (PTSD) are highly comorbid and exhibit strong correlations with one another. We aimed to investigate mechanisms of underlying relationships between PTSD and three kinds of depressive phenotypes, namely, MDD, depressed affect (DAF), and depression (DEP, including both MDD and the broad definition of depression).

METHODS

Genetic correlations between PTSD and the depressive phenotypes were tested using linkage disequilibrium score regression. Polygenic overlap analysis was used to estimate shared and trait-specific causal variants across a pair of traits. Causal relationships between PTSD and the depressive phenotypes were investigated using Mendelian randomization. Shared genomic loci between PTSD and MDD were identified using cross-trait meta-analysis.

RESULTS

Genetic correlations of PTSD with the depressive phenotypes were in the range of 0.71~0.80. The estimated numbers of causal variants were 14,565, 12,965, 10,565, and 4,986 for MDD, DEP, DAF, and PTSD, respectively. In each case, causal variants contributing to PTSD were completely or largely covered by causal variants defining each of the depressive phenotypes. Mendelian randomization analysis indicates that the genetically determined depressive phenotypes confer a causal effect on PTSD (b = 0.21~0.31). Notably, genetically determined PTSD confers a causal effect on DEP (b = 0.14) and DAF (b = 0.15), but not MDD. Cross-trait meta-analysis of MDD and PTSD identifies 47 genomic loci, including 29 loci shared between PTSD and MDD.

CONCLUSION

Evidence from shared genetics suggests that PTSD is a subtype of MDD. This study provides support to the efforts in reducing diagnostic heterogeneity in psychiatric nosology.

FUNDING

The National Key Research and Development Program of China (2018YFC1314300) and the National Natural Science Foundation of China (81471364 and 81971255).

A lead candidate functional single nucleotide polymorphism within the WARS2 gene associated with waist-hip-ratio does not alter RNA stability

We have prioritised a single nucleotide polymorphism (SNP) rs2645294 as one candidate functional SNP in the TBX15-WARS2 waist-hip-ratio locus using posterior probability analysis. This SNP is located in the 3' untranslated region of the WARS2 (tryptophanyl tRNA synthetase 2, mitochondrial) gene with which it has an expression quantitative trait in subcutaneous white adipose tissue. We show that transcripts of the WARS2 gene in a human white adipose cell line, heterozygous for the rs2645294 SNP, showed allelic imbalance. We tested whether the rs2645294 SNP altered WARS2 RNA stability using three different methods: actinomycin-D inhibition and RNA decay, mature and nascent RNA analysis and luciferase reporter assays. We found no evidence of a difference in RNA stability between the rs2645294 alleles indicating that the allelic expression imbalance was likely due to transcriptional regulation.

A genotype imputation method for de-identified haplotype reference information by using recurrent neural network

Genotype imputation estimates the genotypes of unobserved variants using the genotype data of other observed variants based on a collection of haplotypes for thousands of individuals, which is known as a haplotype reference panel. In general, more accurate imputation results were obtained using a larger size of haplotype reference panel. Most of the existing genotype imputation methods explicitly require the haplotype reference panel in precise form, but the accessibility of haplotype data is often limited, due to the requirement of agreements from the donors. Since de-identified information such as summary statistics or model parameters can be used publicly, imputation methods using de-identified haplotype reference information might be useful to enhance the quality of imputation results under the condition where the access of the haplotype data is limited. In this study, we proposed a novel imputation method that handles the reference panel as its model parameters by using bidirectional recurrent neural network (RNN). The model parameters are presented in the form of de-identified information from which the restoration of the genotype data at the individual-level is almost impossible. We demonstrated that the proposed method provides comparable imputation accuracy when compared with the existing imputation methods using haplotype datasets from the 1000 Genomes Project (1KGP) and the Haplotype Reference Consortium. We also considered a scenario where a subset of haplotypes is made available only in de-identified form for the haplotype reference panel. In the evaluation using the 1KGP dataset under the scenario, the imputation accuracy of the proposed method is much higher than that of the existing imputation methods. We therefore conclude that our RNN-based method is quite promising to further promote the data-sharing of sensitive genome data under the recent movement for the protection of individuals’ privacy.

Genotype imputation estimates the genotypes of unobserved variants using the genotype data of other observed variants based on a collection of genome data of a large number of individuals called a reference panel. In general, more accurate imputation results are obtained using a larger size of the reference panel. Although most of the existing imputation methods use the reference panel in an explicit form, the accessibility of genome data is often limited due to the requirement of agreements from the donors. We thus proposed a new imputation method that handles the reference panel as its model parameters by using bidirectional recurrent neural network. Since it is almost impossible to restore genome data at the individual-level from the model parameters, they can be shared publicly as the de-identified information even when the accessibility of the original reference panel is limited. We demonstrate that the proposed method provides comparable imputation accuracy with the existing methods. We also considered a scenario where a part of the genome data is made available only in de-identified form for the reference panel and have shown that the imputation accuracy of the proposed method is much higher than that of the existing methods under the scenario.

Two Synthetic 18-Way Outcrossed Populations of Diploid Budding Yeast with Utility for Complex Trait Dissection

Advanced-generation multiparent populations (MPPs) are a valuable tool for dissecting complex traits, having more power than genome-wide association studies to detect rare variants and higher resolution than F2 linkage mapping. To extend the advantages of MPPs in budding yeast, we describe the creation and characterization of two outbred MPPs derived from 18 genetically diverse founding strains. We carried out de novo assemblies of the genomes of the 18 founder strains, such that virtually all variation segregating between these strains is known, and represented those assemblies as Santa Cruz Genome Browser tracks. We discovered complex patterns of structural variation segregating among the founders, including a large deletion within the vacuolar ATPase VMA1, several different deletions within the osmosensor MSB2, a series of deletions and insertions at PRM7 and the adjacent BSC1, as well as copy number variation at the dehydrogenase ALD2 Resequenced haploid recombinant clones from the two MPPs have a median unrecombined block size of 66 kb, demonstrating that the population is highly recombined. We pool-sequenced the two MPPs to 3270× and 2226× coverage and demonstrated that we can accurately estimate local haplotype frequencies using pooled data. We further downsampled the pool-sequenced data to ∼20-40× and showed that local haplotype frequency estimates remained accurate, with median error rates 0.8 and 0.6% at 20× and 40×, respectively. Haplotypes frequencies are estimated much more accurately than SNP frequencies obtained directly from the same data. Deep sequencing of the two populations revealed that 10 or more founders are present at a detectable frequency for > 98% of the genome, validating the utility of this resource for the exploration of the role of standing variation in the architecture of complex traits.

Multi-tissue transcriptome analyses identify genetic mechanisms underlying neuropsychiatric traits

The genetic architecture of psychiatric disorders is characterized by a large number of small-effect variants1 located primarily in non-coding regions, suggesting that the underlying causal effects may influence disease risk by modulating gene expression2-4. We provide comprehensive analyses using transcriptome data from an unprecedented collection of tissues to gain pathophysiological insights into the role of the brain, neuroendocrine factors (adrenal gland) and gastrointestinal systems (colon) in psychiatric disorders. In each tissue, we perform PrediXcan analysis and identify trait-associated genes for schizophrenia (n associations = 499; n unique genes = 275), bipolar disorder (n associations = 17; n unique genes = 13), attention deficit hyperactivity disorder (n associations = 19; n unique genes = 12) and broad depression (n associations = 41; n unique genes = 31). Importantly, both PrediXcan and summary-data-based Mendelian randomization/heterogeneity in dependent instruments analyses suggest potentially causal genes in non-brain tissues, showing the utility of these tissues for mapping psychiatric disease genetic predisposition. Our analyses further highlight the importance of joint tissue approaches as 76% of the genes were detected only in difficult-to-acquire tissues.

Association score testing for rare variants and binary traits in family data with shared controls

Genome-wide association studies have been an important approach used to localize trait loci, with primary focus on common variants. The multiple rare variant-common disease hypothesis may explain the missing heritability remaining after accounting for identified common variants. Advances of sequencing technologies with their decreasing costs, coupled with methodological advances in the context of association studies in large samples, now make the study of rare variants at a genome-wide scale feasible. The resurgence of family-based association designs because of their advantage in studying rare variants has also stimulated more methods development, mainly based on linear mixed models (LMMs). Other tests such as score tests can have advantages over the LMMs, but to date have mainly been proposed for single-marker association tests. In this article, we extend several score tests (χcorrected2, WQLS, and SKAT) to the multiple variant association framework. We evaluate and compare their statistical performances relative with the LMM. Moreover, we show that three tests can be cast as the difference between marker allele frequencies (AFs) estimated in each of the group of affected and unaffected subjects. We show that these tests are flexible, as they can be based on related, unrelated or both related and unrelated subjects. They also make feasible an increasingly common design that only sequences a subset of affected subjects (related or unrelated) and uses for comparison publicly available AFs estimated in a group of healthy subjects. Finally, we show the great impact of linkage disequilibrium on the performance of all these tests.

Age at first birth in women is genetically associated with increased risk of schizophrenia

Previous studies have shown an increased risk for mental health problems in children born to both younger and older parents compared to children of average-aged parents. We previously used a novel design to reveal a latent mechanism of genetic association between schizophrenia and age at first birth in women (AFB). Here, we use independent data from the UK Biobank (N = 38,892) to replicate the finding of an association between predicted genetic risk of schizophrenia and AFB in women, and to estimate the genetic correlation between schizophrenia and AFB in women stratified into younger and older groups. We find evidence for an association between predicted genetic risk of schizophrenia and AFB in women (P-value = 1.12E-05), and we show genetic heterogeneity between younger and older AFB groups (P-value = 3.45E-03). The genetic correlation between schizophrenia and AFB in the younger AFB group is -0.16 (SE = 0.04) while that between schizophrenia and AFB in the older AFB group is 0.14 (SE = 0.08). Our results suggest that early, and perhaps also late, age at first birth in women is associated with increased genetic risk for schizophrenia in the UK Biobank sample. These findings contribute new insights into factors contributing to the complex bio-social risk architecture underpinning the association between parental age and offspring mental health.

A correction for sample overlap in genome-wide association studies in a polygenic pleiotropy-informed framework

BACKGROUND

There is considerable evidence that many complex traits have a partially shared genetic basis, termed pleiotropy. It is therefore useful to consider integrating genome-wide association study (GWAS) data across several traits, usually at the summary statistic level. A major practical challenge arises when these GWAS have overlapping subjects. This is particularly an issue when estimating pleiotropy using methods that condition the significance of one trait on the signficance of a second, such as the covariate-modulated false discovery rate (cmfdr).

RESULTS

We propose a method for correcting for sample overlap at the summary statistic level. We quantify the expected amount of spurious correlation between the summary statistics from two GWAS due to sample overlap, and use this estimated correlation in a simple linear correction that adjusts the joint distribution of test statistics from the two GWAS. The correction is appropriate for GWAS with case-control or quantitative outcomes. Our simulations and data example show that without correcting for sample overlap, the cmfdr is not properly controlled, leading to an excessive number of false discoveries and an excessive false discovery proportion. Our correction for sample overlap is effective in that it restores proper control of the false discovery rate, at very little loss in power.

CONCLUSIONS

With our proposed correction, it is possible to integrate GWAS summary statistics with overlapping samples in a statistical framework that is dependent on the joint distribution of the two GWAS.

Polygenic Scores for Major Depressive Disorder and Risk of Alcohol Dependence

Importance

Major depressive disorder (MDD) and alcohol dependence (AD) are heritable disorders with significant public health burdens, and they are frequently comorbid. Common genetic factors that influence the co-occurrence of MDD and AD have been sought in family, twin, and adoption studies, and results to date have been promising but inconclusive.

Objective

To examine whether AD and MDD overlap genetically, using a polygenic score approach.

Design, Settings, and Participants

Association analyses were conducted between MDD polygenic risk score (PRS) and AD case-control status in European ancestry samples from 4 independent genome-wide association study (GWAS) data sets: the Collaborative Study on the Genetics of Alcoholism (COGA); the Study of Addiction, Genetics, and Environment (SAGE); the Yale-Penn genetic study of substance dependence; and the National Health and Resilience in Veterans Study (NHRVS). Results from a meta-analysis of MDD (9240 patients with MDD and 9519 controls) from the Psychiatric Genomics Consortium were applied to calculate PRS at thresholds from P < .05 to P ≤ .99 in each AD GWAS data set.

Main Outcomes and Measures

Association between MDD PRS and AD.

Results

Participants analyzed included 788 cases (548 [69.5%] men; mean [SD] age, 38.2 [10.8] years) and 522 controls (151 [28.9.%] men; age [SD], 43.9 [11.6] years) from COGA; 631 cases (333 [52.8%] men; age [SD], 35.0 [7.7] years) and 756 controls (260 [34.4%] male; age [SD] 36.1 [7.7] years) from SAGE; 2135 cases (1375 [64.4%] men; age [SD], 39.4 [11.5] years) and 350 controls (126 [36.0%] men; age [SD], 43.5 [13.9] years) from Yale-Penn; and 317 cases (295 [93.1%] men; age [SD], 59.1 [13.1] years) and 1719 controls (1545 [89.9%] men; age [SD], 64.5 [13.3] years) from NHRVS. Higher MDD PRS was associated with a significantly increased risk of AD in all samples (COGA: best P = 1.7 × 10-6, R2 = 0.026; SAGE: best P = .001, R2 = 0.01; Yale-Penn: best P = .035, R2 = 0.0018; and NHRVS: best P = .004, R2 = 0.0074), with stronger evidence for association after meta-analysis of the 4 samples (best P = 3.3 × 10-9). In analyses adjusted for MDD status in 3 AD GWAS data sets, similar patterns of association were observed (COGA: best P = 7.6 × 10-6, R2 = 0.023; Yale-Penn: best P = .08, R2 = 0.0013; and NHRVS: best P = .006, R2 = 0.0072). After recalculating MDD PRS using MDD GWAS data sets without comorbid MDD-AD cases, significant evidence was observed for an association between the MDD PRS and AD in the meta-analysis of 3 GWAS AD samples without MDD cases (best P = .007).

Conclusions and Relevance

These results suggest that shared genetic susceptibility contributes modestly to MDD and AD comorbidity. Individuals with elevated polygenic risk for MDD may also be at risk for AD.

Reproduction and In-Depth Evaluation of Genome-Wide Association Studies and Genome-Wide Meta-analyses Using Summary Statistics

In line with open-source genetics, we report a novel linear regression technique for genome-wide association studies (GWAS), called Open GWAS algoriTHm (OATH). When individual-level data are not available, OATH can not only completely reproduce reported results from an experimental model, but also recover underreported results from other alternative models with a different combination of nuisance parameters using naïve summary statistics (NSS). OATH can also reliably evaluate all reported results in-depth (e.g., p-value variance analysis), as demonstrated for 42 Arabidopsis phenotypes under three magnesium (Mg) conditions. In addition, OATH can be used for consortium-driven genome-wide association meta-analyses (GWAMA), and can greatly improve the flexibility of GWAMA. A prototype of OATH is available in the Genetic Analysis Repository (https://github.com/gc5k/GEAR).