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.

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

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.

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.

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.

QTL analysis of proteome and transcriptome variations for dissecting the genetic architecture of complex traits in maize.

In this review, we present some studies on genetic analysis of proteome and transcriptome variations, which exemplify new strategies for a better understanding of the molecular and genetic bases of complex traits. A large genetic variability was revealed at the proteome expression level, which raised the possibility to predict phenotypical performance on the basis of gene product variability. This approach yielded limited results, but could be re-newed by extensive identification of proteins now allowed by mass spectrometry. The dissection of the genetic basis of the variation of individual protein amounts proves very powerful to select 'candidate' proteins, physiologically relevant for a given phenotypical trait, as shown by a study on the effect of water stress in maize. In order to investigate factors of grain quality in maize, we selected a regulatory locus known to control the expression of several storage protein genes, Opaque-2, and investigated the relationships between variability in zein amount and composition and the molecular polymorphism at this locus. Moreover, a QTL analysis revealed that the variability in Opaque-2 transcript abundance was controlled by several polymorphic trans-acting regulators unlinked to the Opaque-2 structural gene. Such genetic approaches should represent additional tools for physiological analysis of the huge amounts of data generated by transcritome and proteome projects.

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.

Optimisation of enzyme concentrations for unbranched reaction chains: the concept of combined response coefficient.

In the metabolic control theory, the control coefficient is a key parameter in quantifying the sensitivity of the flux towards an infinitesimal variation of enzyme activity. This concept does not apply just as it is for variations of enzyme concentrations whenever there is spatial, energy or resources limitations in the cell. Due to constraint on total enzyme concentration, the variation of concentration of any given enzyme may affect the concentrations of other enzymes. To take into account these correlations between enzyme concentrations, we propose the concept of "combined response coefficient". Its definition is similar to that of the control coefficient, but its mathematical expression is different. Its range of variation is from -infinity to +1, the null value corresponding to optimum enzyme concentration, i.e. to concentrations that maximise the flux, and the negative values to concentrations beyond the optimum value. A summation property could be derived using a simple weighting of the combined response coefficients, the sum of the weighed coefficient being 0.

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.

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.

Proteomics for genetic and physiological studies in plants.

Proteomics is becoming a necessity in plant biology, as it is in medicine, zoology and microbiology, for deciphering the function and role of the genes that are or will be sequenced. In this review we focus on the various, mainly genetic, applications of the proteomic tools that have been developed in recent years: characterization of individuals or lines, estimation of genetic variability within and between populations, establishment of genetic distances that can be used in phylogenetic studies, characterization of mutants and localization of the genes encoding the revealed proteins. Improvements in specifically devoted software have permitted precise quantification of the variation in amounts of proteins, leading to the concept of "protein quantity loci" which, combined with the "quantitative trait loci" approach, results in testable hypotheses regarding the role of "candidate proteins" in the metabolism or phenotype under study. This new development is exemplified by the reaction of plants to drought, a trait of major agronomic interest. The accumulation of data regarding genomic and cDNA sequencing will be connected to the protein databases currently developed in plants.

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.

Characterization of storage proteins in Daucus carota L.: two novel proteins display zygotic embryo specificity.

Although maturation-related proteins are well known in the endosperm of albuminous seeds, an important question is whether the zygotic embryo possesses its own maturation proteins. We report on the isolation and partial characterization of storage proteins of carrot (Daucus carota L. var Nandor) dry achenes and isolated zygotic embryos, using one- and two-dimensional electrophoresis techniques, HPLC and amino acid sequencing. The presence of a series of abundant polypeptides showing charge heterogeneity, that are rapidly degraded upon germination, was revealed in the endosperm. These proteins consisted of glycoproteins, the most abundant of which displayed a molecular mass (M(r)) of 58,000, albumins of M(r) 42,000 comprising at least one beta-1,3-glucanase, and two globulins of M(r) 90,000 and 50,000-55,000 respectively, the second being an oligomer composed of three subunits of M(r) 13,000, 20,000 and 30,000. None of these storage proteins identified in the endosperm were detected in zygotic embryos. In contrast, two novel proteins were isolated from zygotic embryos, namely a globulin family of M(r) 50,000 and pI 6.3-6.8, which was named "daucin", and a late embryogenesis abundant (LEA) protein family of M(r) 25,000 and pI 6.3-6.6, named "RAB25". Since the latter proteins are apparently absent of the endosperm, these results suggest that the maturation of carrot zygotic embryos requires its own specific set of storage and LEA proteins.

Somatic and zygotic embryos of Daucus carota L. display different protein patterns until conversion to plants.

Total protein patterns of different developmental stages of carrot zygotic and somatic embryos revealed by one- and two-dimensional gel electrophoresis were compared using statistical dissimilarity index matrix, and some major polypeptides were partially sequenced. In spite of similar morphology, the protein patterns of somatic embryos at the torpedo stage were clearly different from those of zygotic embryos. In particular, none of the proteins specific of zygotic embryos required for maturation, previously identified, were accumulated in somatic embryos, namely the daucin (a globulin-type storage protein), the RAB25 protein (a late embryogenesis abundant protein) (Dodeman et al. 1998), as well as a novel globulin of M(r) 30,000, that we proposed to name apiacin. Somatic plantlets and seedlings also showed different patterns. This discrepancy likely reflects culture conditions, since somatic embryos recover a protein pattern close to that of seedlings after conversion to plant and growth on a carbon-free medium.