GWAS of three molecular traits highlights core genes and pathways alongside a highly polygenic background

A scalable and robust variance components method reveals insights into the architecture of gene-environment interactions underlying complex traits

Understanding the contribution of gene-environment interactions (GxE) to complex trait variation can provide insights into disease mechanisms, explain sources of heritability, and improve genetic risk prediction. While large biobanks with genetic and deep phenotypic data hold promise for obtaining novel insights into GxE, our understanding of GxE architecture in complex traits remains limited. We introduce a method to estimate the proportion of trait variance explained by GxE (GxE heritability) and additive genetic effects (additive heritability) across the genome and within specific genomic annotations. We show that our method is accurate in simulations and computationally efficient for biobank-scale datasets. We applied our method to common array SNPs (MAF ≥1%), fifty quantitative traits, and four environmental variables (smoking, sex, age, and statin usage) in unrelated white British individuals in the UK Biobank. We found 68 trait-E pairs with significant genome-wide GxE heritability (p<0.05/200) with a ratio of GxE to additive heritability of ≈6.8% on average. Analyzing ≈8 million imputed SNPs (MAF ≥0.1%), we documented an approximate 28% increase in genome-wide GxE heritability compared to array SNPs. We partitioned GxE heritability across minor allele frequency (MAF) and local linkage disequilibrium (LD) values, revealing that, like additive allelic effects, GxE allelic effects tend to increase with decreasing MAF and LD. Analyzing GxE heritability near genes highly expressed in specific tissues, we find significant brain-specific enrichment for body mass index (BMI) and basal metabolic rate in the context of smoking and adipose-specific enrichment for waist-hip ratio (WHR) in the context of sex.

The pathogenesis of gout: molecular insights from genetic, epigenomic and transcriptomic studies

The pathogenesis of gout involves a series of steps beginning with hyperuricaemia, followed by the deposition of monosodium urate crystal in articular structures and culminating in an innate immune response, mediated by the NLRP3 inflammasome, to the deposited crystals. Large genome-wide association studies (GWAS) of serum urate levels initially identified the genetic variants with the strongest effects, mapping mainly to genes that encode urate transporters in the kidney and gut. Other GWAS highlighted the importance of uncommon genetic variants. More recently, genetic and epigenetic genome-wide studies have revealed new pathways in the inflammatory process of gout, including genetic associations with epigenomic modifiers. Epigenome-wide association studies are also implicating epigenomic remodelling in gout, which perhaps regulates the responsiveness of the innate immune system to monosodium urate crystals. Notably, genes implicated in gout GWAS do not include those encoding components of the NLRP3 inflammasome itself, but instead include genes encoding molecules involved in its regulation. Knowledge of the molecular mechanisms underlying gout has advanced through the translation of genetic associations into specific molecular mechanisms. Notable examples include ABCG2, HNF4A, PDZK1, MAF and IL37. Current genetic studies are dominated by participants of European ancestry; however, studies focusing on other population groups are discovering informative population-specific variants associated with gout.

Discovering mechanisms of human genetic variation and controlling cell states at scale

Population-scale sequencing efforts have catalogued substantial genetic variation in humans such that variant discovery dramatically outpaces interpretation. We discuss how single-cell sequencing is poised to reveal genetic mechanisms at a rate that may soon approach that of variant discovery. The functional genomics toolkit is sufficiently modular to systematically profile almost any type of variation within increasingly diverse contexts and with molecularly comprehensive and unbiased readouts. As a result, we can construct deep phenotypic atlases of variant effects that span the entire regulatory cascade. The same conceptual approach to interpreting genetic variation should be applied to engineering therapeutic cell states. In this way, variant mechanism discovery and cell state engineering will become reciprocating and iterative processes towards genomic medicine.

Multivariate genetic architecture reveals testosterone-driven sexual antagonism in contemporary humans

Sex difference (SD) is ubiquitous in humans despite shared genetic architecture (SGA) between the sexes. A univariate approach, i.e., studying SD in single traits by estimating genetic correlation, does not provide a complete biological overview, because traits are not independent and are genetically correlated. The multivariate genetic architecture between the sexes can be summarized by estimating the additive genetic (co)variance across shared traits, which, apart from the cross-trait and cross-sex covariances, also includes the cross-sex-cross-trait covariances, e.g., between height in males and weight in females. Using such a multivariate approach, we investigated SD in the genetic architecture of 12 anthropometric, fat depositional, and sex-hormonal phenotypes. We uncovered sexual antagonism (SA) in the cross-sex-cross-trait covariances in humans, most prominently between testosterone and the anthropometric traits - a trend similar to phenotypic correlations. 27% of such cross-sex-cross-trait covariances were of opposite sign, contributing to asymmetry in the SGA. Intriguingly, using multivariate evolutionary simulations, we observed that the SGA acts as a genetic constraint to the evolution of SD in humans only when selection is sexually antagonistic and not concordant. Remarkably, we found that the lifetime reproductive success in both the sexes shows a positive genetic correlation with anthropometric traits, but not with testosterone. Moreover, we demonstrated that genetic variance is depleted along multivariate trait combinations in both the sexes but in different directions, suggesting absolute genetic constraint to evolution. Our results indicate that testosterone drives SA in contemporary humans and emphasize the necessity and significance of using a multivariate framework in studying SD.

Genetic influences on testosterone and PTSD

Females are twice as likely to experience PTSD as compared to males. Although sex differences in prevalence are well-established, little is known about why such sex differences occur. Biological factors that vary with sex, including sex hormone production, may contribute to these differences. Considerable evidence links sex hormones, such as testosterone, to PTSD risk though less is known about the shared genetic underpinnings. The objective of the present study was to test for genetic relationships between testosterone and PTSD. To do so, we used summary statistics from large, publicly available genetic consortia to conduct linkage disequilibrium score regression to estimate the genetic correlations between PTSD and testosterone in males and females, and two-sample, bi-directional Mendelian randomization to examine potential causal relationships of testosterone on PTSD and the reverse. Heritability estimates of testosterone were significantly higher in males (0.17, SE = 0.02) than females (0.11, SE = 0.01; z = 2.46, p = 00.01). The correlation between testosterone and PTSD was negative in males (rg = -0.11, SE = 0.02, p = 6.7 x 10-6), but not significant in females (rg = 0.002, SE = 0.03, p = 0.95). MR analyses found no evidence of a causal effect of testosterone on PTSD or the reverse. Findings are consistent with phenotypic literature suggesting a relationship between testosterone and PTSD that may be sex-specific. This work provides early evidence of a relationship between testosterone and PTSD genotypically and suggests an avenue for future research that will enable a better understanding of disparities in PTSD.

Tradeoffs in Modeling Context Dependency in Complex Trait Genetics

Genetic effects on complex traits may depend on context, such as age, sex, environmental exposures or social settings. However, it is often unclear if the extent of context dependency, or Gene-by-Environment interaction (GxE), merits more involved models than the additive model typically used to analyze data from genome-wide association studies (GWAS). Here, we suggest considering the utility of GxE models in GWAS as a tradeoff between bias and variance parameters. In particular, We derive a decision rule for choosing between competing models for the estimation of allelic effects. The rule weighs the increased estimation noise when context is considered against the potential bias when context dependency is ignored. In the empirical example of GxSex in human physiology, the increased noise of context-specific estimation often outweighs the bias reduction, rendering GxE models less useful when variants are considered independently. However, we argue that for complex traits, the joint consideration of context dependency across many variants mitigates both noise and bias. As a result, polygenic GxE models can improve both estimation and trait prediction. Finally, we exemplify (using GxDiet effects on longevity in fruit flies) how analyses based on independently ascertained "top hits" alone can be misleading, and that considering polygenic patterns of GxE can improve interpretation.

Fathers and sons, mothers and daughters: Sex-specific genetic architecture for fetal testosterone in a wild mammal

Testosterone plays a critical role in mediating fitness-related traits in many species. Although it is highly responsive to environmental and social conditions, evidence from several species show a heritable component to its individual variation. Despite the known effects that in utero testosterone exposure have on adult fitness, the heritable component of individual testosterone variation in fetuses is mostly unexplored. Furthermore, testosterone has sex-differential effects on fetal development, i.e., a specific level may be beneficial for male fetuses but detrimental for females, producing sexual conflict. Such sexual conflict may be resolved by the evolution of a sex-specific genetic architecture of the trait. Here, we quantified fetal testosterone levels in a wild species, free-ranging nutrias (Myocastor coypus) using hair-testing and estimated testosterone heritability between parent and offspring from the same and opposite sex. We found that in utero accumulated hair testosterone levels were heritable between parents and offspring of the same sex. Moreover, there was a low additive genetic covariance between the sexes, and a low cross-sex genetic correlation, suggesting a potential for sex-specific trait evolution, expressed early on, in utero.

Obesity-Dependent Association of the rs10454142 PPP1R21 with Breast Cancer

The purpose of this work was to find a link between the breast cancer (BC)-risk effects of sex hormone-binding globulin (SHBG)-associated polymorphisms and obesity. The study was conducted on a sample of 1498 women (358 BC; 1140 controls) who, depending on the presence/absence of obesity, were divided into two groups: obese (119 BC; 253 controls) and non-obese (239 BC; 887 controls). Genotyping of nine SHBG-associated single nucleotide polymorphisms (SNP)-rs17496332 PRMT6, rs780093 GCKR, rs10454142 PPP1R21, rs3779195 BAIAP2L1, rs440837 ZBTB10, rs7910927 JMJD1C, rs4149056 SLCO1B1, rs8023580 NR2F2, and rs12150660 SHBG-was executed, and the BC-risk impact of these loci was analyzed by logistic regression separately in each group of obese/non-obese women. We found that the BC-risk effect correlated by GWAS with the SHBG-level polymorphism rs10454142 PPP1R21 depends on the presence/absence of obesity. The SHBG-lowering allele C rs10454142 PPP1R21 has a risk value for BC in obese women (allelic model: CvsT, OR = 1.52, 95%CI = 1.10-2.11, and pperm = 0.013; additive model: CCvsTCvsTT, OR = 1.71, 95%CI = 1.15-2.62, and pperm = 0.011; dominant model: CC + TCvsTT, OR = 1.95, 95%CI = 1.13-3.37, and pperm = 0.017) and is not associated with the disease in women without obesity. SNP rs10454142 PPP1R21 and 10 proxy SNPs have adipose-specific regulatory effects (epigenetic modifications of promoters/enhancers, DNA interaction with 51 transcription factors, eQTL/sQTL effects on five genes (PPP1R21, RP11-460M2.1, GTF2A1L, STON1-GTF2A1L, and STON1), etc.), can be "likely cancer driver" SNPs, and are involved in cancer-significant pathways. In conclusion, our study detected an obesity-dependent association of the rs10454142 PPP1R21 with BC in women.

Better together against genetic heterogeneity: A sex-combined joint main and interaction analysis of 290 quantitative traits in the UK Biobank

Genetic effects can be sex-specific, particularly for traits such as testosterone, a sex hormone. While sex-stratified analysis provides easily interpretable sex-specific effect size estimates, the presence of sex-differences in SNP effect implies a SNP×sex interaction. This suggests the usage of the often overlooked joint test, testing for an SNP's main and SNP×sex interaction effects simultaneously. Notably, even without individual-level data, the joint test statistic can be derived from sex-stratified summary statistics through an omnibus meta-analysis. Utilizing the available sex-stratified summary statistics of the UK Biobank, we performed such omnibus meta-analyses for 290 quantitative traits. Results revealed that this approach is robust to genetic effect heterogeneity and can outperform the traditional sex-stratified or sex-combined main effect-only tests. Therefore, we advocate using the omnibus meta-analysis that captures both the main and interaction effects. Subsequent sex-stratified analysis should be conducted for sex-specific effect size estimation and interpretation.

Genetic correlations of alcohol consumption and alcohol use disorder with sex hormone levels in females and males

Background

Alcohol consumption behaviors and alcohol use disorder risk and presentation differ by sex, and these complex traits are associated with blood concentrations of the steroid sex hormones, testosterone and estradiol, and their regulatory binding proteins, sex hormone binding globulin (SHBG) and albumin. Genetic variation is associated with alcohol consumption and alcohol use disorder, as well as levels of steroid sex hormones and their binding proteins.

Methods

To assess the contribution of genetic factors to previously described phenotypic associations between alcohol-use traits and sex-hormone levels, we estimated genetic correlations (rg) using summary statistics from prior published, large sample size genome-wide association studies (GWAS) of alcohol consumption, alcohol dependence, testosterone, estradiol, SHBG, and albumin.

Results

For alcohol consumption, we observed positive genetic correlation (i.e. genetic effects in the same direction) with total testosterone in males (rg = 0.084, p = 0.007) and trends toward positive genetic correlation with bioavailable testosterone (rg = 0.060, p = 0.084) and SHBG in males (rg = 0.056, p = 0.086) and with albumin in a sex-combined cohort (rg = 0.082, p = 0.015); however in females, we observed positive genetic correlation with SHBG (rg = 0.089, p = 0.004) and a trend toward negative genetic correlation (i.e. genetic effects in opposite directions) with bioavailable testosterone (rg = -0.064, p = 0.032). For alcohol dependence, we observed a trend toward negative genetic correlation with total testosterone in females (rg = -0.106, p = 0.024) and positive genetic correlation with BMI-adjusted SHBG in males (rg = 0.119, p = 0.017). Several of these genetic correlations differed between females and males and were not in the same direction as the corresponding phenotypic associations.

Conclusions

Findings suggest that shared genetic effects may contribute to positive associations of alcohol consumption with albumin in both sexes, as well as positive associations between alcohol consumption and bioavailable testosterone and between alcohol dependence and SHBG in males. However, relative contributions of heritable and environmental factors to associations between alcohol-use traits and sex-hormone levels may differ by sex, with genetic factors contributing more in males and environmental factors contributing more in females.

Sex-Hormone-Binding Globulin Gene Polymorphisms and Breast Cancer Risk in Caucasian Women of Russia

In our work, the associations of GWAS (genome-wide associative studies) impact for sex-hormone-binding globulin (SHBG)-level SNPs with the risk of breast cancer (BC) in the cohort of Caucasian women of Russia were assessed. The work was performed on a sample of 1498 women (358 BC patients and 1140 control (non BC) subjects). SHBG correlated in previously GWAS nine polymorphisms such as rs780093 GCKR, rs17496332 PRMT6, rs3779195 BAIAP2L1, rs10454142 PPP1R21, rs7910927 JMJD1C, rs4149056 SLCO1B1, rs440837 ZBTB10, rs12150660 SHBG, and rs8023580 NR2F2 have been genotyped. BC risk effects of allelic and non-allelic SHBG-linked gene SNPs interactions were detected by regression analysis. The risk genetic factor for BC developing is an SHBG-lowering allele variant C rs10454142 PPP1R21 ([additive genetic model] OR = 1.31; 95%CI = 1.08-1.65; pperm = 0.024; power = 85.26%), which determines 0.32% of the cancer variance. Eight of the nine studied SHBG-related SNPs have been involved in cancer susceptibility as part of nine different non-allelic gene interaction models, the greatest contribution to which is made by rs10454142 PPP1R21 (included in all nine models, 100%) and four more SNPs-rs7910927 JMJD1C (five models, 55.56%), rs17496332 PRMT6 (four models, 44.44%), rs780093 GCKR (four models, 44.44%), and rs440837 ZBTB10 (four models, 44.44%). For SHBG-related loci, pronounced functionality in the organism (including breast, liver, fibroblasts, etc.) was predicted in silico, having a direct relationship through many pathways with cancer pathophysiology. In conclusion, our results demonstrated the involvement of SHBG-correlated genes polymorphisms in BC risk in Caucasian women in Russia.

Large-scale integrative analysis of juvenile idiopathic arthritis for new insight into its pathogenesis

BACKGROUND

Juvenile idiopathic arthritis (JIA) is one of the most prevalent rheumatic disorders in children and is classified as an autoimmune disease (AID). While a robust genetic contribution to JIA etiology has been established, the exact pathogenesis remains unclear.

METHODS

To prioritize biologically interpretable susceptibility genes and proteins for JIA, we conducted transcriptome-wide and proteome-wide association studies (TWAS/PWAS). Then, to understand the genetic architecture of JIA, we systematically analyzed single-nucleotide polymorphism (SNP)-based heritability, a signature of natural selection, and polygenicity. Next, we conducted HLA typing using multi-ethnicity RNA sequencing data. Additionally, we examined the T cell receptor (TCR) repertoire at a single-cell level to explore the potential links between immunity and JIA risk.

RESULTS

We have identified 19 TWAS genes and two PWAS proteins associated with JIA risks. Furthermore, we observe that the heritability and cell type enrichment analysis of JIA are enriched in T lymphocytes and HLA regions and that JIA shows higher polygenicity compared to other AIDs. In multi-ancestry HLA typing, B*45:01 is more prevalent in African JIA patients than in European JIA patients, whereas DQA1*01:01, DQA1*03:01, and DRB1*04:01 exhibit a higher frequency in European JIA patients. Using single-cell immune repertoire analysis, we identify clonally expanded T cell subpopulations in JIA patients, including CXCL13+BHLHE40+ TH cells which are significantly associated with JIA risks.

CONCLUSION

Our findings shed new light on the pathogenesis of JIA and provide a strong foundation for future mechanistic studies aimed at uncovering the molecular drivers of JIA.

Are novel treatments for brain disorders hiding in plain sight?

Psychiatric diseases are strongly influenced by genetics, but genetically guided treatments have been slow to develop, and precise molecular mechanisms remain mysterious. Although individual locations in the genome tend to not contribute powerfully to psychiatric disease incidence, genome-wide association studies (GWAS) have now successfully linked hundreds of specific genetic loci to psychiatric disorders [1-3]. Here, building upon results from well-powered GWAS of four phenotypes relevant to psychiatry, we motivate an exploratory workflow leading from GWAS screening, through causal testing in animal models using methods such as optogenetics, to new therapies in human beings. We focus on schizophrenia and the dopamine D2 receptor (DRD2), hot flashes and the neurokinin B receptor (TACR3), cigarette smoking and receptors bound by nicotine (CHRNA5, CHRNA3, CHRNB4), and alcohol use and enzymes that help to break down alcohol (ADH1B, ADH1C, ADH7). A single genomic locus may not powerfully determine disease at the level of the population, but the same locus may nevertheless represent a potent treatment target suitable for population-wide therapeutic approaches.

Chargaff's second parity rule lies at the origin of additive genetic interactions in quantitative traits to make omnigenic selection possible

Background

Francis Crick's central dogma provides a residue-by-residue mechanistic explanation of the flow of genetic information in living systems. However, this principle may not be sufficient for explaining how random mutations cause continuous variation of quantitative highly polygenic complex traits. Chargaff's second parity rule (CSPR), also referred to as intrastrand DNA symmetry, defined as near-exact equalities G ≈ C and A ≈ T within a single DNA strand, is a statistical property of cellular genomes. The phenomenon of intrastrand DNA symmetry was discovered more than 50 years ago; at present, it remains unclear what its biological role is, what the mechanisms are that force cellular genomes to comply strictly with CSPR, and why genomes of certain noncellular organisms have broken intrastrand DNA symmetry. The present work is aimed at studying a possible link between intrastrand DNA symmetry and the origin of genetic interactions in quantitative traits.

Methods

Computational analysis of single-nucleotide polymorphisms in human and mouse populations and of nucleotide composition biases at different codon positions in bacterial and human proteomes.

Results

The analysis of mutation spectra inferred from single-nucleotide polymorphisms observed in murine and human populations revealed near-exact equalities of numbers of reverse complementary mutations, indicating that random genetic variations obey CSPR. Furthermore, nucleotide compositions of coding sequences proved to be statistically interwoven via CSPR because pyrimidine bias at the 3rd codon position compensates purine bias at the 1st and 2nd positions.

Conclusions

According to Fisher's infinitesimal model, we propose that accumulation of reverse complementary mutations results in a continuous phenotypic variation due to small additive effects of statistically interwoven genetic variations. Therefore, additive genetic interactions can be inferred as a statistical entanglement of nucleotide compositions of separate genetic loci. CSPR challenges the neutral theory of molecular evolution-because all random mutations participate in variation of a trait-and provides an alternative solution to Haldane's dilemma by making a gene function diffuse. We propose that CSPR is symmetry of Fisher's infinitesimal model and that genetic information can be transferred in an implicit contactless manner.

Polygenic risk score phenome-wide association study reveals an association between endometriosis and testosterone

BACKGROUND

Endometriosis affects 1 in 9 women, yet it is poorly understood with long diagnostic delays, invasive diagnoses, and poor treatment outcomes. Characterised by the presence of endometrial-like tissue outside of the uterus, its main symptoms are pain and infertility. Endometriosis often co-occurs with other conditions, which may provide insights into the origins of endometriosis.

METHODS

Here a polygenic risk score phenome-wide association study of endometriosis was conducted in the UK Biobank to investigate the pleiotropic effects of a genetic liability to endometriosis. The relationship between the polygenic risk score for endometriosis and health conditions, blood and urine biomarkers and reproductive factors were investigated separately in females, males and females without an endometriosis diagnosis. The relationship between endometriosis and the blood and urine biomarkers was further investigated using genetic correlation and Mendelian randomisation approaches to identify causal relationships.

RESULTS

Multiple health conditions, blood and urine biomarkers and reproductive factors were associated with genetic liability to endometriosis in each group, indicating many endometriosis comorbidities are not dependent on the physical manifestation of endometriosis. Differences in the associated traits between males and females highlighted the importance of sex-specific pathways in the overlap of endometriosis with many other traits. Notably, an association of genetic liability to endometriosis with lower testosterone levels was identified. Follow-up analysis utilising Mendelian randomisation approaches suggested lower testosterone may be causal for both endometriosis and clear cell ovarian cancer.

CONCLUSIONS

This study highlights the diversity of the pleiotropic effects of genetic risk to endometriosis irrespective of a diagnosis of endometriosis. A key finding was the identification of a causal effect of the genetic liability to lower testosterone on endometriosis using Mendelian randomisation.

Speos: an ensemble graph representation learning framework to predict core gene candidates for complex diseases

Understanding phenotype-to-genotype relationships is a grand challenge of 21st century biology with translational implications. The recently proposed "omnigenic" model postulates that effects of genetic variation on traits are mediated by core-genes and -proteins whose activities mechanistically influence the phenotype, whereas peripheral genes encode a regulatory network that indirectly affects phenotypes via core gene products. Here, we develop a positive-unlabeled graph representation-learning ensemble-approach based on a nested cross-validation to predict core-like genes for diverse diseases using Mendelian disorder genes for training. Employing mouse knockout phenotypes for external validations, we demonstrate that core-like genes display several key properties of core genes: Mouse knockouts of genes corresponding to our most confident predictions give rise to relevant mouse phenotypes at rates on par with the Mendelian disorder genes, and all candidates exhibit core gene properties like transcriptional deregulation in disease and loss-of-function intolerance. Moreover, as predicted for core genes, our candidates are enriched for drug targets and druggable proteins. In contrast to Mendelian disorder genes the new core-like genes are enriched for druggable yet untargeted gene products, which are therefore attractive targets for drug development. Interpretation of the underlying deep learning model suggests plausible explanations for our core gene predictions in form of molecular mechanisms and physical interactions. Our results demonstrate the potential of graph representation learning for the interpretation of biological complexity and pave the way for studying core gene properties and future drug development.

Factorizing polygenic epistasis improves prediction and uncovers biological pathways in complex traits

Epistasis is central in many domains of biology, but it has not yet been proven useful for understanding the etiology of complex traits. This is partly because complex-trait epistasis involves polygenic interactions that are poorly captured in current models. To address this gap, we developed a model called Epistasis Factor Analysis (EFA). EFA assumes that polygenic epistasis can be factorized into interactions between a few epistasis factors (EFs), which represent latent polygenic components of the observed complex trait. The statistical goals of EFA are to improve polygenic prediction and to increase power to detect epistasis, while the biological goal is to unravel genetic effects into more-homogeneous units. We mathematically characterize EFA and use simulations to show that EFA outperforms current epistasis models when its assumptions approximately hold. Applied to predicting yeast growth rates, EFA outperforms the additive model for several traits with large epistasis heritability and uniformly outperforms the standard epistasis model. We replicate these prediction improvements in a second dataset. We then apply EFA to four previously characterized traits in the UK Biobank and find statistically significant epistasis in all four, including two that are robust to scale transformation. Moreover, we find that the inferred EFs partly recover pre-defined biological pathways for two of the traits. Our results demonstrate that more realistic models can identify biologically and statistically meaningful epistasis in complex traits, indicating that epistasis has potential for precision medicine and characterizing the biology underlying GWAS results.

Structural variants and short tandem repeats impact gene expression and splicing in bovine testis tissue

Structural variants (SVs) and short tandem repeats (STRs) are significant sources of genetic variation. However, the impacts of these variants on gene regulation have not been investigated in cattle. Here, we genotyped and characterized 19,408 SVs and 374,821 STRs in 183 bovine genomes and investigated their impact on molecular phenotypes derived from testis transcriptomes. We found that 71% STRs were multiallelic. The vast majority (95%) of STRs and SVs were in intergenic and intronic regions. Only 37% SVs and 40% STRs were in high linkage disequilibrium (LD) (R2 > 0.8) with surrounding SNPs/insertions and deletions (Indels), indicating that SNP-based association testing and genomic prediction are blind to a nonnegligible portion of genetic variation. We showed that both SVs and STRs were more than 2-fold enriched among expression and splicing QTL (e/sQTL) relative to SNPs/Indels and were often associated with differential expression and splicing of multiple genes. Deletions and duplications had larger impacts on splicing and expression than any other type of SV. Exonic duplications predominantly increased gene expression either through alternative splicing or other mechanisms, whereas expression- and splicing-associated STRs primarily resided in intronic regions and exhibited bimodal effects on the molecular phenotypes investigated. Most e/sQTL resided within 100 kb of the affected genes or splicing junctions. We pinpoint candidate causal STRs and SVs associated with the expression of SLC13A4 and TTC7B and alternative splicing of a lncRNA and CAPP1. We provide a catalog of STRs and SVs for taurine cattle and show that these variants contribute substantially to gene expression and splicing variation.

Sex-biased gene and microRNA expression in the developing mouse brain is associated with neurodevelopmental functions and neurological phenotypes

BACKGROUND

Sex differences pose a challenge and an opportunity in biomedical research. Understanding how sex chromosomes and hormones affect disease-causing mechanisms will shed light on the mechanisms underlying predominantly idiopathic sex-biased neurodevelopmental disorders such as ADHD, schizophrenia, and autism. Gene expression is a crucial conduit for the influence of sex on developmental processes; therefore, this study focused on sex differences in gene expression and the regulation of gene expression. The increasing interest in microRNAs (miRNAs), small, non-coding RNAs, for their contribution to normal and pathological neurodevelopment prompted us to test how miRNA expression differs between the sexes in the developing brain.

METHODS

High-throughput sequencing approaches were used to identify transcripts, including miRNAs, that showed significantly different expression between male and female brains on day 15.5 of development (E15.5).

RESULTS

Robust sex differences were identified for some genes and miRNAs, confirming the influence of biological sex on RNA. Many miRNAs that exhibit the greatest differences between males and females have established roles in neurodevelopment, implying that sex-biased expression may drive sex differences in developmental processes. In addition to highlighting sex differences for individual miRNAs, gene ontology analysis suggested several broad categories in which sex-biased RNAs might act to establish sex differences in the embryonic mouse brain. Finally, mining publicly available SNP data indicated that some sex-biased miRNAs reside near the genomic regions associated with neurodevelopmental disorders.

CONCLUSIONS

Together, these findings reinforce the importance of cataloguing sex differences in molecular biology research and highlight genes, miRNAs, and pathways of interest that may be important for sexual differentiation in the mouse and possibly the human brain.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

A genome-wide genetic screen uncovers determinants of human pigmentation

Skin color, one of the most diverse human traits, is determined by the quantity, type, and distribution of melanin. In this study, we leveraged the light-scattering properties of melanin to conduct a genome-wide screen for regulators of melanogenesis. We identified 169 functionally diverse genes that converge on melanosome biogenesis, endosomal transport, and gene regulation, of which 135 represented previously unknown associations with pigmentation. In agreement with their melanin-promoting function, the majority of screen hits were up-regulated in melanocytes from darkly pigmented individuals. We further unraveled functions of KLF6 as a transcription factor that regulates melanosome maturation and pigmentation in vivo, and of the endosomal trafficking protein COMMD3 in modulating melanosomal pH. Our study reveals a plethora of melanin-promoting genes, with broad implications for human variation, cell biology, and medicine.

The genetics of a "femaleness/maleness" score in cardiometabolic traits in the UK biobank

We recently devised continuous "sex-scores" that sum up multiple quantitative traits, weighted by their respective sex-difference effect sizes, as an approach to estimating polyphenotypic "maleness/femaleness" within each binary sex. To identify the genetic architecture underlying these sex-scores, we conducted sex-specific genome-wide association studies (GWASs) in the UK Biobank cohort (females: n = 161,906; males: n = 141,980). As a control, we also conducted GWASs of sex-specific "sum-scores", simply aggregating the same traits, without weighting by sex differences. Among GWAS-identified genes, while sum-score genes were enriched for genes differentially expressed in the liver in both sexes, sex-score genes were enriched for genes differentially expressed in the cervix and across brain tissues, particularly for females. We then considered single nucleotide polymorphisms with significantly different effects (sdSNPs) between the sexes for sex-scores and sum-scores, mapping to male-dominant and female-dominant genes. Here, we identified brain-related enrichment for sex-scores, especially for male-dominant genes; these findings were present but weaker for sum-scores. Genetic correlation analyses of sex-biased diseases indicated that both sex-scores and sum-scores were associated with cardiometabolic, immune, and psychiatric disorders.

Amplification is the primary mode of gene-by-sex interaction in complex human traits

Sex differences in complex traits are suspected to be in part due to widespread gene-by-sex interactions (GxSex), but empirical evidence has been elusive. Here, we infer the mixture of ways in which polygenic effects on physiological traits covary between males and females. We find that GxSex is pervasive but acts primarily through systematic sex differences in the magnitude of many genetic effects ("amplification") rather than in the identity of causal variants. Amplification patterns account for sex differences in trait variance. In some cases, testosterone may mediate amplification. Finally, we develop a population-genetic test linking GxSex to contemporary natural selection and find evidence of sexually antagonistic selection on variants affecting testosterone levels. Our results suggest that amplification of polygenic effects is a common mode of GxSex that may contribute to sex differences and fuel their evolution.

Sex-Specific Features of the Correlation between GWAS-Noticeable Polymorphisms and Hypertension in Europeans of Russia

The aim of the study was directed at studying the sex-specific features of the correlation between genome-wide association studies (GWAS)-noticeable polymorphisms and hypertension (HTN). In two groups of European subjects of Russia (n = 1405 in total), such as men (n = 821 in total: n = 564 HTN, n = 257 control) and women (n = 584 in total: n = 375 HTN, n = 209 control), the distribution of ten specially selected polymorphisms (they have confirmed associations of GWAS level with blood pressure (BP) parameters and/or HTN in Europeans) has been considered. The list of studied loci was as follows: (PLCE1) rs932764 A > G, (AC026703.1) rs1173771 G > A, (CERS5) rs7302981 G > A, (HFE) rs1799945 C > G, (OBFC1) rs4387287 C > A, (BAG6) rs805303 G > A, (RGL3) rs167479 T > G, (ARHGAP42) rs633185 C > G, (TBX2) rs8068318 T > C, and (ATP2B1) rs2681472 A > G. The contribution of individual loci and their inter-locus interactions to the HTN susceptibility with bioinformatic interpretation of associative links was evaluated separately in men's and women's cohorts. The men-women differences in involvement in the disease of the BP/HTN-associated GWAS SNPs were detected. Among women, the HTN risk has been associated with HFE rs1799945 C > G (genotype GG was risky; OR_GG = 11.15 p_permGG = 0.014) and inter-locus interactions of all 10 examined SNPs as part of 26 intergenic interactions models. In men, the polymorphism BAG6 rs805303 G > A (genotype AA was protective; OR_AA = 0.30 p_permAA = 0.0008) and inter-SNPs interactions of eight loci in only seven models have been founded as HTN-correlated. HTN-linked loci and strongly linked SNPs were characterized by pronounced polyvector functionality in both men and women, but at the same time, signaling pathways of HTN-linked genes/SNPs in women and men were similar and were represented mainly by immune mechanisms. As a result, the present study has demonstrated a more pronounced contribution of BP/HTN-associated GWAS SNPs to the HTN susceptibility (due to weightier intergenic interactions) in European women than in men.

Genetic effects on the timing of parturition and links to fetal birth weight

The timing of parturition is crucial for neonatal survival and infant health. Yet, its genetic basis remains largely unresolved. We present a maternal genome-wide meta-analysis of gestational duration (n = 195,555), identifying 22 associated loci (24 independent variants) and an enrichment in genes differentially expressed during labor. A meta-analysis of preterm delivery (18,797 cases, 260,246 controls) revealed six associated loci and large genetic similarities with gestational duration. Analysis of the parental transmitted and nontransmitted alleles (n = 136,833) shows that 15 of the gestational duration genetic variants act through the maternal genome, whereas 7 act both through the maternal and fetal genomes and 2 act only via the fetal genome. Finally, the maternal effects on gestational duration show signs of antagonistic pleiotropy with the fetal effects on birth weight: maternal alleles that increase gestational duration have negative fetal effects on birth weight. The present study provides insights into the genetic effects on the timing of parturition and the complex maternal-fetal relationship between gestational duration and birth weight.

Effects of protein-coding variants on blood metabolite measurements and clinical biomarkers in the UK Biobank

Genome-wide association studies (GWASs) have established the contribution of common and low-frequency variants to metabolic blood measurements in the UK Biobank (UKB). To complement existing GWAS findings, we assessed the contribution of rare protein-coding variants in relation to 355 metabolic blood measurements-including 325 predominantly lipid-related nuclear magnetic resonance (NMR)-derived blood metabolite measurements (Nightingale Health Plc) and 30 clinical blood biomarkers-using 412,393 exome sequences from four genetically diverse ancestries in the UKB. Gene-level collapsing analyses were conducted to evaluate a diverse range of rare-variant architectures for the metabolic blood measurements. Altogether, we identified significant associations (p < 1 × 10^-8) for 205 distinct genes that involved 1,968 significant relationships for the Nightingale blood metabolite measurements and 331 for the clinical blood biomarkers. These include associations for rare non-synonymous variants in PLIN1 and CREB3L3 with lipid metabolite measurements and SYT7 with creatinine, among others, which may not only provide insights into novel biology but also deepen our understanding of established disease mechanisms. Of the study-wide significant clinical biomarker associations, 40% were not previously detected on analyzing coding variants in a GWAS in the same cohort, reinforcing the importance of studying rare variation to fully understand the genetic architecture of metabolic blood measurements.

Genetic Risk, Adherence to a Healthy Lifestyle, and Hyperuricemia: The TCLSIH Cohort Study

BACKGROUND

Genetic factors have been associated with hyperuricemia in large studies, but the extent to which this can be offset by a healthy lifestyle is unknown. This study aimed to examine whether healthy lifestyle could reduce hyperuricemia risk among individuals with different genetic profiles.

METHODS

We defined a lifestyle score using body mass index, smoking, alcohol consumption, physical activities, and diets in 2796 unrelated individuals from the Tianjin Chronic Low-grade Systemic Inflammation and Health (TCLSIH) cohort study. Polygenic risk scores (PRS) were constructed based on uric acid loci. Associations of combined lifestyle factors and genetic risk and incident hyperuricemia were estimated using Cox proportional hazard regression.

RESULTS

Of 2796 individuals, 747 participants (26.7%) developed hyperuricemia. Genetic risk and lifestyle were predictors of incident events, and they showed an interaction for the outcome. Compared with high PRS, low PRS reduced risk of incident hyperuricemia by 40%, and compared with unhealthy lifestyle, healthy lifestyle reduced risk of incident hyperuricemia by 41%. Compared with unhealthy lifestyle and high genetic risk, adherence to healthy lifestyle was associated with a 68% (95% confidence interval, 44%-81%) lower risk of hyperuricemia among those at a low genetic risk.

CONCLUSIONS

In this prospective cohort study, we observed an interaction between genetics and lifestyle and the risk of hyperuricemia. The public health implication is that a healthy lifestyle is important for hyperuricemia prevention, especially for individuals with high genetic risk scores.

Genetic analyses implicate complex links between adult testosterone levels and health and disease

BACKGROUND

Testosterone levels are linked with diverse characteristics of human health, yet, whether these associations reflect correlation or causation remains debated. Here, we provide a broad perspective on the role of genetically determined testosterone on complex diseases in both sexes.

METHODS

Leveraging genetic and health registry data from the UK Biobank and FinnGen (total N = 625,650), we constructed polygenic scores (PGS) for total testosterone, sex-hormone binding globulin (SHBG) and free testosterone, associating these with 36 endpoints across different disease categories in the FinnGen. These analyses were combined with Mendelian Randomization (MR) and cross-sex PGS analyses to address causality.

RESULTS

We show testosterone and SHBG levels are intricately tied to metabolic health, but report lack of causality behind most associations, including type 2 diabetes (T2D). Across other disease domains, including 13 behavioral and neurological diseases, we similarly find little evidence for a substantial contribution from normal variation in testosterone levels. We nonetheless find genetically predicted testosterone affects many sex-specific traits, with a pronounced impact on female reproductive health, including causal contribution to PCOS-related traits like hirsutism and post-menopausal bleeding (PMB). We also illustrate how testosterone levels associate with antagonistic effects on stroke risk and reproductive endpoints between the sexes.

CONCLUSIONS

Overall, these findings provide insight into how genetically determined testosterone correlates with several health parameters in both sexes. Yet the lack of evidence for a causal contribution to most traits beyond sex-specific health underscores the complexity of the mechanisms linking testosterone levels to disease risk and sex differences.

Probabilistic integration of transcriptome-wide association studies and colocalization analysis identifies key molecular pathways of complex traits

Integrative genetic association methods have shown great promise in post-GWAS (genome-wide association study) analyses, in which one of the most challenging tasks is identifying putative causal genes and uncovering molecular mechanisms of complex traits. Recent studies suggest that prevailing computational approaches, including transcriptome-wide association studies (TWASs) and colocalization analysis, are individually imperfect, but their joint usage can yield robust and powerful inference results. This paper presents INTACT, a computational framework to integrate probabilistic evidence from these distinct types of analyses and implicate putative causal genes. This procedure is flexible and can work with a wide range of existing integrative analysis approaches. It has the unique ability to quantify the uncertainty of implicated genes, enabling rigorous control of false-positive discoveries. Taking advantage of this highly desirable feature, we further propose an efficient algorithm, INTACT-GSE, for gene set enrichment analysis based on the integrated probabilistic evidence. We examine the proposed computational methods and illustrate their improved performance over the existing approaches through simulation studies. We apply the proposed methods to analyze the multi-tissue eQTL data from the GTEx project and eight large-scale complex- and molecular-trait GWAS datasets from multiple consortia and the UK Biobank. Overall, we find that the proposed methods markedly improve the existing putative gene implication methods and are particularly advantageous in evaluating and identifying key gene sets and biological pathways underlying complex traits.

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.

Taking quantitative genomics into the wild

We organized this special issue to highlight new work and review recent advances at the cutting edge of 'wild quantitative genomics'. In this editorial, we will present some history of wild quantitative genetic and genomic studies, before discussing the main themes in the papers published in this special issue and highlighting the future outlook of this dynamic field.

How broad is the selfing syndrome? Insights from convergent evolution of gene expression across species and tissues in the Capsella genus

The shift from outcrossing to selfing is one of the main evolutionary transitions in plants. It is accompanied by profound effects on reproductive traits, the so-called selfing syndrome. Because the transition to selfing also implies deep genomic and ecological changes, one also expects to observe a genomic selfing syndrome. We took advantage of the three independent transitions from outcrossing to selfing in the Capsella genus to characterize the overall impact of mating system change on RNA expression, in flowers but also in leaves and roots. We quantified the extent of both selfing and genomic syndromes, and tested whether changes in expression corresponded to adaptation to selfing or to relaxed selection on traits that were constrained in outcrossers. Mating system change affected gene expression in all three tissues but more so in flowers than in roots and leaves. Gene expression in selfing species tended to converge in flowers but diverged in the two other tissues. Hence, convergent adaptation to selfing dominates in flowers, whereas genetic drift plays a more important role in leaves and roots. The effect of mating system transition is not limited to reproductive tissues and corresponds to both adaptation to selfing and relaxed selection on previously constrained traits.

Multi-omics peripheral and core regions of cancer

Thousands of genes are perturbed by cancer, and these disturbances can be seen in transcriptome, methylation, somatic mutation, and copy number variation omics studies. Understanding their connectivity patterns as an omnigenic neighbourhood in a molecular interaction network (interactome) is a key step towards advancing knowledge of the molecular mechanisms underlying cancers. Here, we introduce a unified connectivity line (CLine) to pinpoint omics-specific omnigenic patterns across 15 curated cancers. Taking advantage of the universality of CLine, we distinguish the peripheral and core genes for each omics aspect. We propose a network-based framework, multi-omics periphery and core (MOPC), to combine peripheral and core genes from different omics into a button-like structure. On the basis of network proximity, we provide evidence that core genes tend to be specifically perturbed in one omics, but the peripheral genes are diversely perturbed in multiple omics. And the core of one omics is regulated by multiple omics peripheries. Finally, we take the MOPC as an omnigenic neighbourhood, describe its characteristics, and explore its relative contribution to network-based mechanisms of cancer. We were able to present how multi-omics perturbations percolate through the human interactome and contribute to an integrated periphery and core.

Sex Hormone Candidate Gene Polymorphisms Are Associated with Endometriosis

The present study was designed to examine whether sex hormone polymorphisms proven by GWAS are associated with endometriosis risk. Unrelated female participants totaling 1376 in number (395 endometriosis patients and 981 controls) were recruited into the study. Nine single-nucleotide polymorphisms (SNPs) which GWAS correlated with circulating levels of sex hormones were genotyped using a TaqMan allelic discrimination assay. FSH-lowering, and LH- and testosterone-heightening polymorphisms of the FSHB promoter (allelic variants A rs11031002 and C rs11031005) exhibit a protective effect for endometriosis (OR = 0.60-0.68). By contrast, the TT haplotype loci that were GWAS correlated with higher FSH levels and lower LH and testosterone concentrations determined an increased risk for endometriosis (OR = 2.03). Endometriosis-involved epistatic interactions were found between eight loci of sex hormone genes (without rs148982377 ZNF789) within twelve genetic simulation models. In silico examination established that 8 disorder-related loci and 80 proxy SNPs are genome variants affecting the expression, splicing, epigenetic and amino acid conformation of the 34 genes which enrich the organic anion transport and secondary carrier transporter pathways. In conclusion, the present study showed that sex hormone polymorphisms proven by GWAS are associated with endometriosis risk and involved in the molecular pathophysiology of the disease due to their functionality.

Leveraging global multi-ancestry meta-analysis in the study of idiopathic pulmonary fibrosis genetics

The research of rare and devastating orphan diseases, such as idiopathic pulmonary fibrosis (IPF) has been limited by the rarity of the disease itself. The prognosis is poor-the prevalence of IPF is only approximately four times the incidence, limiting the recruitment of patients to trials and studies of the underlying biology. Global biobanking efforts can dramatically alter the future of IPF research. We describe a large-scale meta-analysis of IPF, with 8,492 patients and 1,355,819 population controls from 13 biobanks around the globe. Finally, we combine this meta-analysis with the largest available meta-analysis of IPF, reaching 11,160 patients and 1,364,410 population controls. We identify seven novel genome-wide significant loci, only one of which would have been identified if the analysis had been limited to European ancestry individuals. We observe notable pleiotropy across IPF susceptibility and severe COVID-19 infection and note an unexplained sex-heterogeneity effect at the strongest IPF locus MUC5B.

A Shifting Relationship Between Sex Hormone-Binding Globulin and Total Testosterone Across Puberty in Boys

CONTEXT

Sex hormone-binding globulin (SHBG) is associated with levels of total testosterone (total-T), and both total-T and SHBG are associated with obesity.

OBJECTIVE

We aimed to clarify the nature of the relationship between testosterone and SHBG and improve our understanding of their relationships with obesity. We hypothesize that the hypothalamic-pituitary-gonadal axis contributes to the homeostasis of testosterone by increasing the production of gonadal testosterone through a feedback mechanism that might operate differently at different pubertal stages.

METHODS

We investigated the dynamics of the relationship between SHBG, total-T, and body mass index (BMI) throughout puberty (from age 9 to 17) using longitudinal data obtained in 507 males. The directionality of this relationship was explored using polygenic scores of SHBG and total-T, and a two-sample Mendelian Randomization (MR) in male adults.

RESULTS

Consistent with our hypothesis, we found positive relationships between SHBG and total-T at age 15 and 17 but either no relationship or a negative relationship during the earlier time points. Such shifting relationships explained age-related changes in the association between total-T and BMI. Polygenic scores of SHBG and total-T in mediation analyses and the two-sample MR in male adults suggested an effect of SHBG on total-T but also a somewhat weaker effect of total-T on SHBG. Two-sample MR also showed an effect of BMI on SHBG but no effect of SHBG on BMI.

CONCLUSION

These results clarify the nature of the relationship between testosterone and SHBG during puberty and adulthood and shed new light on their possible relationship with obesity.

Integrative analysis of metabolite GWAS illuminates the molecular basis of pleiotropy and genetic correlation

Pleiotropy and genetic correlation are widespread features in genome-wide association studies (GWAS), but they are often difficult to interpret at the molecular level. Here, we perform GWAS of 16 metabolites clustered at the intersection of amino acid catabolism, glycolysis, and ketone body metabolism in a subset of UK Biobank. We utilize the well-documented biochemistry jointly impacting these metabolites to analyze pleiotropic effects in the context of their pathways. Among the 213 lead GWAS hits, we find a strong enrichment for genes encoding pathway-relevant enzymes and transporters. We demonstrate that the effect directions of variants acting on biology between metabolite pairs often contrast with those of upstream or downstream variants as well as the polygenic background. Thus, we find that these outlier variants often reflect biology local to the traits. Finally, we explore the implications for interpreting disease GWAS, underscoring the potential of unifying biochemistry with dense metabolomics data to understand the molecular basis of pleiotropy in complex traits and diseases.

Singling out motor neurons in the age of single-cell transcriptomics

Motor neurons are a remarkably powerful cell type in the central nervous system. They innervate and control the contraction of virtually every muscle in the body and their dysfunction underlies numerous neuromuscular diseases. Some motor neurons seem resistant to degeneration whereas others are vulnerable. The intrinsic heterogeneity of motor neurons in adult organisms has remained elusive. The development of high-throughput single-cell transcriptomics has changed the paradigm, empowering rapid isolation and profiling of motor neuron nuclei, revealing remarkable transcriptional diversity within the skeletal and autonomic nervous systems. Here, we discuss emerging technologies for defining motor neuron heterogeneity in the adult motor system as well as implications for disease and spinal cord injury. We establish a roadmap for future applications of emerging techniques - such as epigenetic profiling, spatial RNA sequencing, and single-cell somatic mutational profiling to adult motor neurons, which will revolutionize our understanding of the healthy and degenerating adult motor system.

Genome-wide meta-analysis of insomnia prioritizes genes associated with metabolic and psychiatric pathways

Insomnia is a heritable, highly prevalent sleep disorder for which no sufficient treatment currently exists. Previous genome-wide association studies with up to 1.3 million subjects identified over 200 associated loci. This extreme polygenicity suggested that many more loci remain to be discovered. The current study almost doubled the sample size to 593,724 cases and 1,771,286 controls, thereby increasing statistical power, and identified 554 risk loci (including 364 novel loci). To capitalize on this large number of loci, we propose a novel strategy to prioritize genes using external biological resources and functional interactions between genes across risk loci. Of all 3,898 genes naively implicated from the risk loci, we prioritize 289 and find brain-tissue expression specificity and enrichment in specific gene sets of synaptic signaling functions and neuronal differentiation. We show that this novel gene prioritization strategy yields specific hypotheses on underlying mechanisms of insomnia that would have been missed by traditional approaches.

Integrated analysis of direct and proxy genome wide association studies highlights polygenicity of Alzheimer's disease outside of the APOE region

Recent meta-analyses combining direct genome-wide association studies (GWAS) with those of family history (GWAX) have indicated very low SNP heritability of Alzheimer's disease (AD). These low estimates may call into question the prospects of continued progress in genetic discovery for AD within the spectrum of common variants. We highlight dramatic downward biases in previous methods, and we validate a novel method for the estimation of SNP heritability via integration of GWAS and GWAX summary data. We apply our method to investigate the genetic architecture of AD using GWAX from UK Biobank and direct case-control GWAS from the International Genomics of Alzheimer's Project (IGAP). We estimate the liability scale common variant SNP heritability of Clinical AD outside of APOE region at ~7-11%, and we project the corresponding estimate for AD pathology to be up to approximately 23%. We estimate that nearly 90% of common variant SNP heritability of Clinical AD exists outside the APOE region. Rare variants not tagged in standard GWAS may account for additional variance. Our results indicate that, while GWAX for AD in UK Biobank may result in greater attenuation of genetic effects beyond that conventionally assumed, it does not introduce appreciable contamination of signal by genetically distinct traits relative to direct case-control GWAS in IGAP. Genetic risk for AD represents a strong effect of APOE superimposed upon a highly polygenic background.

Enrichment analyses identify shared associations for 25 quantitative traits in over 600,000 individuals from seven diverse ancestries

Since 2005, genome-wide association (GWA) datasets have been largely biased toward sampling European ancestry individuals, and recent studies have shown that GWA results estimated from self-identified European individuals are not transferable to non-European individuals because of various confounding challenges. Here, we demonstrate that enrichment analyses that aggregate SNP-level association statistics at multiple genomic scales-from genes to genomic regions and pathways-have been underutilized in the GWA era and can generate biologically interpretable hypotheses regarding the genetic basis of complex trait architecture. We illustrate examples of the robust associations generated by enrichment analyses while studying 25 continuous traits assayed in 566,786 individuals from seven diverse self-identified human ancestries in the UK Biobank and the Biobank Japan as well as 44,348 admixed individuals from the PAGE consortium including cohorts of African American, Hispanic and Latin American, Native Hawaiian, and American Indian/Alaska Native individuals. We identify 1,000 gene-level associations that are genome-wide significant in at least two ancestry cohorts across these 25 traits as well as highly conserved pathway associations with triglyceride levels in European, East Asian, and Native Hawaiian cohorts.

Natura Non Facit Saltus: The Adaptive Significance of Arginine Vasopressin in Human Affect, Cognition, and Behavior

Hormones coordinate internal bodily systems with cognition, affect, and behavior, and thereby influence aspects of social interactions including cooperation, competition, isolation, and loneliness. The adaptive significance and contextuality of oxytocin (OXT) and testosterone (T) have been well-studied, but a unified theory and evolutionary framework for understanding the adaptive functions of arginine vasopressin (AVP) remain undeveloped. We propose and evaluate the hypothesis that AVP mediates adaptive variation in the presence and strength of social and sociosexual salience, attention and behavior specifically in situations that involve combinations of cooperation with conflict or competition. This hypothesis can help to explain the ancestral, original functions of AVP-like peptides, and their continuity with the current roles of AVP, for humans, in male-male competition, male-male reciprocity, male-to-female pair bonding, female-female interactions, social integration, and social attention and anxiety. In this context, social isolation and loneliness may be mediated by reduced abilities or interests in navigation of social opportunities and situations, due in part to low AVP levels or reactivity, and in part to reductions in levels of OXT-mediated social reward.

Computational Genomics in the Era of Precision Medicine: Applications to Variant Analysis and Gene Therapy

Rapid methodological advances in statistical and computational genomics have enabled researchers to better identify and interpret both rare and common variants responsible for complex human diseases. As we continue to see an expansion of these advances in the field, it is now imperative for researchers to understand the resources and methodologies available for various data types and study designs. In this review, we provide an overview of recent methods for identifying rare and common variants and understanding their roles in disease etiology. Additionally, we discuss the strategy, challenge, and promise of gene therapy. As computational and statistical approaches continue to improve, we will have an opportunity to translate human genetic findings into personalized health care.

Bench Research Informed by GWAS Results

Scientifically interesting as well as practically important phenotypes often belong to the realm of complex traits. To the extent that these traits are hereditary, they are usually 'highly polygenic'. The study of such traits presents a challenge for researchers, as the complex genetic architecture of such traits makes it nearly impossible to utilise many of the usual methods of reverse genetics, which often focus on specific genes. In recent years, thousands of genome-wide association studies (GWAS) were undertaken to explore the relationships between complex traits and a large number of genetic factors, most of which are characterised by tiny effects. In this review, we aim to familiarise 'wet biologists' with approaches for the interpretation of GWAS results, to clarify some issues that may seem counterintuitive and to assess the possibility of using GWAS results in experiments on various complex traits.

Integrative analysis of the plasma proteome and polygenic risk of cardiometabolic diseases

Cardiometabolic diseases are frequently polygenic in architecture, comprising a large number of risk alleles with small effects spread across the genome1-3. Polygenic scores (PGS) aggregate these into a metric representing an individual's genetic predisposition to disease. PGS have shown promise for early risk prediction4-7 and there is an open question as to whether PGS can also be used to understand disease biology8. Here, we demonstrate that cardiometabolic disease PGS can be used to elucidate the proteins underlying disease pathogenesis. In 3,087 healthy individuals, we found that PGS for coronary artery disease, type 2 diabetes, chronic kidney disease and ischaemic stroke are associated with the levels of 49 plasma proteins. Associations were polygenic in architecture, largely independent of cis and trans protein quantitative trait loci and present for proteins without quantitative trait loci. Over a follow-up of 7.7 years, 28 of these proteins associated with future myocardial infarction or type 2 diabetes events, 16 of which were mediators between polygenic risk and incident disease. Twelve of these were druggable targets with therapeutic potential. Our results demonstrate the potential for PGS to uncover causal disease biology and targets with therapeutic potential, including those that may be missed by approaches utilizing information at a single locus.

Strain-Specific Liver Metabolite Profiles in Medaka

The relationship between genetic variation and phenotypic traits is often poorly understood since specific genotypes do not always easily translate into associated phenotypes, especially for complex disorders. The genetic background has been shown to affect metabolic pathways and thus contribute to variations in the metabolome. Here, we tested the suitability of NMR metabolomics for comparative analysis of fish lines as a first step towards phenotype-genotype association studies. The Japanese rice fish, medaka (Oryzias latipes), is a widely used genetic vertebrate model with several isogenic inbred laboratory strains. We used liver extracts of medaka iCab and HO5 strains as a paradigm to test the feasibility of distinguishing the metabolome of two different inbred strains. Fifteen metabolites could be detected in uni- and multivariate analyses that showed strain-specific levels. Differences could be assigned to specific metabolic pathways. Our results show that NMR spectroscopy is a suitable method to detect variance of the metabolome caused by subtle genetic differences. Thus, it has the potential to address genotype–phenotype associations in medaka, providing an additional level of phenotypic analysis.

Discovery and implications of polygenicity of common diseases

The sequencing of the human genome has allowed the study of the genetic architecture of common diseases: the number of genomic variants that contribute to risk of disease and their joint frequency and effect size distribution. Common diseases are polygenic, with many loci contributing to phenotype, and the cumulative burden of risk alleles determines individual risk in conjunction with environmental factors. Most risk loci occur in noncoding regions of the genome regulating cell- and context-specific gene expression. Although the effect sizes of most risk alleles are small, their cumulative effects in individuals, quantified as a polygenic (risk) score, can identify people at increased risk of disease, thereby facilitating prevention or early intervention.

From variant to function in human disease genetics

[Figure: see text].

The genetics of testosterone contributes to "femaleness/maleness" of cardiometabolic traits and type 2 diabetes

The genetic architecture of testosterone is highly distinct between sexes. Moreover, obesity is associated with higher testosterone in females but lower testosterone in males. Here, we ask whether male-specific testosterone variants are associated with a male pattern of obesity and type 2 diabetes (T2D) in females, and vice versa. In the UK Biobank, we conducted sex-specific genome-wide association studies and computed polygenic scores for total (PGS_TT) and bioavailable testosterone (PGS_BT). We tested sex-congruent and sex-incongruent associations between sex-specific PGS_Ts and metabolic traits, as well as T2D diagnosis. Female-specific PGS_BT was associated with an elevated cardiometabolic risk and probability of T2D, in both sexes. Male-specific PGS_TT was associated with traits conferring a lower cardiometabolic risk and probability of T2D, in both sexes. We demonstrate the value in considering polygenic testosterone as sex-related continuous traits, in each sex.