Genotype And Phenotype Research Articles

A population-genomic approach for estimating selection on polygenic traits in heterogeneous environments.

Strong selection can cause rapid evolutionary change, but temporal fluctuations in the form, direction and intensity of selection can limit net evolutionary change over longer time periods. Fluctuating selection could affect molecular diversity levels and the evolution of plasticity and ecological specialization. Nonetheless, this phenomenon remains understudied, in part because of analytical limitations and the general difficulty of detecting selection that does not occur in a consistent manner. Herein, I fill this analytical gap by presenting an approximate Bayesian computation (ABC) method to detect and quantify fluctuating selection on polygenic traits from population genomic time-series data. I propose a model for environment-dependent phenotypic selection. The evolutionary genetic consequences of selection are then modelled based on a genotype-phenotype map. Using simulations, I show that the proposed method generates accurate and precise estimates of selection when the generative model for the data is similar to the model assumed by the method. The performance of the method when applied to an evolve-and-resequence study of host adaptation in the cowpea seed beetle (Callosobruchusmaculatus) was more idiosyncratic and depended on specific analytical choices. Despite some limitations, these results suggest the proposed method provides a powerful approach to connect the causes of (variable) selection to traits and genome-wide patterns of evolution. Documentation and open-source computer software (fsabc) implementing this method are available from github (https://github.com/zgompert/fsabc.git).

High-resolution genome-wide association study pinpoints metal transporter and chelator genes involved in the genetic control of element levels in maize grain.

Despite its importance to plant function and human health, the genetics underpinning element levels in maize grain remain largely unknown. Through a genome-wide association study in the maize Ames panel of nearly 2,000 inbred lines that was imputed with ∼7.7 million SNP markers, we investigated the genetic basis of natural variation for the concentration of 11 elements in grain. Novel associations were detected for the metal transporter genes rte2 (rotten ear2) and irt1 (iron-regulated transporter1) with boron and nickel, respectively. We also further resolved loci that were previously found to be associated with one or more of five elements (copper, iron, manganese, molybdenum, and/or zinc), with two metal chelator and five metal transporter candidate causal genes identified. The nas5 (nicotianamine synthase5) gene involved in the synthesis of nicotianamine, a metal chelator, was found associated with both zinc and iron and suggests a common genetic basis controlling the accumulation of these two metals in the grain. Furthermore, moderate predictive abilities were obtained for the 11 elemental grain phenotypes with two whole-genome prediction models: Bayesian Ridge Regression (0.33–0.51) and BayesB (0.33–0.53). Of the two models, BayesB, with its greater emphasis on large-effect loci, showed ∼4–10% higher predictive abilities for nickel, molybdenum, and copper. Altogether, our findings contribute to an improved genotype-phenotype map for grain element accumulation in maize.

Directed Evolution of Microbial Communities.

Directed evolution is a form of artificial selection that has been used for decades to find biomolecules and organisms with new or enhanced functional traits. Directed evolution can be conceptualized as a guided exploration of the genotype-phenotype map, where genetic variants with desirable phenotypes are first selected and then mutagenized to search the genotype space for an even better mutant. In recent years, the idea of applying artificial selection to microbial communities has gained momentum. In this article, we review the main limitations of artificial selection when applied to large and diverse collectives of asexually dividing microbes and discuss how the tools of directed evolution may be deployed to engineer communities from the top down. We conceptualize directed evolution of microbial communities as a guided exploration of an ecological structure-function landscape and propose practical guidelines for navigating these ecological landscapes.

Nuclease-Deficient Clustered Regularly Interspaced Short Palindromic Repeat-Based Approaches for In Vitro and In Vivo Gene Activation.

Clustered regularly interspaced short palindromic repeat (CRISPR)-based technology has been adapted to achieve a wide range of genome modifications, including transcription regulation. The focus of this review is on the application of CRISPR-based platforms such as nuclease-deficient Cas9 and Cas12a, to achieve targeted gene activation. We review studies to date that have used CRISPR-based activation technology for the elucidation of biological mechanism and disease correction, as well as its application in genetic screens as a powerful tool for high-throughput genotype-phenotype mapping. In addition to our synthesis and critical analysis of published studies, we explore key considerations for the potential clinical translation of CRISPR-based activation technology.

On the evolution and development of morphological complexity: A view from gene regulatory networks.

How does morphological complexity evolve? This study suggests that the likelihood of mutations increasing phenotypic complexity becomes smaller when the phenotype itself is complex. In addition, the complexity of the genotype-phenotype map (GPM) also increases with the phenotypic complexity. We show that complex GPMs and the above mutational asymmetry are inevitable consequences of how genes need to be wired in order to build complex and robust phenotypes during development.We randomly wired genes and cell behaviors into networks in EmbryoMaker. EmbryoMaker is a mathematical model of development that can simulate any gene network, all animal cell behaviors (division, adhesion, apoptosis, etc.), cell signaling, cell and tissues biophysics, and the regulation of those behaviors by gene products. Through EmbryoMaker we simulated how each random network regulates development and the resulting morphology (i.e. a specific distribution of cells and gene expression in 3D). This way we obtained a zoo of possible 3D morphologies. Real gene networks are not random, but a random search allows a relatively unbiased exploration of what is needed to develop complex robust morphologies. Compared to the networks leading to simple morphologies, the networks leading to complex morphologies have the following in common: 1) They are rarer; 2) They need to be finely tuned; 3) Mutations in them tend to decrease morphological complexity; 4) They are less robust to noise; and 5) They have more complex GPMs. These results imply that, when complexity evolves, it does so at a progressively decreasing rate over generations. This is because as morphological complexity increases, the likelihood of mutations increasing complexity decreases, morphologies become less robust to noise, and the GPM becomes more complex. We find some properties in common, but also some important differences, with non-developmental GPM models (e.g. RNA, protein and gene networks in single cells).

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.

Genotype and phenotype studies of Lowe syndrome in three families in Taiwan

Genotype and phenotype studies of Lowe syndrome in three families in Taiwan

Prediction of complex phenotypes using the Drosophila melanogaster metabolome.

Understanding the genotype-phenotype map and how variation at different levels of biological organization is associated are central topics in modern biology. Fast developments in sequencing technologies and other molecular omic tools enable researchers to obtain detailed information on variation at DNA level and on intermediate endophenotypes, such as RNA, proteins and metabolites. This can facilitate our understanding of the link between genotypes and molecular and functional organismal phenotypes. Here, we use the Drosophila melanogaster Genetic Reference Panel and nuclear magnetic resonance (NMR) metabolomics to investigate the ability of the metabolome to predict organismal phenotypes. We performed NMR metabolomics on four replicate pools of male flies from each of 170 different isogenic lines. Our results show that metabolite profiles are variable among the investigated lines and that this variation is highly heritable. Second, we identify genes associated with metabolome variation. Third, using the metabolome gave better prediction accuracies than genomic information for four of five quantitative traits analyzed. Our comprehensive characterization of population-scale diversity of metabolomes and its genetic basis illustrates that metabolites have large potential as predictors of organismal phenotypes. This finding is of great importance, e.g., in human medicine, evolutionary biology and animal and plant breeding.

An integrated clinical and molecular study of a cohort of Turkish patients with Marfan syndrome harboring known and novel FBN1 variants.

Marfan syndrome (MFS) is an autosomal dominant genetic condition that mainly affects connective tissue in many parts of the body. Cardinal manifestations involve the ocular, skeletal, and cardiovascular systems. The diagnosis of MFS relies on the revised Ghent criteria, outlined by international expert opinion to facilitate accurate recognition of this syndrome as well as to improve patient management and counseling. However, it may not always be possible to make a definitive diagnosis according to these criteria in each patient and thus molecular confirmation is necessary in subjects with suspected MFS. This debilitating, if not fatal, disorder is caused by mutations in FBN1, which encodes a major constitutive element of extracellular microfibrils. Here, we present a detailed clinical and molecular analysis of 76 Turkish patients with definitive or suspected MFS diagnosed at our center between 2014 and 2019. We were able to identify a total of 51 different FBN1 variants in our cohort, 31 of which have previously been reported in the relevant scientific literature. The remaining 20 variants have not been documented to date. In one patient, we detected a large deletion including the entire FBN1 gene using the array CGH approach. Currently, there are very few studies on the genotype-phenotype correlation of patients with MFS, and no clear genotype-phenotype maps for MFS have been constructed so far, except for some cases. We believe that our findings will make a rich and peculiar contribution to the elusive genotype-phenotype relationship in MFS, especially in this large and populous ethnic group.

The Combined Analysis of Pleiotropy and Epistasis (CAPE).

Epistasis, or gene-gene interaction, contributes substantially to trait variation in organisms ranging from yeast to humans, and modeling epistasis directly is critical to understanding the genotype-phenotype map. However, inference of genetic interactions is challenging compared to inference of individual allele effects due to low statistical power. Furthermore, genetic interactions can appear inconsistent across different quantitative traits, presenting a challenge for the interpretation of detected interactions. Here we present a method called the Combined Analysis of Pleiotropy and Epistasis (CAPE) that combines information across multiple quantitative traits to infer directed epistatic interactions. By combining information across multiple traits, CAPE not only increases power to detect genetic interactions but also interprets these interactions across traits to identify a single interaction that is consistent across all observed data. This method generates informative, interpretable interaction networks that explain how variants interact with each other to influence groups of related traits. This method could potentially be used to link genetic variants to gene expression, physiological endophenotypes, and higher-level disease traits.

Analyzing High-Order Epistasis from Genotype-Phenotype Maps Using 'Epistasis' Package.

Epistasis is the phenomenon about the interactions between genes, leading to complex phenotypic effects. The interactions between three or more mutations called "high-order epistasis" aroused significant interests in recent studies. However, there are still debates for analysis of high-order epistasis due to the non-linear model complexity and statistical artifacts. A recent "epistasis" Python package was therefore developed to characterize high-order epistasis by estimating non-linear scaling for mutation effects to extract high-order epistasis using linear models. This method successfully discovered statistically significant high-order epistasis on several real genotype-phenotype maps. We provided a concise and step-by-step guide to apply the "epistasis" by reproducing the high-order epistasis discoveries on real genotype-phenotype data using the latest API of the package.

Correction to 'Using small samples to estimate neutral component size and robustness in the genotype-phenotype map of RNA secondary structure'.

Correction to 'Using small samples to estimate neutral component size and robustness in the genotype-phenotype map of RNA secondary structure'.

A Multiple-Trait Bayesian Variable Selection Regression Method for Integrating Phenotypic Causal Networks in Genome-Wide Association Studies

Bayesian regression methods that incorporate different mixture priors for marker effects are used in multi-trait genomic prediction. These methods can also be extended to genome-wide association studies (GWAS). In multiple-trait GWAS, incorporating the underlying causal structures among traits is essential for comprehensively understanding the relationship between genotypes and traits of interest. Therefore, we develop a GWAS methodology, SEM-Bayesian alphabet, which, by applying the structural equation model (SEM), can be used to incorporate causal structures into multi-trait Bayesian regression methods. SEM-Bayesian alphabet provides a more comprehensive understanding of the genotype-phenotype mapping than multi-trait GWAS by performing GWAS based on indirect, direct and overall marker effects. The superior performance of SEM-Bayesian alphabet was demonstrated by comparing its GWAS results with other similar multi-trait GWAS methods on real and simulated data. The software tool JWAS offers open-source routines to perform these analyses.

Deep learning of individual aesthetics

Accurate evaluation of human aesthetic preferences represents a major challenge for creative evolutionary and generative systems research. Prior work has tended to focus on feature measures of the artefact, such as symmetry, complexity and coherence. However, research models from psychology suggest that human aesthetic experiences encapsulate factors beyond the artefact, making accurate computational models very difficult to design. The interactive genetic algorithm circumvents the problem through human-in-the-loop, subjective evaluation of aesthetics, but is limited due to user fatigue and small population sizes. In this paper, we look at how recent advances in deep learning can assist in automating personal aesthetic judgement. Using a leading artist’s computer art dataset, we investigate the relationship between image measures, such as complexity, and human aesthetic evaluation. We use dimension reduction methods to visualise both genotype and phenotype space in order to support the exploration of new territory in a generative system. Convolutional neural networks trained on the artist’s prior aesthetic evaluations are used to suggest new possibilities similar or between known high-quality genotype-phenotype mappings. We integrate this classification and discovery system into a software tool for evolving complex generative art and design.

Neutral components show a hierarchical community structure in the genotype-phenotype map of RNA secondary structure.

Genotype-phenotype (GP) maps describe the relationship between biological sequences and structural or functional outcomes. They can be represented as networks in which genotypes are the nodes, and one-point mutations between them are the edges. The genotypes that map to the same phenotype form subnetworks consisting of one or multiple disjoint connected components-so-called neutral components (NCs). For the GP map of RNA secondary structure, the NCs have been found to exhibit distinctive network features that can affect the dynamical processes taking place on them. Here, we focus on the community structure of RNA secondary structure NCs. Building on previous findings, we introduce a method to reveal the hierarchical community structure solely from the sequence constraints and composition of the genotypes that form a given NC. Thereby, we obtain modularity values similar to common community detection algorithms, which are much more complex. From this knowledge, we endorse a sampling method that allows a fast exploration of the different communities of a given NC. Furthermore, we introduce a way to estimate the community structure from genotype samples, which is useful when an exhaustive analysis of the NC is not feasible, as is the case for longer sequence lengths.

Inferring a complete genotype-phenotype map from a small number of measured phenotypes.

Understanding evolution requires detailed knowledge of genotype-phenotype maps; however, it can be a herculean task to measure every phenotype in a combinatorial map. We have developed a computational strategy to predict the missing phenotypes from an incomplete, combinatorial genotype-phenotype map. As a test case, we used an incomplete genotype-phenotype dataset previously generated for the malaria parasite's 'chloroquine resistance transporter' (PfCRT). Wild-type PfCRT (PfCRT3D7) lacks significant chloroquine (CQ) transport activity, but the introduction of the eight mutations present in the 'Dd2' isoform of PfCRT (PfCRTDd2) enables the protein to transport CQ away from its site of antimalarial action. This gain of a transport function imparts CQ resistance to the parasite. A combinatorial map between PfCRT3D7 and PfCRTDd2 consists of 256 genotypes, of which only 52 have had their CQ transport activities measured through expression in the Xenopus laevis oocyte. We trained a statistical model with these 52 measurements to infer the CQ transport activity for the remaining 204 combinatorial genotypes between PfCRT3D7 and PfCRTDd2. Our best-performing model incorporated a binary classifier, a nonlinear scale, and additive effects for each mutation. The addition of specific pairwise- and high-order-epistatic coefficients decreased the predictive power of the model. We evaluated our predictions by experimentally measuring the CQ transport activities of 24 additional PfCRT genotypes. The R2 value between our predicted and newly-measured phenotypes was 0.90. We then used the model to probe the accessibility of evolutionary trajectories through the map. Approximately 1% of the possible trajectories between PfCRT3D7 and PfCRTDd2 are accessible; however, none of the trajectories entailed eight successive increases in CQ transport activity. These results demonstrate that phenotypes can be inferred with known uncertainty from a partial genotype-phenotype dataset. We also validated our approach against a collection of previously published genotype-phenotype maps. The model therefore appears general and should be applicable to a large number of genotype-phenotype maps.

Genotype networks of 80 quantitative Arabidopsis thaliana phenotypes reveal phenotypic evolvability despite pervasive epistasis.

We study the genotype-phenotype maps of 80 quantitative phenotypes in the model plant Arabidopsis thaliana, by representing the genotypes affecting each phenotype as a genotype network. In such a network, each vertex or node corresponds to an individual’s genotype at all those genomic loci that affect a given phenotype. Two vertices are connected by an edge if the associated genotypes differ in exactly one nucleotide. The 80 genotype networks we analyze are based on data from genome-wide association studies of 199 A. thaliana accessions. They form connected graphs whose topography differs substantially among phenotypes. We focus our analysis on the incidence of epistasis (non-additive interactions among mutations) because a high incidence of epistasis can reduce the accessibility of evolutionary paths towards high or low phenotypic values. We find epistatic interactions in 67 phenotypes, and in 51 phenotypes every pairwise mutant interaction is epistatic. Moreover, we find phenotype-specific differences in the fraction of accessible mutational paths to maximum phenotypic values. However, even though epistasis affects the accessibility of maximum phenotypic values, the relationships between genotypic and phenotypic change of our analyzed phenotypes are sufficiently smooth that some evolutionary paths remain accessible for most phenotypes, even where epistasis is pervasive. The genotype network representation we use can complement existing approaches to understand the genetic architecture of polygenic traits in many different organisms.

Genome-wide mapping and prediction of plant architecture in a sorghum nested association mapping population.

Modifying plant architecture is often necessary for yield improvement and climate adaptation, but we lack understanding of the genotype-phenotype map for plant morphology in sorghum. Here, we use a nested association mapping (NAM) population that captures global allelic diversity of sorghum to characterize the genetics of leaf erectness, leaf width (at two stages), and stem diameter. Recombinant inbred lines (n=2200) were phenotyped in multiple environments (35,200 observations) and joint linkage mapping was performed with ∼93,000 markers. Fifty-four QTL of small to large effect were identified for trait BLUPs (9-16 per trait) each explaining 0.4-4% of variation across the NAM population. While some of these QTL colocalize with sorghum homologs of grass genes (e.g., those involved in transcriptional regulation of hormone synthesis [rice SPINDLY] and transcriptional regulation of development [rice Ideal plant architecture1]), most QTL did not colocalize with an a priori candidate gene (92%). Genomic prediction accuracy was generally high in five-fold cross-validation (0.65-0.83), and varied from low to high in leave-one-family-out cross-validation (0.04-0.61). The findings provide a foundation to identify the molecular basis of architecture variation in sorghum and establish genomic-enabled breeding for improved plant architecture.

ABO blood group genotypes and ventilatory dysfunction in patients with allergic and nonallergic asthma.

Aim ABO blood group genotypes are established as a genetic factor in pathophysiology of various diseases, such as cardiovascular disorders, cancers, infectious diseases and there is rising evidence of their involvement in other conditions. The aim of this study was to determine if ventilatory changes of lung function in asthma, measured by biomarkers/parameters, are connected to certain ABO blood group genotypes in Croatia. Methods A case-control study included 149 patients with asthma and 153 healthy individuals (blood donors). ABO genotyping on five main alleles was performed using PCR-SSP method. All patients had spirometry performed and severity of asthma was estimated. Clinical parameters of spirometry (FEV1, FEV/FVC, PEF), biomarkers FeNO, IgE and pO2 were measured. The χ2 test, Fisher's test, Kruskal-Wallis test and Spearman's correlation coefficients with p<0.05 were used as statistically significant. Results There was no determined statistically significant difference in both ABO genotypes and phenotypes between patient and control groups. Comparison of the lung function in different ABO phenotypes in asthmatic patients also did not show any statistically significant differences in FEV1 values, FEV/FVC ratio or PEF. Statistically significant differences in oxygenation between different ABO blood types have not been noticed (p=0.326). Differences in quantitative values of biomarkers (FeNO and IgE) between different ABO blood phenotypes in patients with asthma were not significant, except for IgE that had marginal values (p=0.074). Conclusion No correlation was found between certain ABO blood group genotypes and parameters/biomarkers of ventilatory dysfunction in patients with allergic and nonallergic asthma.

Services and tools for mammalian genotype and phenotype data integration and comparison

Services and tools for mammalian genotype and phenotype data integration and comparison