The many faces of pleiotropy

Maintenance of quantitative genetic variance in complex, multi-trait phenotypes: The contribution of rare, large effect variants in two Drosophila species

The interaction of evolutionary processes to determine quantitative genetic variation has implications for contemporary and future phenotypic evolution, as well as for our ability to detect causal genetic variants. While theoretical studies have provided robust predictions to discriminate among competing models, empirical assessment of these has been limited. In particular, theory highlights the importance of pleiotropy in resolving observations of selection and mutation, but empirical investigations have typically been limited to few traits. Here, we applied high-dimensional Bayesian Sparse Factor Genetic modeling to gene expression datasets in 2 species, Drosophila melanogaster and Drosophila serrata, to explore the distributions of genetic variance across high-dimensional phenotypic space. Surprisingly, most of the heritable trait covariation was due to few lines (genotypes) with extreme [>3 interquartile ranges (IQR) from the median] values. Intriguingly, while genotypes extreme for a multivariate factor also tended to have a higher proportion of individual traits that were extreme, we also observed genotypes that were extreme for multivariate factors but not for any individual trait. We observed other consistent differences between heritable multivariate factors with outlier lines vs those factors without extreme values, including differences in gene functions. We use these observations to identify further data required to advance our understanding of the evolutionary dynamics and nature of standing genetic variation for quantitative traits.

Phenotypic screening of low molecular weight compounds is rich ground for repurposed, on-target drugs

A target-based drug discovery strategy has led to a bias away from low molecular weight (MWT) drug discovery. Analysis of the ACS chemistry registration system shows that most low MWT drugs were first made in the time era before target-based drug discovery. Therapeutic activity among most low MWT drugs was identified in the era of phenotypic drug discovery when drugs were selected based on their phenotypic effects and before in vitro screening, mechanism of action considerations and experiences with fragment screening became known. The common perception that drugs cannot be found among low MWT compounds is incorrect based on both drug discovery history and our own experience with MLR-1023. The greater proportion of low MWT compounds that are commercially available compared to higher MWT compounds is a factor that should facilitate biology study. We posit that low MWT compounds are more suited to identification of new therapeutic activity using phenotypic screens provided that the phenotypic screening method has enough screening capacity. On-target and off-target therapeutic activities are discussed from both a chemistry and biology perspective because of a concern that either phenotypic or low MWT drug discovery might bias towards promiscuous compounds that combine on-target and off-target effects. Among ideal drug repositioning candidates (late-stage pre-clinical or clinically-experience compounds), pleiotropic activity (multiple therapeutic actions) is far more likely due to on-target effects arising where a single target mediates multiple therapeutic benefits, a desirable outcome for drug development purposes compared to the off-target alternative. Our exemplar of a low MWT compound, MLR-1023, discovered by phenotypic screening and subsequently found to have a single mechanism of action would have been overlooked based on current era medicinal chemistry precedent. The diverse therapeutic activities described for this compound by us, and others arise from the same pleiotropic lyn kinase activation molecular target. MLR-1023 serves as a proof-of-principle that potent, on target, low MWT drugs can be discovered by phenotypic screening.

Estimating the causal effect of liability to disease on healthcare costs using Mendelian Randomization

Accurate measurement of the effects of disease status on healthcare costs is important in the pragmatic evaluation of interventions but is complicated by endogeneity bias. Mendelian Randomization, the use of random perturbations in germline genetic variation as instrumental variables, can avoid these limitations. We used a novel Mendelian Randomization analysis to model the causal impact on inpatient hospital costs of liability to six prevalent diseases and health conditions: asthma, eczema, migraine, coronary heart disease, Type 2 diabetes, and depression. We identified genetic variants from replicated genome-wide associations studies and estimated their association with inpatient hospital costs on over 300,000 individuals. There was concordance of findings across varieties of sensitivity analyses, including stratification by sex and methods robust to violations of the exclusion restriction. Results overall were imprecise and we could not rule out large effects of liability to disease on healthcare costs. In particular, genetic liability to coronary heart disease had substantial impacts on costs.

Solving the grand challenge of phenotypic integration: allometry across scales

Phenotypic integration is a concept related to the cascade of trait relationships from the lowest organizational levels, i.e. genes, to the highest, i.e. whole-organism traits. However, the cause-and-effect linkages between traits are notoriously difficult to determine. In particular, we still lack a mathematical framework to model the relationships involved in the integration of phenotypic traits. Here, we argue that allometric models developed in ecology offer testable mathematical equations of trait relationships across scales. We first show that allometric relationships are pervasive in biology at different organizational scales and in different taxa. We then present mechanistic models that explain the origin of allometric relationships. In addition, we emphasized that recent studies showed that natural variation does exist for allometric parameters, suggesting a role for genetic variability, selection and evolution. Consequently, we advocate that it is time to examine the genetic determinism of allometries, as well as to question in more detail the role of genome size in subsequent scaling relationships. More broadly, a possible-but so far neglected-solution to understand phenotypic integration is to examine allometric relationships at different organizational levels (cell, tissue, organ, organism) and in contrasted species.

Comparative chemical genomic profiling across plant-based hydrolysate toxins reveals widespread antagonism in fitness contributions

The budding yeast Saccharomyces cerevisiae has been used extensively in fermentative industrial processes, including biofuel production from sustainable plant-based hydrolysates. Myriad toxins and stressors found in hydrolysates inhibit microbial metabolism and product formation. Overcoming these stresses requires mitigation strategies that include strain engineering. To identify shared and divergent mechanisms of toxicity and to implicate gene targets for genetic engineering, we used a chemical genomic approach to study fitness effects across a library of S. cerevisiae deletion mutants cultured anaerobically in dozens of individual compounds found in different types of hydrolysates. Relationships in chemical genomic profiles identified classes of toxins that provoked similar cellular responses, spanning inhibitor relationships that were not expected from chemical classification. Our results also revealed widespread antagonistic effects across inhibitors, such that the same gene deletions were beneficial for surviving some toxins but detrimental for others. This work presents a rich dataset relating gene function to chemical compounds, which both expands our understanding of plant-based hydrolysates and provides a useful resource to identify engineering targets.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

Pleiotropic effects of trans-regulatory mutations on fitness and gene expression

Variation in gene expression arises from cis- and trans-regulatory mutations, which contribute differentially to expression divergence. We compare the impacts on gene expression and fitness resulting from cis- and trans-regulatory mutations in Saccharomyces cerevisiae, with a focus on the TDH3 gene. We use the effects of cis-regulatory mutations to infer effects of trans-regulatory mutations attributable to impacts beyond the focal gene, revealing a distribution of pleiotropic effects. Cis- and trans-regulatory mutations had different effects on gene expression with pleiotropic effects of trans-regulatory mutants affecting expression of genes both in parallel to and downstream of the focal gene. The more widespread and deleterious effects of trans-regulatory mutations we observed are consistent with their decreasing relative contribution to expression differences over evolutionary time.

Three laws of teleonometrics

We define teleonometrics as the theoretical and empirical study of teleonomy. We propose three laws for teleonometrics. The first law describes the hierarchical organization of teleonomic functions across biological levels from genes to individuals. According to this law, the number of goal-directed functions increases from individuals (one goal, maximizing inclusive fitness) to intermediate levels and to genes and alleles (myriad time-, space- and context-dependent goals, depending upon degrees and patterns of pleiotropy). The second law describes the operation of teleonomic functions under trade-offs, coadaptations and negative and positive pleiotropies, which are universal in biological systems. According to this law, the functions of an allele, gene or trait are described and defined by patterns of antagonistic (trading off) and compatible (coadapted) functions. The third law of teleonometrics is that the major transitions in evolution are driven by the origins of novel, emergent goals associated with functional changes and by the breaking and reshaping of trade-offs, especially by mechanisms involving increases in resources or time, and new divisions of labour or function. We illustrate the application of these laws using data from three empirical vignettes, which help to show the usefulness of teleonometric viewpoints for understanding the interfaces between function, trade-offs and dysfunctions manifest as disease.

PolarMorphism enables discovery of shared genetic variants across multiple traits from GWAS summary statistics

MOTIVATION

Pleiotropic SNPs are associated with multiple traits. Such SNPs can help pinpoint biological processes with an effect on multiple traits or point to a shared etiology between traits. We present PolarMorphism, a new method for the identification of pleiotropic SNPs from genome-wide association studies (GWAS) summary statistics. PolarMorphism can be readily applied to more than two traits or whole trait domains. PolarMorphism makes use of the fact that trait-specific SNP effect sizes can be seen as Cartesian coordinates and can thus be converted to polar coordinates r (distance from the origin) and theta (angle with the Cartesian x-axis, in the case of two traits). r describes the overall effect of a SNP, while theta describes the extent to which a SNP is shared. r and theta are used to determine the significance of SNP sharedness, resulting in a P-value per SNP that can be used for further analysis.

RESULTS

We apply PolarMorphism to a large collection of publicly available GWAS summary statistics enabling the construction of a pleiotropy network that shows the extent to which traits share SNPs. We show how PolarMorphism can be used to gain insight into relationships between traits and trait domains and contrast it with genetic correlation. Furthermore, pathway analysis of the newly discovered pleiotropic SNPs demonstrates that analysis of more than two traits simultaneously yields more biologically relevant results than the combined results of pairwise analysis of the same traits. Finally, we show that PolarMorphism is more efficient and more powerful than previously published methods.

AVAILABILITY AND IMPLEMENTATION

code: https://github.com/UMCUGenetics/PolarMorphism, results: 10.5281/zenodo.5844193.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Open problems in human trait genetics

Genetic studies of human traits have revolutionized our understanding of the variation between individuals, and yet, the genetics of most traits is still poorly understood. In this review, we highlight the major open problems that need to be solved, and by discussing these challenges provide a primer to the field. We cover general issues such as population structure, epistasis and gene-environment interactions, data-related issues such as ancestry diversity and rare genetic variants, and specific challenges related to heritability estimates, genetic association studies, and polygenic risk scores. We emphasize the interconnectedness of these problems and suggest promising avenues to address them.

An Analysis of Different Distance-Linkage Methods for Clustering Gene Expression Data and Observing Pleiotropy: Empirical Study

BACKGROUND

Large amounts of biological data have been generated over the last few decades, encouraging scientists to look for connections between genes that cause various diseases. Clustering illustrates such a relationship between numerous species and genes. Finding an appropriate distance-linkage metric to construct clusters from diverse biological data sets has thus become critical. Pleiotropy is also important for a gene's expression to vary and create varied consequences in living things. Finding the pleiotropy of genes responsible for various diseases has become a major research challenge.

OBJECTIVE

Our goal was to establish the optimal distance-linkage strategy for creating reliable clusters from diverse data sets and identifying the common genes that cause various tumors to observe genes with pleiotropic effect.

METHODS

We considered 4 linking methods-single, complete, average, and ward-and 3 distance metrics-Euclidean, maximum, and Manhattan distance. For assessing the quality of different sets of clusters, we used a fitness function that combines silhouette width and within-cluster distance.

RESULTS

According to our findings, the maximum distance measure produces the highest-quality clusters. Moreover, for medium data set, the average linkage method, and for large data set, the ward linkage method works best. The outcome is not improved by using ensemble clustering. We also discovered genes that cause 3 different cancers and used gene enrichment to confirm our findings.

CONCLUSIONS

Accuracy is crucial in clustering, and we investigated the accuracy of numerous clustering techniques in our research. Other studies may aid related works if the data set is similar to ours.

Pleiotropic constraints promote the evolution of cooperation in cellular groups

The evolution of cooperation in cellular groups is threatened by lineages of cheaters that proliferate at the expense of the group. These cell lineages occur within microbial communities, and multicellular organisms in the form of tumours and cancer. In contrast to an earlier study, here we show how the evolution of pleiotropic genetic architectures-which link the expression of cooperative and private traits-can protect against cheater lineages and allow cooperation to evolve. We develop an age-structured model of cellular groups and show that cooperation breaks down more slowly within groups that tie expression to a private trait than in groups that do not. We then show that this results in group selection for pleiotropy, which strongly promotes cooperation by limiting the emergence of cheater lineages. These results predict that pleiotropy will rapidly evolve, so long as groups persist long enough for cheater lineages to threaten cooperation. Our results hold when pleiotropic links can be undermined by mutations, when pleiotropy is itself costly, and in mixed-genotype groups such as those that occur in microbes. Finally, we consider features of multicellular organisms-a germ line and delayed reproductive maturity-and show that pleiotropy is again predicted to be important for maintaining cooperation. The study of cancer in multicellular organisms provides the best evidence for pleiotropic constraints, where abberant cell proliferation is linked to apoptosis, senescence, and terminal differentiation. Alongside development from a single cell, we propose that the evolution of pleiotropic constraints has been critical for cooperation in many cellular groups.

A synthetic synthesis to explore animal evolution and development

Identifying the general principles by which genotypes are converted into phenotypes remains a challenge in the post-genomic era. We still lack a predictive understanding of how genes shape interactions among cells and tissues in response to signalling and environmental cues, and hence how regulatory networks generate the phenotypic variation required for adaptive evolution. Here, we discuss how techniques borrowed from synthetic biology may facilitate a systematic exploration of evolvability across biological scales. Synthetic approaches permit controlled manipulation of both endogenous and fully engineered systems, providing a flexible platform for investigating causal mechanisms in vivo. Combining synthetic approaches with multi-level phenotyping (phenomics) will supply a detailed, quantitative characterization of how internal and external stimuli shape the morphology and behaviour of living organisms. We advocate integrating high-throughput experimental data with mathematical and computational techniques from a variety of disciplines in order to pursue a comprehensive theory of evolution.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Genetic basis of speciation and adaptation: from loci to causative mutations

Does evolution proceed in small steps or large leaps? How repeatable is evolution? How constrained is the evolutionary process? Answering these long-standing questions in evolutionary biology is indispensable for both understanding how extant biodiversity has evolved and predicting how organisms and ecosystems will respond to changing environments in the future. Understanding the genetic basis of phenotypic diversification and speciation in natural populations is key to properly answering these questions. The leap forward in genome sequencing technologies has made it increasingly easier to not only investigate the genetic architecture but also identify the variant sites underlying adaptation and speciation in natural populations. Furthermore, recent advances in genome editing technologies are making it possible to investigate the functions of each candidate gene in organisms from natural populations. In this article, we discuss how these recent technological advances enable the analysis of causative genes and mutations and how such analysis can help answer long-standing evolutionary biology questions.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Natural variation in Drosophila shows weak pleiotropic effects

Background

Pleiotropy describes the phenomenon in which a gene affects multiple phenotypes. The extent of pleiotropy is still disputed, mainly because of issues of inadequate power of analyses. A further challenge is that empirical tests of pleiotropy are restricted to a small subset of all possible phenotypes. To overcome these limitations, we propose a new measurement of pleiotropy that integrates across many phenotypes and multiple generations to improve power.

Results

We infer pleiotropy from the fitness cost imposed by frequency changes of pleiotropic loci. Mixing Drosophila simulans populations, which adapted independently to the same new environment using different sets of genes, we show that the adaptive frequency changes have been accompanied by measurable fitness costs.

Conclusions

Unlike previous studies characterizing the molecular basis of pleiotropy, we show that many loci, each of weak effect, contribute to genome-wide pleiotropy. We propose that the costs of pleiotropy are reduced by the modular architecture of gene expression, which facilitates adaptive gene expression changes with low impact on other functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02680-4.

The acquisitive-conservative axis of leaf trait variation emerges even in homogeneous environments

BACKGROUND AND AIMS

The acquisitive-conservative axis of plant ecological strategies results in a pattern of leaf trait covariation that captures the balance between leaf construction costs and plant growth potential. Studies evaluating trait covariation within species are scarcer, and have mostly dealt with variation in response to environmental gradients. Little work has been published on intraspecific patterns of leaf trait covariation in the absence of strong environmental variation.

METHODS

We analysed covariation of four leaf functional traits [specific leaf area (SLA) leaf dry matter content (LDMC), force to tear (Ft) and leaf nitrogen content (Nm)] in six Poaceae and four Fabaceae species common in the dry Chaco forest of Central Argentina, growing in the field and in a common garden. We compared intraspecific covariation patterns (slopes, correlation and effect size) of leaf functional traits with global interspecific covariation patterns. Additionally, we checked for possible climatic and edaphic factors that could affect the intraspecific covariation pattern.

KEY RESULTS

We found negative correlations for the LDMC-SLA, Ft-SLA, LDMC-Nm and Ft-Nm trait pairs. This intraspecific covariation pattern found both in the field and in the common garden and not explained by climatic or edaphic variation in the field follows the expected acquisitive-conservative axis. At the same time, we found quantitative differences in slopes among different species, and between these intraspecific patterns and the interspecific ones. Many of these differences seem to be idiosyncratic, but some appear consistent among species (e.g. all the intraspecific LDMC-SLA and LDMC-Nm slopes tend to be shallower than the global pattern).

CONCLUSIONS

Our study indicates that the acquisitive-conservative leaf functional trait covariation pattern occurs at the intraspecific level even in the absence of relevant environmental variation in the field. This suggests a high degree of variation-covariation in leaf functional traits not driven by environmental variables.

riceExplorer: Uncovering the Hidden Potential of a National Genomic Resource Against a Global Database

Agricultural crop breeding programs, particularly at the national level, typically consist of a core panel of elite breeding cultivars alongside a number of local landrace varieties (or other endemic cultivars) that provide additional sources of phenotypic and genomic variation or contribute as experimental materials (e.g., in GWAS studies). Three issues commonly arise. First, focusing primarily on core development accessions may mean that the potential contributions of landraces or other secondary accessions may be overlooked. Second, elite cultivars may accumulate deleterious alleles away from nontarget loci due to the strong effects of artificial selection. Finally, a tendency to focus solely on SNP-based methods may cause incomplete or erroneous identification of functional variants. In practice, integration of local breeding programs with findings from global database projects may be challenging. First, local GWAS experiments may only indicate useful functional variants according to the diversity of the experimental panel, while other potentially useful loci-identifiable at a global level-may remain undiscovered. Second, large-scale experiments such as GWAS may prove prohibitively costly or logistically challenging for some agencies. Here, we present a fully automated bioinformatics pipeline (riceExplorer) that can easily integrate local breeding program sequence data with international database resources, without relying on any phenotypic experimental procedure. It identifies associated functional haplotypes that may prove more robust in determining the genotypic determinants of desirable crop phenotypes. In brief, riceExplorer evaluates a global crop database (IRRI 3000 Rice Genomes) to identify haplotypes that are associated with extreme phenotypic variation at the global level and recorded in the database. It then examines which potentially useful variants are present in the local crop panel, before distinguishing between those that are already incorporated into the elite breeding accessions and those only found among secondary varieties (e.g., landraces). Results highlight the effectiveness of our pipeline, identifying potentially useful functional haplotypes across the genome that are absent from elite cultivars and found among landraces and other secondary varieties in our breeding program. riceExplorer can automatically conduct a full genome analysis and produces annotated graphical output of chromosomal maps, potential global diversity sources, and summary tables.

Genetic correlation, pleiotropy, and causal associations between substance use and psychiatric disorder

BACKGROUND

Substance use occurs at a high rate in persons with a psychiatric disorder. Genetically informative studies have the potential to elucidate the etiology of these phenomena. Recent developments in genome-wide association studies (GWAS) allow new avenues of investigation.

METHOD

Using results of GWAS meta-analyses, we performed a factor analysis of the genetic correlation structure, a genome-wide search of shared loci, and causally informative tests for six substance use phenotypes (four smoking, one alcohol, and one cannabis use) and five psychiatric disorders (ADHD, anorexia, depression, bipolar disorder, and schizophrenia).

RESULTS

Two correlated externalizing and internalizing/psychosis factor were found, although model fit was beneath conventional standards. Of 458 loci reported in previous univariate GWAS of substance use and psychiatric disorders, about 50% (230 loci) were pleiotropic with additional 111 pleiotropic loci not reported from past GWAS. Of the 341 pleiotropic loci, 152 were associated with both substance use and psychiatric disorders, implicating neurodevelopment, cell morphogenesis, biological adhesion pathways, and enrichment in 13 different brain tissues. Seventy-five and 114 pleiotropic loci were specific to either psychiatric disorders or substance use phenotypes, implicating neuronal signaling pathway and clathrin-binding functions/structures, respectively. No consistent evidence for phenotypic causation was found across different Mendelian randomization methods.

CONCLUSIONS

Genetic etiology of substance use and psychiatric disorders is highly pleiotropic and involves shared neurodevelopmental path, neurotransmission, and intracellular trafficking. In aggregate, the patterns are not consistent with vertical pleiotropy, more likely reflecting horizontal pleiotropy or more complex forms of phenotypic causation.

Genetics in the ADHD Clinic: How Can Genetic Testing Support the Current Clinical Practice?

Attention-deficit/hyperactivity disorder (ADHD) is a neurodevelopmental disorder with a childhood prevalence of 5%. In about two-thirds of the cases, ADHD symptoms persist into adulthood and often cause significant functional impairment. Based on the results of family and twin studies, the estimated heritability of ADHD approximates 80%, suggests a significant genetic component in the etiological background of the disorder; however, the potential genetic effects on disease risk, symptom severity, and persistence are unclear. This article provides a brief review of the genome-wide and candidate gene association studies with a focus on the clinical aspects, summarizing findings of ADHD disease risk, ADHD core symptoms as dimensional traits, and other traits frequently associated with ADHD, which may contribute to the susceptibility to other comorbid psychiatric disorders. Furthermore, neuropsychological impairment and measures from neuroimaging and electrophysiological paradigms, emerging as potential biomarkers, also provide a prominent target for molecular genetic studies, since they lie in the pathway from genes to behavior; therefore, they can contribute to the understanding of the underlying neurobiological mechanisms and the interindividual heterogeneity of clinical symptoms. Beyond the aforementioned aspects, throughout the review, we also give a brief summary of the genetic results, including polygenic risk scores that can potentially predict individual response to different treatment options and may offer a possibility for personalized treatment for the therapy of ADHD in the future.

Genetic Risk of Cardiovascular Disease Is Associated with Macular Ganglion Cell-Inner Plexiform Layer Thinning in an Early Glaucoma Cohort

PURPOSE

To evaluate the association between genetic risk for cardiovascular disease and retinal thinning in early glaucoma.

DESIGN

Prospective, observational genetic association study.

PARTICIPANTS

Multicohort study combining a cohort of patients with suspect and early manifest primary open-angle glaucoma (POAG), a cohort of patients with perimetric POAG, and an external normative control cohort.

METHODS

A cardiovascular disease genetic risk score was calculated for 828 participants from the Progression Risk of Glaucoma: Relevant SNPs [single nucleotide polymorphisms] with Significant Association (PROGRESSA) study. Participants were characterized as showing either predominantly macular ganglion cell-inner plexiform layer (GCIPL), predominantly peripapillary retinal nerve fiber layer (pRNFL) or equivalent macular GCIPL and pRNFL spectral-domain OCT thinning. The cardiovascular disease genetic risk scores for these groups were compared to an internal reference group of stable suspected glaucoma and of an external normative population. Replication was undertaken by comparing the phenotypes of participants from the Australia New Zealand Registry of Advanced Glaucoma (ANZRAG) with the normative control group.

MAIN OUTCOME MEASURES

Spectral-domain OCT and Humphrey Visual Field (HVF) change.

RESULTS

After accounting for age, sex, and intraocular pressure (IOP), participants with predominantly macular GCIPL thinning showed a higher cardiovascular disease genetic risk score than reference participants (odds ratio [OR], 1.76/standard deviation [SD]; 95% confidence interval [CI], 1.18-2.62; P = 0.005) and than normative participants (OR, 1.32/SD; 95% CI, 1.12-1.54; P = 0.002). This finding was replicated by comparing ANZRAG participants with predominantly macular GCIPL change with the normative population (OR, 1.39/SD; 95% CI, 1.05-1.83; P = 0.022). Review of HVF data identified that participants with paracentral visual field defects also demonstrated a higher cardiovascular disease genetic risk score than reference participants (OR, 1.85/SD; 95% CI, 1.16-2.97; P = 0.010). Participants with predominantly macular GCIPL thinning exhibited a higher vertical cup-to-disc ratio genetic risk score (OR, 1.48/SD; 95% CI, 1.24-1.76; P < 0.001), but an IOP genetic risk score (OR, 1.12/SD; 95% CI, 0.95-1.33; P = 0.179) comparable with that of the normative population.

CONCLUSIONS

This study highlighted the relationship between cardiovascular disease and retinal thinning in suspect and manifest glaucoma cases.

High-value pleiotropic genes for developing multiple stress-tolerant biofortified crops for 21st-century challenges

The agriculture-based livelihood systems that are already vulnerable due to multiple challenges face immediate risk of increased crop failures due to weather vagaries. As breeders and biotechnologists, our strategy is to advance and innovate breeding for weather-proofing crops. Plant stress tolerance is a genetically complex trait. Additionally, crops rarely face a single type of stress in isolation, and it is difficult for plants to deal with multiple stresses simultaneously. One of the most helpful approaches to creating stress-resilient crops is genome editing and trans- or cis-genesis. Out of hundreds of stress-responsive genes, many have been used to impart tolerance against a particular stress factor, while a few used in combination for gene pyramiding against multiple stresses. However, a better approach would be to use multi-role pleiotropic genes that enable plants to adapt to numerous environmental stresses simultaneously. Herein we attempt to integrate and present the scattered information published in the past three decades about these pleiotropic genes for crop improvement and remodeling future cropping systems. Research articles validating functional roles of genes in transgenic plants were used to create groups of multi-role pleiotropic genes that could be candidate genes for developing weather-proof crop varieties. These biotech crop varieties will help create 'high-value farms' to meet the goal of a sustainable increase in global food productivity and stabilize food prices by ensuring a fluctuation-free assured food supply. It could also help create a gene repository through artificial gene synthesis for 'resilient high-value food production' for the 21st century.

Plasticity, pleiotropy and fitness trade-offs in Arabidopsis genotypes with different telomere lengths

Telomere length has been implicated in the organismal response to stress, but the underlying mechanisms are unknown. Here we examine the impact of telomere length changes on the responses to three contrasting abiotic environments in Arabidopsis, and measure 32 fitness, developmental, physiological and leaf-level anatomical traits. We report that telomere length in wild-type and short-telomere mutants is resistant to abiotic stress, while the elongated telomeres in ku70 mutants are more plastic. We detected significant pleiotropic effects of telomere length on flowering time and key leaf physiological and anatomical traits. Furthermore, our data reveal a significant genotype by environment (G × E) interaction for reproductive fitness, with the benefits and costs to performance depending on the growth conditions. These results imply that life-history trade-offs between flowering time and reproductive fitness are impacted by telomere length variation. We postulate that telomere length in plants is subject to natural selection imposed by different environments.

Regenerative strategies for the consequences of myocardial infarction: Chronological indication and upcoming visions

Heart muscle injury and an elevated troponin level signify myocardial infarction (MI), which may result in defective and uncoordinated segments, reduced cardiac output, and ultimately, death. Physicians apply thrombolytic therapy, coronary artery bypass graft (CABG) surgery, or percutaneous coronary intervention (PCI) to recanalize and restore blood flow to the coronary arteries, albeit they were not convincingly able to solve the heart problems. Thus, researchers aim to introduce novel substitutional therapies for regenerating and functionalizing damaged cardiac tissue based on engineering concepts. Cell-based engineering approaches, utilizing biomaterials, gene, drug, growth factor delivery systems, and tissue engineering are the most leading studies in the field of heart regeneration. Also, understanding the primary cause of MI and thus selecting the most efficient treatment method can be enhanced by preparing microdevices so-called heart-on-a-chip. In this regard, microfluidic approaches can be used as diagnostic platforms or drug screening in cardiac disease treatment. Additionally, bioprinting technique with whole organ 3D printing of human heart with major vessels, cardiomyocytes and endothelial cells can be an ideal goal for cardiac tissue engineering and remarkable achievement in near future. Consequently, this review discusses the different aspects, advancements, and challenges of the mentioned methods with presenting the advantages and disadvantages, chronological indications, and application prospects of various novel therapeutic approaches.

Leveling with Tinbergen: Four levels simplified to causes and consequences

In 1963, Niko Tinbergen published his foundational manuscript identifying the four questions we ask in animal behavior-how does the behavior emerge across the lifespan (development); how does it work (mechanism); how and why did it evolve (evolution); and why is it adaptive (function). Tinbergen clarified that these 'levels of analysis' are complementary, not competing, thereby avoiding many fruitless scientific debates. However, the relationships among the four levels was never established. Here, we propose 'leveling' Tinbergen's questions to a single temporal timescale divided into causes (encompassing mechanism, development, and evolution) and consequences (encompassing function). Scientific advances now seamlessly link evolution, development, and mechanism into a continuum of 'causes'. The cause-consequence distinction separates the processes that precede (and lead to) a behavior, from the processes that come after (and result from) a behavior. Even for past behaviors, the functional outcomes are (historical) consequences of the causes that preceded them.

The population genetics of collateral resistance and sensitivity

Resistance mutations against one drug can elicit collateral sensitivity against other drugs. Multi-drug treatments exploiting such trade-offs can help slow down the evolution of resistance. However, if mutations with diverse collateral effects are available, a treated population may evolve either collateral sensitivity or collateral resistance. How to design treatments robust to such uncertainty is unclear. We show that many resistance mutations in Escherichia coli against various antibiotics indeed have diverse collateral effects. We propose to characterize such diversity with a joint distribution of fitness effects (JDFE) and develop a theory for describing and predicting collateral evolution based on simple statistics of the JDFE. We show how to robustly rank drug pairs to minimize the risk of collateral resistance and how to estimate JDFEs. In addition to practical applications, these results have implications for our understanding of evolution in variable environments.

Drugs known as antibiotics are the main treatment for most serious infections caused by bacteria. However, many bacteria are acquiring genetic mutations that make them resistant to the effects of one or more types of antibiotics, making them harder to eliminate.

One way to tackle drug-resistant bacteria is to develop new types of antibiotics; however, in recent years, the rate at which new antibiotics have become available has been dwindling. Using two or more existing drugs, one after another, can also be an effective way to eliminate resistant bacteria. The success of any such ‘multi-drug’ treatment lies in being able to predict whether mutations that make the bacteria resistant to one drug simultaneously make it sensitive to another, a phenomenon known as collateral sensitivity.

Different resistance mutations may have different collateral effects: some may increase the bacteria’s sensitivity to the second drug, while others might make the bacteria more resistant. However, it is currently unclear how to design robust multi-drug treatments that take this diversity of collateral effects into account.

Here, Ardell and Kryazhimskiy used a concept called JDFE (short for the joint distribution of fitness effects) to describe the diversity of collateral effects in a population of bacteria exposed to a single drug. This information was then used to mathematically model how collateral effects evolved in the population over time.

Ardell and Kryazhimskiy showed that this approach can predict how likely a population is to become collaterally sensitive or collaterally resistant to a second antibiotic. Drug pairs can then be ranked according to the risk of collateral resistance emerging, so long as information on the variety of resistance mutations available to the bacteria are included in the model.

Each year, more than 700,000 people die from infections caused by bacteria that are resistant to one or more antibiotics. The findings of Ardell and Kryazhimskiy may eventually help clinicians design multi-drug treatments that effectively eliminate bacterial infections and help to prevent more bacteria from evolving resistance to antibiotics. However, to achieve this goal, more research is needed to fully understand the range collateral effects caused by resistance mutations.

Does a complex life cycle affect adaptation to environmental change? Genome-informed insights for characterizing selection across complex life cycle

Complex life cycles, in which discrete life stages of the same organism differ in form or function and often occupy different ecological niches, are common in nature. Because stages share the same genome, selective effects on one stage may have cascading consequences through the entire life cycle. Theoretical and empirical studies have not yet generated clear predictions about how life cycle complexity will influence patterns of adaptation in response to rapidly changing environments or tested theoretical predictions for fitness trade-offs (or lack thereof) across life stages. We discuss complex life cycle evolution and outline three hypotheses—ontogenetic decoupling, antagonistic ontogenetic pleiotropy and synergistic ontogenetic pleiotropy—for how selection may operate on organisms with complex life cycles. We suggest a within-generation experimental design that promises significant insight into composite selection across life cycle stages. As part of this design, we conducted simulations to determine the power needed to detect selection across a life cycle using a population genetic framework. This analysis demonstrated that recently published studies reporting within-generation selection were underpowered to detect small allele frequency changes (approx. 0.1). The power analysis indicates challenging but attainable sampling requirements for many systems, though plants and marine invertebrates with high fecundity are excellent systems for exploring how organisms with complex life cycles may adapt to climate change.

Two-Sample Multivariable Mendelian Randomization Analysis Using R

Mendelian randomization is a framework that uses measured variation in genes for assessing and estimating the causal effect of an exposure on an outcome. Multivariable Mendelian randomization is an extension that can assess the causal effect of multiple exposures on an outcome, and can be advantageous when considering a set (>1) of potentially correlated candidate risk factors in evaluating the causal effect of each on a health outcome, accounting for measured pleiotropy. This can be seen, for example, in determining the causal effects of lipids and cholesterol on type 2 diabetes risk, where the correlated risk factors share genetic predictors. Similar to univariate Mendelian randomization, multivariable Mendelian randomization can be conducted using two-sample summary-level data where the gene-exposure and gene-outcome associations are derived from separate samples from the same underlying population. Here, we present a protocol for conducting a two-sample multivariable Mendelian randomization study using the 'MVMR' package in R and summary-level genetic data. We also provide a protocol for searching and obtaining instruments using available data sources in the 'MRInstruments' R package. Finally, we provide general guidelines and discuss the utility of performing a multivariable Mendelian randomization analysis for simultaneously assessing causality of multiple exposures. © 2021 Wiley Periodicals LLC. Basic Protocol: Performing a two-sample multivariable Mendelian randomization analysis using the 'MVMR' package in R and summarized genetic data Support Protocol 1: Installing the 'MVMR' R package Support Protocol 2: Obtaining instruments from the 'MRInstruments' R package.

SNP- and Haplotype-Based GWAS of Flowering-Related Traits in Brassica napus

Traits related to flowering time are the most promising agronomic traits that directly impact the seed yield and oil quality of rapeseed (Brassica napus L.). Developing early flowering and maturity rapeseed varieties is an important breeding objective in B. napus. Many studies have reported on days to flowering, but few have reported on budding, bolting, and the interval between bolting and DTF. Therefore, elucidating the genetic architecture of QTLs and genes regulating flowering time, we presented an integrated investigation on SNP and haplotype-based genome-wide association study of 373 diverse B. napus germplasm, which were genotyped by the 60K SNP array and were phenotyped in the four environments. The results showed that a total of 15 and 37 QTLs were detected from SNP and haplotype-based GWAS, respectively. Among them, seven QTL clusters were identified by haplotype-based GWAS. Moreover, three and eight environmentally stable QTLs were detected by SNP-GWAS and haplotype-based GWAS, respectively. By integrating the above two approaches and by co-localizing the four traits, ten (10) genomic regions were under selection on chromosomes A03, A07, A08, A10, C06, C07, and C08. Interestingly, the genomic regions FT.A07.1, FT.A08, FT.C06, and FT.C07 were identified as novel. In these ten regions, a total of 197 genes controlling FT were detected, of which 14 highly expressed DEGs were orthologous to 13 Arabidopsis thaliana genes after integration with transcriptome results. In a nutshell, the above results uncovered the genetic architecture of important agronomic traits related to flowering time and provided a basis for multiple molecular marker-trait associations in B. napus.

Discordant associations of educational attainment with ASD and ADHD implicate a polygenic form of pleiotropy

Autism Spectrum Disorder (ASD) and Attention-Deficit/Hyperactivity Disorder (ADHD) are complex co-occurring neurodevelopmental conditions. Their genetic architectures reveal striking similarities but also differences, including strong, discordant polygenic associations with educational attainment (EA). To study genetic mechanisms that present as ASD-related positive and ADHD-related negative genetic correlations with EA, we carry out multivariable regression analyses using genome-wide summary statistics (N = 10,610-766,345). Our results show that EA-related genetic variation is shared across ASD and ADHD architectures, involving identical marker alleles. However, the polygenic association profile with EA, across shared marker alleles, is discordant for ASD versus ADHD risk, indicating independent effects. At the single-variant level, our results suggest either biological pleiotropy or co-localisation of different risk variants, implicating MIR19A/19B microRNA mechanisms. At the polygenic level, they point to a polygenic form of pleiotropy that contributes to the detectable genome-wide correlation between ASD and ADHD and is consistent with effect cancellation across EA-related regions.

Genetic Support of A Causal Relationship Between Iron Status and Type 2 Diabetes: A Mendelian Randomization Study

CONTEXT

Iron overload is a known risk factor for type 2 diabetes (T2D); however, iron overload and iron deficiency have both been associated with metabolic disorders in observational studies.

OBJECTIVE

Using mendelian randomization (MR), we assessed how genetically predicted systemic iron status affected T2D risk.

METHODS

A 2-sample MR analysis was used to obtain a causal estimate. We selected genetic variants strongly associated (P < 5 × 10-8) with 4 biomarkers of systemic iron status from a study involving 48 972 individuals performed by the Genetics of Iron Status consortium and applied these biomarkers to the T2D case-control study (74 124 cases and 824 006 controls) performed by the Diabetes Genetics Replication and Meta-analysis consortium. The simple median, weighted median, MR-Egger, MR analysis using mixture-model, weighted allele scores, and MR based on a Bayesian model averaging approaches were used for the sensitivity analysis.

RESULTS

Genetically instrumented serum iron (odds ratio [OR]: 1.07; 95% CI, 1.02-1.12), ferritin (OR: 1.19; 95% CI, 1.08-1.32), and transferrin saturation (OR: 1.06; 95% CI, 1.02-1.09) were positively associated with T2D. In contrast, genetically instrumented transferrin, a marker of reduced iron status, was inversely associated with T2D (OR: 0.91; 95% CI, 0.87-0.96).

CONCLUSION

Genetic evidence supports a causal link between increased systemic iron status and increased T2D risk. Further studies involving various ethnic backgrounds based on individual-level data and studies regarding the underlying mechanism are warranted for reducing the risk of T2D.

White Matter Integrity and Nicotine Dependence: Evaluating Vertical and Horizontal Pleiotropy

Tobacco smoking is an addictive behavior that supports nicotine dependence and is an independent risk factor for cancer and other illnesses. Its neurogenetic mechanisms are not fully understood but may act through alterations in the cerebral white matter (WM). We hypothesized that the vertical pleiotropic pathways, where genetic variants influence a trait that in turn influences another trait, link genetic factors, integrity of cerebral WM, and nicotine addiction. We tested this hypothesis using individual genetic factors, WM integrity measured by fractional anisotropy (FA), and nicotine dependence-related smoking phenotypes, including smoking status (SS) and cigarettes per day (CPDs), in a large epidemiological sample collected by the UK Biobank. We performed a genome-wide association study (GWAS) to identify previously reported loci associated with smoking behavior. Smoking was found to be associated with reduced WM integrity in multiple brain regions. We then evaluated two competing vertical pathways: Genes → WM integrity → Smoking versus Genes → Smoking → WM integrity and a horizontal pleiotropy pathway where genetic factors independently affect both smoking and WM integrity. The causal pathway analysis identified 272 pleiotropic single-nucleotide polymorphisms (SNPs) whose effects on SS were mediated by FA, as well as 22 pleiotropic SNPs whose effects on FA were mediated by CPD. These SNPs were mainly located in important susceptibility genes for smoking-induced diseases NCAM1 and IREB2. Our findings revealed the role of cerebral WM in the maintenance of the complex addiction and provided potential genetic targets for future research in examining how changes in WM integrity contribute to the nicotine effects on the brain.

A pleiotropic ATM variant (rs1800057 C>G) is associated with risk of multiple cancers

ATM (ataxia-telangiectasia mutated) is an important cell-cycle checkpoint kinase required for cellular response to DNA damage. Activated by DNA double strand breaks, ATM regulates the activities of many downstream proteins involved in various carcinogenic events. Therefore, ATM or its genetic variants may have a pleiotropic effect in cancer development. We conducted a pleiotropic analysis to evaluate associations between genetic variants of ATM and risk of multiple cancers. With genotyping data extracted from previously published genome-wide association studies of various cancers, we performed multivariate logistic regression analysis, followed by a meta-analysis for each cancer site, to identify cancer risk-associated single-nucleotide polymorphisms (SNPs). In the ASSET two-sided analysis, we found that two ATM SNPs were significantly associated with risk of multiple cancers. One tagging SNP (rs1800057 C>G) was associated with risk of multiple cancers (two-sided P=5.27×10 -7). Because ATM rs1800057 is a missense variant, we also explored the intermediate phenotypes through which this variant may confer risk of multiple cancers and identified a possible immune-mediated effect of this variant. Our findings indicate that genetic variants of ATM may have a pleiotropic effect on cancer risk and thus provide an important insight into common mechanisms of carcinogenesis.

Investigating genetic correlation and causality between nicotine dependence and ADHD in a broader psychiatric context

People with attention-deficit/hyperactivity disorder (ADHD) or other psychiatric disorders show high rates of nicotine dependence (ND). This comorbidity might be (partly) explained by shared genetic factors. Genetic correlations between ND and ADHD (or other psychiatric disorders) have not yet been estimated. A significant genetic correlation might indicate genetic overlap, but could also reflect a causal relationship. In the present study we investigated the genetic correlation (with LD score regression analyses) between ND and ADHD, as well as between ND and other major psychiatric conditions (major depressive disorder, schizophrenia, anxiety, bipolar disorder, autism spectrum, anorexia nervosa, and antisocial behavior) based on the summary statistics of large Genome Wide Association studies. We explored the causal nature of the relationship between ND and ADHD using two-sample Mendelian randomization. We found a high genetic correlation between ND and ADHD (r_g  = .53, p = 1.85 × 10^-13 ), and to a lesser extent also between ND-major depressive disorder (r_g  = .42, p = 3.6 × 10^-11 ) and ND-schizophrenia (r_g  = .18, p = 1.1 × 10^-3 ). We did not find evidence for a causal relationship from liability for ADHD to ND (which could be due to a lack of power). The strong genetic correlations might reflect different phenotypic manifestations of (partly) shared underlying genetic vulnerabilities. Combined with the lack of evidence for a causal relationship from liability for ADHD to ND, these findings stress the importance to further investigate the underlying genetic vulnerability explaining co-morbidity in psychiatric disorders.

Shared genetic architecture across psychiatric disorders

Psychiatric disorders overlap substantially at the genetic level, with family-based methods long pointing toward transdiagnostic risk pathways. Psychiatric genomics has progressed rapidly in the last decade, shedding light on the biological makeup of cross-disorder risk at multiple levels of analysis. Over a hundred genetic variants have been identified that affect multiple disorders, with many more to be uncovered as sample sizes continue to grow. Cross-disorder mechanistic studies build on these findings to cluster transdiagnostic variants into meaningful categories, including in what tissues or when in development these variants are expressed. At the upper-most level, methods have been developed to estimate the overall shared genetic signal across pairs of traits (i.e. single-nucleotide polymorphism-based genetic correlations) and subsequently model these relationships to identify overarching, genomic risk factors. These factors can subsequently be associated with external traits (e.g. functional imaging phenotypes) to begin to understand the makeup of these transdiagnostic risk factors. As psychiatric genomic efforts continue to expand, we can begin to gain even greater insight by including more fine-grained phenotypes (i.e. symptom-level data) and explicitly considering the environment. The culmination of these efforts will help to inform bottom-up revisions of our current nosology.

Evolution, Chance, and Aging

Aging has provided fruitful challenges for evolutionary theory, and evolutionary theory has deepened our understanding of aging. A great deal of genetic and molecular data now exists concerning mortality regulation and there is a growing body of knowledge concerning the life histories of diverse species. Assimilating all relevant data into a framework for the evolution of aging promises to significantly advance the field. We propose extensions of some key concepts to provide greater precision when applying these concepts to age-structured contexts. Secondary or byproduct effects of mutations are proposed as an important factor affecting survival patterns, including effects that may operate in small populations subject to genetic drift, widening the possibilities for mutation accumulation and pleiotropy. Molecular and genetic studies have indicated a diverse array of mechanisms that can modify aging and mortality rates, while transcriptome data indicate a high level of tissue and species specificity for genes affected by aging. The diversity of mechanisms and gene effects that can contribute to the pattern of aging in different organisms may mirror the complex evolutionary processes behind aging.

Major Depressive Disorder and Lifestyle: Correlated Genetic Effects in Extended Twin Pedigrees

In recent years, evidence has accumulated with regard to the ubiquity of pleiotropy across the genome, and shared genetic etiology is thought to play a large role in the widespread comorbidity among psychiatric disorders and risk factors. Recent methods investigate pleiotropy by estimating genetic correlation from genome-wide association summary statistics. More comprehensive estimates can be derived from the known relatedness between genetic relatives. Analysis of extended twin pedigree data allows for the estimation of genetic correlation for additive and non-additive genetic effects, as well as a shared household effect. Here we conduct a series of bivariate genetic analyses in extended twin pedigree data on lifetime major depressive disorder (MDD) and three indicators of lifestyle, namely smoking behavior, physical inactivity, and obesity, decomposing phenotypic variance and covariance into genetic and environmental components. We analyze lifetime MDD and lifestyle data in a large multigenerational dataset of 19,496 individuals by variance component analysis in the ‘Mendel’ software. We find genetic correlations for MDD and smoking behavior (r_G = 0.249), physical inactivity (r_G = 0.161), body-mass index (r_G = 0.081), and obesity (r_G = 0.155), which were primarily driven by additive genetic effects. These outcomes provide evidence in favor of a shared genetic etiology between MDD and the lifestyle factors.

Application of Genetic Studies to Flow Cytometry Data and Its Impact on Therapeutic Intervention for Autoimmune Disease

In recent years, systematic genome-wide association studies of quantitative immune cell traits, represented by circulating levels of cell subtypes established by flow cytometry, have revealed numerous association signals, a large fraction of which overlap perfectly with genetic signals associated with autoimmune diseases. By identifying further overlaps with association signals influencing gene expression and cell surface protein levels, it has also been possible, in several cases, to identify causal genes and infer candidate proteins affecting immune cell traits linked to autoimmune disease risk. Overall, these results provide a more detailed picture of how genetic variation affects the human immune system and autoimmune disease risk. They also highlight druggable proteins in the pathogenesis of autoimmune diseases; predict the efficacy and side effects of existing therapies; provide new indications for use for some of them; and optimize the research and development of new, more effective and safer treatments for autoimmune diseases. Here we review the genetic-driven approach that couples systematic multi-parametric flow cytometry with high-resolution genetics and transcriptomics to identify endophenotypes of autoimmune diseases for the development of new therapies.

Pleiotropy data resource as a primer for investigating co-morbidities/multi-morbidities and their role in disease

Most current biomedical and protein research focuses only on a small proportion of genes, which results in a lost opportunity to identify new gene-disease associations and explore new opportunities for therapeutic intervention. The International Mouse Phenotyping Consortium (IMPC) focuses on elucidating gene function at scale for poorly characterized and/or under-studied genes. A key component of the IMPC initiative is the implementation of a broad phenotyping pipeline, which is facilitating the discovery of pleiotropy. Characterizing pleiotropy is essential to identify gene-disease associations, and it is of particular importance when elucidating the genetic causes of syndromic disorders. Here we show how the IMPC is effectively uncovering pleiotropy and how the new mouse models and gene function hypotheses generated by the IMPC are increasing our understanding of the mammalian genome, forming the basis of new research and identifying new gene-disease associations.

Genome-wide shifts in climate-related variation underpin responses to selective breeding in a widespread conifer

Temperate trees originating from warmer localities usually grow faster and acclimate to winter later than trees of the same species from colder localities. However, when trees are selected for faster growth, are climatically adaptive genes and traits affected? Our research demonstrates a simple, sensitive, broadly applicable, and elusive approach to integrating complex polygenic variation into applied environmental management. We show that selective breeding increases allele frequencies of many trait-associated genes, and that alleles that increase most display strong genetic linkage and potential trade-offs among traits. Increasing tree growth while maintaining adaptation is essential to reforestation in a changing climate, but this may be challenging because trait-associated genetic variation that underpins climate adaptation is responsive to selection for tree growth.

Locally adapted temperate tree populations exhibit genetic trade-offs among climate-related traits that can be exacerbated by selective breeding and are challenging to manage under climate change. To inform climatically adaptive forest management, we investigated the genetic architecture and impacts of selective breeding on four climate-related traits in 105 natural and 20 selectively bred lodgepole pine populations from western Canada. Growth, cold injury, growth initiation, and growth cessation phenotypes were tested for associations with 18,600 single-nucleotide polymorphisms (SNPs) in natural populations to identify “positive effect alleles” (PEAs). The effects of artificial selection for faster growth on the frequency of PEAs associated with each trait were quantified in breeding populations from different climates. Substantial shifts in PEA proportions and frequencies were observed across many loci after two generations of selective breeding for height, and responses of phenology-associated PEAs differed strongly among climatic regions. Extensive genetic overlap was evident among traits. Alleles most strongly associated with greater height were often associated with greater cold injury and delayed phenology, although it is unclear whether potential trade-offs arose directly from pleiotropy or indirectly via genetic linkage. Modest variation in multilocus PEA frequencies among populations was associated with large phenotypic differences and strong climatic gradients, providing support for assisted gene flow polices. Relationships among genotypes, phenotypes, and climate in natural populations were maintained or strengthened by selective breeding. However, future adaptive phenotypes and assisted gene flow may be compromised if selective breeding further increases the PEA frequencies of SNPs involved in adaptive trade-offs among climate-related traits.

Favorable pleiotropic loci for fiber yield and quality in upland cotton (Gossypium hirsutum)

Upland cotton (Gossypium hirsutum L.) is an important economic crop for renewable textile fibers. However, the simultaneous improvement of yield and fiber quality in cotton is difficult as the linkage drag. Compared with breaking the linkage drag, identification of the favorable pleiotropic loci on the genome level by genome-wide association study (GWAS) provides a new way to improve the yield and fiber quality simultaneously. In our study restriction-site-associated DNA sequencing (RAD-seq) was used to genotype 316 cotton accessions. Eight major traits in three categories including yield, fiber quality and maturation were investigated in nine environments (3 sites × 3 years). 231 SNPs associated with these eight traits (− log₁₀(P) > 5.27) were identified, located in 27 genomic regions respectively by linkage disequilibrium analysis. Further analysis showed that four genomic regions (the region 1, 6, 8 and 23) held favorable pleiotropic loci and 6 candidate genes were identified. Through genotyping, 14 elite accessions carrying the favorable loci on four pleiotropic regions were identified. These favorable pleiotropic loci and elite genotypes identified in this study will be utilized to improve the yield and fiber quality simultaneously in future cotton breeding.

The contribution of mutation and selection to multivariate quantitative genetic variance in an outbred population of Drosophila serrata

Genetic variance is not equal for all multivariate combinations of traits. This inequality, in which some combinations of traits have abundant genetic variation while others have very little, biases the rate and direction of multivariate phenotypic evolution. However, we still understand little about what causes genetic variance to differ among trait combinations. Here, we investigate the relative roles of mutation and selection in determining the genetic variance of multivariate phenotypes. We accumulated mutations in an outbred population of Drosophila serrata and analyzed wing shape and size traits for over 35,000 flies to simultaneously estimate the additive genetic and additive mutational (co)variances. This experimental design allowed us to gain insight into the phenotypic effects of mutation as they arise and come under selection in naturally outbred populations. Multivariate phenotypes associated with more (less) genetic variance were also associated with more (less) mutational variance, suggesting that differences in mutational input contribute to differences in genetic variance. However, mutational correlations between traits were stronger than genetic correlations, and most mutational variance was associated with only one multivariate trait combination, while genetic variance was relatively more equal across multivariate traits. Therefore, selection is implicated in breaking down trait covariance and resulting in a different pattern of genetic variance among multivariate combinations of traits than that predicted by mutation and drift. Overall, while low mutational input might slow evolution of some multivariate phenotypes, stabilizing selection appears to reduce the strength of evolutionary bias introduced by pleiotropic mutation.

Highly pleiotropic variants of human traits are enriched in genomic regions with strong background selection

Recent studies have shown the ubiquity of pleiotropy for variants affecting human complex traits. These studies also show that rare variants tend to be less pleiotropic than common ones, suggesting that purifying natural selection acts against highly pleiotropic variants of large effect. Here, we investigate the mean frequency, effect size and recombination rate associated with pleiotropic variants, and focus particularly on whether highly pleiotropic variants are enriched in regions with putative strong background selection. We evaluate variants for 41 human traits using data from the NHGRI-EBI GWAS Catalog, as well as data from other three studies. Our results show that variants involving a higher degree of pleiotropy tend to be more common, have larger mean effect sizes, and contribute more to heritability than variants with a lower degree of pleiotropy. This is consistent with the fact that variants of large effect and frequency are more likely detected by GWAS. Using data from four different studies, we also show that more pleiotropic variants are enriched in genome regions with stronger background selection than less pleiotropic variants, suggesting that highly pleiotropic variants are subjected to strong purifying selection. From the above results, we hypothesized that a number of highly pleiotropic variants of low effect/frequency may pass undetected by GWAS.

From QTL to gene: C. elegans facilitates discoveries of the genetic mechanisms underlying natural variation

Although many studies have examined quantitative trait variation across many species, only a small number of genes and thereby molecular mechanisms have been discovered. Without these data, we can only speculate about evolutionary processes that underlie trait variation. Here, we review how quantitative and molecular genetics in the nematode Caenorhabditis elegans led to the discovery and validation of 37 quantitative trait genes over the past 15 years. Using these data, we can start to make inferences about evolution from these quantitative trait genes, including the roles that coding versus noncoding variation, gene family expansion, common versus rare variants, pleiotropy, and epistasis play in trait variation across this species.

Molecular evolution of bumble bee vitellogenin and vitellogenin-like genes

Vitellogenin (Vg), a storage protein, has been significantly studied for its egg yolk precursor role in oviparous animals. Recent studies found that vitellogenin and its Vg-like homologs were fundamentally involved in many other biological processes in social insects such as female caste differences and oxidative stress resilience. In this study, we conducted the first large-scale molecular evolutionary analyses of vitellogenin coding genes (Vg) and Vg-like genes of bumble bees, a primitively eusocial insect belonging to the genus Bombus. We obtained sequences for each of the four genes (Vg, Vg-like-A, Vg-like-B, and Vg-like-C) from 27 bumble bee genomes (nine were newly sequenced in this study), and sequences from the two closest clades of Bombus, including five Apis species and five Tetragonula species. Our molecular evolutionary analyses show that in bumble bee, the conventional Vg experienced strong positive selection, while the Vg-like genes showed overall relaxation of purifying selection. In Apis and Tetragonula; however, all four genes were found under purifying selection. Furthermore, the conventional Vg showed signs of strong positive selection in most subgenera in Bombus, apart from the obligate parasitic subgenus Psithyrus which has no caste differentiation. Together, these results indicate that the conventional Vg, a key pleiotropic gene in social insects, is the most rapidly evolving copy, potentially due to its multiple known social functions for both worker and queen castes. This study shows that concerted evolution and purifying selection shaped the evolution of the Vg gene family following their ancient gene duplication and may be the leading forces behind the evolution of new potential protein function enabling functional social pleiotropy.

Genetically Predicted Serum Iron Status Is Associated with Altered Risk of Systemic Lupus Erythematosus among European Populations

BACKGROUND

Observational epidemiological studies have reported an inconsistent relation between iron status and risk of systemic lupus erythematosus (SLE). Moreover, it remains uncertain whether the observed association is causal or due to confounding or reverse causality.

OBJECTIVES

We aimed to investigate the association between serum iron status and risk of SLE using a 2-sample Mendelian randomization (MR) approach.

METHODS

Genetic instruments for iron status including serum iron, log-transformed ferritin, transferrin saturation, and transferrin were identified from a large-scale genome-wide association study (GWAS) performed by the Genetics of Iron Status Consortium among 48,972 individuals of European ancestry (55% female). Three independent single nucleotide polymorphisms (rs1800562, rs1799945, and rs855791) concordantly related with 4 iron status biomarkers were selected as instrumental variables. Summary statistics of SLE were obtained from a publicly available GWAS of 4036 patients with SLE and 6959 controls of European descent. The MR study was conducted using the inverse-variance weighted (IVW) method, supplemented with MR-Egger regression and simple- and weighted-median methods. Leave-one-out analysis was further performed to test the robustness of our findings. ORs with 95% CIs were calculated.

RESULTS

Genetically predicted iron status was associated with altered risk of SLE, with ORs of 0.79 (95% CI: 0.66, 0.94), 0.54 (95% CI: 0.34, 0.85), 0.82 (95% CI: 0.71, 0.94), and 1.36 (95% CI: 1.06, 1.76) per 1-SD increase in iron, log-transformed ferritin, transferrin saturation, and transferrin using the IVW method, respectively. MR-Egger regression did not indicate potential pleiotropic bias. Sensitivity analyses produced similar findings, suggesting the robustness of the association.

CONCLUSIONS

Our study suggested that high iron status may be associated with a reduced risk of SLE among European populations. Further studies are warranted to elucidate the mechanism underlying the protective role of iron against susceptibility to SLE.

Transhumanist Genetic Enhancement: Creation of a 'New Man' Through Technological Innovation

The transhumanist project of reshaping human beings by promoting their improvement through technological innovations has a broad agenda. This study focuses on the enhancement of the human organism through genetic modification techniques. Transhumanism values and a discussion of their philosophical background provide a framework to understand its ideals. Genetics and ethics are employed to assess the claims of the transhumanist program of human enhancement. A succinct description of central concepts in genetics and an explanation of current techniques to edit the human genome serve to assess the capabilities and limitations of editing techniques. Potential benefits and liabilities of human enhancement through genome editing are discussed to appraise its feasibility. Ethical considerations of genome editing inform a reflection on the implications of introducing heritable changes in the genome of individuals. It is concluded that the transhumanist program is underpinned by a large number of hypotheses rather than by sufficient evidence.

Hidden paths to endless forms most wonderful: parasite-blind diversification of host quality

Evolutionary diversification can occur in allopatry or sympatry, can be driven by selection or unselected, and can be phenotypically manifested immediately or remain latent until manifested in a newly encountered environment. Diversification of host-parasite interactions is frequently studied in the context of intrinsically selective coevolution, but the potential for host-parasite interaction phenotypes to diversify latently during parasite-blind host evolution is rarely considered. Here, we use a social bacterium experimentally adapted to several environments in the absence of phage to analyse allopatric diversification of host quality-the degree to which a host population supports a viral epidemic. Phage-blind evolution reduced host quality overall, with some bacteria becoming completely resistant to growth suppression by phage. Selective-environment differences generated only mild divergence in host quality. However, selective environments nonetheless played a major role in shaping evolution by determining the degree of stochastic diversification among replicate populations within treatments. Ancestral motility genotype was also found to strongly shape patterns of latent host-quality evolution and diversification. These outcomes show that (i) adaptive landscapes can differ in how they constrain stochastic diversification of a latent phenotype and (ii) major effects of selection on biological diversification can be missed by focusing on trait means. Collectively, our findings suggest that latent-phenotype evolution should inform host-parasite evolution theory and that diversification should be conceived broadly to include latent phenotypes.

Complex pleiotropic genetic architecture of evolved heat stress and oxidative stress resistance in the nematode Caenorhabditis remanei

The adaptation of complex organisms to changing environments has been a central question in evolutionary quantitative genetics since its inception. The structure of the genotype-phenotype maps is critical because pleiotropic effects can generate widespread correlated responses to selection and potentially restrict the extent of evolutionary change. In this study, we use experimental evolution to dissect the genetic architecture of natural variation for acute heat stress and oxidative stress response in the nematode Caenorhabiditis remanei. Previous work in the classic model nematode Caenorhabiditis elegans has found that abiotic stress response is controlled by a handful of genes of major effect and that mutations in any one of these genes can have widespread pleiotropic effects on multiple stress response traits. Here, we find that acute heat stress response and acute oxidative response in C. remanei are polygenic, complex traits, with hundreds of genomic regions responding to selection. In contrast to expectation from mutation studies, we find that evolved acute heat stress and acute oxidative stress response for the most part display independent genetic bases. This lack of correlation is reflected at the levels of phenotype, gene expression, and in the genomic response to selection. Thus, while these findings support the general view that rapid adaptation can be generated by changes at hundreds to thousands of sites in the genome, the architecture of segregating variation is likely to be determined by the pleiotropic structure of the underlying genetic networks.

Immune cells lacking Y chromosome show dysregulation of autosomal gene expression

Epidemiological investigations show that mosaic loss of chromosome Y (LOY) in leukocytes is associated with earlier mortality and morbidity from many diseases in men. LOY is the most common acquired mutation and is associated with aberrant clonal expansion of cells, yet it remains unclear whether this mosaicism exerts a direct physiological effect. We studied DNA and RNA from leukocytes in sorted- and single-cells in vivo and in vitro. DNA analyses of sorted cells showed that men diagnosed with Alzheimer's disease was primarily affected with LOY in NK cells whereas prostate cancer patients more frequently displayed LOY in CD4 + T cells and granulocytes. Moreover, bulk and single-cell RNA sequencing in leukocytes allowed scoring of LOY from mRNA data and confirmed considerable variation in the rate of LOY across individuals and cell types. LOY-associated transcriptional effect (LATE) was observed in ~ 500 autosomal genes showing dysregulation in leukocytes with LOY. The fraction of LATE genes within specific cell types was substantially larger than the fraction of LATE genes shared between different subsets of leukocytes, suggesting that LOY might have pleiotropic effects. LATE genes are involved in immune functions but also encode proteins with roles in other diverse biological processes. Our findings highlight a surprisingly broad role for chromosome Y, challenging the view of it as a "genetic wasteland", and support the hypothesis that altered immune function in leukocytes could be a mechanism linking LOY to increased risk for disease.