Estimates Of SNP Heritability Research Articles

Functional and molecular characterization of suicidality factors using phenotypic and genome-wide data

Genome-wide association studies (GWAS) of suicidal thoughts and behaviors support the existence of genetic contributions. Continuous measures of psychiatric disorder symptom severity can sometimes model polygenic risk better than binarized definitions. We compared two severity measures of suicidal thoughts and behaviors at the molecular and functional levels using genome-wide data. We used summary association data from GWAS of four traits analyzed in 122,935 individuals of European ancestry: thought life was not worth living (TLNWL), thoughts of self-harm, actual self-harm, and attempted suicide. A new trait for suicidal thoughts and behaviors was constructed first, phenotypically, by aggregating the previous four traits (termed “suicidality”) and second, genetically, by using genomic structural equation modeling (gSEM; termed S-factor). Suicidality and S-factor were compared using SNP-heritability (h2) estimates, genetic correlation (rg), partitioned h2, effect size distribution, transcriptomic correlations (ρGE) in the brain, and cross-population polygenic scoring (PGS). The S-factor had good model fit (χ2 = 0.21, AIC = 16.21, CFI = 1.00, SRMR = 0.024). Suicidality (h2 = 7.6%) had higher h2 than the S-factor (h2 = 2.54, Pdiff = 4.78 × 10−13). Although the S-factor had a larger number of non-null susceptibility loci (πc = 0.010), these loci had small effect sizes compared to those influencing suicidality (πc = 0.005, Pdiff = 0.045). The h2 of both traits was enriched for conserved biological pathways. The rg and ρGE support highly overlapping genetic and transcriptomic features between suicidality and the S-factor. PGS using European-ancestry SNP effect sizes strongly associated with TLNWL in Admixed Americans: Nagelkerke’s R2 = 8.56%, P = 0.009 (PGSsuicidality) and Nagelkerke’s R2 = 7.48%, P = 0.045 (PGSS-factor). An aggregate suicidality phenotype was statistically more heritable than the S-factor across all analyses and may be more informative for future genetic study designs interested in common genetic factors among different suicide related phenotypes.

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

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

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.

Assessing the contribution of genetic nurture to refractive error

Parents pass on both their genes and environment to offspring, prompting debate about the relative importance of nature versus nurture in the inheritance of complex traits. Advances in molecular genetics now make it possible to quantify an individual’s genetic predisposition to a trait via his or her ‘polygenic score’. However, part of the risk captured by an individual’s polygenic score may actually be attributed to the genotype of their parents. In the most well-studied example of this indirect ‘genetic nurture’ effect, about half the genetic contribution to educational attainment was found to be attributed to parental alleles, even if those alleles were not inherited by the child. Refractive errors, such as myopia, are a common cause of visual impairment and pose high economic and quality-of-life costs. Despite strong evidence that refractive errors are highly heritable, the extent to which genetic risk is conferred directly via transmitted risk alleles or indirectly via the environment that parents create for their children is entirely unknown. Here, an instrumental variable analysis in 1944 pairs of adult siblings from the United Kingdom was used to quantify the proportion of the genetic risk (‘single nucleotide polymorphism (SNP) heritability’) of refractive error contributed by genetic nurture. We found no evidence of a contribution from genetic nurture: non-within-family SNP-heritability estimate = 0.213 (95% confidence interval 0.134–0.310) and within-family SNP-heritability estimate = 0.250 (0.152–0.372). Our findings imply the genetic contribution to refractive error is principally an intrinsic effect from alleles transmitted from parents to offspring.

Shared genomic architectures of COVID-19 and antisocial behavior

Little is known about the genetics of norm violation and aggression in relation to coronavirus disease 2019 (COVID-19). To investigate this, we used summary statistics from genome-wide association studies and linkage disequilibrium score regression to calculate a matrix of genetic correlations (rgs) for antisocial behavior (ASB), COVID-19, and various health and behavioral traits. After false-discovery rate correction, ASB was genetically correlated with COVID-19 (rg = 0.51; P = 1.54E-02) and 19 other traits. ASB and COVID-19 were both positively genetically correlated with having a noisy workplace, doing heavy manual labor, chronic obstructive pulmonary disease, and genitourinary diseases. ASB and COVID-19 were both inversely genetically correlated with average income, education years, healthspan, verbal reasoning, lifespan, cheese intake, and being breastfed as a baby. But keep in mind that rgs are not necessarily causal. And, if causal, their prevailing directions of effect (which causes which) are indiscernible from rgs alone. Moreover, the SNP-heritability (h_g^2) estimates for two measures of COVID-19 were very small, restricting the overlap of genetic variance in absolute terms between ASB and COVID-19. Nonetheless, our findings suggest that those with antisocial tendencies possibly have a higher risk of exposure to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) than those without antisocial tendencies. This may have been especially true early in the pandemic before vaccines against SARS-CoV-2 were available and before the emergence of the highly transmissible Omicron variant.

Efficient estimation of SNP heritability using Gaussian predictive process in large scale cohort studies.

With the advent of high throughput genetic data, there have been attempts to estimate heritability from genome-wide SNP data on a cohort of distantly related individuals using linear mixed model (LMM). Fitting such an LMM in a large scale cohort study, however, is tremendously challenging due to its high dimensional linear algebraic operations. In this paper, we propose a new method named PredLMM approximating the aforementioned LMM motivated by the concepts of genetic coalescence and Gaussian predictive process. PredLMM has substantially better computational complexity than most of the existing LMM based methods and thus, provides a fast alternative for estimating heritability in large scale cohort studies. Theoretically, we show that under a model of genetic coalescence, the limiting form of our approximation is the celebrated predictive process approximation of large Gaussian process likelihoods that has well-established accuracy standards. We illustrate our approach with extensive simulation studies and use it to estimate the heritability of multiple quantitative traits from the UK Biobank cohort.

Probabilistic inference of the genetic architecture underlying functional enrichment of complex traits

We develop a Bayesian model (BayesRR-RC) that provides robust SNP-heritability estimation, an alternative to marker discovery, and accurate genomic prediction, taking 22 seconds per iteration to estimate 8.4 million SNP-effects and 78 SNP-heritability parameters in the UK Biobank. We find that only ≤10% of the genetic variation captured for height, body mass index, cardiovascular disease, and type 2 diabetes is attributable to proximal regulatory regions within 10kb upstream of genes, while 12-25% is attributed to coding regions, 32–44% to introns, and 22-28% to distal 10-500kb upstream regions. Up to 24% of all cis and coding regions of each chromosome are associated with each trait, with over 3,100 independent exonic and intronic regions and over 5,400 independent regulatory regions having ≥95% probability of contributing ≥0.001% to the genetic variance of these four traits. Our open-source software (GMRM) provides a scalable alternative to current approaches for biobank data.

Genetic aetiology of self-harm ideation and behaviour

Family studies have identified a heritable component to self-harm that is partially independent from comorbid psychiatric disorders. However, the genetic aetiology of broad sense (non-suicidal and suicidal) self-harm has not been characterised on the molecular level. In addition, controversy exists about the degree to which suicidal and non-suicidal self-harm share a common genetic aetiology. In the present study, we conduct genome-wide association studies (GWAS) on lifetime self-harm ideation and self-harm behaviour (i.e. any lifetime self-harm act regardless of suicidal intent) using data from the UK Biobank (n > 156,000). We also perform genome wide gene-based tests and characterize the SNP heritability and genetic correlations between these traits. Finally, we test whether polygenic risk scores (PRS) for self-harm ideation and self-harm behaviour predict suicide attempt, suicide thoughts and non-suicidal self-harm (NSSH) in an independent target sample of 8,703 Australian adults. Our GWAS results identified one genome-wide significant locus associated with each of the two phenotypes. SNP heritability (hsnp2) estimates were ~10%, and both traits were highly genetically correlated (LDSC rg > 0.8). Gene-based tests identified seven genes associated with self-harm ideation and four with self-harm behaviour. Furthermore, in the target sample, PRS for self-harm ideation were significantly associated with suicide thoughts and NSSH, and PRS for self-harm behaviour predicted suicide thoughts and suicide attempt. Follow up regressions identified a shared genetic aetiology between NSSH and suicide thoughts, and between suicide thoughts and suicide attempt. Evidence for shared genetic aetiology between NSSH and suicide attempt was not statistically significant.

Investigating the genetic and causal relationship between initiation or use of alcohol, caffeine, cannabis and nicotine

BackgroundCaffeine, alcohol, nicotine and cannabis are commonly used psychoactive substances. While the use of these substances has been previously shown to be genetically correlated, causality between these substance use traits remains unclear. We aimed to revisit the genetic relationships among different measures of SU using genome-wide association study (GWAS) summary statistics from the UK Biobank, International Cannabis Consortium, and GWAS & Sequencing Consortium of Alcohol and Nicotine use. MethodsWe obtained GWAS summary statistics from the aforementioned consortia for ten substance use traits including various measures of alcohol consumption, caffeine consumption, cannabis initiation and smoking behaviours. We then conducted SNP-heritability (h2) estimation for individual SU traits, followed by genetic correlation analyses and two-sample Mendelian randomisation (MR) studies between substance use trait pairs. ResultsSNP h2 of the ten traits ranged from 0.03 to 0.11. After multiple testing correction, 29 of the 45 trait pairs showed evidence of being genetically correlated. MR analyses revealed that most SU traits were not causally associated with each other. However, we found evidence for an MR association between regular smoking initiation and caffeine consumption 40.17 mg; 95 % CI: [24.01, 56.33] increase in caffeine intake per doubling of odds in smoking initiation). Our findings were robust against horizontal pleiotropy, SNP-outliers, and the direction of causality was consistent in all MR analyses. ConclusionsMost of the substance traits were genetically correlated but there is little evidence to establish causality apart from the relationship between smoking initiation and caffeine consumption.

Estimation of dynamic SNP-heritability with Bayesian Gaussian process models

MotivationImproved DNA technology has made it practical to estimate single-nucleotide polymorphism (SNP)-heritability among distantly related individuals with unknown relationships. For growth- and development-related traits, it is meaningful to base SNP-heritability estimation on longitudinal data due to the time-dependency of the process. However, only few statistical methods have been developed so far for estimating dynamic SNP-heritability and quantifying its full uncertainty.ResultsWe introduce a completely tuning-free Bayesian Gaussian process (GP)-based approach for estimating dynamic variance components and heritability as their function. For parameter estimation, we use a modern Markov Chain Monte Carlo method which allows full uncertainty quantification. Several datasets are analysed and our results clearly illustrate that the 95% credible intervals of the proposed joint estimation method (which ‘borrows strength’ from adjacent time points) are significantly narrower than of a two-stage baseline method that first estimates the variance components at each time point independently and then performs smoothing. We compare the method with a random regression model using MTG2 and BLUPF90 software and quantitative measures indicate superior performance of our method. Results are presented for simulated and real data with up to 1000 time points. Finally, we demonstrate scalability of the proposed method for simulated data with tens of thousands of individuals.Availability and implementationThe C++ implementation dynBGP and simulated data are available in GitHub: https://github.com/aarjas/dynBGP. The programmes can be run in R. Real datasets are available in QTL archive: https://phenome.jax.org/centers/QTLA.Supplementary information Supplementary data are available at Bioinformatics online.

Genomic Prediction Using Low Density Marker Panels in Aquaculture: Performance Across Species, Traits, and Genotyping Platforms.

Genomic selection increases the rate of genetic gain in breeding programs, which results in significant cumulative improvements in commercially important traits such as disease resistance. Genomic selection currently relies on collecting genome-wide genotype data accross a large number of individuals, which requires substantial economic investment. However, global aquaculture production predominantly occurs in small and medium sized enterprises for whom this technology can be prohibitively expensive. For genomic selection to benefit these aquaculture sectors, more cost-efficient genotyping is necessary. In this study the utility of low and medium density SNP panels (ranging from 100 to 9,000 SNPs) to accurately predict breeding values was tested and compared in four aquaculture datasets with different characteristics (species, genome size, genotyping platform, family number and size, total population size, and target trait). The traits show heritabilities between 0.19–0.49, and genomic prediction accuracies using the full density panel of 0.55–0.87. A consistent pattern of genomic prediction accuracy was observed across species with little or no accuracy reduction until SNP density was reduced below 1,000 SNPs (prediction accuracies of 0.44–0.75). Below this SNP density, heritability estimates and genomic prediction accuracies tended to be lower and more variable (93% of maximum accuracy achieved with 1,000 SNPs, 89% with 500 SNPs, and 70% with 100 SNPs). A notable drop in accuracy was observed between 200 SNP panels (0.44–0.75) and 100 SNP panels (0.39–0.66). Now that a multitude of studies have highlighted the benefits of genomic over pedigree-based prediction of breeding values in aquaculture species, the results of the current study highlight that these benefits can be achieved at lower SNP densities and at lower cost, raising the possibility of a broader application of genetic improvement in smaller and more fragmented aquaculture settings.

The behavioral, cellular and immune mediators of HIV-1 acquisition: New insights from population genetics

Millions are exposed to the human immunodeficiency virus type 1 (HIV-1) every year, but not all acquire the virus, suggesting a potential role for host genetics in the moderation of HIV-1 acquisition. Here, we analyzed summary statistics from the largest genome-wide association study of HIV-1 acquisition to-date, consisting of 6,334 infected patients and 7,247 population controls, to advance our understanding of the genetic mechanisms implicated in this trait. We found that HIV-1 acquisition is polygenic and heritable, with SNP heritability estimates explaining 28–42% of the variance in this trait at a population level. Genetic correlations alongside UK Biobank data revealed associations with smoking, prospective memory and socioeconomic traits. Gene-level enrichment analysis identified EF-hand calcium binding domain 14 as a novel susceptibility gene for HIV–1 acquisition. We also observed that susceptibility variants for HIV-1 acquisition were significantly enriched for genes expressed in T-cells, but also in striatal and hippocampal neurons. Finally, we tested how polygenic risk scores for HIV-1 acquisition influence blood levels of 35 inflammatory markers in 406 HIV-1-negative individuals. We found that higher genetic risk for HIV-1 acquisition was associated with lower levels of C-C motif chemokine ligand 17. Our findings corroborate a complex model for HIV-1 acquisition, whereby susceptibility is partly heritable and moderated by specific behavioral, cellular and immunological parameters.

Novel genome-wide associations for anhedonia, genetic correlation with psychiatric disorders, and polygenic association with brain structure

Anhedonia is a core symptom of several psychiatric disorders but its biological underpinnings are poorly understood. We performed a genome-wide association study of state anhedonia in 375,275 UK Biobank participants and assessed for genetic correlation between anhedonia and neuropsychiatric conditions (major depressive disorder, schizophrenia, bipolar disorder, obsessive compulsive disorder and Parkinson’s Disease). We then used a polygenic risk score approach to test for association between genetic loading for anhedonia and both brain structure and brain function. This included: magnetic resonance imaging (MRI) assessments of total grey matter volume, white matter volume, cerebrospinal fluid volume, and 15 cortical/subcortical regions of interest; diffusion tensor imaging (DTI) measures of white matter tract integrity; and functional MRI activity during an emotion processing task. We identified 11 novel loci associated at genome-wide significance with anhedonia, with a SNP heritability estimate (h2SNP) of 5.6%. Strong positive genetic correlations were found between anhedonia and major depressive disorder, schizophrenia and bipolar disorder; but not with obsessive compulsive disorder or Parkinson’s Disease. Polygenic risk for anhedonia was associated with poorer brain white matter integrity, smaller total grey matter volume, and smaller volumes of brain regions linked to reward and pleasure processing, including orbito-frontal cortex. In summary, the identification of novel anhedonia-associated loci substantially expands our current understanding of the biological basis of state anhedonia and genetic correlations with several psychiatric disorders confirm the utility of this phenotype as a transdiagnostic marker of vulnerability to mental illness. We also provide the first evidence that genetic risk for state anhedonia influences brain structure, including in regions associated with reward and pleasure processing.

Contribution of Common Genetic Variants to Familial Aggregation of Disease and Implications for Sequencing Studies.

Despite genetics being accepted as the primary cause of familial aggregation for most diseases, it is still unclear whether afflicted families are likely to share a single highly penetrant rare variant, many minimally penetrant common variants, or a combination of the two types of variants. We therefore use recent estimates of SNP heritability and the liability threshold model to estimate the proportion of afflicted families likely to carry a rare, causal variant. We then show that Polygenic Risk Scores (PRS) may be useful for identifying families likely to carry such a rare variant and therefore for prioritizing families to include in sequencing studies with that aim. Specifically, we introduce a new statistic that estimates the proportion of individuals carrying causal rare variants based on the family structure, disease pattern, and PRS of genotyped individuals. Finally, we consider data from the MelaNostrum consortium and show that, despite an estimated PRS heritability of only 0.05 for melanoma, families carrying putative causal variants had a statistically significantly lower PRS, supporting the idea that PRS prioritization may be a useful future tool. However, it will be important to evaluate whether the presence of rare mendelian variants are generally associated with the proposed test statistic or lower PRS in future and larger studies.

Accurate estimation of SNP-heritability from biobank-scale data irrespective of genetic architecture.

SNP-heritability is a fundamental quantity in the study of complex traits. Recent works have shown that existing methods to estimate genome-wide SNP-heritability yield biases when their assumptions are violated. While various approaches have been proposed to account for frequency- and LD-dependent genetic architectures, it remains unclear which estimates reported in the literature are reliable. Here we show that genome-wide SNP-heritability can be accurately estimated from biobank-scale data irrespective of genetic architecture, without specifying a heritability model or partitioning SNPs by allele frequency and/or LD. We show analytically and through extensive simulations starting from real genotypes (UK Biobank, N = 337K) that, unlike existing methods, our closed-form estimator is robust across a wide range of architectures. We provide estimates of SNP-heritability for 22 complex traits in the UK Biobank and show that, consistent with our results in simulations, existing biobank-scale methods yield estimates up to 30% different from our theoretically-justified approach.

The genomic basis of mood instability: identification of 46 loci in 363,705 UK Biobank participants, genetic correlation with psychiatric disorders, and association with gene expression and function.

Genome-wide association studies (GWAS) of psychiatric phenotypes have tended to focus on categorical diagnoses, but to understand the biology of mental illness it may be more useful to study traits which cut across traditional boundaries. Here we report the results of a GWAS of mood instability as a trait in a large population cohort (UK Biobank, n = 363,705). We also assess the clinical and biological relevance of the findings, including whether genetic associations show enrichment for nervous system pathways. Forty six unique loci associated with mood instability were identified with a SNP heritability estimate of 9%. Linkage Disequilibrium Score Regression (LDSR) analyses identified genetic correlations with Major Depressive Disorder (MDD), Bipolar Disorder (BD), Schizophrenia, anxiety and Post Traumatic Stress Disorder (PTSD). Gene-level and gene set analyses identified 244 significant genes and 6 enriched gene sets. Tissue expression analysis of the SNP-level data found enrichment in multiple brain regions, and eQTL analyses highlighted an inversion on chromosome 17 plus two brain-specific eQTLs. Additionally, we used a Phenotype Linkage Network (PLN) analysis and community analysis to assess for enrichment of nervous system gene sets using mouse orthologue databases. The PLN analysis found enrichment in nervous system PLNs for a community containing serotonin and melatonin receptors. In summary, this work has identified novel loci, tissues and gene sets contributing to mood instability. These findings may be relevant for the identification of novel trans-diagnostic drug targets and could help to inform future precision medicine innovations in mental health.

Overlapping genetic effects between suicidal ideation and neurocognitive functioning

BackgroundSuicide is the second leading cause of death among adolescents and young adults. Several studies have indicated significant genetic influences on suicide-related phenotypes and mounting evidence from neurobiological research has linked deficits in neurocognitive abilities to suicide phenotypes. The goal of the present study was to estimate the heritability of suicidal ideation (SI) in a large sample of adolescents and determine if SI is genetically correlated with neurocognitive functioning. MethodsGenome-wide data (N = 3564 unrelated individuals of European Ancestry) were drawn from the Philadelphia Neurodevelopment Cohort. Adolescents completed a psychiatric assessment, as well as a computerized neurocognitive battery to assess performance across four domains: memory, executive function, social cognition, and complex cognition. Genomic-relatedness-matrix restricted maximum likelihood (GREML) estimation was used to determine SNP-heritability (h2SNP) of SI and the genetic correlation (rG) between SI and neurocognitive domains. ResultsNearly 17% of adolescents reported SI. The SNP-heritability estimate for SI was marginally significant (h2SNP = 11%, SE = 8%, p = 0.086). Bivariate analyses indicated a significant rG between SI and emotion identification (rG = 0.79, SE = 0.45, p = 0.006; phenotypic correlation r = 0.04, p = 0.017). LimitationsIt is possible that SI may represent a related, but differentially heritable construct from suicide attempts/completion and other comorbid psychopathology. Additionally, though genetic correlations point to shared genetic factors across traits, direct causal mechanisms cannot be deduced. ConclusionsCommon heritable factors contribute to variation in SI and neurocognitive functioning. Genetic factors influencing emotion identification have significant genetic overlap with SI.

66 - GENOME-WIDE ASSOCIATION STUDY OF ANXIETY AND STRESS-RELATED DISORDERS

Background In Europe, over 70 million people are annually suffering from anxiety and stress-related disorders and over 90% of patients with these disorders have another comorbid mental or somatic disorder. Therefore, it is of strong interest to identify the risk factors of these disorders. Heritability estimates vary between 30% and 50%, however first Genome-Wide Association Studies (GWAS) have been limited by their small sample sizes. As anxiety and stress-related disorders are likely to configure various expressions of abnormalities in the basic stress-response system, we aimed to conduct a GWAS aggregating a large number of cases with varying diagnoses of anxiety and stress-related disorders. Methods The study is making use of a Danish nation-wide population-based sample collected by the Lundbeck Foundation Initiative for Integrative Psychiatric Research including 14,711 cases with anxiety and stress-related diagnoses and 17,857 controls. Data were processed using the Ricopili pipeline; imputation is based on the 1000 genomes phase 3 as reference panel. GWAS was conducted using the imputed marker dosages and an additive logistic regression model. Analyses were supplemented by comparing severe mental disorder with comorbid anxiety and stress-related disorder to those without, using a propensity-score-matched design. SNP heritability and genetic correlation with other traits were computed using LD-score regression. Results In case-control design several SNPs located in the PDE4B gene achieved genome-wide significance. The analyses comparing mental disorders with and without comorbid anxiety and stress-related disorders reveled no significant signal. The SNP heritability estimate for anxiety was 14%, 33% for stress-related disorders and 26% for the combined phenotype. We found a positive significant genetic correlation of anxiety and stress-related disorders with neuroticism, depression, bipolar disorder, schizophrenia and smoking; and a negative genetic correlation with college completion, years of schooling, age at menarche and age of first birth. Discussion This is the first GWAS meta-analysis on anxiety and stress-related disorders demonstrating influences of common genetic variation in the etiology of these disorders and indicating PDE4B as a susceptibility locus. Consistent with phenotypic observations, anxiety and stress-related disorders share genetic variants with psychiatric traits in the expected direction.

OCMA: Fast, Memory-Efficient Factorization of Prohibitively Large Relationship Matrices.

Matrices representing genetic relatedness among individuals (i.e., Genomic Relationship Matrices, GRMs) play a central role in genetic analysis. The eigen-decomposition of GRMs (or its alternative that generates fewer top singular values using genotype matrices) is a necessary step for many analyses including estimation of SNP-heritability, Principal Component Analysis (PCA), and genomic prediction. However, the GRMs and genotype matrices provided by modern biobanks are too large to be stored in active memory. To accommodate the current and future “bigger-data”, we develop a disk-based tool, Out-of-Core Matrices Analyzer (OCMA), using state-of-the-art computational techniques that can nimbly perform eigen and Singular Value Decomposition (SVD) analyses. By integrating memory mapping (mmap) and the latest matrix factorization libraries, our tool is fast and memory-efficient. To demonstrate the impressive performance of OCMA, we test it on a personal computer. For full eigen-decomposition, it solves an ordinary GRM (N = 10,000) in 55 sec. For SVD, a commonly used faster alternative of full eigen-decomposition in genomic analyses, OCMA solves the top 200 singular values (SVs) in half an hour, top 2,000 SVs in 0.95 hr, and all 5,000 SVs in 1.77 hr based on a very large genotype matrix (N = 1,000,000, M = 5,000) on the same personal computer. OCMA also supports multi-threading when running in a desktop or HPC cluster. Our OCMA tool can thus alleviate the computing bottleneck of classical analyses on large genomic matrices, and make it possible to scale up current and emerging analytical methods to big genomics data using lightweight computing resources.

Extracting stability increases the SNP heritability of emotional problems in young people

Twin studies have shown that emotional problems (anxiety and depression) in childhood and adolescence are moderately heritable (~20–50%). In contrast, DNA-based ‘SNP heritability’ estimates are generally <15% and non-significant. One notable feature of emotional problems is that they can be somewhat transient, but the moderate stability seen across time and across raters is predominantly influenced by stable genetic influences. This suggests that by capturing what is in common across time and across raters, we might be more likely to tap into any underlying genetic vulnerability. We therefore hypothesised that a phenotype capturing the pervasive stability of emotional problems would show higher heritability. We fitted single-factor latent trait models using 12 emotional problems measures across ages 7, 12 and 16, rated by parents, teachers and children themselves in the Twins Early Development Study sample. Twin and SNP heritability estimates for stable emotional problems (N = 6110 pairs and 6110 unrelated individuals, respectively) were compared to those for individual measures. Twin heritability increased from 45% on average for individual measures to 76% (se = 0.023) by focusing on stable trait variance. SNP heritability rose from 5% on average (n.s.) to 14% (se = 0.049; p = 0.002). Heritability was also higher for stable within-rater composites. Polygenic scores for both adult anxiety and depression significantly explained variance in stable emotional problems (0.4%; p = 0.0001). The variance explained was more than in most individual measures. Stable emotional problems also showed significant genetic correlation with adult depression and anxiety (average = 52%). These results demonstrate the value of examining stable emotional problems in gene-finding and prediction studies.