Disease and Polygenic Architecture: Avoid Trio Design and Appropriately Account for Unscreened Control Subjects for Common Disease

Pourcain B, Robinson EB, Bourgeron T, Baron-Cohen S, Børglum AD, Martin HC, Warrier V. An axis of genetic heterogeneity in autism is indexed by age at diagnosis and is associated with varying developmental and mental health profiles

There is growing recognition that earliest signs of autism need not clearly manifest in the first three years of life. To what extent is this variation in developmental trajectories associated with age at autism diagnosis? Does the genetic profile of autism vary with age at autism diagnosis? Using longitudinal data from four birth cohorts, we demonstrate that two different trajectories of socio-emotional behaviours are associated with age at diagnosis. We further demonstrate that the age at autism diagnosis is partly heritable (h2 SNP = 0.12, s.e.m = 0.01), and is associated with two moderately correlated (rg = 0.38, s.e.m = 0.07) autism polygenic factors. One of these factors is associated with earlier diagnosis of autism, lower social and communication abilities in early childhood. The second factor is associated with later autism diagnosis, increased socio-emotional difficulties in adolescence, and has moderate to high positive genetic correlations with Attention-Deficit/Hyperactivity Disorder, mental health conditions, and trauma. Overall, our research identifies an axis of heterogeneity in autism, indexed by age at diagnosis, which partly explains heterogeneity in autism and the profiles of co-occurring neurodevelopmental and mental health profiles. Our findings have important implications for how we conceptualise autism and provide one model to explain some of the diversity within autism.

The genetic architecture of youth anxiety: a study protocol

BACKGROUND

Anxiety disorders are the most common psychiatric problems among Canadian youth and typically have an onset in childhood or adolescence. They are characterized by high rates of relapse and chronicity, often resulting in substantial impairment across the lifespan. Genetic factors play an important role in the vulnerability toward anxiety disorders. However, genetic contribution to anxiety in youth is not well understood and can change across developmental stages. Large-scale genetic studies of youth are needed with detailed assessments of symptoms of anxiety disorders and their major comorbidities to inform early intervention or preventative strategies and suggest novel targets for therapeutics and personalization of care.

METHODS

The Genetic Architecture of Youth Anxiety (GAYA) study is a Pan-Canadian effort of clinical and genetic experts with specific recruitment sites in Calgary, Halifax, Hamilton, Toronto, and Vancouver. Youth aged 10-19 (n = 13,000) will be recruited from both clinical and community settings and will provide saliva samples, complete online questionnaires on demographics, symptoms of mental health concerns, and behavioural inhibition, and complete neurocognitive tasks. A subset of youth will be offered access to a self-managed Internet-based cognitive behavioral therapy resource. Analyses will focus on the identification of novel genetic risk loci for anxiety disorders in youth and assess how much of the genetic risk for anxiety disorders is unique or shared across the life span.

DISCUSSION

Results will substantially inform early intervention or preventative strategies and suggest novel targets for therapeutics and personalization of care. Given that the GAYA study will be the biggest genomic study of anxiety disorders in youth in Canada, this project will further foster collaborations nationally and across the world.

Distinguishing different psychiatric disorders using DDx-PRS

Despite great progress on methods for case-control polygenic prediction (e.g. schizophrenia vs. control), there remains an unmet need for a method that genetically distinguishes clinically related disorders (e.g. schizophrenia (SCZ) vs. bipolar disorder (BIP) vs. depression (MDD) vs. control); such a method could have important clinical value, especially at disorder onset when differential diagnosis can be challenging. Here, we introduce a method, Differential Diagnosis-Polygenic Risk Score (DDx-PRS), that jointly estimates posterior probabilities of each possible diagnostic category (e.g. SCZ=50%, BIP=25%, MDD=15%, control=10%) by modeling variance/covariance structure across disorders, leveraging case-control polygenic risk scores (PRS) for each disorder (computed using existing methods) and prior clinical probabilities for each diagnostic category. DDx-PRS uses only summary-level training data and does not use tuning data, facilitating implementation in clinical settings. In simulations, DDx-PRS was well-calibrated (whereas a simpler approach that analyzes each disorder marginally was poorly calibrated), and effective in distinguishing each diagnostic category vs. the rest. We then applied DDx-PRS to Psychiatric Genomics Consortium SCZ/BIP/MDD/control data, including summary-level training data from 3 case-control GWAS ( N =41,917-173,140 cases; total N =1,048,683) and held-out test data from different cohorts with equal numbers of each diagnostic category (total N =11,460). DDx-PRS was well-calibrated and well-powered relative to these training sample sizes, attaining AUCs of 0.66 for SCZ vs. rest, 0.64 for BIP vs. rest, 0.59 for MDD vs. rest, and 0.68 for control vs. rest. DDx-PRS produced comparable results to methods that leverage tuning data, confirming that DDx-PRS is an effective method. True diagnosis probabilities in top deciles of predicted diagnosis probabilities were considerably larger than prior baseline probabilities, particularly in projections to larger training sample sizes, implying considerable potential for clinical utility under certain circumstances. In conclusion, DDx-PRS is an effective method for distinguishing clinically related disorders.

Estimating Disorder Probability Based on Polygenic Prediction Using the BPC Approach

Polygenic Scores (PGSs) summarize an individual's genetic propensity for a given trait in a single value, based on SNP effect sizes derived from Genome-Wide Association Study (GWAS) results. Methods have been developed that apply Bayesian approaches to improve the prediction accuracy of PGSs through optimization of estimated effect sizes. While these methods are generally well-calibrated for continuous traits (implying the predicted values are on average equal to the true trait values), they are not well-calibrated for binary disorder traits in ascertained samples. This is a problem because well-calibrated PGSs are needed to reliably compute the absolute disorder probability for an individual to facilitate future clinical implementation. Here we introduce the Bayesian polygenic score Probability Conversion (BPC) approach, which computes an individual's predicted disorder probability using GWAS summary statistics, an existing Bayesian PGS method (e.g. PRScs, SBayesR), the individual's genotype data, and a prior disorder probability. The BPC approach transforms the PGS to its underlying liability scale, computes the variances of the PGS in cases and controls, and applies Bayes' Theorem to compute the absolute disorder probability; it is practical in its application as it does not require a tuning dataset with both genotype and phenotype data. We applied the BPC approach to extensive simulated data and empirical data of nine disorders. The BPC approach yielded well-calibrated results that were consistently better than the results of another recently published approach.

Assortative mating and parental genetic relatedness contribute to the pathogenicity of variably expressive variants

We examined more than 97,000 families from four neurodevelopmental disease cohorts and the UK Biobank to identify phenotypic and genetic patterns in parents contributing to neurodevelopmental disease risk in children. We identified within- and cross-disorder correlations between six phenotypes in parents and children, such as obsessive-compulsive disorder (R = 0.32-0.38, p < 10-126). We also found that measures of sub-clinical autism features in parents are associated with several autism severity measures in children, including biparental mean Social Responsiveness Scale scores and proband Repetitive Behaviors Scale scores (regression coefficient = 0.14, p = 3.38 × 10-4). We further describe patterns of phenotypic similarity between spouses, where spouses show correlations for six neurological and psychiatric phenotypes, including a within-disorder correlation for depression (R = 0.24-0.68, p < 0.001) and a cross-disorder correlation between anxiety and bipolar disorder (R = 0.09-0.22, p < 10-92). Using a simulated population, we also found that assortative mating can lead to increases in disease liability over generations and the appearance of "genetic anticipation" in families carrying rare variants. We identified several families in a neurodevelopmental disease cohort where the proband inherited multiple rare variants in disease-associated genes from each of their affected parents. We further identified parental relatedness as a risk factor for neurodevelopmental disorders through its inverse relationship with variant pathogenicity and propose that parental relatedness modulates disease risk by increasing genome-wide homozygosity in children (R = 0.05-0.26, p < 0.05). Our results highlight the utility of assessing parent phenotypes and genotypes toward predicting features in children who carry rare variably expressive variants and implicate assortative mating as a risk factor for increased disease severity in these families.

Investigating Direct and Indirect Genetic Effects in Attention-Deficit/Hyperactivity Disorder Using Parent-Offspring Trios

BACKGROUND

Attention-deficit/hyperactivity disorder (ADHD) is highly heritable, but little is known about the relative effects of transmitted (i.e., direct) and nontransmitted (i.e., indirect) common variant risks. Using parent-offspring trios, we tested whether polygenic liability for neurodevelopmental and psychiatric disorders and lower cognitive ability is overtransmitted to ADHD probands. We also tested for indirect or genetic nurture effects by examining whether nontransmitted ADHD polygenic liability is elevated. Finally, we examined whether complete trios are representative of the clinical ADHD population.

METHODS

Polygenic risk scores (PRSs) for ADHD, anxiety, autism, bipolar disorder, depression, obsessive-compulsive disorder, schizophrenia, Tourette syndrome, and cognitive ability were calculated in UK control subjects (n = 5081), UK probands with ADHD (n = 857), their biological parents (n = 328 trios), and also a replication sample of 844 ADHD trios.

RESULTS

ADHD PRSs were overtransmitted and cognitive ability and obsessive-compulsive disorder PRSs were undertransmitted. These results were independently replicated. Overtransmission of polygenic liability was not observed for other disorders. Nontransmitted alleles were not enriched for ADHD liability compared with control subjects. Probands from incomplete trios had more hyperactive-impulsive and conduct disorder symptoms, lower IQ, and lower socioeconomic status than complete trios. PRS did not vary by trio status.

CONCLUSIONS

The results support direct transmission of polygenic liability for ADHD and cognitive ability from parents to offspring, but not for other neurodevelopmental/psychiatric disorders. They also suggest that nontransmitted neurodevelopmental/psychiatric parental alleles do not contribute indirectly to ADHD via genetic nurture. Furthermore, ascertainment of complete ADHD trios may be nonrandom, in terms of demographic and clinical factors.

Pervasive Downward Bias in Estimates of Liability-Scale Heritability in Genome-wide Association Study Meta-analysis: A Simple Solution

BACKGROUND

Single nucleotide polymorphism-based heritability is a fundamental quantity in the genetic analysis of complex traits. For case-control phenotypes, for which the continuous distribution of risk in the population is unobserved, observed-scale heritability estimates must be transformed to the more interpretable liability scale. This article describes how the field standard approach incorrectly performs the liability correction in that it does not appropriately account for variation in the proportion of cases across the cohorts comprising the meta-analysis. We propose a simple solution that incorporates cohort-specific ascertainment using the summation of effective sample sizes across cohorts. This solution is applied at the stage of single nucleotide polymorphism-based heritability estimation and does not require generating updated meta-analytic genome-wide association study summary statistics.

METHODS

We began by performing a series of simulations to examine the ability of the standard approach and our proposed approach to recapture liability-scale heritability in the population. We went on to examine the differences in estimates obtained from these 2 approaches for real data for 12 major case-control genome-wide association studies of psychiatric and neurologic traits.

RESULTS

We found that the field standard approach for performing the liability conversion can downwardly bias estimates by as much as approximately 50% in simulation and approximately 30% in real data.

CONCLUSIONS

Prior estimates of liability-scale heritability for genome-wide association study meta-analysis may be drastically underestimated. To this end, we strongly recommend using our proposed approach of using the sum of effective sample sizes across contributing cohorts to obtain unbiased estimates.

Genetic correlates of phenotypic heterogeneity in autism

The substantial phenotypic heterogeneity in autism limits our understanding of its genetic etiology. To address this gap, here we investigated genetic differences between autistic individuals (nmax = 12,893) based on core and associated features of autism, co-occurring developmental disabilities and sex. We conducted a comprehensive factor analysis of core autism features in autistic individuals and identified six factors. Common genetic variants were associated with the core factors, but de novo variants were not. We found that higher autism polygenic scores (PGS) were associated with lower likelihood of co-occurring developmental disabilities in autistic individuals. Furthermore, in autistic individuals without co-occurring intellectual disability (ID), autism PGS are overinherited by autistic females compared to males. Finally, we observed higher SNP heritability for autistic males and for autistic individuals without ID. Deeper phenotypic characterization will be critical in determining how the complex underlying genetics shape cognition, behavior and co-occurring conditions in autism.

Polygenic risk scores of endo-phenotypes identify the effect of genetic background in congenital heart disease

Congenital heart disease (CHD) is a rare structural defect that occurs in ∼1% of live births. Studies on CHD genetic architecture have identified pathogenic single-gene mutations in less than 30% of cases. Single-gene mutations often show incomplete penetrance and variable expressivity. Therefore, we hypothesize that genetic background may play a role in modulating disease expression. Polygenic risk scores (PRSs) aggregate effects of common genetic variants to investigate whether, cumulatively, these variants are associated with disease penetrance or severity. However, the major limitations in this field have been in generating sufficient sample sizes for these studies. Here we used CHD-phenotype matched genome-wide association study (GWAS) summary statistics from the UK Biobank (UKBB) as our base study and whole-genome sequencing data from the CHD cohort (n1 = 711 trios, n2 = 362 European trios) of the Gabriella Miller Kids First dataset as our target study to develop PRSs for CHD. PRSs estimated using a GWAS for heart valve problems and heart murmur explain 2.5% of the variance in case-control status of CHD (all SNVs, p = 7.90 × 10-3; fetal cardiac SNVs, p = 8.00 × 10-3) and 1.8% of the variance in severity of CHD (fetal cardiac SNVs, p = 6.20 × 10-3; all SNVs, p = 0.015). These results show that common variants captured in CHD phenotype-matched GWASs have a modest but significant contribution to phenotypic expression of CHD. Further exploration of the cumulative effect of common variants is necessary for understanding the complex genetic etiology of CHD and other rare diseases.

Integrated analysis of direct and proxy genome wide association studies highlights polygenicity of Alzheimer's disease outside of the APOE region

Recent meta-analyses combining direct genome-wide association studies (GWAS) with those of family history (GWAX) have indicated very low SNP heritability of Alzheimer's disease (AD). These low estimates may call into question the prospects of continued progress in genetic discovery for AD within the spectrum of common variants. We highlight dramatic downward biases in previous methods, and we validate a novel method for the estimation of SNP heritability via integration of GWAS and GWAX summary data. We apply our method to investigate the genetic architecture of AD using GWAX from UK Biobank and direct case-control GWAS from the International Genomics of Alzheimer's Project (IGAP). We estimate the liability scale common variant SNP heritability of Clinical AD outside of APOE region at ~7-11%, and we project the corresponding estimate for AD pathology to be up to approximately 23%. We estimate that nearly 90% of common variant SNP heritability of Clinical AD exists outside the APOE region. Rare variants not tagged in standard GWAS may account for additional variance. Our results indicate that, while GWAX for AD in UK Biobank may result in greater attenuation of genetic effects beyond that conventionally assumed, it does not introduce appreciable contamination of signal by genetically distinct traits relative to direct case-control GWAS in IGAP. Genetic risk for AD represents a strong effect of APOE superimposed upon a highly polygenic background.

How genetic analysis may contribute to the understanding of avoidant/restrictive food intake disorder (ARFID)

Avoidant/restrictive food intake disorder (ARFID) was introduced in the fifth edition of the Diagnostic and Statistical Manual of Mental Disorders (DSM-5). Unlike anorexia nervosa, ARFID is characterised by avoidant or restricted food intake that is not driven by weight or body shape-related concerns. As with other eating disorders, it is expected that ARFID will have a significant genetic risk component; however, sufficiently large-scale genetic investigations are yet to be performed in this group of patients. This narrative review considers the current literature on the diagnosis, presentation, and course of ARFID, including evidence for different presentations, and identifies fundamental questions about how ARFID might fit into the fluid landscape of other eating and mental disorders. In the absence of large ARFID GWAS, we consider genetic research on related conditions to point to possible features or mechanisms relevant to future ARFID investigations, and discuss the theoretical and clinical implications an ARFID GWAS. An argument for a collaborative approach to recruit ARFID participants for genome-wide association study is presented, as understanding the underlying genomic architecture of ARFID will be a key step in clarifying the biological mechanisms involved, and the development of interventions and treatments for this serious, and often debilitating disorder.

Exploring the genetic heterogeneity in major depression across diagnostic criteria

Major depressive disorder (MDD) is defined differently across genetic research studies and this may be a key source of heterogeneity. While previous literature highlights differences between minimal and strict phenotypes, the components contributing to this heterogeneity have not been identified. Using the cardinal symptoms (depressed mood/anhedonia) as a baseline, we build MDD phenotypes using five components-(1) five or more symptoms, (2) episode duration, (3) functional impairment, (4) episode persistence, and (5) episode recurrence-to determine the contributors to such heterogeneity. Thirty-two depression phenotypes which systematically incorporate different combinations of MDD components were created using the mental health questionnaire data within the UK Biobank. SNP-based heritabilities and genetic correlations with three previously defined major depression phenotypes were calculated (Psychiatric Genomics Consortium (PGC) defined depression, 23andMe self-reported depression and broad depression) and differences between estimates analysed. All phenotypes were heritable (h2SNP range: 0.102-0.162) and showed substantial genetic correlations with other major depression phenotypes (Rg range: 0.651-0.895 (PGC); 0.652-0.837 (23andMe); 0.699-0.900 (broad depression)). The strongest effect on SNP-based heritability was from the requirement for five or more symptoms (1.4% average increase) and for a long episode duration (2.7% average decrease). No significant differences were noted between genetic correlations. While there is some variation, the two cardinal symptoms largely reflect the genetic aetiology of phenotypes incorporating more MDD components. These components may index severity, however, their impact on heterogeneity in genetic results is likely to be limited.

Exploring Drugs and Vaccines Associated with Altered Risks and Severity of COVID-19: A UK Biobank Cohort Study of All ATC Level-4 Drug Categories Reveals Repositioning Opportunities

Effective therapies for COVID-19 are still lacking, and drug repositioning is a promising approach to address this problem. Here, we adopted a medical informatics approach to repositioning. We leveraged a large prospective cohort, the UK-Biobank (UKBB, N ~ 397,000), and studied associations of prior use of all level-4 ATC drug categories (N = 819, including vaccines) with COVID-19 diagnosis and severity. Effects of drugs on the risk of infection, disease severity, and mortality were investigated separately. Logistic regression was conducted, controlling for main confounders. We observed strong and highly consistent protective associations with statins. Many top-listed protective drugs were also cardiovascular medications, such as angiotensin-converting enzyme inhibitors (ACEI), angiotensin receptor blockers (ARB), calcium channel blocker (CCB), and beta-blockers. Some other drugs showing protective associations included biguanides (metformin), estrogens, thyroid hormones, proton pump inhibitors, and testosterone-5-alpha reductase inhibitors, among others. We also observed protective associations by influenza, pneumococcal, and several other vaccines. Subgroup and interaction analyses were also conducted, which revealed differences in protective effects in various subgroups. For example, protective effects of flu/pneumococcal vaccines were weaker in obese individuals, while protection by statins was stronger in cardiovascular patients. To conclude, our analysis revealed many drug repositioning candidates, for example several cardiovascular medications. Further studies are required for validation.

Analysis of genetic differences between psychiatric disorders: exploring pathways and cell types/tissues involved and ability to differentiate the disorders by polygenic scores

Although displaying genetic correlations, psychiatric disorders are clinically defined as categorical entities as they each have distinguishing clinical features and may involve different treatments. Identifying differential genetic variations between these disorders may reveal how the disorders differ biologically and help to guide more personalized treatment. Here we presented a statistical framework and comprehensive analysis to identify genetic markers differentially associated with various psychiatric disorders/traits based on GWAS summary statistics, covering 18 psychiatric traits/disorders and 26 comparisons. We also conducted comprehensive analysis to unravel the genes, pathways and SNP functional categories involved, and the cell types and tissues implicated. We also assessed how well one could distinguish between psychiatric disorders by polygenic risk scores (PRS). SNP-based heritabilities (h2snp) were significantly larger than zero for most comparisons. Based on current GWAS data, PRS have mostly modest power to distinguish between psychiatric disorders. For example, we estimated that AUC for distinguishing schizophrenia from major depressive disorder (MDD), bipolar disorder (BPD) from MDD and schizophrenia from BPD were 0.694, 0.602 and 0.618, respectively, while the maximum AUC (based on h2snp) were 0.763, 0.749 and 0.726, respectively. We also uncovered differences in each pair of studied traits in terms of their differences in genetic correlation with comorbid traits. For example, clinically defined MDD appeared to more strongly genetically correlated with other psychiatric disorders and heart disease, when compared to non-clinically defined depression in UK Biobank. Our findings highlight genetic differences between psychiatric disorders and the mechanisms involved. PRS may help differential diagnosis of selected psychiatric disorders in the future with larger GWAS samples.

Risk in Relatives, Heritability, SNP-Based Heritability, and Genetic Correlations in Psychiatric Disorders: A Review

The genetic contribution to psychiatric disorders is observed through the increased rates of disorders in the relatives of those diagnosed with disorders. These increased rates are observed to be nonspecific; for example, children of those with schizophrenia have increased rates of schizophrenia but also a broad range of other psychiatric diagnoses. While many factors contribute to risk, epidemiological evidence suggests that the genetic contribution carries the highest risk burden. The patterns of inheritance are consistent with a polygenic architecture of many contributing risk loci. The genetic studies of the past decade have provided empirical evidence identifying thousands of DNA variants associated with psychiatric disorders. Here, we describe how these latest results are consistent with observations from epidemiology. We provide an R tool (CHARRGe) to calculate genetic parameters from epidemiological parameters and vice versa. We discuss how the single nucleotide polymorphism-based estimates of heritability and genetic correlation relate to those estimated from family records.

Common genetic risk variants identified in the SPARK cohort support DDHD2 as a candidate risk gene for autism

Autism spectrum disorder (ASD) is a highly heritable neurodevelopmental disorder. Large genetically informative cohorts of individuals with ASD have led to the identification of a limited number of common genome-wide significant (GWS) risk loci to date. However, many more common genetic variants are expected to contribute to ASD risk given the high heritability. Here, we performed a genome-wide association study (GWAS) on 6222 case-pseudocontrol pairs from the Simons Foundation Powering Autism Research for Knowledge (SPARK) dataset to identify additional common genetic risk factors and molecular mechanisms underlying risk for ASD. We identified one novel GWS locus from the SPARK GWAS and four significant loci, including an additional novel locus from meta-analysis with a previous GWAS. We replicated the previous observation of significant enrichment of ASD heritability within regulatory regions of the developing cortex, indicating that disruption of gene regulation during neurodevelopment is critical for ASD risk. We further employed a massively parallel reporter assay (MPRA) and identified a putative causal variant at the novel locus from SPARK GWAS with strong impacts on gene regulation (rs7001340). Expression quantitative trait loci data demonstrated an association between the risk allele and decreased expression of DDHD2 (DDHD domain containing 2) in both adult and prenatal brains. In conclusion, by integrating genetic association data with multi-omic gene regulatory annotations and experimental validation, we fine-mapped a causal risk variant and demonstrated that DDHD2 is a novel gene associated with ASD risk.

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.

Large-scale GWAS reveals insights into the genetic architecture of same-sex sexual behavior

Twin and family studies have shown that same-sex sexual behavior is partly genetically influenced, but previous searches for specific genes involved have been underpowered. We performed a genome-wide association study (GWAS) on 477,522 individuals, revealing five loci significantly associated with same-sex sexual behavior. In aggregate, all tested genetic variants accounted for 8 to 25% of variation in same-sex sexual behavior, only partially overlapped between males and females, and do not allow meaningful prediction of an individual's sexual behavior. Comparing these GWAS results with those for the proportion of same-sex to total number of sexual partners among nonheterosexuals suggests that there is no single continuum from opposite-sex to same-sex sexual behavior. Overall, our findings provide insights into the genetics underlying same-sex sexual behavior and underscore the complexity of sexuality.

Psychiatric Polygenic Risk Scores as Predictor for Attention Deficit/Hyperactivity Disorder and Autism Spectrum Disorder in a Clinical Child and Adolescent Sample

Neurodevelopmental disorders such as attention deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) are highly heritable and influenced by many single nucleotide polymorphisms (SNPs). SNPs can be used to calculate individual polygenic risk scores (PRS) for a disorder. We aim to explore the association between the PRS for ADHD, ASD and for Schizophrenia (SCZ), and ADHD and ASD diagnoses in a clinical child and adolescent population. Based on the most recent genome wide association studies of ADHD, ASD and SCZ, PRS of each disorder were calculated for individuals of a clinical child and adolescent target sample (N = 688) and for adult controls (N = 943). We tested with logistic regression analyses for an association with (1) a single diagnosis of ADHD (N = 280), (2) a single diagnosis of ASD (N = 295), and (3) combining the two diagnoses, thus subjects with either ASD, ADHD or both (N = 688). Our results showed a significant association of the ADHD PRS with ADHD status (OR 1.6, P = 1.39 × 10⁻⁰⁷) and with the combined ADHD/ASD status (OR 1.36, P = 1.211 × 10⁻⁰⁵), but not with ASD status (OR 1.14, P = 1). No associations for the ASD and SCZ PRS were observed. In sum, the PRS of ADHD is significantly associated with the combined ADHD/ASD status. Yet, this association is primarily driven by ADHD status, suggesting disorder specific genetic effects of the ADHD PRS.

PolyGEE: a generalized estimating equation approach to the efficient and robust estimation of polygenic effects in large-scale association studies

To quantify polygenic effects, i.e. undetected genetic effects, in large-scale association studies, we propose a generalized estimating equation (GEE) based estimation framework. We develop a marginal model for single-variant association test statistics of complex diseases that generalizes existing approaches such as LD Score regression and that is applicable to population-based designs, to family-based designs or to arbitrary combinations of both. We extend the standard GEE approach so that the parameters of the proposed marginal model can be estimated based on working-correlation/linkage-disequilibrium (LD) matrices from external reference panels. Our method achieves substantial efficiency gains over standard approaches, while it is robust against misspecification of the LD structure, i.e. the LD structure of the reference panel can differ substantially from the true LD structure in the study population. In simulation studies and in applications to population-based and family-based studies, we illustrate the features of the proposed GEE framework. Our results suggest that our approach can be up to 100% more efficient than existing methodology.

Genetic influences on eight psychiatric disorders based on family data of 4 408 646 full and half-siblings, and genetic data of 333 748 cases and controls

BACKGROUND

Most studies underline the contribution of heritable factors for psychiatric disorders. However, heritability estimates depend on the population under study, diagnostic instruments, and study designs that each has its inherent assumptions, strengths, and biases. We aim to test the homogeneity in heritability estimates between two powerful, and state of the art study designs for eight psychiatric disorders.

METHODS

We assessed heritability based on data of Swedish siblings (N = 4 408 646 full and maternal half-siblings), and based on summary data of eight samples with measured genotypes (N = 125 533 cases and 208 215 controls). All data were based on standard diagnostic criteria. Eight psychiatric disorders were studied: (1) alcohol dependence (AD), (2) anorexia nervosa, (3) attention deficit/hyperactivity disorder (ADHD), (4) autism spectrum disorder, (5) bipolar disorder, (6) major depressive disorder, (7) obsessive-compulsive disorder (OCD), and (8) schizophrenia.

RESULTS

Heritability estimates from sibling data varied from 0.30 for Major Depression to 0.80 for ADHD. The estimates based on the measured genotypes were lower, ranging from 0.10 for AD to 0.28 for OCD, but were significant, and correlated positively (0.19) with national sibling-based estimates. When removing OCD from the data the correlation increased to 0.50.

CONCLUSIONS

Given the unique character of each study design, the convergent findings for these eight psychiatric conditions suggest that heritability estimates are robust across different methods. The findings also highlight large differences in genetic and environmental influences between psychiatric disorders, providing future directions for etiological psychiatric research.

A simple approximation to the bias of gene-environment interactions in case-control studies with silent disease

One of the most important research areas in case-control Genome-Wide Association Studies is to determine how the effect of a genotype varies across the environment or to measure the gene-environment interaction (G × E). We consider the scenario when some of the "healthy" controls actually have the disease and when the frequency of these latent cases varies by the environmental variable of interest. In this scenario, performing logistic regression with the clinically diagnosed disease status as an outcome variable and will result in biased estimates of G × E interaction. Here, we derive a general theoretical approximation to the bias in the estimates of the G × E interaction and show, through extensive simulation, that this approximation is accurate in finite samples. Moreover, we apply this approximation to evaluate the bias in the effect estimates of the genetic variants related to mitochondrial proteins a large-scale prostate cancer study.

Does Childhood Trauma Moderate Polygenic Risk for Depression? A Meta-analysis of 5765 Subjects From the Psychiatric Genomics Consortium

BACKGROUND

The heterogeneity of genetic effects on major depressive disorder (MDD) may be partly attributable to moderation of genetic effects by environment, such as exposure to childhood trauma (CT). Indeed, previous findings in two independent cohorts showed evidence for interaction between polygenic risk scores (PRSs) and CT, albeit in opposing directions. This study aims to meta-analyze MDD-PRS × CT interaction results across these two and other cohorts, while applying more accurate PRSs based on a larger discovery sample.

METHODS

Data were combined from 3024 MDD cases and 2741 control subjects from nine cohorts contributing to the MDD Working Group of the Psychiatric Genomics Consortium. MDD-PRS were based on a discovery sample of ∼110,000 independent individuals. CT was assessed as exposure to sexual or physical abuse during childhood. In a subset of 1957 cases and 2002 control subjects, a more detailed five-domain measure additionally included emotional abuse, physical neglect, and emotional neglect.

RESULTS

MDD was associated with the MDD-PRS (odds ratio [OR] = 1.24, p = 3.6 × 10-5, R2 = 1.18%) and with CT (OR = 2.63, p = 3.5 × 10-18 and OR = 2.62, p = 1.4 ×10-5 for the two- and five-domain measures, respectively). No interaction was found between MDD-PRS and the two-domain and five-domain CT measure (OR = 1.00, p = .89 and OR = 1.05, p = .66).

CONCLUSIONS

No meta-analytic evidence for interaction between MDD-PRS and CT was found. This suggests that the previously reported interaction effects, although both statistically significant, can best be interpreted as chance findings. Further research is required, but this study suggests that the genetic heterogeneity of MDD is not attributable to genome-wide moderation of genetic effects by CT.

Analysis of shared heritability in common disorders of the brain

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

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 GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia

BACKGROUND

Over the past decade genome-wide association studies (GWAS) have been applied to aid in the understanding of the biology of traits. The success of this approach is governed by the underlying effect sizes carried by the true risk variants and the corresponding statistical power to observe such effects given the study design and sample size under investigation. Previous ASD GWAS have identified genome-wide significant (GWS) risk loci; however, these studies were of only of low statistical power to identify GWS loci at the lower effect sizes (odds ratio (OR) <1.15).

METHODS

We conducted a large-scale coordinated international collaboration to combine independent genotyping data to improve the statistical power and aid in robust discovery of GWS loci. This study uses genome-wide genotyping data from a discovery sample (7387 ASD cases and 8567 controls) followed by meta-analysis of summary statistics from two replication sets (7783 ASD cases and 11359 controls; and 1369 ASD cases and 137308 controls).

RESULTS

We observe a GWS locus at 10q24.32 that overlaps several genes including PITX3, which encodes a transcription factor identified as playing a role in neuronal differentiation and CUEDC2 previously reported to be associated with social skills in an independent population cohort. We also observe overlap with regions previously implicated in schizophrenia which was further supported by a strong genetic correlation between these disorders (Rg = 0.23; P = 9 × 10^-6). We further combined these Psychiatric Genomics Consortium (PGC) ASD GWAS data with the recent PGC schizophrenia GWAS to identify additional regions which may be important in a common neurodevelopmental phenotype and identified 12 novel GWS loci. These include loci previously implicated in ASD such as FOXP1 at 3p13, ATP2B2 at 3p25.3, and a 'neurodevelopmental hub' on chromosome 8p11.23.

CONCLUSIONS

This study is an important step in the ongoing endeavour to identify the loci which underpin the common variant signal in ASD. In addition to novel GWS loci, we have identified a significant genetic correlation with schizophrenia and association of ASD with several neurodevelopmental-related genes such as EXT1, ASTN2, MACROD2, and HDAC4.