Phylogenomic inference of two widespread European leaf miner species complexes suggests mechanisms for sympatric speciation (Lepidoptera: Nepticulidae: Ectoedemia)

Leafmining insects have extraordinarily intimate relationships with their host plants and are therefore prime candidates for potential ecological speciation in sympatry, but how commonly this mode of speciation occurs in any group of life remains debated. Prior research on the pygmy moth (Nepticulidae) genus Ectoedemia using morphology, DNA barcoding, and multi-marker phylogenetic approaches left two unresolved species complexes, each with each four described species. In this study, we thoroughly sampled each complex across a range of host plants and localities across Europe. We used a double digest restriction-site associated DNA (ddRAD) sequencing approach to clarify species boundaries. In the E. rubivora complex, ddRAD data resolved all four species, contrary to morphological and COI data, which supports a potential scenario of host plant-driven speciation where the host plant specialization provides an ecological barrier to hybridization. However, we found no indication of host race formation within the oligophagous E. atricollis (Stainton). In the E. subbimaculella complex, SNP data only partly distinguishes between the parapatric E. subbimaculella (Haworth) and E. heringi (Toll), but with some statistical overlap, suggesting incomplete lineage sorting which may represent early phases of host-based ecological speciation, or admixture following a period of isolation.     

Prominent intraspecific genetic divergence within Anopheles gambiae sibling species triggered by habitat discontinuities across a riverine landscape

The Anopheles gambiae complex of mosquitoes includes malaria vectors at different stages of speciation, whose study enables a better understanding of how adaptation to divergent environmental conditions leads to evolution of reproductive isolation. We investigated the population genetic structure of closely related sympatric taxa that have recently been proposed as separate species (An. coluzzii and An. gambiae), sampled from diverse habitats along the Gambia river in West Africa. We characterized putatively neutral microsatellite loci as well as chromosomal inversion polymorphisms known to be associated with ecological adaptation. The results revealed strong ecologically associated population subdivisions within both species. Microsatellite loci on chromosome‐3L revealed clear differentiation between coastal and inland populations, which in An. coluzzii is reinforced by a unusual inversion polymorphism pattern, supporting the hypothesis of genetic divergence driven by adaptation to the coastal habitat. A strong reduction of gene flow was observed between An. gambiae populations west and east of an extensively rice‐cultivated region apparently colonized exclusively by An. coluzzii. Notably, this ‘intraspecific’ differentiation is higher than that observed between the two species and involves also the centromeric region of chromosome‐X which has previously been considered a marker of speciation within this complex, possibly suggesting that the two populations may be at an advanced stage of differentiation triggered by human‐made habitat fragmentation. These results confirm ongoing ecological speciation within these most important Afro‐tropical malaria vectors and raise new questions on the possible effect of this process in malaria transmission.     

Integrative taxonomy tests possible hybridisation between Central Asian cerambycids (Coleoptera)

Hybridisation can lead to the formation of new evolutionary lineages and some described insect species may in fact be first-generation hybrids of their parental species. Anoplistes forticornis and A. galusoi are two closely related cerambycid taxa endemic to Central Asia that exhibit a peculiar geographical distribution. Although the common diagnostic elytral pattern makes these beetles superficially easily distinguishable, there are marked resemblances in their morphology. By applying an integrative taxonomy approach, we verify the taxonomic status of A. forticornis and A. galusoi and test the possibility of hybridisation. The results of our comprehensive morphological examination and molecular analyses targeting mitochondrial (COI, 16S rRNA) and nuclear (arginine kinase, CAD, 28S rRNA) genes support the specific status of these taxa. In mtDNA, the putative hybrid specimen is more closely related to A. forticornis, which indicates that a female of this species bred with a male of A. galusoi. The supposed ecotypic variability (manifested in elytral pattern) is apparently not related to the topography or spatial structure of habitat. The phylogenetic hypothesis supports a peripatric scenario in which populations of the common ancestor had been separated for a significant period of time, first by the mountain ranges of Dzungarian Alatau (~11–10 Ma) and ultimately by the Ili River. The diverged taxa came into recent secondary contact following these isolation events. The methodology presented herein can be widely applied to identify first-generation hybridisation in Coleoptera.     

Genomic footprints of an old affair: Single nucleotide polymorphism data reveal historical hybridization and the subsequent evolution of reproductive barriers in two recently diverged grasshoppers with partly overlapping distributions

Secondary contact in close relatives can result in hybridization and the admixture of previously isolated gene pools. However, after an initial period of hybridization, reproductive isolation can evolve through different processes and lead to the interruption of gene flow and the completion of the speciation process. Omocestus minutissimus and O. uhagonii are two closely related grasshoppers with partially overlapping distributions in the Central System mountains of the Iberian Peninsula. To analyse spatial patterns of historical and/or contemporary hybridization between these two taxa and understand how species boundaries are maintained in the region of secondary contact, we sampled sympatric and allopatric populations of the two species and obtained genome‐wide single nucleotide polymorphism data using a restriction site‐associated DNA sequencing approach. We used Bayesian clustering analyses to test the hypothesis of contemporary hybridization in sympatric populations and employed a suite of phylogenomic approaches and a coalescent‐based simulation framework to evaluate alternative hypothetical scenarios of interspecific gene flow. Our analyses rejected the hypothesis of contemporary hybridization but revealed past introgression in the area where the distributions of the two species overlap. Overall, these results point to a scenario of historical gene flow after secondary contact followed by the evolution of reproductive isolation that currently prevents hybridization among sympatric populations.     

Deep mitochondrial introgression and hybridization among ecologically divergent vole species

The completion of speciation is typically difficult to ascertain in rapidly diverging taxa but the amount of hybridization and gene flow in sympatry or parapatry contains important information about the level of reproductive isolation achieved. Here, we examined the progress in speciation between the Mediterranean (Microtus duodecimcostatus) and the Lusitanian pine vole (M. lusitanicus), which are part of the most rapid radiation of species known in mammals. These two Iberian pine voles are classified as separate species because of differences in morphology and ecology, but relatively many ambiguous individuals can be found in sympatric conditions. Our phylogenetic analyses of rangewide data from the mitochondrial cytochrome b gene (mtDNA) demonstrated high levels of diversity and a basal separation in two parapatric lineages. However, mtDNA affiliation was at odds with morphological classification or geographical distribution of the taxa. In contrast, statistical analyses of microsatellites (nucDNA) showed two clear genetic clusters in allopatry and sympatry generally matching morphological classification. This cytonuclear discordance over a large geographic area suggests historical introgression of mtDNA from M. duodecimcostatus to M. lusitanicus. There was statistical evidence for at least two recent hybrids in the sympatry zone but gene flow is apparently low given clear‐cut differences in nucDNA. Our results indicate a relatively advanced speciation process in these Iberian pine voles without fully established reproductive isolation. This situation enables use of combined population genomic and experimental approaches for the separation of patterns and mechanisms in the ongoing explosive diversification of these and other Arvicoline rodents in the future.     

Speciation,population structure,and demographic history of the Mojave Fringe‐toed Lizard (Uma scoparia), a species of conservation concern

The North American deserts were impacted by both Neogene plate tectonics and Quaternary climatic fluctuations, yet it remains unclear how these events influenced speciation in this region. We tested published hypotheses regarding the timing and mode of speciation, population structure, and demographic history of the Mojave Fringe‐toed Lizard (Uma scoparia), a sand dune specialist endemic to the Mojave Desert of California and Arizona. We sampled 109 individual lizards representing 22 insular dune localities, obtained DNA sequences for 14 nuclear loci, and found that U. scoparia has low genetic diversity relative to the U. notata species complex, comparable to that of chimpanzees and southern elephant seals. Analyses of genotypes using Bayesian clustering algorithms did not identify discrete populations within U. scoparia. Using isolation‐with‐migration (IM) models and a novel coalescent‐based hypothesis testing approach, we estimated that U. scoparia diverged from U. notata in the Pleistocene epoch. The likelihood ratio test and the Akaike Information Criterion consistently rejected nested speciation models that included parameters for migration and population growth of U. scoparia. We reject the Neogene vicariance hypothesis for the speciation of U. scoparia and define this species as a single evolutionarily significant unit for conservation purposes.     

Evolutionary and demographic history of the Californian scrub white oak species complex: an integrative approach

Understanding the factors promoting species formation is a major task in evolutionary research. Here, we employ an integrative approach to study the evolutionary history of the Californian scrub white oak species complex (genus Quercus). To infer the relative importance of geographical isolation and ecological divergence in driving the speciation process, we (i) analysed inter‐ and intraspecific patterns of genetic differentiation and employed an approximate Bayesian computation (ABC) framework to evaluate different plausible scenarios of species divergence. In a second step, we (ii) linked the inferred divergence pathways with current and past species distribution models (SDMs) and (iii) tested for niche differentiation and phylogenetic niche conservatism across taxa. ABC analyses showed that the most plausible scenario is the one considering the divergence of two main lineages followed by a more recent pulse of speciation. Genotypic data in conjunction with SDMs and niche differentiation analyses support that different factors (geography vs. environment) and modes of speciation (parapatry, allopatry and maybe sympatry) have played a role in the divergence process within this complex. We found no significant relationship between genetic differentiation and niche overlap, which probably reflects niche lability and/or that multiple factors, have contributed to speciation. Our study shows that different mechanisms can drive divergence even among closely related taxa representing early stages of species formation and exemplifies the importance of adopting integrative approaches to get a better understanding of the speciation process.     

Divergent selection and drift shape the genomes of two avian sister species spanning a saline–freshwater ecotone

The role of species divergence due to ecologically based divergent selection—or ecological speciation—in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline–freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes—the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)—to evaluate the influence of selective and demographic processes in shaping genome‐wide patterns of divergence. Genome‐wide divergence between these two recently diverged sister species was notably high (genome‐wide F_ST = 0.32). Against a background of high genome‐wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism—all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.     

Stable climate corridors promote gene flow in the Cape sand snake species complex (Psammophiidae)

Climate shifts during the Quaternary Period have driven changes in regional range dynamics for many species, influencing population structure of species and in some cases promoting speciation. Within southern Africa, the psammophine snakes Psammophis trinasalis and P. namibensis were historically considered subspecies of P. leightoni but were elevated to species rank based on ecological differences. Preliminary phylogenetic analyses suggested intraspecific, not interspecific genetic variation between these taxa, but this finding was based on very limited data and could not be confirmed. To assess the level of genetic differentiation within the P. leightoni species complex, we explored the evolutionary history of these snakes by combining phylogenetic analyses, species distribution modelling and an examination of morphology. We generated a comprehensive, multi-gene phylogeny for Psammophis that included wider geographic sampling of the three species in the complex. Using this phylogeny, Bayesian and distance-based species delimitation analyses showed intraspecific, not interspecific divergences between taxa in the complex, suggesting that they collectively represent a single taxon. Furthermore, non-metric multidimensional scaling analysis of scalation characters showed no differences between the species. Moreover, palaeo-modelling at three time periods since the last interglacial period suggest that there have been varying levels of connectivity between these taxa, which has likely facilitated gene flow between them. Given the evidence, we propose that the P. leightoni complex represents a single species and therefore formally synonymise the three species.     

Ecological speciation in anemone‐associated snapping shrimps (Alpheus armatus species complex)

Divergent natural selection driven by competition for limited resources can promote speciation, even in the presence of gene flow. Reproductive isolation is more likely to result from divergent selection when the partitioned resource is closely linked to mating. Obligate symbiosis and host fidelity (mating on or near the host) can provide this link, creating ideal conditions for speciation in the absence of physical barriers to dispersal. Symbiotic organisms often experience competition for hosts, and host fidelity ensures that divergent selection for a specific host or host habitat can lead to speciation and strengthen pre‐existing reproductive barriers. Here, we present evidence that diversification of a sympatric species complex occurred despite the potential for gene flow and that partitioning of host resources (both by species and by host habitat) has contributed to this diversification. Four species of snapping shrimps (Alpheus armatus, A. immaculatus, A. polystictus and A. roquensis) are distributed mainly sympatrically in the Caribbean, while the fifth species (A. rudolphi) is restricted to Brazil. All five species are obligate commensals of sea anemones with a high degree of fidelity and ecological specificity for host species and habitat. We analysed sequence data from 10 nuclear genes and the mitochondrial COI gene in 11–16 individuals from each of the Caribbean taxa and from the only available specimen of the Brazilian taxon. Phylogenetic analyses support morphology‐based species assignments and a well‐supported Caribbean clade. The Brazilian A. rudolphi is recovered as an outgroup to the Caribbean taxa. Isolation–migration coalescent analysis provides evidence for historical gene flow among sympatric sister species. Our data suggest that both selection for a novel host and selection for host microhabitat may have promoted diversification of this complex despite gene flow.     

Floral scent and species divergence in a pair of sexually deceptive orchids

Speciation is typically accompanied by the formation of isolation barriers between lineages. Commonly, reproductive barriers are separated into pre‐ and post‐zygotic mechanisms that can evolve with different speed. In this study, we measured the strength of different reproductive barriers in two closely related, sympatric orchids of the Ophrys insectifera group, namely Ophrys insectifera and Ophrys aymoninii to infer possible mechanisms of speciation. We quantified pre‐ and post‐pollination barriers through observation of pollen flow, by performing artificial inter‐ and intraspecific crosses and analyzing scent bouquets. Additionally, we investigated differences in mycorrhizal fungi as a potential extrinsic factor of post‐zygotic isolation. Our results show that floral isolation mediated by the attraction of different pollinators acts apparently as the sole reproductive barrier between the two orchid species, with later‐acting intrinsic barriers seemingly absent. Also, the two orchids share most of their fungal mycorrhizal partners in sympatry, suggesting little or no importance of mycorrhizal symbiosis in reproductive isolation. Key traits underlying floral isolation were two alkenes and wax ester, present predominantly in the floral scent of O. aymoninii. These compounds, when applied to flowers of O. insectifera, triggered attraction and a copulation attempt of the bee pollinator of O. aymoninii and thus led to the (partial) breakdown of floral isolation. Based on our results, we suggest that adaptation to different pollinators, mediated by floral scent, underlies species isolation in this plant group. Pollinator switches may be promoted by low pollination success of individuals in dense patches of plants, an assumption that we also confirmed in our study.     

Comparative phylogeography of two related plant species with overlapping ranges in Europe, and the potential effects of climate change on their intraspecific genetic diversity

Background  

The aim of the present study was to use a combined phylogeographic and species distribution modelling approach to compare the glacial histories of two plant species with overlapping distributions, Orthilia secunda (one-sided wintergreen) and Monotropa hypopitys (yellow bird's nest). Phylogeographic analysis was carried out to determine the distribution of genetic variation across the range of each species and to test whether both correspond to the "classic" model of high diversity in the south, with decreasing diversity at higher latitudes, or whether the cold-adapted O. secunda might retain more genetic variation in northern populations. In addition, projected species distributions based on a future climate scenario were modelled to assess how changes in the species ranges might impact on total intraspecific diversity in both cases.     

Phylogenomics unravels Quaternary vicariance and allopatric speciation patterns in temperate‐montane plant species: A case study on the Ranunculus auricomus species complex

The time frame and geographical patterns of diversification processes in European temperate‐montane herbs are still not well understood. We used the sexual species of the Ranunculus auricomus complex as a model system to understand how vicariance versus dispersal processes in the context of Pleistocene climatic fluctuations have triggered speciation in temperate‐montane plant species. We used target enrichment sequence data from about 600 nuclear genes and coalescent‐based species tree inference methods to resolve phylogenetic relationships among the sexual taxa of the complex. We estimated absolute divergence times and, using ancestral range reconstruction, we tested if speciation was enhanced by vicariance or by dispersal processes. Phylogenetic relationships among taxa were fully resolved with some incongruence in the position of the tetraploid R. marsicus. Speciation events took place in a very short time at the end of the Mid‐Pleistocene Transition (830–580 thousand years ago [ka]). A second wave of intraspecific geographical differentiation occurred at the end of the Riss glaciation or during the Eemian interglacial between 200 and 100 ka. Ancestral range reconstruction suggests a widespread European ancestor of the R. auricomus complex. Vicariance has triggered allopatric speciation in temperate‐montane plant species during the climatic deterioration that occurred during the last phase of the Mid‐Pleistocene Transition. Vegetation restructuring from forest into tundra could have confined these forest species into isolated glacial macro‐ and microrefugia. During subsequent warming periods, range expansions of these species could have been hampered by apomictic derivatives and by other congeneric competitors in the same habitat.     

Early events in speciation: Cryptic species of Drosophila aldrichi

Understanding the earliest events in speciation remains a major challenge in evolutionary biology. Thus identifying species whose populations are beginning to diverge can provide useful systems to study the process of speciation. Drosophila aldrichi, a cactophilic fruit fly species with a broad distribution in North America, has long been assumed to be a single species owing to its morphological uniformity. While previous reports either of genetic divergence or reproductive isolation among different D. aldrichi strains have hinted at the existence of cryptic species, the evolutionary relationships of this species across its range have not been thoroughly investigated. Here we show that D. aldrichi actually is paraphyletic with respect to its closest relative, Drosophila wheeleri, and that divergent D. aldrichi lineages show complete hybrid male sterility when crossed. Our data support the interpretation that there are at least two species of D. aldrichi, making these flies particularly attractive for studies of speciation in an ecological and geographical context.     

Integrative analyses on Western Palearctic Lasiommata reveal a mosaic of nascent butterfly species

Satyrinae butterflies occurring in the Mediterranean apparently have reduced gene flow over sea straits, and for several species, recent wide-scale biodiversity surveys indicate the existence of divergent mitochondrial lineages. Here, we apply an integrative approach and examine the phylogeography of the genus Lasiommata in the Western Palearctic. Our research comprised molecular analyses (mitochondrial and nuclear DNA) and geometric morphometrics (wings and genitalia) for two main species groups, and a comparative GMYC analysis, based on COI, of all the tribes within Satyrinae from this region. The GMYC approach revealed a particularly fast coalescence rate in the Parargina subtribe. The Lasiommata group was divided into 12 evolutionary significant units: six clades for the L. maera species group, five for the L. megera species group, and one for L. petropolitana, with divergences of about 1%. The patterns of COI were mirrored by ITS2 in L. maera, but the two markers were generally inconsistent in L. megera. On the contrary, morphological differences were coherent with the results of COI for L. megera, but less clearly so for L. maera. L. paramegaera and L. meadewaldoi were considerably differentiated for all the analyzed markers and likely proceeded faster in the process of speciation because of geographic isolation and reduced effective population size, rendering the rest paraphyletic. Our study illustrates the continuous nature of speciation and the difficulties of delimiting species. In Lasiommata, the recognition of taxa as diverging lineages or distinct, possibly paraphyletic species, mostly depends on the criteria adopted by different species concepts.     

Evaluation of species boundaries in sympatric and parapatric populations of Mesoamerican toads

Approaches that integrate multiple independent, yet complimentary, lines of evidence have been effectively utilized to identify and evaluate species diversity. Integrative approaches are especially useful in taxa that exhibit cryptic diversity and are highly morphologically conserved, as well as organisms whose distributions may be sympatric or parapatric. The Incilius coccifer complex in Honduras is comprised of three putative taxa: I. coccifer, I. ibarrai and I. porteri. The taxonomy of the I. coccifer complex has been a source of debate among specialists, with some recognizing three species, while others choose to recognize one widespread taxon. To assess species boundaries and evaluate the taxonomic structure for the I. coccifer complex, we utilized a combination of comprehensive field sampling, molecular phylogenetics and macroecological modelling. Using 58 samples representing all three putative taxa, we generated sequence data from the mitochondrial loci 16S and COI in order to assess genetic diversity and phylogenetic relationships, and tested putative species boundaries using General Mixed Yule‐Coalescent models. To evaluate macroecological differences in the distribution of putative taxa, we utilized maximum entropy modelling and identified areas of suitable and non‐suitable habitat, as well as identifying potential areas of overlap between species habitats. We recovered three clades that broadly correspond to the three named taxa that, while being monophyletic, are separated by relatively small genetic distances. Species distribution models revealed that I. coccifer is macroecologically different than the other two taxa, but that I. ibarrai and I. porteri are highly similar. We uncovered cases of sympatry between pairs of species in at least three localities in Honduras, suggesting the potential for hybridization in these closely related lineages.     

Following the cold: geographical differentiation between interglacial refugia and speciation in the arcto‐alpine species complex Bombus monticola (Hymenoptera: Apidae)

Cold‐adapted species are expected to have reached their largest distribution range during a part of the Ice Ages whereas postglacial warming has led to their range contracting toward high‐latitude and high‐altitude areas. This has resulted in an extant allopatric distribution of populations and possibly to trait differentiations (selected or not) or even speciation. Assessing inter‐refugium differentiation or speciation remains challenging for such organisms because of sampling difficulties (several allopatric populations) and disagreements on species concept. In the present study, we assessed postglacial inter‐refugia differentiation and potential speciation among populations of one of the most common arcto‐alpine bumblebee species in European mountains, Bombus monticola Smith, 1849. Based on mitochondrial DNA/nuclear DNA markers and eco‐chemical traits, we performed integrative taxonomic analysis to evaluate alternative species delimitation hypotheses and to assess geographical differentiation between interglacial refugia and speciation in arcto‐alpine species. Our results show that trait differentiations occurred between most Southern European mountains (i.e. Alps, Balkan, Pyrenees, and Apennines) and Arctic regions. We suggest that the monticola complex actually includes three species: B. konradini   stat.n. status distributed in Italy (Central Apennine mountains), B. monticola with five subspecies, including B. monticola mathildis   ssp.n. distributed in the North Apennine mountains ; and B. lapponicus . Our results support the hypothesis that post‐Ice Age periods can lead to speciation in cold‐adapted species through distribution range contraction. We underline the importance of an integrative taxonomic approach for rigorous species delimitation, and for evolutionary study and conservation of taxonomically challenging taxa.     

Fuzzy species limits in Mediterranean gorgonians (Cnidaria,Octocorallia): inferences on speciation processes

The study of the interplay between speciation and hybridization is of primary importance in evolutionary biology. Octocorals are ecologically important species whose shallow phylogenetic relationships often remain to be studied. In the Mediterranean Sea, three congeneric octocorals can be observed in sympatry: Eunicella verrucosa, Eunicella cavolini and Eunicella singularis. They display morphological differences and E. singularis hosts photosynthetic Symbiodinium, contrary to the two other species. Two nuclear sequence markers were used to study speciation and gene flow between these species, through network analysis and Approximate Bayesian Computation (ABC). Shared sequences indicated the possibility of hybridization or incomplete lineage sorting. According to ABC, a scenario of gene flow through secondary contact was the best model to explain these results. At the intraspecific level, neither geographical nor ecological isolation corresponded to distinct genetic lineages in E. cavolini. These results are discussed in the light of the potential role of ecology and genetic incompatibilities in the persistence of species limits.     

Effects of environmental factors on pollen production in anemophilous woody species

The aim of this study was to estimate the amount of pollen produced by anemophilous woody taxa with allergenic properties and with considerable contribution in the concentration of pollen in the air of a Mediterranean city (Thessaloniki, Greece). The taxa selected are Corylus avellana, Cupressus sempervirens var. horizontalis and var. pyramidalis, Olea europaea and Platanus orientalis; each was studied in more than one sampling stations differing in elevation, direction or both. O. europaea produced the highest number of pollen grains per flower (1.3 × 10⁵ ± 0.1 × 10⁵) and P. orientalis the highest per inflorescence (3.3 × 10⁶ ± 0.2 × 10⁶). At the level of crown, pollen grains produced were of the order of 10⁹ per surface/volume unit for O. europaea and the two C. sempervirens varieties; for the other two taxa, they were of the order of 10⁶. Pollen production was lower at higher elevation and northern direction and depended on the size of the floral unit sampled (flower for O. europaea, inflorescence for all other species): the bigger the floral unit, the more pollen it contained. Our results and reports from other areas, where C. sempervirens and O. europaea grow, show that these two Mediterranean species produce comparable amounts of pollen at the levels of inflorescence or flower, respectively, wherever they occur.     

Population genetic structure of Oryza rufipogon and Oryza nivara: implications for the origin of O. nivara

Ecological speciation plays a primary role in driving species divergence and adaptation. Oryza rufipogon and Oryza nivara are two incipient species at the early stage of speciation with distinct differences in morphology, life history traits and habitat preference, and therefore provide a unique model for the study of ecological speciation. However, the population genetic structure of the ancestral O. rufipogon has been controversial despite substantial study, and the origin of the derivative O. nivara remains unclear. Here, based on sequences of 10 nuclear and two chloroplast loci from 26 wild populations across the entire geographic ranges of the two species, we conducted comprehensive analyses using population genetics, phylogeography and species distribution modelling (SDM) approaches. In addition to supporting the two previously reported major subdivisions, we detected four genetically distinct groups within O. rufipogon and found no correlation between the genetic groups and either species identity or geographical regions. The SDM clearly showed substantial change in the distribution range of O. rufipogon in history, demonstrating that the repeated extinction and colonization of local populations due to multiple glacial–interglacial cycles during the Quaternary was most likely the main factor shaping the confounding population genetic structure of O. rufipogon. Moreover, we found significant differences between the two species in climate preferences, suggestive of an important role for climatic factors in the adaptation, persistence and expansion of O. nivara. Finally, based on the genetic pattern and dynamics of the O. nivara populations, we hypothesize that O. nivara might have independently originated multiple times from different O. rufipogon populations.