Embracing polygenicity: a review of methods and tools for psychiatric genetics research.

The availability of genome-wide genetic data on hundreds of thousands of people has led to an equally rapid growth in methodologies available to analyse these data. While the motivation for undertaking genome-wide association studies (GWAS) is identification of genetic markers associated with complex traits, once generated these data can be used for many other analyses. GWAS have demonstrated that complex traits exhibit a highly polygenic genetic architecture, often with shared genetic risk factors across traits. New methods to analyse data from GWAS are increasingly being used to address a diverse set of questions about the aetiology of complex traits and diseases, including psychiatric disorders. Here, we give an overview of some of these methods and present examples of how they have contributed to our understanding of psychiatric disorders. We consider: (i) estimation of the extent of genetic influence on traits, (ii) uncovering of shared genetic control between traits, (iii) predictions of genetic risk for individuals, (iv) uncovering of causal relationships between traits, (v) identifying causal single-nucleotide polymorphisms and genes or (vi) the detection of genetic heterogeneity. This classification helps organise the large number of recently developed methods, although some could be placed in more than one category. While some methods require GWAS data on individual people, others simply use GWAS summary statistics data, allowing novel well-powered analyses to be conducted at a low computational burden

Meta-analysis of genome-wide association studies of anxiety disorders.

Anxiety disorders (ADs), namely generalized AD, panic disorder and phobias, are common, etiologically complex conditions with a partially genetic basis. Despite differing on diagnostic definitions based on clinical presentation, ADs likely represent various expressions of an underlying common diathesis of abnormal regulation of basic threat-response systems. We conducted genome-wide association analyses in nine samples of European ancestry from seven large, independent studies. To identify genetic variants contributing to genetic susceptibility shared across interview-generated DSM-based ADs, we applied two phenotypic approaches: (1) comparisons between categorical AD cases and supernormal controls, and (2) quantitative phenotypic factor scores (FS) derived from a multivariate analysis combining information across the clinical phenotypes. We used logistic and linear regression, respectively, to analyze the association between these phenotypes and genome-wide single nucleotide polymorphisms. Meta-analysis for each phenotype combined results across the nine samples for over 18 000 unrelated individuals. Each meta-analysis identified a different genome-wide significant region, with the following markers showing the strongest association: for case-control contrasts, rs1709393 located in an uncharacterized non-coding RNA locus on chromosomal band 3q12.3 (P=1.65 × 10(-8)); for FS, rs1067327 within CAMKMT encoding the calmodulin-lysine N-methyltransferase on chromosomal band 2p21 (P=2.86 × 10(-9)). Independent replication and further exploration of these findings are needed to more fully understand the role of these variants in risk and expression of ADs.Molecular Psychiatry advance online publication, 12 January 2016; doi:10.1038/mp.2015.197

Genetic overlap between diagnostic subtypes of ischemic stroke

Background and Purpose: Despite moderate heritability, the phenotypic heterogeneity of ischemic stroke has hampered gene discovery, motivating analyses of diagnostic subtypes with reduced sample sizes. We assessed evidence for a shared genetic basis among the 3 major subtypes: large artery atherosclerosis (LAA), cardioembolism, and small vessel disease (SVD), to inform potential cross-subtype analyses. Methods: Analyses used genome-wide summary data for 12 389 ischemic stroke cases (including 2167 LAA, 2405 cardioembolism, and 1854 SVD) and 62 004 controls from the Metastroke consortium. For 4561 cases and 7094 controls, individual-level genotype data were also available. Genetic correlations between subtypes were estimated using linear mixed models and polygenic profile scores. Meta-analysis of a combined LAA-SVD phenotype (4021 cases and 51 976 controls) was performed to identify shared risk alleles. Results: High genetic correlation was identified between LAA and SVD using linear mixed models (rg=0.96, SE=0.47, P=9×10-4) and profile scores (rg=0.72; 95% confid

Dissection of genetic variation and evidence for pleiotropy in male pattern baldness.

Male pattern baldness (MPB) is a sex-limited, age-related, complex trait. We study MPB genetics in 205,327 European males from the UK Biobank. Here we show that MPB is strongly heritable and polygenic, with pedigree-heritability of 0.62 (SE = 0.03) estimated from close relatives, and SNP-heritability of 0.39 (SE = 0.01) from conventionally-unrelated males. We detect 624 near-independent genome-wide loci, contributing SNP-heritability of 0.25 (SE = 0.01), of which 26 X-chromosome loci explain 11.6%. Autosomal genetic variance is enriched for common variants and regions of lower linkage disequilibrium. We identify plausible genetic correlations between MPB and multiple sex-limited markers of earlier puberty, increased bone mineral density (r <sub>g</sub>  = 0.15) and pancreatic β-cell function (r <sub>g</sub>  = 0.12). Correlations with reproductive traits imply an effect on fitness, consistent with an estimated linear selection gradient of -0.018 per MPB standard deviation. Overall, we provide genetic insights into MPB: a phenotype of interest in its own right, with value as a model sex-limited, complex trait

Evaluating the Impact of Nonrandom Mating: Psychiatric Outcomes Among the Offspring of Pairs Diagnosed With Schizophrenia and Bipolar Disorder

Background: Nonrandom mating has been shown for psychiatric diagnoses, with hypothesized—but not quantified—implications for offspring liability. This national cohort study enumerated the incidence of major psychiatric disorders among the offspring of parent pairs affected with schizophrenia (SCZ) and/or bipolar disorder (BIP) (i.e., dual-affected pairs). Methods: Participants were all Swedish residents alive or born between 1968 and 2013 (n = 4,255,196 unique pairs and 8,343,951 offspring). Offspring with dual-affected, single-affected, and unaffected parents were followed (1973–2013) for incidence of broad psychiatric disorders. Primary outcomes included hazard ratio (HR) and cumulative incidence for SCZ and BIP in the offspring. Additional outcomes included any neuropsychiatric, anxiety, depressive, personality, or substance use disorders. Cumulative incidences of SCZ and BIP were used to inform heritability models for these disorders. Results: Hazards were highest within disorder (e.g., offspring of dual-SCZ pairs had sharply raised hazards for SCZ [HR = 55.3]); however, they were significantly raised for all diagnoses (HR range = 2.89–11.84). Incidences were significantly higher for the majority of outcomes, with 43.4% to 48.5% diagnosed with “any” disorder over follow-up. Risks were retained, with modest attenuations, for the offspring of heterotypic pairs. The estimated heritability of liability for SCZ (h2 = 0.62, 95% confidence interval = 0.55–0.70) and BIP (h2 = 0.52, 95% confidence interval = 0.46–0.58) did not differ significantly from estimates derived from single-affected parents. Conclusions: Risks for a broad spectrum of psychiatric diagnoses are significantly raised in the offspring of dual-affected parents, in line with expectations from a polygenic model of liability to disease risk. How these risks may contribute to population maintenance of these disorders is considered

mBAT-combo: A more powerful test to detect gene-trait associations from GWAS data

Gene-based association tests aggregate multiple SNP-trait associations into sets defined by gene boundaries and are widely used in post-GWAS analysis. A common approach for gene-based tests is to combine SNPs associations by computing the sum of χ2 statistics. However, this strategy ignores the directions of SNP effects, which could result in a loss of power for SNPs with masking effects, e.g., when the product of two SNP effects and the linkage disequilibrium (LD) correlation is negative. Here, we introduce “mBAT-combo,” a set-based test that is better powered than other methods to detect multi-SNP associations in the context of masking effects. We validate the method through simulations and applications to real data. We find that of 35 blood and urine biomarker traits in the UK Biobank, 34 traits show evidence for masking effects in a total of 4,273 gene-trait pairs, indicating that masking effects is common in complex traits. We further validate the improved power of our method in height, body mass index, and schizophrenia with different GWAS sample sizes and show that on average 95.7% of the genes detected only by mBAT-combo with smaller sample sizes can be identified by the single-SNP approach with a 1.7-fold increase in sample sizes. Eleven genes significant only in mBAT-combo for schizophrenia are confirmed by functionally informed fine-mapping or Mendelian randomization integrating gene expression data. The framework of mBAT-combo can be applied to any set of SNPs to refine trait-association signals hidden in genomic regions with complex LD structures

Polygenic risk prediction: why and when out-of-sample prediction R2 can exceed SNP-based heritability

In polygenic score (PGS) analysis, the coefficient of determination (R2) is a key statistic to evaluate efficacy. R2 is the proportion of phenotypic variance explained by the PGS, calculated in a cohort that is independent of the genome-wide association study (GWAS) that provided estimates of allelic effect sizes. The SNP-based heritability (hSNP2, the proportion of total phenotypic variances attributable to all common SNPs) is the theoretical upper limit of the out-of-sample prediction R2. However, in real data analyses R2 has been reported to exceed hSNP2, which occurs in parallel with the observation that hSNP2 estimates tend to decline as the number of cohorts being meta-analyzed increases. Here, we quantify why and when these observations are expected. Using theory and simulation, we show that if heterogeneities in cohort-specific hSNP2 exist, or if genetic correlations between cohorts are less than one, hSNP2 estimates can decrease as the number of cohorts being meta-analyzed increases. We derive conditions when the out-of-sample prediction R2 will be greater than hSNP2 and show the validity of our derivations with real data from a binary trait (major depression) and a continuous trait (educational attainment). Our research calls for a better approach to integrating information from multiple cohorts to address issues of between-cohort heterogeneity

Cohort profile: the Australian genetics of depression study

PURPOSE:Depression is the most common psychiatric disorder and the largest contributor to global disability. The Australian Genetics of Depression study was established to recruit a large cohort of individuals who have been diagnosed with depression at some point in their lifetime. The purpose of establishing this cohort is to investigate genetic and environmental risk factors for depression and response to commonly prescribed antidepressants. PARTICIPANTS:A total of 20 689 participants were recruited through the Australian Department of Human Services and a media campaign, 75% of whom were female. The average age of participants was 43 years±15 years. Participants completed an online questionnaire that consisted of a compulsory module that assessed self-reported psychiatric history, clinical depression using the Composite Interview Diagnostic Interview Short Form and experiences of using commonly prescribed antidepressants. Further voluntary modules assessed a wide range of traits of relevance to psychopathology. Participants who reported they were willing to provide a DNA sample (75%) were sent a saliva kit in the mail. FINDINGS TO DATE:95% of participants reported being given a diagnosis of depression by a medical practitioner and 88% met the criteria for a lifetime depressive episode. 68% of the sample report having been diagnosed with another psychiatric disorder in addition to depression. In line with findings from clinical trials, only 33% of the sample report responding well to the first antidepressant they were prescribed. FUTURE PLANS:A number of analyses to investigate the genetic architecture of depression and common comorbidities will be conducted. The cohort will contribute to the global effort to identify genetic variants that increase risk to depression. Furthermore, a thorough investigation of genetic and psychosocial predictors of antidepressant response and side effects is planned.Enda M Byrne, Katherine M Kirk, Sarah E Medland, John J McGrath, Lucia Colodro-Conde, Richard Parker ... et al

Quantifying Between-Cohort and Between-Sex Genetic Heterogeneity in Major Depressive Disorder

Major depressive disorder (MDD) is clinically heterogeneous with prevalence rates twice as high in women as in men. There are many possible sources of heterogeneity in MDD most of which are not measured in a sufficiently comparable way across study samples. Here, we assess genetic heterogeneity based on two fundamental measures, between-cohort and between-sex heterogeneity. First, we used genome-wide association study (GWAS) summary statistics to investigate between-cohort genetic heterogeneity using the 29 research cohorts of the Psychiatric Genomics Consortium (PGC; N cases = 16,823, N controls = 25,632) and found that some of the cohort heterogeneity can be attributed to ascertainment differences (such as recruitment of cases from hospital vs community sources). Second, we evaluated between-sex genetic heterogeneity using GWAS summary statistics from the PGC, Kaiser Permanente GERA, UK Biobank and the Danish iPSYCH studies but did not find convincing evidence for genetic differences between the sexes. We conclude that there is no evidence that the heterogeneity between MDD data sets and between sexes reflects genetic heterogeneity. Larger sample sizes with detailed phenotypic records and genomic data remain the key to overcome heterogeneity inherent in assessment of MDD