Ecological and geographic modes of species divergence in wild tomatoes

Insights from ecological niche modeling on the taxonomic distinction and niche differentiation between the black-spotted and red-spotted tokay geckoes (Gekko gecko)

The black-spotted tokay and the red-spotted tokay are morphologically distinct and have largely allopatric distributions. The black-spotted tokay is characterized by a small body size and dark skin with sundry spots, while the red-spotted tokay has a relatively large body size and red spots. Based on morphological, karyotypic, genetic, and distribution differences, recent studies suggested their species status; however, their classifications remain controversial, and additional data such as ecological niches are necessary to establish firm hypotheses regarding their taxonomic status. We reconstructed their ecological niches models using climatic and geographic data. We then performed niche similarity tests (niche identity and background tests) and point-based analyses to explore whether ecological differentiation has occurred, and whether such differences are sufficient to explain the maintenance of their separate segments of environmental ranges. We found that both niche models of the black- and the red-spotted tokay had a good fit and a robust performance, as indicated by the high area under the curve (AUC) values (“black” = 0.982, SD = ± 0.002, “red” = 0.966 ± 0.02). Significant ecological differentiation across the entire geographic range was found, indicating that the involvement of ecological differentiation is important for species differentiation. Divergence along the environmental axes is highly associated with climatic conditions, with isothermality being important for the “black” form, while temperature seasonality, precipitation of warmest quarter, and annual temperature range together being important for the “red” form. These factors are likely important factors in niche differentiation between the two forms, which result in morphological replacement. Overall, beside morphological and genetic differentiation information, our results contribute to additional insights into taxonomic distinction and niche differentiation between the black- and the red-spotted tokay.

Morphological and anatomical determinants of mesophyll conductance in wild relatives of tomato (Solanum sect. Lycopersicon, sect. Lycopersicoides; Solanaceae)

Natural selection on photosynthetic performance is a primary factor determining leaf phenotypes. The complex CO2 diffusion path from substomatal cavities to the chloroplasts - the mesophyll conductance (g(m)) - limits photosynthetic rate in many species and hence shapes variation in leaf morphology and anatomy. Among sclerophyllous and succulent taxa, structural investment in leaves, measured as the leaf dry mass per area (LMA), has been implicated in decreased gm . However, in herbaceous taxa with high g(m), it is less certain how LMA impacts CO2 diffusion and whether it significantly affects photosynthetic performance. We addressed these questions in the context of understanding the ecophysiological significance of leaf trait variation in wild tomatoes, a closely related group of herbaceous perennials. Although g(m) was high in wild tomatoes, variation in g(m) significantly affected photosynthesis. Even in these tender-leaved herbaceous species, greater LMA led to reduced g(m). This relationship between g(m) and LMA is partially mediated by cell packing and leaf thickness, although amphistomy (equal distribution of stomata on both sides of the leaf) mitigates the effect of leaf thickness. Understanding the costs of increased LMA will inform future work on the adaptive significance of leaf trait variation across ecological gradients in wild tomatoes and other systems.

Genetic structure of the four wild tomato species in the Solanum peruvianum s.l. species complex

The most diverse wild tomato species Solanum peruvianum sensu lato (s.l.) has been reclassified into four separate species: Solanum peruvianum sensu stricto (s.s.), Solanum corneliomuelleri, Solanum huaylasense, and Solanum arcanum. However, reproductive barriers among the species are incomplete and this can lead to discrepancies regarding genetic identity of germplasm. We used genotyping by sequencing (GBS) of S. peruvianum s.l., Solanum neorickii, and Solanum chmielewskii to develop tens of thousands of mapped single nucleotide polymorphisms (SNPs) to analyze genetic relationships within and among species. The data set was condensed to 14,043 SNPs with no missing data across 46 sampled plants. Origins of accessions were mapped using geographical information systems (GIS). Isolation by distance, pairwise genetic distances, and number of clusters were estimated using population genetics approaches. Isolation by distance was strongly supported, especially between interspecific pairs. Eriopersicon (S. peruvianum s.s., S. corneliomuelleri, S. huaylasense) and Arcanum (S. arcanum, S. neorickii, S. chmielewskii) species groups were genetically distinct, except for S. huaylasense which showed 50% membership proportions in each group. Solanum peruvianum and S. corneliomuelleri were not significantly differentiated from each other. Many thousands of SNP markers were identified that could potentially be used to distinguish pairs of species, including S. peruvianum versus S. corneliomuelleri, if they are verified on larger numbers of samples. Diagnostic markers will be valuable for delimiting morphologically similar and interfertile species in germplasm management. Approximately 12% of the SNPs rejected a genome-wide test of selective neutrality based on differentiation among species of S. peruvianum s.l. These are candidates for more comprehensive studies of microevolutionary processes within this species complex.

The geography and ecology of plant speciation: range overlap and niche divergence in sister species

A goal of evolutionary biology is to understand the roles of geography and ecology in speciation. The recent shared ancestry of sister species can leave a major imprint on their geographical and ecological attributes, possibly revealing processes involved in speciation. We examined how ecological similarity, range overlap and range asymmetry are related to time since divergence of 71 sister species pairs in the California Floristic Province (CFP). We found that plants exhibit strikingly different age-range correlation patterns from those found for animals; the latter broadly support allopatric speciation as the primary mode of speciation. By contrast, plant sisters in the CFP were sympatric in 80% of cases and range sizes of sisters differed by a mean of 10-fold. Range overlap and range asymmetry were greatest in younger sisters. These results suggest that speciation mechanisms broadly grouped under 'budding' speciation, in which a larger ranged progenitor gives rise to a smaller ranged derivative species, are probably common. The ecological and reproductive similarity of sisters was significantly greater than that of sister-non-sister congeners for every trait assessed. However, shifts in at least one trait were present in 93% of the sister pairs; habitat and soil shifts were especially common. Ecological divergence did not increase with range overlap contrary to expectations under character displacement in sympatry. Our results suggest that vicariant speciation is more ubiquitous in animals than plants, perhaps owing to the sensitivity of plants to fine-scale environmental heterogeneity. Despite high levels of range overlap, ecological shifts in the process of budding speciation may result in low rates of fine-scale spatial co-occurrence. These results have implications for ecological studies of trait evolution and community assembly; despite high levels of sympatry, sister taxa and potentially other close relatives, may be missing from local communities.

Restoring pistil-side self-incompatibility factors recapitulates an interspecific reproductive barrier between tomato species

Interspecific reproductive barriers are poorly understood, but are central to the biological species concept. The pre-zygotic barriers between red- and green-fruited species in the tomato clade of the genus Solanum provide a model to better understand these barriers in plants. Compatibility usually follows the SI x SC rule: pollen from self-compatible (SC) red-fruited species is rejected on pistils of the predominantly self-incompatible (SI) green-fruited species, but the reciprocal crosses are compatible. This suggests that the interspecific reproductive barrier may be linked to the intraspecific SI mechanism. However, pollen from the SC red-fruited species is also rejected by SC accessions of green-fruited species that lack S-RNase, a key protein expressed in pistils of SI Solanum species. Thus, multiple mechanisms may contribute to the barrier between red- and green-fruited species. We tested whether an S-RNase-dependent barrier is sufficient for rejection of pollen from red-fruited species by introducing functional S-RNase, HT-A and HT-B genes from SI species into Solanum lycopersicum (cultivated tomato). We found that expressing S-RNase in combination with either HT-A or HT-B in the pistil is sufficient to cause rejection of pollen from all four red-fruited species. Thus, redundant mechanisms must operate side by side to prevent crosses between red- and green-fruited species in the clade, underlining the complexity of interspecific pollination barriers. Our results also have implications for mating system transitions. We suggest that these transitions must occur in a specific sequence, and that the transition from SI to SC also affects interspecific compatibility.

Niche overlap of congeneric invaders supports a single-species hypothesis and provides insight into future invasion risk: implications for global management of the Bactrocera dorsalis complex

Background

The invasive fruit fly, Bactrocera invadens, has expanded its range rapidly over the past 10 years. Here we aimed to determine if the recent range expansion of Bactrocera invadens into southern Africa can be better understood through niche exploration tools, ecological niche models (ENMs), and through incorporating information about Bactrocera dorsalis s.s., a putative conspecific species from Asia. We test for niche overlap of environmental variables between Bactrocera invadens and Bactrocera dorsalis s.s. as well as two other putative conspecific species, Bactrocera philippinensis and B. papayae. We examine overlap and similarity in the geographical expression of each species’ realised niche through reciprocal distribution models between Africa and Asia. We explore different geographical backgrounds, environmental variables and model complexity with multiple and single Bactrocera species hypotheses in an attempt to predict the recent range expansion of B. invadens into northern parts of South Africa.

Principal Findings

Bactrocera invadens has a high degree of niche overlap with B. dorsalis s.s. (and B. philippinensis and B. papayae). Ecological niche models built for Bactrocera dorsalis s.s. have high transferability to describe the range of B. invadens, and B. invadens is able to project to the core range of B. dorsalis s.s. The ENMs of both Bactrocera dorsalis and B. dorsalis combined with B. philipenesis and B. papayae have significantly higher predictive ability to capture the distribution points in South Africa than for B. invadens alone.

Conclusions/Significance

Consistent with other studies proposing these Bactrocera species as conspecific, niche similarity and overlap between these species is high. Considering these other Bactrocera dorsalis complex species simultaneously better describes the range expansion and invasion potential of B. invadens in South Africa. We suggest that these species should be considered the same–at least functionally–and global quarantine and management strategies applied equally to these Bactrocera species.

Designing optimized multi-species monitoring networks to detect range shifts driven by climate change: a case study with bats in the North of Portugal

Here we develop a framework to design multi-species monitoring networks using species distribution models and conservation planning tools to optimize the location of monitoring stations to detect potential range shifts driven by climate change. For this study, we focused on seven bat species in Northern Portugal (Western Europe). Maximum entropy modelling was used to predict the likely occurrence of those species under present and future climatic conditions. By comparing present and future predicted distributions, we identified areas where each species is likely to gain, lose or maintain suitable climatic space. We then used a decision support tool (the Marxan software) to design three optimized monitoring networks considering: a) changes in species likely occurrence, b) species conservation status, and c) level of volunteer commitment. For present climatic conditions, species distribution models revealed that areas suitable for most species occur in the north-eastern part of the region. However, areas predicted to become climatically suitable in the future shifted towards west. The three simulated monitoring networks, adaptable for an unpredictable volunteer commitment, included 28, 54 and 110 sampling locations respectively, distributed across the study area and covering the potential full range of conditions where species range shifts may occur. Our results show that our framework outperforms the traditional approach that only considers current species ranges, in allocating monitoring stations distributed across different categories of predicted shifts in species distributions. This study presents a straightforward framework to design monitoring schemes aimed specifically at testing hypotheses about where and when species ranges may shift with climatic changes, while also ensuring surveillance of general population trends.

Evidence for large-scale effects of competition: niche displacement in Canada lynx and bobcat

Determining the patterns, causes and consequences of character displacement is central to our understanding of competition in ecological communities. However, the majority of competition research has occurred over small spatial extents or focused on fine-scale differences in morphology or behaviour. The effects of competition on broad-scale distribution and niche characteristics of species remain poorly understood but critically important. Using range-wide species distribution models, we evaluated whether Canada lynx (Lynx canadensis) or bobcat (Lynx rufus) were displaced in regions of sympatry. Consistent with our prediction, we found that lynx niches were less similar to those of bobcat in areas of sympatry versus allopatry, with a stronger reliance on snow cover driving lynx niche divergence in the sympatric zone. By contrast, bobcat increased niche breadth in zones of sympatry, and bobcat niches were equally similar to those of lynx in zones of sympatry and allopatry. These findings suggest that competitively disadvantaged species avoid competition at large scales by restricting their niche to highly suitable conditions, while superior competitors expand the diversity of environments used. Our results indicate that competition can manifest within climatic niche space across species' ranges, highlighting the importance of biotic interactions occurring at large spatial scales on niche dynamics.

Unveiling current Guanaco distribution in chile based upon niche structure of phylogeographic lineages: Andean puna to subpolar forests

Niche description and differentiation at broad geographic scales have been recent major topics in ecology and evolution. Describing the environmental niche structure of sister taxa with known evolutionary trajectories stands out as a useful exercise in understanding niche requirements. Here we model the environmental niche structure and distribution of the recently resolved phylogeography of guanaco (Lama guanicoe) lineages on the western slope of the southern Andes. Using a maximum entropy framework, field data, and information on climate, topography, human density, and vegetation cover, we identify differences between the two subspecies (L.g.cacsilensis, L.g.guanicoe) and their intermediate-hybrid lineage, that most likely determine the distribution of this species. While aridity seems to be a major factor influencing the distribution at the species-level (annual precipitation <900 mm), we also document important differences in niche specificity for each subspecies, where distribution of Northern lineage is explained mainly by elevation (mean = 3,413 m) and precipitation seasonality (mean = 161 mm), hybrid lineage by annual precipitation (mean = 139 mm), and Southern subspecies by annual precipitation (mean = 553 mm), precipitation seasonality (mean = 21 mm) and grass cover (mean = 8.2%). Among lineages, we detected low levels of niche overlap: I (Similarity Index) = 0.06 and D (Schoener's Similarity Index) = 0.01; and higher levels when comparing Northern and Southern subspecies with hybrids lineage ( I = 0.32-0.10 and D = 0.12-0.03, respectively). This suggests that important ecological and/or evolutionary processes are shaping the niche of guanacos in Chile, producing discrepancies when comparing range distribution at the species-level (81,756 km(2)) with lineages-level (65,321 km(2)). The subspecies-specific description of niche structure is provided here based upon detailed spatial distribution of the lineages of guanacos in Chile. Such description provides a scientific tool to further develop large scale plans for habitat conservation and preservation of intraspecific genetic variability for this far ranging South American camelid, which inhabits a diversity of ecoregion types from Andean puna to subpolar forests.

Sequence evolution and expression regulation of stress-responsive genes in natural populations of wild tomato

The wild tomato species Solanum chilense and S. peruvianum are a valuable non-model system for studying plant adaptation since they grow in diverse environments facing many abiotic constraints. Here we investigate the sequence evolution of regulatory regions of drought and cold responsive genes and their expression regulation. The coding regions of these genes were previously shown to exhibit signatures of positive selection. Expression profiles and sequence evolution of regulatory regions of members of the Asr (ABA/water stress/ripening induced) gene family and the dehydrin gene pLC30-15 were analyzed in wild tomato populations from contrasting environments. For S. chilense, we found that Asr4 and pLC30-15 appear to respond much faster to drought conditions in accessions from very dry environments than accessions from more mesic locations. Sequence analysis suggests that the promoter of Asr2 and the downstream region of pLC30-15 are under positive selection in some local populations of S. chilense. By investigating gene expression differences at the population level we provide further support of our previous conclusions that Asr2, Asr4, and pLC30-15 are promising candidates for functional studies of adaptation. Our analysis also demonstrates the power of the candidate gene approach in evolutionary biology research and highlights the importance of wild Solanum species as a genetic resource for their cultivated relatives.

Chilling-induced water stress: variation in shoot turgor maintenance among wild tomato species from diverse habitats

PREMISE OF THE STUDY

Cultivated tomato, Solanum lycopersicum, suffers chilling induced wilting because water movement through its roots decreases with declining soil temperatures. Certain wild tomato species exhibit resistance to chilling-induced wilting, but the extent of this chilling tolerance in wild tomatoes is not known. •

METHODS

We measured shoot wilting during root chilling in wild Solanum accessions from habitats differing in elevation, temperature, and precipitation. We also measured shoot wilting during root chilling in introgression lines (ILs) with chromosome 9 segments collinear to the shoot turgor maintenance QTL stm9 region from chilling-tolerant S. habrochaites, chilling and drought-tolerant S. lycopersicoides, or drought-tolerant S. pennellii. •

KEY RESULTS

Wild tomato species, which experience chilling temperatures (<10°C) in their native habitat, maintain shoot turgor under root chilling. Among accessions of S. lycopersicum var. cerasiforme, a typically chilling sensitive species, shoot turgor maintenance during root chilling was correlated with the precipitation of the native habitat. By contrast, S. pennellii, a species that is typically drought adapted, did not maintain turgor under root chilling. Grafted plants with roots containing S. habrochaites and S. lycopersicoides introgressions improved shoot turgor maintenance under root chilling. •

CONCLUSIONS

Resistance to chilling-induced water stress is an important adaptation to chilling temperatures in wild tomatoes. There is some overlap in adaptation to drought and chilling stress in some tomato species. Root-based resistance to chilling-induced water stress in wild tomatoes may involve orthologous gene(s).

Differences in insect resistance between tomato species endemic to the Galapagos Islands

BACKGROUND

The Galapagos Islands constitute a highly diverse ecosystem and a unique source of variation in the form of endemic species. There are two endemic tomato species, Solanum galapagense and S. cheesmaniae and two introduced tomato species, S. pimpinellifolium and S. lycopersicum. Morphologically the two endemic tomato species of the Galapagos Islands are clearly distinct, but molecular marker analysis showed no clear separation. Tomatoes on the Galapagos are affected by both native and exotic herbivores. Bemisia tabaci is an important introduced insect species that feeds on a wide range of plants. In this article, we address the question whether the differentiation between S. galapagense and S. cheesmaniae may be related to differences in susceptibility towards phloem-feeders and used B. tabaci as a model to evaluate this.

RESULTS

We have characterized 12 accessions of S. galapagense, 22 of S. cheesmaniae, and one of S. lycopersicum as reference for whitefly resistance using no-choice experiments. Whitefly resistance was found in S. galapagense only and was associated with the presence of relatively high levels of acyl sugars and the presence of glandular trichomes of type I and IV. Genetic fingerprinting using 3316 SNP markers did not show a clear differentiation between the two endemic species. Acyl sugar accumulation as well as the climatic and geographical conditions at the collection sites of the accessions did not follow the morphological species boundaries.

CONCLUSION

Our results suggest that S. galapagense and S. cheesmaniae might be morphotypes rather than two species and that their co-existence is likely the result of selective pressure.

Comparative transcriptomics reveals patterns of selection in domesticated and wild tomato

Although applied over extremely short timescales, artificial selection has dramatically altered the form, physiology, and life history of cultivated plants. We have used RNAseq to define both gene sequence and expression divergence between cultivated tomato and five related wild species. Based on sequence differences, we detect footprints of positive selection in over 50 genes. We also document thousands of shifts in gene-expression level, many of which resulted from changes in selection pressure. These rapidly evolving genes are commonly associated with environmental response and stress tolerance. The importance of environmental inputs during evolution of gene expression is further highlighted by large-scale alteration of the light response coexpression network between wild and cultivated accessions. Human manipulation of the genome has heavily impacted the tomato transcriptome through directed admixture and by indirectly favoring nonsynonymous over synonymous substitutions. Taken together, our results shed light on the pervasive effects artificial and natural selection have had on the transcriptomes of tomato and its wild relatives.

Evolution of a complex locus for terpene biosynthesis in solanum

Functional gene clusters, containing two or more genes encoding different enzymes for the same pathway, are sometimes observed in plant genomes, most often when the genes specify the synthesis of specialized defensive metabolites. Here, we show that a cluster of genes in tomato (Solanum lycopersicum; Solanaceae) contains genes for terpene synthases (TPSs) that specify the synthesis of monoterpenes and diterpenes from cis-prenyl diphosphates, substrates that are synthesized by enzymes encoded by cis-prenyl transferase (CPT) genes also located within the same cluster. The monoterpene synthase genes in the cluster likely evolved from a diterpene synthase gene in the cluster by duplication and divergence. In the orthologous cluster in Solanum habrochaites, a new sesquiterpene synthase gene was created by a duplication event of a monoterpene synthase followed by a localized gene conversion event directed by a diterpene synthase gene. The TPS genes in the Solanum cluster encoding cis-prenyl diphosphate-utilizing enzymes are closely related to a tobacco (Nicotiana tabacum; Solanaceae) diterpene synthase encoding Z-abienol synthase (Nt-ABS). Nt-ABS uses the substrate copal-8-ol diphosphate, which is made from the all-trans geranylgeranyl diphosphate by copal-8-ol diphosphate synthase (Nt-CPS2). The Solanum gene cluster also contains an ortholog of Nt-CPS2, but it appears to encode a nonfunctional protein. Thus, the Solanum functional gene cluster evolved by duplication and divergence of TPS genes, together with alterations in substrate specificity to utilize cis-prenyl diphosphates and through the acquisition of CPT genes.

Bioclimatic, ecological, and phenotypic intermediacy and high genetic admixture in a natural hybrid of octoploid strawberries

PREMISE OF THE STUDY

Hybrid zones provide "natural laboratories" for understanding the processes of selection, reinforcement, and speciation. We sought to gain insight into the degree of introgression and the extent of ecological-phenotypic intermediacy in the natural hybrid strawberry, Fragaria × ananassa subsp. cuneifolia. •

METHODS

We used whole-plastome sequencing to identify parental species-specific (Fragaria chiloensis and F. virginiana) chloroplast single-nucleotide polymorphisms and combined the use of these with nuclear microsatellite markers to genetically characterize the hybrid zone. We assessed the potential role of selection in the observed geographic patterns by bioclimatically characterizing the niche of the hybrid populations and phenotypically characterizing hybrid individuals of known genomic constitution. •

KEY RESULTS

Significant admixture and little overall maternal bias in chloroplast or nuclear genomes suggest a high degree of interfertility among the parental and hybrid species and point to a long history of backcrossing and genetic mixing in the hybrid zone. Even though hybrids were phenotypically intermediate to the parental species, there was a discernible fingerprint of the parental genotype within hybrid individuals. Thus, although the pattern of introgression observed suggests geographic limitations to gene flow, it may be reinforced by selection for specific parental traits in the bioclimatically intermediate habitat occupied by the hybrid. •

CONCLUSIONS

This work uncovered the genetic complexity underlying the hybrid zone of the wild relatives of the cultivated strawberry. It lays the foundation for experimental dissection of the causes of genomic introgression and nuclear-cytoplasmic disassociation, and for understanding other parts of Fragaria evolutionary history.

Reconsidering the specialist-generalist paradigm in niche breadth dynamics: resource gradient selection by Canada lynx and bobcat

The long-standing view in ecology is that disparity in overall resource selection is the basis for identifying niche breadth patterns, with species having narrow selection being classified "specialists" and those with broader selection being "generalists". The standard model of niche breadth characterizes generalists and specialists as having comparable levels of overall total resource exploitation, with specialists exploiting resources at a higher level of performance over a narrower range of conditions. This view has gone largely unchallenged. An alternate model predicts total resource use being lower for the specialized species with both peaking at a comparable level of performance over a particular resource gradient. To reconcile the niche breadth paradigm we contrasted both models by developing range-wide species distribution models for Canada lynx, Lynx canadensis, and bobcat, Lynx rufus. Using a suite of environmental factors to define each species' niche, we determined that Canada lynx demonstrated higher total performance over a restricted set of variables, specifically those related to snow and altitude, while bobcat had higher total performance across most variables. Unlike predictions generated by the standard model, bobcat level of exploitation was not compromised by the trade-off with peak performance, and Canada lynx were not restricted to exploiting a narrower range of conditions. Instead, the emergent pattern was that specialist species have a higher total resource utilization and peak performance value within a smaller number of resources or environmental axes than generalists. Our results also indicate that relative differences in niche breadth are strongly dependent on the variable under consideration, implying that the appropriate model describing niche breadth dynamics between specialists and generalists may be more complex than either the traditional heuristic or our modified version. Our results demonstrate a need to re-evaluate traditional, but largely untested, assumptions regarding resource utilization in species with broad and narrow niches.

Using morphological, molecular and climatic data to delimitate yews along the Hindu Kush-Himalaya and adjacent regions

Background

Despite the availability of several studies to clarify taxonomic problems on the highly threatened yews of the Hindu Kush-Himalaya (HKH) and adjacent regions, the total number of species and their exact distribution ranges remains controversial. We explored the use of comprehensive sets of morphological, molecular and climatic data to clarify taxonomy and distributions of yews in this region.

Methodology/Principal Findings

A total of 743 samples from 46 populations of wild yew and 47 representative herbarium specimens were analyzed. Principle component analyses on 27 morphological characters and 15 bioclimatic variables plus altitude and maximum parsimony analysis on molecular ITS and trnL-F sequences indicated the existence of three distinct species occurring in different ecological (climatic) and altitudinal gradients along the HKH and adjacent regions Taxus contorta from eastern Afghanistan to the eastern end of Central Nepal, T. wallichiana from the western end of Central Nepal to Northwest China, and the first report of the South China low to mid-elevation species T. mairei in Nepal, Bhutan, Northeast India, Myanmar and South Vietnam.

Conclusion/Significance

The detailed sampling and combination of different data sets allowed us to identify three clearly delineated species and their precise distribution ranges in the HKH and adjacent regions, which showed no overlap or no distinct hybrid zone. This might be due to differences in the ecological (climatic) requirements of the species. The analyses further provided the selection of diagnostic morphological characters for the identification of yews occurring in the HKH and adjacent regions. Our work demonstrates that extensive sampling combined with the analysis of diverse data sets can reliably address the taxonomy of morphologically challenging plant taxa.

The tomato genome: implications for plant breeding, genomics and evolution

The genome sequence of tomato (Solanum lycopersicum), one of the most important vegetable crops, has recently been decoded. We address implications of the tomato genome for plant breeding, genomics and evolutionary studies, and its potential to fuel future crop biology research.

Physiological adaptation along environmental gradients and replicated hybrid zone structure in swordtails (Teleostei: Xiphophorus)

Local adaptation is often invoked to explain hybrid zone structure, but empirical evidence of this is generally rare. Hybrid zones between two poeciliid fishes, Xiphophorus birchmanni and X. malinche, occur in multiple tributaries with independent replication of upstream-to-downstream gradients in morphology and allele frequencies. Ecological niche modelling revealed that temperature is a central predictive factor in the spatial distribution of pure parental species and their hybrids and explains spatial and temporal variation in the frequency of neutral genetic markers in hybrid populations. Among populations of parentals and hybrids, both thermal tolerance and heat-shock protein expression vary strongly, indicating that spatial and temporal structure is likely driven by adaptation to local thermal environments. Therefore, hybrid zone structure is strongly influenced by interspecific differences in physiological mechanisms for coping with the thermal environment.

Balancing selection at the tomato RCR3 Guardee gene family maintains variation in strength of pathogen defense

Coevolution between hosts and pathogens is thought to occur between interacting molecules of both species. This results in the maintenance of genetic diversity at pathogen antigens (or so-called effectors) and host resistance genes such as the major histocompatibility complex (MHC) in mammals or resistance (R) genes in plants. In plant-pathogen interactions, the current paradigm posits that a specific defense response is activated upon recognition of pathogen effectors via interaction with their corresponding R proteins. According to the "Guard-Hypothesis," R proteins (the "guards") can sense modification of target molecules in the host (the "guardees") by pathogen effectors and subsequently trigger the defense response. Multiple studies have reported high genetic diversity at R genes maintained by balancing selection. In contrast, little is known about the evolutionary mechanisms shaping the guardee, which may be subject to contrasting evolutionary forces. Here we show that the evolution of the guardee RCR3 is characterized by gene duplication, frequent gene conversion, and balancing selection in the wild tomato species Solanum peruvianum. Investigating the functional characteristics of 54 natural variants through in vitro and in planta assays, we detected differences in recognition of the pathogen effector through interaction with the guardee, as well as substantial variation in the strength of the defense response. This variation is maintained by balancing selection at each copy of the RCR3 gene. Our analyses pinpoint three amino acid polymorphisms with key functional consequences for the coevolution between the guardee (RCR3) and its guard (Cf-2). We conclude that, in addition to coevolution at the "guardee-effector" interface for pathogen recognition, natural selection acts on the "guard-guardee" interface. Guardee evolution may be governed by a counterbalance between improved activation in the presence and prevention of auto-immune responses in the absence of the corresponding pathogen.

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.

Little ecological divergence associated with speciation in two African rain forest tree genera

BACKGROUND

The tropical rain forests (TRF) of Africa are the second largest block of this biome after the Amazon and exhibit high levels of plant endemism and diversity. Two main hypotheses have been advanced to explain speciation processes that have led to this high level of biodiversity: allopatric speciation linked to geographic isolation and ecological speciation linked to ecological gradients. Both these hypotheses rely on ecology: in the former conservation of ecological niches through time is implied, while in the latter adaptation via selection to alternative ecological niches would be a prerequisite. Here, we investigate the role of ecology in explaining present day species diversity in African TRF using a species level phylogeny and ecological niche modeling of two predominantly restricted TRF tree genera, Isolona and Monodora (Annonaceae). Both these genera, with 20 and 14 species, respectively, are widely distributed in African TRFs, with a few species occurring in slightly less humid regions such as in East Africa.

RESULTS

A total of 11 sister species pairs were identified most of them occurring in allopatry or with little geographical overlap. Our results provide a mixed answer on the role of ecology in speciation. Although no sister species have identical niches, just under half of the tests suggest that sister species do have more similar niches than expected by chance. PCA analyses also support little ecological differences between sister species. Most speciation events within both genera predate the Pleistocene, occurring during the Late Miocene and Pliocene periods.

CONCLUSIONS

Ecology is almost always involved in speciation, however, it would seem to have had a little role in species generation within Isolona and Monodora at the scale analyzed here. This is consistent with the geographical speciation model for TRF diversification. These results contrast to other studies for non-TRF plant species where ecological speciation was found to be an important factor of diversification. The Pliocene period appears to be a vital time in the generation of African TRF diversity, whereas Pleistocene climatic fluctuations have had a smaller role on speciation than previously thought.Ecological niche modeling, species level phylogeny, ecological speciation, African tropics, Isolona, Monodora, Annonaceae.

Inference of seed bank parameters in two wild tomato species using ecological and genetic data

Seed and egg dormancy is a prevalent life-history trait in plants and invertebrates whose storage effect buffers against environmental variability, modulates species extinction in fragmented habitats, and increases genetic variation. Experimental evidence for reliable differences in dormancy over evolutionary scales (e.g., differences in seed banks between sister species) is scarce because complex ecological experiments in the field are needed to measure them. To cope with these difficulties, we developed an approximate Bayesian computation (ABC) framework that integrates ecological information on population census sizes in the priors of the parameters, along with a coalescent model accounting simultaneously for seed banks and spatial genetic structuring of populations. We collected SNP data at seven nuclear loci (over 300 SNPs) using a combination of three spatial sampling schemes: population, pooled, and species-wide samples. We provide evidence for the existence of a seed bank in two wild tomato species (Solanum chilense and Solanum peruvianum) found in western South America. Although accounting for uncertainties in ecological data, we infer for each species (i) the past demography and (ii) ecological parameters, such as the germination rate, migration rates, and minimum number of demes in the metapopulation. The inferred difference in germination rate between the two species may reflect divergent seed dormancy adaptations, in agreement with previous population genetic analyses and the ecology of these two sister species: Seeds spend, on average, a shorter time in the soil in the specialist species (S. chilense) than in the generalist species (S. peruvianum).

Distinct niche divergence characterizes the homoploid hybrid speciation of Pinus densata on the Tibetan plateau

Ecological divergence and selection for novel adaptations to new habitats have been theoretically proposed to play important roles in promoting homoploid hybrid speciation (HHS). The successful establishment of Pinus densata on the Tibetan Plateau is one of the few known examples of HHS. In this study, we carried out extensive field expeditions to obtain representative coverage of occurrence sites of P. densata and its two putative parents. We then applied a series of geographic information system-based analyses to define the patterns of environmental variation within and among the three pine species, to remove potentially confounding effects of spatial autocorrelation in the environmental data due to allopatric ranges, and to build species distribution models. All results consistently indicated that the ecological preferences of P. densata and its parental species have diverged, and they identified candidate ecological factors associated with habitat-specific adaptation. Projections from niche modeling indicated that P. densata could extend across a vast range along the parallel valley systems of the southeastern Tibetan Plateau. Our findings provide evidence of a distinct niche shift in P. densata and support the hypothesis that local adaptation and geographic isolation help maintain and reinforce between-species differences and reproductive isolation in the species complex.

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.

Ecological niche modelling and nDNA sequencing support a new, morphologically cryptic beetle species unveiled by DNA barcoding

BACKGROUND

DNA sequencing techniques used to estimate biodiversity, such as DNA barcoding, may reveal cryptic species. However, disagreements between barcoding and morphological data have already led to controversy. Species delimitation should therefore not be based on mtDNA alone. Here, we explore the use of nDNA and bioclimatic modelling in a new species of aquatic beetle revealed by mtDNA sequence data.

METHODOLOGY/PRINCIPAL FINDINGS

The aquatic beetle fauna of Australia is characterised by high degrees of endemism, including local radiations such as the genus Antiporus. Antiporus femoralis was previously considered to exist in two disjunct, but morphologically indistinguishable populations in south-western and south-eastern Australia. We constructed a phylogeny of Antiporus and detected a deep split between these populations. Diagnostic characters from the highly variable nuclear protein encoding arginine kinase gene confirmed the presence of two isolated populations. We then used ecological niche modelling to examine the climatic niche characteristics of the two populations. All results support the status of the two populations as distinct species. We describe the south-western species as Antiporus occidentalis sp.n.

CONCLUSION/SIGNIFICANCE

In addition to nDNA sequence data and extended use of mitochondrial sequences, ecological niche modelling has great potential for delineating morphologically cryptic species.

Fitness effects of derived deleterious mutations in four closely related wild tomato species with spatial structure

A key issue in evolutionary biology is an improved understanding of the genetic mechanisms by which species adapt to various environments. Using DNA sequence data, it is possible to quantify the number of adaptive and deleterious mutations, and the distribution of fitness effects of new mutations (its mean and variance) by simultaneously taking into account the demography of a given species. We investigated how selection functions at eight housekeeping genes of four closely related, outcrossing species of wild tomatoes that are native to diverse environments in western South America (Solanum arcanum, S. chilense, S. habrochaites and S. peruvianum). We found little evidence for adaptive mutations but pervasive evidence for strong purifying selection in coding regions of the four species. In contrast, the strength of purifying selection seems to vary among the four species in non-coding (NC) regions (introns). Using F(ST)-based measures of fixation in subdivided populations, we suggest that weak purifying selection has affected the NC regions of S. habrochaites, S. chilense and S. peruvianum. In contrast, NC regions in S. arcanum show a distribution of fitness effects with mutations being either nearly neutral or very strongly deleterious. These results suggest that closely related species with similar genetic backgrounds but experiencing contrasting environments differ in the variance of deleterious fitness effects.

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.