The Minnesota Center for Twin and Family Research genome-wide association study

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity

STUDY QUESTION

Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins?

SUMMARY ANSWER

We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3.

WHAT IS KNOWN ALREADY

The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures.

STUDY DESIGN, SIZE, DURATION

We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation.

MAIN RESULTS AND THE ROLE OF CHANCE

This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle.

LARGE SCALE DATA

The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/).

LIMITATIONS, REASONS FOR CAUTION

Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility.

WIDER IMPLICATIONS OF THE FINDINGS

About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility.

STUDY FUNDING/COMPETING INTEREST(S)

Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Skłodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests.

TRIAL REGISTRATION NUMBER

N/A.

Evaluating substance use outcomes of recreational cannabis legalization using a unique co-twin control design

Background: As more states pass recreational cannabis legalization (RCL), we must understand how RCL affects substance use.Objectives: The current study aims to examine the effect of RCL on lifetime and past-year use of cannabis, alcohol, tobacco, and other drugs, frequency of cannabis, alcohol, and tobacco use, co-use of cannabis with alcohol and tobacco, and consequences from cannabis and alcohol use.Methods: We used a unique, co-twin control design of twin pairs who were discordant for living in a state with RCL between 2018 and 2021. The sample consisted of 3,830 adult twins (41% male), including 232 twin pairs discordant for RCL. Problems from alcohol and cannabis use were assessed via the Brief Marijuana Consequences Questionnaire and the Brief Young Adult Alcohol Consequences Questionnaire.Results: Results indicated that the twin living in an RCL state was more likely to endorse past-year cannabis use (OR = 1.56, p = .009), greater number of cannabis use days in the past 6 months (β = 0.47, p = .019), but not more negative consequences from cannabis use (β = 0.21, p = .456) compared to their co-twin in a non-RCL state. There were no differences within-twin pairs in frequency of alcohol use (β=-0.05, p = .601), but the RCL twin reported fewer negative consequences from alcohol use (β=-0.29, p = .016) compared to their co-twin in a non-RCL state. We did not observe any other differences within-twin pairs on other outcomes.Conclusion: These results suggest that living in an RCL state is associated with greater cannabis frequency but not more negative consequences from cannabis use than living in a non-RCL state.

Exploring Relationships Between Internalizing Problems and Risky Sexual Behavior: A Twin Study

Previous research links risky sexual behavior (RSB) to externalizing problems and to substance use, but little research has been conducted on relationships between internalizing problems (INT) and RSB. The current study addresses that literature gap, using both a twin sample from Colorado (N = 2567) and a second twin sample from Minnesota (N = 1131) in attempt to replicate initial results. We explored the hypothesis that the latent variable INT would be more strongly associated with the latent variable RSB for females than for males, examining relationships between INT and RSB via phenotypic confirmatory factor analysis and multivariate twin analyses. We found a small but significant phenotypic association between the latent variables. However, despite using two large twin samples, limited power restricted our ability to identify the genetic and environmental mechanisms underlying this association. Our sex differences hypothesis was not fully supported in either sample and requires further investigation. Our findings illustrate the complexity of the relationship between internalizing problems and risky sexual behavior.

Simultaneous selection of multiple important single nucleotide polymorphisms in familial genome wide association studies data

We propose a resampling-based fast variable selection technique for detecting relevant single nucleotide polymorphisms (SNP) in a multi-marker mixed effect model. Due to computational complexity, current practice primarily involves testing the effect of one SNP at a time, commonly termed as 'single SNP association analysis'. Joint modeling of genetic variants within a gene or pathway may have better power to detect associated genetic variants, especially the ones with weak effects. In this paper, we propose a computationally efficient model selection approach-based on the e-values framework-for single SNP detection in families while utilizing information on multiple SNPs simultaneously. To overcome computational bottleneck of traditional model selection methods, our method trains one single model, and utilizes a fast and scalable bootstrap procedure. We illustrate through numerical studies that our proposed method is more effective in detecting SNPs associated with a trait than either single-marker analysis using family data or model selection methods that ignore the familial dependency structure. Further, we perform gene-level analysis in Minnesota Center for Twin and Family Research (MCTFR) dataset using our method to detect several SNPs using this that have been implicated to be associated with alcohol consumption.

An Efficient Bayesian Method for Estimating the Degree of the Skewness of X Chromosome Inactivation Based on the Mixture of General Pedigrees and Unrelated Females

Skewed X chromosome inactivation (XCI-S) has been reported to be associated with some X-linked diseases. Several methods have been proposed to estimate the degree of XCI-S (denoted as γ) for quantitative and qualitative traits based on unrelated females. However, there is no method available for estimating γ based on general pedigrees. Therefore, in this paper, we propose a Bayesian method to obtain the point estimate and the credible interval of γ based on the mixture of general pedigrees and unrelated females (called mixed data for brevity), which is also suitable for only general pedigrees. We consider the truncated normal prior and the uniform prior for γ. Further, we apply the eigenvalue decomposition and Cholesky decomposition to our proposed methods to accelerate the computation speed. We conduct extensive simulation studies to compare the performances of our proposed methods and two existing Bayesian methods which are only applicable to unrelated females. The simulation results show that the incorporation of general pedigrees can improve the efficiency of the point estimation and the precision and the accuracy of the interval estimation of γ. Finally, we apply the proposed methods to the Minnesota Center for Twin and Family Research data for their practical use.

A Century of Behavioral Genetics at the University of Minnesota

The University of Minnesota has played an important role in the resurgence and eventual mainstreaming of human behavioral genetics in psychology and psychiatry. We describe this history in the context of three major movements in behavioral genetics: (1) radical eugenics in the early 20th century, (2) resurgence of human behavioral genetics in the 1960s, largely using twin and adoption designs to obtain more precise estimates of genetic and environmental influences on individual differences in behavior; and (3) use of measured genotypes to understand behavior. University of Minnesota scientists made significant contributions especially in (2) and (3) in the domains of cognitive ability, drug abuse and mental health, and endophenotypes. These contributions are illustrated through a historical perspective of major figures and events in behavioral genetics.

Uncovering the genetic architecture of broad antisocial behavior through a genome-wide association study meta-analysis

Despite the substantial heritability of antisocial behavior (ASB), specific genetic variants robustly associated with the trait have not been identified. The present study by the Broad Antisocial Behavior Consortium (BroadABC) meta-analyzed data from 28 discovery samples (N = 85,359) and five independent replication samples (N = 8058) with genotypic data and broad measures of ASB. We identified the first significant genetic associations with broad ASB, involving common intronic variants in the forkhead box protein P2 (FOXP2) gene (lead SNP rs12536335, p = 6.32 × 10-10). Furthermore, we observed intronic variation in Foxp2 and one of its targets (Cntnap2) distinguishing a mouse model of pathological aggression (BALB/cJ strain) from controls (BALB/cByJ strain). Polygenic risk score (PRS) analyses in independent samples revealed that the genetic risk for ASB was associated with several antisocial outcomes across the lifespan, including diagnosis of conduct disorder, official criminal convictions, and trajectories of antisocial development. We found substantial genetic correlations of ASB with mental health (depression rg = 0.63, insomnia rg = 0.47), physical health (overweight rg = 0.19, waist-to-hip ratio rg = 0.32), smoking (rg = 0.54), cognitive ability (intelligence rg = -0.40), educational attainment (years of schooling rg = -0.46) and reproductive traits (age at first birth rg = -0.58, father's age at death rg = -0.54). Our findings provide a starting point toward identifying critical biosocial risk mechanisms for the development of ASB.

Longitudinal effects and environmental moderation of ALDH2 and ADH1B gene variants on substance use from age 14 to 40

Alcohol use and dependence are strongly affected by variation in aldehyde dehydrogenase (ALDH2) and, to a lesser extent, alcohol dehydrogenase (ADH1B) genes. We use this genetic variation with an adoption design to test the causal role of alcohol use on other drug use, as well as the moderating role of adoptive parent, sibling, and peer alcohol use. Longitudinal models were run on 412 genotyped adopted individuals of East Asian ancestry with multiple assessments between ages 14 and 40. We found robust associations between alcohol frequency, quantity, and maximum drinks and ALDH2, but not ADH1B, status. The magnitude of the ALDH2 protective effect increased with age, particularly for maximum drinks, though estimates were smaller than previously reported in ancestrally similar individuals in East/North-East Asian countries. These results suggest that sociocultural factors in Minnesota may reduce the protective effects of ALDH2. We found that peer alcohol use, but not parent or sibling use, predicted adopted offspring's use, and that these environmental influences did not vary by ALDH2 status. Finally, we did not find strong evidence of associations between ALDH2 status and tobacco, marijuana, or illegal drug use, contrary to expectation if alcohol serves as a gateway to use of other drugs.

Rare variant association study of veteran twin whole-genomes links severe depression with a nonsynonymous change in the neuronal gene BHLHE22

OBJECTIVES

Major Depressive Disorder (MDD) is a complex neuropsychiatric disease with known genetic associations, but without known links to rare variation in the human genome. Here we aim to identify rare genetic variants associated with MDD using deep whole-genome sequencing data in an independent population.

METHODS

We report the sequencing of 1,688 whole genomes in a large sample of male-male Veteran twins. Depression status was classified based on a structured diagnostic interview according to DSM-III-R diagnostic criteria. Searching only rare variants in genomic regions from recent GWAS on MDD, we used the optimised sequence kernel association test and Fisher's Exact test to fine map loci associated with severe depression.

RESULTS

Our analysis identified one gene associated with severe depression, basic helix loop helix e22 (P_Adjusted = 0.03) via SKAT-O test between unrelated severely depressed cases compared to unrelated non-depressed controls. The same gene BHLHE22 had a non-silent variant rs13279074 (P_Adjusted = 0.032) based on a single variant Fisher's Exact test between unrelated severely depressed cases compared to unrelated non-depressed controls.

CONCLUSION

The gene BHLHE22 shows compelling genetic evidence of directly impacting the severe depression phenotype. Together these results advance understanding of the genetic contribution to major depressive disorder in a new cohort and link a rare variant to severe forms of the disorder.

Genetic and environmental contributions to IQ in adoptive and biological families with 30-year-old offspring

While adoption studies have provided key insights into the influence of the familial environment on IQ scores of adolescents and children, few have followed adopted offspring long past the time spent living in the family home. To improve confidence about the extent to which shared environment exerts enduring effects on IQ, we estimated genetic and environmental effects on adulthood IQ in a unique sample of 486 biological and adoptive families. These families, tested previously on measures of IQ when offspring averaged age 15, were assessed a second time nearly two decades later ( M offspring age = 32 years). We estimated the proportions of the variance in IQ attributable to environmentally mediated effects of parental IQs, sibling-specific shared environment, and gene-environment covariance to be .01 [95% CI .00, .02], .04 [95% CI .00, .15], and .03 [95% CI .00, .07] respectively; these components jointly accounted for 8 percent of the IQ variance in adulthood. The heritability was estimated to be .42 [95% CI .21, .64]. Together, these findings provide further evidence for the predominance of genetic influences on adult intelligence over any other systematic source of variation.

The role of parental genotype in predicting offspring years of education: evidence for genetic nurture

Similarities between parent and offspring are widespread in psychology; however, shared genetic variants often confound causal inference for offspring outcomes. A polygenic score (PGS) derived from genome-wide association studies (GWAS) can be used to test for the presence of parental influence that controls for genetic variants shared across generations. We use a PGS for educational attainment (EA3; N ≈ 750 thousand) to predict offspring years of education in a sample of 2517 twins and both parents. We find that within families, the dizygotic twin with the higher PGS is more likely to attain higher education (unstandardized β = 0.32; p < 0.001). Additionally, however, we find an effect of parental genotype on offspring outcome that is independent of the offspring's own genotype; this raises the variance explained in offspring years of education from 9.3 to 11.1% (∆R² = 0.018, p < 0.001). Controlling for parental IQ or socioeconomic status substantially attenuated or eliminated this effect of parental genotype. These findings suggest a role of environmental factors affected by heritable characteristics of the parents in fostering offspring years of education.

A statistical measure for the skewness of X chromosome inactivation for quantitative traits and its application to the MCTFR data

Background

X chromosome inactivation (XCI) is that one of two chromosomes in mammalian females is silenced during early development of embryos. There has been a statistical measure for the degree of the skewness of XCI for qualitative traits. However, no method is available for such task at quantitative trait loci.

Results

In this article, we extend the existing statistical measure for the skewness of XCI for qualitative traits, and the likelihood ratio, Fieller’s and delta methods for constructing the corresponding confidence intervals, and make them accommodate quantitative traits. The proposed measure is a ratio of two linear regression coefficients when association exists. Noting that XCI may cause variance heterogeneity of the traits across different genotypes in females, we obtain the point estimate and confidence intervals of the measure by incorporating such information. The hypothesis testing of the proposed methods is also investigated. We conduct extensive simulation studies to assess the performance of the proposed methods. Simulation results demonstrate that the median of the point estimates of the measure is very close to the pre-specified true value. The likelihood ratio and Fieller’s methods control the size well, and have the similar test power and accurate coverage probability, which perform better than the delta method. So far, we are not aware of any association study for the X-chromosomal loci in the Minnesota Center for Twin and Family Research data. So, we apply our proposed methods to these data for their practical use and find that only the rs792959 locus, which is simultaneously associated with the illicit drug composite score and behavioral disinhibition composite score, may undergo XCI skewing. However, this needs to be confirmed by molecular genetics.

Conclusions

We recommend the Fieller’s method in practical use because it is a non-iterative procedure and has the similar performance to the likelihood ratio method.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00978-z.

Genome-Wide SNP Analysis for Milk Performance Traits in Indigenous Sheep: A Case Study in the Egyptian Barki Sheep

Simple Summary

The Barki sheep is one of the three main breeds in Egypt, which is spread mainly throughout the northwestern coastal zone, which has harsh conditions. Considering the harsh, semi-arid habitat of this breed, milk performance traits such as milk yield and milk composition have a very important role in the feeding of newborn lambs and affect their growth during the early stage of life. In this study, rare milk performance data and genomic information of Barki sheep were used to uncover diversified genomic regions that could explain the variability of milk yield and milk quality traits in the studied population of Barki ewes. Genome-wide analysis identified genomic regions harboring interesting candidate genes such as SLC5A8, NUB1, TBC1D1, KLF3 and ABHD5 for milk yield and PPARA and FBLN1 genes for milk quality traits. The findings offer valuable information for obtaining a better understanding of the genetics of milk performance traits and contribute to the genetic improvement of these traits in Barki sheep.

Abstract

Sheep milk yield and milk composition traits play an important role in supplying newborn lambs with essential components such as amino acids, energy, vitamins and immune antibodies and are also of interest in terms of the nutritional value of the milk for human consumption. The aim of this study was to identify genomic regions and candidate genes for milk yield and milk composition traits through genome-wide SNP analyses between high and low performing ewes of the Egyptian Barki sheep breed, which is well adapted to the harsh conditions of North-East Africa. Therefore, out of a herd of 111 ewes of the Egyptian Barki sheep breed (IBD = 0.08), ewes representing extremes in milk yield and milk quality traits (n = 25 for each group of animals) were genotyped using the Illumina OvineSNP50 V2 BeadChip. The fixation index (F_ST) for each SNP was calculated between the diversified groups. F_ST values were Z-transformed and used to identify putative SNPs for further analysis (Z(F_ST) > 10). Genome-wide SNP analysis revealed genomic regions covering promising candidate genes related to milk performance traits such as SLC5A8, NUB1, TBC1D1, KLF3 and ABHD5 for milk yield and PPARA and FBLN1 genes for milk quality trait. The results of this study may contribute to the genetic improvement of milk performance traits in Barki sheep breed and to the general understanding of the genetic contribution to variability in milk yield and quality traits.

Mechanisms of parent-child transmission of tobacco and alcohol use with polygenic risk scores: Evidence for a genetic nurture effect

Parent-child similarity is a function of genetic and environmental transmission. In addition, genetic effects not transmitted to offspring may drive parental behavior, thereby affecting the rearing environment of the child. Measuring genetic proclivity directly, through polygenic risk scores (PRSs), provides a way to test for the effect of nontransmitted parental genotype, on offspring outcome, termed a genetic nurture effect-in other words, if and how parental genomes might affect their children through the environment. The current study used polygenic risk scores for smoking initiation, cigarettes per day, and drinks per week to predict substance use in a sample of 3,008 twins, assessed prospectively from age 17-29, and their parents, from the Minnesota Center for Twin and Family Research. Mixed-effects models were used to test for a genetic nurture effect whereby parental PRSs predict offspring tobacco and alcohol use after statistically adjusting for offspring's own PRS. Parental smoking initiation PRS predicted offspring cigarettes per day at age 24 (β = .103, 95% CI [.03, .17]) and alcohol use at age 17 (β = .091, 95% CI [.04, .14]) independent of shared genetics. There was also a suggestive independent association between the parent PRS and offspring smoking at age 17 (β = .096; 95% CI [.02, .17]). Mediation analyses provided some evidence for environmental effects of parental smoking, alcohol use, and family socioeconomic status. These findings, and more broadly the molecular genetic method used, have implications on the identification of environmental effects on developmental outcomes such as substance use. (PsycInfo Database Record (c) 2021 APA, all rights reserved).

Using multivariate endophenotypes to identify psychophysiological mechanisms associated with polygenic scores for substance use, schizophrenia, and education attainment

BACKGROUND

To better characterize brain-based mechanisms of polygenic liability for psychopathology and psychological traits, we extended our previous report (Liu et al. Psychophysiological endophenotypes to characterize mechanisms of known schizophrenia genetic loci. Psychological Medicine, 2017), focused solely on schizophrenia, to test the association between multivariate psychophysiological candidate endophenotypes (including novel measures of θ/δ oscillatory activity) and a range of polygenic scores (PGSs), namely alcohol/cannabis/nicotine use, an updated schizophrenia PGS (containing 52 more genome-wide significant loci than the PGS used in our previous report) and educational attainment.

METHOD

A large community-based twin/family sample (N = 4893) was genome-wide genotyped and imputed. PGSs were constructed for alcohol use, regular smoking initiation, lifetime cannabis use, schizophrenia, and educational attainment. Eleven endophenotypes were assessed: visual oddball task event-related electroencephalogram (EEG) measures (target-related parietal P3 amplitude, frontal θ, and parietal δ energy/inter-trial phase clustering), band-limited resting-state EEG power, antisaccade error rate. Principal component analysis exploited covariation among endophenotypes to extract a smaller number of meaningful dimensions/components for statistical analysis.

RESULTS

Endophenotypes were heritable. PGSs showed expected intercorrelations (e.g. schizophrenia PGS correlated positively with alcohol/nicotine/cannabis PGSs). Schizophrenia PGS was negatively associated with an event-related P3/δ component [β = -0.032, nonparametric bootstrap 95% confidence interval (CI) -0.059 to -0.003]. A prefrontal control component (event-related θ/antisaccade errors) was negatively associated with alcohol (β = -0.034, 95% CI -0.063 to -0.006) and regular smoking PGSs (β = -0.032, 95% CI -0.061 to -0.005) and positively associated with educational attainment PGS (β = 0.031, 95% CI 0.003-0.058).

CONCLUSIONS

Evidence suggests that multivariate endophenotypes of decision-making (P3/δ) and cognitive/attentional control (θ/antisaccade error) relate to alcohol/nicotine, schizophrenia, and educational attainment PGSs and represent promising targets for future research.

A robust and unified framework for estimating heritability in twin studies using generalized estimating equations

The 'heritability' of a phenotype measures the proportion of trait variance due to genetic factors in a population. In the past 50 years, studies with monozygotic and dizygotic twins have estimated heritability for 17,804 traits;¹ thus twin studies are popular for estimating heritability. Researchers are often interested in estimating heritability for non-normally distributed outcomes such as binary, counts, skewed or heavy-tailed continuous traits. In these settings, the traditional normal ACE model (NACE) and Falconer's method can produce poor coverage of the true heritability. Therefore, we propose a robust generalized estimating equations (GEE2) framework for estimating the heritability of non-normally distributed outcomes. The traditional NACE and Falconer's method are derived within this unified GEE2 framework, which additionally provides robust standard errors. Although the traditional Falconer's method cannot adjust for covariates, the corresponding 'GEE2-Falconer' can incorporate mean and variance-level covariate effects (e.g. let heritability vary by sex or age). Given a non-normally distributed outcome, the GEE2 models are shown to attain better coverage of the true heritability compared to traditional methods. Finally, a scenario is demonstrated where NACE produces biased estimates of heritability while Falconer remains unbiased. Therefore, we recommend GEE2-Falconer for estimating the heritability of non-normally distributed outcomes in twin studies.

Modeling the Dependence Structure in Genome Wide Association Studies of Binary Phenotypes in Family Data

Genome-wide association studies (GWASs) are a popular tool for detecting association between genetic variants or single nucleotide polymorphisms (SNPs) and complex traits. Family data introduce complexity due to the non-independence of the family members. Methods for non-independent data are well established, but when the GWAS contains distinct family types, explicit modeling of between-family-type differences in the dependence structure comes at the cost of significantly increased computational burden. The situation is exacerbated with binary traits. In this paper, we perform several simulation studies to compare multiple candidate methods to perform single SNP association analysis with binary traits. We consider generalized estimating equations (GEE), generalized linear mixed models (GLMMs), or generalized least square (GLS) approaches. We study the influence of different working correlation structures for GEE on the GWAS findings and also the performance of different analysis method(s) to conduct a GWAS with binary trait data in families. We discuss the merits of each approach with attention to their applicability in a GWAS. We also compare the performances of the methods on the alcoholism data from the Minnesota Center for Twin and Family Research (MCTFR) study.

The Contribution of Cognitive and Noncognitive Skills to Intergenerational Social Mobility

We investigated intergenerational educational and occupational mobility in a sample of 2,594 adult offspring and 2,530 of their parents. Participants completed assessments of general cognitive ability and five noncognitive factors related to social achievement; 88% were also genotyped, allowing computation of educational-attainment polygenic scores. Most offspring were socially mobile. Offspring who scored at least 1 standard deviation higher than their parents on both cognitive and noncognitive measures rarely moved down and frequently moved up. Polygenic scores were also associated with social mobility. Inheritance of a favorable subset of parent alleles was associated with moving up, and inheritance of an unfavorable subset was associated with moving down. Parents' education did not moderate the association of offspring's skill with mobility, suggesting that low-skilled offspring from advantaged homes were not protected from downward mobility. These data suggest that cognitive and noncognitive skills as well as genetic factors contribute to the reordering of social standing that takes place across generations.

Elucidating the role of maternal environmental exposures on offspring health and disease using two-sample Mendelian randomization

BACKGROUND

There is considerable interest in estimating the causal effect of a range of maternal environmental exposures on offspring health-related outcomes. Previous attempts to do this using Mendelian randomization methodologies have been hampered by the paucity of epidemiological cohorts with large numbers of genotyped mother-offspring pairs.

METHODS

We describe a new statistical model that we have created which can be used to estimate the effect of maternal genotypes on offspring outcomes conditional on offspring genotype, using both individual-level and summary-results data, even when the extent of sample overlap is unknown.

RESULTS

We describe how the estimates obtained from our method can subsequently be used in large-scale two-sample Mendelian randomization studies to investigate the causal effect of maternal environmental exposures on offspring outcomes. This includes studies that aim to assess the causal effect of in utero exposures related to fetal growth restriction on future risk of disease in offspring. We illustrate our framework using examples related to offspring birthweight and cardiometabolic disease, although the general principles we espouse are relevant for many other offspring phenotypes.

CONCLUSIONS

We advocate for the establishment of large-scale international genetics consortia that are focused on the identification of maternal genetic effects and committed to the public sharing of genome-wide summary-results data from such efforts. This information will facilitate the application of powerful two-sample Mendelian randomization studies of maternal exposures and offspring outcomes.

Adolescent Externalizing Psychopathology and Its Prospective Relationship to Marijuana Use Development from Age 14 to 30: Replication Across Independent Longitudinal Twin Samples

Externalizing psychopathology in early adolescence is a highly heritable risk factor for drug use, yet how it relates to marijuana use development is not well-characterized. We evaluate this issue in independent twin samples from Colorado (N = 2608) and Minnesota (N = 3630), assessed from adolescence to early adulthood. We used a biometric latent growth model of marijuana use frequency with data from up to five waves of assessment from ages 14 to 30, to examine change in marijuana use and its relationship with a factor model of adolescent externalizing psychopathology. The factor structure of adolescent externalizing psychopathology was similar across samples, as was the association between that common factor and early marijuana use (Minnesota r = 0.67 [0.60, 0.75]; Colorado r = 0.69 [0.59, 0.78]), and increase in use (Minnesota r = 0.18 [0.10, 0.26]; Colorado r = 0.20 [0.07, 0.34]). Early use was moderately heritable in both samples (Minnesota h² = 0.57 [0.37, 0.79]; Colorado h² = 0.42 [0.14, 0.73]). Increase in use was highly heritable in Minnesota (h² = 0.82 [0.72, 0.88]), less so in Colorado (h² = 0.22 [0.01, 0.66]), and shared environmental effects were larger in Colorado (c² = 0.55 [0.14, 0.83]) than Minnesota (c² = 0 [0, 0.06]). We found moderate genetic correlations between externalizing psychopathology and early use in both samples. Finally, additional analyses in the Minnesota sample indicated that marijuana use decreased during the late 20s. This decline is strongly heritable (h² = 0.73 [0.49, 0.91]) and moderately negatively correlated with adolescent externalizing psychopathology (r = - 0.41 [- 0.54, - 0.28]). Adolescent externalizing psychopathology is genetically correlated with change in late adolescent marijuana use (late teens, early 20s), as well as maintenance of use in early adulthood (late 20 s) even after controlling for the effects of early use.

Minnesota Center for Twin and Family Research

The Minnesota Center for Twin and Family Research (MCTFR) comprises multiple longitudinal, community-representative investigations of twin and adoptive families that focus on psychological adjustment, personality, cognitive ability and brain function, with a special emphasis on substance use and related psychopathology. The MCTFR includes the Minnesota Twin Registry (MTR), a cohort of twins who have completed assessments in middle and older adulthood; the Minnesota Twin Family Study (MTFS) of twins assessed from childhood and adolescence into middle adulthood; the Enrichment Study (ES) of twins oversampled for high risk for substance-use disorders assessed from childhood into young adulthood; the Adolescent Brain (AdBrain) study, a neuroimaging study of adolescent twins; and the Siblings Interaction and Behavior Study (SIBS), a study of adoptive and nonadoptive families assessed from adolescence into young adulthood. Here we provide a brief overview of key features of these established studies and describe new MCTFR investigations that follow up and expand upon existing studies or recruit and assess new samples, including the MTR Study of Relationships, Personality, and Health (MTR-RPH); the Colorado-Minnesota (COMN) Marijuana Study; the Adolescent Brain Cognitive Development (ABCD) study; the Colorado Online Twins (CoTwins) study and the Children of Twins (CoT) study.

The causal influence of brain size on human intelligence: Evidence from within-family phenotypic associations and GWAS modeling

There exists a moderate correlation between MRI-measured brain size and the general factor of IQ performance (g), but the question of whether the association reflects a theoretically important causal relationship or spurious confounding remains somewhat open. Previous small studies (n < 100) looking for the persistence of this correlation within families failed to find a tendency for the sibling with the larger brain to obtain a higher test score. We studied the within-family relationship between brain volume and intelligence in the much larger sample provided by the Human Connectome Project (n = 1,022) and found a highly significant correlation (disattenuated ρ = 0.18, p < .001). We replicated this result in the Minnesota Center for Twin and Family Research (n = 2,698), finding a highly significant within-family correlation between head circumference and intelligence (disattenuated ρ = 0.19, p < .001). We also employed novel methods of causal inference relying on summary statistics from genome-wide association studies (GWAS) of head size (n ≈ 10,000) and measures of cognition (257,000 < n < 767,000). Using bivariate LD Score regression, we found a genetic correlation between intracranial volume (ICV) and years of education (EduYears) of 0.41 (p < .001). Using the Latent Causal Variable method, we found a genetic causality proportion of 0.72 (p < .001); thus the genetic correlation arises from an asymmetric pattern, extending to sub-significant loci, of genetic variants associated with ICV also being associated with EduYears but many genetic variants associated with EduYears not being associated with ICV. This is the pattern of genetic results expected from a causal effect of brain size on intelligence. These findings give reason to take up the hypothesis that the dramatic increase in brain volume over the course of human evolution has been the result of natural selection favoring general intelligence.

The accuracy of LD Score regression as an estimator of confounding and genetic correlations in genome-wide association studies

To infer that a single-nucleotide polymorphism (SNP) either affects a phenotype or is linkage disequilibrium with a causal site, we must have some assurance that any SNP-phenotype correlation is not the result of confounding with environmental variables that also affect the trait. In this study, we study the properties of linkage disequilibrium (LD) Score regression, a recently developed method for using summary statistics from genome-wide association studies to ensure that confounding does not inflate the number of false positives. We do not treat the effects of genetic variation as a random variable and thus are able to obtain results about the unbiasedness of this method. We demonstrate that LD Score regression can produce estimates of confounding at null SNPs that are unbiased or conservative under fairly general conditions. This robustness holds in the case of the parent genotype affecting the offspring phenotype through some environmental mechanism, despite the resulting correlation over SNPs between LD Scores and the degree of confounding. Additionally, we demonstrate that LD Score regression can produce reasonably robust estimates of the genetic correlation, even when its estimates of the genetic covariance and the two univariate heritabilities are substantially biased.

Extremely rare variants reveal patterns of germline mutation rate heterogeneity in humans

A detailed understanding of the genome-wide variability of single-nucleotide germline mutation rates is essential to studying human genome evolution. Here, we use ~36 million singleton variants from 3560 whole-genome sequences to infer fine-scale patterns of mutation rate heterogeneity. Mutability is jointly affected by adjacent nucleotide context and diverse genomic features of the surrounding region, including histone modifications, replication timing, and recombination rate, sometimes suggesting specific mutagenic mechanisms. Remarkably, GC content, DNase hypersensitivity, CpG islands, and H3K36 trimethylation are associated with both increased and decreased mutation rates depending on nucleotide context. We validate these estimated effects in an independent dataset of ~46,000 de novo mutations, and confirm our estimates are more accurate than previously published results based on ancestrally older variants without considering genomic features. Our results thus provide the most refined portrait to date of the factors contributing to genome-wide variability of the human germline mutation rate.

A linear mixed model framework for gene-based gene-environment interaction tests in twin studies

Interaction between genes and environments (G×E) can be well investigated in families due to the shared genes and environment among family members. However, the majority of the current tests of G×E interaction between a set of variants and an environment are only suitable for studies with unrelated subjects. In this paper, we extend several G×E interaction tests to a linear mixed model framework to study interaction between a set of correlated environments and a candidate gene in families. The correlated environments can either be modeled separately or jointly in one model. We demonstrate theoretically that the tests developed by modeling correlated environments separately are valid and present a computationally fast alternative to detect G×E interaction in families. For either strategy, we propose treating the genetic main effects as a random effect to reduce the number of main-effect parameters and thus improve the power to detect interactions. Additionally, we propose a generalization of a test of interaction that adaptively sums the interactions using a sequential algorithm. This generalized set of tests, referred to as the sequential algorithm for the sum of powered score (Seq-SPU) family of tests, can be expressed as a weighted version of the SPU. We find that the adaptive version of our test, Seq-aSPU, can outperform aSPU in cases where the interactions effects are in opposite directions. We applied these methods to the Minnesota Center for Twin and Family Research data set and found one significant gene in interaction with four psychosocial environmental factors affecting the alcohol consumption among the twins.

Proper conditional analysis in the presence of missing data: Application to large scale meta-analysis of tobacco use phenotypes

Meta-analysis of genetic association studies increases sample size and the power for mapping complex traits. Existing methods are mostly developed for datasets without missing values, i.e. the summary association statistics are measured for all variants in contributing studies. In practice, genotype imputation is not always effective. This may be the case when targeted genotyping/sequencing assays are used or when the un-typed genetic variant is rare. Therefore, contributed summary statistics often contain missing values. Existing methods for imputing missing summary association statistics and using imputed values in meta-analysis, approximate conditional analysis, or simple strategies such as complete case analysis all have theoretical limitations. Applying these approaches can bias genetic effect estimates and lead to seriously inflated type-I or type-II errors in conditional analysis, which is a critical tool for identifying independently associated variants. To address this challenge and complement imputation methods, we developed a method to combine summary statistics across participating studies and consistently estimate joint effects, even when the contributed summary statistics contain large amounts of missing values. Based on this estimator, we proposed a score statistic called PCBS (partial correlation based score statistic) for conditional analysis of single-variant and gene-level associations. Through extensive analysis of simulated and real data, we showed that the new method produces well-calibrated type-I errors and is substantially more powerful than existing approaches. We applied the proposed approach to one of the largest meta-analyses to date for the cigarettes-per-day phenotype. Using the new method, we identified multiple novel independently associated variants at known loci for tobacco use, which were otherwise missed by alternative methods. Together, the phenotypic variance explained by these variants was 1.1%, improving that of previously reported associations by 71%. These findings illustrate the extent of locus allelic heterogeneity and can help pinpoint causal variants.

It is of great interest to estimate the joint effects of multiple variants from large scale meta-analyses, in order to fine-map causal variants and understand the genetic architecture for complex traits. The summary association statistics from participating studies in a meta-analysis often contain missing values at some variant sites, as the imputation methods may not work well and the variants with low imputation quality will be filtered out. Missingness is especially likely when the underlying genetic variant is rare or the participating studies use targeted genotyping array that is not suitable for imputation. Existing methods for conditional meta-analysis do not properly handle missing data, and can incorrectly estimate correlations between score statistics. As a result, they can produce highly inflated type-I errors for conditional analysis, which will result in overestimated phenotypic variance explained and incorrect identification of causal variants. We systematically evaluated this bias and proposed a novel partial correlation based score statistic. The new statistic has valid type-I errors for conditional analysis and much higher power than the existing methods, even when the contributed summary statistics contain a large fraction of missing values. We expect this method to be highly useful in the sequencing age for complex trait genetics.

A sibling method for identifying vQTLs

The propensity of a trait to vary within a population may have evolutionary, ecological, or clinical significance. In the present study we deploy sibling models to offer a novel and unbiased way to ascertain loci associated with the extent to which phenotypes vary (variance-controlling quantitative trait loci, or vQTLs). Previous methods for vQTL-mapping either exclude genetically related individuals or treat genetic relatedness among individuals as a complicating factor addressed by adjusting estimates for non-independence in phenotypes. The present method uses genetic relatedness as a tool to obtain unbiased estimates of variance effects rather than as a nuisance. The family-based approach, which utilizes random variation between siblings in minor allele counts at a locus, also allows controls for parental genotype, mean effects, and non-linear (dominance) effects that may spuriously appear to generate variation. Simulations show that the approach performs equally well as two existing methods (squared Z-score and DGLM) in controlling type I error rates when there is no unobserved confounding, and performs significantly better than these methods in the presence of small degrees of confounding. Using height and BMI as empirical applications, we investigate SNPs that alter within-family variation in height and BMI, as well as pathways that appear to be enriched. One significant SNP for BMI variability, in the MAST4 gene, replicated. Pathway analysis revealed one gene set, encoding members of several signaling pathways related to gap junction function, which appears significantly enriched for associations with within-family height variation in both datasets (while not enriched in analysis of mean levels). We recommend approximating laboratory random assignment of genotype using family data and more careful attention to the possible conflation of mean and variance effects.

Adaptive SNP-Set Association Testing in Generalized Linear Mixed Models with Application to Family Studies

In genome-wide association studies (GWAS), it has been increasingly recognized that, as a complementary approach to standard single SNP analyses, it may be beneficial to analyze a group of functionally related SNPs together. Among the existent population-based SNP-set association tests, two adaptive tests, the aSPU test and the aSPUpath test, offer a powerful and general approach at the gene- and pathway-levels by data-adaptively combining the results across multiple SNPs (and genes) such that high statistical power can be maintained across a wide range of scenarios. We extend the aSPU and the aSPUpath test to familial data under the framework of the generalized linear mixed models (GLMMs), which can take account of both subject relatedness and possible population structure. As in population-based GWAS, the proposed aSPU and aSPUpath tests require only fitting a single and common GLMM (under the null hypothesis) for all the SNPs, thus are computationally efficient and feasible for large GWAS data. We illustrate our approaches in identifying genes and pathways associated with alcohol dependence in the Minnesota Twin Family Study. The aSPU test detected a gene associated with the trait, in contrast to none by the standard single SNP analysis. Our aSPU test also controlled Type I errors satisfactorily in a small simulation study. We provide R code to conduct the aSPU and aSPUpath tests for familial and other correlated data.

Resampling-based tests for Lasso in genome-wide association studies

Background

Genome-wide association studies involve detecting association between millions of genetic variants and a trait, which typically use univariate regression to test association between each single variant and the phenotype. Alternatively, Lasso penalized regression allows one to jointly model the relationship between all genetic variants and the phenotype. However, it is unclear how to best conduct inference on the individual Lasso coefficients, especially in high-dimensional settings.

Methods

We consider six methods for testing the Lasso coefficients: two permutation (Lasso-Ayers, Lasso-PL) and one analytic approach (Lasso-AL) to select the penalty parameter for type-1-error control, residual bootstrap (Lasso-RB), modified residual bootstrap (Lasso-MRB), and a permutation test (Lasso-PT). Methods are compared via simulations and application to the Minnesota Center for Twins and Family Study.

Results

We show that for finite sample sizes with increasing number of null predictors, Lasso-RB, Lasso-MRB, and Lasso-PT fail to be viable methods of inference. However, Lasso-PL and Lasso-AL remain fast and powerful tools for conducting inference with the Lasso, even in high-dimensions.

Conclusion

Our results suggest that the proposed permutation selection procedure (Lasso-PL) and the analytic selection method (Lasso-AL) are fast and powerful alternatives to the standard univariate analysis in genome-wide association studies.

Genetic tests in major psychiatric disorders-integrating molecular medicine with clinical psychiatry-why is it so difficult?

With the advent of post-genomic era, new technologies create extraordinary possibilities for diagnostics and personalized therapy, transforming todays' medicine. Rooted in both medical genetics and clinical psychiatry, the paper is designed as an integrated source of information of the current and potential future application of emerging genomic technologies as diagnostic tools in psychiatry, moving beyond the classical concept of patient approach. Selected approaches are presented, starting from currently used technologies (next-generation sequencing (NGS) and microarrays), followed by newer options (reverse phenotyping). Next, we describe an old concept in a new light (endophenotypes), subsequently coming up with a sophisticated and complex approach (gene networks) ending by a nascent field (computational psychiatry). The challenges and barriers that exist to translate genomic research to real-world patient assessment are further discussed. We emphasize the view that only a paradigm shift can bring a fundamental change in psychiatric practice, allowing to disentangle the intricacies of mental diseases. All the diagnostic methods, as described, are directed at uncovering the integrity of the system including many types of relations within a complex structure. The integrative system approach offers new opportunity to connect genetic background with specific diseases entities, or concurrently, with symptoms regardless of a diagnosis. To advance the field, we propose concerted cross-disciplinary effort to provide a diagnostic platform operating at the general level of genetic pathogenesis of complex-trait psychiatric disorders rather than at the individual level of a specific disease.

Psychophysiological endophenotypes to characterize mechanisms of known schizophrenia genetic loci

BACKGROUND

Endophenotypes are laboratory-based measures hypothesized to lie in the causal chain between genes and clinical disorder, and to serve as a more powerful way to identify genes associated with the disorder. One promise of endophenotypes is that they may assist in elucidating the neurobehavioral mechanisms by which an associated genetic polymorphism affects disorder risk in complex traits. We evaluated this promise by testing the extent to which variants discovered to be associated with schizophrenia through large-scale meta-analysis show associations with psychophysiological endophenotypes.

METHOD

We genome-wide genotyped and imputed 4905 individuals. Of these, 1837 were whole-genome-sequenced at 11× depth. In a community-based sample, we conducted targeted tests of variants within schizophrenia-associated loci, as well as genome-wide polygenic tests of association, with 17 psychophysiological endophenotypes including acoustic startle response and affective startle modulation, antisaccade, multiple frequencies of resting electroencephalogram (EEG), electrodermal activity and P300 event-related potential.

RESULTS

Using single variant tests and gene-based tests we found suggestive evidence for an association between contactin 4 (CNTN4) and antisaccade and P300. We were unable to find any other variant or gene within the 108 schizophrenia loci significantly associated with any of our 17 endophenotypes. Polygenic risk scores indexing genetic vulnerability to schizophrenia were not related to any of the psychophysiological endophenotypes after correction for multiple testing.

CONCLUSIONS

The results indicate significant difficulty in using psychophysiological endophenotypes to characterize the genetically influenced neurobehavioral mechanisms by which risk loci identified in genome-wide association studies affect disorder risk.

Deep Sequencing of 71 Candidate Genes to Characterize Variation Associated with Alcohol Dependence

BACKGROUND

Previous genomewide association studies (GWASs) have identified a number of putative risk loci for alcohol dependence (AD). However, only a few loci have replicated and these replicated variants only explain a small proportion of AD risk. Using an innovative approach, the goal of this study was to generate hypotheses about potentially causal variants for AD that can be explored further through functional studies.

METHODS

We employed targeted capture of 71 candidate loci and flanking regions followed by next-generation deep sequencing (mean coverage 78X) in 806 European Americans. Regions included in our targeted capture library were genes identified through published GWAS of alcohol, all human alcohol and aldehyde dehydrogenases, reward system genes including dopaminergic and opioid receptors, prioritized candidate genes based on previous associations, and genes involved in the absorption, distribution, metabolism, and excretion of drugs. We performed single-locus tests to determine if any single variant was associated with AD symptom count. Sets of variants that overlapped with biologically meaningful annotations were tested for association in aggregate.

RESULTS

No single, common variant was significantly associated with AD in our study. We did, however, find evidence for association with several variant sets. Two variant sets were significant at the q-value <0.10 level: a genic enhancer for ADHFE1 (p = 1.47 × 10^-5 ; q = 0.019), an alcohol dehydrogenase, and ADORA1 (p = 5.29 × 10^-5 ; q = 0.035), an adenosine receptor that belongs to a G-protein-coupled receptor gene family.

CONCLUSIONS

To our knowledge, this is the first sequencing study of AD to examine variants in entire genes, including flanking and regulatory regions. We found that in addition to protein coding variant sets, regulatory variant sets may play a role in AD. From these findings, we have generated initial functional hypotheses about how these sets may influence AD.

What can time-frequency and phase coherence measures tell us about the genetic basis of P3 amplitude?

In a recent comprehensive investigation, we largely failed to identify significant genetic markers associated with P3 amplitude or to corroborate previous associations between P3 and specific single nucleotide polymorphisms (SNPs) or genes. In the present study we extended this line of investigation to examine time-frequency (TF) activity and intertrial phase coherence (ITPC) in the P3 time window, both of which are associated with P3 amplitude. Previous genome-wide research has reported associations between P3-related theta and delta activity and individual genetic variants. A large, population-based sample of 4211 subjects, comprising male and female adolescent twins and their parents, was genotyped for 527,828 single nucleotide polymorphisms (SNPs), from which over six million SNPs were accurately imputed. Heritability estimates were greater for TF energy than ITPC, whether based on biometric models or the combined influence of all measured SNPs (derived from genome-wide complex trait analysis). The magnitude of overlap in the specific SNPs associated with delta energy and ITPC and P3 amplitude was significant. A genome-wide analysis of all SNPs, accompanied by an analysis of approximately 17,600 genes, indicated a region of chromosome 2 around TEKT4 that was significantly associated with theta ITPC. Analysis of candidate SNPs and genes previously reported to be associated with P3 or related phenotypes yielded one association surviving correction for multiple tests: between theta energy and CRHR1. However, we did not obtain significant associations for SNPs implicated in previous genome-wide studies of TF measures. Identifying specific genetic variants associated with P3 amplitude remains a challenge.

A Test-Replicate Approach to Candidate Gene Research on Addiction and Externalizing Disorders: A Collaboration Across Five Longitudinal Studies

This study presents results from a collaboration across five longitudinal studies seeking to test and replicate models of gene-environment interplay in the development of substance use and externalizing disorders (SUDs, EXT). We describe an overview of our conceptual models, plan for gene-environment interplay analyses, and present main effects results evaluating six candidate genes potentially relevant to SUDs and EXT (MAOA, 5-HTTLPR, COMT, DRD2, DAT1, and DRD4). All samples included rich longitudinal and phenotypic measurements from childhood/adolescence (ages 5-13) through early adulthood (ages 25-33); sample sizes ranged from 3487 in the test sample, to ~600-1000 in the replication samples. Phenotypes included lifetime symptom counts of SUDs (nicotine, alcohol and cannabis), adult antisocial behavior, and an aggregate externalizing disorder composite. Covariates included the first 10 ancestral principal components computed using all autosomal markers in subjects across the data sets, and age at the most recent assessment. Sex, ancestry, and exposure effects were thoroughly evaluated. After correcting for multiple testing, only one significant main effect was found in the test sample, but it was not replicated. Implications for subsequent gene-environment interplay analyses are discussed.

Identification of Common Genetic Variants Influencing Spontaneous Dizygotic Twinning and Female Fertility

Spontaneous dizygotic (DZ) twinning occurs in 1%-4% of women, with familial clustering and unknown physiological pathways and genetic origin. DZ twinning might index increased fertility and has distinct health implications for mother and child. We performed a GWAS in 1,980 mothers of spontaneous DZ twins and 12,953 control subjects. Findings were replicated in a large Icelandic cohort and tested for association across a broad range of fertility traits in women. Two SNPs were identified (rs11031006 near FSHB, p = 1.54 × 10(-9), and rs17293443 in SMAD3, p = 1.57 × 10(-8)) and replicated (p = 3 × 10(-3) and p = 1.44 × 10(-4), respectively). Based on ∼90,000 births in Iceland, the risk of a mother delivering twins increased by 18% for each copy of allele rs11031006-G and 9% for rs17293443-C. A higher polygenic risk score (PRS) for DZ twinning, calculated based on the results of the DZ twinning GWAS, was significantly associated with DZ twinning in Iceland (p = 0.001). A higher PRS was also associated with having children (p = 0.01), greater lifetime parity (p = 0.03), and earlier age at first child (p = 0.02). Allele rs11031006-G was associated with higher serum FSH levels, earlier age at menarche, earlier age at first child, higher lifetime parity, lower PCOS risk, and earlier age at menopause. Conversely, rs17293443-C was associated with later age at last child. We identified robust genetic risk variants for DZ twinning: one near FSHB and a second within SMAD3, the product of which plays an important role in gonadal responsiveness to FSH. These loci contribute to crucial aspects of reproductive capacity and health.

Association of the OPRM1 Variant rs1799971 (A118G) with Non-Specific Liability to Substance Dependence in a Collaborative de novo Meta-Analysis of European-Ancestry Cohorts

The mu1 opioid receptor gene, OPRM1, has long been a high-priority candidate for human genetic studies of addiction. Because of its potential functional significance, the non-synonymous variant rs1799971 (A118G, Asn40Asp) in OPRM1 has been extensively studied, yet its role in addiction has remained unclear, with conflicting association findings. To resolve the question of what effect, if any, rs1799971 has on substance dependence risk, we conducted collaborative meta-analyses of 25 datasets with over 28,000 European-ancestry subjects. We investigated non-specific risk for "general" substance dependence, comparing cases dependent on any substance to controls who were non-dependent on all assessed substances. We also examined five specific substance dependence diagnoses: DSM-IV alcohol, opioid, cannabis, and cocaine dependence, and nicotine dependence defined by the proxy of heavy/light smoking (cigarettes-per-day >20 vs. ≤ 10). The G allele showed a modest protective effect on general substance dependence (OR = 0.90, 95% C.I. [0.83-0.97], p value = 0.0095, N = 16,908). We observed similar effects for each individual substance, although these were not statistically significant, likely because of reduced sample sizes. We conclude that rs1799971 contributes to mechanisms of addiction liability that are shared across different addictive substances. This project highlights the benefits of examining addictive behaviors collectively and the power of collaborative data sharing and meta-analyses.

Heritability and molecular-genetic basis of the P3 event-related brain potential: a genome-wide association study

P3 amplitude is a candidate endophenotype for disinhibitory psychopathology, psychosis, and other disorders. The present study is a comprehensive analysis of the behavioral- and molecular-genetic basis of P3 amplitude and a P3 genetic factor score in a large community sample (N = 4,211) of adolescent twins and their parents, genotyped for 527,829 single nucleotide polymorphisms (SNPs). Biometric models indicated that as much as 65% of the variance in each measure was due to additive genes. All SNPs in aggregate accounted for approximately 40% to 50% of the heritable variance. However, analyses of individual SNPs did not yield any significant associations. Analyses of individual genes did not confirm previous associations between P3 amplitude and candidate genes but did yield a novel association with myelin expression factor 2 (MYEF2). Main effects of individual variants may be too small to be detected by GWAS without larger samples.

Heritability and molecular-genetic basis of resting EEG activity: a genome-wide association study

Several EEG parameters are potential endophenotypes for different psychiatric disorders. The present study consists of a comprehensive behavioral- and molecular-genetic analysis of such parameters in a large community sample (N = 4,026) of adolescent twins and their parents, genotyped for 527,829 single nucleotide polymorphisms (SNPs). Biometric heritability estimates ranged from .49 to .85, with a median of .78. The additive effect of all SNPs (SNP heritability) varied across electrodes. Although individual SNPs were not significantly associated with EEG parameters, several genes were associated with delta power. We also obtained an association between the GABRA2 gene and beta power (p < .014), consistent with findings reported by others, although this did not survive Bonferroni correction. If EEG parameters conform to a largely polygenic model of inheritance, larger sample sizes will be required to detect individual variants reliably.

Heritability and molecular genetic basis of electrodermal activity: a genome-wide association study

The molecular genetic basis of electrodermal activity (EDA) was analyzed using 527,829 single nucleotide polymorphisms (SNPs) in a large population-representative sample of twins and parents (N = 4,424) in relation to various EDA indices. Biometric analyses suggested that approximately 50% or more of variance in all EDA indices was heritable. The combined effect of all SNPs together accounted for a significant amount of variance in each index, affirming their polygenic basis and heritability. However, none of the SNPs were genome-wide significant for any EDA index. Previously reported SNP associations with disorders such as substance dependence or schizophrenia, which have been linked to EDA abnormalities, were not significant; nor were associations between EDA and genes in specific neurotransmitter systems. These results suggest that EDA is influenced by multiple genes rather than by polymorphisms with large effects.

Premorbid risk factors for major depressive disorder: are they associated with early onset and recurrent course?

Premorbid risk for major depressive disorder (MDD) and predictors of an earlier onset and recurrent course were examined in two studies in a large, community-based sample of parents and offspring, prospectively assessed from late childhood into adulthood. In Study 1 (N = 2,764 offspring and their parents), parental psychiatric status, offspring personality at age 11, and age 11 offspring internalizing and externalizing symptoms predicted the subsequent development of MDD, as did poor quality parent-child relationships, poor academic functioning, early pubertal development, and childhood maltreatment by age 11. Parental MDD and adult antisocial behavior, offspring negative emotionality and disconstraint, externalizing symptoms, and childhood maltreatment predicted an earlier onset of MDD, after accounting for course; lower positive emotionality, trait anxiety, and childhood maltreatment predicted recurrent MDD, after accounting for age of onset. In Study 2 (N = 7,146), we examined molecular genetic risk for MDD by extending recent reports of associations with glutamatergic system genes. We failed to confirm associations with MDD using either individual single nucleotide polymorphism based tests or gene-based analyses. Overall, results speak to the pervasiveness of risk for MDD, as well as specific risk for early onset MDD; risk for recurrent MDD appears to be largely a function of its often earlier onset.

Genome-Wide Association Study of Behavioral Disinhibition in a Selected Adolescent Sample

Behavioral disinhibition (BD) is a quantitative measure designed to capture the heritable variation encompassing risky and impulsive behaviors. As a result, BD represents an ideal target for discovering genetic loci that predispose individuals to a wide range of antisocial behaviors and substance misuse that together represent a large cost to society as a whole. Published genome-wide association studies (GWAS) have examined specific phenotypes that fall under the umbrella of BD (e.g. alcohol dependence, conduct disorder); however no GWAS has specifically examined the overall BD construct. We conducted a GWAS of BD using a sample of 1,901 adolescents over-selected for characteristics that define high BD, such as substance and antisocial behavior problems, finding no individual locus that surpassed genome-wide significance. Although no single SNP was significantly associated with BD, restricted maximum likelihood analysis estimated that 49.3 % of the variance in BD within the Caucasian sub-sample was accounted for by the genotyped SNPs (p = 0.06). Gene-based tests identified seven genes associated with BD (p ≤ 2.0 × 10(-6)). Although the current study was unable to identify specific SNPs or pathways with replicable effects on BD, the substantial sample variance that could be explained by all genotyped SNPs suggests that larger studies could successfully identify common variants associated with BD.

Determination of Zygosity in Adult Chinese Twins Using the 450K Methylation Array versus Questionnaire Data

Previous studies have shown that both single nucleotide polymorphisms (SNPs) and questionnaires-based method can be used for twin zygosity determination, but few validation studies have been conducted using Chinese populations. In the current study, we recruited 192 same sex Chinese adult twin pairs to evaluate the validity of using genetic markers-based method and questionnaire-based method in zygosity determination. We considered the relatedness analysis based on more than 0.6 million SNPs genotyping as the golden standards for zygosity determination. After quality control, qualified twins were left for relatedness analysis based on identical by descent calculation. Then those same sex twin pairs were included in the zygosity questionnaire validation analysis. Logistic regression model was applied to assess the discriminant ability of age, sex and the three questions in zygosity determination. Leave one out cross-validation was used as a measurement of internal validation. The results of zygosity determination based on 65 SNPs in 450k methylation array were all consistent with genotyping. Age, gender, questions of appearance confused by strangers and previously perceived zygosity consisted of the most predictable model with a consistency rate of 0.8698, cross validation predictive error of 0.1347. For twin studies with genotyping and\or 450k methylation array, there would be no need to conduct other zygosity testing for the sake of costs consideration.

Heritability and molecular genetic basis of antisaccade eye tracking error rate: a genome-wide association study

Antisaccade deficits reflect abnormalities in executive function linked to various disorders including schizophrenia, externalizing psychopathology, and neurological conditions. We examined the genetic bases of antisaccade error in a sample of community-based twins and parents (N = 4,469). Biometric models showed that about half of the variance in the antisaccade response was due to genetic factors and half due to nonshared environmental factors. Molecular genetic analyses supported these results, showing that the heritability accounted for by common molecular genetic variants approximated biometric estimates. Genome-wide analyses revealed several SNPs as well as two genes-B3GNT7 and NCL-on Chromosome 2 associated with antisaccade error. SNPs and genes hypothesized to be associated with antisaccade error based on prior work, although generating some suggestive findings for MIR137, GRM8, and CACNG2, could not be confirmed.

Genome-wide scans of genetic variants for psychophysiological endophenotypes: a methodological overview

This article provides an introductory overview of the investigative strategy employed to evaluate the genetic basis of 17 endophenotypes examined as part of a 20-year data collection effort from the Minnesota Center for Twin and Family Research. Included are characterization of the study samples, descriptive statistics for key properties of the psychophysiological measures, and rationale behind the steps taken in the molecular genetic study design. The statistical approach included (a) biometric analysis of twin and family data, (b) heritability analysis using 527,829 single nucleotide polymorphisms (SNPs), (c) genome-wide association analysis of these SNPs and 17,601 autosomal genes, (d) follow-up analyses of candidate SNPs and genes hypothesized to have an association with each endophenotype, (e) rare variant analysis of nonsynonymous SNPs in the exome, and (f) whole genome sequencing association analysis using 27 million genetic variants. These methods were used in the accompanying empirical articles comprising this special issue, Genome-Wide Scans of Genetic Variants for Psychophysiological Endophenotypes.

Heritability and molecular genetic basis of acoustic startle eye blink and affectively modulated startle response: a genome-wide association study

Acoustic startle responses have been studied extensively in relation to individual differences and psychopathology. We examined three indices of the blink response in a picture-viewing paradigm-overall startle magnitude across all picture types, and aversive and pleasant modulation scores-in 3,323 twins and parents. Biometric models and molecular genetic analyses showed that half the variance in overall startle was due to additive genetic effects. No single nucleotide polymorphism was genome-wide significant, but GRIK3 produced a significant effect when examined as part of a candidate gene set. In contrast, emotion modulation scores showed little evidence of heritability in either biometric or molecular genetic analyses. However, in a genome-wide scan, PARP14 produced a significant effect for aversive modulation. We conclude that, although overall startle retains potential as an endophenotype, emotion-modulated startle does not.

Genetic associations of nonsynonymous exonic variants with psychophysiological endophenotypes

We mapped ∼85,000 rare nonsynonymous exonic single nucleotide polymorphisms (SNPs) to 17 psychophysiological endophenotypes in 4,905 individuals, including antisaccade eye movements, resting EEG, P300 amplitude, electrodermal activity, affect-modulated startle eye blink. Nonsynonymous SNPs are predicted to directly change or disrupt proteins encoded by genes and are expected to have significant biological consequences. Most such variants are rare, and new technologies can efficiently assay them on a large scale. We assayed 247,870 mostly rare SNPs on an Illumina exome array. Approximately 85,000 of the SNPs were polymorphic, rare (MAF < .05), and nonsynonymous. Single variant association tests identified a SNP in the PARD3 gene associated with theta resting EEG power. The sequence kernel association test, a gene-based test, identified a gene PNPLA7 associated with pleasant difference startle, the difference in startle magnitude between pleasant and neutral images. No other single nonsynonymous variant, or gene-based group of variants, was strongly associated with any endophenotype.

In search of rare variants: preliminary results from whole genome sequencing of 1,325 individuals with psychophysiological endophenotypes

Whole genome sequencing was completed on 1,325 individuals from 602 families, identifying 27 million autosomal variants. Genetic association tests were conducted for those individuals who had been assessed for one or more of 17 endophenotypes (N range = 802-1,185). No significant associations were found. These 27 million variants were then imputed into the full sample of individuals with psychophysiological data (N range = 3,088-4,469) and again tested for associations with the 17 endophenotypes. No association was significant. Using a gene-based variable threshold burden test of nonsynonymous variants, we obtained five significant associations. These findings are preliminary and call for additional analysis of this rich sample. We argue that larger samples, alternative study designs, and additional bioinformatics approaches will be necessary to discover associations between these endophenotypes and genomic variation.

Results of a "GWAS plus:" general cognitive ability is substantially heritable and massively polygenic

We carried out a genome-wide association study (GWAS) for general cognitive ability (GCA) plus three other analyses of GWAS data that aggregate the effects of multiple single-nucleotide polymorphisms (SNPs) in various ways. Our multigenerational sample comprised 7,100 Caucasian participants, drawn from two longitudinal family studies, who had been assessed with an age-appropriate IQ test and had provided DNA samples passing quality screens. We conducted the GWAS across ∼2.5 million SNPs (both typed and imputed), using a generalized least-squares method appropriate for the different family structures present in our sample, and subsequently conducted gene-based association tests. We also conducted polygenic prediction analyses under five-fold cross-validation, using two different schemes of weighting SNPs. Using parametric bootstrapping, we assessed the performance of this prediction procedure under the null. Finally, we estimated the proportion of variance attributable to all genotyped SNPs as random effects with software GCTA. The study is limited chiefly by its power to detect realistic single-SNP or single-gene effects, none of which reached genome-wide significance, though some genomic inflation was evident from the GWAS. Unit SNP weights performed about as well as least-squares regression weights under cross-validation, but the performance of both increased as more SNPs were included in calculating the polygenic score. Estimates from GCTA were 35% of phenotypic variance at the recommended biological-relatedness ceiling. Taken together, our results concur with other recent studies: they support a substantial heritability of GCA, arising from a very large number of causal SNPs, each of very small effect. We place our study in the context of the literature–both contemporary and historical–and provide accessible explication of our statistical methods.