Ecological and evolutionary genomics in the wild tomatoes (Solanum sect. Lycopersicon)

The promise of genomics in the study of plant-pollinator interactions

Flowers exist in exceedingly complex fitness landscapes, in which subtle variation in each trait can affect the pollinators, herbivores and pleiotropically linked traits in other plant tissues. A whole-genome approach to flower evolution will help our understanding of plant-pollinator interactions.

Identification of novel loci regulating interspecific variation in root morphology and cellular development in tomato

While the Arabidopsis (Arabidopsis thaliana) root has been elegantly characterized with respect to specification of cell identity, its development is missing a number of cellular features present in other species. We have characterized the root development of a wild and a domesticated tomato species, Solanum pennellii and Solanum lycopersicum 'M82.' We found extensive differences between these species for root morphology and cellular development including root length, a novel gravity set point angle, differences in cortical cell layer patterning, stem cell niche structure, and radial cell division. Using an introgression line population between these two species, we identified numerous loci that regulate these distinct aspects of development. Specifically we comprehensively identified loci that regulate (1) root length by distinct mechanisms including regulation of cell production within the meristem and the balance between cell division and expansion, (2) the gravity set point angle, and (3) radial cell division or expansion either in specific cell types or generally across multiple cell types. Our findings provide a novel perspective on the regulation of root growth and development between species. These loci have exciting implications with respect to regulation of drought resistance or salinity tolerance and regulation of root development in a family that has undergone domestication.

Role of an esterase in flavor volatile variation within the tomato clade

Tomato flavor is dependent upon a complex mixture of volatiles including multiple acetate esters. Red-fruited species of the tomato clade accumulate a relatively low content of acetate esters in comparison with the green-fruited species. We show that the difference in volatile ester content between the red- and green-fruited species is associated with insertion of a retrotransposon adjacent to the most enzymatically active member of a family of esterases. This insertion causes higher expression of the esterase, resulting in the reduced levels of multiple esters that are negatively correlated with human preferences for tomato. The insertion was evolutionarily fixed in the red-fruited species, suggesting that high expression of the esterase and consequent low ester content may provide an adaptive advantage in the ancestor of the red-fruited species. These results illustrate at a molecular level how closely related species exhibit major differences in volatile production by altering a volatile-associated catabolic activity.

What is needed for next-generation ecological and evolutionary genomics?

Ecological and evolutionary genomics (EEG) aims to link gene functions and genomic features to phenotypes and ecological factors. Although the rapid development of technologies allows central questions to be addressed at an unprecedented level of molecular detail, they do not alleviate one of the major challenges of EEG, which is that a large fraction of genes remains without any annotation. Here, we propose two solutions to this challenge. The first solution is in the form of a database that regroups associations between genes, organismal attributes and abiotic and biotic conditions. This database would result in an ecological annotation of genes by allowing cross-referencing across studies and taxa. Our second solution is to use new functional techniques to characterize genes implicated in the response to ecological challenges.

Structural homology in the Solanaceae: analysis of genomic regions in support of synteny studies in tomato, potato and pepper

We have analysed the structural homology in euchromatin regions of tomato, potato and pepper with special attention for the long arm of chromosome 2 (2L). Molecular organization and colinear junctions were delineated using multi-color BAC FISH analysis and comparative sequence alignment. We found large-scale rearrangements including inversions and segmental translocations that were not reported in previous comparative studies. Some of the structural rearrangements are specific for the tomato clade, and differentiate tomato from potato, pepper and other Solanaceous species. Although local gene vicinity is largely preserved, there are many small-scale synteny perturbations. Gene adjacency in the aligned segments was frequently disrupted for 47% of the ortholog pairs as a result of gene and LTR retrotransposon insertions, and occasionally by single gene inversions and translocations. Our data also suggests that long distance intra-chromosomal rearrangements and local gene rearrangements have evolved frequently during speciation in the Solanum genus, and that small changes are more prevalent than large-scale differences. The occurrence of sonata and harbinger transposable elements and other repeats near or at junction breaks is considered in the light of repeat-mediated rearrangements and a reconstruction scenario for an ancestral 2L topology is discussed.

Combined transcriptome, genetic diversity and metabolite profiling in tomato fruit reveals that the ethylene response factor SlERF6 plays an important role in ripening and carotenoid accumulation

Solanum lycopersicum (tomato) and its wild relatives harbor genetic diversity that yields heritable variation in fruit chemistry that could be exploited to identify genes regulating their synthesis and accumulation. Carotenoids, for example, are essential in plant and animal nutrition, and are the visual indicators of ripening for many fruits, including tomato. Whereas carotenoid synthesis is well characterized, factors regulating flux through the pathway are poorly understood at the molecular level. To exploit the impact of tomato genetic diversity on carotenoids, Solanum pennellii introgression lines were used as a source of defined natural variation and as a resource for the identification of candidate regulatory genes. Ripe fruits were analyzed for numerous fruit metabolites and transcriptome profiles generated using a 12,000 unigene oligoarray. Correlation analysis between carotenoid content and gene expression profiles revealed 953 carotenoid-correlated genes. To narrow the pool, subnetwork analysis of carotenoid-correlated transcription revealed 38 candidates. One candidate for impact on trans-lycopene and β-carotene accumulation was functionally charaterized, SlERF6, revealing that it indeed influences carotenoid biosynthesis and additional ripening phenotypes. Reduced expression of SlERF6 by RNAi enhanced both carotenoid and ethylene levels during fruit ripening, demonstrating an important role for SlERF6 in ripening, integrating the ethylene and carotenoid synthesis pathways.

The fruit cuticles of wild tomato species exhibit architectural and chemical diversity, providing a new model for studying the evolution of cuticle function

The cuticle covers the aerial epidermis of land plants and plays a primary role in water regulation and protection from external stresses. Remarkable species diversity in the structure and composition of its components, cutin and wax, have been catalogued, but few functional or genetic correlations have emerged. Tomato (Solanum lycopersicum) is part of a complex of closely related wild species endemic to the northern Andes and the Galapagos Islands (Solanum Sect. Lycopersicon). Although sharing an ancestor <7 million years ago, these species are found in diverse environments and are subject to unique selective pressures. Furthermore, they are genetically tractable, since they can be crossed with S. lycopersicum, which has a sequenced genome. With the aim of evaluating the relationships between evolution, structure and function of the cuticle, we characterized the morphological and chemical diversity of fruit cuticles of seven species from Solanum Sect. Lycopersicon. Striking differences in cuticular architecture and quantities of cutin and waxes were observed, with the wax coverage of wild species exceeding that of S. lycopersicum by up to seven fold. Wax composition varied in the occurrence of wax esters and triterpenoid isomers. Using a Solanum habrochaites introgression line population, we mapped triterpenoid differences to a genomic region that includes two S. lycopersicum triterpene synthases. Based on known metabolic pathways for acyl wax compounds, hypotheses are discussed to explain the appearance of wax esters with atypical chain lengths. These results establish a model system for understanding the ecological and evolutionary functional genomics of plant cuticles.

Native environment modulates leaf size and response to simulated foliar shade across wild tomato species

The laminae of leaves optimize photosynthetic rates by serving as a platform for both light capture and gas exchange, while minimizing water losses associated with thermoregulation and transpiration. Many have speculated that plants maximize photosynthetic output and minimize associated costs through leaf size, complexity, and shape, but a unifying theory linking the plethora of observed leaf forms with the environment remains elusive. Additionally, the leaf itself is a plastic structure, responsive to its surroundings, further complicating the relationship. Despite extensive knowledge of the genetic mechanisms underlying angiosperm leaf development, little is known about how phenotypic plasticity and selective pressures converge to create the diversity of leaf shapes and sizes across lineages. Here, we use wild tomato accessions, collected from locales with diverse levels of foliar shade, temperature, and precipitation, as a model to assay the extent of shade avoidance in leaf traits and the degree to which these leaf traits correlate with environmental factors. We find that leaf size is correlated with measures of foliar shade across the wild tomato species sampled and that leaf size and serration correlate in a species-dependent fashion with temperature and precipitation. We use far-red induced changes in leaf length as a proxy measure of the shade avoidance response, and find that shade avoidance in leaves negatively correlates with the level of foliar shade recorded at the point of origin of an accession. The direction and magnitude of these correlations varies across the leaf series, suggesting that heterochronic and/or ontogenic programs are a mechanism by which selective pressures can alter leaf size and form. This study highlights the value of wild tomato accessions for studies of both morphological and light-regulated development of compound leaves, and promises to be useful in the future identification of genes regulating potentially adaptive plastic leaf traits.

Adaptation to drought in two wild tomato species: the evolution of the Asr gene family

Wild tomato species are a valuable system in which to study local adaptation to drought: they grow in diverse environments ranging from mesic to extremely arid conditions. Here, we investigate the evolution of members of the Asr (ABA/water stress/ripening induced) gene family, which have been reported to be involved in the water stress response. We analysed molecular variation in the Asr gene family in populations of two closely related species, Solanum chilense and Solanum peruvianum. We concluded that Asr1 has evolved under strong purifying selection. In contrast to previous reports, we did not detect evidence for positive selection at Asr2. However, Asr4 shows patterns consistent with local adaptation in an S. chilense population that lives in an extremely dry environment. We also discovered a new member of the gene family, Asr5. Our results show that the Asr genes constitute a dynamic gene family and provide an excellent example of tandemly arrayed genes that are of importance in adaptation. Taking the potential distribution of the species into account, it appears that S. peruvianum can cope with a great variety of environmental conditions without undergoing local adaptation, whereas S. chilense undergoes local adaptation more frequently.

MTML-msBayes: approximate Bayesian comparative phylogeographic inference from multiple taxa and multiple loci with rate heterogeneity

Background

MTML-msBayes uses hierarchical approximate Bayesian computation (HABC) under a coalescent model to infer temporal patterns of divergence and gene flow across codistributed taxon-pairs. Under a model of multiple codistributed taxa that diverge into taxon-pairs with subsequent gene flow or isolation, one can estimate hyper-parameters that quantify the mean and variability in divergence times or test models of migration and isolation. The software uses multi-locus DNA sequence data collected from multiple taxon-pairs and allows variation across taxa in demographic parameters as well as heterogeneity in DNA mutation rates across loci. The method also allows a flexible sampling scheme: different numbers of loci of varying length can be sampled from different taxon-pairs.

Results

Simulation tests reveal increasing power with increasing numbers of loci when attempting to distinguish temporal congruence from incongruence in divergence times across taxon-pairs. These results are robust to DNA mutation rate heterogeneity. Estimating mean divergence times and testing simultaneous divergence was less accurate with migration, but improved if one specified the correct migration model. Simulation validation tests demonstrated that one can detect the correct migration or isolation model with high probability, and that this HABC model testing procedure was greatly improved by incorporating a summary statistic originally developed for this task (Wakeley's Ψ_W). The method is applied to an empirical data set of three Australian avian taxon-pairs and a result of simultaneous divergence with some subsequent gene flow is inferred.

Conclusions

To retain flexibility and compatibility with existing bioinformatics tools, MTML-msBayes is a pipeline software package consisting of Perl, C and R programs that are executed via the command line. Source code and binaries are available for download at http://msbayes.sourceforge.net/ under an open source license (GNU Public License).

Evolutionary and Ecological Genomics of Non-Model Plants

Dissecting evolutionary dynamics of ecologically important traits is a long-term challenge for biologists. Attempts to understand natural variation and molecular mechanisms have motivated a move from laboratory model systems to non-model systems in diverse natural environments. Next generation sequencing methods, along with an expansion of genomic resources and tools, have fostered new links between diverse disciplines, including molecular biology, evolution, and ecology, and genomics. Great progress has been made in a few non-model wild plants, such as Arabidopsis relatives, monkey flowers, and wild sunflowers. Until recently, the lack of comprehensive genomic information has limited evolutionary and ecological studies to larger QTL regions rather than single gene resolution, and has hindered recognition of general patterns of natural variation and local adaptation. Further efforts in accumulating genomic data and developing bioinformatic and biostatistical tools are now poised to move this field forward. Integrative national and international collaborations and research communities are needed to facilitate development in the field of evolutionary and ecological genomics.

Functional gametophytic self-incompatibility in a peripheral population of Solanum peruvianum (Solanaceae)

The transition from self-incompatibility to self-compatibility is a common transition in angiosperms often reported in populations at the edge of species range limits. Geographically distinct populations of wild tomato species (Solanum section Lycopersicon (Solanaceae)) have been described as polymorphic for mating system with both self-incompatible and self-compatible populations. Using controlled pollinations and sequencing of the S-RNase mating system gene, we test the compatibility status of a population of S. peruvianum located near its southern range limit. Pollinations among plants of known genotypes revealed strong self-incompatibility; fruit set following compatible pollinations was significantly higher than following incompatible pollinations for all tested individuals. Sequencing of the S-RNase gene in parents and progeny arrays was also as predicted under self-incompatibility. Molecular variation at the S-RNase locus revealed a diverse set of alleles, and heterozygosity in over 500 genotyped individuals. We used controlled crosses to test the specificity of sequences recovered in this study; in all cases, results were consistent with a unique allelic specificity for each tested sequence, including two alleles sharing 92% amino-acid similarity. Site-specific patterns of selection at the S-RNase gene indicate positive selection in regions of the gene associated with allelic specificity determination and purifying selection in previously characterized conserved regions. Further, there is broad convergence between the present and previous studies in specific amino-acid positions inferred to be evolving under positive selection.

Nucleotide diversity patterns of local adaptation at drought-related candidate genes in wild tomatoes

We surveyed nucleotide diversity at two candidate genes LeNCED1 and pLC30-15, involved in an ABA (abscisic acid) signalling pathway, in two closely related tomato species Solanum peruvianum and Solanum chilense. Our six population samples (three for each species) cover a range of mesic to very dry habitats. The ABA pathway plays an important role in the plants' response to drought stress. LeNCED1 is an upstream gene involved in ABA biosynthesis, and pLC30-15 is a dehydrin gene positioned downstream in the pathway. The two genes show very different patterns of nucleotide variation. LeNCED1 exhibits very low nucleotide diversity relative to the eight neutral reference loci that were previously surveyed in these populations. This suggests that strong purifying selection has been acting on this gene. In contrast, pLC30-15 exhibits higher levels of nucleotide diversity and, in particular in S. chilense, higher genetic differentiation between populations than the reference loci, which is indicative of local adaptation. In the more drought-tolerant species S. chilense, one population (from Quicacha) shows a significant haplotype structure, which appears to be the result of positive (diversifying) selection.

Reciprocal insights into adaptation from agricultural and evolutionary studies in tomato

Although traditionally separated by different aims and methodologies, research on agricultural and evolutionary problems shares a common goal of understanding the mechanisms underlying functionally important traits. As such, research in both fields offers potential complementary and reciprocal insights. Here, we discuss adaptive stress responses (specifically to water stress) as an example of potentially fruitful research reciprocity, where agricultural research has clearly produced advances that could benefit evolutionary studies, while evolutionary studies offer approaches and insights underexplored in crop studies. We focus on research on Solanum species that include the domesticated tomato and its wild relatives. Integrated approaches to understanding ecological adaptation are particularly attractive in tomato and its wild relatives: many presumptively adaptive phenotypic differences characterize wild species, and the physiological and mechanistic basis of many relevant traits and environmental responses has already been examined in the context of cultivated tomato and some wild species. We highlight four specific instances where these reciprocal insights can be combined to better address questions that are fundamental both to agriculture and evolution.

Sexual and apomictic plant reproduction in the genomics era: exploring the mechanisms potentially useful in crop plants

Arabidopsis, Mimulus and tomato have emerged as model plants in researching genetic and molecular basis of differences in mating systems. Variations in floral traits and loss of self-incompatibility have been associated with mating system differences in crops. Genomics research has advanced considerably, both in model and crop plants, which may provide opportunities to modify breeding systems as evidenced in Arabidopsis and tomato. Mating system, however, not recombination per se, has greater effect on the level of polymorphism. Generating targeted recombination remains one of the most important factors for crop genetic enhancement. Asexual reproduction through seeds or apomixis, by producing maternal clones, presents a tremendous potential for agriculture. Although believed to be under simple genetic control, recent research has revealed that apomixis results as a consequence of the deregulation of the timing of sexual events rather than being the product of specific apomixis genes. Further, forward genetic studies in Arabidopsis have permitted the isolation of novel genes reported to control meiosis I and II entry. Mutations in these genes trigger the production of unreduced or apomeiotic megagametes and are an important step toward understanding and engineering apomixis.

SolEST database: a "one-stop shop" approach to the study of Solanaceae transcriptomes

BACKGROUND

Since no genome sequences of solanaceous plants have yet been completed, expressed sequence tag (EST) collections represent a reliable tool for broad sampling of Solanaceae transcriptomes, an attractive route for understanding Solanaceae genome functionality and a powerful reference for the structural annotation of emerging Solanaceae genome sequences.

DESCRIPTION

We describe the SolEST database http://biosrv.cab.unina.it/solestdb which integrates different EST datasets from both cultivated and wild Solanaceae species and from two species of the genus Coffea. Background as well as processed data contained in the database, extensively linked to external related resources, represent an invaluable source of information for these plant families. Two novel features differentiate SolEST from other resources: i) the option of accessing and then visualizing Solanaceae EST/TC alignments along the emerging tomato and potato genome sequences; ii) the opportunity to compare different Solanaceae assemblies generated by diverse research groups in the attempt to address a common complaint in the SOL community.

CONCLUSION

Different databases have been established worldwide for collecting Solanaceae ESTs and are related in concept, content and utility to the one presented herein. However, the SolEST database has several distinguishing features that make it appealing for the research community and facilitates a "one-stop shop" for the study of Solanaceae transcriptomes.

Linking genotypes to phenotypes and fitness: how mechanistic biology can inform molecular ecology

The accessibility of new genomic resources, high-throughput molecular technologies and analytical approaches such as genome scans have made finding genes contributing to fitness variation in natural populations an increasingly feasible task. Once candidate genes are identified, we argue that it is necessary to take a mechanistic approach and work up through the levels of biological organization to fully understand the impacts of genetic variation at these candidate genes. We demonstrate how this approach provides testable hypotheses about the causal links among levels of biological organization, and assists in designing relevant experiments to test the effects of genetic variation on phenotype, whole-organism performance capabilities and fitness. We review some of the research programs that have incorporated mechanistic approaches when examining naturally occurring genetic and phenotypic variation and use these examples to highlight the value of developing a comprehensive understanding of the relationship between genotype and fitness. We give suggestions to guide future research aimed at uncovering and understanding the genetic basis of adaptation and argue that further integration of mechanistic approaches will help molecular ecologists better understand the evolution of natural populations.

Antagonistic epistasis for ecophysiological trait differences between Solanum species

Epistasis, the nonadditive interaction between loci, is thought to play a role in many fundamental evolutionary processes, including adaptive differentiation and speciation. Focusing on species differences in ecophysiological traits, we examined the strength and direction of pairwise epistatic interactions between target chromosomal regions from one species, when co-introgressed into the genetic background of a foreign species. A full diallel cross was performed using 15 near-isogenic lines (NILs) constructed between two tomato species (Solanum habrochaites and Solanum lycopersicum) to compare the phenotypic effects of each chromosomal region singly and in combination with each other region. We detected main effect quantitative trait loci (QTLs) for two of our three focal traits. Epistatic effects accounted for c. 25% of detected effects on trait means, depending on the trait. Strikingly, all but two interactions were antagonistic, with the combined effect of chromosomal regions acting in the opposite direction from that of one or both individual chromosomal regions. Our study is one of the few to systematically examine pairwise epistatic effects in a nonmicrobial system. Our results suggest that epistatic interactions can contribute substantially to the genetic basis of traits involved in adaptive species differentiation, especially highly complex, multivariate traits.