Declining autozygosity over time: An exploration in over 1 million individuals from three diverse cohorts.

Previous studies have hypothesized that autozygosity is decreasing over generational time. However, these studies were limited to relatively small samples (n < 11,000) lacking in diversity, which may limit the generalizability of their findings. We present data that partially support this hypothesis from three large cohorts of diverse ancestries, two from the US (All of Us, n = 82,474; the Million Veteran Program, n = 622,497) and one from the UK (UK Biobank, n = 380,899). Our results from a mixed-effect meta-analysis demonstrate an overall trend of decreasing autozygosity over generational time (meta-analyzed slope = -0.029, SE = 0.009, p = 6.03e-4). On the basis of our estimates, we would predict FROH to decline 0.29% for every 20-year increase in birth year. We determined that a model including an ancestry-by-country interaction term fit the data best, indicating that ancestry differences in this trend differ by country. We found further evidence to suggest a difference between the US and UK cohorts by meta-analyzing within country, observing a significant negative estimate in the US cohorts (meta-analyzed slope = -0.058, SE = 0.015, p = 1.50e-4) but a non-significant estimate in the UK (meta-analyzed slope = -0.001, SE = 0.008, p = 0.945). The association between autozygosity and birth year was substantially attenuated when accounting for educational attainment and income (meta-analyzed slope = -0.011, SE = 0.008, p = 0.167), suggesting they may partially account for decreasing autozygosity over time. Overall, we demonstrate decreasing autozygosity over time in a large, modern sample and speculate that this trend can be attributed to increases in urbanization and panmixia and differences in sociodemographic processes lead to country-specific differences in the rate of decline.

Genomic partitioning of inbreeding depression in humans.

Across species, offspring of related individuals often exhibit significant reduction in fitness-related traits, known as inbreeding depression (ID), yet the genetic and molecular basis for ID remains elusive. Here, we develop a method to quantify enrichment of ID within specific genomic annotations and apply it to human data. We analyzed the phenomes and genomes of ∼350,000 unrelated participants of the UK Biobank and found, on average of over 11 traits, significant enrichment of ID within genomic regions with high recombination rates (>21-fold; p < 10-5), with conserved function across species (>19-fold; p < 10-4), and within regulatory elements such as DNase I hypersensitive sites (∼5-fold; p = 8.9 × 10-7). We also quantified enrichment of ID within trait-associated regions and found suggestive evidence that genomic regions contributing to additive genetic variance in the population are enriched for ID signal. We find strong correlations between functional enrichment of SNP-based heritability and that of ID (r = 0.8, standard error: 0.1). These findings provide empirical evidence that ID is most likely due to many partially recessive deleterious alleles in low linkage disequilibrium regions of the genome. Our study suggests that functional characterization of ID may further elucidate the genetic architectures and biological mechanisms underlying complex traits and diseases.

Predicting causal genes from psychiatric genome-wide association studies using high-level etiological knowledge.

Genome-wide association studies have discovered hundreds of genomic loci associated with psychiatric traits, but the causal genes underlying these associations are often unclear, a research gap that has hindered clinical translation. Here, we present a Psychiatric Omnilocus Prioritization Score (PsyOPS) derived from just three binary features encapsulating high-level assumptions about psychiatric disease etiology - namely, that causal psychiatric disease genes are likely to be mutationally constrained, be specifically expressed in the brain, and overlap with known neurodevelopmental disease genes. To our knowledge, PsyOPS is the first method specifically tailored to prioritizing causal genes at psychiatric GWAS loci. We show that, despite its extreme simplicity, PsyOPS achieves state-of-the-art performance at this task, comparable to a prior domain-agnostic approach relying on tens of thousands of features. Genes prioritized by PsyOPS are substantially more likely than other genes at the same loci to have convergent evidence of direct regulation by the GWAS variant according to both DNA looping assays and expression or splicing quantitative trait locus (QTL) maps. We provide examples of genes hundreds of kilobases away from the lead variant, like GABBR1 for schizophrenia, that are prioritized by all three of PsyOPS, DNA looping and QTLs. Our results underscore the power of incorporating high-level knowledge of trait etiology into causal gene prediction at GWAS loci, and comprise a resource for researchers interested in experimentally characterizing psychiatric gene candidates.

Disentangling the influences of parental genetics on offspring's cognition, education, and psychopathology via genetic and phenotypic pathways.

BACKGROUND: Specific pathways of intergenerational transmission of behavioral traits remain unclear. Here, we aim to investigate how parental genetics influence offspring cognition, educational attainment, and psychopathology in youth. METHODS: Participants for the discovery sample were 2,189 offspring (aged 6-14 years), 1898 mothers and 1,017 fathers who underwent genotyping, psychiatric, and cognitive assessments. We calculated polygenic scores (PGS) for cognition, educational attainment, attention-deficit hyperactivity disorder (ADHD), and schizophrenia for the trios. Phenotypes studied included educational and cognitive measures, ADHD and psychotic symptoms. We used a stepwise approach and multiple mediation models to analyze the effect of parental PGS on offspring traits via offspring PGS and parental phenotype. Significant results were replicated in a sample of 1,029 adolescents, 363 mothers, and 307 fathers. RESULTS: Maternal and paternal PGS for cognition influenced offspring general intelligence and executive function via offspring PGS (genetic pathway) and parental education (phenotypic pathway). Similar results were found for parental PGS for educational attainment and offspring reading and writing skills. These pathways fully explained associations between parental PGS and offspring phenotypes, without residual direct association. Associations with maternal, but not paternal, PGS were replicated. No associations were found between parental PGS for psychopathology and offspring specific symptoms. CONCLUSIONS: Our findings indicate that parental genetics influences offspring cognition and educational attainment by genetic and phenotypic pathways, suggesting the expression of parental phenotypes partially explain the association between parental genetic risk and offspring outcomes. Multiple mediations might represent an effective approach to disentangle distinct pathways for intergenerational transmission of behavioral traits.

The challenges of sociogenomics make it more, not less, worthy of careful and innovative investigation.

Influences on social traits involve a tangled interplay of genetic, social, and environmental factors. Moreover, there is increasing awareness that gene-environment correlations are real and potentially measurable. Such gene-environment correlations can mislead if they are uncontrolled and genetic associations are interpreted as being purely because of direct genetic effects. This complexity is cause for more and better investigation, not a reason to refrain from researching one of the potentially important factors (genetics) influencing trait variation.

Novel approach for parallelizing pairwise comparison problems as applied to detecting segments identical by decent in whole-genome data.

MOTIVATION: Pairwise comparison problems arise in many areas of science. In genomics, datasets are already large and getting larger, and so operations that require pairwise comparisons-either on pairs of SNPs or pairs of individuals-are extremely computationally challenging. We propose a generic algorithm for addressing pairwise comparison problems that breaks a large problem (of order n2 comparisons) into multiple smaller ones (each of order n comparisons), allowing for massive parallelization. RESULTS: We demonstrated that this approach is very efficient for calling identical by descent (IBD) segments between all pairs of individuals in the UK Biobank dataset, with a 250-fold savings in time and 750-fold savings in memory over the standard approach to detecting such segments across the full dataset. This efficiency should extend to other methods of IBD calling and, more generally, to other pairwise comparison tasks in genomics or other areas of science. AVAILABILITY AND IMPLEMENTATION: A GitHub page is available at https://github.com/emmanuelsapin with the code to generate data needed for the implementation.

Exploring the Relationships Between Autozygosity, Educational Attainment, and Cognitive Ability in a Contemporary, Trans-Ancestral American Sample.

Previous studies have found significant associations between estimated autozygosity - the proportion of an individual's genome contained in homozygous segments due to distant inbreeding - and multiple traits, including educational attainment (EA) and cognitive ability. In one study, estimated autozygosity showed a stronger association with parental EA than the subject's own EA. This was likely driven by parental EA's association with mobility: more educated parents tended to migrate further from their hometown, and because of the strong correlation between ancestry and geography in the Netherlands, these individuals chose partners farther from their ancestry and therefore more different from them genetically. We examined the associations between estimated autozygosity, cognitive ability, and parental EA in a contemporary sub-sample of adolescents from the Adolescent Brain Cognitive Development Study℠ (ABCD Study®) (analytic N = 6,504). We found a negative association between autozygosity and child cognitive ability consistent with previous studies, while the associations between autozygosity and parental EA were in the expected direction of effect (with greater levels of autozygosity being associated with lower EA) but the effect sizes were significantly weaker than those estimated in previous work. We also found a lower mean level of autozygosity in the ABCD sample compared to previous autozygosity studies, which may reflect overall decreasing levels of autozygosity over generations. Variation in spousal similarities in ancestral background in the ABCD study compared to other studies may explain the pattern of associations between estimated autozygosity, EA, and cognitive ability in the current study.

The estimation of environmental and genetic parental influences.

Parents share half of their genes with their children, but they also share background social factors and actively help shape their child's environment - making it difficult to disentangle genetic and environmental causes of parent-offspring similarity. While adoption and extended twin family designs have been extremely useful for distinguishing genetic and nongenetic parental influences, these designs entail stringent assumptions about phenotypic similarity between relatives and require samples that are difficult to collect and therefore are typically small and not publicly shared. Here, we describe these traditional designs, as well as modern approaches that use large, publicly available genome-wide data sets to estimate parental effects. We focus in particular on an approach we recently developed, structural equation modeling (SEM)-polygenic score (PGS), that instantiates the logic of modern PGS-based methods within the flexible SEM framework used in traditional designs. Genetically informative designs such as SEM-PGS rely on different and, in some cases, less rigid assumptions than traditional approaches; thus, they allow researchers to capitalize on new data sources and answer questions that could not previously be investigated. We believe that SEM-PGS and similar approaches can lead to improved insight into how nature and nurture combine to create the incredible diversity underlying human behavior.

Estimation of Parental Effects Using Polygenic Scores.

Offspring resemble their parents for both genetic and environmental reasons. Understanding the relative magnitude of these alternatives has long been a core interest in behavioral genetics research, but traditional designs, which compare phenotypic covariances to make inferences about unmeasured genetic and environmental factors, have struggled to disentangle them. Recently, Kong et al. (2018) showed that by correlating offspring phenotypic values with the measured polygenic score of parents' nontransmitted alleles, one can estimate the effect of "genetic nurture"-a type of passive gene-environment covariation that arises when heritable parental traits directly influence offspring traits. Here, we instantiate this basic idea in a set of causal models that provide novel insights into the estimation of parental influences on offspring. Most importantly, we show how jointly modeling the parental polygenic scores and the offspring phenotypes can provide an unbiased estimate of the variation attributable to the environmental influence of parents on offspring, even when the polygenic score accounts for a small fraction of trait heritability. This model can be further extended to (a) account for the influence of different types of assortative mating, (b) estimate the total variation due to additive genetic effects and their covariance with the familial environment (i.e., the full genetic nurture effect), and (c) model situations where a parental trait influences a different offspring trait. By utilizing structural equation modeling techniques developed for extended twin family designs, our approach provides a general framework for modeling polygenic scores in family studies and allows for various model extensions that can be used to answer old questions about familial influences in new ways.

Bias and Precision of Parameter Estimates from Models Using Polygenic Scores to Estimate Environmental and Genetic Parental Influences.

In a companion paper Balbona et al. (Behav Genet, in press), we introduced a series of causal models that use polygenic scores from transmitted and nontransmitted alleles, the offspring trait, and parental traits to estimate the variation due to the environmental influences the parental trait has on the offspring trait (vertical transmission) as well as additive genetic effects. These models also estimate and account for the gene-gene and gene-environment covariation that arises from assortative mating and vertical transmission respectively. In the current study, we simulated polygenic scores and phenotypes of parents and offspring under genetic and vertical transmission scenarios, assuming two types of assortative mating. We instantiated the models from our companion paper in the OpenMx software, and compared the true values of parameters to maximum likelihood estimates from models fitted on the simulated data to quantify the bias and precision of estimates. We show that parameter estimates from these models are unbiased when assumptions are met, but as expected, they are biased to the degree that assumptions are unmet. Standard errors of the estimated variances due to vertical transmission and to genetic effects decrease with increasing sample sizes and with increasing [Formula: see text] values of the polygenic score. Even when the polygenic score explains a modest amount of trait variation ([Formula: see text]), standard errors of these standardized estimates are reasonable ([Formula: see text]) for [Formula: see text] trios, and can even be reasonable for smaller sample sizes (e.g., down to 4K) when the polygenic score is more predictive. These causal models offer a novel approach for understanding how parents influence their offspring, but their use requires polygenic scores on relevant traits that are modestly predictive (e.g., [Formula: see text] as well as datasets with genomic and phenotypic information on parents and offspring. The utility of polygenic scores for elucidating parental influences should thus serve as additional motivation for large genomic biobanks to perform GWAS's on traits that may be relevant to parenting and to oversample close relatives, particularly parents and offspring.

Relationships between estimated autozygosity and complex traits in the UK Biobank.

Inbreeding increases the risk of certain Mendelian disorders in humans but may also reduce fitness through its effects on complex traits and diseases. Such inbreeding depression is thought to occur due to increased homozygosity at causal variants that are recessive with respect to fitness. Until recently it has been difficult to amass large enough sample sizes to investigate the effects of inbreeding depression on complex traits using genome-wide single nucleotide polymorphism (SNP) data in population-based samples. Further, it is difficult to infer causation in analyses that relate degree of inbreeding to complex traits because confounding variables (e.g., education) may influence both the likelihood for parents to outbreed and offspring trait values. The present study used runs of homozygosity in genome-wide SNP data in up to 400,000 individuals in the UK Biobank to estimate the proportion of the autosome that exists in autozygous tracts-stretches of the genome which are identical due to a shared common ancestor. After multiple testing corrections and controlling for possible sociodemographic confounders, we found significant relationships in the predicted direction between estimated autozygosity and three of the 26 traits we investigated: age at first sexual intercourse, fluid intelligence, and forced expiratory volume in 1 second. Our findings corroborate those of several published studies. These results may imply that these traits have been associated with Darwinian fitness over evolutionary time. However, some of the autozygosity-trait relationships were attenuated after controlling for background sociodemographic characteristics, suggesting that alternative explanations for these associations have not been eliminated. Care needs to be taken in the design and interpretation of ROH studies in order to glean reliable information about the genetic architecture and evolutionary history of complex traits.

No Evidence for Social Genetic Effects or Genetic Similarity Among Friends Beyond that Due to Population Stratification: A Reappraisal of Domingue et al (2018).

Using data from 5500 adolescents from the National Longitudinal Study of Adolescent to Adult Health, Domingue et al. (Proc Natl Acad Sci 25:256., 2018) claimed to show that friends are genetically more similar to one another than randomly selected peers, beyond the confounding effects of population stratification by ancestry. The authors also claimed to show 'social-genetic' effects, whereby individuals' educational attainment (EA) is influenced by their friends' genes. We argue that neither claim is justified by the data. Mathematically we show that (1) the genetic similarity reported between friends is far larger than theoretically possible if it was caused by phenotypic assortment as the authors claim; uncontrolled population stratification is a likely reason for the genetic similarity they observed, and (2) significant association between individuals' EA and their friends' polygenic scores for EA is a necessary consequence of EA similarity among friends, and does not provide evidence for social-genetic effects. Going forward, we urge caution in the analysis and interpretation of data at the intersection of human genetics and the social sciences.

The Role of A Priori-Identified Addiction and Smoking Gene Sets in Smoking Behaviors.

INTRODUCTION: Smoking is a leading cause of death, and genetic variation contributes to smoking behaviors. Identifying genes and sets of genes that contribute to risk for addiction is necessary to prioritize targets for functional characterization and for personalized medicine. METHODS: We performed a gene set-based association and heritable enrichment study of two addiction-related gene sets, those on the Smokescreen Genotyping Array and the nicotinic acetylcholine receptors, using the largest available GWAS summary statistics. We assessed smoking initiation, cigarettes per day, smoking cessation, and age of smoking initiation. RESULTS: Individual genes within each gene set were significantly associated with smoking behaviors. Both sets of genes were significantly associated with cigarettes per day, smoking initiation, and smoking cessation. Age of initiation was only associated with the Smokescreen gene set. Although both sets of genes were enriched for trait heritability, each accounts for only a small proportion of the single nucleotide polymorphism-based heritability (2%-12%). CONCLUSIONS: These two gene sets are associated with smoking behaviors, but collectively account for a limited amount of the genetic and phenotypic variation of these complex traits, consistent with high polygenicity. IMPLICATIONS: We evaluated evidence for the association and heritable contribution of expert-curated and bioinformatically identified sets of genes related to smoking. Although they impact smoking behaviors, these specifically targeted genes do not account for much of the heritability in smoking and will be of limited use for predictive purposes. Advanced genome-wide approaches and integration of other 'omics data will be needed to fully account for the genetic variation in smoking phenotypes.

Nick Martin as a Mentor - A Perspective.

Nick Martin has had an outsized influence on the field of behavioral genetics. Much of this influence stems from his mentorship of young scientists. I describe Nick's mentorship, and what makes it special, from a personal perspective.

Shared additive genetic variation for alcohol dependence among subjects of African and European ancestry.

Alcohol dependence (AD) affects individuals from all racial/ethnic groups, and previous research suggests that there is considerable variation in AD risk between and among various ancestrally defined groups in the United States. Although the reasons for these differences are likely due in part to contributions of complex sociocultural factors, limited research has attempted to examine whether similar genetic variation plays a role across ancestral groups. Using a pooled sample of individuals of African and European ancestry (AA/EA) obtained through data shared within the Database for Genotypes and Phenotypes, we estimated the extent to which additive genetic similarity for AD between AA and EAs using common single nucleotide polymorphisms overlapped across the two populations. AD was represented as a factor score by using Diagnostic and Statistical Manual dependence criteria, and genetic data were imputed by using the 1000 Genomes Reference Panel. Analyses revealed a significant single nucleotide polymorphism-based heritability of 17 percent (SE = 5) in EAs and 24 percent (SE = 15) in AAs. Further, a significant genetic correlation of 0.77 (SE = 0.46) suggests that the allelic architecture influencing the AD factor for EAs and AAs is largely similar across the two populations. Analyses indicated that investigating the genetic underpinnings of alcohol dependence in different ethnic groups may serve to highlight core etiological factors common to both groups and unique etiological factors specific to each ethnic group.

No Reliable Association between Runs of Homozygosity and Schizophrenia in a Well-Powered Replication Study.

It is well known that inbreeding increases the risk of recessive monogenic diseases, but it is less certain whether it contributes to the etiology of complex diseases such as schizophrenia. One way to estimate the effects of inbreeding is to examine the association between disease diagnosis and genome-wide autozygosity estimated using runs of homozygosity (ROH) in genome-wide single nucleotide polymorphism arrays. Using data for schizophrenia from the Psychiatric Genomics Consortium (n = 21,868), Keller et al. (2012) estimated that the odds of developing schizophrenia increased by approximately 17% for every additional percent of the genome that is autozygous (ß = 16.1, CI(ß) = [6.93, 25.7], Z = 3.44, p = 0.0006). Here we describe replication results from 22 independent schizophrenia case-control datasets from the Psychiatric Genomics Consortium (n = 39,830). Using the same ROH calling thresholds and procedures as Keller et al. (2012), we were unable to replicate the significant association between ROH burden and schizophrenia in the independent PGC phase II data, although the effect was in the predicted direction, and the combined (original + replication) dataset yielded an attenuated but significant relationship between Froh and schizophrenia (ß = 4.86,CI(ß) = [0.90,8.83],Z = 2.40,p = 0.02). Since Keller et al. (2012), several studies reported inconsistent association of ROH burden with complex traits, particularly in case-control data. These conflicting results might suggest that the effects of autozygosity are confounded by various factors, such as socioeconomic status, education, urbanicity, and religiosity, which may be associated with both real inbreeding and the outcome measures of interest.

GeneEvolve: a fast and memory efficient forward-time simulator of realistic whole-genome sequence and SNP data.

MOTIVATION: Computer simulations are excellent tools for understanding the evolutionary and genetic consequences of complex processes that cannot be analytically predicted and for creating realistic genetic data. There are many software packages that simulate genetic data, but they are typically not fast or memory efficient enough to simulate realistic, individual-level genome-wide SNP/sequence data. RESULTS: GeneEvolve is a user-friendly and efficient population genetics simulator that handles complex evolutionary and life history scenarios and generates individual-level phenotypes and realistic whole-genome sequence or SNP data. GeneEvolve runs forward-in-time, which allows it to provide a wide range of scenarios for mating systems, selection, population size and structure, migration, recombination and environmental effects. The software is designed to use as input data from real or previously simulated phased haplotypes, allowing it to mimic very closely the properties of real genomic data. AVAILABILITY AND IMPLEMENTATION: GeneEvolve is freely available at https://github.com/rtahmasbi/GeneEvolve CONTACT: Rasool.Tahmasbi@Colorado.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Narrow-sense heritability estimation of complex traits using identity-by-descent information.

Heritability is a fundamental parameter in genetics. Traditional estimates based on family or twin studies can be biased due to shared environmental or non-additive genetic variance. Alternatively, those based on genotyped or imputed variants typically underestimate narrow-sense heritability contributed by rare or otherwise poorly tagged causal variants. Identical-by-descent (IBD) segments of the genome share all variants between pairs of chromosomes except new mutations that have arisen since the last common ancestor. Therefore, relating phenotypic similarity to degree of IBD sharing among classically unrelated individuals is an appealing approach to estimating the near full additive genetic variance while possibly avoiding biases that can occur when modeling close relatives. We applied an IBD-based approach (GREML-IBD) to estimate heritability in unrelated individuals using phenotypic simulation with thousands of whole-genome sequences across a range of stratification, polygenicity levels, and the minor allele frequencies of causal variants (CVs). In simulations, the IBD-based approach produced unbiased heritability estimates, even when CVs were extremely rare, although precision was low. However, population stratification and non-genetic familial environmental effects shared across generations led to strong biases in IBD-based heritability. We used data on two traits in ~120,000 people from the UK Biobank to demonstrate that, depending on the trait and possible confounding environmental effects, GREML-IBD can be applied to very large genetic datasets to infer the contribution of very rare variants lost using other methods. However, we observed apparent biases in these real data, suggesting that more work may be required to understand and mitigate factors that influence IBD-based heritability estimates.