Cospeciation vs host-shift speciation: methods for testing, evidence from natural associations and relation to coevolution.

Hosts and their symbionts are involved in intimate physiological and ecological interactions. The impact of these interactions on the evolution of each partner depends on the time-scale considered. Short-term dynamics - 'coevolution' in the narrow sense - has been reviewed elsewhere. We focus here on the long-term evolutionary dynamics of cospeciation and speciation following host shifts. Whether hosts and their symbionts speciate in parallel, by cospeciation, or through host shifts, is a key issue in host-symbiont evolution. In this review, we first outline approaches to compare divergence between pairwise associated groups of species, their advantages and pitfalls. We then consider recent insights into the long-term evolution of host-parasite and host-mutualist associations by critically reviewing the literature. We show that convincing cases of cospeciation are rare (7%) and that cophylogenetic methods overestimate the occurrence of such events. Finally, we examine the relationships between short-term coevolutionary dynamics and long-term patterns of diversification in host-symbiont associations. We review theoretical and experimental studies showing that short-term dynamics can foster parasite specialization, but that these events can occur following host shifts and do not necessarily involve cospeciation. Overall, there is now substantial evidence to suggest that coevolutionary dynamics of hosts and parasites do not favor long-term cospeciation.

The genotype-phenotype relationship and evolutionary genetics in the light of the Metabolic Control Analysis.

Metabolic control analysis has long been used as a systemic model of the genotype-phenotype (GP) relationship. By considering kinetic parameters and enzyme concentrations as reflecting the genotype level and metabolic fluxes or pools as phenotypes related to fitness, MCA has given a biological basis to the relationship between these two levels. The non-linear and concave relationship between enzymes and fluxes can account for common genetic effects that reductionist approaches have been powerless to explain, such as the dominance of active alleles over less active alleles, the various types of epistasis and heterosis, and reveals the structural links between these genetic effects. The summation property of the flux control coefficients accounts for the L-shaped distribution of Quantitative Trait Locus (QTL) effects, irrespective of other possible causes. Metabolic models of response to selection results in evolutionary scenarios that are markedly different from those derived from the classical infinitesimal model of quantitative genetics. In particular, evolution towards selective neutrality appears to be a consequence of the diminishing return of the flux-enzyme relationship. In this paper, we survey the historical and recent achievements of MCA in genetics, quantitative genetics and evolution, focusing on epistasis and the evolution of flux in relation to enzyme concentrations.

The decoupling between genetic structure and metabolic phenotypes in Escherichia coli leads to continuous phenotypic diversity.

To assess the extent of intra-species diversity and the links between phylogeny, lifestyle (habitat and pathogenicity) and phenotype, we assayed the growth yield on 95 carbon sources of 168 Escherichia strains. We also correlated the growth capacities of 14 E. coli strains with the presence/absence of enzyme-coding genes. Globally, we found that the genetic distance, based on multilocus sequence typing data, was a weak indicator of the metabolic phenotypic distance. Besides, lifestyle and phylogroup had almost no impact on the growth yield of non-Shigella E. coli strains. In these strains, the presence/absence of the metabolic pathways, which was linked to the phylogeny, explained most of the growth capacities. However, few discrepancies blurred the link between metabolic phenotypic distance and metabolic pathway distance. This study shows that a prokaryotic species structured into well-defined genetic and lifestyle groups can yet exhibit continuous phenotypic diversity, possibly caused by gene regulatory effects.

Sex in Penicillium: combined phylogenetic and experimental approaches.

We studied the mode of reproduction and its evolution in the fungal subgenus Penicillium Biverticillium using phylogenetic and experimental approaches. We sequenced mating type (MAT) genes and nuclear DNA fragments in sexual and putatively asexual species. Examination of the concordance between individual trees supported the recognition of the morphological species. MAT genes were detected in two putatively asexual species and were found to evolve mostly under purifying selection, although high substitution rates were detected at some sites in some clades. The first steps of sexual reproduction could be induced under controlled conditions in one of the two species, although no mature cleistothecia were produced. Altogether, these findings suggest that the asexual Penicillium species may have lost sex only very recently and/or that the MAT genes are involved in other functions. An ancestral state reconstruction analysis indicated several events of putative sex loss in the genus. Alternatively, it is possible that the supposedly asexual Penicillium species may have retained a cryptic sexual stage.

Phylogenetic determinants of potential host shifts in fungal pathogens.

Understanding what determines the host range of pathogens and the potential for host shifts is of critical importance to controlling their introductions into new environments. The phylogeny of the hosts has been shown to be important: pathogens are more likely to be infectious on hosts closely related to their host-of-origin because of the similar host environments that is shared by descent. The importance of pathogen phylogenies for predicting host range has never been investigated, although a pathogen should also be able to exploit a new host that its close relative can infect. We performed cross-inoculations using a plant-fungal association and showed that both host and pathogen phylogenies were significant predictors of host range, with at least partly independent effects. Furthermore, we showed that some pathogens were better at infecting novel hosts. Our results should have implications in the context of biological invasions and emergences of new diseases due to globalization.

Evidence for autotetraploidy associated with reproductive isolation in Saccharomyces cerevisiae: towards a new domesticated species.

Partial or whole-genome duplications have played a major role in the evolution of new species. We have investigated the variation of ploidy level in a panel of domesticated strains of Saccharomyces cerevisiae coming from different geographical origins. Segregation studies and crosses with tester strains of different ploidy levels showed that part of the strains were well-balanced autotetraploids displaying tetrasomic inheritance. The presence of up to four different alleles for various loci is consistent with a polyploidization mechanism relying on the fusion of two nonreduced meiospores coming from two S. cerevisiae strains. Autotetraploidy was also in accordance with karyotype and flow cytometry analyses. Interestingly, most bakery strains were tetraploids, suggesting a link between ploidy level and human use. The null or drastically reduced fertility of the hybrids between tetraploid and diploid strains indicated that domesticated S. cerevisiae strains are composed of two groups isolated by post-zygotic reproductive barriers.

Hybrid sterility and inviability in the parasitic fungal species complex Microbotryum.

Microbotryum violaceum, the anther-smut fungus, forms a complex of sibling species which specialize on different plants. Previous studies have shown the presence of partial ecological isolation and F1 inviability, but did not detect assortative mating apart from a high selfing rate. We investigated other post-mating barriers and show that F1 hybrid sterility, the inability of gametes to mate, increased gradually with the increasing genetic distance between the parents. F2 hybrids showed a reduced ability to infect the plants that was also correlated with the genetic distance. The host on which the F2 hybrids were passaged caused a selection for alleles derived from the pathogen species originally isolated from that host, but this effect was not detectable for the most closely related species. The post-mating barriers thus remain weak among the closest species pairs, suggesting that premating barriers are sufficient to initiate divergence in this system.

Patterns of recombination and MLH1 foci density along mouse chromosomes: modeling effects of interference and obligate chiasma.

Crossover interference in meiosis is often modeled via stationary renewal processes. Here we consider a new model to incorporate the known biological feature of "obligate chiasma" whereby in most organisms each bivalent almost always has at least one crossover. The initial crossover is modeled as uniformly distributed along the chromosome, and starting from its position, subsequent crossovers are placed with forward and backward stationary renewal processes using a chi-square distribution of intercrossover distances. We used our model as well as the standard chi-square model to simulate the patterns of crossover densities along bivalents or chromatids for those having zero, one, two, or three or more crossovers; indeed, such patterns depend on the number of crossovers. With both models, simulated patterns compare very well to those found experimentally in mice, both for MLH1 foci on bivalents and for crossovers on genetic maps. However, our model provides a better fit to experimental data as compared to the standard chi-square model, particularly regarding the distribution of numbers of crossovers per chromosome. Finally, our model predicts an enhancement of the recombination rate near the extremities, which, however, explains only a part of the pattern observed in mouse.

Assessment of pleiotropic effects of a gene substitution in pea by two-dimensional polyacrylamide gel electrophoresis.

We examined, by two-dimensional polyacrylamide gel electrophoresis (2D-PAGE), near-isogenic lines of the r-gene in pea (Pisum sativum) which determines round (RR) vs. wrinkled (rr) seed. The study was undertaken to assess the number of protein changes resulting from a single gene substitution as a means of quantifying pleiotropic effects. A total of 636 to 770 resolvable polypeptides were identical in all respects between RR and rr for roots, shoots, leaflets, stipules, young ovaries, and young embryos. A single difference between the lines became evident about 21-23 days after anthesis in the embryos. Mature seeds of the two lines showed 62 spot differences in addition to differences in four clusters of spots, representing about 10% of the total number of spots visible on the gels. The protein differences are presumably involved in the many known physiological differences of the two seed types. 2-D PAGE analyses of near-isogenic lines are likely to be valuable in a number of quantitative and developmental genetic contexts.

Fluxes and metabolic pools as model traits for quantitative genetics. I. The L-shaped distribution of gene effects.

The fluxes through metabolic pathways can be considered as model quantitative traits, whose QTL are the polymorphic loci controlling the activity or quantity of the enzymes. Relying on metabolic control theory, we investigated the relationships between the variations of enzyme activity along metabolic pathways and the variations of the flux in a population with biallelic QTL. Two kinds of variations were taken into account, the variation of the average enzyme activity across the loci, and the variation of the activity of each enzyme of the pathway among the individuals of the population. We proposed analytical approximations for the flux mean and variance in the population as well as for the additive and dominance variances of the individual QTL. Monte Carlo simulations based on these approximations showed that an L-shaped distribution of the contributions of individual QTL to the flux variance (R(2)) is consistently expected in an F(2) progeny. This result could partly account for the classically observed L-shaped distribution of QTL effects for quantitative traits. The high correlation we found between R(2) value and flux control coefficients variance suggests that such a distribution is an intrinsic property of metabolic pathways due to the summation property of control coefficients.

Genetic and nongenetic bases for the L-shaped distribution of quantitative trait loci effects.

The L-shaped distribution of estimated QTL effects (R(2)) has long been reported. We recently showed that a metabolic mechanism could account for this phenomenon. But other nonexclusive genetic or nongenetic causes may contribute to generate such a distribution. Using analysis and simulations of an additive genetic model, we show that linkage disequilibrium between QTL, low heritability, and small population size may also be involved, regardless of the gene effect distribution. In addition, a comparison of the additive and metabolic genetic models revealed that estimates of the QTL effects for traits proportional to metabolic flux are far less robust than for additive traits. However, in both models the highest R(2)'s repeatedly correspond to the same set of QTL.

Quantitative trait loci underlying gene product variation: a novel perspective for analyzing regulation of genome expression.

A methodology to dissect the genetic architecture of quantitative variation of numerous gene products simultaneously is proposed. For each individual of a segregating progeny, proteins extracted from a given organ are separated using two-dimensional electrophoresis, and their amounts are estimated with a computer-assisted system for spot quantification. Provided a complete genetic map is available, statistical procedures allow determination of the number, effects and chromosomal locations of factors controlling the amounts of individual proteins. This approach was applied to anonymous proteins of etiolated coleoptiles of maize, in an F2 progeny between two distant lines. The genetic map included both restriction fragment length polymorphism and protein markers. Minimum estimates of one to five unlinked regulatory factors were found for 42 of the 72 proteins analyzed, with a large diversity of effects. Dominance and epistasis interactions were involved in the control of 38% and 14% of the 72 proteins, respectively. Such a methodology might help understanding the architecture of regulatory networks and the possible adaptive or phenotypic significance of the polymorphism of the genes involved.

Protein changes in response to progressive water deficit in maize . Quantitative variation and polypeptide identification

Three-week-old plants of two unrelated lines of maize (Zea mays L.) and their hybrid were submitted to progressive water stress for 10 d. Changes induced in leaf proteins were studied by two-dimensional electrophoresis and quantitatively analyzed using image analysis. Seventy-eight proteins out of a total of 413 showed a significant quantitative variation (increase or decrease), with 38 of them exhibiting a different expression in the two genotypes. Eleven proteins that increased by a factor of 1.3 to 5 in stressed plants and 8 proteins detected only in stressed plants were selected for internal amino acid microsequencing, and by similarity search 16 were found to be closely related to previously reported proteins. In addition to proteins already known to be involved in the response to water stress (e.g. RAB17 [Responsive to ABA]), several enzymes involved in basic metabolic cellular pathways such as glycolysis and the Krebs cycle (e.g. enolase and triose phosphate isomerase) were identified, as well as several others, including caffeate O-methyltransferase, the induction of which could be related to lignification.

PERMANENT GENETIC RESOURCES: Isolation of 60 polymorphic microsatellite loci in EST libraries of four sibling species of the phytopathogenic fungal complex Microbotryum.

We report the development of 60 microsatellite markers on four species of the fungal complex Microbotryum, causing anther smut of the Caryophyllaceae. Microsatellites were found in four expressed sequence tag (EST) libraries, built from isolates of M. lychnis-dioicae, M. violaceum sensus stricto, M. lagerheimii and M. dianthorum, collected, respectively, from the plants Silene latifolia, S. nutans, S. vulgaris and Dianthus carthusianorum. Intrapopulation polymorphism was investigated using 24 isolates, and cross-amplification was explored using 23 isolates belonging to at least 10 different Microbotryum species. This study provides numerous microsatellite markers for population genetics and mapping studies.

When can host shifts produce congruent host and parasite phylogenies? A simulation approach.

Congruence between host and parasite phylogenies is often taken as evidence for cospeciation. However, 'pseudocospeciation', resulting from host-switches followed by parasite speciation, may also generate congruent trees. To investigate this process and the conditions favouring its appearance, we here simulated the adaptive radiation of a parasite onto a new range of hosts. A very high congruence between the host tree and the resulting parasite trees was obtained when parasites switched between closely related hosts. Setting a shorter time lag for speciation after switches between distantly related hosts further increased the degree of congruence. The shape of the host tree, however, had a strong impact, as no congruence could be obtained when starting with highly unbalanced host trees. The strong congruences obtained were erroneously interpreted as the result of cospeciations by commonly used phylogenetic software packages despite the fact that all speciations resulted from host-switches in our model. These results highlight the importance of estimating the age of nodes in host and parasite phylogenies when testing for cospeciation and also demonstrate that the results obtained with software packages simulating evolutionary events must be interpreted with caution.

Maize Sh2 gene is constrained by natural selection but escaped domestication.

In Zea mays L., we studied the molecular evolution of Shrunken2 (Sh2), a gene that encodes the large subunits of a major enzyme in endosperm starch biosynthesis, ADP-glucose pyrophosphorylase. We compared 4669 bp of the Sh2 coding region on 50 accessions of maize and teosinte. Very few nucleotide polymorphisms were found when compared with other genes in Z. mays, revealing an effect of purifying selection in the whole species that predates domestication. Additionally, the comparison of Sh2 sequences in all Z. mays subspecies and outgroups Z. diploperennis and Tripsacum dactyloides suggests the occurrence of an ancient selective sweep in the Sh2 3' region. The amount and nature of nucleotide diversity are similar in both maize and teosinte, confirming previous results that suggested that Sh2 has not been involved in maize domestication. The very low level of nucleotide diversity as well as the highly conserved protein sequence suggest that natural selection retained effective Sh2 allele(s) long before agriculture started, making human selection inefficient on this gene.

Silver staining of proteins: standardized procedure for two-dimensional gels bound to polyester sheets.

A system for the simultaneous silver staining of numerous two-dimensional gels bound to polyester sheets is proposed. Under controlled conditions, a very good reproducibility of staining can be obtained, even between different experiments. Procedures to avoid artifacts due to polyester sheets (background, silver mirror, and vertical streaking) are described.

Biochemical and genetic properties of oligomeric structures: a general approach.

The oligomers constituted by association of different subunits can exist under multiple forms. In the case of the genetically variable proteins, such a multiplicity leads to numerous questions (i) on the enumerations: what is the number of active forms when a given subunit can make the oligomer inactive, or when the subunits are encoded by s alleles; (ii) on the subunit effects on biochemical properties: how to estimate these effects, are they equal, are there interactions between subunits, etc. Theoretical methods for the study of such oligomeric structures are developed, which mainly rely on linear model techniques. Peculiar properties examined are Vmax and Km, but also the quantities of the various oligomers, which depend on their association law. This approach is extended to the oligomers composed of different sets of subunits, as are for example some enzymes. These aspects are discussed from numerous bibliographic examples, with special reference to molecular interactions (protein complementation or molecular heterosis). Otherwise the genetic application of this theoretical approach is presented: it is possible to consider a genotype as an oligomer of alleles, and thus to study their effects and their interactions, in the one-locus case as well as in the several-loci case. The relevance of this generalization is discussed in connection with two other concepts, the "sequence space" used in molecular evolution and the regression of the genotypic values on the number of alleles used in quantitative genetics.

Proteomics: a link between genomics, genetics and physiology.

Thanks to spectacular advances in the techniques for identifying proteins separated by two-dimensional electrophoresis and in methods for large-scale analysis of proteome variations, proteomics is becoming an essential methodology in various fields of plant biology. In the study of pleiotropic effects of mutants and in the analysis of responses to hormones and to environmental changes, the identification of involved metabolic pathways can be deduced from the function of affected proteins. In molecular quantitative genetics, proteomics can be used to map translated genes and loci controlling their expression, which can be used to identify proteins accounting for the variation of complex phenotypic traits. Linking gene expression to cell metabolism on the one hand and to genetic maps on the other, proteomics is a central tool for functional genomics.

Differential Two-Dimensional Protein Patterns as Related to Tissue Specificity and Water Conditions in Brassica napus var oleifera Root System.

Differential two-dimensional protein patterns as related to tissue specificity and water conditions were investigated within Brassica napus var oleifera root system. The different parts of the root system (tap root, lateral roots, and drought-induced short roots) were analyzed under various moisture regimes (regular watering at field capacity, progressive drought stress, and rewatering). Tissue specificity was evident from 25 differences in protein patterns (qualitative and quantitative) between well-watered lateral and tap roots. Twice as many polypeptides (52) were drought-affected and the response to the water stress was shown to be similar in both root types. In addition, more than half of the polypeptides detected as organ-specific were affected by drought. Based upon the trend of variation observed under drought and rehydration, three categories of polypeptides could be defined that might be differently involved in drought susceptibility or tolerance. A highly differentiated protein pattern characterized the drought-induced short roots. This pattern appeared as far from the watered as from the water-stressed normal roots. In particular, 13 unique polypeptides were detected which could be relevant to their adaptive morphogenesis and/or their specific drought tolerance induction. Upon rehydration, their polypeptide pattern and their specific morphology returned to a normal well-watered lateral root type.