Assortative mating biases marker-based heritability estimators

The genetic architecture of dog ownership: large-scale genome-wide association study in 97,552 European-ancestry individuals

Dog ownership has been associated with several complex traits, and there is evidence of genetic influence. We performed a genome-wide association study of dog ownership through a meta-analysis of 31,566 Swedish twins in 5 discovery cohorts and an additional 65,986 European-ancestry individuals in 3 replication cohorts from Sweden, Norway, and the United Kingdom. Association tests with >7.4 million single-nucleotide polymorphisms were meta-analyzed using a fixed effect model after controlling for population structure and relatedness. We identified 2 suggestive loci using discovery cohorts, which did not reach genome-wide significance after meta-analysis with replication cohorts. Single-nucleotide polymorphism-based heritability of dog ownership using linkage disequilibrium score regression was estimated at 0.123 (CI 0.038-0.207) using the discovery cohorts and 0.018 (CI -0.002 to 0.039) when adding in replication cohorts. Negative genetic correlation with complex traits including type 2 diabetes, depression, neuroticism, and asthma was only found using discovery summary data. Furthermore, we did not identify any genes/gene-sets reaching even a suggestive level of significance. This genome-wide association study does not, by itself, provide clear evidence on common genetic variants that influence dog ownership among European-ancestry individuals.

The impact of PM2.5 on lung function and chronic respiratory diseases: insights from genetic evidence

BACKGROUND

PM2.5 has been associated with various adverse health effects, particularly affecting lung function and chronic respiratory diseases. However, the genetic causality relationship between PM2.5 exposure and lung function as well as chronic respiratory diseases remains poorly understood.

METHOD

We conducted a two-sample Mendelian randomization analysis to investigate the causal impact of PM2.5 on lung function and chronic respiratory diseases. Instrumental variables were carefully selected, with significance thresholds (P < 5 × 10- 8), and linkage disequilibrium with an r2 value below 0.001. Additionally, SNPs with an F-statistic exceeding 10 were included to mitigate potential bias stemming from weak instrumental variables. The primary analytical approach employed the Inverse Variance Weighted method, supplemented by the Weighted Median, MR-Egger, Simple Model, and Weighted Model. Furthermore, pleiotropy and heterogeneity were evaluated through the MR-Egger intercept test and Cochrane's Q test, with a sensitivity analysis conducted using the leave-one-out method.

RESULTS

Eight SNPs significantly associated with PM2.5 exposure were identified as Instrumental variables. Mendelian randomization analysis revealed a significant causal association between PM2.5 exposure and lung function (FEV), with an OR of 0.7284 (95% CI: 0.5799-0.9150). Similarly, PM2.5 exposure demonstrated a substantial causal effect on asthma, with an OR of 1.5280 (95% CI: 1.0470-2.2299). However, no causal association was observed between PM2.5 exposure and chronic obstructive pulmonary disease, with an OR of 1.5176 (95% CI: 0.8294-2.7768).

CONCLUSION

These findings emphasize the necessity for continued research efforts in environmental health to develop effective strategies for the prevention and management of chronic respiratory diseases.

Phenotypic and ancestry-related assortative mating in autism

BACKGROUND

Positive assortative mating (AM) in several neuropsychiatric traits, including autism, has been noted. However, it is unknown whether the pattern of AM is different in phenotypically defined autism subgroups [e.g., autism with and without intellectually disability (ID)]. It is also unclear what proportion of the phenotypic AM can be explained by the genetic similarity between parents of children with an autism diagnosis, and the consequences of AM on the genetic structure of the population.

METHODS

To address these questions, we analyzed two family-based autism collections: the Simons Foundation Powering Autism Research for Knowledge (SPARK) (1575 families) and the Simons Simplex Collection (SSC) (2283 families).

RESULTS

We found a similar degree of phenotypic and ancestry-related AM in parents of children with an autism diagnosis regardless of the presence of ID. We did not find evidence of AM for autism based on autism polygenic scores (PGS) (at a threshold of |r|> 0.1). The adjustment of ancestry-related AM or autism PGS accounted for only 0.3-4% of the fractional change in the estimate of the phenotypic AM. The ancestry-related AM introduced higher long-range linkage disequilibrium (LD) between single nucleotide polymorphisms (SNPs) on different chromosomes that are highly ancestry-informative compared to SNPs that are less ancestry-informative (D2 on the order of 1 × 10-5).

LIMITATIONS

We only analyzed participants of European ancestry, limiting the generalizability of our results to individuals of non-European ancestry. SPARK and SSC were both multicenter studies. Therefore, there could be ancestry-related AM in SPARK and SSC due to geographic stratification. The study participants from each site were unknown, so we were unable to evaluate for geographic stratification.

CONCLUSIONS

This study showed similar patterns of AM in autism with and without ID, and demonstrated that the common genetic influences of autism are likely relevant to both autism groups. The adjustment of ancestry-related AM and autism PGS accounted for < 5% of the fractional change in the estimate of the phenotypic AM. Future studies are needed to evaluate if the small increase of long-range LD induced by ancestry-related AM has impact on the downstream analysis.

Genetic similarity between relatives provides evidence on the presence and history of assortative mating

Assortative mating - the non-random mating of individuals with similar traits - is known to increase trait-specific genetic variance and genetic similarity between relatives. However, empirical evidence is limited for many traits, and the implications hinge on whether assortative mating has started recently or many generations ago. Here we show theoretically and empirically that genetic similarity between relatives can provide evidence on the presence and history of assortative mating. First, we employed path analysis to understand how assortative mating affects genetic similarity between family members across generations, finding that similarity between distant relatives is more affected than close relatives. Next, we correlated polygenic indices of 47,135 co-parents from the Norwegian Mother, Father, and Child Cohort Study (MoBa) and found genetic evidence of assortative mating in nine out of sixteen examined traits. The same traits showed elevated similarity between relatives, especially distant relatives. Six of the nine traits, including educational attainment, showed greater genetic variance among offspring, which is inconsistent with stable assortative mating over many generations. These results suggest an ongoing increase in familial similarity for these traits. The implications of this research extend to genetic methodology and the understanding of social and economic disparities.

Investigating the causal associations between five anthropometric indicators and nonalcoholic fatty liver disease: Mendelian randomization study

BACKGROUND

Although the etiology of nonalcoholic fatty liver disease (NAFLD) has not been thoroughly understood, the emerging roles of anthropometric indicators in assessing and predicting the risk of NAFLD have been highlighted by accumulating evidence.

AIM

To evaluate the causal relationships between five anthropometric indicators and NAFLD employing Mendelian randomization (MR) design.

METHODS

The Anthropometric Consortium provided genetic exposure data for five anthropometric indicators, including hip circumference (HC), waist circumference (WC), waist-to-hip ratio (WHR), body mass index (BMI), and body fat percentage (BF). Genetic outcome data for NAFLD were obtained from the United Kingdom Biobank and FinnGen Consortium. Genome-wide significant single nucleotide polymorphisms were chosen as instrumental variables. Univariable MR (UVMR) and multivariable MR (MVMR) designs with analytical approaches, including inverse variance weighted (IVW), MR-Egger, weighted median (WM), and weighted mode methods, were used to assess the causal relationships between anthropometric indicators and NAFLD.

RESULTS

Causal relationships were revealed by UVMR, indicating that a higher risk of NAFLD was associated with a per-unit increase in WC [IVW: odds ratio (OR) = 2.67, 95%CI: 1.42-5.02, P = 2.25 × 10-3], and BF was causally associated with an increased risk of NAFLD (WM: OR = 2.23, 95%CI: 1.07-4.66, P = 0.033). The presence of causal effects of WC on the decreased risk of NAFLD was supported by MVMR after adjusting for BMI and smoking. However, no causal association between BF and NAFLD was observed. In addition, other causal relationships of HC, WHR (BMI adjusted), and BMI with the risk of NAFLD were not retained after FDR correction.

CONCLUSION

This study establishes a causal relationship, indicating that an increase in WC is associated with a higher risk of NAFLD. This demonstrates that a suitable decrease in WC is advantageous for preventing NAFLD.

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

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

An evolutionary perspective on complex neuropsychiatric disease

The forces of evolution-mutation, selection, migration, and genetic drift-shape the genetic architecture of human traits, including the genetic architecture of complex neuropsychiatric illnesses. Studying these illnesses in populations that are diverse in genetic ancestry, historical demography, and cultural history can reveal how evolutionary forces have guided adaptation over time and place. A fundamental truth of shared human biology is that an allele responsible for a disease in anyone, anywhere, reveals a gene critical to the normal biology underlying that condition in everyone, everywhere. Understanding the genetic causes of neuropsychiatric disease in the widest possible range of human populations thus yields the greatest possible range of insight into genes critical to human brain development. In this perspective, we explore some of the relationships between genes, adaptation, and history that can be illuminated by an evolutionary perspective on studies of complex neuropsychiatric disease in diverse populations.

Causal Relationships between Air Pollutant Exposure and Bone Mineral Density and the Risk of Bone Fractures: Evidence from a Two-Stage Mendelian Randomization Analysis

A number of studies from the literature have suggested that exposure to air pollutants is associated with a declined bone mineral density (BMD), and increased risks of osteoporosis (OP) and bone fractures. This study was performed to systemically assess the genetically causal associations of air pollutants with site-/age-specific BMD and risk of bone fractures with the implementation of two-sample Mendelian randomization (TSMR) and multivariate Mendelian randomization (MVMR). The TSMR analysis was implemented to infer the causal associations between air pollutants and BMD and the risk of bone fractures, additional MVMR analysis was used to further estimate the direct causal effects between air pollutants and BMD, the occurrence of OP, and bone fractures. The results showed that NOx exposure contributed to lower femoral neck BMD (FN-BMD) (β = -0.71, 95%CI: -1.22, -0.20, p = 0.006) and total body BMD (TB-BMD) (β = -0.55, 95%CI: -0.90, -0.21, p = 0.002). Additionally, exposure to PM10 was found to be associated with a decreased TB-BMD (B β = -0.42, 95%CI: -0.66, -0.18, p = 0.001), further age-specific subgroup analysis demonstrated the causal effect of PM10 exposure on the decreased TB-BMD in a subgroup aged 45 to 60 years (β = -0.70, 95%CI: -1.12, -0.29, p = 0.001). Moreover, the findings of the MVMR analysis implied that there was a direct causal effect between PM10 exposure and the decreased TB-BMD (45 < age < 60), after adjusting for PM2.5 and PM2.5 -10 exposure. Our study provides additional evidence to support the causal associations of higher concentrations of air pollutant exposure with decreased BMD, especially in those populations aged between 45 to 60 years, suggesting that early intervention measures and public policy should be considered to improve public health awareness and promote bone health.

Inferring disease architecture and predictive ability with LDpred2-auto

LDpred2 is a widely used Bayesian method for building polygenic scores (PGSs). LDpred2-auto can infer the two parameters from the LDpred model, the SNP heritability h2 and polygenicity p, so that it does not require an additional validation dataset to choose best-performing parameters. The main aim of this paper is to properly validate the use of LDpred2-auto for inferring multiple genetic parameters. Here, we present a new version of LDpred2-auto that adds an optional third parameter α to its model, for modeling negative selection. We then validate the inference of these three parameters (or two, when using the previous model). We also show that LDpred2-auto provides per-variant probabilities of being causal that are well calibrated and can therefore be used for fine-mapping purposes. We also introduce a formula to infer the out-of-sample predictive performance r2 of the resulting PGS directly from the Gibbs sampler of LDpred2-auto. Finally, we extend the set of HapMap3 variants recommended to use with LDpred2 with 37% more variants to improve the coverage of this set, and we show that this new set of variants captures 12% more heritability and provides 6% more predictive performance, on average, in UK Biobank analyses.

Causal interpretations of family GWAS in the presence of heterogeneous effects

Family-based genome-wide association studies (GWAS) have emerged as a gold standard for assessing causal effects of alleles and polygenic scores. Notably, family studies are often claimed to provide an unbiased estimate of the average causal effect (or average treatment effect; ATE) of an allele, on the basis of an analogy between the random transmission of alleles from parents to children and a randomized controlled trial. Here, we show that this interpretation does not hold in general. Because Mendelian segregation only randomizes alleles among children of heterozygotes, the effects of alleles in the children of homozygotes are not observable. Consequently, if an allele has different average effects in the children of homozygotes and heterozygotes, as can arise in the presence of gene-by-environment interactions, gene-by-gene interactions, or differences in LD patterns, family studies provide a biased estimate of the average effect in the sample. At a single locus, family-based association studies can be thought of as providing an unbiased estimate of the average effect in the children of heterozygotes (i.e., a local average treatment effect; LATE). This interpretation does not extend to polygenic scores, however, because different sets of SNPs are heterozygous in each family. Therefore, other than under specific conditions, the within-family regression slope of a PGS cannot be assumed to provide an unbiased estimate for any subset or weighted average of families. Instead, family-based studies can be reinterpreted as enabling an unbiased estimate of the extent to which Mendelian segregation at loci in the PGS contributes to the population-level variance in the trait. Because this estimate does not include the between-family variance, however, this interpretation applies to only (roughly) half of the sample PGS variance. In practice, the potential biases of a family-based GWAS are likely smaller than those arising from confounding in a standard, population-based GWAS, and so family studies remain important for the dissection of genetic contributions to phenotypic variation. Nonetheless, the causal interpretation of family-based GWAS estimates is less straightforward than has been widely appreciated.

Revealing polygenic pleiotropy using genetic risk scores for asthma

In this study we examined how genetic risk for asthma associates with different features of the disease and with other medical conditions and traits. Using summary statistics from two multi-ancestry genome-wide association studies of asthma, we modeled polygenic risk scores (PRSs) and validated their predictive performance in the UK Biobank. We then performed phenome-wide association studies of the asthma PRSs with 371 heritable traits in the UK Biobank. We identified 228 total significant associations across a variety of organ systems, including associations that varied by PRS model, sex, age of asthma onset, ancestry, and human leukocyte antigen region alleles. Our results highlight pervasive pleiotropy between asthma and numerous other traits and conditions and elucidate pathways that contribute to asthma and its comorbidities.

A General Approach to Adjusting Genetic Studies for Assortative Mating

The effects of assortative mating (AM) on estimates from genetic studies has been receiving increasing attention in recent years. We extend existing AM theory to more general models of sorting and conclude that correct theory-based AM adjustments require knowledge of complicated, unknown historical sorting patterns. We propose a simple, general-purpose approach using polygenic indexes (PGIs). Our approach can estimate the fraction of genetic variance and genetic correlation that is driven by AM. Our approach is less effective when applied to Mendelian randomization (MR) studies for two reasons: AM can induce a form of selection bias in MR studies that remains after our adjustment; and, in the MR context, the adjustment is particularly sensitive to PGI estimation error. Using data from the UK Biobank, we find that AM inflates genetic correlation estimates between health traits and education by 14% on average. Our results suggest caution in interpreting genetic correlations or MR estimates for traits subject to AM.

Evidence of correlations between human partners based on systematic reviews and meta-analyses of 22 traits and UK Biobank analysis of 133 traits

Positive correlations between mates can increase trait variation and prevalence, as well as bias estimates from genetically informed study designs. While past studies of similarity between human mating partners have largely found evidence of positive correlations, to our knowledge, no formal meta-analysis has examined human partner correlations across multiple categories of traits. Thus, we conducted systematic reviews and random-effects meta-analyses of human male-female partner correlations across 22 traits commonly studied by psychologists, economists, sociologists, anthropologists, epidemiologists and geneticists. Using ScienceDirect, PubMed and Google Scholar, we incorporated 480 partner correlations from 199 peer-reviewed studies of co-parents, engaged pairs, married pairs and/or cohabitating pairs that were published on or before 16 August 2022. We also calculated 133 trait correlations using up to 79,074 male-female couples in the UK Biobank (UKB). Estimates of the 22 mean meta-analysed correlations ranged from rmeta = 0.08 (adjusted 95% CI = 0.03, 0.13) for extraversion to rmeta = 0.58 (adjusted 95% CI = 0.50, 0.64) for political values, with funnel plots showing little evidence of publication bias across traits. The 133 UKB correlations ranged from rUKB = -0.18 (adjusted 95% CI = -0.20, -0.16) for chronotype (being a 'morning' or 'evening' person) to rUKB = 0.87 (adjusted 95% CI = 0.86, 0.87) for birth year. Across analyses, political and religious attitudes, educational attainment and some substance use traits showed the highest correlations, while psychological (that is, psychiatric/personality) and anthropometric traits generally yielded lower but positive correlations. We observed high levels of between-sample heterogeneity for most meta-analysed traits, probably because of both systematic differences between samples and true differences in partner correlations across populations.

Extremely sparse models of linkage disequilibrium in ancestrally diverse association studies

Linkage disequilibrium (LD) is the correlation among nearby genetic variants. In genetic association studies, LD is often modeled using large correlation matrices, but this approach is inefficient, especially in ancestrally diverse studies. In the present study, we introduce LD graphical models (LDGMs), which are an extremely sparse and efficient representation of LD. LDGMs are derived from genome-wide genealogies; statistical relationships among alleles in the LDGM correspond to genealogical relationships among haplotypes. We published LDGMs and ancestry-specific LDGM precision matrices for 18 million common variants (minor allele frequency >1%) in five ancestry groups, validated their accuracy and demonstrated order-of-magnitude improvements in runtime for commonly used LD matrix computations. We implemented an extremely fast multiancestry polygenic prediction method, BLUPx-ldgm, which performs better than a similar method based on the reference LD correlation matrix. LDGMs will enable sophisticated methods that scale to ancestrally diverse genetic association data across millions of variants and individuals.

rBahadur: efficient simulation of structured high-dimensional genotype data with applications to assortative mating

Existing methods for generating synthetic genotype data are ill-suited for replicating the effects of assortative mating (AM). We propose rb_dplr, a novel and computationally efficient algorithm for generating high-dimensional binary random variates that effectively recapitulates AM-induced genetic architectures using the Bahadur order-2 approximation of the multivariate Bernoulli distribution. The rBahadur R library is available through the Comprehensive R Archive Network at https://CRAN.R-project.org/package=rBahadur .

Overcoming attenuation bias in regressions using polygenic indices

Measurement error in polygenic indices (PGIs) attenuates the estimation of their effects in regression models. We analyze and compare two approaches addressing this attenuation bias: Obviously Related Instrumental Variables (ORIV) and the PGI Repository Correction (PGI-RC). Through simulations, we show that the PGI-RC performs slightly better than ORIV, unless the prediction sample is very small (N < 1000) or when there is considerable assortative mating. Within families, ORIV is the best choice since the PGI-RC correction factor is generally not available. We verify the empirical validity of the simulations by predicting educational attainment and height in a sample of siblings from the UK Biobank. We show that applying ORIV between families increases the standardized effect of the PGI by 12% (height) and by 22% (educational attainment) compared to a meta-analysis-based PGI, yet estimates remain slightly below the PGI-RC estimates. Furthermore, within-family ORIV regression provides the tightest lower bound for the direct genetic effect, increasing the lower bound for the standardized direct genetic effect on educational attainment from 0.14 to 0.18 (+29%), and for height from 0.54 to 0.61 (+13%) compared to a meta-analysis-based PGI.

Participation bias in the UK Biobank distorts genetic associations and downstream analyses

While volunteer-based studies such as the UK Biobank have become the cornerstone of genetic epidemiology, the participating individuals are rarely representative of their target population. To evaluate the impact of selective participation, here we derived UK Biobank participation probabilities on the basis of 14 variables harmonized across the UK Biobank and a representative sample. We then conducted weighted genome-wide association analyses on 19 traits. Comparing the output from weighted genome-wide association analyses (neffective = 94,643 to 102,215) with that from standard genome-wide association analyses (n = 263,464 to 283,749), we found that increasing representativeness led to changes in SNP effect sizes and identified novel SNP associations for 12 traits. While heritability estimates were less impacted by weighting (maximum change in h2, 5%), we found substantial discrepancies for genetic correlations (maximum change in rg, 0.31) and Mendelian randomization estimates (maximum change in βSTD, 0.15) for socio-behavioural traits. We urge the field to increase representativeness in biobank samples, especially when studying genetic correlates of behaviour, lifestyles and social outcomes.

The genetic basis of major depressive disorder

The genetic dissection of major depressive disorder (MDD) ranks as one of the success stories of psychiatric genetics, with genome-wide association studies (GWAS) identifying 178 genetic risk loci and proposing more than 200 candidate genes. However, the GWAS results derive from the analysis of cohorts in which most cases are diagnosed by minimal phenotyping, a method that has low specificity. I review data indicating that there is a large genetic component unique to MDD that remains inaccessible to minimal phenotyping strategies and that the majority of genetic risk loci identified with minimal phenotyping approaches are unlikely to be MDD risk loci. I show that inventive uses of biobank data, novel imputation methods, combined with more interviewer diagnosed cases, can identify loci that contribute to the episodic severe shifts of mood, and neurovegetative and cognitive changes that are central to MDD. Furthermore, new theories about the nature and causes of MDD, drawing upon advances in neuroscience and psychology, can provide handles on how best to interpret and exploit genetic mapping results.

Transcriptome-wide gene-gene interaction associations elucidate pathways and functional enrichment of complex traits

It remains unknown to what extent gene-gene interactions contribute to complex traits. Here, we introduce a new approach using predicted gene expression to perform exhaustive transcriptome-wide interaction studies (TWISs) for multiple traits across all pairs of genes expressed in several tissue types. Using imputed transcriptomes, we simultaneously reduce the computational challenge and improve interpretability and statistical power. We discover (in the UK Biobank) and replicate (in independent cohorts) several interaction associations, and find several hub genes with numerous interactions. We also demonstrate that TWIS can identify novel associated genes because genes with many or strong interactions have smaller single-locus model effect sizes. Finally, we develop a method to test gene set enrichment of TWIS associations (E-TWIS), finding numerous pathways and networks enriched in interaction associations. Epistasis is may be widespread, and our procedure represents a tractable framework for beginning to explore gene interactions and identify novel genomic targets.

Partner choice, confounding and trait convergence all contribute to phenotypic partner similarity

Partners are often similar in terms of their physical and behavioural traits, such as their education, political affiliation and height. However, it is currently unclear what exactly causes this similarity-partner choice, partner influence increasing similarity over time or confounding factors such as shared environment or indirect assortment. Here, we applied Mendelian randomization to the data of 51,664 couples in the UK Biobank and investigated partner similarity in 118 traits. We found evidence of partner choice for 64 traits, 40 of which had larger phenotypic correlation than causal effect. This suggests that confounders contribute to trait similarity, among which household income, overall health rating and education accounted for 29.8, 14.1 and 11.6% of correlations between partners, respectively. Finally, mediation analysis revealed that most causal associations between different traits in the two partners are indirect. In summary, our results show the mechanisms through which indirect assortment increases the observed partner similarity.

Celebrating a Century of Research in Behavioral Genetics

A century after the first twin and adoption studies of behavior in the 1920s, this review looks back on the journey and celebrates milestones in behavioral genetic research. After a whistle-stop tour of early quantitative genetic research and the parallel journey of molecular genetics, the travelogue focuses on the last fifty years. Just as quantitative genetic discoveries were beginning to slow down in the 1990s, molecular genetics made it possible to assess DNA variation directly. From a rocky start with candidate gene association research, by 2005 the technological advance of DNA microarrays enabled genome-wide association studies, which have successfully identified some of the DNA variants that contribute to the ubiquitous heritability of behavioral traits. The ability to aggregate the effects of thousands of DNA variants in polygenic scores has created a DNA revolution in the behavioral sciences by making it possible to use DNA to predict individual differences in behavior from early in life.

Human height: a model common complex trait

CONTEXT

Like other complex phenotypes, human height reflects a combination of environmental and genetic factors, but is notable for being exceptionally easy to measure. Height has therefore been commonly used to make observations later generalised to other phenotypes though the appropriateness of such generalisations is not always considered.

OBJECTIVES

We aimed to assess height's suitability as a model for other complex phenotypes and review recent advances in height genetics with regard to their implications for complex phenotypes more broadly.

METHODS

We conducted a comprehensive literature search in PubMed and Google Scholar for articles relevant to the genetics of height and its comparatibility to other phenotypes.

RESULTS

Height is broadly similar to other phenotypes apart from its high heritability and ease of measurment. Recent genome-wide association studies (GWAS) have identified over 12,000 independent signals associated with height and saturated height's common single nucleotide polymorphism based heritability of height within a subset of the genome in individuals similar to European reference populations.

CONCLUSIONS

Given the similarity of height to other complex traits, the saturation of GWAS's ability to discover additional height-associated variants signals potential limitations to the omnigenic model of complex-phenotype inheritance, indicating the likely future power of polygenic scores and risk scores, and highlights the increasing need for large-scale variant-to-gene mapping efforts.

The genetic architecture of the human skeletal form

The human skeletal form underlies our ability to walk on two legs, but unlike standing height, the genetic basis of limb lengths and skeletal proportions is less well understood. Here we applied a deep learning model to 31,221 whole body dual-energy X-ray absorptiometry (DXA) images from the UK Biobank (UKB) to extract 23 different image-derived phenotypes (IDPs) that include all long bone lengths as well as hip and shoulder width, which we analyzed while controlling for height. All skeletal proportions are highly heritable (∼40-50%), and genome-wide association studies (GWAS) of these traits identified 179 independent loci, of which 102 loci were not associated with height. These loci are enriched in genes regulating skeletal development as well as associated with rare human skeletal diseases and abnormal mouse skeletal phenotypes. Genetic correlation and genomic structural equation modeling indicated that limb proportions exhibited strong genetic sharing but were genetically independent of width and torso proportions. Phenotypic and polygenic risk score analyses identified specific associations between osteoarthritis (OA) of the hip and knee, the leading causes of adult disability in the United States, and skeletal proportions of the corresponding regions. We also found genomic evidence of evolutionary change in arm-to-leg and hip-width proportions in humans consistent with striking anatomical changes in these skeletal proportions in the hominin fossil record. In contrast to cardiovascular, auto-immune, metabolic, and other categories of traits, loci associated with these skeletal proportions are significantly enriched in human accelerated regions (HARs), and regulatory elements of genes differentially expressed through development between humans and the great apes. Taken together, our work validates the use of deep learning models on DXA images to identify novel and specific genetic variants affecting the human skeletal form and ties a major evolutionary facet of human anatomical change to pathogenesis.

Decomposing heritability and genetic covariance by direct and indirect effect paths

Estimation of heritability and genetic covariance is crucial for quantifying and understanding complex trait genetic architecture and is employed in almost all recent genome-wide association studies (GWAS). However, many existing approaches for heritability estimation and almost all methods for estimating genetic correlation ignore the presence of indirect genetic effects, i.e., genotype-phenotype associations confounded by the parental genome and family environment, and may thus lead to incorrect interpretation especially for human sociobehavioral phenotypes. In this work, we introduce a statistical framework to decompose heritability and genetic covariance into multiple components representing direct and indirect effect paths. Applied to five traits in UK Biobank, we found substantial involvement of indirect genetic components in shared genetic architecture across traits. These results demonstrate the effectiveness of our approach and highlight the importance of accounting for indirect effects in variance component analysis of complex traits.

Cross-trait assortative mating is widespread and inflates genetic correlation estimates

The observation of genetic correlations between disparate human traits has been interpreted as evidence of widespread pleiotropy. Here, we introduce cross-trait assortative mating (xAM) as an alternative explanation. We observe that xAM affects many phenotypes and that phenotypic cross-mate correlation estimates are strongly associated with genetic correlation estimates (R²=74%). We demonstrate that existing xAM plausibly accounts for substantial fractions of genetic correlation estimates and that previously reported genetic correlation estimates between some pairs of psychiatric disorders are congruent with xAM alone. Finally, we provide evidence for a history of xAM at the genetic level using cross-trait even/odd chromosome polygenic score correlations. Together, our results demonstrate that previous reports have likely overestimated the true genetic similarity between many phenotypes.

Genetic assortative mating for schizophrenia and bipolar disorder

BACKGROUND

Psychiatric disorders are highly polygenic and show patterns of partner resemblance. Partner resemblance has direct population-level genetic implications if it is caused by assortative mating, but not if it is caused by convergence or social homogamy. Using genetics may help distinguish these different mechanisms. Here, we investigated whether partner resemblance for schizophrenia and bipolar disorder is influenced by assortative mating using polygenic risk scores (PRSs).

METHODS

PRSs from The Danish High-Risk and Resilience Study-VIA 7 were compared between parents in three subsamples: population-based control parent pairs (N=198), parent pairs where at least one parent had schizophrenia (N=193), and parent pairs where at least one parent had bipolar disorder (N=115).

RESULTS

The PRS for schizophrenia was predictive of schizophrenia in the full sample and showed a significant correlation between parent pairs (r=0.121, p=0.0440), indicative of assortative mating. The PRS for bipolar disorder was also correlated between parent pairs (r=0.162, p=0.0067), but it was not predictive of bipolar disorder in the full sample, limiting the interpretation.

CONCLUSIONS

Our study provides genetic evidence for assortative mating for schizophrenia, with important implications for our understanding of the genetics of schizophrenia.

Modeling assortative mating and genetic similarities between partners, siblings, and in-laws

Assortative mating on heritable traits can have implications for the genetic resemblance between siblings and in-laws in succeeding generations. We studied polygenic scores and phenotypic data from pairs of partners (n = 26,681), siblings (n = 2,170), siblings-in-law (n = 3,905), and co-siblings-in-law (n = 1,763) in the Norwegian Mother, Father and Child Cohort Study. Using structural equation models, we estimated associations between measurement error-free latent genetic and phenotypic variables. We found evidence of genetic similarity between partners for educational attainment (r_g = 0.37), height (r_g = 0.13), and depression (r_g = 0.08). Common genetic variants associated with educational attainment correlated between siblings above 0.50 (r_g = 0.68) and between siblings-in-law (r_g = 0.25) and co-siblings-in-law (r_g = 0.09). Indirect assortment on secondary traits accounted for partner similarity in education and depression, but not in height. Comparisons between the genetic similarities of partners and siblings indicated that genetic variances were in intergenerational equilibrium. This study shows genetic similarities between extended family members and that assortative mating has taken place for several generations.