Detecting directional epistasis and dominance from cross-line analyses in alpine populations of Arabidopsis thaliana.

The contribution of non-additive genetic effects to the genetic architecture of fitness, and to the evolutionary potential of populations, has been a topic of theoretical and empirical interest for a long time. Yet, the empirical study of these effects in natural populations remains scarce, perhaps because measuring dominance and epistasis relies heavily on experimental line crosses. In this study, we explored the contribution of dominance and epistasis in natural alpine populations of Arabidopsis thaliana, for two fitness traits, the dry biomass and the estimated number of siliques, measured in a greenhouse. We found that, on average, crosses between inbred lines of A. thaliana led to mid-parent heterosis for dry biomass, but outbreeding depression for estimated number of siliques. While heterosis for dry biomass was due to dominance, we found that outbreeding depression for estimated number of siliques could be attributed to the breakdown of beneficial epistatic interactions. We simulated and discussed the implication of these results for the adaptive potential of the studied populations, as well as the use of line-cross analyses to detect non-additive genetic effects.

A population genetics theory for piRNA-regulated transposable elements.

Transposable elements (TEs) are self-reproducing selfish DNA sequences that can invade the genome of virtually all living species. Population genetics models have shown that TE copy numbers generally reach a limit, either because the transposition rate decreases with the number of copies (transposition regulation) or because TE copies are deleterious, and thus purged by natural selection. Yet, recent empirical discoveries suggest that TE regulation may mostly rely on piRNAs, which require a specific mutational event (the insertion of a TE copy in a piRNA cluster) to be activated - the so-called TE regulation "trap model". We derived new population genetics models accounting for this trap mechanism, and showed that the resulting equilibria differ substantially from previous expectations based on a transposition-selection equilibrium. We proposed three sub-models, depending on whether or not genomic TE copies and piRNA cluster TE copies are selectively neutral or deleterious, and we provide analytical expressions for maximum and equilibrium copy numbers, as well as cluster frequencies for all of them. In the full neutral model, the equilibrium is achieved when transposition is completely silenced, and this equilibrium does not depend on the transposition rate. When genomic TE copies are deleterious but not cluster TE copies, no long-term equilibrium is possible, and active TEs are eventually eliminated after an active incomplete invasion stage. When all TE copies are deleterious, a transposition-selection equilibrium exists, but the invasion dynamics is not monotonic, and the copy number peaks before decreasing. Mathematical predictions were in good agreement with numerical simulations, except when genetic drift and/or linkage disequilibrium dominates. Overall, the trap-model dynamics appeared to be substantially more stochastic and less repeatable than traditional regulation models.

Using phenotypic plasticity to understand the structure and evolution of the genotype-phenotype map.

Deciphering the genotype-phenotype map necessitates relating variation at the genetic level to variation at the phenotypic level. This endeavour is inherently limited by the availability of standing genetic variation, the rate of spontaneous mutation to novo genetic variants, and possible biases associated with induced mutagenesis. An interesting alternative is to instead rely on the environment as a source of variation. Many phenotypic traits change plastically in response to the environment, and these changes are generally underlain by changes in gene expression. Relating gene expression plasticity to the phenotypic plasticity of more integrated organismal traits thus provides useful information about which genes influence the development and expression of which traits, even in the absence of genetic variation. We here appraise the prospects and limits of such an environment-for-gene substitution for investigating the genotype-phenotype map. We review models of gene regulatory networks, and discuss the different ways in which they can incorporate the environment to mechanistically model phenotypic plasticity and its evolution. We suggest that substantial progress can be made in deciphering this genotype-environment-phenotype map, by connecting theory on gene regulatory network to empirical patterns of gene co-expression, and by more explicitly relating gene expression to the expression and development of phenotypes, both theoretically and empirically.

Optimization of sampling designs for pedigrees and association studies.

In many studies, related individuals are phenotyped in order to infer how their genotype contributes to their phenotype, through the estimation of parameters such as breeding values or locus effects. When it is not possible to phenotype all the individuals, it is important to properly sample the population to improve the precision of the statistical analysis. This article studies how to optimize such sampling designs for pedigrees and association studies. Two sampling methods are developed, stratified sampling and D optimality. It is found that it is important to take account of mutation when sampling pedigrees with many generations: as the size of mutation effects increases, optimized designs sample more individuals in late generations. Optimized designs for association studies tend to improve the joint estimation of breeding values and locus effects, all the more as sample size is low and the genetic architecture of the trait is simple. When the trait is determined by few loci, they are reminiscent of classical experimental designs for regression models and tend to select homozygous individuals. When the trait is determined by many loci, locus effects may be difficult to estimate, even if an optimized design is used.

Quantitative trait loci involved in the reproductive success of a parasitoid wasp.

Dissecting the genetic basis of intraspecific variations in life history traits is essential to understand their evolution, notably for potential biocontrol agents. Such variations are observed in the endoparasitoid Cotesia typhae (Hymenoptera: Braconidae), specialized on the pest Sesamia nonagrioides (Lepidoptera: Noctuidae). Previously, we identified two strains of C. typhae that differed significantly for life history traits on an allopatric host population. To investigate the genetic basis underlying these phenotypic differences, we used a quantitative trait locus (QTL) approach based on restriction site-associated DNA markers. The characteristic of C. typhae reproduction allowed us generating sisters sharing almost the same genetic content, named clonal sibship. Crosses between individuals from the two strains were performed to generate F2 and F8 recombinant CSS. The genotypes of 181 clonal sibships were determined as well as the phenotypes of the corresponding 4,000 females. Informative markers were then used to build a high-quality genetic map. These 465 markers spanned a total length of 1,300 cM and were organized in 10 linkage groups which corresponded to the number of C. typhae chromosomes. Three QTLs were detected for parasitism success and two for offspring number, while none were identified for sex ratio. The QTLs explained, respectively, 27.7% and 24.5% of the phenotypic variation observed. The gene content of the genomic intervals was investigated based on the genome of C. congregata and revealed 67 interesting candidates, as potentially involved in the studied traits, including components of the venom and of the symbiotic virus (bracovirus) shown to be necessary for parasitism success in related wasps.

Symbiont-Driven Male Mating Success in the Neotropical Drosophila paulistorum Superspecies.

Microbial symbionts are ubiquitous associates of living organisms but their role in mediating reproductive isolation (RI) remains controversial. We addressed this knowledge gap by employing the Drosophila paulistorum-Wolbachia model system. Semispecies in the D. paulistorum species complex exhibit strong RI between each other and knockdown of obligate mutualistic Wolbachia bacteria in female D. paulistorum flies triggers loss of assortative mating behavior against males carrying incompatible Wolbachia strains. Here we set out to determine whether de novo RI can be introduced by Wolbachia-knockdown in D. paulistorum males. We show that Wolbachia-knockdown D. paulistorum males (i) are rejected as mates by wild type females, (ii) express altered sexual pheromone profiles, and (iii) are devoid of the endosymbiont in pheromone producing cells. Our findings suggest that changes in Wolbachia titer and tissue tropism can induce de novo premating isolation by directly or indirectly modulating sexual behavior of their native D. paulistorum hosts.

Experimental evolution reveals hyperparasitic interactions among transposable elements.

Transposable elements (TEs) are repeated DNA sequences that can constitute a substantial part of genomes. Studying TEs' activity, interactions, and accumulation dynamics is thus of major interest to understand genome evolution. Here, we describe the transposition dynamics of cut-and-paste mariner elements during experimental (short- and longer-term) evolution in Drosophila melanogaster Flies with autonomous and nonautonomous mariner copies were introduced in populations containing no active mariner, and TE accumulation was tracked by quantitative PCR for up to 100 generations. Our results demonstrate that (i) active mariner elements are highly invasive and characterized by an elevated transposition rate, confirming their capacity to spread in populations, as predicted by the "selfish-DNA" mechanism; (ii) nonautonomous copies act as parasites of autonomous mariner elements by hijacking the transposition machinery produced by active mariner, which can be considered as a case of hyperparasitism; (iii) this behavior resulted in a failure of active copies to amplify which systematically drove the whole family to extinction in less than 100 generations. This study nicely illustrates how the presence of transposition-competitive variants can deeply impair TE dynamics and gives clues to the extraordinary diversity of TE evolutionary histories observed in genomes.

Modelling the influence of parental effects on gene-network evolution.

Understanding the importance of nongenetic heredity in the evolutionary process is a major topic in modern evolutionary biology. We modified a classical gene-network model by allowing parental transmission of gene expression and studied its evolutionary properties through individual-based simulations. We identified ontogenetic time (i.e. the time gene networks have to stabilize before being submitted to natural selection) as a crucial factor in determining the evolutionary impact of this phenotypic inheritance. Indeed, fast-developing organisms display enhanced adaptation and greater robustness to mutations when evolving in presence of nongenetic inheritance (NGI). In contrast, in our model, long development reduces the influence of the inherited state of the gene network. NGI thus had a negligible effect on the evolution of gene networks when the speed at which transcription levels reach equilibrium is not constrained. Nevertheless, simulations show that intergenerational transmission of the gene-network state negatively affects the evolution of robustness to environmental disturbances for either fast- or slow-developing organisms. Therefore, these results suggest that the evolutionary consequences of NGI might not be sought only in the way species respond to selection, but also on the evolution of emergent properties (such as environmental and genetic canalization) in complex genetic architectures.

Fatty acid synthase cooperates with glyoxalase 1 to protect against sugar toxicity.

Fatty acid (FA) metabolism is deregulated in several human diseases including metabolic syndrome, type 2 diabetes and cancers. Therefore, FA-metabolic enzymes are potential targets for drug therapy, although the consequence of these treatments must be precisely evaluated at the organismal and cellular levels. In healthy organism, synthesis of triacylglycerols (TAGs)-composed of three FA units esterified to a glycerol backbone-is increased in response to dietary sugar. Saturation in the storage and synthesis capacity of TAGs is associated with type 2 diabetes progression. Sugar toxicity likely depends on advanced-glycation-end-products (AGEs) that form through covalent bounding between amine groups and carbonyl groups of sugar or their derivatives α-oxoaldehydes. Methylglyoxal (MG) is a highly reactive α-oxoaldehyde that is derived from glycolysis through a non-enzymatic reaction. Glyoxalase 1 (Glo1) works to neutralize MG, reducing its deleterious effects. Here, we have used the power of Drosophila genetics to generate Fatty acid synthase (FASN) mutants, allowing us to investigate the consequence of this deficiency upon sugar-supplemented diets. We found that FASN mutants are lethal but can be rescued by an appropriate lipid diet. Rescued animals do not exhibit insulin resistance, are dramatically sensitive to dietary sugar and accumulate AGEs. We show that FASN and Glo1 cooperate at systemic and cell-autonomous levels to protect against sugar toxicity. We observed that the size of FASN mutant cells decreases as dietary sucrose increases. Genetic interactions at the cell-autonomous level, where glycolytic enzymes or Glo1 were manipulated in FASN mutant cells, revealed that this sugar-dependent size reduction is a direct consequence of MG-derived-AGE accumulation. In summary, our findings indicate that FASN is dispensable for cell growth if extracellular lipids are available. In contrast, FA-synthesis appears to be required to limit a cell-autonomous accumulation of MG-derived-AGEs, supporting the notion that MG is the most deleterious α-oxoaldehyde at the intracellular level.

VHICA, a New Method to Discriminate between Vertical and Horizontal Transposon Transfer: Application to the Mariner Family within Drosophila.

Transposable elements (TEs) are genomic repeated sequences that display complex evolutionary patterns. They are usually inherited vertically, but can occasionally be transmitted between sexually independent species, through so-called horizontal transposon transfers (HTTs). Recurrent HTTs are supposed to be essential in life cycle of TEs, which are otherwise destined for eventual decay. HTTs also impact the host genome evolution. However, the extent of HTTs in eukaryotes is largely unknown, due to the lack of efficient, statistically supported methods that can be applied to multiple species sequence data sets. Here, we developed a new automated method available as a R package "vhica" that discriminates whether a given TE family was vertically or horizontally transferred, and potentially infers donor and receptor species. The method is well suited for TE sequences extracted from complete genomes, and applicable to multiple TEs and species at the same time. We first validated our method using Drosophila TE families with well-known evolutionary histories, displaying both HTTs and vertical transmission. We then tested 26 different lineages of mariner elements recently characterized in 20 Drosophila genomes, and found HTTs in 24 of them. Furthermore, several independent HTT events could often be detected within the same mariner lineage. The VHICA (Vertical and Horizontal Inheritance Consistence Analysis) method thus appears as a valuable tool to analyze the evolutionary history of TEs across a large range of species.

Why and how genetic canalization evolves in gene regulatory networks.

BACKGROUND: Genetic canalization reflects the capacity of an organism's phenotype to remain unchanged in spite of mutations. As selection on genetic canalization is weak and indirect, whether or not genetic canalization can reasonably evolve in complex genetic architectures is still an open question. In this paper, we use a quantitative model of gene regulatory network to describe the conditions in which substantial canalization is expected to emerge in a stable environment. RESULTS: Through an individual-based simulation framework, we confirmed that most parameters associated with the network topology (complexity and size of the network) have less influence than mutational parameters (rate and size of mutations) on the evolution of genetic canalization. We also established that selecting for extreme phenotypic optima (nil or full gene expression) leads to much higher canalization levels than selecting for intermediate expression levels. Overall, constrained networks evolve less canalization than networks in which some genes could evolve freely (i.e. without direct stabilizing selection pressure on gene expression). CONCLUSIONS: Taken together, these results lead us to propose a two-fold mechanism involved in the evolution of genetic canalization in gene regulatory networks: the shrinkage of mutational target (useless genes are virtually removed from the network) and redundancy in gene regulation (so that some regulatory factors can be lost without affecting gene expression).

Epistasis-Induced Evolutionary Plateaus in Selection Responses.

Understanding and predicting evolution is a central challenge in both population and quantitative genetics. The amount of genetic variance for quantitative traits available in a population conditions the particular way in which this population will (or will not) evolve under natural or artificial selection. Here, we explore the potential of gene-gene interactions (epistasis) to induce evolutionary plateaus at which evolutionary change virtually collapses for a number of generations, followed by the release of previously cryptic genetic variation. First, we demonstrate theoretically that a wide range of epistatic interactions has the potential to generate temporary decelerations in the course of response to selection. Second, we perform simulations to show that such microevolutionary plateaus may occur in selection responses under empirically based assumptions. Finally, we show that such events can be traced in artificial selection experiments, thus providing further empirical evidence for this phenomenon.

Evolutionary time-series analysis reveals the signature of frequency-dependent selection on a female mating polymorphism.

A major challenge in evolutionary biology is understanding how stochastic and deterministic factors interact and influence macroevolutionary dynamics in natural populations. One classical approach is to record frequency changes of heritable and visible genetic polymorphisms over multiple generations. Here, we combined this approach with a maximum likelihood-based population-genetic model with the aim of understanding and quantifying the evolutionary processes operating on a female mating polymorphism in the blue-tailed damselfly Ischnura elegans. Previous studies on this color-polymorphic species have suggested that males form a search image for females, which leads to excessive mating harassment of common female morphs. We analyzed a large temporally and spatially replicated data set of between-generation morph frequency changes in I. elegans. Morph frequencies were more stable than expected from genetic drift alone, suggesting the presence of selection toward a stable equilibrium that prevents local loss or fixation of morphs. This can be interpreted as the signature of negative frequency-dependent selection maintaining the phenotypic stasis and genetic diversity in these populations. Our novel analytical approach allows the estimation of the strength of frequency-dependent selection from the morph frequency fluctuations around their inferred long-term equilibria. This approach can be extended and applied to other polymorphic organisms for which time-series data across multiple generations are available.

Male manipulation impinges on social-dependent tumor suppression in Drosophila melanogaster females.

Physiological status can influence social behavior, which in turn can affect physiology and health. Previously, we reported that tumor growth in Drosophila virgin females depends on the social context, but did not investigate the underlying physiological mechanisms. Here, we sought to characterize the signal perceived between tumorous flies, ultimately discovering that the tumor suppressive effect varies depending on reproductive status. Firstly, we show that the tumor suppressive effect is neither dependent on remnant pheromone-like products nor on the microbiota. Transcriptome analysis of the heads of these tumorous flies reveals social-dependent gene-expression changes related to nervous-system activity, suggesting that a cognitive-like relay might mediate the tumor suppressive effect. The transcriptome also reveals changes in the expression of genes related to mating behavior. Surprisingly, we observed that this social-dependent tumor-suppressive effect is lost in fertilized females. After mating, Drosophila females change their behavior-favoring offspring survival-in response to peptides transferred via the male ejaculate, a phenomenon called "male manipulation". Remarkably, the social-dependent tumor suppressive effect is restored in females mated by sex-peptide deficient males. Since male manipulation has likely been selected to favor male gene transmission, our findings indicate that this evolutionary trait impedes social-dependent tumor growth slowdown.

Genomic parasites or symbionts? Modeling the effects of environmental pressure on transposition activity in asexual populations.

Transposable elements are DNA segments capable of persisting in host genomes by self-replication in spite of deleterious mutagenic effects. The theoretical dynamics of these elements within genomes has been studied extensively, and population genetic models predict that they can invade and maintain as a result of both intra-genomic and inter-individual selection in sexual species. In asexuals, the success of selfish DNA is more difficult to explain. However, most theoretical work assumes constant environment. Here, we analyze the impact of environmental change on the dynamics of transposition activity when horizontal DNA exchange is absent, based on a stochastic computational model of transposable element proliferation. We argue that repeated changes in the phenotypic optimum in a multidimensional fitness landscape may induce explosive bursts of transposition activity associated with faster adaptation. However, long-term maintenance of transposition activity is unlikely. This could contribute to the significant variation in the transposable element copy number among closely related species.

Correlated stabilizing selection shapes the topology of gene regulatory networks.

The evolution of gene expression is constrained by the topology of gene regulatory networks, as co-expressed genes are likely to have their expressions affected together by mutations. Conversely, co-expression can also be an advantage when genes are under joint selection. Here, we assessed theoretically whether correlated selection (selection for a combination of traits) was able to affect the pattern of correlated gene expressions and the underlying gene regulatory networks. We ran individual-based simulations, applying a stabilizing correlated fitness function to three genetic architectures: a quantitative genetics (multilinear) model featuring epistasis and pleiotropy, a quantitative genetics model where each genes has an independent mutational structure, and a gene regulatory network model, mimicking the mechanisms of gene expression regulation. Simulations showed that correlated mutational effects evolved in the three genetic architectures as a response to correlated selection, but the response in gene networks was specific. The intensity of gene co-expression was mostly explained by the regulatory distance between genes (largest correlations being associated to genes directly interacting with each other), and the sign of co-expression was associated with the nature of the regulation (transcription activation or inhibition). These results concur to the idea that gene network topologies could partly reflect past selection patterns on gene expression.

Modelling of genetic interactions improves prediction of hybrid patterns--a case study in domestic fowl.

A major challenge in complex trait genetics is to unravel how multiple loci and environmental factors together cause phenotypic diversity. Both first (F(1)) and second (F(2)) generation hybrids often display phenotypes that deviate from what is expected under intermediate inheritance. We have here studied two chicken F(2) populations generated by crossing divergent chicken lines to assess how epistatic loci, identified in earlier quantitative trait locus (QTL) studies, contribute to hybrid deviations from the mid-parent phenotype. Empirical evidence suggests that the average phenotypes of the intercross birds tend to be lower than the midpoint between the parental means in both crosses. Our results confirm that epistatic interactions, despite a relatively small contribution to the phenotypic variance, play an important role in the deviation of hybrid phenotypes from the mid-parent values (i.e. multi-locus hybrid genotypes lead to lower rather than higher body weights). To a lesser extent, dominance also appears to contribute to the mid-parent deviation, at least in one of the crosses. This observation coincides with the hypothesis that hybridization tends to break up co-adapted gene complexes, i.e. generate Bateson-Dobzhansky-Muller incompatibilities.

Gene network simulations provide testable predictions for the molecular domestication syndrome.

The domestication of plant species leads to repeatable morphological evolution, often referred to as the phenotypic domestication syndrome. Domestication is also associated with important genomic changes, such as the loss of genetic diversity compared with adequately large wild populations, and modifications of gene expression patterns. Here, we explored theoretically the effect of a domestication-like scenario on the evolution of gene regulatory networks. We ran population genetics simulations in which individuals were featured by their genotype (an interaction matrix encoding a gene regulatory network) and their gene expressions, representing the phenotypic level. Our domestication scenario included a population bottleneck and a selection switch mimicking human-mediated directional and canalizing selection, i.e., change in the optimal gene expression level and selection toward more stable expression across environments. We showed that domestication profoundly alters genetic architectures. Based on four examples of plant domestication scenarios, our simulations predict (1) a drop in neutral allelic diversity; (2) a change in gene expression variance that depends upon the domestication scenario; (3) transient maladaptive plasticity; (4) a deep rewiring of the gene regulatory networks, with a trend toward gain of regulatory interactions; and (5) a global increase in the genetic correlations among gene expressions, with a loss of modularity in the resulting coexpression patterns and in the underlying networks. We provide empirically testable predictions on the differences of genetic architectures between wild and domesticated forms. The characterization of such systematic evolutionary changes in the genetic architecture of traits contributes to define a molecular domestication syndrome.

Strong and consistent natural selection associated with armour reduction in sticklebacks.

Measuring the strength of natural selection is tremendously important in evolutionary biology, but remains a challenging task. In this work, we analyse the characteristics of selection for a morphological change (lateral-plate reduction) in the threespine stickleback Gasterosteus aculeatus. Adaptation to freshwater, leading with the reduction or loss of the bony lateral armour, has occurred in parallel on numerous occasions in this species. Completely-plated and low-plated sticklebacks were introduced into a pond, and the phenotypic changes were tracked for 20 years. Fish from the last generation were genotyped for the Ectodysplasin-A (Eda) locus, the major gene involved in armour development. We found a strong fitness advantage for the freshwater-type fish (on average, 20% fitness advantage for the freshwater morph, and 92% for the freshwater genotype). The trend is best explained by assuming that this fitness advantage is maximum at the beginning of the invasion and decreases with time. Such fitness differences provide a quantifiable example of rapid selection-driven phenotypic evolution associated with environmental change in a natural population.