THE BIOLOGY OF SPECIATION

Local adaptation, recombination, and the fate of neopolyploids

Polyploidy is widely recognized as an important speciation mechanism because it isolates tetraploids from their diploid progenitors. Polyploidy also provides new genetic material that may facilitate adaptive evolution. However, new mutations are more likely to arise after a neopolyploid has already successfully invaded a population. Thus, the role of adaptive forces in establishing a polyploid remains unclear. One solution to this apparent paradox may lie in the capacity of polyploids to suppress recombination among preexisting locally adapted alleles. The local adaptation mechanism requires that spatially heterogeneous selection acts on multiple loci and that gene flow introduces maladapted alleles to the population where the polyploid forms. The mechanism requires neither strong genetic drift nor any intrinsic benefit of genome doubling and can accommodate any mode of gene action. A unique prediction of the mechanism is that adaptive alleles should predate polyploidization, a pattern consistent with observations from a few well-studied polyploids. The mechanism is also consistent with the coexistence of both diploid and tetraploid cytotypes, fitness heterogeneity among independently derived polyploids, and the prevalence of outcrossing among older polyploids. The local adaptation mechanism also makes novel predictions about circumstances favoring polyploid invasions that can be tested using molecular genetic or comparative approaches.

Ice age-driven range shifts of diploids and expanding autotetraploids of Biscutella laevigata within a conserved niche

Early studies of the textbook mixed-ploidy system Biscutella laevigata highlighted diploids restricted to never-glaciated lowlands and tetraploids at high elevations across the European Alps, promoting the hypothesis that whole-genome duplication (WGD) is advantageous under environmental changes. Here we addressed long-held hypotheses on the role of hybridisation at the origin of the tetraploids, their single vs multiple origins, and whether a shift in climatic niche accompanied WGD. Climatic niche modelling together with spatial genetics and coalescent modelling based on ddRAD-seq genotyping of 17 diploid and 19 tetraploid populations was used to revisit the evolution of this species complex in space and time. Diploids differentiated into four genetic lineages corresponding to allopatric glacial refugia at the onset of the last ice age, whereas tetraploids displaying tetrasomic inheritance formed a uniform group that originated from southern diploids before the last glacial maximum. Derived from diploids occurring at high elevation, autotetraploids likely inherited their adaptation to high elevation rather than having evolved it through or after WGD. They further presented considerable postglacial expansion across the Alps and underwent admixture with diploids. Although the underpinnings of the successful expansion of autotetraploids remain elusive, differentiation in B. laevigata was chiefly driven by the glacial history of the Alps.

The Ecology of Hybrid Incompatibilities

Ecologically mediated selection against hybrids, caused by hybrid phenotypes fitting poorly into available niches, is typically viewed as distinct from selection caused by epistatic Dobzhansky-Muller hybrid incompatibilities. Here, we show how selection against transgressive phenotypes in hybrids manifests as incompatibility. After outlining our logic, we summarize current approaches for studying ecology-based selection on hybrids. We then quantitatively review QTL-mapping studies and find traits differing between parent taxa are typically polygenic. Next, we describe how verbal models of selection on hybrids translate to phenotypic and genetic fitness landscapes, highlighting emerging approaches for detecting polygenic incompatibilities. Finally, in a synthesis of published data, we report that trait transgression-and thus possibly extrinsic hybrid incompatibility in hybrids-escalates with the phenotypic divergence between parents. We discuss conceptual implications and conclude that studying the ecological basis of hybrid incompatibility will facilitate new discoveries about mechanisms of speciation.

The Cretaceous Diatom Database: A tool for investigating early diatom evolution

The Cretaceous period is the time of the first appearance of the diatoms in the fossil record. These fossils give us direct evidence of the age and early evolution of the diatom lineage. The fossil record, however, is incomplete and therefore often extrapolated through time-calibrated phylogenies. These two approaches offer different perspectives on the early evolution of diatoms, which is still poorly understood. We compiled the first comprehensive Cretaceous Diatom Database, a tool to investigate the taxonomy, diversity, and occurrence of the earliest known diatom lineages. To further aid the integration and use of the oldest diatom fossils in molecular clock analyses, we present a set of well-documented Cretaceous fossils that can be placed onto molecular phylogenetic trees of extant and extinct species, making them ideal candidates for the calibration of molecular clocks. The analysis of the fossil record and the Cretaceous Diatom Database revealed Cretaceous diversity is substantially greater than previously thought, yet considerable taxonomic work is still needed. The Cretaceous Diatom Database and the list of Cretaceous fossils for calibrating molecular clocks represent valuable resources for future evolutionary and taxonomic studies of modern and fossil diatoms.

Geographical isolation and hyperendemicity of Hepatozoon felis: Epidemiological scenario in Skopelos, Greece, and phylogenetic analysis

Feline hepatozoonosis is a vector-borne disease caused by different species of the genus Hepatozoon, i.e. Hepatozoon felis, Hepatozoon silvestris and Hepatozoon canis. Knowledge on the biology, epidemiology and taxonomy of Hepatozoon spp. is still limited, despite the fact that the number of documented Hepatozoon spp. infections in domestic cats increased in recent years in different countries. This study was carried out to evaluate the prevalence and the genetic profile of Hepatozoon spp. in cats living on the island of Skopelos, Greece. Individual blood samples were collected from 54 owned cats and were subjected to Giemsa-stained blood smear examination to investigate the presence of Hepatozoon spp. gamonts and to a specific PCR protocol targeting the 18S rRNA gene of Hepatozoon. A total of 45 cats (83.3%) were found infected by Hepatozoon spp. by at least one of the methods applied. In particular, 43 (79.6%) of the cats were PCR-positive, and in 6 (11.1%) cats gamonts of Hepatozoon spp. were found in the blood smears. A total of 26 H. felis sequences were obtained and the presence of three undescribed single nucleotide polymorphisms were detected. The present results indicate that H. felis species complex may be hyperendemic in isolated/confined areas. In such contexts, geographical isolation may favor the origin of new genotypes or haplotypes or even new species.

Museum genomics approach to study the taxonomy and evolution of Woolly-necked storks using historic specimens

The accessibility of genomic tools in evolutionary biology has allowed for a thorough exploration of various evolutionary processes associated with adaptation and speciation. However, genomic studies in natural systems present numerous challenges, reflecting the inherent complexities of studying organisms in their native habitats. The utilization of museum specimens for genomics research has received increased attention in recent times, facilitated by advancements in ancient DNA techniques. In this study, we have utilized a museum genomics approach to analyze historic specimens of Woolly-necked storks (Ciconia spp.) and examine their genetic composition and taxonomic status and explore the evolutionary and adaptive trajectories of populations over the years. The Woolly-necked storks are distributed in Asia and Africa with a taxonomic classification that has been a matter of ambiguity. Asian and African Woollynecks were recently recognized as different species based on their morphological differences; however, their genomic validation was lacking. In this study, we have used ∼70-year-old museum samples for whole-genome population-scale sequencing. Our study has revealed that Asian and African Woollynecks are genetically distinct, consistent with the current taxonomic classification based on morphological features. However, we also found a high genetic divergence between the Asian subspecies Ciconia episcopus neglecta and Ciconia episcopus episcopus, suggesting this classification requires a detailed examination to explore processes of ongoing speciation. Because taxonomic classification directly impacts conservation efforts, and there is evidence of declining populations of Asian Woollynecks in Southeast Asia, our results highlight that population-scale studies are urgent to determine the genetic, ecological, and phylogenetic diversity of these birds.

Phenotypic plasticity rather than ecotypic differentiation explains the broad realized niche of a Neotropical orchid species

Local adaptation is common in plant species, and knowing whether a population is locally adapted has fundamental and applied relevance. However, local adaptation in tropical plants remains largely less studied, and covering this gap is not simple since reciprocal transplantation - the gold standard for detecting local adaptation - is not feasible for most species. Here, we combined genetic, climatic and phenotypic data to investigate ecotypic differentiation, an important aspect of local adaptation, in coastal and inland populations of the orchid Epidendrum fulgens Brongn., a long-lived tropical plant for which reciprocal transplantation would not be feasible. We used nine microsatellite markers to estimate genetic divergence between inland and coastal populations. Moreover, occurrence data and climate data were used to test for differences in the realized niche of those populations. Finally, we assessed saturated water content, leaf specific area, height, and stomatal density in common garden and in situ to investigate the effects of ecotypic differentiation and plasticity on the phenotype. Coastal and inland groups' niches do not overlap, the former occupying a wetter and warmer area. However, this differentiation does not seem to be driven by ecotypic differentiation since there was no positive correlation between genetic structure and climate dissimilarity. Moreover, specific leaf area and leaf saturated water content, which are important phenotypic traits related to soil fertility and drought stress, were rather plastic. We conclude that ecotypic differentiation is absent, since phenotypic plasticity is an important mechanism explaining the niche broadness of this species.

Adaptive divergence, historical population dynamics, and simulation of suitable distributions for Picea Meyeri and P. Mongolica at the whole-genome level

The taxonomic classification of Picea meyeri and P. mongolica has long been controversial. To investigate the genetic relatedness, evolutionary history, and population history dynamics of these species, genotyping-by-sequencing (GBS) technology was utilized to acquire whole-genome single nucleotide polymorphism (SNP) markers, which were subsequently used to assess population structure, population dynamics, and adaptive differentiation. Phylogenetic and population structural analyses at the genomic level indicated that although the ancestor of P. mongolica was a hybrid of P. meyeri and P. koraiensis, P. mongolica is an independent Picea species. Additionally, P. mongolica is more closely related to P. meyeri than to P. koraiensis, which is consistent with its geographic distribution. There were up to eight instances of interspecific and intraspecific gene flow between P. meyeri and P. mongolica. The P. meyeri and P. mongolica effective population sizes generally decreased, and Maxent modeling revealed that from the Last Glacial Maximum (LGM) to the present, their habitat areas decreased initially and then increased. However, under future climate scenarios, the habitat areas of both species were projected to decrease, especially under high-emission scenarios, which would place P. mongolica at risk of extinction and in urgent need of protection. Local adaptation has promoted differentiation between P. meyeri and P. mongolica. Genotype‒environment association analysis revealed 96,543 SNPs associated with environmental factors, mainly related to plant adaptations to moisture and temperature. Selective sweeps revealed that the selected genes among P. meyeri, P. mongolica and P. koraiensis are primarily associated in vascular plants with flowering, fruit development, and stress resistance. This research enhances our understanding of Picea species classification and provides a basis for future genetic improvement and species conservation efforts.

Rapid vertebrate speciation via isolation, bottlenecks, and drift

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish (Cyprinodon tularosa) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global FST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose.

Divergent dynamics of sexual and habitat isolation at the transition between stick insect populations and species

Speciation is often viewed as a continuum along which populations diverge until they become reproductively-isolated species. However, such divergence may be heterogeneous, proceeding in fits and bursts, rather than being uniform and gradual. We show in Timema stick insects that one component of reproductive isolation evolves non-uniformly across this continuum, whereas another does not. Specifically, we use thousands of host-preference and mating trials to study habitat and sexual isolation among 42 pairs of taxa spanning a range of genomic differentiation and divergence time. We find that habitat isolation is uncoupled from genomic differentiation within species, but accumulates linearly with it between species. In contrast, sexual isolation accumulates linearly across the speciation continuum, and thus exhibits similar dynamics to morphological traits not implicated in reproductive isolation. The results show different evolutionary dynamics for different components of reproductive isolation and highlight a special relevance for species status in the process of speciation.

Chromosomal evolution, environmental heterogeneity, and migration drive spatial patterns of species richness in Calochortus (Liliaceae)

We used nuclear genomic data and statistical models to evaluate the ecological and evolutionary processes shaping spatial variation in species richness in Calochortus (Liliaceae, 74 spp.). Calochortus occupies diverse habitats in the western United States and Mexico and has a center of diversity in the California Floristic Province, marked by multiple orogenies, winter rainfall, and highly divergent climates and substrates (including serpentine). We used sequences of 294 low-copy nuclear loci to produce a time-calibrated phylogeny, estimate historical biogeography, and test hypotheses regarding drivers of present-day spatial patterns in species number. Speciation and species coexistence require reproductive isolation and ecological divergence, so we examined the roles of chromosome number, environmental heterogeneity, and migration in shaping local species richness. Six major clades-inhabiting different geographic/climatic areas, and often marked by different base chromosome numbers (n = 6 to 10)-began diverging from each other ~10.3 Mya. As predicted, local species number increased significantly with local heterogeneity in chromosome number, elevation, soil characteristics, and serpentine presence. Species richness is greatest in the Transverse/Peninsular Ranges where clades with different chromosome numbers overlap, topographic complexity provides diverse conditions over short distances, and several physiographic provinces meet allowing immigration by several clades. Recently diverged sister-species pairs generally have peri-patric distributions, and maximum geographic overlap between species increases over the first million years since divergence, suggesting that chromosomal evolution, genetic divergence leading to gametic isolation or hybrid inviability/sterility, and/or ecological divergence over small spatial scales may permit species co-occurrence.

Demographic dynamics and molecular evolution of the rare and endangered subsect. Gerardianae of Pinus: insights from chloroplast genomes and mitochondrial DNA markers

MAIN CONCLUSION

The divergence of subsect. Gerardianae was likely triggered by the uplift of the Qinghai-Tibetan Plateau and adjacent mountains. Pinus bungeana might have probably experienced expansion since Last Interglacial period. Historical geological and climatic oscillations have profoundly affected patterns of nucleotide variability, evolutionary history, and species divergence in numerous plants of the Northern Hemisphere. However, how long-lived conifers responded to geological and climatic fluctuations in East Asia remain poorly understood. Here, based on paternally inherited chloroplast genomes and maternally inherited mitochondrial DNA markers, we investigated the population demographic history and molecular evolution of subsect. Gerardianae (only including three species, Pinus bungeana, P. gerardiana, and P. squamata) of Pinus. A low level of nucleotide diversity was found in P. bungeana (π was 0.00016 in chloroplast DNA sequences, and 0.00304 in mitochondrial DNAs). The haplotype-based phylogenetic topology and unimodal distributions of demographic analysis suggested that P. bungeana probably originated in the southern Qinling Mountains and experienced rapid population expansion since Last Interglacial period. Phylogenetic analysis revealed that P. gerardiana and P. squamata had closer genetic relationship. The species divergence of subsect. Gerardianae occurred about 27.18 million years ago (Mya) during the middle to late Oligocene, which was significantly associated with the uplift of the Qinghai-Tibetan Plateau and adjacent mountains from the Eocene to the mid-Pliocene. The molecular evolutionary analysis showed that two chloroplast genes (psaI and ycf1) were under positive selection, the genetic lineages of P. bungeana exhibited higher transition and nonsynonymous mutations, which were involved with the strongly environmental adaptation. These findings shed light on the population evolutionary history of white pine species and provide striking insights for comprehension of their species divergence and molecular evolution.

Temporal Variation in Selection Influences Microgeographic Local Adaptation

AbstractEcological heterogeneity can lead to local adaptation when populations exhibit fitness trade-offs among habitats. However, the degree to which local adaptation is affected by the spatial and temporal scale of environmental variation is poorly understood. A multiyear reciprocal transplant experiment was performed with populations of the annual plant Leptosiphon parviflorus living on adjacent serpentine and nonserpentine soil. Local adaptation over this small geographic scale was observed, but there were differences in the temporal variability of selection across habitats. On serpentine soil, the local population had a consistently large survival advantage, presumably as a result of the temporal stability in selection imposed by soil cation content. In contrast, a fecundity advantage was observed for the sandstone population on its native soil type but only in the two study years with the highest rainfall. A manipulative greenhouse experiment demonstrated that the fitness advantage of the sandstone population in its native soil type depends critically on water availability. The temporal variability in local adaptation driven by variation in precipitation suggests that continued drought conditions have the potential to erode local adaptation in these populations. These results show how different selective factors can influence spatial and temporal patterns of variation in fitness trade-offs.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

Chemical signaling glands are unlinked to species diversification in lizards

Sexual selection has long been thought to increase species diversification. Sexually selected traits, such as sexual signals that contribute to reproductive isolation, were thought to promote diversification. However, studies exploring links between sexually selected traits and species diversification have thus far primarily focused on visual or acoustic signals. Many animals often employ chemical signals (i.e., pheromones) for sexual communications, but large-scale analyses on the role of chemical communications in driving species diversification have been missing. Here, for the first time, we investigate whether traits associated with chemical communications-the presence of follicular epidermal glands-promote diversification across 6,672 lizard species. In most analyses, we found no strong association between the presence of follicular epidermal glands and species diversification rates, either across all lizard species or at lower phylogenetic scales. Previous studies suggest that follicular gland secretions act as species recognition signals that prevent hybridization during speciation in lizards. However, we show that geographic range overlap was no different in sibling species pairs with and without follicular epidermal glands. Together, these results imply that either follicular epidermal glands do not primarily function in sexual communications or sexually selected traits in general (here chemical communication) have a limited effect on species diversification. In our additional analysis accounting for sex-specific differences in glands, we again found no detectable effect of follicular epidermal glands on species diversification rates. Thus, our study challenges the general role of sexually selected traits in broad-scale species diversification patterns.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Geographic-genomic and geographic-phenotypic differentiation of the Aquilegia viridiflora complex

How species diverge into different lineages is a central issue in evolutionary biology. Despite the increasing evidence indicating that such divergences do not need geographic isolation, the correlation between lineage divergence and the adaptive ecological divergence of phenotype corresponding to distribution is still unknown. In addition, gene flow has been widely detected during and through such diverging processes. We used one widely distributed Aquilegia viridiflora complex as a model system to examine genomic differentiation and corresponding phenotypic variations along geographic gradients. Our phenotypic analyses of 20 populations from northwest to northeast China identified two phenotypic groups along the geographic cline. All examined traits are distinct from each other, although a few intermediate individuals occur in their contacting regions. We further sequenced the genomes of representative individuals of each population. However, four distinct genetic lineages were detected based on nuclear genomes. In particular, we recovered numerous genetic hybrids in the contact regions of four lineages. Gene flow is widespread and continuous between four lineages but much higher between contacting lineages than geographically isolated lineages. Gene flow and natural selection might result in inconsistency between heredity and phenotype. Moreover, many genes with fast lineage-specific mutations were identified to be involved in local adaptation. Our results suggest that both geographic isolation and local selection exerted by the environment and pollinators may together create geographic distributions of phenotypic variations as well as the underlying genomic divergences in numerous lineages.

Multiple pre- and postzygotic components of reproductive isolation between two co-occurring Lysimachia species

Genetic divergence between species depends on reproductive isolation (RI) due to traits that reduce interspecific mating (prezygotic isolation) or are due to reduced hybrid fitness (postzygotic isolation). Previous research found that prezygotic barriers tend to be stronger than postzygotic barriers, but most studies are based on the evaluation of F₁ hybrid fitness in early life cycle stages. We combined field and experimental data to determine the strength of 17 prezygotic and postzygotic reproductive barriers between two Lysimachia species that often co-occur and share pollinators. We assessed postzygotic barriers up to F₂ hybrids and backcrosses. The two species showed near complete RI due to the cumulative effect of multiple barriers, with an uneven and asymmetric contribution to isolation. In allopatry, prezygotic barriers contributed more to reduce gene flow than postzygotic barriers, but their contributions were more similar in sympatry. The strength of postzygotic RI was up to three times lower for F₁ progeny than for F₂ or backcrossed progenies, and RI was only complete when late F₁ stages and either F₂ or backcrosses were accounted for. Our results thus suggest that the relative strength of postzygotic RI may be underestimated when its effects on late stages of the life cycle are disregarded.

Unrecognized diversity, genetic structuring, and phylogeography of the genus Triplophysa (Cypriniformes: Nemacheilidae) sheds light on two opposite colonization routes during Quaternary glaciation that occurred in the Qilian Mountains

In recent years, the influence of historical geological and climatic events on the evolution of flora and fauna in the Tibetan Plateau has been a hot research topic. The Qilian Mountain region is one of the most important sources of biodiversity on the Qinghai-Tibet Plateau. Many species existed in the region during the Pleistocene glacial oscillation, and the complex geographical environment provided suitable conditions for the survival of local species. The shrinkage, expansion, and transfer of the distribution range and population size of species have significant effects on genetic diversity and intraspecific differentiation. To reveal the effects of geological uplift and climate oscillation on the evolution of fish populations in the Qilian Mountains, we investigated the genetic structure, phylogenetic relationship, and phylogeographical characteristics of genus Triplophysa species in the Qilian Mountains using the mitochondrial DNA gene (COI), three nuclear genes (RAG1, sRH, and Myh6) and 11 pairs of nuclear microsatellite markers. We collected 11 species of genus Triplophysa living in the Qilian Mountains, among which Triplophysa hsutschouensis and Triplophysa papillosolabiata are widely distributed in the rivers on the northern slope of the Qilian Mountains. There was a high degree of lineage differentiation among species, and the genetic diversity of endemic species was low. The different geographical groups of T. papillosolabiata presented some allogeneic adaptation and differentiation, which was closely related to the changes in the river system. Except for the population expansion event of T. hsutschouensis during the last glacial period of the Qinghai-Tibet Plateau (0.025 MYA), the population sizes of other plateau loach species remained stable without significant population expansion. Starting from the east and west sides of the Qilian Mountains, T. hsutschouensis, and T. papillosolabiata showed two species colonization routes in opposite directions. The geological events of the uplift of the Qinghai-Tibet Plateau and the climatic oscillation of the Quaternary glaciation had a great influence on the genetic structure of the plateau loach in the Qilian Mountains, which promoted the genetic differentiation of the plateau loach and formed some unique new species. The results of this study have important guiding significance for fish habitat protection in the Qilian Mountains.

Genomic evidence reveals intraspecific divergence of the hot-spring snake (Thermophis baileyi), an endangered reptile endemic to the Qinghai-Tibet plateau

Understanding how and why species evolve requires knowledge on intraspecific divergence. In this study, we examined intraspecific divergence in the endangered hot-spring snake (Thermophis baileyi), an endemic species on the Qinghai-Tibet Plateau (QTP). Whole-genome resequencing of 58 sampled individuals from 15 populations was performed to identify the drivers of intraspecific divergence and explore the potential roles of genes under selection. Our analyses resolved three groups, with major intergroup admixture occurring in regions of group contact. Divergence probably occurred during the Pleistocene as a result of glacial climatic oscillations, Yadong-Gulu rift, and geothermal fields differentiation, while complex gene flow between group pairs reflected a unique intraspecific divergence pattern on the QTP. Intergroup fixed loci involved selected genes functionally related to divergence and local adaptation, especially adaptation to hot spring microenvironments in different geothermal fields. Analysis of structural variants, genetic diversity, inbreeding, and genetic load indicated that the hot-spring snake population has declined to a low level with decreased diversity, which is important for the conservation management of this endangered species. Our study demonstrated that the integration of demographic history, gene flow, genomic divergence genes, and other information is necessary to distinguish the evolutionary processes involved in speciation.

Spatially and temporally varying selection influence species boundaries in two sympatric Mimulus

Spatially and temporally varying selection can maintain genetic variation within and between populations, but it is less well known how these forces influence divergence between closely related species. We identify the interaction of temporal and spatial variation in selection and their role in either reinforcing or eroding divergence between two closely related Mimulus species. Using repeated reciprocal transplant experiments with advanced generation hybrids, we compare the strength of selection on quantitative traits involved in adaptation and reproductive isolation in Mimulus guttatus and Mimulus laciniatus between two years with dramatically different water availability. We found strong divergent habitat-mediated selection on traits in the direction of species differences during a drought in 2013, suggesting that spatially varying selection maintains species divergence. However, a relaxation in divergent selection on most traits in an unusually wet year (2019), including flowering time, which is involved in pre-zygotic isolation, suggests that temporal variation in selection may weaken species differences. Therefore, we find evidence that temporally and spatially varying selection may have opposing roles in mediating species boundaries. Given our changing climate, future growing seasons are expected to be more similar to the dry year, suggesting that in this system climate change may actually increase species divergence.

Initiation of speciation across multiple dimensions in a rock-restricted, tropical lizard

Population isolation and concomitant genetic divergence, resulting in strong phylogeographical structure, is a core aspect of speciation initiation. If and how speciation then proceeds and ultimately completes depends on multiple factors that mediate reproductive isolation, including divergence in genomes, ecology and mating traits. Here we explored these multiple dimensions in two young (Plio-Pleistocene) species complexes of gekkonid lizards (Heteronotia) from the Kimberley-Victoria River regions of tropical Australia. Using mitochondrial DNA screening and exon capture phylogenomics, we show that the rock-restricted Heteronotia planiceps exhibits exceptional fine-scale phylogeographical structure compared to the codistributed habitat generalist Heteronotia binoei. This indicates pervasive population isolation and persistence in the rock-specialist, and thus a high rate of speciation initiation across this geographically complex region, with levels of genomic divergence spanning the "grey zone" of speciation. Proximal lineages of H. planiceps were often separated by different rock substrates, suggesting a potential role for ecological isolation; however, phylogenetic incongruence and historical introgression were inferred between one such pair. Ecomorphological divergence among lineages within both H. planiceps and H. binoei was limited, except that limestone-restricted lineages of H. planiceps tended to be larger than rock-generalists. By contrast, among-lineage divergence in the chemical composition of epidermal pore secretions (putative mating trait) exceeded ecomorphology in both complexes, but with less trait overlap among lineages in H. planiceps. This system-particularly the rock-specialist H. planiceps-highlights the role of multidimensional divergence during incipient speciation, with divergence in genomes, ecomorphology and chemical signals all at play at very fine spatial scales.

Climate influences the genetic structure and niche differentiation among populations of the olive field mouse Abrothrix olivacea (Cricetidae: Abrotrichini)

Even when environmental variation over time and space is commonly considered as an important driver of population divergence, few evaluations of intraspecific genetic variation explicitly assess whether observed structure has been caused by or is correlated with landscape heterogeneity. Several phylogeographic studies have characterized the mitochondrial diversity of Abrothrix olivacea, but none has incorporated landscape genetics analyses and ecological niche modeling, leaving a gap in the understanding of the species evolutionary history. Here, these aspects were addressed based on 186 single nucleotide polymorphisms, extracted from sequences of 801 bp of Cytb gene, gathered from 416 individuals collected at 103 localities in Argentina and Chile. Employing multivariate statistical analyses (gPCA, Mantel and Partial Mantel Tests, Procrustes Analysis, and RDA), associations between genetic differences and geographic and climatic distances were evaluated. Presence data was employed to estimate the potential geographic distribution of this species during historical and contemporary climatic scenarios, and to address differences among the climatic niches of their main mitochondrial lineages. The significant influence of landscape features in structuring mitochondrial variability was evidenced at different spatial scales, as well as the role of past climatic dynamics in driving geographic range shifts, mostly associated to Quaternary glaciations. Overall, these results suggest that throughout geographic range gene flow is unevenly influenced by climatic dissimilarity and the geographic distancing, and that studied lineages do not exhibit distributional signals of climatic niche conservatism. Additionally, genetic differentiation occurred by more complex evolutionary processes than mere disruption of gene flow or drift.

Ecology and the origin of non-ephemeral species

AbstractResearch over the past three decades has shown that ecology-based extrinsic reproductive barriers can rapidly arise to generate incipient species-but such barriers can also rapidly dissolve when environments change, resulting in incipient species collapse. Understanding the evolution of unconditional, "intrinsic" reproductive barriers is therefore important for understanding the longer-term buildup of biodiversity. In this article, we consider ecology's role in the evolution of intrinsic reproductive isolation. We suggest that this topic has fallen into a gap between disciplines: while evolutionary ecologists have traditionally focused on the rapid evolution of extrinsic isolation between co-occurring ecotypes, speciation geneticists studying intrinsic isolation in other taxa have devoted little attention to the ecological context in which it evolves. We argue that for evolutionary ecology to close this gap, the field will have to expand its focus beyond rapid adaptation and its traditional model systems. Synthesizing data from several subfields, we present circumstantial evidence for and against different forms of ecological adaptation as promoters of intrinsic isolation and discuss alternative forces that may be significant. We conclude by outlining complementary approaches that can better address the role of ecology in the evolution of nonephemeral reproductive barriers and, by extension, less ephemeral species.

Eco-Geography and Phenology Are the Major Drivers of Reproductive Isolation in the Royal Irises, a Species Complex in the Course of Speciation

The continuous nature of speciation implies that different species are found at different stages of divergence, from no- to complete reproductive isolation. This process and its underlying mechanisms are best viewed in incipient species. Moreover, the species complex can offer unique insight into how reproductive isolation (RI) has evolved. The royal irises (Iris section Oncocyclus) are a young group of species in the course of speciation, providing an ideal system for speciation study. We quantified pre- and post-zygotic reproductive barriers between the eight Israeli species of this complex and estimated the total RI among them. We tested for both pre-pollination and post-pollination reproductive barriers. Pre-pollination barriers, i.e., eco-geographic divergence and phenological differentiation were the major contributors to RI among the Iris species. On the other hand, post-pollination barriers, namely pollen-stigma interactions, fruit set, and seed viability had negligible contributions to total RI. The strength of RI was not uniform across the species complex, suggesting that species may have diverged at different rates. Overall, this study in a young, recently diverged group of species provides insight into the first steps of speciation, suggesting a crucial role of the pre-zygotic barriers.

Barriers to gene flow in the deepest ocean ecosystems: Evidence from global population genomics of a cosmopolitan amphipod

The deepest marine ecosystem, the hadal zone, hosts endemic biodiversity resulting from geographic isolation and environmental selection pressures. However, the pan-ocean distribution of some fauna challenges the concept that the hadal zone is a series of isolated island-like habitats. Whether this remains true at the population genomic level is untested. We investigated phylogeographic patterns of the amphipod, Bathycallisoma schellenbergi, from 12 hadal features across the Pacific, Atlantic, Indian, and Southern oceans and analyzed genome-wide single-nucleotide polymorphism markers and two mitochondrial regions. Despite a cosmopolitan distribution, populations were highly restricted to individual features with only limited gene flow between topographically connected features. This lack of connectivity suggests that populations are on separate evolutionary trajectories, with evidence of potential cryptic speciation at the Atacama Trench. Together, this global study demonstrates that the shallower ocean floor separating hadal features poses strong barriers to dispersal, driving genetic isolation and creating pockets of diversity to conserve.

The strength of reproductive isolating barriers in seed plants: Insights from studies quantifying premating and postmating reproductive barriers over the past 15 years

Speciation is driven by the evolution of reproductive isolating barriers that reduce, and ultimately prevent, substantial gene flow between lineages. Despite its central role in evolutionary biology, the process can be difficult to study because it proceeds differently among groups and may occur over long timescales. Due to this complexity, we typically rely on generalizations of empirical data to describe and understand the process. Previous reviews of reproductive isolation (RI) in flowering plants have suggested that prezygotic or extrinsic barriers generally have a stronger effect on reducing gene flow compared to postzygotic or intrinsic barriers. Past conclusions have rested on relatively few empirical estimates of RI; however, RI data have become increasingly abundant over the past 15 years. We analyzed data from recent studies quantifying multiple pre- and postmating barriers in plants and compared the strengths of isolating barriers across 89 taxa pairs using standardized RI metrics. Individual prezygotic barriers were on average stronger than individual postzygotic barriers, and the total strength of prezygotic RI was approximately twice that of postzygotic RI. These findings corroborate that ecological divergence and extrinsic factors, as opposed to solely the accumulation of genetic incompatibilities, are important to speciation and the maintenance of species boundaries in plants. Despite an emphasis in the literature on asymmetric postmating and postzygotic RI, we found that prezygotic barriers acted equally asymmetrically. Overall, substantial variability in the strengths of 12 isolating barriers highlights the great diversity of mechanisms that contribute to plant diversification.

Distribution of Biodiversity of Wild Beet Species (Genus Beta L.) in Armenia under Ongoing Climate Change Conditions

The reported annual temperature increase and significant precipitation drop in Armenia impact the country’s ecosystems and biodiversity. The present study surveyed the geographical distribution of the local wild beet species under the ongoing climate change conditions. We showed that B. lomatogona, B. corolliflora and B. macrorhiza are sensitive to climate change and were affected to various degrees, depending on their location. The most affected species was B. lomatogona, which is at the verge of extinction. Migration for ca. 90 and 200–300 m up the mountain belt was recorded for B. lomatogona and B. macrorhiza, respectively. B. corolliflora was found at 100–150 m lower altitudes than in the 1980s. A general reduction in the beet’s population size in the native habitats was observed, with an increased number of plants within the populations, recorded for B. corolliflora and B. macrorhiza. A new natural hybrid Beta x intermedium Aloyan between B. corolliflora and B. macrorhiza was described and confirmed using chloroplast DNA trnL-trnF intergenic spacer (LF) and partially sequenced alcohol dehydrogenase (adh) of nuclear DNA. An overview of the wild beets reported in Armenia with the taxonomic background, morphological features, and distribution is provided. Conservation measures for preservation of these genetic resources are presented.

Flexible drought deciduousness in a neotropical understory herb

PREMISE

Adaptive divergence across environmental gradients is a key driver of speciation. Precipitation seasonality gradients are common in the tropics, yet drought adaptation is nearly unexplored in neotropical understory herbs. Here, we examined two recently diverged neotropical spiral gingers, one adapted to seasonal drought and one reliant on perennial water, to uncover the basis of drought adaptation.

METHODS

We combined ecophysiological trait measurements in the field and greenhouse with experimental and observational assessments of real-time drought response to determine how Costus villosissimus (Costaceae) differs from C. allenii to achieve drought adaptation.

RESULTS

We found that drought-adapted C. villosissimus has several characteristics indicating flexible dehydration avoidance via semi-drought-deciduousness and a fast economic strategy. Although the two species do not differ in water-use efficiency, C. villosissimus has a more rapid growth rate, lower leaf mass per area, lower stem density, higher leaf nitrogen, and a strong trend of greater light-saturated photosynthetic rates. These fast economic strategy traits align with both field-based observations and experimental dry-down results. During drought, C. villosissimus displays facultative drought-deciduousness, losing lower leaves during the dry season and rapidly growing new leaves in the wet season.

CONCLUSIONS

We revealed a drought adaptation strategy that has not, to our knowledge, previously been documented in tropical herbs. This divergent drought adaptation evolved recently and is an important component of reproductive isolation between C. villosissimus and C. allenii, indicating that adaptive shifts to survive seasonal drought may be an underappreciated axis of neotropical understory plant diversification.

Three problems in the genetics of speciation by selection

Speciation is the process by which barriers to gene flow evolve between populations. Although we now know that speciation is largely driven by natural selection, knowledge of the agents of selection and the genetic and genomic mechanisms that facilitate divergence is required for a satisfactory theory of speciation. In this essay, we highlight three advances/problems in our understanding of speciation that have arisen from studies of the genes and genomic regions that underlie the evolution of reproductive isolation. First, we describe how the identification of “speciation” genes makes it possible to identify the agents of selection causing the evolution of reproductive isolation, while also noting that the link between the genetics of phenotypic divergence and intrinsic postzygotic reproductive barriers remains tenuous. Second, we discuss the important role of recombination suppressors in facilitating speciation with gene flow, but point out that the means and timing by which reproductive barriers become associated with recombination cold spots remains uncertain. Third, we establish the importance of ancient genetic variation in speciation, although we argue that the focus of speciation studies on evolutionarily young groups may bias conclusions in favor of ancient variation relative to new mutations.

Linking large-scale genetic structure of three Argynnini butterfly species to geography and environment

Understanding which factors and processes are associated with genetic differentiation within and among species remains a major goal in evolutionary biology. To explore differences and similarities in genetic structure and its association with geographical and climatic factors in sympatric sister species, we conducted a large‐scale (>32° latitude and >36° longitude) comparative phylogeographical study on three Argynnini butterfly species (Speyeria aglaja, Fabriciana adippe and F. niobe) that have similar life histories, but differ in ecological generalism and dispersal abilities. Analyses of nuclear (ddRAD‐sequencing derived SNP markers) and mitochondrial (COI sequences) data revealed differences between species in genetic structure and how genetic differentiation was associated with climatic factors (temperature, solar radiation, precipitation, wind speed). Geographical proximity accounted for much of the variation in nuclear and mitochondrial structure and evolutionary relationships in F. adippe and F. niobe, but only explained the pattern observed in the nuclear data in S. aglaja, for which mitonuclear discordance was documented. In all species, Iberian and Balkan individuals formed genetic clusters, suggesting isolation in glacial refugia and limited postglacial expansion. Solar radiation and precipitation were associated with the genetic structure on a regional scale in all species, but the specific combinations of environmental and geographical factors linked to variation within species were unique, pointing to species‐specific responses to common environments. Our findings show that the species share similar colonization histories, and that the same ecological factors, such as niche breadth and dispersal capacity, covary with genetic differentiation within these species to some extent, thereby highlighting the importance of comparative phylogeographical studies in sympatric sister species.

Reproductive Isolation Among Three Nocturnal Moth-Pollinated Sympatric Habenaria Species (Orchidaceae)

Comparison and quantification of multiple pre- and post-pollination barriers to interspecific hybridization are important to understand the factors promoting reproductive isolation. Such isolating factors have been studied recently in many flowering plant species which seek after the general roles and relative strengths of different pre- and post-pollination barriers. In this study, we quantified six isolating factors (ecogeographic isolation, phenological isolation, pollinator isolation, pollinia-pistil interactions, fruit production, and seed development) that could possibly be acting as reproductive barriers at different stages among three sympatric Habenaria species (H. limprichtii, H. davidii, and H. delavayi). These three species overlap geographically but occupy different microhabitats varying in soil water content. They were isolated through pollinator interactions both ethologically (pollinator preference) and mechanically (pollinia attachment site), but to a variable degree for different species pairs. Interspecific crosses between H. limprichtii and H. davidii result in high fruit set, and embryo development suggested weak post-pollination barriers, whereas bidirectional crosses of H. delavayi with either of the other two species fail to produce fruits. Our results revealed that pollinators were the most important isolating barrier including both ethological and mechanical mechanisms, to maintain the boundaries among these three sympatric Habenaria species. Our study also highlights the importance of a combination of pre-and post-pollination barriers for species co-existence in Orchidaceae.

Phylogenomic analysis does not support a classic but controversial hypothesis of progenitor-derivative origins for the serpentine endemic Clarkia franciscana

Budding speciation involves isolation of marginal populations at the periphery of a species range and is thought to be a prominent mode of speciation in organisms with low dispersal and/or strong local adaptation among populations. Budding speciation is typically evidenced by abutting, asymmetric ranges of ecologically divergent sister species and low genetic diversity in putative budded species. Yet these indirect patterns may be unreliable, instead caused by postspeciation processes such as range or demographic shifts. Nested phylogenetic relationships provide the most conclusive evidence of budding speciation. A putative case of budding speciation in the serpentine endemic Clarkia franciscana and two closely related widespread congeners was studied by Harlan Lewis, Peter Raven, Leslie Gottlieb, and others over a 20-year period, yet the origin of C. franciscana remains controversial. Here, we reinvestigate this system with phylogenomic analyses to determine whether C. franciscana is a recently derived budded species, phylogenetically nested within one of the other two putative progenitor species. In contrast to the hypothesized pattern of relatedness among the three Clarkia species, we find no evidence for recent budding speciation. Instead, the data suggest the three species diverged simultaneously. We urge caution in using contemporary range patterns to infer geographic modes of speciation.

Population and Landscape Genetics Provide Insights Into Species Conservation of Two Evergreen Oaks in Qinghai-Tibet Plateau and Adjacent Regions

The combination of population and landscape genetics can facilitate the understanding of conservation strategy under the changing climate. Here, we focused on the two most diverse and ecologically important evergreen oaks: Quercus aquifolioides and Quercus spinosa in Qinghai-Tibetan Plateau (QTP), which is considered as world's biodiversity hotspot. We genotyped 1,657 individuals of 106 populations at 15 nuclear microsatellite loci throughout the species distribution range. Spatial patterns of genetic diversity were identified by mapping the allelic richness (AR) and locally common alleles (LCA) according to the circular neighborhood methodology. Migration routes from QTP were detected by historical gene flow estimation. The response pattern of genetic variation to environmental gradient was assessed by the genotype-environment association (GEA) analysis. The overall genetic structure showed a high level of intra-species genetic divergence of a strong west-east pattern. The West-to-East migration route indicated the complex demographic history of two oak species. We found evidence of isolation by the environment in Q. aqu-East and Q. spi-West lineage but not in Q. aqu-West and Q. spi-East lineage. Furthermore, priority for conservation should be given to populations that retain higher spatial genetic diversity or isolated at the edge of the distribution range. Our findings indicate that knowledge of spatial diversity and migration route can provide valuable information for the conservation of existing populations. This study provides an important guide for species conservation for two oak species by the integration of population and landscape genetic methods.

Dark-colored Oncocyclus irises in Israel analyzed by AFLP, whole chloroplast genome sequencing and species distribution modeling

The Haynei is one of seven species aggregates (clusters of species having similar flower morphology) recognized in section Oncocyclus of genus Iris. This aggregate, characterized by dark-colored flowers, is represented by six species in Israel and adjacent Jordan. There is, however, no knowledge of the genetic relationship of these species making verification of their taxonomic status impossible. We investigated genetic variation in this group using analysis of whole chloroplast genomes and amplified fragment length polymorphism (AFLP). We also used species distribution modeling (SDM) to predict species ranges under current climatic conditions. We found some population groups within the currently recognized species of section Oncocyclus to represent dramatically different genetic entities which devaluates a general trend of merging many previously recognized species of section Oncocyclus based on their flower morphology. Despite the importance of homoploid hybridization in this group’s evolution and some apparently sporadically happening inter-specific gene flow, the main evolutionary forces in Oncocyclus appear to be vicariance and spatial isolation. Our findings suggest that some of the currently recognized species in section Oncocyclus need revision. A revision must be based on genetic analyses allowing the reconstruction of ancestry and recognition of the importance of vicariance and spatial isolation in the evolution of this group. The implications of the present findings for conservation are discussed.

A snapshot of progenitor-derivative speciation in Iberodes (Boraginaceae)

Traditional classification of speciation modes has focused on physical barriers to gene flow. Allopatric speciation with complete reproductive isolation is viewed as the most common mechanism of speciation. Parapatry and sympatry, by contrast, entail speciation in the face of ongoing gene flow, making them more difficult to detect. The genus Iberodes (Boraginaceae, NW Europe) comprises five species with contrasting morphological traits, habitats, and species distributions. Based on the predominance of narrow and geographically distant endemic species, we hypothesized that geographic barriers were responsible for most speciation events in Iberodes. We undertook an integrative study including: (i) phylogenomics through restriction-site associated DNA sequencing, (ii) genetic structure analyses, (iii) demographic modeling, (iv) morphometrics, and (v) climatic niche modeling and niche overlap analysis. Results revealed a history of recurrent progenitor-derivative speciation manifested by a paraphyletic pattern of nested species differentiation. Budding speciation mediated by ecological differentiation is suggested for the coastal lineage, deriving from the inland widespread I. linifolia during Late Pliocene. Meanwhile, geographic isolation followed by niche shifts are suggested for the more recent differentiation of the coastland taxa. Our work provides a model for distinguishing speciation via ecological differentiation of peripheral, narrowly endemic I. kuzinskyanae and I. littoralis from a widespread extant ancestor, I. linifolia. Ultimately, our results illustrate a case of Pliocene speciation in the probable absence of geographic barriers and get away from the traditional cladistic perspective of speciation as producing two species from an extinct ancestor, thus reminding us that phylogenetic trees tell only part of the story.

Comparative Ecological Perspectives of Two Ancient Lineages of Gray Wolves: Woolly Wolf (Canis lupus chanco) and Indian Wolf (Canis lupus pallipes)

Geographical isolation can often lead to speciation, and two disconnected populations of the same species living in drastically different bioclimatic regions provide an opportunity to understand the process of speciation. The Woolly wolf is found in the cold-arid, Trans-Himalayan landscape, while the Indian wolf inhabits the semi-arid grasslands of Central India. Both the lineages of wolves from India have generated scientific debate on their taxonomic status in recent years. In this study, we collected data and reviewed published literature to document the ecological and behavioral differences between the Woolly wolf and the Indian wolf. Most studies have used genetic data; hence we discuss variation in spatial ecology, habitat preferences, vocalization, diet diversity and cranial measurements of these two subspecies. The spatial ecology of two lineages was compared from the data on three Woolly and ten Indian wolves tagged with GPS collars. The telemetry data shows that there has been no difference in the day-night movement of Woolly wolves, whereas Indian wolves show significant high displacement during the night. The BBMM method indicated that Woolly wolf home ranges were three times larger than the Indian wolf. The Woolly wolf diet is comprised of 20 different types of food items, whereas the Indian wolf diet consists of 17 types. The Woolly and Indian wolf largely depend upon domestic prey base, i.e., 48.44 and 40.34%, respectively. We found no differences in the howling parameters of these subspecies. Moreover, the Woolly wolf skull was significantly longer and broader than the Indian wolf. Wolves of India are ancient and diverged from the main clade about 200,000–1,000,000 years ago. Their genetic and ecological evolution in different bioclimatic zones has resulted in considerable differences as distinct subspecies. The present study is a step in understanding ecological differences between two important, genetically unique subspecies of wolves.

Reproductive Isolation in the Cryptic Species Complex of a Key Pest: Analysis of Mating and Rejection Behaviour of Onion Thrips (Thrips tabaci Lindeman)

Simple Summary

The onion thrips (Thrips tabaci), which is a key pest worldwide, includes three characteristic, distinct groups (i.e., lineages) under the same species name. In the current report, we addressed the question of whether individuals from these lineages recognize and assess each other as potential mating partners. We demonstrated that at least two of these lineages do not belong to the same species, since under our laboratory conditions no mating occurred between them. Moreover, specimens from these cross pairings often exhibited an escape response upon contact with the other thrips, while most of the pairs from the same lineages readily mated at their first interaction. The behaviour of males clearly indicated that they can assess the mating status of a female and usually only try to copulate with virgins. Our findings are important not only because in agriculture it is vital to know whether morphologically similar entities belong to the same species, but also because knowledge about the communication of insects and the possible role of the cues involved can help the development of new plant protection techniques.

Abstract

Thrips tabaci Lindeman is a serious pest of various cultivated plants, with three, distinct lineages within a cryptic species complex. Despite the well-known significance of this pest, many attributes of these lineages are not yet fully understood, including their reproductive behaviour. We performed no-choice-design cross-mating experiments under a controlled laboratory environment with virgin adult individuals from all three lineages. The behaviour of thrips was recorded with a camera mounted on a stereomicroscope, and the recordings were analysed in detail. We found that the so-called leek-associated lineages of this cryptic species complex are reproductively isolated from the tobacco-associated lineage; therefore, they represent different species. Divergence in the behaviour of conspecific and heterospecific pairs became evident only after contact. There were no marked differences between the lineages in their precopulatory and copulatory behaviour, except in the duration of the latter. We confirmed mating between thelytokous females and arrhenotokous males; however, we assume some form of loss of function in the sexual traits of asexual females. The post-mating behaviour of males indicated the presence and role of an anti-aphrodisiac pheromone. We also demonstrated differences between lineages regarding their activity and their propensity for exhibiting an escape response upon interaction with heterospecific thrips.

Extensive introgression at late stages of species formation: Insights from grasshopper hybrid zones

The process of species formation is characterised by the accumulation of multiple reproductive barriers. The evolution of hybrid male sterility, or Haldane's rule, typically characterises later stages of species formation, when reproductive isolation is strongest. Yet, understanding how quickly reproductive barriers evolve and their consequences for maintaining genetic boundaries between emerging species remains a challenging task because it requires studying taxa that hybridise in nature. Here, we address these questions using the meadow grasshopper Pseudochorthippus parallelus, where populations that show multiple reproductive barriers, including hybrid male sterility, hybridise in two natural hybrid zones. Using mitochondrial data, we infer that such populations have diverged some 100,000 years ago, at the beginning of the last glacial cycle in Europe. Nuclear data shows that contractions at multiple glacial refugia, and post-glacial expansions have facilitated genetic differentiation between lineages that today interact in hybrid zones. We find extensive introgression throughout the sampled species range, irrespective of current strength of reproductive isolation. Populations exhibiting hybrid male sterility in two hybrid zones show repeatable patterns of genomic differentiation, consistent with shared genomic constraints affecting ancestral divergence or with the role of those regions in reproductive isolation. Together, our results suggest that reproductive barriers that characterise late stages of species formation can evolve relatively quickly, particularly when associated with strong demographic changes. Moreover, we show that such barriers persist in the face of extensive gene flow, allowing future studies to identify associated genomic regions.

Reproductive character displacement allows two sexually deceptive orchids to coexist and attract the same specific pollinator

An increased divergence in characters between species in secondary contact can be shaped by selection against competition for a common resource (ecological character displacement, ECD) or against maladapted hybridization (reproductive character displacement, RCD). These selective pressures can act between incipient species (reinforcement) or well-separated species that already completed the speciation process, but that can still hybridize and produce maladapted hybrids. Here, we investigated two well-separated sexually deceptive orchid species that, unusually, share their specific pollinator. Sympatric individuals of these species are more divergent than allopatric ones in floral characters involved in a mechanical isolating barrier, a pattern suggestive of RCD. To experimentally test this scenario, we built an artificial sympatric population with allopatric individuals. We measured flower characters, genotyped the offspring in natural and artificial sympatry and estimated fertility of hybrids. Different from naturally sympatric individuals, allopatric individuals in artificial sympatry hybridized widely. Hybrids showed lower pollination success and seed viability than parentals. Character displacement did not affect plant pollination success. These findings suggest that RCD evolved between these species to avoid hybridization and that selection on reinforcement may be very strong even in plants with highly specialized pollination.

Ecological Opportunity and Necessity: Biotic and Abiotic Drivers Interact During Diversification of Digital Host-Parasite Communities

Most of Earth’s diversity has been produced in rounds of adaptive radiation, but the ecological drivers of diversification, such as abiotic complexity (i.e., ecological opportunity) or predation and parasitism (i.e., ecological necessity), are hard to disentangle. However, most of these radiations occurred hundreds of thousands if not millions of years ago, and the mechanisms promoting contemporary coexistence are not necessarily the same mechanisms that drove diversification in the first place. Experimental evolution has been one fruitful approach used to understand how different ecological mechanisms promote diversification in simple microbial microcosms, but these microbial systems come with their own limitations. To test how ecological necessity and opportunity interact, we use an unusual system of self-replicating computer programs that diversify to fill niches in a virtual environment. These organisms are subject to ecological pressures just like their natural counterparts. They experience biotic interactions from digital parasites, which steal host resources to replicate their own code and spread in the population. With the control afforded by experimenting with computational ecologies, we begin to unweave the complex interplay between ecological drivers of diversification. In particular, we find that the complexity of the abiotic environment and the size of the phenotypic space in which organisms are able to interact play different roles depending on the ecological driver of diversification. We find that in some situations, both ecological opportunity and necessity drive similar levels of diversity. However, the phenotypes that hosts uncover while coevolving with parasites are dramatically more complex than hosts evolving alone.

The Build-Up of Population Genetic Divergence along the Speciation Continuum during a Recent Adaptive Radiation of Rhagoletis Flies

New species form through the evolution of genetic barriers to gene flow between previously interbreeding populations. The understanding of how speciation proceeds is hampered by our inability to follow cases of incipient speciation through time. Comparative approaches examining different diverging taxa may offer limited inferences, unless they fulfill criteria that make the comparisons relevant. Here, we test for those criteria in a recent adaptive radiation of the Rhagoletis pomonella species group (RPSG) hypothesized to have diverged in sympatry via adaptation to different host fruits. We use a large-scale population genetic survey of 1568 flies across 33 populations to: (1) detect on-going hybridization, (2) determine whether the RPSG is derived from the same proximate ancestor, and (3) examine patterns of clustering and differentiation among sympatric populations. We find that divergence of each in-group RPSG taxon is occurring under current gene flow, that the derived members are nested within the large pool of genetic variation present in hawthorn-infesting populations of R. pomonella, and that sympatric population pairs differ markedly in their degree of genotypic clustering and differentiation across loci. We conclude that the RPSG provides a particularly robust opportunity to make direct comparisons to test hypotheses about how ecological speciation proceeds despite on-going gene flow.

Geological and climatic influences on population differentiation of the Phrynocephalus vlangalii species complex (Sauria: Agamidae) in the northern Qinghai-Tibet Plateau

Extremely heterogeneous topography and complex paleoclimate history of the Qinghai-Tibet Plateau (QTP) have a key role in promoting genetic divergence and lineage/species formation. Here, we sequenced one nuclear and three mitochondrial markers of 532 individuals from the entire range of the Phrynocephalus vlangalii species complex including two species, P. putjatai and P. vlangalii, endemic to the northern QTP. We integrated multilocus phylogeny, demographic analysis and geographic barrier detection to evaluate the population structure and dynamics. We found a new mitochondrial clade (PV-I) in the Gonghe County population of P. vlangalii, partial mitochondrial DNA replacement within P. vlangalii and complete mitochondrial DNA replacement between P. putjatai and P. vlangalii. Neutrality test, mismatch distribution analysis and Extended Bayesian Skyline Plot (EBSP) analysis all supported a significant expansion of the Qaidam Basin population of P. vlangalii (PV-II-2) from 0.091 to 0.026 Ma after Penultimate Glaciation. The uplift of the Arjin and Anyemanqen Mountains during the Kunhuang Movement (∼1.2 Ma) split populations of P. vlangalii in Akesai, Qaidam Basin and source of the Yellow River. The uplift of the Elashan Mountains during the second phase of the Qingzang Movement (∼2.5 Ma) contributed to the divergence of the Gonghe County population of P. vlangalii from other conspecific populations. The third phase of the Qingzang Movement (∼1.7 Ma) contributed to the divergence of the Xinghai population of P. vlangalii from P. putjatai and to the divergence of the northern populations of P. putjatai from the southern conspecific populations. Our data support the idea that the geological and climatic changes following the orogeny of the QTP may have promoted population differentiation and shaped the current population patterns of the P. vlangalii species complex in the northeastern QTP.

Misinterpretation of Genomic Data Matters for Endangered Species Listing: The Sub-specific Status of the Peñasco Least Chipmunk (Neotamias minimus atristriatus)

We provide a response to a recently published evaluation of the subspecies status of the Peñasco least chipmunk (Neotamias minimus atristriatus). The work we discuss used exon capture genomic approaches and concluded that their results did not support the distinction of this taxon as a subspecies, with recommendation that it be synonymized with N. m. operarius. We refute the interpretations, conclusions, and taxonomic recommendations of this study, and explain in clearer terms how to interpret genomic analyses for applied management. We identify four broad conceptual issues that led to errant recommendations: (1) interpretation of subspecies and diagnosability, (2) inappropriate use of reciprocal monophyly as a criterion for subspecies, (3) importance of geographic isolation, and (4) error in hypothesis testing and misinterpretation of results. We conclude that the data from this genomic appraisal add to information from prior studies providing strong support for recognition of N. m. atristriatus as a subspecies. Our conclusions have important and immediate implications for the proposed listing of N. m. atristriatus as an endangered species under the U.S. Endangered Species Act.