Disentangling potential causal effects of educational duration on well-being, and mental and physical health outcomes

BACKGROUND

Extensive research has focused on the potential benefits of education on various mental and physical health outcomes. However, whether the associations reflect a causal effect is harder to establish.

METHODS

To examine associations between educational duration and specific aspects of well-being, anxiety and mood disorders, and cardiovascular health in a sample of European Ancestry UK Biobank participants born in England and Wales, we apply four different causal inference methods (a natural policy experiment leveraging the minimum school-leaving age, a sibling-control design, Mendelian randomization [MR], and within-family MR), and assess if the methods converge on the same conclusion.

RESULTS

A comparison of results across the four methods reveals that associations between educational duration and these outcomes appears predominantly to be the result of confounding or bias rather than a true causal effect of education on well-being and health outcomes. Although we do consistently find no associations between educational duration and happiness, family satisfaction, work satisfaction, meaning in life, anxiety, and bipolar disorder, we do not find consistent significant associations across all methods for the other phenotypes (health satisfaction, depression, financial satisfaction, friendship satisfaction, neuroticism, and cardiovascular outcomes).

CONCLUSIONS

We discuss inconsistencies in results across methods considering their respective limitations and biases, and additionally discuss the generalizability of our findings in light of the sample and phenotype limitations. Overall, this study strengthens the idea that triangulation across different methods is necessary to enhance our understanding of the causal consequences of educational duration.

Evaluating the causal relationship between educational attainment and mental health

We investigate the causal relationship between educational attainment (EA) and mental health using two research designs. First, we compare the relationship between EA and 18 psychiatric diagnoses within sibship in Dutch national registry data (N=1.7 million), thereby controlling for unmeasured familial factors. Second, we apply two-sample Mendelian Randomization, which uses genetic variants related to EA or psychiatric diagnosis as instrumental variables, to test whether there is a causal relation in either direction. Our results suggest that lower levels of EA causally increase the risk of MDD, ADHD, alcohol dependence, GAD and PTSD diagnoses. We also find evidence of a causal effect of ADHD on EA. For schizophrenia, anorexia nervosa, OCD, and bipolar disorder, results were inconsistent across the different approaches, highlighting the importance of using multiple research designs to understand complex relationships such as between EA and mental health.

Using Polygenic Scores for Circadian Rhythms to Predict Wellbeing, Depressive Symptoms, Chronotype, and Health

The association between circadian rhythms and diseases has been well established, while the association with mental health is less explored. Given the heritable nature of circadian rhythms, this study aimed to investigate the relationship between genes underlying circadian rhythms and mental health outcomes, as well as a possible gene-environment correlation for circadian rhythms. Polygenic scores (PGSs) represent the genetic predisposition to develop a certain trait or disease. In a sample from the Netherlands Twin Register (N = 14,021), PGSs were calculated for two circadian rhythm measures: morningness and relative amplitude (RA). The PGSs were used to predict mental health outcomes such as subjective happiness, quality of life, and depressive symptoms. In addition, we performed the same prediction analysis in a within-family design in a subset of dizygotic twins. The PGS for morningness significantly predicted morningness (R2 = 1.55%) and depressive symptoms (R2 = 0.22%). The PGS for RA significantly predicted general health (R2 = 0.12%) and depressive symptoms (R2 = 0.20%). Item analysis of the depressive symptoms showed that 4 out of 14 items were significantly associated with the PGSs. Overall, the results showed that people with a genetic predisposition of being a morning person or with a high RA are likely to have fewer depressive symptoms. The four associated depressive symptoms described symptoms related to decision-making, energy, and feeling worthless or inferior, rather than sleep. Based on our findings future research should include a substantial role for circadian rhythms in depression research and should further explore the gene-environment correlation in circadian rhythms.

Beyond the factor indeterminacy problem using genome-wide association data

Latent factors, such as general intelligence, depression and risk tolerance, are invoked in nearly all social science research where a construct is measured via aggregation of symptoms, question responses or other measurements. Because latent factors cannot be directly observed, they are inferred by fitting a specific model to empirical patterns of correlations among measured variables. A long-standing critique of latent factor theories is that the correlations used to infer latent factors can be produced by alternative data-generating mechanisms that do not include latent factors. This is referred to as the factor indeterminacy problem. Researchers have recently begun to overcome this problem by using information on the associations between individual genetic variants and measured variables. We review historical work on the factor indeterminacy problem and describe recent efforts in genomics to rigorously test the validity of latent factors, advancing the understanding of behavioural science constructs.

More than nature and nurture, indirect genetic effects on children's academic achievement are consequences of dynastic social processes

Families transmit genes and environments across generations. When parents' genetics affect their children's environments, these two modes of inheritance can produce an 'indirect genetic effect'. Such indirect genetic effects may account for up to half of the estimated genetic variance in educational attainment. Here we tested if indirect genetic effects reflect within-nuclear-family transmission ('genetic nurture') or instead a multi-generational process of social stratification ('dynastic effects'). We analysed indirect genetic effects on children's academic achievement in their fifth to ninth years of schooling in N = 37,117 parent-offspring trios in the Norwegian Mother, Father, and Child Cohort Study (MoBa). We used pairs of genetically related families (parents were siblings, children were cousins; N = 10,913) to distinguish within-nuclear-family genetic-nurture effects from dynastic effects shared by cousins in different nuclear families. We found that indirect genetic effects on children's academic achievement cannot be explained by processes that operate exclusively within the nuclear family.

The shared genetic architecture and evolution of human language and musical rhythm

Rhythm and language-related traits are phenotypically correlated, but their genetic overlap is largely unknown. Here, we leveraged two large-scale genome-wide association studies performed to shed light on the shared genetics of rhythm (N=606,825) and dyslexia (N=1,138,870). Our results reveal an intricate shared genetic and neurobiological architecture, and lay groundwork for resolving longstanding debates about the potential co-evolution of human language and musical traits.

Genetic contributions of noncognitive skills to academic development

Noncognitive skills such as motivation and self-regulation, are partly heritable and predict academic achievement beyond cognitive skills. However, how the relationship between noncognitive skills and academic achievement changes over development is unclear. The current study examined how cognitive and noncognitive skills contribute to academic achievement from ages 7 to 16 in a sample of over 10,000 children from England and Wales. Noncognitive skills were increasingly predictive of academic achievement across development. Twin and polygenic scores analyses found that the contribution of noncognitive genetics to academic achievement became stronger over the school years. Results from within-family analyses indicated that associations with noncognitive genetics could not simply be attributed to confounding by environmental differences between nuclear families and are consistent with a possible role for evocative/active gene-environment correlations. By studying genetic effects through a developmental lens, we provide novel insights into the role of noncognitive skills in academic development.

Genetic analysis of selection bias in a natural experiment: Investigating in-utero famine effects on elevated body mass index in the Dutch Hunger Winter Families Study

Natural-experiment designs that compare survivors of in-utero famine exposure to unaffected controls suggest that in-utero undernutrition predisposes to development of obesity. However, birth rates drop dramatically during famines. Selection bias could arise if factors that contribute to obesity also protect fertility and/or fetal survival under famine conditions. We investigated this hypothesis using genetic analysis of a famine-exposed birth cohort. We genotyped participants in the Dutch Hunger Winter Families Study (DHWFS, N=950; 45% male), of whom 51% were exposed to the 1944-1945 Dutch Famine during gestation and 49% were their unexposed same-sex siblings or "time controls" born before or after the famine in the same hospitals. We computed body-mass index (BMI) polygenic indices (PGIs) in DHWFS participants and compared BMI PGIs between famine-exposed and control groups. Participants with higher polygenic risk had higher BMIs (Pearson r=0.42, p<0.001). However, differences between BMI PGIs of famine-exposed participants and controls were small and not statistically different from zero across specifications (Cohen's d=0.10, p>0.092). Our findings did not indicate selection bias, supporting the validity of the natural-experiment design within DHWFS. In summary, our study outlines a novel approach to explore the presence of selection bias in famine and other natural experiment studies.

Educational attainment, health outcomes and mortality: a within-sibship Mendelian randomization study

BACKGROUND

Previous Mendelian randomization (MR) studies using population samples (population MR) have provided evidence for beneficial effects of educational attainment on health outcomes in adulthood. However, estimates from these studies may have been susceptible to bias from population stratification, assortative mating and indirect genetic effects due to unadjusted parental genotypes. MR using genetic association estimates derived from within-sibship models (within-sibship MR) can avoid these potential biases because genetic differences between siblings are due to random segregation at meiosis.

METHODS

Applying both population and within-sibship MR, we estimated the effects of genetic liability to educational attainment on body mass index (BMI), cigarette smoking, systolic blood pressure (SBP) and all-cause mortality. MR analyses used individual-level data on 72 932 siblings from UK Biobank and the Norwegian HUNT study, and summary-level data from a within-sibship Genome-wide Association Study including >140 000 individuals.

RESULTS

Both population and within-sibship MR estimates provided evidence that educational attainment decreased BMI, cigarette smoking and SBP. Genetic variant-outcome associations attenuated in the within-sibship model, but genetic variant-educational attainment associations also attenuated to a similar extent. Thus, within-sibship and population MR estimates were largely consistent. The within-sibship MR estimate of education on mortality was imprecise but consistent with a putative effect.

CONCLUSIONS

These results provide evidence of beneficial individual-level effects of education (or liability to education) on adulthood health, independently of potential demographic and family-level confounders.

Direct and Indirect Genetic Effects on Aggression

Background

Family members resemble each other in their propensity for aggression. In twin studies, approximately 50% of the variance in aggression can be explained by genetic influences. However, if there are genotype-environment correlation mechanisms, such as environmental manifestations of parental and sibling genotypes, genetic influences may partly reflect environmental influences. In this study, we investigated the importance of indirect polygenic score (PGS) effects on aggression.

Methods

We modeled the effect of PGSs based on 3 genome-wide association studies: early-life aggression, educational attainment, and attention-deficit/hyperactivity disorder (ADHD). The associations with aggression were tested in a within- and between-family design (37,796 measures from 7740 individuals, ages 3-86 years [mean = 14.20 years, SE = 12.03], from 3107 families, 55% female) and in a transmitted/nontransmitted PGS design (42,649 measures from 6653 individuals, ages 3-61 years [mean = 11.81 years, SE = 8.68], from 3024 families, 55% female). All participants are enrolled in the Netherlands Twin Register.

Results

We found no evidence for contributions of indirect PGS effects on aggression in either a within- and between-family design or a transmitted/nontransmitted PGS design. Results indicate significant direct effects on aggression for the PGSs based on early-life aggression, educational attainment, and ADHD, although explained variance was low (within- and between-family: early-life aggression R2 = 0.3%, early-life ADHD R2 = 0.6%, educational attainment R2 = 0.7%; transmitted/nontransmitted PGSs: early-life aggression R2 = 0.2%, early-life ADHD R2 = 0.9%, educational attainment R2 = 0.5%).

Conclusions

PGSs included in the current study had a direct (but no indirect) effect on aggression, consistent with results of previous twin and family studies. Further research involving other PGSs for aggression and related phenotypes is needed to determine whether this conclusion generalizes to overall genetic influences on aggression.

Patterns of item nonresponse behaviour to survey questionnaires are systematic and associated with genetic loci

Response to survey questionnaires is vital for social and behavioural research, and most analyses assume full and accurate response by participants. However, nonresponse is common and impedes proper interpretation and generalizability of results. We examined item nonresponse behaviour across 109 questionnaire items in the UK Biobank (N = 360,628). Phenotypic factor scores for two participant-selected nonresponse answers, 'Prefer not to answer' (PNA) and 'I don't know' (IDK), each predicted participant nonresponse in follow-up surveys (incremental pseudo-R2 = 0.056), even when controlling for education and self-reported health (incremental pseudo-R2 = 0.046). After performing genome-wide association studies of our factors, PNA and IDK were highly genetically correlated with one another (rg = 0.73 (s.e. = 0.03)) and with education (rg,PNA = -0.51 (s.e. = 0.03); rg,IDK = -0.38 (s.e. = 0.02)), health (rg,PNA = 0.51 (s.e. = 0.03); rg,IDK = 0.49 (s.e. = 0.02)) and income (rg,PNA = -0.57 (s.e. = 0.04); rg,IDK = -0.46 (s.e. = 0.02)), with additional unique genetic associations observed for both PNA and IDK (P < 5 × 10-8). We discuss how these associations may bias studies of traits correlated with item nonresponse and demonstrate how this bias may substantially affect genome-wide association studies. While the UK Biobank data are deidentified, we further protected participant privacy by avoiding exploring non-response behaviour to single questions, assuring that no information can be used to associate results with any particular respondents.

Genetic structure of major depression symptoms across clinical and community cohorts

Diagnostic criteria for major depressive disorder allow for heterogeneous symptom profiles but genetic analysis of major depressive symptoms has the potential to identify clinical and aetiological subtypes. There are several challenges to integrating symptom data from genetically-informative cohorts, such as sample size differences between clinical and community cohorts and various patterns of missing data. We conducted genome-wide association studies of major depressive symptoms in three clinical cohorts that were enriched for affected participants (Psychiatric Genomics Consortium, Australian Genetics of Depression Study, Generation Scotland) and three community cohorts (Avon Longitudinal Study of Parents and Children, Estonian Biobank, and UK Biobank). We fit a series of confirmatory factor models with factors that accounted for how symptom data was sampled and then compared alternative models with different symptom factors. The best fitting model had a distinct factor for Appetite/Weight symptoms and an additional measurement factor that accounted for missing data patterns in the community cohorts (use of Depression and Anhedonia as gating symptoms). The results show the importance of assessing the directionality of symptoms (such as hypersomnia versus insomnia) and of accounting for study and measurement design when meta-analysing genetic association data.

Genetic contributions of noncognitive skills to academic development

Noncognitive skills such as motivation and self-regulation, predict academic achievement beyond cognitive skills. However, the role of genetic and environmental factors and of their interplay in these developmental associations remains unclear. We provide a comprehensive account of how cognitive and noncognitive skills contribute to academic achievement from ages 7 to 16 in a sample of >10,000 children from England and Wales. Results indicated that noncognitive skills become increasingly predictive of academic achievement across development. Triangulating genetic methods, including twin analyses and polygenic scores (PGS), we found that the contribution of noncognitive genetics to academic achievement becomes stronger over development. The PGS for noncognitive skills predicted academic achievement developmentally, with prediction nearly doubling by age 16, pointing to gene-environment correlation (rGE). Within-family analyses indicated both passive and active/evocative rGE processes driven by noncognitive genetics. By studying genetic effects through a developmental lens, we provide novel insights into the role of noncognitive skills in academic development.

The Relationship of Attention-Deficit/Hyperactivity Disorder With Posttraumatic Stress Disorder: A Two-Sample Mendelian Randomization and Population-Based Sibling Comparison Study

BACKGROUND

Attention-deficit/hyperactivity disorder (ADHD) and posttraumatic stress disorder (PTSD) are associated, but it is unclear if this is a causal relationship or confounding. We used genetic analyses and sibling comparisons to clarify the direction of this relationship.

METHODS

Linkage disequilibrium score regression and 2-sample Mendelian randomization were used to test for genetic correlation (rg) and bidirectional causal effects using European ancestry genome-wide association studies of ADHD (20,183 cases and 35,191 controls) and 6 PTSD definitions (up to 320,369 individuals). Several additional variables were included in the analysis to verify the independence of the ADHD-PTSD relationship. In a population-based sibling comparison (N = 2,082,118 individuals), Cox regression models were fitted to account for time at risk, a range of sociodemographic factors, and unmeasured familial confounders (via sibling comparisons).

RESULTS

ADHD and PTSD had consistent rg (rg range, 0.43-0.52; p <  .001). ADHD genetic liability was causally linked with increased risk for PTSD (β = 0.367; 95% CI, 0.186-0.552; p = 7.68 × 10-5). This result was not affected by heterogeneity, horizontal pleiotropy (Mendelian randomization Egger intercept = 4.34 × 10-4, p = .961), or other phenotypes and was consistent across PTSD datasets. However, we found no consistent associations between PTSD genetic liability and ADHD risk. Individuals diagnosed with ADHD were at a higher risk for developing PTSD than their undiagnosed sibling (hazard ratio = 2.37; 95% CI, 1.98-3.53).

CONCLUSIONS

Our findings add novel evidence supporting the need for early and effective treatment of ADHD, as patients with this diagnosis are at significantly higher risk to develop PTSD later in life.

Plasma biomarkers predict amyloid pathology in cognitively normal monozygotic twins after 10 years

Blood-based biomarkers could prove useful to predict Alzheimer's disease core pathologies in advance of clinical symptoms. Implementation of such biomarkers requires a solid understanding of their long-term dynamics and the contribution of confounding to their association with Alzheimer's disease pathology. Here we assess the value of plasma amyloid-β_1-42/1-40, phosphorylated-tau181 and glial fibrillary acidic protein to detect early Alzheimer's disease pathology, accounting for confounding by genetic and early environmental factors. Participants were 200 monozygotic twins, aged ≥60 years with normal cognition from the european medical information framework for Alzheimer's disease study. All twins had amyloid-β status and plasma samples available at study enrolment. For 80 twins, additional plasma samples were available that had been collected approximately 10 years prior to amyloid-β status assessment. Single-molecule array assays were applied to measure amyloid-β_1-42/1-40, phosphorylated-tau181 and glial fibrillary acidic protein. Predictive value of and longitudinal change in these biomarkers were assessed using receiver operating characteristic curve analysis and linear mixed models. Amyloid pathology could be predicted using blood-based biomarkers obtained at the time of amyloid status assessment (amyloid-β_1-42/1-40: area under the curve = 0.65, P = 0.01; phosphorylated-tau181: area under the curve = 0.84, P < 0.001; glial fibrillary acidic protein: area under the curve = 0.74, P < 0.001), as well as using those obtained 10 years prior to amyloid status assessment (amyloid-β_1-42/1-40: area under the curve = 0.69, P = 0.03; phosphorylated-tau181: area under the curve = 0.92, P < 0.001; glial fibrillary acidic protein: area under the curve = 0.84, P < 0.001). Longitudinally, amyloid-β_1-42/1-40 levels decreased [β (SE) = -0.12 (0.01), P < 0.001] and phosphorylated-tau181 levels increased [β (SE) = 0.02 (0.01), P = 0.004]. Amyloid-β-positive individuals showed a steeper increase in phosphorylated-tau181 compared with amyloid-β-negative individuals [β (SE) = 0.06 (0.02), P = 0.004]. Also amyloid-β-positive individuals tended to show a steeper increase in glial fibrillary acidic protein [β (SE) = 0.04 (0.02), P = 0.07]. Within monozygotic twin pairs, those with higher plasma phosphorylated-tau181 and lower amyloid-β_1-42/1-40 levels were more likely to be amyloid-β positive [β (SE) = 0.95 (0.26), P < 0.001; β (SE) = -0.28 (0.14), P < 0.05] indicating minimal contribution of confounding by genetic and early environmental factors. Our data support the use of amyloid-β_1-42/1-40, phosphorylated-tau181 and glial fibrillary acidic protein as screening tools for Alzheimer's disease pathology in the normal aging population, which is of importance for enrolment of high-risk subjects in secondary, or even primary, prevention trials. Furthermore, these markers show potential as low-invasive monitoring tool of disease progression and possibly treatment effects in clinical trials.

Multivariate analyses of molecular genetic associations between childhood psychopathology and adult mood disorders and related traits

Ubiquitous associations have been detected between different types of childhood psychopathology and polygenic risk scores based on adult psychiatric disorders and related adult outcomes, indicating that genetic factors partly explain the association between childhood psychopathology and adult outcomes. However, these analyses in general do not take into account the correlations between the adult trait and disorder polygenic risk scores. This study aimed to further clarify the influence of genetic factors on associations between childhood psychopathology and adult outcomes by accounting for these correlations. Using a multivariate multivariable regression, we analyzed associations of childhood attention-deficit/hyperactivity disorder (ADHD), internalizing, and social problems, with polygenic scores (PGS) of adult disorders and traits including major depression, bipolar disorder, subjective well-being, neuroticism, insomnia, educational attainment, and body mass index (BMI), derived for 20,539 children aged 8.5-10.5 years. After correcting for correlations between the adult phenotypes, major depression PGS were associated with all three childhood traits, that is, ADHD, internalizing, and social problems. In addition, BMI PGS were associated with ADHD symptoms and social problems, while neuroticism PGS were only associated with internalizing problems and educational attainment PGS were only associated with ADHD symptoms. PGS of bipolar disorder, subjective well-being, and insomnia were not associated with any childhood traits. Our findings suggest that associations between childhood psychopathology and adult traits like insomnia and subjective well-being may be primarily driven by genetic factors that influence adult major depression. Additionally, specific childhood phenotypes are genetically associated with educational attainment, BMI and neuroticism.

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.

Transcriptome-wide and stratified genomic structural equation modeling identify neurobiological pathways shared across diverse cognitive traits

Functional genomic methods are needed that consider multiple genetically correlated traits. Here we develop and validate Transcriptome-wide Structural Equation Modeling (T-SEM), a multivariate method for studying the effects of tissue-specific gene expression across genetically overlapping traits. T-SEM allows for modeling effects on broad dimensions spanning constellations of traits, while safeguarding against false positives that can arise when effects of gene expression are specific to a subset of traits. We apply T-SEM to investigate the biological mechanisms shared across seven distinct cognitive traits (N = 11,263-331,679), as indexed by a general dimension of genetic sharing (g). We identify 184 genes whose tissue-specific expression is associated with g, including 10 genes not identified in univariate analysis for the individual cognitive traits for any tissue type, and three genes whose expression explained a significant portion of the genetic sharing across g and different subclusters of psychiatric disorders. We go on to apply Stratified Genomic SEM to identify enrichment for g within 28 functional categories. This includes categories indexing the intersection of protein-truncating variant intolerant (PI) genes and specific neuronal cell types, which we also find to be enriched for the genetic covariance between g and a psychotic disorders factor.

Ultra-rare and common genetic variant analysis converge to implicate negative selection and neuronal processes in the aetiology of schizophrenia

Both common and rare genetic variants (minor allele frequency >1% and <0.1% respectively) have been implicated in the aetiology of schizophrenia. In this study, we integrate single-cell gene expression data with publicly available Genome-Wide Association Study (GWAS) and exome sequenced data in order to investigate in parallel, the enrichment of common and (ultra-)rare variants related to schizophrenia in several functionally relevant gene-sets. Four types of gene-sets were constructed 1) protein-truncating variant (PTV)-intolerant (PI) genes 2) genes expressed in brain cell types and neurons ascertained from mouse and human brain tissue 3) genes defined by synaptic function and location and 4) intersection genes, i.e., PI genes that are expressed in the human and mouse brain cell gene-sets. We show that common as well as ultra-rare schizophrenia-associated variants are overrepresented in PI genes, in excitatory neurons from the prefrontal cortex and hippocampus, medium spiny neurons, and genes enriched for synaptic processes. We also observed stronger enrichment in the intersection genes. Our findings suggest that across the allele frequency spectrum, genes and genetic variants likely to be under stringent selection, and those expressed in particular brain cell types, are involved in the same biological pathways influencing the risk for schizophrenia.

Genome-wide Association Meta-analysis of Childhood and Adolescent Internalizing Symptoms

OBJECTIVE

To investigate the genetic architecture of internalizing symptoms in childhood and adolescence.

METHOD

In 22 cohorts, multiple univariate genome-wide association studies (GWASs) were performed using repeated assessments of internalizing symptoms, in a total of 64,561 children and adolescents between 3 and 18 years of age. Results were aggregated in meta-analyses that accounted for sample overlap, first using all available data, and then using subsets of measurements grouped by rater, age, and instrument.

RESULTS

The meta-analysis of overall internalizing symptoms (INToverall) detected no genome-wide significant hits and showed low single nucleotide polymorphism (SNP) heritability (1.66%, 95% CI = 0.84-2.48%, neffective = 132,260). Stratified analyses indicated rater-based heterogeneity in genetic effects, with self-reported internalizing symptoms showing the highest heritability (5.63%, 95% CI = 3.08%-8.18%). The contribution of additive genetic effects on internalizing symptoms appeared to be stable over age, with overlapping estimates of SNP heritability from early childhood to adolescence. Genetic correlations were observed with adult anxiety, depression, and the well-being spectrum (|rg| > 0.70), as well as with insomnia, loneliness, attention-deficit/hyperactivity disorder, autism, and childhood aggression (range |rg| = 0.42-0.60), whereas there were no robust associations with schizophrenia, bipolar disorder, obsessive-compulsive disorder, or anorexia nervosa.

CONCLUSION

Genetic correlations indicate that childhood and adolescent internalizing symptoms share substantial genetic vulnerabilities with adult internalizing disorders and other childhood psychiatric traits, which could partially explain both the persistence of internalizing symptoms over time and the high comorbidity among childhood psychiatric traits. Reducing phenotypic heterogeneity in childhood samples will be key in paving the way to future GWAS success.

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.

Genetic architecture of 11 major psychiatric disorders at biobehavioral, functional genomic and molecular genetic levels of analysis

We interrogate the joint genetic architecture of 11 major psychiatric disorders at biobehavioral, functional genomic and molecular genetic levels of analysis. We identify four broad factors (neurodevelopmental, compulsive, psychotic and internalizing) that underlie genetic correlations among the disorders and test whether these factors adequately explain their genetic correlations with biobehavioral traits. We introduce stratified genomic structural equation modeling, which we use to identify gene sets that disproportionately contribute to genetic risk sharing. This includes protein-truncating variant-intolerant genes expressed in excitatory and GABAergic brain cells that are enriched for genetic overlap across disorders with psychotic features. Multivariate association analyses detect 152 (20 new) independent loci that act on the individual factors and identify nine loci that act heterogeneously across disorders within a factor. Despite moderate-to-high genetic correlations across all 11 disorders, we find little utility of a single dimension of genetic risk across psychiatric disorders either at the level of biobehavioral correlates or at the level of individual variants.

Within-sibship genome-wide association analyses decrease bias in estimates of direct genetic effects

Estimates from genome-wide association studies (GWAS) of unrelated individuals capture effects of inherited variation (direct effects), demography (population stratification, assortative mating) and relatives (indirect genetic effects). Family-based GWAS designs can control for demographic and indirect genetic effects, but large-scale family datasets have been lacking. We combined data from 178,086 siblings from 19 cohorts to generate population (between-family) and within-sibship (within-family) GWAS estimates for 25 phenotypes. Within-sibship GWAS estimates were smaller than population estimates for height, educational attainment, age at first birth, number of children, cognitive ability, depressive symptoms and smoking. Some differences were observed in downstream SNP heritability, genetic correlations and Mendelian randomization analyses. For example, the within-sibship genetic correlation between educational attainment and body mass index attenuated towards zero. In contrast, analyses of most molecular phenotypes (for example, low-density lipoprotein-cholesterol) were generally consistent. We also found within-sibship evidence of polygenic adaptation on taller height. Here, we illustrate the importance of family-based GWAS data for phenotypes influenced by demographic and indirect genetic effects.

Identifying the Common Genetic Basis of Antidepressant Response

Background

Antidepressants are a first-line treatment for depression. However, only a third of individuals experience remission after the first treatment. Common genetic variation, in part, likely regulates antidepressant response, yet the success of previous genome-wide association studies has been limited by sample size. This study performs the largest genetic analysis of prospectively assessed antidepressant response in major depressive disorder to gain insight into the underlying biology and enable out-of-sample prediction.

Methods

Genome-wide analysis of remission (n remit = 1852, n nonremit = 3299) and percentage improvement (n = 5218) was performed. Single nucleotide polymorphism-based heritability was estimated using genome-wide complex trait analysis. Genetic covariance with eight mental health phenotypes was estimated using polygenic scores/AVENGEME. Out-of-sample prediction of antidepressant response polygenic scores was assessed. Gene-level association analysis was performed using MAGMA and transcriptome-wide association study. Tissue, pathway, and drug binding enrichment were estimated using MAGMA.

Results

Neither genome-wide association study identified genome-wide significant associations. Single nucleotide polymorphism-based heritability was significantly different from zero for remission (h 2 = 0.132, SE = 0.056) but not for percentage improvement (h 2 = -0.018, SE = 0.032). Better antidepressant response was negatively associated with genetic risk for schizophrenia and positively associated with genetic propensity for educational attainment. Leave-one-out validation of antidepressant response polygenic scores demonstrated significant evidence of out-of-sample prediction, though results varied in external cohorts. Gene-based analyses identified ETV4 and DHX8 as significantly associated with antidepressant response.

Conclusions

This study demonstrates that antidepressant response is influenced by common genetic variation, has a genetic overlap schizophrenia and educational attainment, and provides a useful resource for future research. Larger sample sizes are required to attain the potential of genetics for understanding and predicting antidepressant response.

No effects of siblings and twin testosterone transfer on autistic traits

Background

Having twin and non-twin siblings might influence autistic traits both prenatally and postnatally. The twin testosterone transfer hypothesis suggests that girls with a twin brother are exposed to higher levels of prenatal testosterone than girls with a twin sister, and that increased testosterone exposure masculinizes neural development and increases autistic traits. Postnatally, siblings may provide example behaviour, which could reduce autistic traits.

Methods

We studied pre- and postnatal influences of twin and non-twin siblings on mother and teacher-reported autistic traits in 7714 dizygotic twins. We examined the effect of sex of the proband child and of the siblings. We fitted regression models (for boys and girls separately) with sex of co-twin and having older and/or younger siblings of each sex as predictors.

Results

Girls' mother-reported autistic traits were slightly lower for those with a twin brother than those with a twin sister (β = -.08, p = .001, Cohen's d = -.13). This difference was not replicated in teacher-reported autistic traits (β = .01, p = .734). Boys' (mother and teacher-reported) autistic traits were not related to the sex of their co-twin (ps > 0.50). Teacher-reported autistic traits were slightly higher if girls had an older brother (β = .07, p = .013, Cohen's d = .12). Other than this small effect, we found no effect of non-twin siblings on autistic traits in either girls or boys (ps > .18).

Conclusions

We did not find increased autistic traits in girls with a twin brother compared to girls with a twin sister. This finding contributes to a body of literature that rejects the twin testosterone transfer hypothesis. In addition, we found little evidence for pre- and postnatal sibling influences. Our findings align with high heritability and absence of shared-environmental influences in ASD.

Item-Level Genome-Wide Association Study of the Alcohol Use Disorders Identification Test in Three Population-Based Cohorts

OBJECTIVE

Genome-wide association studies (GWASs) of the Alcohol Use Disorders Identification Test (AUDIT), a 10-item screen for alcohol use disorder (AUD), have elucidated novel loci for alcohol consumption and misuse. However, these studies also revealed that GWASs can be influenced by numerous biases (e.g., measurement error, selection bias), which may have led to inconsistent genetic correlations between alcohol involvement and AUD, as well as paradoxically negative genetic correlations between alcohol involvement and psychiatric disorders and/or medical conditions. The authors used genomic structural equation modeling to elucidate the genetics of alcohol consumption and problematic consequences of alcohol use as measured by AUDIT.

METHODS

To explore these unexpected differences in genetic correlations, the authors conducted the first item-level and the largest GWAS of AUDIT items (N=160,824) and applied a multivariate framework to mitigate previous biases.

RESULTS

The authors identified novel patterns of similarity (and dissimilarity) among the AUDIT items and found evidence of a correlated two-factor structure at the genetic level ("consumption" and "problems," rg=0.80). Moreover, by applying empirically derived weights to each of the AUDIT items, the authors constructed an aggregate measure of alcohol consumption that was strongly associated with alcohol dependence (rg=0.67), moderately associated with several other psychiatric disorders, and no longer positively associated with health and positive socioeconomic outcomes. Lastly, by conducting polygenic analyses in three independent cohorts that differed in their ascertainment and prevalence of AUD, the authors identified novel genetic associations between alcohol consumption, alcohol misuse, and health.

CONCLUSIONS

This work further emphasizes the value of AUDIT for both clinical and genetic studies of AUD and the importance of using multivariate methods to study genetic associations that are more closely related to AUD.

[What can genetic research teach us about the overlap and differences between psychiatric disorders?]

Genetic research has made continuous progress in the past years, as a results of the rapidly growing sample sizes and development of novel analytical tools. <br> AIM: To present the latest developments in research into the genetic overlap and differences between psychiatric disorders.<br> METHOD: Description of findings from recent literature.<br> RESULTS: First of all, psychiatric disorders differ in their heritability, and consequently impact of environmental factors. Second, some disorders show high concordance in direction of genetic effects, while other disorders have rather different effects. Using these genetic (dis)similarities, subgroups of psychiatric disorders can also be defined. <br> CONCLUSION: Genetic research is developing at rapid pace, and much more will be learned about the genetic overlap and differences between psychiatric disorders in the next years. The big promise -and challenge- will be to integrate genetic research with research of environmental factors, other biological measures and more detailed phenotypic information.<br>.

Estimating direct and indirect genetic effects on offspring phenotypes using genome-wide summary results data

Estimation of direct and indirect (i.e. parental and/or sibling) genetic effects on phenotypes is becoming increasingly important. We compare several multivariate methods that utilize summary results statistics from genome-wide association studies to determine how well they estimate direct and indirect genetic effects. Using data from the UK Biobank, we contrast point estimates and standard errors at individual loci compared to those obtained using individual level data. We show that Genomic structural equation modelling (SEM) outperforms the other methods in accurately estimating conditional genetic effects and their standard errors. We apply Genomic SEM to fertility data in the UK Biobank and partition the genetic effect into female and male fertility and a sibling specific effect. We identify a novel locus for fertility and genetic correlations between fertility and educational attainment, risk taking behaviour, autism and subjective well-being. We recommend Genomic SEM be used to partition genetic effects into direct and indirect components when using summary results from genome-wide association studies.

Continuity of Genetic Risk for Aggressive Behavior Across the Life-Course

We test whether genetic influences that explain individual differences in aggression in early life also explain individual differences across the life-course. In two cohorts from The Netherlands (N = 13,471) and Australia (N = 5628), polygenic scores (PGSs) were computed based on a genome-wide meta-analysis of childhood/adolescence aggression. In a novel analytic approach, we ran a mixed effects model for each age (Netherlands: 12-70 years, Australia: 16-73 years), with observations at the focus age weighted as 1, and decaying weights for ages further away. We call this approach a 'rolling weights' model. In The Netherlands, the estimated effect of the PGS was relatively similar from age 12 to age 41, and decreased from age 41-70. In Australia, there was a peak in the effect of the PGS around age 40 years. These results are a first indication from a molecular genetics perspective that genetic influences on aggressive behavior that are expressed in childhood continue to play a role later in life.

Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression

Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.

Genetic correlates of socio-economic status influence the pattern of shared heritability across mental health traits

Epidemiological studies show high comorbidity between different mental health problems, indicating that individuals with a diagnosis of one disorder are more likely to develop other mental health problems. Genetic studies reveal substantial sharing of genetic factors across mental health traits. However, mental health is also genetically correlated with socio-economic status (SES), and it is therefore important to investigate and disentangle the genetic relationship between mental health and SES. We used summary statistics from large genome-wide association studies (average N ~ 160,000) to estimate the genetic overlap across nine psychiatric disorders and seven substance use traits and explored the genetic influence of three different indicators of SES. Using genomic structural equation modelling, we show significant changes in patterns of genetic correlations after partialling out SES-associated genetic variation. Our approach allows the separation of disease-specific genetic variation and genetic variation shared with SES, thereby improving our understanding of the genetic architecture of mental health.

Genetic association study of childhood aggression across raters, instruments, and age

Childhood aggressive behavior (AGG) has a substantial heritability of around 50%. Here we present a genome-wide association meta-analysis (GWAMA) of childhood AGG, in which all phenotype measures across childhood ages from multiple assessors were included. We analyzed phenotype assessments for a total of 328 935 observations from 87 485 children aged between 1.5 and 18 years, while accounting for sample overlap. We also meta-analyzed within subsets of the data, i.e., within rater, instrument and age. SNP-heritability for the overall meta-analysis (AGGoverall) was 3.31% (SE = 0.0038). We found no genome-wide significant SNPs for AGGoverall. The gene-based analysis returned three significant genes: ST3GAL3 (P = 1.6E-06), PCDH7 (P = 2.0E-06), and IPO13 (P = 2.5E-06). All three genes have previously been associated with educational traits. Polygenic scores based on our GWAMA significantly predicted aggression in a holdout sample of children (variance explained = 0.44%) and in retrospectively assessed childhood aggression (variance explained = 0.20%). Genetic correlations (rg) among rater-specific assessment of AGG ranged from rg = 0.46 between self- and teacher-assessment to rg = 0.81 between mother- and teacher-assessment. We obtained moderate-to-strong rgs with selected phenotypes from multiple domains, but hardly with any of the classical biomarkers thought to be associated with AGG. Significant genetic correlations were observed with most psychiatric and psychological traits (range [Formula: see text]: 0.19-1.00), except for obsessive-compulsive disorder. Aggression had a negative genetic correlation (rg = ~-0.5) with cognitive traits and age at first birth. Aggression was strongly genetically correlated with smoking phenotypes (range [Formula: see text]: 0.46-0.60). The genetic correlations between aggression and psychiatric disorders were weaker for teacher-reported AGG than for mother- and self-reported AGG. The current GWAMA of childhood aggression provides a powerful tool to interrogate the rater-specific genetic etiology of AGG.

Genetic analyses identify widespread sex-differential participation bias

Genetic association results are often interpreted with the assumption that study participation does not affect downstream analyses. Understanding the genetic basis of participation bias is challenging since it requires the genotypes of unseen individuals. Here we demonstrate that it is possible to estimate comparative biases by performing a genome-wide association study contrasting one subgroup versus another. For example, we showed that sex exhibits artifactual autosomal heritability in the presence of sex-differential participation bias. By performing a genome-wide association study of sex in approximately 3.3 million males and females, we identified over 158 autosomal loci spuriously associated with sex and highlighted complex traits underpinning differences in study participation between the sexes. For example, the body mass index-increasing allele at FTO was observed at higher frequency in males compared to females (odds ratio = 1.02, P = 4.4 × 10^-36). Finally, we demonstrated how these biases can potentially lead to incorrect inferences in downstream analyses and propose a conceptual framework for addressing such biases. Our findings highlight a new challenge that genetic studies may face as sample sizes continue to grow.

Genetic analyses identify widespread sex-differential participation bias

Genetic association results are often interpreted with the assumption that study participation does not affect downstream analyses. Understanding the genetic basis of participation bias is challenging since it requires the genotypes of unseen individuals. Here we demonstrate that it is possible to estimate comparative biases by performing a genome-wide association study contrasting one subgroup versus another. For example, we showed that sex exhibits artifactual autosomal heritability in the presence of sex-differential participation bias. By performing a genome-wide association study of sex in approximately 3.3 million males and females, we identified over 158 autosomal loci spuriously associated with sex and highlighted complex traits underpinning differences in study participation between the sexes. For example, the body mass index-increasing allele at FTO was observed at higher frequency in males compared to females (odds ratio = 1.02, P = 4.4 × 10^-36). Finally, we demonstrated how these biases can potentially lead to incorrect inferences in downstream analyses and propose a conceptual framework for addressing such biases. Our findings highlight a new challenge that genetic studies may face as sample sizes continue to grow.

Response to Comment on "Large-scale GWAS reveals insights into the genetic architecture of same-sex sexual behavior"

Hamer et al argue that the variable "ever versus never had a same-sex partner" does not capture the complexity of human sexuality. We agree and said so in our paper. But Hamer et al neglect to mention that we also reported follow-up analyses showing substantial overlap of the genetic influences on our main variable and on more nuanced measures of sexual behavior, attraction, and identity.

Onset of Preclinical Alzheimer Disease in Monozygotic Twins

Objective

The present work was undertaken to study the genetic contribution to the start of Alzheimer's disease (AD) with amyloid and tau biomarkers in cognitively intact older identical twins.

Methods

We studied in 96 monozygotic twin‐pairs relationships between amyloid‐beta (Aβ) aggregation as measured by the Aβ1–42/1–40 ratio in cerebrospinal fluid (CSF; n = 126) and positron emission tomography (PET, n = 194), and CSF markers for Aβ production (beta‐secretase 1, Aβ1–40, and Aβ1–38) and CSF tau. Associations among markers were tested with generalized estimating equations including a random effect for twin status, adjusted for age, gender, and apolipoprotein E ε4 genotype. We used twin analyses to determine relative contributions of genetic and/or environmental factors to AD pathophysiological processes.

Results

Twenty‐seven individuals (14%) had an abnormal amyloid PET, and 14 twin‐pairs (15%) showed discordant amyloid PET scans. Within twin‐pairs, Aβ production markers and total‐tau (t‐tau) levels strongly correlated (r range = 0.73–0.86, all p < 0.0001), and Aβ aggregation markers and 181‐phosphorylated‐tau (p‐tau) levels correlated moderately strongly (r range = 0.50–0.64, all p < 0.0001). Cross‐twin cross‐trait analysis showed that Aβ1–38 in one twin correlated with Aβ1–42/1–40 ratios, and t‐tau and p‐tau levels in their cotwins (r range = −0.28 to 0.58, all p < .007). Within‐pair differences in Aβ production markers related to differences in tau levels (r range = 0.49–0.61, all p < 0.0001). Twin discordance analyses suggest that Aβ production and tau levels show coordinated increases in very early AD.

Interpretation

Our results suggest a substantial genetic/shared environmental background contributes to both Aβ and tau increases, suggesting that modulation of environmental risk factors may aid in delaying the onset of AD pathophysiological processes. ANN NEUROL 2021;89:987–1000

Investigating the genetic architecture of noncognitive skills using GWAS-by-subtraction

Little is known about the genetic architecture of traits affecting educational attainment other than cognitive ability. We used genomic structural equation modeling and prior genome-wide association studies (GWASs) of educational attainment (n = 1,131,881) and cognitive test performance (n = 257,841) to estimate SNP associations with educational attainment variation that is independent of cognitive ability. We identified 157 genome-wide-significant loci and a polygenic architecture accounting for 57% of genetic variance in educational attainment. Noncognitive genetics were enriched in the same brain tissues and cell types as cognitive performance, but showed different associations with gray-matter brain volumes. Noncognitive genetics were further distinguished by associations with personality traits, less risky behavior and increased risk for certain psychiatric disorders. For socioeconomic success and longevity, noncognitive and cognitive-performance genetics demonstrated associations of similar magnitude. By conducting a GWAS of a phenotype that was not directly measured, we offer a view of genetic architecture of noncognitive skills influencing educational success.

Content, diagnostic, correlational, and genetic similarities between common measures of childhood aggressive behaviors and related psychiatric traits

BACKGROUND

Given the role of childhood aggressive behavior (AGG) in everyday child development, precise and accurate measurement is critical in clinical practice and research. This study aims to quantify agreement among widely used measures of childhood AGG regarding item content, clinical concordance, correlation, and underlying genetic construct.

METHODS

We analyzed data from 1254 Dutch twin pairs (age 8-10 years, 51.1% boys) from a general population sample for whom both parents completed the A-TAC, CBCL, and SDQ at the same occasion.

RESULTS

There was substantial variation in item content among AGG measures, ranging from .00 (i.e., mutually exclusive) to .50 (moderate agreement). Clinical concordance (i.e., do the same children score above a clinical threshold among AGG measures) was very weak to moderate with estimates ranging between .01 and .43 for mother-reports and between .12 and .42 for father-reports. Correlations among scales were weak to strong, ranging from .32 to .70 for mother-reports and from .32 to .64 for father-reports. We found weak to very strong genetic correlations among the measures, with estimates between .65 and .84 for mother-reports and between .30 and .87 for father-reports.

CONCLUSIONS

Our results demonstrated that degree of agreement between measures of AGG depends on the type (i.e., item content, clinical concordance, correlation, genetic correlation) of agreement considered. Because agreement was higher for correlations compared to clinical concordance (i.e., above or below a clinical cutoff), we propose the use of continuous scores to assess AGG, especially for combining data with different measures. Although item content can be different and agreement among observed measures may not be high, the genetic correlations indicate that the underlying genetic liability for childhood AGG is consistent across measures.

Avoiding dynastic, assortative mating, and population stratification biases in Mendelian randomization through within-family analyses

Estimates from Mendelian randomization studies of unrelated individuals can be biased due to uncontrolled confounding from familial effects. Here we describe methods for within-family Mendelian randomization analyses and use simulation studies to show that family-based analyses can reduce such biases. We illustrate empirically how familial effects can affect estimates using data from 61,008 siblings from the Nord-Trøndelag Health Study and UK Biobank and replicated our findings using 222,368 siblings from 23andMe. Both Mendelian randomization estimates using unrelated individuals and within family methods reproduced established effects of lower BMI reducing risk of diabetes and high blood pressure. However, while Mendelian randomization estimates from samples of unrelated individuals suggested that taller height and lower BMI increase educational attainment, these effects were strongly attenuated in within-family Mendelian randomization analyses. Our findings indicate the necessity of controlling for population structure and familial effects in Mendelian randomization studies.

Genetic Associations Between Childhood Psychopathology and Adult Depression and Associated Traits in 42 998 Individuals: A Meta-analysis

Importance

Adult mood disorders are often preceded by behavioral and emotional problems in childhood. It is yet unclear what explains the associations between childhood psychopathology and adult traits.

Objective

To investigate whether genetic risk for adult mood disorders and associated traits is associated with childhood disorders.

Design, Setting, and Participants

This meta-analysis examined data from 7 ongoing longitudinal birth and childhood cohorts from the UK, the Netherlands, Sweden, Norway, and Finland. Starting points of data collection ranged from July 1985 to April 2002. Participants were repeatedly assessed for childhood psychopathology from ages 6 to 17 years. Data analysis occurred from September 2017 to May 2019.

Exposures

Individual polygenic scores (PGS) were constructed in children based on genome-wide association studies of adult major depression, bipolar disorder, subjective well-being, neuroticism, insomnia, educational attainment, and body mass index (BMI).

Main Outcomes and Measures

Regression meta-analyses were used to test associations between PGS and attention-deficit/hyperactivity disorder (ADHD) symptoms and internalizing and social problems measured repeatedly across childhood and adolescence and whether these associations depended on childhood phenotype, age, and rater.

Results

The sample included 42 998 participants aged 6 to 17 years. Male participants varied from 43.0% (1040 of 2417 participants) to 53.1% (2434 of 4583 participants) by age and across all cohorts. The PGS of adult major depression, neuroticism, BMI, and insomnia were positively associated with childhood psychopathology (β estimate range, 0.023-0.042 [95% CI, 0.017-0.049]), while associations with PGS of subjective well-being and educational attainment were negative (β, -0.026 to -0.046 [95% CI, -0.020 to -0.057]). There was no moderation of age, type of childhood phenotype, or rater with the associations. The exceptions were stronger associations between educational attainment PGS and ADHD compared with internalizing problems (Δβ, 0.0561 [Δ95% CI, 0.0318-0.0804]; ΔSE, 0.0124) and social problems (Δβ, 0.0528 [Δ95% CI, 0.0282-0.0775]; ΔSE, 0.0126), and between BMI PGS and ADHD and social problems (Δβ, -0.0001 [Δ95% CI, -0.0102 to 0.0100]; ΔSE, 0.0052), compared with internalizing problems (Δβ, -0.0310 [Δ95% CI, -0.0456 to -0.0164]; ΔSE, 0.0074). Furthermore, the association between educational attainment PGS and ADHD increased with age (Δβ, -0.0032 [Δ 95% CI, -0.0048 to -0.0017]; ΔSE, 0.0008).

Conclusions and Relevance

Results from this study suggest the existence of a set of genetic factors influencing a range of traits across the life span with stable associations present throughout childhood. Knowledge of underlying mechanisms may affect treatment and long-term outcomes of individuals with psychopathology.

Integration of epidemiologic, pharmacologic, genetic and gut microbiome data in a drug-metabolite atlas

Progress in high-throughput metabolic profiling provides unprecedented opportunities to obtain insights into the effects of drugs on human metabolism. The Biobanking BioMolecular Research Infrastructure of the Netherlands has constructed an atlas of drug-metabolite associations for 87 commonly prescribed drugs and 150 clinically relevant plasma-based metabolites assessed by proton nuclear magnetic resonance. The atlas includes a meta-analysis of ten cohorts (18,873 persons) and uncovers 1,071 drug-metabolite associations after evaluation of confounders including co-treatment. We show that the effect estimates of statins on metabolites from the cross-sectional study are comparable to those from intervention and genetic observational studies. Further data integration links proton pump inhibitors to circulating metabolites, liver function, hepatic steatosis and the gut microbiome. Our atlas provides a tool for targeted experimental pharmaceutical research and clinical trials to improve drug efficacy, safety and repurposing. We provide a web-based resource for visualization of the atlas (http://bbmri.researchlumc.nl/atlas/).

Phenome-wide investigation of health outcomes associated with genetic predisposition to loneliness

Humans are social animals that experience intense suffering when they perceive a lack of social connection. Modern societies are experiencing an epidemic of loneliness. Although the experience of loneliness is universally human, some people report experiencing greater loneliness than others. Loneliness is more strongly associated with mortality than obesity, emphasizing the need to understand the nature of the relationship between loneliness and health. Although it is intuitive that circumstantial factors such as marital status and age influence loneliness, there is also compelling evidence of a genetic predisposition toward loneliness. To better understand the genetic architecture of loneliness and its relationship with associated outcomes, we extended the genome-wide association study meta-analysis of loneliness to 511 280 subjects, and detect 19 significant genetic variants from 16 loci, including four novel loci, as well as 58 significantly associated genes. We investigated the genetic overlap with a wide range of physical and mental health traits by computing genetic correlations and by building loneliness polygenic scores in an independent sample of 18 498 individuals with EHR data to conduct a PheWAS with. A genetic predisposition toward loneliness was associated with cardiovascular, psychiatric, and metabolic disorders and triglycerides and high-density lipoproteins. Mendelian randomization analyses showed evidence of a causal, increasing, the effect of both BMI and body fat on loneliness. Our results provide a framework for future studies of the genetic basis of loneliness and its relationship to mental and physical health.

Genetic correlates of social stratification in Great Britain

Human DNA polymorphisms vary across geographic regions, with the most commonly observed variation reflecting distant ancestry differences. Here we investigate the geographic clustering of common genetic variants that influence complex traits in a sample of ~450,000 individuals from Great Britain. Of 33 traits analysed, 21 showed significant geographic clustering at the genetic level after controlling for ancestry, probably reflecting migration driven by socioeconomic status (SES). Alleles associated with educational attainment (EA) showed the most clustering, with EA-decreasing alleles clustering in lower SES areas such as coal mining areas. Individuals who leave coal mining areas carry more EA-increasing alleles on average than those in the rest of Great Britain. The level of geographic clustering is correlated with genetic associations between complex traits and regional measures of SES, health and cultural outcomes. Our results are consistent with the hypothesis that social stratification leaves visible marks in geographic arrangements of common allele frequencies and gene-environment correlations.