The many faces of pleiotropy

Biological constraint, evolutionary spandrels and antagonistic pleiotropy

Maximum lifespan differs greatly between species, indicating that the process of senescence is largely genetically determined. Senescence evolves in part due to antagonistic pleiotropy (AP), where selection favors gene variants that increase fitness earlier in life but promote pathology later. Identifying the biological mechanisms by which AP causes senescence is key to understanding the endogenous causes of aging and its attendant diseases. Here we argue that the frequent occurrence of AP as a property of genes reflects the presence of constraint in the biological systems that they specify. This arises particularly because the functionally interconnected nature of biological systems constrains the simultaneous optimization of coupled traits (interconnection constraints), or because individual traits cannot evolve (impossibility constraints). We present an account of aging that integrates AP and biological constraint with recent programmatic aging concepts, including costly programs, quasi-programs, hyperfunction and hypofunction. We argue that AP mechanisms of costly programs and triggered quasi-programs are consequences of constraint, in which costs resulting from hyperfunction or hypofunction cause senescent pathology. Impossibility constraint can also cause hypofunction independently of AP. We also describe how AP corresponds to Stephen Jay Gould's constraint-based concept of evolutionary spandrels, and argue that pathologies arising from AP are bad spandrels. Biological constraint is a conceptual missing link between ultimate and proximate causes of senescence, including diseases of aging.

Common Genetic Influence on the Relationship Between Gaming Addiction and Attention Deficit Hyperactivity Disorder in Young Adults: A Twin Study

Although the relationship between gaming addiction (GA) and attention deficit hyperactivity disorder (ADHD) is well established, the causal mechanism of this relationship remains ambiguous. We aimed to investigate whether common genetic and/or environmental factors explain the GA-ADHD relationship. We recruited 1413 South Korean adult twins (837 monozygotic [MZ], 326 same-sex dizygotic [DZ], and 250 opposite-sex DZ twins; mean age = 23.1 ± 2.8 years) who completed an online survey on GA and related traits. Correlational analysis and bivariate model-fitting analysis were conducted. Phenotypic correlation between GA and ADHD in the present sample was 0.55 (95% CI [0.51, 0.59]). Bivariate model-fitting analysis revealed that genetic variances were 69% (95% CI [64%, 73%]) and 68% (95% CI [63%, 72%]) for ADHD and GA respectively. The remaining variances (ADHD: 31%; GA: 32%) were associated with nonshared environmental variances, including measurement error. Genetic and nonshared environmental correlations between ADHD and GA were 0.68 (95% CI [0.62, 0.74]) and 0.22 (95% CI [0.13, 0.30]) respectively, which indicates that shared genes can explain 82% of the phenotypic correlation between ADHD and GA. Our study demonstrated that the ADHD-GA association was largely due to shared genetic vulnerability.

Understanding developmental system drift

Developmental system drift (DSD) occurs when the genetic basis for homologous traits diverges over time despite conservation of the phenotype. In this Review, we examine the key ideas, evidence and open problems arising from studies of DSD. Recent work suggests that DSD may be pervasive, having been detected across a range of different organisms and developmental processes. Although developmental research remains heavily reliant on model organisms, extrapolation of findings to non-model organisms can be error-prone if the lineages have undergone DSD. We suggest how existing data and modelling approaches may be used to detect DSD and estimate its frequency. More direct study of DSD, we propose, can inform null hypotheses for how much genetic divergence to expect on the basis of phylogenetic distance, while also contributing to principles of gene regulatory evolution.

Epistasis and pleiotropy-induced variation for plant breeding

Epistasis refers to nonallelic interaction between genes that cause bias in estimates of genetic parameters for a phenotype with interactions of two or more genes affecting the same trait. Partitioning of epistatic effects allows true estimation of the genetic parameters affecting phenotypes. Multigenic variation plays a central role in the evolution of complex characteristics, among which pleiotropy, where a single gene affects several phenotypic characters, has a large influence. While pleiotropic interactions provide functional specificity, they increase the challenge of gene discovery and functional analysis. Overcoming pleiotropy-based phenotypic trade-offs offers potential for assisting breeding for complex traits. Modelling higher order nonallelic epistatic interaction, pleiotropy and non-pleiotropy-induced variation, and genotype × environment interaction in genomic selection may provide new paths to increase the productivity and stress tolerance for next generation of crop cultivars. Advances in statistical models, software and algorithm developments, and genomic research have facilitated dissecting the nature and extent of pleiotropy and epistasis. We overview emerging approaches to exploit positive (and avoid negative) epistatic and pleiotropic interactions in a plant breeding context, including developing avenues of artificial intelligence, novel exploitation of large-scale genomics and phenomics data, and involvement of genes with minor effects to analyse epistatic interactions and pleiotropic quantitative trait loci, including missing heritability.

Halofilins as emerging bactofilin families of archaeal cell shape plasticity orchestrators

Bactofilins are rigid, nonpolar bacterial cytoskeletal filaments that link cellular processes to specific curvatures of the cytoplasmic membrane. Although homologs of bactofilins have been identified in archaea and eukaryotes, functional studies have remained confined to bacterial systems. Here, we characterize representatives of two families of archaeal bactofilins from the pleomorphic archaeon Haloferax volcanii, halofilin A (HalA) and halofilin B (HalB). HalA and HalB polymerize in vitro, assembling into straight bundles. HalA polymers are highly dynamic and accumulate at positive membrane curvatures in vivo, whereas HalB forms more static foci that localize in areas of local negative curvatures on the outer cell surface. Gene deletions and live-cell imaging show that halofilins are critical in maintaining morphological integrity during shape transition from disk (sessile) to rod (motile). Morphological defects in ΔhalA result in accumulation of highly positive curvatures in rods but not in disks. Conversely, disk-shaped cells are exclusively affected by halB deletion, resulting in flatter cells. Furthermore, while ΔhalA and ΔhalB cells imprecisely determine the future division plane, defects arise predominantly during the disk-to-rod shape remodeling. The deletion of halA in the haloarchaeon Halobacterium salinarum, whose cells are consistently rod-shaped, impacted morphogenesis but not cell division. Increased levels of halofilins enforced drastic deformations in cells devoid of the S-layer, suggesting that HalB polymers are more stable at defective S-layer lattice regions. Our results suggest that halofilins might play a significant mechanical scaffolding role in addition to possibly directing envelope synthesis.

Distinguishing mutants that resist drugs via different mechanisms by examining fitness tradeoffs

There is growing interest in designing multidrug therapies that leverage tradeoffs to combat resistance. Tradeoffs are common in evolution and occur when, for example, resistance to one drug results in sensitivity to another. Major questions remain about the extent to which tradeoffs are reliable, specifically, whether the mutants that provide resistance to a given drug all suffer similar tradeoffs. This question is difficult because the drug-resistant mutants observed in the clinic, and even those evolved in controlled laboratory settings, are often biased towards those that provide large fitness benefits. Thus, the mutations (and mechanisms) that provide drug resistance may be more diverse than current data suggests. Here, we perform evolution experiments utilizing lineage-tracking to capture a fuller spectrum of mutations that give yeast cells a fitness advantage in fluconazole, a common antifungal drug. We then quantify fitness tradeoffs for each of 774 evolved mutants across 12 environments, finding these mutants group into classes with characteristically different tradeoffs. Their unique tradeoffs may imply that each group of mutants affects fitness through different underlying mechanisms. Some of the groupings we find are surprising. For example, we find some mutants that resist single drugs do not resist their combination, while others do. And some mutants to the same gene have different tradeoffs than others. These findings, on one hand, demonstrate the difficulty in relying on consistent or intuitive tradeoffs when designing multidrug treatments. On the other hand, by demonstrating that hundreds of adaptive mutations can be reduced to a few groups with characteristic tradeoffs, our findings may yet empower multidrug strategies that leverage tradeoffs to combat resistance. More generally speaking, by grouping mutants that likely affect fitness through similar underlying mechanisms, our work guides efforts to map the phenotypic effects of mutation.

Most Pleiotropic Effects of Gene Knockouts Are Evolutionarily Transient in Yeasts

Pleiotropy, the phenomenon in which a single gene influences multiple traits, is a fundamental concept in genetics. However, the evolutionary mechanisms underlying pleiotropy require further investigation. In this study, we conducted parallel gene knockouts targeting 100 transcription factors in 2 strains of Saccharomyces cerevisiae. We systematically examined and quantified the pleiotropic effects of these knockouts on gene expression levels for each transcription factor. Our results showed that the knockout of a single gene generally affected the expression levels of multiple genes in both strains, indicating various degrees of pleiotropic effects. Strikingly, the pleiotropic effects of the knockouts change rapidly between strains in different genetic backgrounds, and ∼85% of them were nonconserved. Further analysis revealed that the conserved effects tended to be functionally associated with the deleted transcription factors, while the nonconserved effects appeared to be more ad hoc responses. In addition, we measured 184 yeast cell morphological traits in these knockouts and found consistent patterns. In order to investigate the evolutionary processes underlying pleiotropy, we examined the pleiotropic effects of standing genetic variations in a population consisting of ∼1,000 hybrid progenies of the 2 strains. We observed that newly evolved expression quantitative trait loci impacted the expression of a greater number of genes than did old expression quantitative trait loci, suggesting that natural selection is gradually eliminating maladaptive or slightly deleterious pleiotropic responses. Overall, our results show that, although being prevalent for new mutations, the majority of pleiotropic effects observed are evolutionarily transient, which explains how evolution proceeds despite complicated pleiotropic effects.

FLT4 causes developmental disorders of the cardiovascular and lymphovascular systems via pleiotropic molecular mechanisms

AIMS

Rare, deleterious genetic variants in FLT4 are associated with Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease. The distinct genetic variants in FLT4 are also an established cause of Milroy disease, the most prevalent form of primary hereditary lymphoedema. The phenotypic features of these two conditions are non-overlapping, implying pleiotropic cellular mechanisms during development.

METHODS AND RESULTS

In this study, we show that FLT4 variants identified in patients with TOF, when expressed in primary human endothelial cells, cause aggregation of FLT4 protein in the perinuclear endoplasmic reticulum, activating proteostatic and metabolic signalling, whereas lymphoedema-associated FLT4 variants and wild-type (WT) FLT4 do not. FLT4 TOF variants display characteristic gene expression profiles in key developmental signalling pathways, revealing a role for FLT4 in cardiogenesis distinct from its role in lymphatic development. Inhibition of proteostatic signalling abrogates these effects, identifying potential avenues for therapeutic intervention. Depletion of flt4 in zebrafish caused cardiac phenotypes of reduced heart size and altered heart looping. These phenotypes were rescued with coinjection of WT human FLT4 mRNA, but incompletely or not at all by mRNA harbouring FLT4 TOF variants.

CONCLUSION

Taken together, we identify a pathogenic mechanism for FLT4 variants predisposing to TOF that is distinct from the known dominant negative mechanism of Milroy-causative variants. FLT4 variants give rise to conditions of the two circulatory subdivisions of the vascular system via distinct developmental pleiotropic molecular mechanisms.

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

Drosophila melanogaster pigmentation demonstrates adaptive phenotypic parallelism but genomic unpredictability over multiple timescales

Populations are capable of responding to environmental change over ecological timescales via adaptive tracking. However, the translation from patterns of allele frequency change to rapid adaptation of complex traits remains unresolved. We used abdominal pigmentation in Drosophila melanogaster as a model phenotype to address the nature, genetic architecture, and repeatability of rapid adaptation in the field. We show that D. melanogaster pigmentation evolves as a highly parallel and deterministic response to shared environmental gradients across latitude and season in natural North American populations. We then experimentally evolved replicate, genetically diverse fly populations in field mesocosms to remove any confounding effects of demography and/or cryptic structure that may drive patterns in wild populations; we show that pigmentation rapidly responds, in parallel, in fewer than ten generations. Thus, pigmentation evolves concordantly in response to spatial and temporal climatic gradients. We next examined whether phenotypic differentiation was associated with allele frequency change at loci with established links to genetic variance in pigmentation in natural populations. We found that across all spatial and temporal scales, phenotypic patterns were associated with variation at pigmentation-related loci, and the sets of genes we identified in each context were largely nonoverlapping. Therefore, our findings suggest that parallel phenotypic evolution is associated with an unpredictable genomic response, with distinct components of the polygenic architecture shifting across each environmental gradient to produce redundant adaptive patterns.

Significance Statement

Shifts in global climate conditions have heightened our need to understand the dynamics and pace of adaptation in natural populations. In order to anticipate the population-level response to rapidly changing environmental conditions, we need to understand whether trait evolution is predictable over short timescales, and whether the genetic basis of adaptation is shared or distinct across multiple timescales. Here, we explored parallelism in the adaptive response of a complex phenotype, D. melanogaster pigmentation, to shared conditions that varied over multiple spatiotemporal scales. Our results demonstrate that while phenotypic adaptation proceeds as a predictable response to environmental gradients, even over short timescales, the genetic basis of the adaptive response is variable and nuanced across spatial and temporal contexts.

Diabetes, glycemic profile and risk of vitiligo: A Mendelian randomization study

BACKGROUD

Previous observational studies have shown that vitiligo usually co-manifests with a variety of dysglycemic diseases, such as Type 1 diabetes mellitus (T1DM) and Type 2 diabetes mellitus (T2DM). Mendelian randomization (MR) analysis was performed to further evaluate the causal association between fasting plasma glucose, glycosylated hemoglobin (HbA1c), T1DM, T2DM and vitiligo.

MATERIALS AND METHODS

We used aggregated genome-wide association data from the Integrative Epidemiology Unit (IEU) online database of European adults vitiligo; HbA1c data were from IEU. Fasting blood glucose data were obtained from the European Bioinformatics Institute (EBI). T1DM and T2DM data were from FinnGen. We used bidirectional two-sample and multivariate MR analyses to test whether dysglycemic measures (fasting blood glucose, HbA1c), diabetes-related measures (T1DM, T2DM) are causatively associated with vitiligo. Inverse variance weighting (IVW) method was used as the main test method, MR-Egger, Weighted mode and Weighted median were used as supplementary methods.

RESULTS

We found no statistically significant evidence to support a causal association between dysglycemic traits and vitiligo, but in the correlation analysis of diabetic traits, our data supported a positive causal association between T1DM and vitiligo (p = 0.018). In the follow-up multivariate MR analysis, our results still supported this conclusion (p = 0.016), and suggested that HbA1c was not a mediator of T1DM affecting the pathogenesis of vitiligo. No reverse causality was found in any of the reverse MR Analyses of dysglycemic traits and diabetic traits.

CONCLUSIONS

Our findings support that T1DM is a risk factor for the development of vitiligo, and this conclusion may explain why the co-presentation of T1DM and vitiligo is often seen in observational studies. Clinical use of measures related to T1DM may be a new idea for the prevention or treatment of vitiligo.

Polygenic selection to a changing optimum under self-fertilisation

Many traits are polygenic, affected by multiple genetic variants throughout the genome. Selection acting on these traits involves co-ordinated allele-frequency changes at these underlying variants, and this process has been extensively studied in random-mating populations. Yet many species self-fertilise to some degree, which incurs changes to genetic diversity, recombination and genome segregation. These factors cumulatively influence how polygenic selection is realised in nature. Here, we use analytical modelling and stochastic simulations to investigate to what extent self-fertilisation affects polygenic adaptation to a new environment. Our analytical solutions show that while selfing can increase adaptation to an optimum, it incurs linkage disequilibrium that can slow down the initial spread of favoured mutations due to selection interference, and favours the fixation of alleles with opposing trait effects. Simulations show that while selection interference is present, high levels of selfing (at least 90%) aids adaptation to a new optimum, showing a higher long-term fitness. If mutations are pleiotropic then only a few major-effect variants fix along with many neutral hitchhikers, with a transient increase in linkage disequilibrium. These results show potential advantages to self-fertilisation when adapting to a new environment, and how the mating system affects the genetic composition of polygenic selection.

Parallel evolution of integrated craniofacial traits in trophic specialist pupfishes

Populations may adapt to similar environments via parallel or non-parallel genetic changes, but the frequency of these alternative mechanisms and underlying contributing factors are still poorly understood outside model systems. We used QTL mapping to investigate the genetic basis of highly divergent craniofacial traits between the scale-eater (Cyprinodon desquamator) and molluscivore (C. brontotheroides) pupfish adapting to two different hypersaline lake environments on San Salvador Island, Bahamas. We lab-reared F2 scale-eater x molluscivore intercrosses from two different lake populations, estimated linkage maps, scanned for significant QTL for 29 skeletal and craniofacial traits, female mate preference, and sex. We compared the location of QTL between lakes to quantify parallel and non-parallel genetic changes. We detected significant QTL for six craniofacial traits in at least one lake. However, nearly all shared QTL loci were associated with a different craniofacial trait within each lake. Therefore, our estimate of parallel evolution of craniofacial genetic architecture could range from one out of six identical trait QTL (low parallelism) to five out of six integrated trait QTL (high parallelism). We suggest that pleiotropy and trait integration can affect estimates of parallel evolution, particularly within rapid radiations. We also observed increased adaptive introgression in shared QTL regions, suggesting that gene flow contributed to parallel evolution. Overall, our results suggest that the same genomic regions may contribute to parallel adaptation across integrated suites of craniofacial traits, rather than specific traits, and highlight the need for a more expansive definition of parallel evolution.

Distinguishing mutants that resist drugs via different mechanisms by examining fitness tradeoffs

There is growing interest in designing multidrug therapies that leverage tradeoffs to combat resistance. Tradeoffs are common in evolution and occur when, for example, resistance to one drug results in sensitivity to another. Major questions remain about the extent to which tradeoffs are reliable, specifically, whether the mutants that provide resistance to a given drug all suffer similar tradeoffs. This question is difficult because the drug-resistant mutants observed in the clinic, and even those evolved in controlled laboratory settings, are often biased towards those that provide large fitness benefits. Thus, the mutations (and mechanisms) that provide drug resistance may be more diverse than current data suggests. Here, we perform evolution experiments utilizing lineage-tracking to capture a fuller spectrum of mutations that give yeast cells a fitness advantage in fluconazole, a common antifungal drug. We then quantify fitness tradeoffs for each of 774 evolved mutants across 12 environments, finding these mutants group into 6 classes with characteristically different tradeoffs. Their unique tradeoffs may imply that each group of mutants affects fitness through different underlying mechanisms. Some of the groupings we find are surprising. For example, we find some mutants that resist single drugs do not resist their combination, while others do. And some mutants to the same gene have different tradeoffs than others. These findings, on one hand, demonstrate the difficulty in relying on consistent or intuitive tradeoffs when designing multidrug treatments. On the other hand, by demonstrating that hundreds of adaptive mutations can be reduced to a few groups with characteristic tradeoffs, our findings may yet empower multidrug strategies that leverage tradeoffs to combat resistance. More generally speaking, by grouping mutants that likely affect fitness through similar underlying mechanisms, our work guides efforts to map the phenotypic effects of mutation.

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits

BACKGROUND

Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS).

RESULTS

In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits.

CONCLUSIONS

The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

Genomic regions, candidate genes, and pleiotropic variants associated with physiological and anatomical indicators of heat stress response in lactating sows

BACKGROUND

Heat stress (HS) poses significant threats to the sustainability of livestock production. Genetically improving heat tolerance could enhance animal welfare and minimize production losses during HS events. Measuring phenotypic indicators of HS response and understanding their genetic background are crucial steps to optimize breeding schemes for improved climatic resilience. The identification of genomic regions and candidate genes influencing the traits of interest, including variants with pleiotropic effects, enables the refinement of genotyping panels used to perform genomic prediction of breeding values and contributes to unraveling the biological mechanisms influencing heat stress response. Therefore, the main objectives of this study were to identify genomic regions, candidate genes, and potential pleiotropic variants significantly associated with indicators of HS response in lactating sows using imputed whole-genome sequence (WGS) data. Phenotypic records for 18 traits and genomic information from 1,645 lactating sows were available for the study. The genotypes from the PorcineSNP50K panel containing 50,703 single nucleotide polymorphisms (SNPs) were imputed to WGS and after quality control, 1,622 animals and 7,065,922 SNPs were included in the analyses.

RESULTS

A total of 1,388 unique SNPs located on sixteen chromosomes were found to be associated with 11 traits. Twenty gene ontology terms and 11 biological pathways were shown to be associated with variability in ear skin temperature, shoulder skin temperature, rump skin temperature, tail skin temperature, respiration rate, panting score, vaginal temperature automatically measured every 10 min, vaginal temperature measured at 0800 h, hair density score, body condition score, and ear area. Seven, five, six, two, seven, 15, and 14 genes with potential pleiotropic effects were identified for indicators of skin temperature, vaginal temperature, animal temperature, respiration rate, thermoregulatory traits, anatomical traits, and all traits, respectively.

CONCLUSIONS

Physiological and anatomical indicators of HS response in lactating sows are heritable but highly polygenic. The candidate genes found are associated with important gene ontology terms and biological pathways related to heat shock protein activities, immune response, and cellular oxidative stress. Many of the candidate genes with pleiotropic effects are involved in catalytic activities to reduce cell damage from oxidative stress and cellular mechanisms related to immune response.

Comparative Analysis of Maternal Gene Expression Patterns Unravels Evolutionary Signatures Across Reproductive Modes

Maternal genes have a pivotal role in regulating metazoan early development. As such their functions have been extensively studied since the dawn of developmental biology. The temporal and spatial dynamics of their transcripts have been thoroughly described in model organisms and their functions have been undergoing heavy investigations. Yet, less is known about the evolutionary changes shaping their presence within diverse oocytes. Due to their unique maternal inheritance pattern, a high degree is predicted to be present when it comes to their expression. Insofar only limited and conflicting results have emerged around it. Here, we set out to elucidate which evolutionary changes could be detected in the maternal gene expression patterns using phylogenetic comparative methods on RNAseq data from 43 species. Using normalized gene expression values and fold change information throughout early development we set out to find the best-fitting evolutionary model. Through modeling, we find evidence supporting both the high degree of divergence and constraint on gene expression values, together with their temporal dynamics. Furthermore, we find that maternal gene expression alone can be used to explain the reproductive modes of different species. Together, these results suggest a highly dynamic evolutionary landscape of maternal gene expression. We also propose a possible functional dichotomy of maternal genes which is influenced by the reproductive strategy undertaken by examined species.

Rapid evolutionary change in trait correlations of single proteins

Many organismal traits are genetically determined and covary in evolving populations. The resulting trait correlations can either help or hinder evolvability - the ability to bring forth new and adaptive phenotypes. The evolution of evolvability requires that trait correlations themselves must be able to evolve, but we know little about this ability. To learn more about it, we here study two evolvable systems, a yellow fluorescent protein and the antibiotic resistance protein VIM-2 metallo beta-lactamase. We consider two traits in the fluorescent protein, namely the ability to emit yellow and green light, and three traits in our enzyme, namely the resistance against ampicillin, cefotaxime, and meropenem. We show that correlations between these traits can evolve rapidly through both mutation and selection on short evolutionary time scales. In addition, we show that these correlations are driven by a protein's ability to fold, because single mutations that alter foldability can dramatically change trait correlations. Since foldability is important for most proteins and their traits, mutations affecting protein folding may alter trait correlations mediated by many other proteins. Thus, mutations that affect protein foldability may also help shape the correlations of complex traits that are affected by hundreds of proteins.

Pleiotropic effects of melanin pigmentation: haemoparasite infection intensity but not telomere length is associated with plumage morph in black sparrowhawks

There is increasing recognition of the potential pleiotropic effects of melanin pigmentation, particularly on immunity, with reports of variation in haemoparasite infection intensity and immune responses between the morphs of colour-polymorphic bird species. In a population of the black sparrowhawk (Accipiter melanoleucus) in western South Africa, light morphs have a higher haemoparasite infection intensity, but no physiological effects of this are apparent. Here, we investigate the possible effects of haemoparasite infection on telomere length in this species and explore whether relative telomere length is associated with either plumage morph or sex. Using quantitative polymerase chain reaction analysis, we confirmed that dark morphs had a lower haemoparasite infection intensity than light morphs. However, we found no differences in telomere length associated with either the haemoparasite infection status or morph in adults, although males have longer telomeres than females. While differences in haemoparasite intensity between morphs are consistent with pleiotropic effects of melanin pigmentation in the black sparrowhawk, we found no evidence that telomere length was associated with haemoparasite infection. Further work is needed to investigate the implications of possible pleiotropic effects of plumage morph and their potential role in the maintenance of colour polymorphism in this species.

Decoding and reconstructing disease relations between dry eye and depression: a multimodal investigation comprising meta-analysis, genetic pathways and Mendelian randomization

INTRODUCTION

The clinical presentations of dry eye disease (DED) and depression (DEP) often comanifest. However, the robustness and the mechanisms underlying this association were undetermined.

OBJECTIVES

To this end, we set up a three-segment study that employed multimodality results (meta-analysis, genome-wide association study [GWAS] and Mendelian randomization [MR]) to elucidate the association, common pathways and causality between DED and DEP.

METHODS

A meta-analysis comprising 26 case-control studies was first conducted to confirm the DED-DEP association. Next, we performed a linkage disequilibrium (LD)-adjusted GWAS and targeted phenotype association study (PheWAS) in East Asian TW Biobank (TWB) and European UK Biobank (UKB) populations. Single-nucleotide polymorphisms (SNPs) were further screened for molecular interactions and common pathways at the functional gene level. To further elucidate the activated pathways in DED and DEP, a systemic transcriptome review was conducted on RNA sequencing samples from the Gene Expression Omnibus. Finally, 48 MR experiments were implemented to examine the bidirectional causation between DED and DEP.

RESULTS

Our meta-analysis showed that DED patients are associated with an increased DEP prevalence (OR = 1.83), while DEP patients have a concurrent higher risk of DED (OR = 2.34). Notably, cross-disease GWAS analysis revealed that similar genetic architecture (rG = 0.19) and pleiotropic functional genes contributed to phenotypes in both diseases. Through protein-protein interaction and ontology convergence, we summarized the pleiotropic functional genes under the ontology of immune activation, which was further validated by a transcriptome systemic review. Importantly, the inverse variance-weighted (IVW)-MR experiments in both TWB and UKB populations (p value <0.001) supported the bidirectional exposure-outcome causation for DED-to-DEP and DEP-to-DED. Despite stringent LD-corrected instrumental variable re-selection, the bidirectional causation between DED and DEP remained.

CONCLUSION

With the multi-modal evidence combined, we consolidated the association and causation between DED and DEP.

Reducing Epistasis and Pleiotropy Can Avoid the Survival of the Flattest Tragedy

This study investigates whether reducing epistasis and pleiotropy enhances mutational robustness in evolutionary adaptation, utilizing an indirect encoded model within the "survival of the flattest" (SoF) fitness landscape. By simulating genetic variations and their phenotypic consequences, we explore organisms' adaptive mechanisms to maintain positions on higher, narrower evolutionary peaks amidst environmental and genetic pressures. Our results reveal that organisms can indeed sustain their advantageous positions by minimizing the complexity of genetic interactions-specifically, by reducing the levels of epistasis and pleiotropy. This finding suggests a counterintuitive strategy for evolutionary stability: simpler genetic architectures, characterized by fewer gene interactions and multifunctional genes, confer a survival advantage by enhancing mutational robustness. This study contributes to our understanding of the genetic underpinnings of adaptability and robustness, challenging traditional views that equate complexity with fitness in dynamic environments.

Competition Increases Risk of Species Extinction during Extreme Warming

AbstractTemperature and interspecific competition are fundamental drivers of community structure in natural systems and can interact to affect many measures of species performance. However, surprisingly little is known about the extent to which competition affects extinction temperatures during extreme warming. This information is important for evaluating future threats to species from extreme high-temperature events and heat waves, which are rising in frequency and severity around the world. Using experimental freshwater communities of rotifers and ciliates, this study shows that interspecific competition can lower the threshold temperature at which local extinction occurs, reducing time to extinction during periods of sustained warming by as much as 2 weeks. Competitors may lower extinction temperatures by altering biochemical characteristics of the natural environment that affect temperature tolerance (e.g., levels of dissolved oxygen, nutrients, and metabolic wastes) or by accelerating population decline through traditional effects of resource depletion on life history parameters that affect population growth rates. The results suggest that changes in community structure in space and time could drive variability in upper thermal limits.

Genetic influences on depression and selection into adverse life experiences

Genome-wide association studies find that a large number of genetic variants jointly influence the risk of depression, which is summarized by polygenic indices (PGIs) of depressive symptoms and major depression. But PGIs by design remain agnostic about the causal mechanisms linking genes to depression. Meanwhile, the role of adverse life experiences in shaping depression risk is well-documented, including via gene-environment correlation. Building on theoretical work on dynamic and contingent genetic selection, we suggest that genetic influences may lead to differential selection into negative life experiences, forging gene-environment correlations that manifest in various permutations of depressive behaviors and environmental adversities. We also examine the extent to which apparent genetic influences may reflect spurious associations due to factors such as indirect genetic effects. Using data from two large surveys of middle-aged and older US adults, we investigate to what extent a PGI of depression predicts the risk of 27 different adversities. Further, to glean insights about the kinds of processes that might lead to gene-environment correlation, we augment these analyses with data from an original preregistered survey to measure cultural understandings of the behavioral dependence of various adversities. We find that the PGI predicts the risk of majority of adversities, net of class background and prior depression, and that the selection risk is greater for adversities typically perceived as being dependent on peoples' own behaviors. Taken together, our findings suggest that the PGI of depression largely picks up the risk of behaviorally-influenced adversities, but to a lesser degree also captures other environmental influences. The results invite further exploration into the behavioral and interactional processes that lie along the pathways intervening between genetic differences and wellbeing.

Exploring the costs of phenotypic plasticity for evolvable digital organisms

Phenotypic plasticity is usually defined as a property of individual genotypes to produce different phenotypes when exposed to different environmental conditions. While the benefits of plasticity for adaptation are well established, the costs associated with plasticity remain somewhat obscure. Understanding both why and how these costs arise could help us explain and predict the behavior of living creatures as well as allow the design of more adaptable robotic systems. One of the challenges of conducting such investigations concerns the difficulty of isolating the effects of different types of costs and the lack of control over environmental conditions. The present study addresses these challenges by using virtual worlds (software) to investigate the environmentally regulated phenotypic plasticity of digital organisms. The experimental setup guarantees that potential genetic costs of plasticity are isolated from other plasticity-related costs. Multiple populations of organisms endowed with and without phenotypic plasticity in either the body or the brain are evolved in simulation, and organisms must cope with different environmental conditions. The traits and fitness of the emergent organisms are compared, demonstrating cases in which plasticity is beneficial and cases in which it is neutral. The hypothesis put forward here is that the potential benefits of plasticity might be undermined by the genetic costs related to plasticity itself. The results suggest that this hypothesis is true, while further research is needed to guarantee that the observed effects unequivocally derive from genetic costs and not from some other (unforeseen) mechanism related to plasticity.

Causal relationship between asthma and ulcerative colitis and the mediating role of interleukin-18: a bidirectional Mendelian study and mediation analysis

Objective

Numerous observational investigations have documented a correlation between asthma and ulcerative colitis(UC). In this Mendelian Randomization (MR) study, we utilized extensive summary data from Genome-Wide Association Studies (GWAS) to further estimate the association between adult-onset asthma and the risk of UC, and to investigate the role of Interleukin-18 (IL-18) as a potential mediator.

Materials and methods

A two-step, two-sample MR study was conducted through mediation analysis. For this study, we employed a two-sample MR analysis using the inverse variance-weighted (IVW), weighted median, weighted mode, and MR-Egger regression techniques. We utilized publicly accessible summary statistics from a GWAS meta-analysis of adult-onset asthma in the UK Biobank (n=327,253; cases=26,582; controls=300,671) as the exposure factor. The outcomes were derived from GWAS data of individuals with European ancestry (n=26,405; cases=6,687; controls=19,718). GWAS data for IL-18 were obtained from individuals of European ancestry (n=9,785,222; cases=3,636; controls=9,781,586).

Results

The MR analysis indicates that adult-onset asthma is associated with an increased risk of UC, with an odds ratio (OR) of 1.019 (95% CI 1.001-1.045, P=0.006). However, there is no strong evidence to suggest that UC significantly impacts the risk of adult-onset asthma. IL-18 may act as a potential mediator in the causal relationship between adult-onset asthma and UC, with a mediation proportion of 3.9% (95% CI, 0.6%-6.9%).

Conclusion

In summary, our study established a causal relationship between asthma and UC, in which IL-18 contributes to a small extent. However, the primary factors underlying the influence of asthma on UC remain unclear. Future research should focus on identifying other potential mediators. In clinical practice, it is important to pay greater attention to intestinal lesions in patients with asthma.

Natural variation in the SVP contributes to the pleiotropic adaption of Arabidopsis thaliana across contrasted habitats

Coordinated plant adaptation involves the interplay of multiple traits driven by habitat-specific selection pressures. Pleiotropic effects, wherein genetic variants of a single gene control multiple traits, can expedite such adaptations. Until present, only a limited number of genes have been reported to exhibit pleiotropy. Here, we create a recombinant inbred line (RIL) population derived from two Arabidopsis thaliana (A. thaliana) ecotypes originating from divergent habitats. Using this RIL population, we identify an allelic variation in a MADS-box transcription factor, SHORT VEGETATIVE PHASE (SVP), which exerts a pleiotropic effect on leaf size and drought-versus-humidity tolerance. Further investigation reveals that a natural null variant of the SVP protein disrupts its normal regulatory interactions with target genes, including GRF3, CYP707A1/3, and AtBG1, leading to increased leaf size, enhanced tolerance to humid conditions, and changes in flowering time of humid conditions in A. thaliana. Remarkably, polymorphic variations in this gene have been traced back to early A. thaliana populations, providing a genetic foundation and plasticity for subsequent colonization of diverse habitats by influencing multiple traits. These findings advance our understanding of how plants rapidly adapt to changing environments by virtue of the pleiotropic effects of individual genes on multiple trait alterations.

Casual association between childhood body mass index and risk of psoriasis: A Mendelian randomization study

BACKGROUND

Observational studies have suggested that childhood body mass index (BMI) is associated with the risk of psoriasis. However, their causal relationship remains unclear. In this investigation, we aimed to determine whether an association exists between childhood BMI and psoriasis.

METHODS

Using summary statistics for childhood BMI of European descent from publicly available GWAS meta-analyses (n = 39 620), we conducted Mendelian randomization (MR) research using the inverse variance weighting (IVW), weighted median, and MR-Egger regression techniques. The outcome was a genome-wide association studies (GWAS) for the self-reported non-cancer disease classification psoriasis in the UK Biobank population (total n = 337 159; case = 3871; control = 333 288).

RESULTS

We selected instrumental variables from 16 single-molecule polymorphisms that attained genome-wide significance in GWAS on childhood BMI. Using the IVW method, our findings supported a causal relationship between childhood BMI and psoriasis (beta = 0.003, standard error [SE] = 0.001, p = 0.006). Using MR-Egger regression analysis, we evaluated the potential for directional pleiotropy to bias our results (intercept = 0.00039, p-value = 0.247) and found no causal relationship between childhood BMI and psoriasis (beta = -0.002, SE = 0.004, p = 0.625). The weighted median method, however, provided proof of a causal relationship (beta = 0.003, SE = 0.001, p = 0.029). Cochran's Q test and the funnel plot revealed little proof of heterogeneity or asymmetry, indicating the lack of directional pleiotropy.

CONCLUSION

According to the findings of the MR analysis, an increased childhood BMI may be linked to a higher likelihood of psoriasis.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

Single-cell, whole-embryo phenotyping of mammalian developmental disorders

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

Patterns and evolutionary consequences of pleiotropy

Pleiotropy refers to the phenomenon of one gene or one mutation affecting multiple phenotypic traits. While the concept of pleiotropy is as old as Mendelian genetics, functional genomics has finally allowed the first glimpses of the extent of pleiotropy for a large fraction of genes in a genome. After describing conceptual and operational difficulties in quantifying pleiotropy and the pros and cons of various methods for measuring pleiotropy, I review empirical data on pleiotropy, which generally show an L-shaped distribution of the degree of pleiotropy (i.e., the number of traits affected) with most genes having low pleiotropy. I then review the current understanding of the molecular basis of pleiotropy. The rest of the review discusses evolutionary consequences of pleiotropy, focusing on advances in topics including the cost of complexity, regulatory vs. coding evolution, environmental pleiotropy and adaptation, evolution of ageing and other seemingly harmful traits, and evolutionary resolution of pleiotropy.

Wnt9 directs zebrafish heart tube assembly via a combination of canonical and non-canonical pathway signaling

During zebrafish heart formation, cardiac progenitor cells converge at the embryonic midline where they form the cardiac cone. Subsequently, this structure transforms into a heart tube. Little is known about the molecular mechanisms that control these morphogenetic processes. Here, we use light-sheet microscopy and combine genetic, molecular biological and pharmacological tools to show that the paralogous genes wnt9a/b are required for the assembly of the nascent heart tube. In wnt9a/b double mutants, cardiomyocyte progenitor cells are delayed in their convergence towards the embryonic midline, the formation of the heart cone is impaired and the transformation into an elongated heart tube fails. The same cardiac phenotype occurs when both canonical and non-canonical Wnt signaling pathways are simultaneously blocked by pharmacological inhibition. This demonstrates that Wnt9a/b and canonical and non-canonical Wnt signaling regulate the migration of cardiomyocyte progenitor cells and control the formation of the cardiac tube. This can be partly attributed to their regulation of the timing of cardiac progenitor cell differentiation. Our study demonstrates how these morphogens activate a combination of downstream pathways to direct cardiac morphogenesis.

Inferring sparse structure in genotype-phenotype maps

Correlation among multiple phenotypes across related individuals may reflect some pattern of shared genetic architecture: individual genetic loci affect multiple phenotypes (an effect known as pleiotropy), creating observable relationships between phenotypes. A natural hypothesis is that pleiotropic effects reflect a relatively small set of common "core" cellular processes: each genetic locus affects one or a few core processes, and these core processes in turn determine the observed phenotypes. Here, we propose a method to infer such structure in genotype-phenotype data. Our approach, sparse structure discovery (SSD) is based on a penalized matrix decomposition designed to identify latent structure that is low-dimensional (many fewer core processes than phenotypes and genetic loci), locus-sparse (each locus affects few core processes), and/or phenotype-sparse (each phenotype is influenced by few core processes). Our use of sparsity as a guide in the matrix decomposition is motivated by the results of a novel empirical test indicating evidence of sparse structure in several recent genotype-phenotype datasets. First, we use synthetic data to show that our SSD approach can accurately recover core processes if each genetic locus affects few core processes or if each phenotype is affected by few core processes. Next, we apply the method to three datasets spanning adaptive mutations in yeast, genotoxin robustness assay in human cell lines, and genetic loci identified from a yeast cross, and evaluate the biological plausibility of the core process identified. More generally, we propose sparsity as a guiding prior for resolving latent structure in empirical genotype-phenotype maps.

Trade-off between local replication and long-distance dissemination during experimental evolution of a satellite RNA

Satellite RNAs (satRNAs) are molecular parasites that depend on their non-homologous helper viruses (HVs) for essential biological functions. While there are multiple molecular and phylogenetic studies on satRNAs, there is no experimental evolution study on how satRNAs may evolve in common infection conditions. In this study, we serially passaged the Bamboo mosaic virus (BaMV) associated-satRNA (satBaMV) under conditions in which satBaMV either coinfects an uninfected host plant, Nicotiana benthamiana, with BaMV or superinfects a transgenic N. benthamiana expressing the full-length BaMV genome. Single-nucleotide polymorphisms (SNPs) of satBaMV populations were analyzed by deep sequencing. Forty-eight SNPs were identified across four different experimental treatments. Most SNPs are treatment-specific, and some are also ephemeral. However, mutations at positions 30, 34, 63, and 82, all located at the 5' untranslated region (UTR), are universal in all treatments. These universal SNPs are configured into several haplotypes and follow different population dynamics. We constructed isogenic satBaMV strains only differing at positions 30 and 82 and conducted competition experiments in protoplasts and host plants. We found that the haplotype that reached high frequency in protoplasts and inoculation leaves also exhibited poor dissemination to systemic leaves and vice versa, thus suggesting an apparent trade-off between local replication and long-distance dissemination. We posit that the trade-off is likely caused by antagonistic pleiotropy at the 5' UTR. Our findings revealed a hitherto under-explored connection between satRNA genome replication and movement within a host plant. The significance of such a connection during satRNA evolution warrants a more thorough investigation.

Pleiotropy drives evolutionary repair of the responsiveness of polarized cell growth to environmental cues

The ability of cells to translate different extracellular cues into different intracellular responses is vital for their survival in unpredictable environments. In Saccharomyces cerevisiae, cell polarity is modulated in response to environmental signals which allows cells to adopt varying morphologies in different external conditions. The responsiveness of cell polarity to extracellular cues depends on the integration of the molecular network that regulates polarity establishment with networks that signal environmental changes. The coupling of molecular networks often leads to pleiotropic interactions that can make it difficult to determine whether the ability to respond to external signals emerges as an evolutionary response to environmental challenges or as a result of pleiotropic interactions between traits. Here, we study how the propensity of the polarity network of S. cerevisiae to evolve toward a state that is responsive to extracellular cues depends on the complexity of the environment. We show that the deletion of two genes, BEM3 and NRP1, disrupts the ability of the polarity network to respond to cues that signal the onset of the diauxic shift. By combining experimental evolution with whole-genome sequencing, we find that the restoration of the responsiveness to these cues correlates with mutations in genes involved in the sphingolipid synthesis pathway and that these mutations frequently settle in evolving populations irrespective of the complexity of the selective environment. We conclude that pleiotropic interactions make a significant contribution to the evolution of networks that are responsive to extracellular cues.

A genome-wide search for pleiotropy in more than 100,000 harmonized longitudinal cognitive domain scores

BACKGROUND

More than 75 common variant loci account for only a portion of the heritability for Alzheimer's disease (AD). A more complete understanding of the genetic basis of AD can be deduced by exploring associations with AD-related endophenotypes.

METHODS

We conducted genome-wide scans for cognitive domain performance using harmonized and co-calibrated scores derived by confirmatory factor analyses for executive function, language, and memory. We analyzed 103,796 longitudinal observations from 23,066 members of community-based (FHS, ACT, and ROSMAP) and clinic-based (ADRCs and ADNI) cohorts using generalized linear mixed models including terms for SNP, age, SNP × age interaction, sex, education, and five ancestry principal components. Significance was determined based on a joint test of the SNP's main effect and interaction with age. Results across datasets were combined using inverse-variance meta-analysis. Genome-wide tests of pleiotropy for each domain pair as the outcome were performed using PLACO software.

RESULTS

Individual domain and pleiotropy analyses revealed genome-wide significant (GWS) associations with five established loci for AD and AD-related disorders (BIN1, CR1, GRN, MS4A6A, and APOE) and eight novel loci. ULK2 was associated with executive function in the community-based cohorts (rs157405, P = 2.19 × 10-9). GWS associations for language were identified with CDK14 in the clinic-based cohorts (rs705353, P = 1.73 × 10-8) and LINC02712 in the total sample (rs145012974, P = 3.66 × 10-8). GRN (rs5848, P = 4.21 × 10-8) and PURG (rs117523305, P = 1.73 × 10-8) were associated with memory in the total and community-based cohorts, respectively. GWS pleiotropy was observed for language and memory with LOC107984373 (rs73005629, P = 3.12 × 10-8) in the clinic-based cohorts, and with NCALD (rs56162098, P = 1.23 × 10-9) and PTPRD (rs145989094, P = 8.34 × 10-9) in the community-based cohorts. GWS pleiotropy was also found for executive function and memory with OSGIN1 (rs12447050, P = 4.09 × 10-8) and PTPRD (rs145989094, P = 3.85 × 10-8) in the community-based cohorts. Functional studies have previously linked AD to ULK2, NCALD, and PTPRD.

CONCLUSION

Our results provide some insight into biological pathways underlying processes leading to domain-specific cognitive impairment and AD, as well as a conduit toward a syndrome-specific precision medicine approach to AD. Increasing the number of participants with harmonized cognitive domain scores will enhance the discovery of additional genetic factors of cognitive decline leading to AD and related dementias.

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.

Evolutionary constraints mediate extinction risk under climate change

Mounting evidence suggests that rapid evolutionary adaptation may rescue some organisms from the impacts of climate change. However, evolutionary constraints might hinder this process, especially when different aspects of environmental change generate antagonistic selection on genetically correlated traits. Here, we use individual-based simulations to explore how genetic correlations underlying the thermal physiology of ectotherms might influence their responses to the two major components of climate change-increases in mean temperature and thermal variability. We found that genetic correlations can influence population dynamics under climate change, with declines in population size varying three-fold depending on the type of correlation present. Surprisingly, populations whose thermal performance curves were constrained by genetic correlations often declined less rapidly than unconstrained populations. Our results suggest that accurate forecasts of the impact of climate change on ectotherms will require an understanding of the genetic architecture of the traits under selection.

Expression inheritance and constraints on cis- and trans-regulatory mutations underlying lotus color variation

Both cis- and trans-regulatory mutations drive changes in gene expression that underpin plant phenotypic evolution. However, how and why these two major types of regulatory mutations arise in different genes and how gene expression is inherited and associated with these regulatory changes are unclear. Here, by studying allele-specific expression in F1 hybrids of pink-flowered sacred lotus (Nelumbo nucifera) and yellow-flowered American lotus (N. lutea), we reveal the relative contributions of cis- and trans-regulatory changes to interspecific expression rewiring underlying petal color change and how the expression is inherited in hybrids. Although cis-only variants influenced slightly more genes, trans-only variants had a stronger impact on expression differences between species. In F1 hybrids, genes under cis-only and trans-only regulatory effects showed a propensity toward additive and dominant inheritance, respectively, whereas transgressive inheritance was observed in genes carrying both cis- and trans-variants acting in opposite directions. By investigating anthocyanin and carotenoid coexpression networks in petals, we found that the same category of regulatory mutations, particularly trans-variants, tend to rewire hub genes in coexpression modules underpinning flower color differentiation between species; we identified 45 known genes with cis- and trans-regulatory variants significantly correlated with flower coloration, such as ANTHOCYANIN 5-AROMATIC ACYLTRANSFERASE (ACT), GLUTATHIONE S-TRANSFERASE F11 (GSTF11), and LYCOPENE Ε-CYCLASE (LCYE). Notably, the relative abundance of genes in different categories of regulatory divergence was associated with the inferred magnitude of constraints like expression level and breadth. Overall, our study suggests distinct selective constraints and modes of gene expression inheritance among different regulatory mutations underlying lotus petal color divergence.

Elucidating the patterns of pleiotropy and its biological relevance in maize

Pleiotropy-when a single gene controls two or more seemingly unrelated traits-has been shown to impact genes with effects on flowering time, leaf architecture, and inflorescence morphology in maize. However, the genome-wide impact of biological pleiotropy across all maize phenotypes is largely unknown. Here, we investigate the extent to which biological pleiotropy impacts phenotypes within maize using GWAS summary statistics reanalyzed from previously published metabolite, field, and expression phenotypes across the Nested Association Mapping population and Goodman Association Panel. Through phenotypic saturation of 120,597 traits, we obtain over 480 million significant quantitative trait nucleotides. We estimate that only 1.56-32.3% of intervals show some degree of pleiotropy. We then assess the relationship between pleiotropy and various biological features such as gene expression, chromatin accessibility, sequence conservation, and enrichment for gene ontology terms. We find very little relationship between pleiotropy and these variables when compared to permuted pleiotropy. We hypothesize that biological pleiotropy of common alleles is not widespread in maize and is highly impacted by nuisance terms such as population structure and linkage disequilibrium. Natural selection on large standing natural variation in maize populations may target wide and large effect variants, leaving the prevalence of detectable pleiotropy relatively low.

Bidirectional Causal Associations Between Same-Sex Attraction and Psychological Distress: Testing Moderation and Mediation Effects

Only one study has examined bidirectional causality between sexual minority status (having same-sex attraction) and psychological distress. We combined twin and genomic data from 8700 to 9700 participants in the UK Twins Early Development Study cohort at ≈21 years to replicate and extend these bidirectional causal effects using separate unidirectional Mendelian Randomization-Direction of Causation models. We further modified these models to separately investigate sex differences, moderation by childhood factors (retrospectively-assessed early-life adversity and prospectively-assessed childhood gender nonconformity), and mediation by victimization. All analyses were carried out in OpenMx in R. Same-sex attraction causally influenced psychological distress with significant reverse causation (beta = 0.19 and 0.17; 95% CIs = 0.09, 0.29 and 0.08, 0.25 respectively) and no significant sex differences. The same-sex attraction → psychological distress causal path was partly mediated by victimization (12.5%) while the reverse causal path was attenuated by higher childhood gender nonconformity (moderation coefficient = -0.09, 95% CI: -0.13, -0.04).

A pleiotropic variant in DNAJB4 is associated with multiple myeloma risk

Pleiotropy, which consists of a single gene or allelic variant affecting multiple unrelated traits, is common across cancers, with evidence for genome-wide significant loci shared across cancer and noncancer traits. This feature is particularly relevant in multiple myeloma (MM) because several susceptibility loci that have been identified to date are pleiotropic. Therefore, the aim of this study was to identify novel pleiotropic variants involved in MM risk using 28 684 independent single nucleotide polymorphisms (SNPs) from GWAS Catalog that reached a significant association (P < 5 × 10^-8 ) with their respective trait. The selected SNPs were analyzed in 2434 MM cases and 3446 controls from the International Lymphoma Epidemiology Consortium (InterLymph). The 10 SNPs showing the strongest associations with MM risk in InterLymph were selected for replication in an independent set of 1955 MM cases and 1549 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and 418 MM cases and 147 282 controls from the FinnGen project. The combined analysis of the three studies identified an association between DNAJB4-rs34517439-A and an increased risk of developing MM (OR = 1.22, 95%CI 1.13-1.32, P = 4.81 × 10^-7 ). rs34517439-A is associated with a modified expression of the FUBP1 gene, which encodes a multifunctional DNA and RNA-binding protein that it was observed to influence the regulation of various genes involved in cell cycle regulation, among which various oncogenes and oncosuppressors. In conclusion, with a pleiotropic scan approach we identified DNAJB4-rs34517439 as a potentially novel MM risk locus.

Phage-resistant mutations impact bacteria susceptibility to future phage infections and antibiotic response

Bacteriophage (phage) therapy in combination with antibiotic treatment serves as a potential strategy to overcome the continued rise in antibiotic resistance across bacterial pathogens. Understanding the impacts of evolutionary and ecological processes to the phage-antibiotic-resistance dynamic could advance the development of such combinatorial therapy. We tested whether the acquisition of mutations conferring phage resistance may have antagonistically pleiotropic consequences for antibiotic resistance. First, to determine the robustness of phage resistance across different phage strains, we infected resistant Escherichia coli cultures with phage that were not previously encountered. We found that phage-resistant E. coli mutants that gained resistance to a single phage strain maintain resistance to other phages with overlapping adsorption methods. Mutations underlying the phage-resistant phenotype affects lipopolysaccharide (LPS) structure and/or synthesis. Because LPS is implicated in both phage infection and antibiotic response, we then determined whether phage-resistant trade-offs exist when challenged with different classes of antibiotics. We found that only 1 out of the 4 phage-resistant E. coli mutants yielded trade-offs between phage and antibiotic resistance. Surprisingly, when challenged with novobiocin, we uncovered evidence of synergistic pleiotropy for some mutants allowing for greater antibiotic resistance, even though antibiotic resistance was never selected for. Our results highlight the importance of understanding the role of selective pressures and pleiotropic interactions in the bacterial response to phage-antibiotic combinatorial therapy.

Network Topology Can Explain Differences in Pleiotropy Between Cis- and Trans-regulatory Mutations

A mutation's degree of pleiotropy (i.e., the number of traits it alters) is predicted to impact the probability of the mutation being detrimental to fitness. For mutations that impact gene expression, mutations acting in cis have been hypothesized to generally be less pleiotropic than mutations affecting the same gene's expression in trans, suggesting that cis-regulatory mutations should be less deleterious and more likely to fix over evolutionary time. Here, we use expression and fitness data from Saccharomyces cerevisiae gene deletion strains to test these hypotheses. By treating deletion of each gene as a cis-regulatory mutation affecting its own expression and deletions of other genes affecting expression of this focal gene as trans-regulatory mutations, we find that cis-acting mutations do indeed tend to be less pleiotropic than trans-acting mutations affecting expression of the same gene. This pattern was observed for the vast majority of genes in the data set and could be explained by the topology of the regulatory network controlling gene expression. Comparing the fitness of cis- and trans-acting mutations affecting expression of the same gene also confirmed that trans-acting deletions tend to be more deleterious. These findings provide strong support for pleiotropy playing a role in the preferential fixation of cis-regulatory alleles over evolutionary time.

Novel genetic loci associated with osteoarthritis in multi-ancestry analyses in the Million Veteran Program and UK Biobank

Osteoarthritis is a common progressive joint disease. As no effective medical interventions are available, osteoarthritis often progresses to the end stage, in which only surgical options such as total joint replacement are available. A more thorough understanding of genetic influences of osteoarthritis is essential to develop targeted personalized approaches to treatment, ideally long before the end stage is reached. To date, there have been no large multiancestry genetic studies of osteoarthritis. Here, we leveraged the unique resources of 484,374 participants in the Million Veteran Program and UK Biobank to address this gap. Analyses included participants of European, African, Asian and Hispanic descent. We discovered osteoarthritis-associated genetic variation at 10 loci and replicated findings from previous osteoarthritis studies. We also present evidence that some osteoarthritis-associated regions are robust to population ancestry. Drug repurposing analyses revealed enrichment of targets of several medication classes and provide potential insight into the etiology of beneficial effects of antiepileptics on osteoarthritis pain.

Unintended Consequences of Plant Domestication

Human selection on wild populations mostly favored a common set of plant traits during domestication. This process of direct selection also altered other independent traits that were not directly perceived or desired during crop domestication and improvement. A deeper knowledge of the inadvertent and undesirable phenotypic effects and their underlying genetic causes can help design strategies to mitigate their effects and improve genetic gain in crop plants. We review different factors explaining the negative consequences of plant domestication at the phenotypic and genomic levels. We further describe the genetic causes of undesirable effects that originate from the selection of favorable alleles during plant domestication. In addition, we propose strategies that could be useful in attenuating such effects for crop improvement. With novel -omics and genome-editing tools, it is relatively approachable to understand and manipulate the genetic and biochemical mechanisms responsible for the undesirable phenotypes in domesticated plants.

Hidden paths to endless forms most wonderful: ecology latently shapes evolution of multicellular development in predatory bacteria

Ecological causes of developmental evolution, for example from predation, remain much investigated, but the potential importance of latent phenotypes in eco-evo-devo has received little attention. Using the predatory bacterium Myxococcus xanthus, which undergoes aggregative fruiting body development upon starvation, we tested whether adaptation to distinct growth environments that do not induce development latently alters developmental phenotypes under starvation conditions that do induce development. In an evolution experiment named MyxoEE-3, growing M. xanthus populations swarmed across agar surfaces while adapting to conditions varying at factors such as surface stiffness or prey identity. Such ecological variation during growth was found to greatly impact the latent evolution of development, including fruiting body morphology, the degree of morphological trait correlation, reaction norms, degrees of developmental plasticity and stochastic diversification. For example, some prey environments promoted retention of developmental proficiency whereas others led to its systematic loss. Our results have implications for understanding evolutionary interactions among predation, development and motility in myxobacterial life cycles, and, more broadly, how ecology can profoundly shape the evolution of developmental systems latently rather than by direct selection on developmental features.

Individual Genetic Heterogeneity

Genetic variation has been widely covered in literature, however, not from the perspective of an individual in any species. Here, a synthesis of genetic concepts and variations relevant for individual genetic constitution is provided. All the different levels of genetic information and variation are covered, ranging from whether an organism is unmixed or hybrid, has variations in genome, chromosomes, and more locally in DNA regions, to epigenetic variants or alterations in selfish genetic elements. Genetic constitution and heterogeneity of microbiota are highly relevant for health and wellbeing of an individual. Mutation rates vary widely for variation types, e.g., due to the sequence context. Genetic information guides numerous aspects in organisms. Types of inheritance, whether Mendelian or non-Mendelian, zygosity, sexual reproduction, and sex determination are covered. Functions of DNA and functional effects of variations are introduced, along with mechanism that reduce and modulate functional effects, including TARAR countermeasures and intraindividual genetic conflict. TARAR countermeasures for tolerance, avoidance, repair, attenuation, and resistance are essential for life, integrity of genetic information, and gene expression. The genetic composition, effects of variations, and their expression are considered also in diseases and personalized medicine. The text synthesizes knowledge and insight on individual genetic heterogeneity and organizes and systematizes the central concepts.

What snakes and caecilians have in common? Molecular interaction units and the independent origins of similar morphotypes in Tetrapoda

Developmental pathways encompass transcription factors and cis-regulatory elements that interact as transcription factor-regulatory element (TF-RE) units. Independent origins of similar phenotypes likely involve changes in different parts of these units, a hypothesis promisingly tested addressing the evolution of the rib-associated lumbar (RAL) morphotype that characterizes emblematic animals such as snakes and elephants. Previous investigation in these lineages identified a polymorphism in the Homology region 1 [H1] enhancer of the Myogenic factor-5 [Myf5], which interacts with HOX10 proteins to modulate rib development. Here we address the evolution of TF-RE units focusing on independent origins of RAL morphotypes. We compiled an extensive database for H1-Myf5 and HOX10 sequences with two goals: (i) evaluate if the enhancer polymorphism is present in amphibians exhibiting the RAL morphotype and (ii) test a hypothesis of enhanced evolutionary flexibility mediated by TF-RE units, according to which independent origins of the RAL morphotype might involve changes in either component of the interaction unit. We identified the H1-Myf5 polymorphism in lineages that diverged around 340 Ma, including Lissamphibia. Independent origins of the RAL morphotype in Tetrapoda involved sequence variation in either component of the TF-RE unit, confirming that different changes may similarly affect the phenotypic outcome of a given developmental pathway.

Genetics behind Cerebral Disease with Ocular Comorbidity: Finding Parallels between the Brain and Eye Molecular Pathology

Cerebral visual impairments (CVIs) is an umbrella term that categorizes miscellaneous visual defects with parallel genetic brain disorders. While the manifestations of CVIs are diverse and ambiguous, molecular diagnostics stand out as a powerful approach for understanding pathomechanisms in CVIs. Nevertheless, the characterization of CVI disease cohorts has been fragmented and lacks integration. By revisiting the genome-wide and phenome-wide association studies (GWAS and PheWAS), we clustered a handful of renowned CVIs into five ontology groups, namely ciliopathies (Joubert syndrome, Bardet-Biedl syndrome, Alstrom syndrome), demyelination diseases (multiple sclerosis, Alexander disease, Pelizaeus-Merzbacher disease), transcriptional deregulation diseases (Mowat-Wilson disease, Pitt-Hopkins disease, Rett syndrome, Cockayne syndrome, X-linked alpha-thalassaemia mental retardation), compromised peroxisome disorders (Zellweger spectrum disorder, Refsum disease), and channelopathies (neuromyelitis optica spectrum disorder), and reviewed several mutation hotspots currently found to be associated with the CVIs. Moreover, we discussed the common manifestations in the brain and the eye, and collated animal study findings to discuss plausible gene editing strategies for future CVI correction.