Conserved pleiotropy of an ancient plant homeobox gene uncovered by cis-regulatory dissection.

Divergence of gene function is a hallmark of evolution, but assessing functional divergence over deep time is not trivial. The few alleles available for cross-species studies often fail to expose the entire functional spectrum of genes, potentially obscuring deeply conserved pleiotropic roles. Here, we explore the functional divergence of WUSCHEL HOMEOBOX9 (WOX9), suggested to have species-specific roles in embryo and inflorescence development. Using a cis-regulatory editing drive system, we generate a comprehensive allelic series in tomato, which revealed hidden pleiotropic roles for WOX9. Analysis of accessible chromatin and conserved cis-regulatory sequences identifies the regions responsible for this pleiotropic activity, the functions of which are conserved in groundcherry, a tomato relative. Mimicking these alleles in Arabidopsis, distantly related to tomato and groundcherry, reveals new inflorescence phenotypes, exposing a deeply conserved pleiotropy. We suggest that targeted cis-regulatory mutations can uncover conserved gene functions and reduce undesirable effects in crop improvement.

Genomic Relationships, Novel Loci, and Pleiotropic Mechanisms across Eight Psychiatric Disorders.

Genetic influences on psychiatric disorders transcend diagnostic boundaries, suggesting substantial pleiotropy of contributing loci. However, the nature and mechanisms of these pleiotropic effects remain unclear. We performed analyses of 232,964 cases and 494,162 controls from genome-wide studies of anorexia nervosa, attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depression, obsessive-compulsive disorder, schizophrenia, and Tourette syndrome. Genetic correlation analyses revealed a meaningful structure within the eight disorders, identifying three groups of inter-related disorders. Meta-analysis across these eight disorders detected 109 loci associated with at least two psychiatric disorders, including 23 loci with pleiotropic effects on four or more disorders and 11 loci with antagonistic effects on multiple disorders. The pleiotropic loci are located within genes that show heightened expression in the brain throughout the lifespan, beginning prenatally in the second trimester, and play prominent roles in neurodevelopmental processes. These findings have important implications for psychiatric nosology, drug development, and risk prediction.

Decoupling the Functional Pleiotropy of Stem Cell Factor by Tuning c-Kit Signaling.

Most secreted growth factors and cytokines are functionally pleiotropic because their receptors are expressed on diverse cell types. While important for normal mammalian physiology, pleiotropy limits the efficacy of cytokines and growth factors as therapeutics. Stem cell factor (SCF) is a growth factor that acts through the c-Kit receptor tyrosine kinase to elicit hematopoietic progenitor expansion but can be toxic when administered in vivo because it concurrently activates mast cells. We engineered a mechanism-based SCF partial agonist that impaired c-Kit dimerization, truncating downstream signaling amplitude. This SCF variant elicited biased activation of hematopoietic progenitors over mast cells in vitro and in vivo. Mouse models of SCF-mediated anaphylaxis, radioprotection, and hematopoietic expansion revealed that this SCF partial agonist retained therapeutic efficacy while exhibiting virtually no anaphylactic off-target effects. The approach of biasing cell activation by tuning signaling thresholds and outputs has applications to many dimeric receptor-ligand systems.

Comparative physiological anthropogeny: exploring molecular underpinnings of distinctly human phenotypes.

Anthropogeny is a classic term encompassing transdisciplinary investigations of the origins of the human species. Comparative anthropogeny is a systematic comparison of humans and other living nonhuman hominids (so-called "great apes"), aiming to identify distinctly human features in health and disease, with the overall goal of explaining human origins. We begin with a historical perspective, briefly describing how the field progressed from the earliest evolutionary insights to the current emphasis on in-depth molecular and genomic investigations of "human-specific" biology and an increased appreciation for cultural impacts on human biology. While many such genetic differences between humans and other hominids have been revealed over the last two decades, this information remains insufficient to explain the most distinctive phenotypic traits distinguishing humans from other living hominids. Here we undertake a complementary approach of "comparative physiological anthropogeny," along the lines of the preclinical medical curriculum, i.e., beginning with anatomy and considering each physiological system and in each case considering genetic and molecular components that are relevant. What is ultimately needed is a systematic comparative approach at all levels from molecular to physiological to sociocultural, building networks of related information, drawing inferences, and generating testable hypotheses. The concluding section will touch on distinctive considerations in the study of human evolution, including the importance of gene-culture interactions.

Examining the shared etiology of psychopathology with genome-wide association studies.

Genome-wide association studies (GWASs) have ushered in a new era of reproducible discovery in psychiatric genetics. The field has now identified hundreds of common genetic variants that are associated with mental disorders, and many of them influence more than one disorder. By advancing the understanding of causal biology underlying psychopathology, GWAS results are poised to inform the development of novel therapeutics, stratification of at-risk patients, and perhaps even the revision of top-down classification systems in psychiatry. Here, we provide a concise review of GWAS findings with an emphasis on findings that have elucidated the shared genetic etiology of psychopathology, summarizing insights at three levels of analysis: 1) genome-wide architecture; 2) networks, pathways, and gene sets; and 3) individual variants/genes. Three themes emerge from these efforts. First, all psychiatric phenotypes are heritable, highly polygenic, and influenced by many pleiotropic variants with incomplete penetrance. Second, GWAS results highlight the broad etiological roles of neuronal biology, system-wide effects over localized effects, and early neurodevelopment as a critical period. Third, many loci that are robustly associated with multiple forms of psychopathology harbor genes that are involved in synaptic structure and function. Finally, we conclude our review by discussing the implications that GWAS results hold for the field of psychiatry, as well as expected challenges and future directions in the next stage of psychiatric genetics.

Recent and future perspectives on engineering interferons and other cytokines as therapeutics.

As crucial mediators and regulators of our immune system, cytokines are involved in a broad range of biological processes and are implicated in various disease pathologies. The field of cytokine therapeutics has gained much momentum from the maturation of conventional protein engineering methodologies such as structure-based designs and/or directed evolution, which is further aided by the advent of in silico protein designs and characterization. Just within the past 5 years, there has been an explosion of proof-of-concept, preclinical, and clinical studies that utilize an armory of protein engineering methods to develop cytokine-based drugs. Here, we highlight the key engineering strategies undertaken by recent studies that aim to improve the pharmacodynamic and pharmacokinetic profile of interferons and other cytokines as therapeutics.

A Bayesian approach to Mendelian randomization using summary statistics in the univariable and multivariable settings with correlated pleiotropy.

Mendelian randomization uses genetic variants as instrumental variables to make causal inferences on the effect of an exposure on an outcome. Due to the recent abundance of high-powered genome-wide association studies, many putative causal exposures of interest have large numbers of independent genetic variants with which they associate, each representing a potential instrument for use in a Mendelian randomization analysis. Such polygenic analyses increase the power of the study design to detect causal effects; however, they also increase the potential for bias due to instrument invalidity. Recent attention has been given to dealing with bias caused by correlated pleiotropy, which results from violation of the "instrument strength independent of direct effect" assumption. Although methods have been proposed that can account for this bias, a number of restrictive conditions remain in many commonly used techniques. In this paper, we propose a Bayesian framework for Mendelian randomization that provides valid causal inference under very general settings. We propose the methods MR-Horse and MVMR-Horse, which can be performed without access to individual-level data, using only summary statistics of the type commonly published by genome-wide association studies, and can account for both correlated and uncorrelated pleiotropy. In simulation studies, we show that the approach retains type I error rates below nominal levels even in high-pleiotropy scenarios. We demonstrate the proposed approaches in applied examples in both univariable and multivariable settings, some with very weak instruments.

shaPRS: Leveraging shared genetic effects across traits or ancestries improves accuracy of polygenic scores.

We present shaPRS, a method that leverages widespread pleiotropy between traits or shared genetic effects across ancestries, to improve the accuracy of polygenic scores. The method uses genome-wide summary statistics from two diseases or ancestries to improve the genetic effect estimate and standard error at SNPs where there is homogeneity of effect between the two datasets. When there is significant evidence of heterogeneity, the genetic effect from the disease or population closest to the target population is maintained. We show via simulation and a series of real-world examples that shaPRS substantially enhances the accuracy of polygenic risk scores (PRSs) for complex diseases and greatly improves PRS performance across ancestries. shaPRS is a PRS pre-processing method that is agnostic to the actual PRS generation method, and as a result, it can be integrated into existing PRS generation pipelines and continue to be applied as more performant PRS methods are developed over time.

JASPER: Fast, powerful, multitrait association testing in structured samples gives insight on pleiotropy in gene expression.

Joint association analysis of multiple traits with multiple genetic variants can provide insight into genetic architecture and pleiotropy, improve trait prediction, and increase power for detecting association. Furthermore, some traits are naturally high-dimensional, e.g., images, networks, or longitudinally measured traits. Assessing significance for multitrait genetic association can be challenging, especially when the sample has population sub-structure and/or related individuals. Failure to adequately adjust for sample structure can lead to power loss and inflated type 1 error, and commonly used methods for assessing significance can work poorly with a large number of traits or be computationally slow. We developed JASPER, a fast, powerful, robust method for assessing significance of multitrait association with a set of genetic variants, in samples that have population sub-structure, admixture, and/or relatedness. In simulations, JASPER has higher power, better type 1 error control, and faster computation than existing methods, with the power and speed advantage of JASPER increasing with the number of traits. JASPER is potentially applicable to a wide range of association testing applications, including for multiple disease traits, expression traits, image-derived traits, and microbiome abundances. It allows for covariates, ascertainment, and rare variants and is robust to phenotype model misspecification. We apply JASPER to analyze gene expression in the Framingham Heart Study, where, compared to alternative approaches, JASPER finds more significant associations, including several that indicate pleiotropic effects, most of which replicate previous results, while others have not previously been reported. Our results demonstrate the promise of JASPER for powerful multitrait analysis in structured samples.

The foraging Gene and Its Behavioral Effects: Pleiotropy and Plasticity.

The Drosophila melanogaster foraging (for) gene is a well-established example of a gene with major effects on behavior and natural variation. This gene is best known for underlying the behavioral strategies of rover and sitter foraging larvae, having been mapped and named for this phenotype. Nevertheless, in the last three decades an extensive array of studies describing for's role as a modifier of behavior in a wide range of phenotypes, in both Drosophila and other organisms, has emerged. Furthermore, recent work reveals new insights into the genetic and molecular underpinnings of how for affects these phenotypes. In this article, we discuss the history of the for gene and its role in natural variation in behavior, plasticity, and behavioral pleiotropy, with special attention to recent findings on the molecular structure and transcriptional regulation of this gene.

Pleiotropy of positive selection in ancient ACE2 suggests an alternative hypothesis for bat-specific adaptations to host coronaviruses.

Angiotensin-convertingenzyme 2 (ACE2) has dual functions, regulating cardiovascular physiology and serving as the receptor for coronaviruses. Bats, the only true flying mammals and natural viral reservoirs, have evolved positive alterations in traits related to both functions of ACE2. This suggests significant evolutionary changes in ACE2 during bat evolution. To test this hypothesis, we examine the selection pressure in ACE2 along the ancestral branch of all bats (AncBat-ACE2), where powered flight and bat-coronavirus coevolution occurred, and detect a positive selection signature. To assess the functional effects of positive selection, we resurrect AncBat-ACE2 and its mutant (AncBat-ACE2-mut) created by replacing the positively selected sites. Compared to AncBat-ACE2-mut, AncBat-ACE2 exhibits stronger enzymatic activity, enhances mice's performance in exercise fatigue, and shows lower affinity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our findings indicate the functional pleiotropy of positive selection in the ancient ACE2 of bats, providing an alternative hypothesis for the evolutionary origin of bats' defense against coronaviruses.

A large-effect fitness trade-off across environments is explained by a single mutation affecting cold acclimation.

Identifying the genetic basis of local adaptation and fitness trade-offs across environments is a central goal of evolutionary biology. Cold acclimation is an adaptive plastic response for surviving seasonal freezing, and costs of acclimation may be a general mechanism for fitness trade-offs across environments in temperate zone species. Starting with locally adapted ecotypes of Arabidopsis thaliana from Italy and Sweden, we examined the fitness consequences of a naturally occurring functional polymorphism in CBF2. This gene encodes a transcription factor that is a major regulator of cold-acclimated freezing tolerance and resides within a locus responsible for a genetic trade-off for long-term mean fitness. We estimated the consequences of alternate genotypes of CBF2 on 5-y mean fitness and fitness components at the native field sites by comparing near-isogenic lines with alternate genotypes of CBF2 to their genetic background ecotypes. The effects of CBF2 were validated at the nucleotide level using gene-edited lines in the native genetic backgrounds grown in simulated parental environments. The foreign CBF2 genotype in the local genetic background reduced long-term mean fitness in Sweden by more than 10%, primarily via effects on survival. In Italy, fitness was reduced by more than 20%, primarily via effects on fecundity. At both sites, the effects were temporally variable and much stronger in some years. The gene-edited lines confirmed that CBF2 encodes the causal variant underlying this genetic trade-off. Additionally, we demonstrated a substantial fitness cost of cold acclimation, which has broad implications for potential maladaptive responses to climate change.

Spatiotemporally programming cytokine immunotherapies through protein engineering.

Cytokines have long been considered promising cancer immunotherapy agents due to their endogenous role in activating and proliferating lymphocytes. However, since the initial FDA approvals of Interleukin-2 (IL-2) and Interferon-ɑ (IFNɑ) for oncology over 30 years ago, cytokines have achieved little success in the clinic due to narrow therapeutic windows and dose-limiting toxicities. This is attributable to the discrepancy between the localized, regulated manner in which cytokines are deployed endogenously versus the systemic, untargeted administration used to date in most exogenous cytokine therapies. Furthermore, cytokines' ability to stimulate multiple cell types, often with paradoxical effects, may present significant challenges for their translation into effective therapies. Recently, protein engineering has emerged as a tool to address the shortcomings of first-generation cytokine therapies. In this perspective, we contextualize cytokine engineering strategies such as partial agonism, conditional activation and intratumoral retention through the lens of spatiotemporal regulation. By controlling the time, place, specificity, and duration of cytokine signaling, protein engineering can allow exogenous cytokine therapies to more closely approach their endogenous exposure profile, ultimately moving us closer to unlocking their full therapeutic potential.

Hox, homology, and parsimony: An organismal perspective.

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

A scalable approach to characterize pleiotropy across thousands of human diseases and complex traits using GWAS summary statistics.

Genome-wide association studies (GWASs) across thousands of traits have revealed the pervasive pleiotropy of trait-associated genetic variants. While methods have been proposed to characterize pleiotropic components across groups of phenotypes, scaling these approaches to ultra-large-scale biobanks has been challenging. Here, we propose FactorGo, a scalable variational factor analysis model to identify and characterize pleiotropic components using biobank GWAS summary data. In extensive simulations, we observe that FactorGo outperforms the state-of-the-art (model-free) approach tSVD in capturing latent pleiotropic factors across phenotypes while maintaining a similar computational cost. We apply FactorGo to estimate 100 latent pleiotropic factors from GWAS summary data of 2,483 phenotypes measured in European-ancestry Pan-UK BioBank individuals (N = 420,531). Next, we find that factors from FactorGo are more enriched with relevant tissue-specific annotations than those identified by tSVD (p = 2.58E-10) and validate our approach by recapitulating brain-specific enrichment for BMI and the height-related connection between reproductive system and muscular-skeletal growth. Finally, our analyses suggest shared etiologies between rheumatoid arthritis and periodontal condition in addition to alkaline phosphatase as a candidate prognostic biomarker for prostate cancer. Overall, FactorGo improves our biological understanding of shared etiologies across thousands of GWASs.

Robust multivariable Mendelian randomization based on constrained maximum likelihood.

Mendelian randomization (MR) is a powerful tool for causal inference with observational genome-wide association study (GWAS) summary data. Compared to the more commonly used univariable MR (UVMR), multivariable MR (MVMR) not only is more robust to the notorious problem of genetic (horizontal) pleiotropy but also estimates the direct effect of each exposure on the outcome after accounting for possible mediating effects of other exposures. Despite promising applications, there is a lack of studies on MVMR's theoretical properties and robustness in applications. In this work, we propose an efficient and robust MVMR method based on constrained maximum likelihood (cML), called MVMR-cML, with strong theoretical support. Extensive simulations demonstrate that MVMR-cML performs better than other existing MVMR methods while possessing the above two advantages over its univariable counterpart. An application to several large-scale GWAS summary datasets to infer causal relationships between eight cardiometabolic risk factors and coronary artery disease (CAD) highlights the usefulness and some advantages of the proposed method. For example, after accounting for possible pleiotropic and mediating effects, triglyceride (TG), low-density lipoprotein cholesterol (LDL), and systolic blood pressure (SBP) had direct effects on CAD; in contrast, the effects of high-density lipoprotein cholesterol (HDL), diastolic blood pressure (DBP), and body height diminished after accounting for other risk factors.

Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets.

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

The impact of 22q11.2 copy-number variants on human traits in the general population.

While extensively studied in clinical cohorts, the phenotypic consequences of 22q11.2 copy-number variants (CNVs) in the general population remain understudied. To address this gap, we performed a phenome-wide association scan in 405,324 unrelated UK Biobank (UKBB) participants by using CNV calls from genotyping array. We mapped 236 Human Phenotype Ontology terms linked to any of the 90 genes encompassed by the region to 170 UKBB traits and assessed the association between these traits and the copy-number state of 504 genotyping array probes in the region. We found significant associations for eight continuous and nine binary traits associated under different models (duplication-only, deletion-only, U-shape, and mirror models). The causal effect of the expression level of 22q11.2 genes on associated traits was assessed through transcriptome-wide Mendelian randomization (TWMR), revealing that increased expression of ARVCF increased BMI. Similarly, increased DGCR6 expression causally reduced mean platelet volume, in line with the corresponding CNV effect. Furthermore, cross-trait multivariable Mendelian randomization (MVMR) suggested a predominant role of genuine (horizontal) pleiotropy in the CNV region. Our findings show that within the general population, 22q11.2 CNVs are associated with traits previously linked to genes in the region, and duplications and deletions act upon traits in different fashions. We also showed that gain or loss of distinct segments within 22q11.2 may impact a trait under different association models. Our results have provided new insights to help further the understanding of the complex 22q11.2 region.

Contraception ends the genetic maintenance of human same-sex sexual behavior.

Because human same-sex sexual behavior (SSB) is heritable and leads to fewer offspring, it is puzzling why SSB-associated alleles have not been selectively purged. Current evidence supports the antagonistic pleiotropy hypothesis that SSB-associated alleles benefit individuals exclusively performing opposite-sex sexual behavior by increasing their number of sexual partners and consequently their number of offspring. However, by analyzing the UK Biobank, here, we show that having more sexual partners no longer predicts more offspring since the availability of oral contraceptives in the 1960s and that SSB is now genetically negatively correlated with the number of offspring, suggesting a loss of SSB's genetic maintenance in modern societies.

The paradox of adaptive trait clines with nonclinal patterns in the underlying genes.

Multivariate climate change presents an urgent need to understand how species adapt to complex environments. Population genetic theory predicts that loci under selection will form monotonic allele frequency clines with their selective environment, which has led to the wide use of genotype-environment associations (GEAs). This study used a set of simulations to elucidate the conditions under which allele frequency clines are more or less likely to evolve as multiple quantitative traits adapt to multivariate environments. Phenotypic clines evolved with nonmonotonic (i.e., nonclinal) patterns in allele frequencies under conditions that promoted unique combinations of mutations to achieve the multivariate optimum in different parts of the landscape. Such conditions resulted from interactions among landscape, demography, pleiotropy, and genetic architecture. GEA methods failed to accurately infer the genetic basis of adaptation under a range of scenarios due to first principles (clinal patterns did not evolve) or statistical issues (clinal patterns evolved but were not detected due to overcorrection for structure). Despite the limitations of GEAs, this study shows that a back-transformation of multivariate ordination can accurately predict individual multivariate traits from genotype and environmental data regardless of whether inference from GEAs was accurate. In addition, frameworks are introduced that can be used by empiricists to quantify the importance of clinal alleles in adaptation. This research highlights that multivariate trait prediction from genotype and environmental data can lead to accurate inference regardless of whether the underlying loci display clinal or nonmonotonic patterns.