Genetic architecture of quantitative traits in mice, flies, and humans

Epistasis and pleiotropy-induced variation for plant breeding

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

Evolutionary genomics of socially polymorphic populations of Pogonomyrmex californicus

BACKGROUND

Social insects vary considerably in their social organization both between and within species. In the California harvester ant, Pogonomyrmex californicus (Buckley 1867), colonies are commonly founded and headed by a single queen (haplometrosis, primary monogyny). However, in some populations in California (USA), unrelated queens cooperate not only during founding (pleometrosis) but also throughout the life of the colony (primary polygyny). The genetic architecture and evolutionary dynamics of this complex social niche polymorphism (haplometrosis vs pleometrosis) have remained unknown.

RESULTS

We provide a first analysis of its genomic basis and evolutionary history using population genomics comparing individuals from a haplometrotic population to those from a pleometrotic population. We discovered a recently evolved (< 200 k years), 8-Mb non-recombining region segregating with the observed social niche polymorphism. This region shares several characteristics with supergenes underlying social polymorphisms in other socially polymorphic ant species. However, we also find remarkable differences from previously described social supergenes. Particularly, four additional genomic regions not in linkage with the supergene show signatures of a selective sweep in the pleometrotic population. Within these regions, we find for example genes crucial for epigenetic regulation via histone modification (chameau) and DNA methylation (Dnmt1).

CONCLUSIONS

Altogether, our results suggest that social morph in this species is a polygenic trait involving a potential young supergene. Further studies targeting haplo- and pleometrotic individuals from a single population are however required to conclusively resolve whether these genetic differences underlie the alternative social phenotypes or have emerged through genetic drift.

Genetic, Environmental, and Stochastic Components of Lifespan Variability: The Drosophila Paradigm

Lifespan is a complex quantitative trait involving genetic and non-genetic factors as well as the peculiarities of ontogenesis. As with all quantitative traits, lifespan shows considerable variation within populations and between individuals. Drosophila, a favourite object of geneticists, has greatly advanced our understanding of how different forms of variability affect lifespan. This review considers the role of heritable genetic variability, phenotypic plasticity and stochastic variability in controlling lifespan in Drosophila melanogaster. We discuss the major historical milestones in the development of the genetic approach to study lifespan, the breeding of long-lived lines, advances in lifespan QTL mapping, the environmental factors that have the greatest influence on lifespan in laboratory maintained flies, and the mechanisms, by which individual development affects longevity. The interplay between approaches to study ageing and lifespan limitation will also be discussed. Particular attention will be paid to the interaction of different types of variability in the control of lifespan.

Pleiotropy, epistasis and the genetic architecture of quantitative traits

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

The evolution of decision-making mechanisms under competing demands

Animals in nature are constantly managing multiple demands, and decisions about how to adjust behavior in response to ecologically relevant demands is critical for fitness. Evidence for behavioral correlations across functional contexts (behavioral syndromes) and growing appreciation for shared proximate substrates of behavior prompts novel questions about the existence of distinct neural, molecular, and genetic mechanisms involved in decision-making. Those proximate mechanisms are likely to be an important target of selection, but little is known about how they evolve, their evolutionary history, or where they harbor genetic variation. Herein I provide a conceptual framework for understanding the evolution of mechanisms for decision-making, highlighting insights on decision-making in humans and model organisms, and sketch an emerging synthesis.

A Functional Single-Nucleotide Polymorphism Upstream of the Collagen Type III Gene Is Associated with Catastrophic Fracture Risk in Thoroughbred Horses

Fractures caused by bone overloading are a leading cause of euthanasia in Thoroughbred racehorses. The risk of fatal fracture has been shown to be influenced by both environmental and genetic factors but, to date, no specific genetic mechanisms underpinning fractures have been identified. In this study, we utilised a genome-wide polygenic risk score to establish an in vitro cell system to study bone gene regulation in horses at high and low genetic risk of fracture. Candidate gene expression analysis revealed differential expression of COL3A1 and STAT1 genes in osteoblasts derived from high- and low-risk horses. Whole-genome sequencing of two fracture cases and two control horses revealed a single-nucleotide polymorphism (SNP) upstream of COL3A1 that was confirmed in a larger cohort to be significantly associated with fractures. Bioinformatics tools predicted that this SNP may impact the binding of the transcription factor SOX11. Gene modulation demonstrated SOX11 is upstream of COL3A1, and the region binds to nuclear proteins. Furthermore, luciferase assays demonstrated that the region containing the SNP has promoter activity. However, the specific effect of the SNP depends on the broader genetic background of the cells and suggests other factors may also be involved in regulating COL3A1 expression. In conclusion, we have identified a novel SNP that is significantly associated with fracture risk and provide new insights into the regulation of the COL3A1 gene.

Human genomic data have different statistical properties than the data of randomised controlled trials

Madole & Harden argue that the Mendelian reshuffling of genes and genomes is analogous to randomised controlled trials. We are not convinced by their arguments. First, their recipe for meeting the demands on randomised experiments is inherently inconsistent. Second, disequilibrium across chromosomes conflicts with their assumption of statistical independence. Third, the genome-wide association study (GWAS) method has many pitfalls, including low repeatability.

The long-term effects of genomic selection: 2. Changes in allele frequencies of causal loci and new mutations

Genetic selection has been applied for many generations in animal, plant, and experimental populations. Selection changes the allelic architecture of traits to create genetic gain. It remains unknown whether the changes in allelic architecture are different for the recently introduced technique of genomic selection compared to traditional selection methods and whether they depend on the genetic architectures of traits. Here, we investigate the allele frequency changes of old and new causal loci under 50 generations of phenotypic, pedigree, and genomic selection, for a trait controlled by either additive, additive and dominance, or additive, dominance, and epistatic effects. Genomic selection resulted in slightly larger and faster changes in allele frequencies of causal loci than pedigree selection. For each locus, allele frequency change per generation was not only influenced by its statistical additive effect but also to a large extent by the linkage phase with other loci and its allele frequency. Selection fixed a large number of loci, and 5 times more unfavorable alleles became fixed with genomic and pedigree selection than with phenotypic selection. For pedigree selection, this was mainly a result of increased genetic drift, while genetic hitchhiking had a larger effect on genomic selection. When epistasis was present, the average allele frequency change was smaller (∼15% lower), and a lower number of loci became fixed for all selection methods. We conclude that for long-term genetic improvement using genomic selection, it is important to consider hitchhiking and to limit the loss of favorable alleles.

Leveraging the genetic correlation between traits improves the detection of epistasis in genome-wide association studies

Epistasis, commonly defined as the interaction between genetic loci, is known to play an important role in the phenotypic variation of complex traits. As a result, many statistical methods have been developed to identify genetic variants that are involved in epistasis, and nearly all of these approaches carry out this task by focusing on analyzing one trait at a time. Previous studies have shown that jointly modeling multiple phenotypes can often dramatically increase statistical power for association mapping. In this study, we present the "multivariate MArginal ePIstasis Test" (mvMAPIT)-a multioutcome generalization of a recently proposed epistatic detection method which seeks to detect marginal epistasis or the combined pairwise interaction effects between a given variant and all other variants. By searching for marginal epistatic effects, one can identify genetic variants that are involved in epistasis without the need to identify the exact partners with which the variants interact-thus, potentially alleviating much of the statistical and computational burden associated with conventional explicit search-based methods. Our proposed mvMAPIT builds upon this strategy by taking advantage of correlation structure between traits to improve the identification of variants involved in epistasis. We formulate mvMAPIT as a multivariate linear mixed model and develop a multitrait variance component estimation algorithm for efficient parameter inference and P-value computation. Together with reasonable model approximations, our proposed approach is scalable to moderately sized genome-wide association studies. With simulations, we illustrate the benefits of mvMAPIT over univariate (or single-trait) epistatic mapping strategies. We also apply mvMAPIT framework to protein sequence data from two broadly neutralizing anti-influenza antibodies and approximately 2,000 heterogeneous stock of mice from the Wellcome Trust Centre for Human Genetics. The mvMAPIT R package can be downloaded at https://github.com/lcrawlab/mvMAPIT.

The importance of reaction norms in dietary restriction and ageing research

Ageing research has progressed rapidly through our ability to modulate the ageing process. Pharmacological and dietary treatments can increase lifespan and have been instrumental in our understanding of the mechanisms of ageing. Recently, several studies have reported genetic variance in response to these anti-ageing interventions, questioning their universal application and making a case for personalised medicine in our field. As an extension of these findings the response to dietary restriction was found to not be repeatable when the same genetic mouse lines were retested. We show here that this effect is more widespread with the response to dietary restriction also showing low repeatability across genetic lines in the fly (Drosophila melanogaster). We further argue that variation in reaction norms, the relationship between dose and response, can explain such conflicting findings in our field. We simulate genetic variance in reaction norms and show that such variation can: 1) lead to over- or under-estimation of treatment responses, 2) dampen the response measured if a genetically heterogeneous population is studied, and 3) illustrate that genotype-by-dose-by-environment interactions can lead to low repeatability of DR and potentially other anti-ageing interventions. We suggest that putting experimental biology and personalised geroscience in a reaction norm framework will aid progress in ageing research.

Pierre CL, Bimschleger H, Gao J, Cheng R, Mohammadi P, Chen H, Palmer AA, Polesskaya O, Jablonski MM. Genome-wide association study finds multiple loci associated with intraocular pressure in HS rats

Elevated intraocular pressure (IOP) is influenced by environmental and genetic factors. Increased IOP is a major risk factor for most types of glaucoma, including primary open angle glaucoma (POAG). Investigating the genetic basis of IOP may lead to a better understanding of the molecular mechanisms of POAG. The goal of this study was to identify genetic loci involved in regulating IOP using outbred heterogeneous stock (HS) rats. HS rats are a multigenerational outbred population derived from eight inbred strains that have been fully sequenced. This population is ideal for a genome-wide association study (GWAS) owing to the accumulated recombinations among well-defined haplotypes, the relatively high allele frequencies, the accessibility to a large collection of tissue samples, and the large allelic effect size compared to human studies. Both male and female HS rats (N = 1,812) were used in the study. Genotyping-by-sequencing was used to obtain ∼3.5 million single nucleotide polymorphisms (SNP) from each individual. SNP heritability for IOP in HS rats was 0.32, which agrees with other studies. We performed a GWAS for the IOP phenotype using a linear mixed model and used permutation to determine a genome-wide significance threshold. We identified three genome-wide significant loci for IOP on chromosomes 1, 5, and 16. Next, we sequenced the mRNA of 51 whole eye samples to find cis-eQTLs to aid in identification of candidate genes. We report 5 candidate genes within those loci: Tyr, Ctsc, Plekhf2, Ndufaf6 and Angpt2. Tyr, Ndufaf6 and Angpt2 genes have been previously implicated by human GWAS of IOP-related conditions. Ctsc and Plekhf2 genes represent novel findings that may provide new insight into the molecular basis of IOP. This study highlights the efficacy of HS rats for investigating the genetics of elevated IOP and identifying potential candidate genes for future functional testing.

Genetic architecture of natural variation of cardiac performance from flies to humans

Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP). Genome-wide associations studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict transcription factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart-specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans.

Genomic changes underlying repeated niche shifts in an adaptive radiation

In adaptive radiations, single lineages rapidly diversify by adapting to many new niches. Little is known yet about the genomic mechanisms involved, that is, the source of genetic variation or genomic architecture facilitating or constraining adaptive radiation. Here, we investigate genomic changes associated with repeated invasion of many different freshwater niches by threespine stickleback in the Haida Gwaii archipelago, Canada, by resequencing single genomes from one marine and 28 freshwater populations. We find 89 likely targets of parallel selection in the genome that are enriched for old standing genetic variation. In contrast to theoretical expectations, their genomic architecture is highly dispersed with little clustering. Candidate genes and genotype-environment correlations match the three major environmental axes predation regime, light environment, and ecosystem size. In a niche space with these three dimensions, we find that the more divergent a new niche from the ancestral marine habitat, the more loci show signatures of parallel selection. Our findings suggest that the genomic architecture of parallel adaptation in adaptive radiation depends on the steepness of ecological gradients and the dimensionality of the niche space.

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’.

Feasibility and application of polygenic score analysis to the morphology of human-induced pluripotent stem cells

Genome-wide association studies have identified thousands of significant associations between genetic variants and complex traits. Inferring biological insights from these associations has been challenging. One approach attempted has been to examine the effects of individual variants in cellular models. Here, I demonstrate the feasibility of examining the aggregate effect of many variants on cellular phenotypes. I examine the effects of polygenic scores for cross-psychiatric disorder risk, schizophrenia, body mass index and height on cellular morphology, using 1.5 million induced pluripotent stem cells (iPSC) from 60 European-ancestry donors from the Human iPSC Initiative dataset. I show that measuring multiple cells per donor provides sufficient power for polygenic score analyses, and that cross-psychiatric disorder risk is associated with cell area (p = 0.004). Combined with emerging methods of high-throughput iPSC phenotyping, cellular polygenic scoring is a promising method for understanding potential biological effects of the polygenic component of complex traits.

Challenging Sustainable and Innovative Technologies in Cheese Production: A Review

It is well known that cheese yield and quality are affected by animal genetics, milk quality (chemical, physical, and microbiological), production technology, and the type of rennet and dairy cultures used in production. Major differences in the same type of cheese (i.e., hard cheese) are caused by the rennet and dairy cultures, which affect the ripening process. This review aims to explore current technological advancements in animal genetics, methods for the isolation and production of rennet and dairy cultures, along with possible applications of microencapsulation in rennet and dairy culture production, as well as the challenge posed to current dairy technologies by the preservation of biodiversity. Based on the reviewed scientific literature, it can be concluded that innovative approaches and the described techniques can significantly improve cheese production.

The long-term effects of genomic selection: 1. Response to selection, additive genetic variance, and genetic architecture

Background

Genomic selection has revolutionized genetic improvement in animals and plants, but little is known about its long-term effects. Here, we investigated the long-term effects of genomic selection on response to selection, genetic variance, and the genetic architecture of traits using stochastic simulations. We defined the genetic architecture as the set of causal loci underlying each trait, their allele frequencies, and their statistical additive effects. We simulated a livestock population under 50 generations of phenotypic, pedigree, or genomic selection for a single trait, controlled by either only additive, additive and dominance, or additive, dominance, and epistatic effects. The simulated epistasis was based on yeast data.

Results

Short-term response was always greatest with genomic selection, while response after 50 generations was greater with phenotypic selection than with genomic selection when epistasis was present, and was always greater than with pedigree selection. This was mainly because loss of genetic variance and of segregating loci was much greater with genomic and pedigree selection than with phenotypic selection. Compared to pedigree selection, selection response was always greater with genomic selection. Pedigree and genomic selection lost a similar amount of genetic variance after 50 generations of selection, but genomic selection maintained more segregating loci, which on average had lower minor allele frequencies than with pedigree selection. Based on this result, genomic selection is expected to better maintain genetic gain after 50 generations than pedigree selection. The amount of change in the genetic architecture of traits was considerable across generations and was similar for genomic and pedigree selection, but slightly less for phenotypic selection. Presence of epistasis resulted in smaller changes in allele frequencies and less fixation of causal loci, but resulted in substantial changes in statistical additive effects across generations.

Conclusions

Our results show that genomic selection outperforms pedigree selection in terms of long-term genetic gain, but results in a similar reduction of genetic variance. The genetic architecture of traits changed considerably across generations, especially under selection and when non-additive effects were present. In conclusion, non-additive effects had a substantial impact on the accuracy of selection and long-term response to selection, especially when selection was accurate.

Multiple epistatic DNA variants in a single gene affect gene expression in trans

DNA variants that alter gene expression in trans are important sources of phenotypic variation. Nevertheless, the identity of trans-acting variants remains poorly understood. Single causal variants in several genes have been reported to affect the expression of numerous distant genes in trans. Whether these simple molecular architectures are representative of trans-acting variation is unknown. Here, we studied the large RAS signaling regulator gene IRA2, which contains variants with extensive trans-acting effects on gene expression in the yeast Saccharomyces cerevisiae. We used systematic CRISPR-based genome engineering and a sensitive phenotyping strategy to dissect causal variants to the nucleotide level. In contrast to the simple molecular architectures known so far, IRA2 contained at least seven causal nonsynonymous variants. The effects of these variants were modulated by nonadditive, epistatic interactions. Two variants at the 5'-end affected gene expression and growth only when combined with a third variant that also had no effect in isolation. Our findings indicate that the molecular basis of trans-acting genetic variation may be considerably more complex than previously appreciated.

Genome-Wide Association Study of Soybean Germplasm Derived From Canadian × Chinese Crosses to Mine for Novel Alleles to Improve Seed Yield and Seed Quality Traits

Genome-wide association study (GWAS) has emerged in the past decade as a viable tool for identifying beneficial alleles from a genomic diversity panel. In an ongoing effort to improve soybean [Glycine max (L.) Merr.], which is the third largest field crop in Canada, a GWAS was conducted to identify novel alleles underlying seed yield and seed quality and agronomic traits. The genomic panel consisted of 200 genotypes including lines derived from several generations of bi-parental crosses between modern Canadian × Chinese cultivars (CD-CH). The genomic diversity panel was field evaluated at two field locations in Ontario in 2019 and 2020. Genotyping-by-sequencing (GBS) was conducted and yielded almost 32 K high-quality SNPs. GWAS was conducted using Fixed and random model Circulating Probability Unification (FarmCPU) model on the following traits: seed yield, seed protein concentration, seed oil concentration, plant height, 100 seed weight, days to maturity, and lodging score that allowed to identify five QTL regions controlling seed yield and seed oil and protein content. A candidate gene search identified a putative gene for each of the three traits. The results of this GWAS study provide insight into potentially valuable genetic resources residing in Chinese modern cultivars that breeders may use to further improve soybean seed yield and seed quality traits.

Genetic Bases of Complex Traits: From Quantitative Trait Loci to Prediction

Conceived as a general introduction to the book, this chapter is a reminder of the core concepts of genetic mapping and molecular marker-based prediction. It provides an overview of the principles and the evolution of methods for mapping the variation of complex traits, and methods for QTL-based prediction of human disease risk and animal and plant breeding value. The principles of linkage-based and linkage disequilibrium-based QTL mapping methods are described in the context of the simplest, single-marker, methods. Methodological evolutions are analysed in relation with their ability to account for the complexity of the genotype-phenotype relations. Main characteristics of the genetic architecture of complex traits, drawn from QTL mapping works using large populations of unrelated individuals, are presented. Methods combining marker-QTL association data into polygenic risk score that captures part of an individual's susceptibility to complex diseases are reviewed. Principles of best linear mixed model-based prediction of breeding value in animal- and plant-breeding programs using phenotypic and pedigree data, are summarized and methods for moving from BLUP to marker-QTL BLUP are presented. Factors influencing the additional genetic progress achieved by using molecular data and rules for their optimization are discussed.

CHIT: an allele-specific method for testing the association between molecular quantitative traits and phenotype-genotype interaction

MOTIVATION

Allele-specific differences in molecular traits can be obtained from next-generation sequencing data and could potentially improve testing power, but such information is usually overlooked in association studies. Furthermore, the variation of molecular quantitative traits (e.g. gene expression) could result from the interaction effect of genotypes and phenotypes, but it is challenging to identify such interaction signals in complex disease studies in humans due to small genetic effect sizes and/or small sample sizes.

RESULTS

We develop a novel statistical method, the combined haplotype interaction test (CHIT), which tests for association between molecular quantitative traits and phenotype-genotype interactions by modeling the total read counts and allele-specific reads in a target region. CHIT can be used as a supplementary analysis to the regular linear interaction regression. In our simulations, CHIT obtains non-inflated type I error rates, and it has higher power than a standard interaction quantitative trait locus approach based on linear regression models. Finally, we illustrate CHIT by testing associations between gene expression obtained by RNA-seq and the interaction of SNPs and atopy status from a study of childhood asthma in Puerto Ricans, and results demonstrate that CHIT could be more powerful than a standard linear interaction expression quantitative trait loci approach.

AVAILABILITY AND IMPLEMENTATION

The CHIT algorithm has been implemented in Python. The source code and documentation are available and can be downloaded from https://github.com/QiYanPitt/CHIT.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pleiotropy or linkage? Their relative contributions to the genetic correlation of quantitative traits and detection by multitrait GWA studies

Genetic correlations between traits may cause correlated responses to selection. Previous models described the conditions under which genetic correlations are expected to be maintained. Selection, mutation, and migration are all proposed to affect genetic correlations, regardless of whether the underlying genetic architecture consists of pleiotropic or tightly linked loci affecting the traits. Here, we investigate the conditions under which pleiotropy and linkage have different effects on the genetic correlations between traits by explicitly modeling multiple genetic architectures to look at the effects of selection strength, degree of correlational selection, mutation rate, mutational variance, recombination rate, and migration rate. We show that at mutation-selection(-migration) balance, mutation rates differentially affect the equilibrium levels of genetic correlation when architectures are composed of pairs of physically linked loci compared to architectures of pleiotropic loci. Even when there is perfect linkage (no recombination within pairs of linked loci), a lower genetic correlation is maintained than with pleiotropy, with a lower mutation rate leading to a larger decrease. These results imply that the detection of causal loci in multitrait association studies will be affected by the type of underlying architectures, whereby pleiotropic variants are more likely to be underlying multiple detected associations. We also confirm that tighter linkage between nonpleiotropic causal loci maintains higher genetic correlations at the traits and leads to a greater proportion of false positives in association analyses.

DNA methylation in genes associated with the evolution of ageing and disease: A critical review

Ageing is characterised by a physical decline in biological functioning which results in a progressive risk of mortality with time. As a biological phenomenon, it is underpinned by the dysregulation of a myriad of complex processes. Recently, however, ever-increasing evidence has associated epigenetic mechanisms, such as DNA methylation (DNAm) with age-onset pathologies, including cancer, cardiovascular disease, and Alzheimer's disease. These diseases compromise healthspan. Consequently, there is a medical imperative to understand the link between epigenetic ageing, and healthspan. Evolutionary theory provides a unique way to gain new insights into epigenetic ageing and health. This review will: (1) provide a brief overview of the main evolutionary theories of ageing; (2) discuss recent genetic evidence which has revealed alleles that have pleiotropic effects on fitness at different ages in humans; (3) consider the effects of DNAm on pleiotropic alleles, which are associated with age related disease; (4) discuss how age related DNAm changes resonate with the mutation accumulation, disposable soma and programmed theories of ageing; (5) discuss how DNAm changes associated with caloric restriction intersect with the evolution of ageing; and (6) conclude by discussing how evolutionary theory can be used to inform investigations which quantify age-related DNAm changes which are linked to age onset pathology.

Experimental evaluation of genomic selection prediction for rust resistance in sugarcane

The total sugarcane (Saccharum L.) production has increased worldwide; however, the rate of growth is lower compared with other major crops, mainly due to a plateauing of genetic gain. Genomic selection (GS) has proven to substantially increase the rate of genetic gain in many crops. To investigate the utility of GS in future sugarcane breeding, a field trial was conducted using 432 sugarcane clones using an augmented design with two replications. Two major diseases in sugarcane, brown and orange rust (BR and OR), were screened artificially using whorl inoculation method in the field over two crop cycles. The genotypic data were generated through target enrichment sequencing technologies. After filtering, a set of 8,825 single nucleotide polymorphic markers were used to assess the prediction accuracy of multiple GS models. Using fivefold cross-validation, we observed GS prediction accuracies for BR and OR that ranged from 0.28 to 0.43 and 0.13 to 0.29, respectively, across two crop cycles and combined cycles. The prediction ability further improved by including a known major gene for resistance to BR as a fixed effect in the GS model. It also substantially reduced the minimum number of markers and training population size required for GS. The nonparametric GS models outperformed the parametric GS suggesting that nonadditive genetic effects could contribute genomic sources underlying BR and OR. This study demonstrated that GS could potentially predict the genomic estimated breeding value for selecting the desired germplasm for sugarcane breeding for disease resistance.

Comparative Genomic Analysis of Quantitative Trait Loci Associated With Micronutrient Contents, Grain Quality, and Agronomic Traits in Wheat (Triticum aestivum L.)

Comparative genomics and meta-quantitative trait loci (MQTLs) analysis are important tools for the identification of reliable and stable QTLs and functional genes controlling quantitative traits. We conducted a meta-analysis to identify the most stable QTLs for grain yield (GY), grain quality traits, and micronutrient contents in wheat. A total of 735 QTLs retrieved from 27 independent mapping populations reported in the last 13 years were used for the meta-analysis. The results showed that 449 QTLs were successfully projected onto the genetic consensus map which condensed to 100 MQTLs distributed on wheat chromosomes. This consolidation of MQTLs resulted in a three-fold reduction in the confidence interval (CI) compared with the CI for the initial QTLs. Projection of QTLs revealed that the majority of QTLs and MQTLs were in the non-telomeric regions of chromosomes. The majority of micronutrient MQTLs were located on the A and D genomes. The QTLs of thousand kernel weight (TKW) were frequently associated with QTLs for GY and grain protein content (GPC) with co-localization occurring at 55 and 63%, respectively. The co- localization of QTLs for GY and grain Fe was found to be 52% and for QTLs of grain Fe and Zn, it was found to be 66%. The genomic collinearity within Poaceae allowed us to identify 16 orthologous MQTLs (OrMQTLs) in wheat, rice, and maize. Annotation of promising candidate genes (CGs) located in the genomic intervals of the stable MQTLs indicated that several CGs (e.g., TraesCS2A02G141400, TraesCS3B02G040900, TraesCS4D02G323700, TraesCS3B02G077100, and TraesCS4D02G290900) had effects on micronutrients contents, yield, and yield-related traits. The mapping refinements leading to the identification of these CGs provide an opportunity to understand the genetic mechanisms driving quantitative variation for these traits and apply this information for crop improvement programs.

Genomic analysis of ionome-related QTLs in Arabidopsis thaliana

Ionome contributes to maintain cell integrity and acts as cofactors for catalyzing regulatory pathways. Identifying ionome contributing genomic regions provides a practical framework to dissect the genetic architecture of ionomic traits for use in biofortification. Meta-QTL (MQTL) analysis is a robust method to discover stable genomic regions for traits regardless of the genetic background. This study used information of 483 QTLs for ionomic traits identified from 12 populations for MQTL analysis in Arabidopsis thaliana. The selected QTLs were projected onto the newly constructed genetic consensus map and 33 MQTLs distributed on A. thaliana chromosomes were identified. The average confidence interval (CI) of the drafted MQTLs was 1.30 cM, reduced eight folds from a mean CI of 10.88 cM for the original QTLs. Four MQTLs were considered as stable MQTLs over different genetic backgrounds and environments. In parallel to the gene density over the A. thaliana genome, the genomic distribution of MQTLs over the genetic and physical maps indicated the highest density at non- and sub-telomeric chromosomal regions, respectively. Several candidate genes identified in the MQTLs intervals were associated with ion transportation, tolerance, and homeostasis. The genomic context of the identified MQTLs suggested nine chromosomal regions for Zn, Mn, and Fe control. The QTLs for potassium (K) and phosphorus (P) were the most frequently co-located with Zn (78.3%), Mn (76.2%), and Fe (88.2% and 70.6%) QTLs. The current MQTL analysis demonstrates that meta-QTL analysis is cheaper than, and as informative as genome-wide association study (GWAS) in refining the known QTLs.

Epistasis for head morphology in Drosophila melanogaster

Epistasis-gene-gene interaction-is common for mutations with large phenotypic effects in humans and model organisms. Epistasis impacts quantitative genetic models of speciation, response to natural and artificial selection, genetic mapping, and personalized medicine. However, the existence and magnitude of epistasis between alleles with small quantitative phenotypic effects are controversial and difficult to assess. Here, we use the Drosophila melanogaster Genetic Reference Panel of sequenced inbred lines to evaluate the magnitude of naturally occurring epistasis modifying the effects of mutations in jing and inv, two transcription factors that have subtle quantitative effects on head morphology as homozygotes. We find significant epistasis for both mutations and performed single marker genome-wide association analyses to map candidate modifier variants and loci affecting head morphology. A subset of these loci was significantly enriched for a known genetic interaction network, and mutations of the candidate epistatic modifier loci also affect head morphology.

Genome-wide association studies reveal the coordinated regulatory networks underlying photosynthesis and wood formation in Populus

Photosynthesis and wood formation underlie the ability of trees to provide renewable resources and perform ecological functions; however, the genetic basis and regulatory pathways coordinating these two linked processes remain unclear. Here, we used a systems genetics strategy, integrating genome-wide association studies, transcriptomic analyses, and transgenic experiments, to investigate the genetic architecture of photosynthesis and wood properties among 435 unrelated individuals of Populus tomentosa, and unravel the coordinated regulatory networks resulting in two trait categories. We detected 222 significant single-nucleotide polymorphisms, annotated to 177 candidate genes, for 10 traits of photosynthesis and wood properties. Epistasis uncovered 74 epistatic interactions for phenotypes. Strikingly, we deciphered the coordinated regulation patterns of pleiotropic genes underlying phenotypic variations for two trait categories. Furthermore, expression quantitative trait nucleotide mapping and coexpression analysis were integrated to unravel the potential transcriptional regulatory networks of candidate genes coordinating photosynthesis and wood properties. Finally, heterologous expression of two pleiotropic genes, PtoMYB62 and PtoMYB80, in Arabidopsis thaliana demonstrated that they control regulatory networks balancing photosynthesis and stem secondary cell wall components, respectively. Our study provides insights into the regulatory mechanisms coordinating photosynthesis and wood formation in poplar, and should facilitate genetic breeding in trees via molecular design.

The power and promise of genetic mapping from Plasmodium falciparum crosses utilizing human liver-chimeric mice

Genetic crosses are most powerful for linkage analysis when progeny numbers are high, parental alleles segregate evenly and numbers of inbred progeny are minimized. We previously developed a novel genetic crossing platform for the human malaria parasite Plasmodium falciparum, an obligately sexual, hermaphroditic protozoan, using mice carrying human hepatocytes (the human liver-chimeric FRG NOD huHep mouse) as the vertebrate host. We report on two genetic crosses-(1) an allopatric cross between a laboratory-adapted parasite (NF54) of African origin and a recently patient-derived Asian parasite, and (2) a sympatric cross between two recently patient-derived Asian parasites. We generated 144 unique recombinant clones from the two crosses, doubling the number of unique recombinant progeny generated in the previous 30 years. The allopatric African/Asian cross has minimal levels of inbreeding and extreme segregation distortion, while in the sympatric Asian cross, inbred progeny predominate and parental alleles segregate evenly. Using simulations, we demonstrate that these progeny provide the power to map small-effect mutations and epistatic interactions. The segregation distortion in the allopatric cross slightly erodes power to detect linkage in several genome regions. We greatly increase the power and the precision to map biomedically important traits with these new large progeny panels.

Identification of genetic modifiers of lifespan on a high sugar diet in the Drosophila Genetic Reference Panel

Genome-wide association studies (GWAS) have become beneficial in identifying genetic variants underlying susceptibility to various complex diseases and conditions, including obesity. Utilizing the Drosophila Genetic Reference Panel (DGRP), we performed a GWAS of lifespan of 193 genetically distinct lines on a high sugar diet (HSD). The DGRP analysis pipeline determined the most significant lifespan associated polymorphisms were within loci of genes involved in: neural processes, behavior, development, and apoptosis, among other functions. Next, based on the relevance to obesity pathology, and the availability of transgenic RNAi lines targeting the genes we identified, whole-body in vivo knockdown of several candidate genes was performed. We utilized the GAL4-UAS binary expression system to independently validate the impacts of these loci on Drosophila lifespan during HSD. These loci were largely confirmed to affect lifespan in that HSD setting, as well as a normal diet setting. However, we also detected unexpected dietary effects of the HSD, including inconsistent diet effects on lifespan relative to a normal diet and a strong downregulation of feeding quantity.

Positive Selection in Gene Regulatory Factors Suggests Adaptive Pleiotropic Changes During Human Evolution

Gene regulatory factors (GRFs), such as transcription factors, co-factors and histone-modifying enzymes, play many important roles in modifying gene expression in biological processes. They have also been proposed to underlie speciation and adaptation. To investigate potential contributions of GRFs to primate evolution, we analyzed GRF genes in 27 publicly available primate genomes. Genes coding for zinc finger (ZNF) proteins, especially ZNFs with a Krüppel-associated box (KRAB) domain were the most abundant TFs in all genomes. Gene numbers per TF family differed between all species. To detect signs of positive selection in GRF genes we investigated more than 3,000 human GRFs with their more than 70,000 orthologs in 26 non-human primates. We implemented two independent tests for positive selection, the branch-site-model of the PAML suite and aBSREL of the HyPhy suite, focusing on the human and great ape branch. Our workflow included rigorous procedures to reduce the number of false positives: excluding distantly similar orthologs, manual corrections of alignments, and considering only genes and sites detected by both tests for positive selection. Furthermore, we verified the candidate sites for selection by investigating their variation within human and non-human great ape population data. In order to approximately assign a date to positively selected sites in the human lineage, we analyzed archaic human genomes. Our work revealed with high confidence five GRFs that have been positively selected on the human lineage and one GRF that has been positively selected on the great ape lineage. These GRFs are scattered on different chromosomes and have been previously linked to diverse functions. For some of them a role in speciation and/or adaptation can be proposed based on the expression pattern or association with human diseases, but it seems that they all contributed independently to human evolution. Four of the positively selected GRFs are KRAB-ZNF proteins, that induce changes in target genes co-expression and/or through arms race with transposable elements. Since each positively selected GRF contains several sites with evidence for positive selection, we suggest that these GRFs participated pleiotropically to phenotypic adaptations in humans.

Repeatability of Adaptive Radiation Depends on Spatial Scale: Regional Versus Global Replicates of Stickleback in Lake Versus Stream Habitats

The repeatability of adaptive radiation is expected to be scale-dependent, with determinism decreasing as greater spatial separation among "replicates" leads to their increased genetic and ecological independence. Threespine stickleback (Gasterosteus aculeatus) provide an opportunity to test whether this expectation holds for the early stages of adaptive radiation-their diversification in freshwater ecosystems has been replicated many times. To better understand the repeatability of that adaptive radiation, we examined the influence of geographic scale on levels of parallel evolution by quantifying phenotypic and genetic divergence between lake and stream stickleback pairs sampled at regional (Vancouver Island) and global (North America and Europe) scales. We measured phenotypes known to show lake-stream divergence and used reduced representation genome-wide sequencing to estimate genetic divergence. We assessed the scale dependence of parallel evolution by comparing effect sizes from multivariate models and also the direction and magnitude of lake-stream divergence vectors. At the phenotypic level, parallelism was greater at the regional than the global scale. At the genetic level, putative selected loci showed greater lake-stream parallelism at the regional than the global scale. Generally, the level of parallel evolution was low at both scales, except for some key univariate traits. Divergence vectors were often orthogonal, highlighting possible ecological and genetic constraints on parallel evolution at both scales. Overall, our results confirm that the repeatability of adaptive radiation decreases at increasing spatial scales. We suggest that greater environmental heterogeneity at larger scales imposes different selection regimes, thus generating lower repeatability of adaptive radiation at larger spatial scales.

Correlational selection in the age of genomics

Ecologists and evolutionary biologists are well aware that natural and sexual selection do not operate on traits in isolation, but instead act on combinations of traits. This long-recognized and pervasive phenomenon is known as multivariate selection, or-in the particular case where it favours correlations between interacting traits-correlational selection. Despite broad acknowledgement of correlational selection, the relevant theory has often been overlooked in genomic research. Here, we discuss theory and empirical findings from ecological, quantitative genetic and genomic research, linking key insights from different fields. Correlational selection can operate on both discrete trait combinations and quantitative characters, with profound implications for genomic architecture, linkage, pleiotropy, evolvability, modularity, phenotypic integration and phenotypic plasticity. We synthesize current knowledge and discuss promising research approaches that will enable us to understand how correlational selection shapes genomic architecture, thereby linking quantitative genetic approaches with emerging genomic methods. We suggest that research on correlational selection has great potential to integrate multiple fields in evolutionary biology, including developmental and functional biology, ecology, quantitative genetics, phenotypic polymorphisms, hybrid zones and speciation processes.

MeSCoT: The tool for quantitative trait simulation through the mechanistic modelling of genes' regulatory interactions

This work represents a novel mechanistic approach to simulate and study genomic networks with accompanying regulatory interactions and complex mechanisms of quantitative trait formation. The approach implemented in MeSCoT software is conceptually based on the omnigenic genetic model of quantitative (complex) trait, and closely imitates the basic in vivo mechanisms of quantitative trait realization. The software provides a framework to study molecular mechanisms of gene-by-gene and gene-by-environment interactions underlying quantitative trait’s realization and allows detailed mechanistic studies of impact of genetic and phenotypic variance on gene regulation. MeSCoT performs a detailed simulation of genes’ regulatory interactions for variable genomic architectures and generates complete set of transcriptional and translational data together with simulated quantitative trait values. Such data provide opportunities to study, for example, verification of novel statistical methods aiming to integrate intermediate phenotypes together with final phenotype in quantitative genetic analyses or to investigate novel approaches for exploiting gene-by-gene and gene-by-environment interactions.

Selection, drift, and constraint in cypridinid luciferases and the diversification of bioluminescent signals in sea fireflies

Understanding the genetic causes of evolutionary diversification is challenging because differences across species are complex, often involving many genes. However, cases where single or few genetic loci affect a trait that varies dramatically across a radiation of species provide tractable opportunities to understand the genetics of diversification. Here, we begin to explore how diversification of bioluminescent signals across species of cypridinid ostracods ("sea fireflies") was influenced by evolution of a single gene, cypridinid-luciferase. In addition to emission spectra ("colour") of bioluminescence from 21 cypridinid species, we report 13 new c-luciferase genes from de novo transcriptomes, including in vitro assays to confirm function of four of those genes. Our comparative analyses suggest some amino acid sites in c-luciferase evolved under episodic diversifying selection and may be associated with changes in both enzyme kinetics and colour, two enzymatic functions that directly impact the phenotype of bioluminescent signals. The analyses also suggest multiple other amino acid positions in c-luciferase evolved neutrally or under purifying selection, and may have impacted the variation of colour of bioluminescent signals across genera. Previous mutagenesis studies at candidate sites show epistatic interactions, which could constrain the evolution of c-luciferase function. This work provides important steps toward understanding the genetic basis of diversification of behavioural signals across multiple species, suggesting different evolutionary processes act at different times during a radiation of species. These results set the stage for additional mutagenesis studies that could explicitly link selection, drift, and constraint to the evolution of phenotypic diversification.

Cis-regulatory variation in the shavenbaby gene underlies intraspecific phenotypic variation, mirroring interspecific divergence in the same trait

Despite considerable progress in recent decades in dissecting the genetic causes of natural morphological variation, there is limited understanding of how variation within species ultimately contributes to species differences. We have studied patterning of the non-sensory hairs, commonly known as "trichomes," on the dorsal cuticle of first-instar larvae of Drosophila. Most Drosophila species produce a dense lawn of dorsal trichomes, but a subset of these trichomes were lost in D. sechellia and D. ezoana due entirely to regulatory evolution of the shavenbaby (svb) gene. Here, we describe intraspecific variation in dorsal trichome patterns of first-instar larvae of D. virilis that is similar to the trichome pattern variation identified previously between species. We found that a single large effect QTL, which includes svb, explains most of the trichome number difference between two D. virilis strains and that svb expression correlates with the trichome difference between strains. This QTL does not explain the entire difference between strains, implying that additional loci contribute to variation in trichome numbers. Thus, the genetic architecture of intraspecific variation exhibits similarities and differences with interspecific variation that may reflect differences in long-term and short-term evolutionary processes.

Interspecies Comparative Analyses Reveal Distinct Carbohydrate-Responsive Systems among Drosophila Species

During evolution, organisms have acquired variable feeding habits. Some species are nutritional generalists that adapt to various food resources, while others are specialists, feeding on specific resources. However, much remains to be discovered about how generalists adapt to diversified diets. We find that larvae of the generalists Drosophila melanogaster and D. simulans develop on three diets with different nutrient balances, whereas specialists D. sechellia and D. elegans cannot develop on carbohydrate-rich diets. The generalist D. melanogaster downregulates the expression of diverse metabolic genes systemically by transforming growth factor β (TGF-β)/Activin signaling, maintains metabolic homeostasis, and successfully adapts to the diets. In contrast, the specialist D. sechellia expresses those metabolic genes at higher levels and accumulates various metabolites on the carbohydrate-rich diet, culminating in reduced adaptation. Phenotypic similarities and differences strongly suggest that the robust carbohydrate-responsive regulatory systems are evolutionarily retained through genome-environment interactions in the generalists and contribute to their nutritional adaptabilities.

Functional Captures of Multiple Human Quantitative Trait Loci Regulating Blood Pressure with the Use of Orthologs in Genetically Defined Rat Models

BACKGROUND

Most signals from human genome-wide association studies (GWAS) for blood pressure (BP) are single-nucleotide polymorphisms (SNPs). It was unknown if such SNPs can functionally affect BP. Because BP is similar between humans and rodents, unraveling basic mechanisms from rodents can reveal the same BP-modulating mechanisms in humans originating from their common ancestors while overcoming limitations in human epidemiology.

METHODS

For the first time, we used quantitative trait loci (QTLs) from Dahl salt-sensitive (DSS) rats as functional surrogates to capture human BP QTLs.

RESULTS

A total of 107 human GWAS genes may be classified into 2 common pathways of hypertension pathogeneses. Among them, 4 DSS BP QTLs correspond to 4 human GWAS genes. Each of them independently showed a major impact on BP in vivo and thus functional redundancy. BP was altered by each of these 4 QTLs, but human GWAS SNPs marking these QTLs do not exist in the rat. They cannot be responsible for physiological changes in BP caused by these QTLs and are genome signposts marking positions of the QTLs nearby, rather than being QTLs themselves. These SNPs appeared during primate evolution, independently of BP regulation. Because the functional dosage of QTLs, not their gene dose, determined hypertension pathogenesis, a role for the noncoding GWAS SNPs in BP via regulating gene expressions can be discounted.

CONCLUSIONS

The human QTLs may function in a common pathway, with each involved in a different step in the pathway leading to BP control. These results may be conceptually paradigm shifting.

Opinion: Is gene mapping in wild populations useful for understanding and predicting adaptation to global change?

Changing environmental conditions will inevitably alter selection pressures. Over the long term, populations have to adapt to these altered conditions by evolutionary change to avoid extinction. Quantifying the 'evolutionary potential' of populations to predict whether they will be able to adapt fast enough to forecasted changes is crucial to fully assess the threat for biodiversity posed by climate change. Technological advances in sequencing and high-throughput genotyping have now made genomic studies possible in a wide range of species. Such studies, in theory, allow an unprecedented understanding of the genomics of ecologically relevant traits and thereby a detailed assessment of the population's evolutionary potential. Aimed at a wider audience than only evolutionary geneticists, this paper gives an overview of how gene-mapping studies have contributed to our understanding and prediction of evolutionary adaptations to climate change, identifies potential reasons why their contribution to understanding adaptation to climate change may remain limited, and highlights approaches to study and predict climate change adaptation that may be more promising, at least in the medium term.

Identification and epistasis analysis of quantitative trait loci for zeaxanthin concentration in maize kernel across different generations and environments

Zeaxanthin, a natural fat-soluble pigment, not only increases plant resistance, but also has vital significance for human health. However, quantitative trait loci (QTL) and the epistatic effects of zeaxanthin concentration in maize kernel have not been well studied. To identify QTLs and analyse the epistatic effects of zeaxanthin concentration in maize kernel, two sets of segregating generations derived from the cross between HuangC (a high zeaxanthin concentration inbred line) and Rezi1 (a low zeaxanthin concentration inbred line) were evaluated in three different environments. One major-effect QTL, qZea6a, explains 41.4-71.4% of the phenotypic variation and two QTLs, qZea4a and qZea3a, show LOD > 3 for zeaxanthin concentration detected over two generations and three different environments. Four of the ten QTL pairs show epistatic effects, explaining 7.34-14.3% of the phenotypic variance. Furthermore, additivity was the major allelic action at zeaxanthin concentration QTLs located in F₂ and F_2:3 populations and plants with homozygous HuangC alleles have a strong genetic ability in enhancing zeaxanthin concentration in maize kernel. These results will contribute to understanding these complex loci better and provide awareness about zeaxanthin concentration to maize breeders and scientists involved in maize research.

Modularity/non-cumulativity of quantitative trait loci on blood pressure

Large numbers of quantitative trait loci (QTLs) for blood pressure (BP) exist and have long been thought to function by accumulating their individual miniscule effects. Recent experimental evidence in the functional biology of BP control has tested this intuitive assumption. A new paradigm has emerged that BP is biologically determined in modularity by multiple QTLs. Functionally, when a master regulator is taken out, distinct epistatic modules organize biological 'blocks' into a genetic architecture, and serve as basic functional cores from which numerous QTLs act together to physiologically formulate BP. An epistatic module refers to the grouping of QTLs that perform their functions epistatically to one another and influence BP as a group. The modularity mechanism framework indicates that BP as a quantitatively-measured trait is not cumulatively determined and implies that the QTLs in the same epistatic module may participate in the same pathway leading to the BP control, and the QTLs from separate epistatic modules may act in divergent but parallel pathways. This mechanistic conceptualization and subsequent validations synergize with anticipated demands from current human epidemiological studies, since the outcome from them primarily implicates single nucleotide polymorphisms with unknown functions. Eventually, functional understandings of the human results have to be realized by their pathogenic directionality and mechanisms biologically controlling BP.

Pleiotropy facilitates local adaptation to distant optima in common ragweed (Ambrosia artemisiifolia)

Pleiotropy, the control of multiple phenotypes by a single locus, is expected to slow the rate of adaptation by increasing the chance that beneficial alleles also have deleterious effects. However, a prediction arising from classical theory of quantitative trait evolution states that pleiotropic alleles may have a selective advantage when phenotypes are distant from their selective optima. We examine the role of pleiotropy in regulating adaptive differentiation among populations of common ragweed (Ambrosia artemisiifolia); a species that has recently expanded its North American range due to human-mediated habitat change. We employ a phenotype-free approach by using connectivity in gene networks as a proxy for pleiotropy. First, we identify loci bearing footprints of local adaptation, and then use genotype-expression mapping and co-expression networks to infer the connectivity of the genes. Our results indicate that the putatively adaptive loci are highly pleiotropic, as they are more likely than expected to affect the expression of other genes, and they reside in central positions within the gene networks. We propose that the conditionally advantageous alleles at these loci avoid the cost of pleiotropy by having large phenotypic effects that are beneficial when populations are far from their selective optima. We further use evolutionary simulations to show that these patterns are in agreement with a model where populations face novel selective pressures, as expected during a range expansion. Overall, our results suggest that highly connected genes may be targets of positive selection during environmental change, even though they likely experience strong purifying selection in stable selective environments.

Quantitative analyses of adiposity dynamics in zebrafish

Adipose tissues often exhibit subtle, quantitative differences between individuals, leading to a graded series of adiposity phenotypes at the population level. Robust, quantitative analyses are vital for studying these differences. In this Commentary we highlight two articles from our lab that employ sensitive new methods in zebrafish capable of delineating complex and quantitative adiposity phenotypes. In the first article, we utilized in vivo imaging to systematically quantify zebrafish adipose tissues. We identified 34 regionally distinct zebrafish adipose tissues and developed statistical models to predict the size and variance of each adipose tissue over the course of zebrafish growth. We then employed these models to identify effects of strain and diet on adipose tissue growth. In the second article, we employed deep phenotyping to study complex disease-related adiposity traits. Using this methodology, we identified that adipose tissues have unique capacities to re-deposit lipid following food restriction and re-feeding. These distinct re-deposition potentials led to widespread fat distribution changes following re-feeding. We discuss how these novel findings may provide relevance to health conditions such as anorexia nervosa. Together, the strategies described in these two articles can be used as unbiased and quantitative methods to uncover new relationships between genotype, diet and adiposity.

Functional models in genome-wide selection

The development of sequencing technologies has enabled the discovery of markers that are abundantly distributed over the whole genome. Knowledge about the marker locations in reference genomes provides further insights in the search for causal regions and the prediction of genomic values. The present study proposes a Bayesian functional approach for incorporating the marker locations into genomic analysis using stochastic methods to search causal regions and predict genotypic values. For this, three scenarios were analyzed: F2 population with 300 individuals and three different heritability levels (0.2, 0.5, and 0.8), along with 12,150 SNP markers that were distributed through ten linkage groups; F∞ populations with 320 individuals and three different heritability levels (0.2, 0.5, and 0.8), along with 10,020 SNP markers that were distributed through ten linkage groups; and data related to Eucalyptus spp. to measure the model performance in a real LD setting, with 611 individuals whose phenotypes were simulated from QTLs distributed through a panel of 36,812 SNPs with known positions. The performance of the proposed method was compared with those of other genome selection models, namely, RR-BLUP, Bayes B and Bayesian Lasso. The Bayesian functional model presented higher or similar predictive ability when compared with those classical regressions methods in simulated and real scenarios on different LD structures. In general, the Bayesian functional model also achieved higher computational efficiency, using 12 SNPs per MCMC round. The model was efficient in the identification of causal regions and showed high flexibility of analysis, as it is easily adaptable to any genomic selection model.

Widespread Cumulative Influence of Small Effect Size Mutations on Yeast Quantitative Traits

Quantitative traits are influenced by pathways that have traditionally been defined through genes that have a large loss- or gain-of-function effect. However, in theory, a large number of small effect size genes could cumulatively play a substantial role in pathway function. Here, we determine the number, strength, and identity of all non-essential test genes that affect two quantitative galactose-responsive traits in addition to re-analyzing two previously screened quantitative traits. We find that over a quarter of assayed genes have a detectable, quantitative effect on phenotype. Despite their ubiquity, these genes are enriched in core cellular processes in a trait-specific manner. In a simulated population with 50% frequency of all-or-none alleles, we show that small effect size alleles are capable of contributing more to trait variation than alleles in a canonical, large effect size pathway. In total, by demonstrating that the genes effecting quantitative traits can be highly distributed and interconnected, this work challenges the concept of pathways as modular and independent.

Transcription factors involved in the regulatory networks governing the Calvin-Benson-Bassham cycle

Transcription factors (TFs) play crucial roles in the regulation of photosynthesis; elucidating these roles will facilitate our understanding of photosynthesis and thus accelerate its improvement for enhancing crop yield. Promoter analysis of 52 nuclear-encoded Populus tomentosa Carr. genes involved in the Calvin-Benson-Bassham (CBB) cycle revealed 706 motifs and 326 potentially interacting TFs. A backward elimination random forest (BWERF) algorithm reduced the number of TFs to 40, involved in a three-layer gene regulatory network (GRN) including 46 photosynthesis genes (bottom layer), 25 TFs (second layer) and 15 TFs (top layer). Phenotype-genotype association identified 248 single-nucleotide polymorphisms (SNPs) within 72 genes associated with 11 photosynthesis traits. Of the regulatory pairs identified by the BWERF (202 pairs), 77 TF-target combinations harbored SNPs associated with the same trait, supporting similar mechanisms of phenotype modulation. We used expression quantitative trait nucleotide (eQTN) analysis to identify causal SNPs affecting gene expression, identifying 1851 eQTN signals for 50 eGenes (genes whose expressions are regulated by eQTNs). Distribution patterns identified 14 eQTNs from seven TFs associated with eight expression levels of their downstream targets (defined in the GRN), whereas seven TF-target pairs were also identified by phenotype-genotype associations. To further validate the roles of TFs at the metabolic level, we selected 6764 SNPs from 55 genes (identified by GRN-association or GRN-eQTN pairs or both) for metabolic association, identifying variants within 10 TFs affecting metabolic processes underlying the CBB cycle. Our study provides new insights into the photosynthesis pathway in poplar and may facilitate understanding of processes underlying photosynthesis improvement.

A Polygenic Genetic Architecture of Flowering Time in the Worldwide Arabidopsis thaliana Population

Here, we report an empirical study of the polygenic basis underlying the evolution of complex traits. Flowering time variation measured at 10 and 16°C in the 1,001-genomes Arabidopsis thaliana collection of natural accessions were used as a model. The polygenic architecture of flowering time was defined as the 48 loci that were significantly associated with flowering time-at 10 and/or 16°C and/or their difference-in this population. Contributions from alleles at flowering time associated loci to global and local adaptation were explored by evaluating their distribution across genetically and geographically defined subpopulations across the native range of the species. The dynamics in the genetic architecture of flowering time in response to temperature was evaluated by estimating how the effects of these loci on flowering changed with growth temperature. Overall, the genetic basis of flowering time was stable-about 2/3 of the flowering time loci had similar effects at 10°C and 16°C-but many loci were involved in gene by temperature interactions. Globally present alleles, mostly of moderate effect, contributed to the differences in flowering times between the subpopulations via subtle changes in allele frequencies. More extreme local adaptations were, on several occasions, due to regional alleles with relatively large effects, and their linkage disequilibrium-patterns suggest coevolution of functionally connected alleles within local populations. Overall, these findings provide a significant contribution to our understanding about the possible modes of global and local evolution of a complex adaptive trait in A. thaliana.

Assessing the Architecture of Drosophila mojavensis Locomotor Evolution with Bulk Segregant Analysis

Behavior is frequently predicted to be especially important for evolution in novel environments. If these predictions are accurate, there might be particular patterns of genetic architecture associated with recently diverged behaviors. Specifically, it has been predicted that behaviors linked to population divergence should be underpinned by a few genes of relatively large effect, compared to architectures of intrapopulation behavioral variation, which is considered to be highly polygenic. More mapping studies of behavioral variation between recently diverged populations are needed to continue assessing the generality of these predictions. Here, we used a bulk segregant mapping approach to dissect the genetic architecture of a locomotor trait that has evolved between two populations of the cactophilic fly Drosophila mojavensis We created an F8 mapping population of 1,500 individuals from advanced intercross lines and sequenced the 10% of individuals with the highest and lowest levels of locomotor activity. Using three alternative statistical approaches, we found strong evidence for two relatively large-effect QTL that is localized in a region homologous to a region of densely packed behavior loci in Drosophila melanogaster, suggesting that clustering of behavior genes may display relatively deep evolutionary conservation. Broadly, our data are most consistent with a polygenic architecture, though with several loci explaining a high proportion of variation in comparison to similar behavioral traits. We further note the presence of several antagonistic QTL linked to locomotion and discuss these results in light of theories regarding behavioral evolution and the effect size and direction of QTL for diverging traits in general.

Interactions Between Mitochondrial Haplotype and Dietary Macronutrient Ratios Confer Sex-Specific Effects on Longevity in Drosophila melanogaster

Recent studies have demonstrated that modifications to the ratio of dietary macronutrients affect longevity in a diverse range of species. However, the degree to which levels of natural genotypic variation shape these dietary effects on longevity remains unclear. The mitochondria have long been linked to the aging process. The mitochondria possess their own genome, and previous studies have shown that mitochondrial genetic variation affects longevity in insects. Furthermore, the mitochondria are the sites in which dietary nutrients are oxidized to produce adenosine triphosphate, suggesting a capacity for dietary quality to mediate the link between mitochondrial genotype and longevity. Here, we measured longevity of male and female fruit flies, across a panel of genetic strains of Drosophila melanogaster, which vary only in their mitochondrial haplotype, when fed one of the two isocaloric diets that differed in their protein-to-carbohydrate ratio. The mitochondrial haplotype affected the longevity of flies, but the pattern of these effects differed across the two diets in males, but not in females. We discuss the implications of these results in relation to an evolutionary theory linking maternal inheritance of mitochondria to the accumulation of male-harming mitochondrial mutations, and to the theory exploring the evolution of phenotypic plasticity to novel environments.