Ecological and geographic modes of species divergence in wild tomatoes

Niche comparisons reveal significant divergence despite narrow endemism in Leavenworthia, a genus of rare plants

BACKGROUND AND AIMS

Quantifying niche similarity among closely related species offers myriad insights into evolutionary history and ecology. In this study, our aim was to explore the interplay of geographical and niche space for rare, endemic plant species and to determine whether endemic habitats were environmentally similar or unique.

METHODS

We characterized the niche of all Leavenworthia species, a genus of rare plants endemic to rocky glades in the eastern USA, using WorldClim data, surface geology, elevation and slope. We calculated the area of range overlap and estimated niche similarity between pairs of species in their total occupied niche space and the subset of niche space shared by both species. We used linear discriminant analyses to determine which niche dimensions differed the most between species. We used niche dimensions with consistently high discriminatory power to perform a random forest classification analysis and principal component analysis. Using a linear model, we related geographical distance to distance in niche space.

KEY RESULTS

Most species comparisons concluded that species niches had diverged, with niche similarity increasing linearly with range overlap. Temperature variation, precipitation amount and seasonality, and surface geology were the most divergent niche dimensions among all species comparisons. Geographical distance explained 42 % of the variation in niche space distance. Sites that were closer in niche space than expected were oriented east-west owing to the strong correlation between latitude and scores on the first principal component.

CONCLUSIONS

Despite being endemic seemingly to very similar habitat, niche similarity is low among Leavenworthia species. Low niche similarity, combined with low geographical overlap, suggests that this lineage of rare plants potentially diversified in isolation but across a very small geographical area. The correlation between geographical space and niche space has received considerable attention, but our results suggest that geographical distance is a weak predictor of distance in niche space.

Evolution of gene networks underlying adaptation to drought stress in the wild tomato Solanum chilense

Drought stress is a key limitation for plant growth and colonization of arid habitats. We study the evolution of gene expression response to drought stress in a wild tomato, Solanum chilense, naturally occurring in dry habitats in South America. We conduct a transcriptome analysis under standard and drought experimental conditions to identify drought-responsive gene networks and estimate the age of the involved genes. We identify two main regulatory networks corresponding to two typical drought-responsive strategies: cell cycle and fundamental metabolic processes. The metabolic network exhibits a more recent evolutionary origin and a more variable transcriptome response than the cell cycle network (with ancestral origin and higher conservation of the transcriptional response). We also integrate population genomics analyses to reveal positive selection signals acting at the genes of both networks, revealing that genes exhibiting selective sweeps of older age also exhibit greater connectivity in the networks. These findings suggest that adaptive changes first occur at core genes of drought response networks, driving significant network re-wiring, which likely underpins species divergence and further spread into drier habitats. Combining transcriptomics and population genomics approaches, we decipher the timing of gene network evolution for drought stress response in arid habitats.

Prediction of suitable regions of wild tomato provides insights on domesticated tomato cultivation in China

Climate change is one of the biggest challenges to the world at present. Tomato is also suffered from devastating yield loss due to climate change. The domesticated tomato (Solanum lycopersicum) is presumed to be originated from the wild tomato (S. pimpinellifolium). In this study, we compared the climate data of S. pimpinellifollium with the domesticated tomato, predicted the suitable regions of S. pimpinellifollium in China using MaxEnt model and assessed their tolerance to drought stress. We found that the predicted suitable regions of wild tomato are highly consistent with the current cultivated regions of domesticated tomato, suggesting that the habitat demand of domesticated tomato descended largely from its ancestor, hence the habitat information of wild tomato could provide a reference for tomato cultivation. We further predicted suitable regions of wild tomato in the future in China. Finally, we found that while average drought tolerance between wild and domesticated tomato accessions shows no difference, tolerance levels among wild tomato accessions exhibit higher variation, which could be used for future breeding to improve drought resistance. To summarize, our study shows that suitable regions of wild tomato provide insights into domesticated tomato cultivation in China.

A chromosome-level genome assembly of Solanum chilense, a tomato wild relative associated with resistance to salinity and drought

Introduction

Solanum chilense is a wild relative of tomato reported to exhibit resistance to biotic and abiotic stresses. There is potential to improve tomato cultivars via breeding with wild relatives, a process greatly accelerated by suitable genomic and genetic resources.

Methods

In this study we generated a high-quality, chromosome-level, de novo assembly for the S. chilense accession LA1972 using a hybrid assembly strategy with ~180 Gbp of Illumina short reads and ~50 Gbp long PacBio reads. Further scaffolding was performed using Bionano optical maps and 10x Chromium reads.

Results

The resulting sequences were arranged into 12 pseudomolecules using Hi-C sequencing. This resulted in a 901 Mbp assembly, with a completeness of 95%, as determined by Benchmarking with Universal Single-Copy Orthologs (BUSCO). Sequencing of RNA from multiple tissues resulting in ~219 Gbp of reads was used to annotate the genome assembly with an RNA-Seq guided gene prediction, and for a de novo transcriptome assembly. This chromosome-level, high-quality reference genome for S. chilense accession LA1972 will support future breeding efforts for more sustainable tomato production.

Discussion

Gene sequences related to drought and salt resistance were compared between S. chilense and S. lycopersicum to identify amino acid variations with high potential for functional impact. These variants were subsequently analysed in 84 resequenced tomato lines across 12 different related species to explore the variant distributions. We identified a set of 7 putative impactful amino acid variants some of which may also impact on fruit development for example the ethylene-responsive transcription factor WIN1 and ethylene-insensitive protein 2. These variants could be tested for their ability to confer functional phenotypes to cultivars that have lost these variants.

Fluctuation of ecological niches and geographic range shifts along chile pepper's domestication gradient

Domestication is an ongoing well-described process. However, while many have studied the changes domestication causes in plant genetics, few have explored its impact on the portion of the geographic landscape in which the plants exist. Therefore, the goal of this study was to understand how the process of domestication changed the geographic space suitable for chile pepper (Capsicum annuum) in its center of origin (domestication). C. annuum is a major crop species globally whose center of domestication, Mexico, has been well-studied. It provides a unique opportunity to explore the degree to which ranges of different domestication classes diverged and how these ranges might be altered by climate change. To this end, we created ecological niche models for four domestication classes (wild, semiwild, landrace, modern cultivar) based on present climate and future climate scenarios for 2050, 2070, and 2090. Considering present environment, we found substantial overlap in the geographic niches of all the domestication classes. Yet, environmental and geographic aspects of the current ranges did vary among classes. Wild and commercial varieties could grow in desert conditions, while landraces could not. With projections into the future, habitat was lost asymmetrically, with wild, semiwild, and landraces at greater risk of territorial declines than modern cultivars. Further, we identified areas where future suitability overlap between landraces and wilds is expected to be lost. While range expansion is widely associated with domestication, we found little support of a constant niche expansion (either in environmental or geographical space) throughout the domestication gradient in chile peppers in Mexico. Instead, particular domestication transitions resulted in loss, followed by capturing or recapturing environmental or geographic space. The differences in environmental characterization among domestication gradient classes and their future potential range shifts increase the need for conservation efforts to preserve landraces and semiwild genotypes.

Variation in ripe fruit volatiles across the tomato clade: An evolutionary framework for studying fruit scent diversity in a crop wild relative

PREMISE

The scents of volatile organic compounds (VOCs) are an important component of ripe fleshy fruit attractiveness, yet their variation across closely related wild species is poorly understood. Phylogenetic patterns in these compounds and their biosynthetic pathways offer insight into the evolutionary drivers of fruit diversity, including whether scent can communicate an honest signal of nutrient content to animal dispersers. We assessed ripe fruit VOC content across the tomato clade (Solanum sect. Lycopersicon), with implications for crop improvement since these compounds are key components of tomato flavor.

METHODS

We analyzed ripe fruit volatiles from 13 species of wild tomato grown in a common garden. Interspecific variations in 66 compounds and their biosynthetic pathways were assessed in 32 accessions, and an accession-level phylogeny was constructed to account for relatedness.

RESULTS

Wild tomato species can be differentiated by their VOCs, with Solanum pennellii notably distinct. Phylogenetic conservatism exists to a limited extent. Major cladewide patterns corresponded to divergence of the five brightly colored-fruited species from the nine green-fruited species, particularly for nitrogen-containing compounds (higher in colored-fruited) and esters (higher in green-fruited), the latter appearing to signal a sugar reward.

CONCLUSIONS

We established a framework for fruit scent evolution studies in a crop wild relative system, showing that each species in the tomato clade has a unique VOC profile. Differences between color groups align with fruit syndromes that could be driven by selection from frugivores. The evolution of colored fruits was accompanied by changes in biosynthetic pathways for esters and nitrogen-containing compounds, volatiles important to tomato flavor.

Geographic Variation in Genomic Signals of Admixture Between Two Closely Related European Sepsid Fly Species

The extent of interspecific gene flow and its consequences for the initiation, maintenance, and breakdown of species barriers in natural systems remain poorly understood. Interspecific gene flow by hybridization may weaken adaptive divergence, but can be overcome by selection against hybrids, which may ultimately promote reinforcement. An informative step towards understanding the role of gene flow during speciation is to describe patterns of past gene flow among extant species. We investigate signals of admixture between allopatric and sympatric populations of the two closely related European dung fly species Sepsis cynipsea and S. neocynipsea (Diptera: Sepsidae). Based on microsatellite genotypes, we first inferred a baseline demographic history using Approximate Bayesian Computation. We then used genomic data from pooled DNA of natural and laboratory populations to test for past interspecific gene flow based on allelic configurations discordant with the inferred population tree (ABBA-BABA test with D-statistic). Comparing the detected signals of gene flow with the contemporary geographic relationship among interspecific pairs of populations (sympatric vs. allopatric), we made two contrasting observations. At one site in the French Cevennes, we detected an excess of past interspecific gene flow, while at two sites in Switzerland we observed lower signals of past microsatellite genotypes gene flow among populations in sympatry compared to allopatric populations. These results suggest that the species boundaries between these two species depend on the past and/or present eco-geographic context in Europe, which indicates that there is no uniform link between contemporary geographic proximity and past interspecific gene flow in natural populations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11692-023-09612-5.

Understanding the Geographic Patterns of Closely-Related Species of Paspalum (Poaceae) Using Distribution Modelling and Seed Germination Traits

The sexual species of the Dilatata complex (Paspalum dasypleurum, P. flavescens, P. plurinerve, P. vacarianum, and P. urvillei) are closely related phylogenetically and show allopatric distributions, except P. urvillei. These species show microhabitat similarities and differences in germination traits. We integrated species distribution models (SDMs) and seed germination assays to determine whether germination divergences explain their biogeographic pattern. We trained SDMs in South America using species' presence-absence data and environmental variables. Additionally, populations sampled from highly favourable areas in the SDMs of these species were grown together, and their seeds germinated at different temperatures and dormancy-breaking conditions. Differences among species in seed dormancy and germination niche breadth were tested, and linear regressions between seed dormancy and climatic variables were explored. SDMs correctly classified both the observed presences and absences. Spatial factors and anthropogenic activities were the main factors explaining these distributions. Both SDMs and germination analyses confirmed that the niche of P. urvillei was broader than the other species which showed restricted distributions, narrower germination niches, and high correlations between seed dormancy and precipitation regimes. Both approaches provided evidence about the generalist-specialist status of each species. Divergences in seed dormancy between the specialist species could explain these allopatric distributions.

Selective sweeps linked to the colonization of novel habitats and climatic changes in a wild tomato species

Positive selection is the driving force underpinning local adaptation and leaves footprints of selective sweeps on the underlying major genes. Quantifying the timing of selection and revealing the genetic bases of adaptation in plant species occurring in steep and varying environmental gradients are crucial to predict a species' ability to colonize new niches. We use whole-genome sequence data from six populations across three different habitats of the wild tomato species Solanum chilense to infer the past demographic history and search for genes under strong positive selection. We then correlate current and past climatic projections with the demographic history, allele frequencies, the age of selection events and distribution shifts. Several selective sweeps occur at regulatory networks involved in root-hair development in low altitude and response to photoperiod and vernalization in high-altitude populations. These sweeps appear to occur in a concerted fashion in a given regulatory gene network at particular periods of substantial climatic change. Using a unique combination of genome scans and modelling of past climatic data, we quantify the timing of selection at genes likely underpinning local adaptation to semiarid habitats.

Increased ACS Enzyme Dosage Causes Initiation of Climacteric Ethylene Production in Tomato

Fruits of wild tomato species show different ethylene-dependent ripening characteristics, such as variations in fruit color and whether they exhibit a climacteric or nonclimacteric ripening transition. 1-Aminocyclopropane-1-carboxylic acid (ACC) synthase (ACS) and ACC oxidase (ACO) are key enzymes in the ethylene biosynthetic pathway encoded by multigene families. Gene duplication is a primary driver of plant diversification and angiosperm evolution. Here, interspecific variations in the molecular regulation of ethylene biosynthesis and perception during fruit ripening in domesticated and wild tomatoes were investigated. Results showed that the activated ACS genes were increased in number in red-ripe tomato fruits than in green-ripe tomato fruits; therefore, elevated dosage of ACS enzyme promoted ripening ethylene production. Results showed that the expression of three ACS isogenes ACS1A, ACS2, and ACS4, which are involved in autocatalytic ethylene production, was higher in red-ripe tomato fruits than in green-ripe tomato fruits. Elevated ACS enzyme dosage promoted ethylene production, which corresponded to the climacteric response of red-ripe tomato fruits. The data suggest that autoinhibitory ethylene production is common to all tomato species, while autocatalytic ethylene production is specific to red-ripe species. The essential regulators Non-ripening (NOR) and Ripening-Inhibitor (RIN) have experienced gene activation and overlapped with increasing ACS enzyme dosage. These complex levels of transcript regulation link higher ethylene production with spatiotemporal modulation of gene expression in red-ripe tomato species. Taken together, this study shows that bursts in ethylene production that accompany fruit color changes in red-ripe tomatoes are likely to be an evolutionary adaptation for seed dispersal.

Auxin-driven ecophysiological diversification of leaves in domesticated tomato

Heterobaric leaves have bundle sheath extensions (BSEs) that compartmentalize the parenchyma, whereas homobaric leaves do not. The presence of BSEs affects leaf hydraulics and photosynthetic rate. The tomato (Solanum lycopersicum) obscuravenosa (obv) mutant lacks BSEs. Here, we identify the obv gene and the causative mutation, a nonsynonymous amino acid change that disrupts a C2H2 zinc finger motif in a putative transcription factor. This mutation exists as a polymorphism in the natural range of wild tomatoes but has increased in frequency in domesticated tomatoes, suggesting that the latter diversified into heterobaric and homobaric leaf types. The obv mutant displays reduced vein density, leaf hydraulic conductance and photosynthetic assimilation rate. We show that these and other pleiotropic effects on plant development, including changes in leaf insertion angle, leaf margin serration, minor vein density, and fruit shape, are controlled by OBV via changes in auxin signaling. Loss of function of the transcriptional regulator AUXIN RESPONSE FACTOR 4 (ARF4) also results in defective BSE development, revealing an additional component of a genetic module controlling aspects of leaf development important for ecological adaptation and subject to breeding selection.

Genetic diversity and demography of Bufo japonicus and B. torrenticola (Amphibia: Anura: Bufonidae) influenced by the Quaternary climate

The Quaternary climate affected the present species richness and geographic distribution patterns of amphibians by limiting their activities during the glacial period. The present study examined the phylogenetic relationships of Japanese toads (Bufo japonicus and B. torrenticola) and the demography of each lineage from the past to the present based on mitochondrial sequences and ecological niche models. Japanese toads are a monophyletic group with two main clades (clades A and B). Clade A represents B. j. formosus, including three clades (clades A1, A2, and A3). Clade B contains three clades, two of which corresponded to B. j. japonicus (clades B1 and B2) and the other to B. torrenticola. Clade B2 and B. torrenticola made a sister group, and, thus, B. j. japonicus is paraphyletic. Clades A and B diverged in the late Miocene 5.7 million years ago (Mya) during the period when the Japanese archipelago was constructed. The earliest divergence between the three clades of clade A was estimated at 1.8 Mya. Clades A1 and A2 may have diverged at 0.8 Mya, resulting from the isolation in the multiple different refugia; however, the effects of the glacial climate on the divergence events of clade A3 are unclear. Divergences within clade B occurred from the late Pliocene to the early Pleistocene (3.2-2.2 Mya). Niche similarity between the parapatric clade in clade B (clades B1 and B2) indicated their allopatric divergence. It was suggested that niche segregation between B. japonicus and B. torrenticola contributed to a rapid adaptation of B. torrenticola for lotic breeding. All clade of Japanese toads retreated to each refugium at a low elevation in the glacial period, and effective population sizes increased to construct the current populations after the Last Glacial Maximum. Furthermore, we highlight the areas of climate stability from the last glacial maximum to the present that have served as the refugia of Japanese toads and, thus, affected their present distribution patterns.

Adaptive Divergence without Distinct Species Relationships Indicate Early Stage Ecological Speciation in Species of the Rhododendronpseudochrysanthum Complex Endemic to Taiwan

The testing association of environmental variables with genetic and epigenetic variation could be crucial to deciphering the effects of environmental factors playing roles as selective drivers in ecological speciation. Although ecological speciation may occur in closely related species, species boundaries may not be established over a short evolutionary timescale. Here, we investigated the genetic and epigenetic variations using amplified fragment length polymorphism (AFLP) and methylation-sensitive amplification polymorphism (MSAP), respectively, and tested their associations with environmental variables in populations of four closely related species in the R. pseudochrysanthum complex. No distinctive species relationships were found using genetic clustering analyses, neighbor-joining tree, and neighbor-net tree based on the total AFLP variation, which is suggestive of the incomplete lineage sorting of ancestral variation. Nonetheless, strong isolation-by-environment and adaptive divergence were revealed, despite the significant isolation-by-distance. Annual mean temperature, elevation, normalized difference vegetation index, and annual total potential evapotranspiration were found to be the most important environmental variables explaining outlier genetic and epigenetic variations. Our results suggest that the four closely related species of the R. pseudochrysanthum complex share the polymorphism of their ancestor, but reproductive isolation due to ecological speciation can occur if local environmental divergence persists over time.

The Halophyte Species Solanum chilense Dun. Maintains Its Reproduction despite Sodium Accumulation in Its Floral Organs

Salinity is a growing global concern that affects the yield of crop species, including tomato (Solanum lycopersicum). Its wild relative Solanum chilense was reported to have halophyte properties. We compared salt resistance of both species during the reproductive phase, with a special focus on sodium localization in the flowers. Plants were exposed to NaCl from the seedling stage. Salinity decreased the number of inflorescences in both species but the number of flowers per inflorescence and sepal length only in S. lycopersicum. External salt supply decreased the stamen length in S. chilense, and it was associated with a decrease in pollen production and an increase in pollen viability. Although the fruit set was not affected by salinity, fruit weight and size decreased in S. lycopersicum. Concentrations and localization of Na, K, Mg, and Ca differed in reproductive structures of both species. Inflorescences and fruits of S. chilense accumulated more Na than S. lycopersicum. Sodium was mainly located in male floral organs of S. chilense but in non-reproductive floral organs in S. lycopersicum. The expression of Na transporter genes differed in flowers of both species. Overall, our results indicated that S. chilense was more salt-resistant than S. lycopersicum during the reproductive phase and that differences could be partly related to dissimilarities in element distribution and transport in flowers.

Edaphoclimatic Descriptors of Wild Tomato Species (Solanum Sect. Lycopersicon) and Closely Related Species (Solanum Sect. Juglandifolia and Sect. Lycopersicoides) in South America

Wild species related to cultivated tomato are essential genetic resources in breeding programs focused on food security to face future challenges. The ecogeographic analysis allows identifying the species adaptive ranges and most relevant environmental variables explaining their patterns of actual distribution. The objective of this research was to identify the diversity, ecological descriptors, and statistical relationship of 35 edaphoclimatic variables (20 climatic, 1 geographic and 14 edaphic variables) from 4,649 accessions of 12 wild tomato species and 4 closely related species classified in Solanum sect. Lycopersicon and clustered into four phylogenetic groups, namely “Lycopersicon group” (S. pimpinellifolium, S. cheesmaniae, and S. galapagense), “Arcanum group” (S. arcanum, S. chmielewskii, and S. neorickii), “Eriopersicon group” (S. habrochaites, S. huaylasense, S. corneliomulleri, S. peruvianum, and S. chilense), “Neolycopersicon group” (S. pennellii); and two phylogenetically related groups in Solanum sect. Juglandifolia (S. juglandifolium and S. ochranthum), and section Lycopersicoides (S. lycopersicoides and S. sitiens). The relationship between the climate and edaphic variables were determined by the canonical correlation analysis, reaching 89.2% of variation with the first three canonical correlations. The most significant climatic variables were related to humidity (annual evapotranspiration, annual precipitation, and precipitation of driest month) and physicochemical soil characteristics (bulk density, pH, and base saturation percentage). In all groups, ecological descriptors and diversity patterns were consistent with previous reports. Regarding edaphoclimatic diversity, 12 climate types and 17 soil units were identified among all species. This approach has promissory applications for biodiversity conservation and uses valuable genetic resources related to a leading crop.

A Cysteine-Rich Protein, SpDIR1L, Implicated in S-RNase-Independent Pollen Rejection in the Tomato (Solanum Section Lycopersicon) Clade

Tomato clade species (Solanum sect. Lycopersicon) display multiple interspecific reproductive barriers (IRBs). Some IRBs conform to the SI x SC rule, which describes unilateral incompatibility (UI) where pollen from SC species is rejected on SI species’ pistils, but reciprocal pollinations are successful. However, SC x SC UI also exists, offering opportunities to identify factors that contribute to S-RNase-independent IRBs. For instance, SC Solanum pennellii LA0716 pistils only permit SC Solanum lycopersicum pollen tubes to penetrate to the top third of the pistil, while S. pennellii pollen penetrates to S. lycopersicum ovaries. We identified candidate S. pennellii LA0716 pistil barrier genes based on expression profiles and published results. CRISPR/Cas9 mutants were created in eight candidate genes, and mutants were assessed for changes in S. lycopersicum pollen tube growth. Mutants in a gene designated Defective in Induced Resistance 1-like (SpDIR1L), which encodes a small cysteine-rich protein, permitted S. lycopersicum pollen tubes to grow to the bottom third of the style. We show that SpDIR1L protein accumulation correlates with IRB strength and that species with weak or no IRBs toward S. lycopersicum pollen share a 150 bp deletion in the upstream region of SpDIR1L. These results suggest that SpDIR1L contributes to an S-RNase-independent IRB.

Genetic structure and ecological niche segregation of Indian gray mongoose (Urva edwardsii) in Iran

Combining genetic data with ecological niche models is an effective approach for exploring climatic and nonclimatic environmental variables affecting spatial patterns of intraspecific genetic variation. Here, we adopted this combined approach to evaluate genetic structure and ecological niche of the Indian gray mongoose (Urva edwardsii) in Iran, as the most western part of the species range. Using mtDNA, we confirmed the presence of two highly differentiated clades. Then, we incorporated ensemble of small models (ESMs) using climatic and nonclimatic variables with genetic data to assess whether genetic differentiation among clades was coupled with their ecological niche. Climate niche divergence was also examined based on a principal component analysis on climatic factors only. The relative habitat suitability values predicted by the ESMs for both clades revealed their niche separation. Between-clade climate only niche comparison revealed that climate space occupied by clades is similar to some extent, but the niches that they utilize differ between the distribution ranges of clades. We found that in the absence of evidence for recent genetic exchanges, distribution models suggest the species occurs in different niches and that there are apparent areas of disconnection across the species range. The estimated divergence time between the two Iranian clades (4.9 Mya) coincides with the uplifting of the Zagros Mountains during the Early Pliocene. The Zagros mountain-building event seems to have prevented the distribution of U. edwardsii populations between the western and eastern parts of the mountains as a result of vicariance events. Our findings indicated that the two U. edwardsii genetic clades in Iran can be considered as two conservation units and can be utilized to develop habitat-specific and climate change-integrated management strategies.

Influence of Pleistocene climatic oscillations on the phylogeography and demographic history of endemic vulnerable trees (section Magnolia) of the Tropical Montane Cloud Forest in Mexico

The Tropical Montane Cloud Forest (TMCF) is a highly dynamic ecosystem that has undergone frequent spatial changes in response to the interglacial-glacial cycles of the Pleistocene. These climatic fluctuations between cold and warm cycles have led to species range shifts and contractions-expansions, resulting in complex patterns of genetic structure and lineage divergence in forest tree species. In this study, we sequenced four regions of the chloroplast DNA (trnT-trnL, trnK5-matk, rpl32-trnL, trnS-trnG) for 20 populations and 96 individuals to evaluate the phylogeography, historical demography, and paleodistributions of vulnerable endemic TMCF trees in Mexico: Magnolia pedrazae (north-region), M. schiedeana (central-region), and M. schiedeana population Oaxaca (south-region). Our data recovered 49 haplotypes that showed a significant phylogeographic structure in three regions: north, central, and south. Bayesian Phylogeographic and Ecological Clustering (BPEC) analysis also supported the divergence in three lineages and highlighted the role of environmental factors (temperature and precipitation) in genetic differentiation. Our historical demography analyses revealed demographic expansions predating the Last Interglacial (LIG, ~125,000 years ago), while Approximate Bayesian Computation (ABC) simulations equally supported two contrasting demographic scenarios. The BPEC and haplotype network analyses suggested that ancestral haplotypes were geographically found in central Veracruz. Our paleodistributions modeling showed evidence of range shifts and expansions-contractions from the LIG to the present, which suggested the complex evolutionary dynamics associated to the climatic oscillations of the Pleistocene. Habitat management of remnant forest fragments where large and genetically diverse populations occur in the three TMCF regions analyzed would be key for the conservation of these magnolia populations.

Coupling coordination relationship between ecosystem services and water-land resources for the Daguhe River Basin, China

Water and land resource utilization is an important driving force of changes in ecosystem services; therefore, research on multi-parameter coupling systems that consider “ecosystem services, water resources, and land resources” together has key significance for river basins. This study aims to reveal the interaction and mutual influence of ecosystem services and water and land resources in the Daguhe River Basin, China, based on the coupling coordination degree model. The results showed that during the period from 2000 to 2010, the coupling coordination degree values for the years 2000, 2005, and 2010 were 0.6005, 0.7292, and 0.8037. The corresponding coupling coordination classifications were categorized as “primary coordinated development”, “intermediate coordinated development,” and “well-coordinated development”, respectively. These results reflected the fact that the relationship between water and land resource utilization and the environment tends to evolve in the direction of coordinated development (an improvement in one part corresponds to an improvement in another part) with variation in water and land utilization types, and eventually pushes the whole resource, as well as ecological and environmental systems, from low to high levels of coupling coordination degrees as observed in case of the Daguhe River Basin, China. Our research provides an overview of the interaction between ecosystem services and water and land resources in the Daguhe Basin and even in the Shandong Province. With our results, we offer new perspectives on river basin management and for planning future eco-environmental policies (the policy is specifically designed for the ecological environment) by combining water and land resource utilization.

Comparison of Drought and Heat Resistance Strategies among Six Populations of Solanum chilense and Two Cultivars of Solanum lycopersicum

Within the tomato clade, Solanum chilense is considered one of the most promising sources of genes for tomato (S. lycopersicum) selection to biotic and abiotic stresses. In this study, we compared the effects of drought, high temperature, and their combination in two cultivars of S. lycopersicum and six populations of S. chilense, differing in their local habitat. Plants were grown at 21/19 °C or 28/26 °C under well-watered and water-stressed conditions. Plant growth, physiological responses, and expression of stress-responsive genes were investigated. Our results demonstrated strong variability among accessions. Differences in plant growth parameters were even higher among S. chilense populations than between species. The effects of water stress, high temperature, and their combination also differed according to the accession, suggesting differences in stress resistance between species and populations. Overall, water stress affected plants more negatively than temperature from a morpho-physiological point of view, while the expression of stress-responsive genes was more affected by temperature than by water stress. Accessions clustered in two groups regarding resistance to water stress and high temperature. The sensitive group included the S. lycopersicum cultivars and the S. chilense populations LA2931 and LA1930, and the resistant group included the S. chilense populations LA1958, LA2880, LA2765, and LA4107. Our results suggested that resistance traits were not particularly related to the environmental conditions in the natural habitat of the populations. The expression of stress-responsive genes was more stable in resistant accessions than in sensitive ones in response to water stress and high temperature. Altogether, our results suggest that water stress and high temperature resistance in S. chilense did not depend on single traits but on a combination of morphological, physiological, and genetic traits.

Brassica rapa Domestication: Untangling Wild and Feral Forms and Convergence of Crop Morphotypes

The study of domestication contributes to our knowledge of evolution and crop genetic resources. Human selection has shaped wild Brassica rapa into diverse turnip, leafy, and oilseed crops. Despite its worldwide economic importance and potential as a model for understanding diversification under domestication, insights into the number of domestication events and initial crop(s) domesticated in B. rapa have been limited due to a lack of clarity about the wild or feral status of conspecific noncrop relatives. To address this gap and reconstruct the domestication history of B. rapa, we analyzed 68,468 genotyping-by-sequencing-derived single nucleotide polymorphisms for 416 samples in the largest diversity panel of domesticated and weedy B. rapa to date. To further understand the center of origin, we modeled the potential range of wild B. rapa during the mid-Holocene. Our analyses of genetic diversity across B. rapa morphotypes suggest that noncrop samples from the Caucasus, Siberia, and Italy may be truly wild, whereas those occurring in the Americas and much of Europe are feral. Clustering, tree-based analyses, and parameterized demographic inference further indicate that turnips were likely the first crop type domesticated, from which leafy types in East Asia and Europe were selected from distinct lineages. These findings clarify the domestication history and nature of wild crop genetic resources for B. rapa, which provides the first step toward investigating cases of possible parallel selection, the domestication and feralization syndrome, and novel germplasm for Brassica crop improvement.

The impact of climate change on western Plethodon salamanders' distribution

AIM

Given that salamanders have experienced large shifts in their distributions over time, we determined how each species of Plethodon in the Pacific Northwest would respond to climate change. We incorporated several greenhouse scenarios both on a species-by-species basis, and also using phylogenetic groups, with the aim to determine the best course of action in managing land area to conserve diversity in this group.

LOCATION

Pacific Northwest of the United States (northern CA, OR, WA, ID, and MT).

MAJOR TAXA STUDIED

Western Plethodon salamanders.

METHODS

Species distribution models were estimated using MaxEnt for the current time period and for several future climate scenarios using bioclimatic data layers. We used several methods to quantify the change in habitat suitability over time from the models. We explored aspects of the climate layers to determine whether we can expect a concerted response to climate change due to similarity in ecological niche or independent responses that could be harder to manage.

RESULTS

The distribution of western Plethodon salamander species is strongly influenced by precipitation and less so by temperature. Species responses to climate change resulted in both increases and decreases in predicted suitable habitat, though most species ranges do not contract, especially when taken as a phylogenetic group.

MAIN CONCLUSIONS

While some established habitats may become more or less climatically suitable, the overall distribution of species in this group is unlikely to be significantly affected. Clades of Plethodon species are unlikely to be in danger of extirpation despite the possibility that individual species may be threatened as a result of limited distributions. Grouping species into lineages with similar geographic ranges can be a viable method of determining conservation needs. More biotic and dispersal information is needed to determine the true impact that changes in climate will have on the distribution of Plethodon species.

Genome size evolution is associated with climate seasonality and glucosinolates, but not life history, soil nutrients or range size, across a clade of mustards

BACKGROUND AND AIMS

We investigate patterns of evolution of genome size across a morphologically and ecologically diverse clade of Brassicaceae, in relation to ecological and life history traits. While numerous hypotheses have been put forward regarding autecological and environmental factors that could favour small vs. large genomes, a challenge in understanding genome size evolution in plants is that many hypothesized selective agents are intercorrelated.

METHODS

We contribute genome size estimates for 47 species of Streptanthus Nutt. and close relatives, and take advantage of many data collections for this group to assemble data on climate, life history, soil affinity and composition, geographic range and plant secondary chemistry to identify simultaneous correlates of variation in genome size in an evolutionary framework. We assess models of evolution across clades and use phylogenetically informed analyses as well as model selection and information criteria approaches to identify variables that can best explain genome size variation in this clade.

KEY RESULTS

We find differences in genome size and heterogeneity in its rate of evolution across subclades of Streptanthus and close relatives. We show that clade-wide genome size is positively associated with climate seasonality and glucosinolate compounds. Model selection and information criteria approaches identify a best model that includes temperature seasonality and fraction of aliphatic glucosinolates, suggesting a possible role for genome size in climatic adaptation or a role for biotic interactions in shaping the evolution of genome size. We find no evidence supporting hypotheses of life history, range size or soil nutrients as forces shaping genome size in this system.

CONCLUSIONS

Our findings suggest climate seasonality and biotic interactions as potential forces shaping the evolution of genome size and highlight the importance of evaluating multiple factors in the context of phylogeny to understand the effect of possible selective agents on genome size.

Climatic Diversity and Ecological Descriptors of Wild Tomato Species (Solanum sect. Lycopersicon) and Close Related Species (Solanum sect. Juglandifolia y sect. Lycopersicoides) in Latin America

Conservation and sustainable use of species diversity require a description of the environment where they develop. The objectives were to determine ecological descriptors and climatic diversity of areas along the distribution range of 12 species of wild tomatoes (Solanum sect. Lycopersicon) and four wild species of phylogenetically related groups (Solanum sect. Juglandifolia and sect. Lycopersicoides), as well as their ecological similarity in Latin America. With 4228 selected tomato accessions and an environmental information system (EIS) composed of 21 climatic variables, diversity patterns of the distribution areas were identified for each species, as well as ecological descriptors through the use of geographic information systems (GIS). The contribution of climatic variables to the species geographical distribution was identified by principal component analysis (PCA), and similarity in species distribution as a function of the variables identified with cluster analysis (CA). Climatic characteristics and the environmental amplitude of wild tomatoes and related species along their distributional range were satisfactorily determined by ecological descriptors. Eleven climate types were identified, predominantly BSk (arid, steppe, cold), BWh (arid, desert, hot), and Cfb (temperate, no dry season, warm summer). PCA determined 10 most important variables were the most important for the geographical distribution. Six groups of species were identified according to CA and climatic distribution similarity. This approach has shown promissory applications for biodiversity conservation of valuable genetic resources for tomato crop breeding.

Reliable genomic strategies for species classification of plant genetic resources

BACKGROUND

To address the need for easy and reliable species classification in plant genetic resources collections, we assessed the potential of five classifiers (Random Forest, Neighbour-Joining, 1-Nearest Neighbour, a conservative variety of 3-Nearest Neighbours and Naive Bayes) We investigated the effects of the number of accessions per species and misclassification rate on classification success, and validated theirs generic value results with three complete datasets.

RESULTS

We found the conservative variety of 3-Nearest Neighbours to be the most reliable classifier when varying species representation and misclassification rate. Through the analysis of the three complete datasets, this finding showed generic value. Additionally, we present various options for marker selection for classification taks such as these.

CONCLUSIONS

Large-scale genomic data are increasingly being produced for genetic resources collections. These data are useful to address species classification issues regarding crop wild relatives, and improve genebank documentation. Implementation of a classification method that can improve the quality of bad datasets without gold standard training data is considered an innovative and efficient method to improve gene bank documentation.

Estimating hybridization as a consequence of climatic fluctuations: A case study of two geometridae species

In some cases, the phylogenetic analysis based on the mitochondrial DNA (mtDNA) and the nuclear DNA (ncDNA) are discordant. There are three major causes of the discordance within insects, including hybridization, incomplete lineage sorting (ILS) and infection by Wolbachia. In this study, we used a combination of multilocus and coalescent analyses to explore these processes occurred during the evolutionary history of Limbatochlamys rosthorni Rothschild, 1894 and Limbatochlamys pararosthorni Han and Xue, 2005. The ncDNA phylogenetic tree supported two reciprocally monophyletic species, whereas the mtDNA results failed to reveal such a structure and revealed an extensive level of admixture between two species. Because of very low Wolbachia infection rates (<20%), we firstly excluded this reason for the discordance. The fixed nucleotide differences and large genetic distances (1.5-2.5%) at the ncDNA genes suggested that the lineage sorting process between these two species is nearly complete and two species have experienced a prolonged period of independent evolution. Thus, we secondly excluded ILS. Sharing haplotypes, mtDNA gene flow occurring and the transitional samples with morphological features supported hybridization. The distribution contraction during glaciations and postglacial distribution expansion might have facilitated hybridization. Taken together, our study indicates that the current genetic structure of L. rosthorni and L. pararosthorni is the results of contraction and fragmentation into separated refugia during glaciations, followed by postglacial expansion and admixture.

Biosystematic studies on the status of Solanum chilense

PREMISE

Common taxonomic practices, which condition species' descriptions on diagnostic morphological traits, may systematically lump outcrossing species and unduly split selfing species. Specifically, higher effective population sizes and genetic diversity of obligate outcrossers are expected to result less reliable phenotypic diagnoses. Wild tomatoes, members of Solanum sect. Lycopersicum, are commonly used as a source of exotic germplasm for improvement of the cultivated tomato, and are increasingly employed in basic research. Although the section experienced significant early work, which continues presently, the taxonomic status of many wild species has undergone a number of significant revisions and remains uncertain. Species in this section vary in their breeding systems, notably the expression of self-incompatibility, which determines individual propensity for outcrossing METHODS: Here, we examine the taxonomic status of obligately outcrossing Chilean wild tomato (Solanum chilense) using reduced-representation sequencing (RAD-seq), a range of phylogenetic and population genetic analyses, as well as analyses of crossing and morphological data.

RESULTS

Overall, each of our analyses provides a considerable weight of evidence that the Pacific coastal populations and Andean inland populations of the currently described Solanum chilense represent separately evolving populations, and conceal at least one undescribed cryptic species.

CONCLUSIONS

Despite its vast economic importance, Solanum sect. Lycopersicon still exhibits considerable taxonomic instability. A pattern of under-recognition of outcrossing species may be common, not only in tomatoes, but across flowering plants. We discuss the possible causes and implications of this observation, with a focus on macroevolutionary inference.

Subspecies divergence and pronounced phylogenetic incongruence in the East-Asia-endemic shrub Magnolia sieboldii

BACKGROUND AND AIMS

The biogeographic patterns of the East-Asia-endemic shrub Magnolia sieboldii, in which the range of the subsp. sieboldii is interposed with the disjunct distribution of subsp. japonica, implies a complex evolutionary history, involving rapid speciation and hybridization. Here, we aim to reveal the evolutionary and phylogeographic histories of the species with a particular focus on the time of subspecies divergence, the hypothesis of secondary hybridization and the Pleistocene survival of each subspecies, using a combination of genetic analyses and ecological niche modelling.

METHODS

Genetic variation, genetic structures and phylogenetic relationships were elucidated based on nuclear low-copy genes, chloroplast DNA, and nuclear simple sequence repeats (SSRs). A scenario selection analysis and divergence time estimation were performed using coalescent simulation in DIYABC and *BEAST. Ecological niche modelling and a test of niche differentiation were performed using Maxent and ENMTools.

KEY RESULTS

All marker types showed deep, but pronouncedly incongruent, west-east genetic divergences, with the subspecies being delineated only by the nuclear low-copy genes. Phylogenetic tree topologies suggested that ancient hybridization and introgression were likely to have occurred; however, this scenario did not receive significant support in the DIYABC analysis. The subspecies differentiated their niches, but both showed a dependence on high humidity and were predicted to have persisted during the last glacial cycle by maintaining a stable latitudinal distribution via migration to lower altitudes.

CONCLUSIONS

We found a deep genetic divergence and a pronounced phylogenetic incongruence among the two subspecies of M. sieboldii, which may have been driven by major paleogeographic and paleoclimatic events that have occurred since the Neogene in East Asia, including global cooling, climate oscillations and the formation of land bridges. Both subspecies were, however, considered to persist in situ in stable climatic conditions during the late Pleistocene.

Population Genomics of the "Arcanum" Species Group in Wild Tomatoes: Evidence for Separate Origins of Two Self-Compatible Lineages

Given their diverse mating systems and recent divergence, wild tomatoes (Solanum section Lycopersicon) have become an attractive model system to study ecological divergence, the build-up of reproductive barriers, and the causes and consequences of the breakdown of self-incompatibility. Here we report on a lesser-studied group of species known as the "Arcanum" group, comprising the nominal species Solanum arcanum, Solanum chmielewskii, and Solanum neorickii. The latter two taxa are self-compatible but are thought to self-fertilize at different rates, given their distinct manifestations of the morphological "selfing syndrome." Based on experimental crossings and transcriptome sequencing of a total of 39 different genotypes from as many accessions representing each species' geographic range, we provide compelling evidence for deep genealogical divisions within S. arcanum; only the self-incompatible lineage known as "var. marañón" has close genealogical ties to the two self-compatible species. Moreover, there is evidence under multiple inference schemes for different geographic subsets of S. arcanum var. marañón being closest to S. chmielewskii and S. neorickii, respectively. To broadly characterize the population-genomic consequences of these recent mating-system transitions and their associated speciation events, we fit demographic models indicating strong reductions in effective population size, congruent with reduced nucleotide and S-locus diversity in the two independently derived self-compatible species.

Ecological Factors Generally Not Altitude Related Played Main Roles in Driving Potential Adaptive Evolution at Elevational Range Margin Populations of Taiwan Incense Cedar (Calocedrus formosana)

Population diversification can be shaped by a combination of environmental factors as well as geographic isolation interacting with gene flow. We surveyed genetic variation of 243 samples from 12 populations of Calocedrus formosana using amplified fragment length polymorphism (AFLP) and scored a total of 437 AFLP fragments using 11 selective amplification primer pairs. The AFLP variation was used to assess the role of gene flow on the pattern of genetic diversity and to test environments in driving population adaptive evolution. This study found the relatively lower level of genetic diversity and the higher level of population differentiation in C. formosana compared with those estimated in previous studies of conifers including Cunninghamia konishii, Keteleeria davidiana var. formosana, and Taiwania cryptomerioides occurring in Taiwan. BAYESCAN detected 26 F_ST outlier loci that were found to be associated strongly with various environmental variables using multiple univariate logistic regression, latent factor mixed model, and Bayesian logistic regression. We found several environmentally dependent adaptive loci with high frequencies in low- or high-elevation populations, suggesting their involvement in local adaptation. Ecological factors, including relative humidity and sunshine hours, that are generally not altitude related could have been the most important selective drivers for population divergent evolution in C. formosana. The present study provides fundamental information in relation to adaptive evolution and can be useful for assisted migration program of C. formosana in the future conservation of this species.

Toward Unifying Global Hotspots of Wild and Domesticated Biodiversity

Global biodiversity hotspots are areas containing high levels of species richness, endemism and threat. Similarly, regions of agriculturally relevant diversity have been identified where many domesticated plants and animals originated, and co-occurred with their wild ancestors and relatives. The agro-biodiversity in these regions has, likewise, often been considered threatened. Biodiversity and agro-biodiversity hotspots partly overlap, but their geographic intricacies have rarely been investigated together. Here we review the history of these two concepts and explore their geographic relationship by analysing global distribution and human use data for all plants, and for major crops and associated wild relatives. We highlight a geographic continuum between agro-biodiversity hotspots that contain high richness in species that are intensively used and well known by humanity (i.e., major crops and most viewed species on Wikipedia) and biodiversity hotspots encompassing species that are less heavily used and documented (i.e., crop wild relatives and species lacking information on Wikipedia). Our contribution highlights the key considerations needed for further developing a unifying concept of agro-biodiversity hotspots that encompasses multiple facets of diversity (including genetic and phylogenetic) and the linkage with overall biodiversity. This integration will ultimately enhance our understanding of the geography of human-plant interactions and help guide the preservation of nature and its contributions to people.

Regional differences in the abiotic environment contribute to genomic divergence within a wild tomato species

The wild currant tomato Solanum pimpinellifolium inhabits a wide range of abiotic habitats across its native range of Ecuador and Peru. Although it has served as a key genetic resource for the improvement of domestic cultivars, little is known about the genetic basis of traits underlying local adaptation in this species, nor what abiotic variables are most important for driving differentiation. Here we use redundancy analysis (RDA) and other multivariate statistical methods (structural equation modelling [SEM] and generalized dissimilarity modelling [GDM]) to quantify the relationship of genomic variation (6,830 single nucleotide polymorphisms [SNPs]) with climate and geography, among 140 wild accessions. RDA, SEM and GDM each identified environment as explaining more genomic variation than geography, suggesting that local adaptation to heterogeneous abiotic habitats may be an important source of genetic diversity in this species. Environmental factors describing temporal variation in precipitation and evaporative demand explained the most SNP variation among accessions, indicating that these forces may represent key selective agents. Lastly, by studying how SNP-environment associations vary throughout the genome (44,064 SNPs), we mapped the location and investigated the functions of loci putatively contributing to climatic adaptations. Together, our findings indicate an important role for selection imposed by the abiotic environment in driving genomic differentiation between populations.

Assessing biological factors affecting postspeciation introgression

An increasing number of phylogenomic studies have documented a clear “footprint” of postspeciation introgression among closely related species. Nonetheless, systematic genome‐wide studies of factors that determine the likelihood of introgression remain rare. Here, we propose an a priori hypothesis‐testing framework that uses introgression statistics—including a new metric of estimated introgression, D_p—to evaluate general patterns of introgression prevalence and direction across multiple closely related species. We demonstrate this approach using whole genome sequences from 32 lineages in 11 wild tomato species to assess the effect of three factors on introgression—genetic relatedness, geographical proximity, and mating system differences—based on multiple trios within the “ABBA–BABA” test. Our analyses suggest each factor affects the prevalence of introgression, although our power to detect these is limited by the number of comparisons currently available. We find that of 14 species pairs with geographically “proximate” versus “distant” population comparisons, 13 showed evidence of introgression; in 10 of these cases, this was more prevalent between geographically closer populations. We also find modest evidence that introgression declines with increasing genetic divergence between lineages, is more prevalent between lineages that share the same mating system, and—when it does occur between mating systems—tends to involve gene flow from more inbreeding to more outbreeding lineages. Although our analysis indicates that recent postspeciation introgression is frequent in this group—detected in 15 of 17 tested trios—estimated levels of genetic exchange are modest (0.2–2.5% of the genome), so the relative importance of hybridization in shaping the evolutionary trajectories of these species could be limited. Regardless, similar clade‐wide analyses of genomic introgression would be valuable for disentangling the major ecological, reproductive, and historical determinants of postspeciation gene flow, and for assessing the relative contribution of introgression as a source of genetic variation.

Niche modeling reveals life history shifts in birds at La Brea over the last twenty millennia

A species presence at a particular site can change over time, resulting in temporally dynamic species pools. Ecological niche models provide estimates of species presence at different time intervals. The avifauna of La Brea includes approximately 120 species dating to approximately 15,000 years ago. Niche models predicted presence at the Last Glacial Maximum for over 90% of 89 landbird species. This confirms that niche modeling produces sensible range estimates at the Last Glacial Maximum. For 97 currently local species that are as yet undocumented at La Brea over 90% were predicted to occur; absence is due to insufficient study, lack of the ecological niche, transient occurrence or a behavioral ability to avoid entrapment. Our 366 niche models provide a prospective checklist of the landbird fauna of La Brea. The models indicate fluidity in life history strategies and a higher proportion of resident birds at the LGM (88% to 60%). We evaluated a subset of 103 species in breeding and winter periods using two climate models (MIROC-ESM, CCSM4) with a variety of differing parameters, finding differences in 5% of the niche models. Niche breadths in bark-foraging birds changed little between the present and LGM, suggesting that greater species diversity at the LGM was due to greater niche availability rather than contractions of niche breadths (i.e., niche partitioning).

Leaps and bounds: geographical and ecological distance constrained the colonisation of the Afrotemperate by Erica

Background

The coincidence of long distance dispersal (LDD) and biome shift is assumed to be the result of a multifaceted interplay between geographical distance and ecological suitability of source and sink areas. Here, we test the influence of these factors on the dispersal history of the flowering plant genus Erica (Ericaceae) across the Afrotemperate. We quantify similarity of Erica climate niches per biogeographic area using direct observations of species, and test various colonisation scenarios while estimating ancestral areas for the Erica clade using parametric biogeographic model testing.

Results

We infer that the overall dispersal history of Erica across the Afrotemperate is the result of infrequent colonisation limited by geographic proximity and niche similarity. However, the Drakensberg Mountains represent a colonisation sink, rather than acting as a “stepping stone” between more distant and ecologically dissimilar Cape and Tropical African regions. Strikingly, the most dramatic examples of species radiations in Erica were the result of single unique dispersals over longer distances between ecologically dissimilar areas, contradicting the rule of phylogenetic biome conservatism.

Conclusions

These results highlight the roles of geographical and ecological distance in limiting LDD, but also the importance of rare biome shifts, in which a unique dispersal event fuels evolutionary radiation.

Do selfing species have greater niche breadth? Support from ecological niche modeling

We explore the relationship between plant mating system (selfing or outcrossing) and niche breadth to gain new insights into processes that drive species distributions. Using a comparative approach with highly selfing versus highly outcrossing sister species, we test the extent to which: (1) species pairs have evolved significant niche divergence and less niche overlap, (2) selfers have wider niche breadths than outcrossers or vice versa, and (3) niches of selfers and outcrossers are defined by significant differences in environmental variables. We applied predictive ecological niche modeling approaches to estimate and contrast niche divergence, overlap and breadth, and to identify key environmental variables associated with each species' niche for seven sister species with divergent mating systems. Data from 4862 geo-referenced herbarium occurrence records were compiled for 14 species in Collinsia and Tonella (Plantaginaceae) and 19 environmental variables associated with each record. We found sister species display significant niche divergence, though not as a function of divergence time, and overall, selfers have significantly wider niche breadths compared to their outcrossing sisters. Our results suggest that a selfing mating system likely contributes to the greater capacity to reach, reproduce, establish, and adapt to new habitats, which increases niche breadth of selfers.

Testing the Effect of Mountain Ranges as a Physical Barrier to Current Gene Flow and Environmentally Dependent Adaptive Divergence in Cunninghamia konishii (Cupressaceae)

Populations can be genetically isolated by differences in their ecology or environment that hampered efficient migration, or they may be isolated solely by geographic distance. Moreover, mountain ranges across a species’ distribution area might have acted as barriers to gene flow. Genetic variation was quantified using amplified fragment length polymorphism (AFLP) and 13 selective amplification primer combinations used generated a total of 482 fragments. Here, we tested the barrier effects of mountains on gene flow and environmentally dependent local adaptation of Cunninghamia konishii occur in Taiwan. A pattern of genetic isolation by distance was not found and variation partitioning revealed that environment explained a relatively larger proportion of genetic variation than geography. The effect of mountains as barriers to genetic exchange, despite low population differentiation indicating a high rate of gene flow, was found within the distribution range of C. konishii. Twelve AFLP loci were identified as potential selective outliers using genome-scan methods (BAYESCAN and DFDIST) and strongly associated with environmental variables using regression approaches (LFMM, Samβada, and rstanarm) demonstrating adaptive divergence underlying local adaptation. Annual mean temperature, annual precipitation, and slope could be the most important environmental factors causally associated with adaptive genetic variation in C. konishii. The study revealed the existence of physical barriers to current gene flow and environmentally dependent adaptive divergence, and a significant proportion of the rate of gene flow may represent a reflection of demographic history.

Ecological niche modeling and a lack of phylogeographic structure in Vriesea incurvata suggest historically stable areas in the southern Atlantic Forest

PREMISE

Despite the efforts to understand the processes that shape neotropical biodiversity, the complexity of certain biomes, such as the Atlantic Forest (AF), prevents the generalization of patterns. Initially, ecological niche modeling (ENM), with phylogeographic studies, identified past stable areas in the central and northern portions of the AF, while the southern portion was thought to be highly fragmented. Here, we examined the phylogeography, historical patterns, genetic diversity, and population structure of Vriesea incurvata, an endemic species of the southern portion of the AF.

METHODS

We evaluated 149 individuals using two plastid DNA regions (cpDNA) and 13 nuclear microsatellite markers (nuSSRs) to access the historical patterns, genetic diversity, and structure of V. incurvata populations. We also conducted historical demography and ENM analyses.

RESULTS

We found moderate to high genetic diversity and low population structure for both genomes. The cpDNA network revealed high haplotype sharing. The ENM suggested no drastic changes in suitable areas for V. incurvata occurrence, corroborating the finding of no phylogeographic structure.

CONCLUSIONS

Contrary to some studies, our results indicate that the southern AF was a historically stable climate region for V. incurvata occupation after southward colonization by the species. Past climatic changes probably did not cause structuring among its populations.

Ecological niche differentiation between Acanthodactylus micropholis and A. khamirensis (Sauria: Lacertidae) in southern Iran

Acanthodactylusmicropholis Heidari, Rastegar-Pouyani, Rastegar-Pouyani &amp; Rajabizadeh, 2013 and A.khamirensis Blanford, 1874 are genetically and morphologically distinct, but their ecological differentiation has not previously been evaluated. The ecological niche models of these two sister species Acanthodactylus were reconstructed using climate and geographical data. Species distribution modeling for A.micropholis and A.khamirensis was used to make predictions and showed that most parts of southern and southeastern Iran are suitable for the distribution of both species. Habitat suitability was mostly dependent upon minimum temperature of the coldest month and seasonal precipitation for A.micropholis and A.khamirensis, respectively. Niche similarity tests (niche overlap and identity tests) were performed to evaluate species differentiation based on the ecological species criterion. Our results indicate that both species have different ecological niches and are significantly separated from each other. Therefore, our study corroborates previous analyses based on molecular and morphological evidences that suggested that A.micropholis and A.khamirensis were valid species.

Quaternary climatic fluctuations and resulting climatically suitable areas for Eurasian owlets

AIM

The nested pattern in the geographical distribution of three Indian owlets, resulting in a gradient of endemicity, is hypothesized to be an impact of historical climate change. In current time, the Forest Owlet Athene blewitti is endemic to central India, and its range is encompassed within the ranges of the Jungle Owlet Glaucidium radiatum (distributed through South Asia) and Spotted Owlet Athene brama (distributed through Iran, South and Southeast Asia). Another phylogenetically close species, Little Owl Athene noctua, which is largely Palearctic in distribution, is hypothesized to have undergone severe range reduction during the Last Glacial Maximum, showing a postglacial expansion. The present study tests hypotheses on the possible role of Quaternary climatic fluctuations in shaping geographical ranges of owlets.

METHODS

We used primary field observations, open access data, and climatic niche modeling to construct climatic niches of four owlets for four periods, the Last Interglacial (~120-140 Ka), Last Glacial Maximum (~22 Ka), Mid-Holocene (~6 Ka), and Current (1960-1990). We performed climatic niche extent, breadth, and overlap analyses and tested if climatically suitable areas for owlets are nested in a relatively stable climate.

RESULTS

Climatically suitable areas for all owlets examined underwent cycles of expansion and reduction or a gradual expansion or reduction since the Last Interglacial. The Indian owlets show significant climatic niche overlap in the current period. Climatically suitable areas for Little Owl shifted southwards during the Last Glacial Maximum and expanded northwards in the postglaciation period. For each owlet, the modeled climatic niches were nested in climatically stable areas.

MAIN CONCLUSIONS

The study highlights the impact of Quaternary climate change in shaping the present distribution of owlets. This is relevant to the current scenario of climate change and global warming and can help inform conservation strategies, especially for the extremely range-restricted Forest Owlet.

Phylogeography of the red-bellied lizard, Darevskia parvula in Turkey

The groups of red-bellied lizards had a small distribution area in the Pontic zone. The several studies performed on these lizard groups are based on taxonomy and systematics. Although there were several taxonomic or systematic researches on some species of this group, the phylogeographical pattern and species disturbing boundaries of this group is still not clear. In the present study, we aimed to resolve the taxonomic and phylogenetic relationships of the red-bellied lizards in Turkey, based on two combined mitochondrial gene fragments and one protein-coding nuclear gene (rag1). Also, we evaluated ecological niches differentiations among red-bellied lizard groups. The mitochondrial DNA genes were found to be highly polymorphic in this group. One hundred and one variable nucleotide sites were detected on the combined gene sequences. According to phylogenetic trees based on the maximum likelihood (ML) and Bayesian inference (BI), the red-bellied lizards group have three species groups; Darevskia parvula, D. adjarica and unnamed Darevskia sp. (candidate species for Darevskia genus). This situation was supported by high bootstrap and posterior probability values in the trees of mitochondrial DNA gene fragments. However, no genetic variation was detected according to nuclear DNA (rag1) sequence. Because the species groups have no overlaps in terms of their ecological niches, ecological niche modelling (ENM) results revealed differences among the groups of D. parvula, D. adjarica, and unnamed Darevskia sp. Besides, we detected no geographical overlaps among three species groups, since there were geographical isolation zones among the species groups of red-bellied lizard.