Within-island speciation with an exceptional case of distinct separation between two sibling lizard species divided by a narrow stream

Range-Wide Phylogeography and Ecological Niche Modeling Provide Insights into the Evolutionary History of the Mongolian Racerunner (Eremias argus) in Northeast Asia

The Mongolian racerunner, Eremias argus, is a small lizard endemic to Northeast Asia that can serve as an excellent model for investigating how geography and past climate change have jointly influenced the evolution of biodiversity in this region. To elucidate the processes underlying its diversification and demography, we reconstructed the range-wide phylogeographic pattern and evolutionary trajectory, using phylogenetic, population genetic, landscape genetic, Bayesian phylogeographic reconstruction and ecological niche modeling approaches. Phylogenetic analyses of the mtDNA cyt b gene revealed eight lineages that were unbounded by geographic region. The genetic structure of E. argus was mainly determined by geographic distance. Divergence dating indicated that E. argus and E. brenchleyi diverged during the Mid-Pliocene Warm Period. E. argus was estimated to have coalesced at~0.4351 Ma (Marine Isotope Stage 19). Bayesian phylogeographic diffusion analysis revealed out-of-Inner Mongolia and rapid colonization events from the end of the Last Interglacial to the Last Glacial Maximum, which is consistent with the expanded suitable range of the Last Glacial Maximum. Pre-Last Glacial Maximum growth of population is presented for most lineages of E. argus. The Glacial Maximum contraction model and the previous multiple glacial refugia hypotheses are rejected. This may be due to an increase in the amount of climatically favorable habitats in Northeast Asia. Furthermore, E. argus barbouri most likely represents an invalid taxon. The present study is the first to report a range-wide phylogeography of reptiles over such a large region in Northeast Asia. Our results make a significant contribution towards understanding the biogeography of the entire Northeast Asia.

Integrative species delimitation reveals fine-scale allopatric speciation in a good-flying insect: a case study on

Alpha taxonomy is fundamental for many biological fields. Delineation of the species boundary, however, can be challenging in a species complex, where different species share a similar morphology and diagnostic characters may not be available. In this context, integrative approaches that incorporate molecular and morphological data sets, and account for speciation history can be helpful to alpha taxonomy. Different approaches to species delimitation based on different assumptions are complementary and by integrating the results from multiple approaches we can generate a more reliable and objective taxonomic decision. In this study, we applied three molecular approaches to species delimitation and inferred the demographic history based on an isolation with migration model to test a morphologically based taxonomic hypothesis for the Cylindera pseudocylindriformis complex. We discuss the association between genetic divergence and microhabitat specialisation, and further corroborate that C. subtilis sp. nov. is a valid new species by integrating the results from model-based species delimitation and the genealogical divergence index. We argue that genetic endemism can occur at a small geographic scale, even in a winged insect like tiger beetles. Our results also indicated that there may still be undocumented species diversity of Taiwanese Cylindera remaining to be discovered. ZooBank LSID: urn:lsid:zoobank.org:pub:9DEC1432-365C-4872-8D06-73B95F30624F

Incipient speciation and its impact on taxonomic decision: a case study using a sky island sister-species pair of stag beetles (Lucanidae: Lucanus)

Species delimitation can be difficult when the divergence between focal taxa is in the incipient stage of speciation, because conflicting results are expected among different data sets, and the species limits can differ depending on the species concept applied. We studied speciation history and investigated the impact on taxonomic decision-making when using different types of data in a Taiwanese endemic sister-species pair of stag beetles, Lucanus miwai and Lucanus yulaoensis, from sky island habitats. We showed that the two geographical taxa can be diagnosed by male mandibular shape. We found two mitochondrial co1 lineages with pairwise sequence divergence > 3%; however, L. miwai might not be monophyletic. The result of our multispecies coalescent-based species delimitation using five nuclear loci supported the evolutionary independence of the two sister species, but the calculated values of the genealogical divergence index (gdi) corresponded to the ambiguous zone of species delimitation. We also showed that post-divergence gene flow is unlikely. Our study demonstrates challenges in the delineation of incipient species, but shows the importance of understanding the speciation history and adopting integrative approaches to reconcile seemingly conflicting results before making evolutionarily relevant taxonomic decisions.

Asymmetric acoustic signal recognition led to asymmetric gene flow between two parapatric frogs

Correct discrimination between courtship signals could help to maintain genetic integrity between closely related species. However, asymmetric usage of signals might cause asymmetric gene flow across the contact zone. Buergeria choui and B. otai are sibling-species with a parapatric distribution pattern in Taiwan, having two narrow contact zones on the east and west sides of the island. Combining behavioural experiments with genome-wide RAD-seq analyses, we test whether the ability of signal recognition influences genetic introgression across their species boundary. The playback experiments show that all B. choui populations respond strongest to their own ‘cricket’ trills, while the western population of B. otai have evolved a strong level of reproductive character displacement by showing the inclusive usage of the unique ‘chicken’ signals. In contrast, the eastern B. otai population uses both ‘chicken’ and ‘cricket’ trills, and has a stronger preference for the latter. The weak reproductive character displacement in the eastern population has led to asymmetry genetic introgression from B. choui toward B. otai. Our results support the prediction that a more specialized signal-user, compared to its sibling, generalized signal-user, might have a higher probability of maintaining their genetic integrity in the secondary contact region.

Ancestral Haplotype Retention and Population Expansion Determine the Complicated Population Genetic Structure of the Hilly Lineage of Neolucanus swinhoei Complex (Coleoptera, Lucanidae) on the Subtropical Taiwan Island

The present study demonstrates that the complicated genetic structure of the hilly lineage of the Neolucanus swinhoei complex was driven by its biological features and habitat requirements as well as hindrance by the CMR during periodical Pleistocene glaciations. The results revealed a tendency of geographical differentiation and major and sub- lineage divergences before and after the Riss glaciation, followed by stable population growth during Würm glaciation. At least four refugia were inferred for N. swinhoei during the Riss-Würm glaciations. The ancestral haplotype retention in the cytochrome oxidase subunit I (COI) gene and compensated substitution in 16S rRNA gene is a possible evolutionary scenario resulting in the inconsistent evolution pattern between COI and 16S rRNA gene coupled with the long-distance dispersal of N. swinhoei. Although the CMR did hinder the dispersal of N. swinhoei, its ancestors may have dispersed to eastern Taiwan through the northern and southern low mountains of the CMR before the Riss glaciation. Our finding suggests that the population growth in the Würm glaciation led a dispersal back to western Taiwan, which is contrast to the more common dispersal scenario from western Taiwan to eastern populations proposed in other studies.

Long-term monitoring reveals invariant clutch size and unequal reproductive costs between sexes in a subtropical lacertid lizard

Based on 20,000 records representing c. 11,000 individuals from an 8-year capture-mark-recapture (CMR) study, we tested and confirmed a new case of invariant clutch size (ICS) in a sexually dichromatic lacertid lizard, Takydromus viridipunctatus. In the grassland habitat of the early succession stage, females showed strictly low and invariant clutch size, multiple clutches in a breeding season, high reproductive potential, and annual breeding cycles that correspond to the emergence of male courtship coloration. The hatchlings mature quickly, and join the adult cohort for breeding within a few months, whereas adults show low survival rates and a short lifespan, such that most die within one year. Mortality increased in both sexes during the breeding season, especially in females, indicating an unequal cost of reproduction in survival. These life history characters may be explained by two non-exclusive hypotheses of ICS-arboreal hypothesis and predation hypothesis-within the ecological context of their habitat. Our study highlights a confirmed case of ICS, which adapts well to this r-selected grassland habitat that experiences seasonal fluctuation and frequent disturbance.

A check list and population trends of invasive amphibians and reptiles in Taiwan

Invasive species have impacted biodiversity all around the world. Among various ecosystems, islands are most vulnerable to these impacts due to their high ratio of endemism, highly specialized adaptation, and isolated and unique fauna. As with other subtropical islands, Taiwan faces constant risk of biological invasions and is currently ranked as one of the countries most affected by invasive amphibians and reptiles. In this paper, a comprehensive checklist of all known exotic amphibians and reptiles is provided, including twelve species which have successfully colonized Taiwan and six species with a controversial status. We provide an update on the knowledge of all these species including their distribution, colonization history, threats to native animals, and population trends based on literature records, fauna surveys, and data collected during invasive species eradication and control programs. A list of species with high invasive potentials is also provided. This study reports, for the first time, a comprehensive survey of invasive herpetofauna in Taiwan, which should provide a valuable reference to other regions which might suffer from similar invasion risk.

Differential genetic responses to the stress revealed the mutation-order adaptive divergence between two sympatric ginger species

BACKGROUND

Divergent genetic responses to the same environmental pressures may lead sympatric ecological speciation possible. Such speciation process possibly explains rapid sympatric speciation of island species. Two island endemic ginger species Zingiber kawagoii and Z. shuanglongensis was suggested to be independently originated from inland ancestors, but their island endemism and similar morphologies and habitats lead another hypothesis of in situ ecological speciation. For understanding when and how these two species diverged, intraspecific variation was estimated from three chloroplast DNA fragments (cpDNA) and interspecific genome-wide SNPs and expression differences after saline treatment were examined by transcriptomic analyses.

RESULTS

Extremely low intraspecific genetic variation was estimated by cpDNA sequences in both species: nucleotide diversity π = 0.00002 in Z. kawagoii and no nucleotide substitution but only indels found in Z. shuanglongensis. Nonsignificant inter-population genetic differentiation suggests homogenized genetic variation within species. Based on 53,683 SNPs from 13,842 polymorphic transcripts, in which 10,693 SNPs are fixed between species, Z. kawagoii and Z. shuanglongensis were estimated to be diverged since 218~ 238 thousand generations ago (complete divergence since 41.5~ 43.5 thousand generations ago). This time is more recent than the time of Taiwan Island formation. In addition, high proportion of differential expression genes (DEGs) is non-polymorphic or non-positively selected, suggesting key roles of plastic genetic divergence in broaden the selectability in incipient speciation. While some positive selected DEGs were mainly the biotic and abiotic stress-resistance genes, emphasizing the importance of adaptive divergence of stress-related genes in sympatric ecological speciation. Furthermore, the higher proportional expression of functional classes in Z. kawagoii than in Z. shuanglongensis explains the more widespread distribution of Z. kawagoii in Taiwan.

CONCLUSIONS

Our results contradict the previous hypothesis of independent origination of these two island endemic ginger species from SE China and SW China. Adaptive divergent responses to the stress explain how these gingers maintain genetic differentiation in sympatry. However, the recent speciation and rapid expansion make extremely low intraspecific genetic variation in these two species. This study arise a more probable speciation hypothesis of sympatric speciation within an island via the mutation-order mechanism underlying the same environmental pressure.

Geological data indicate that the interpretation for the age-calibrated phylogeny for the Kurixalus-genus frogs of South, South-east and East Asia (Lv et al., 2018) needs to be rethought

Recently, Lv et al. (2018) published an age-calibrated phylogenetic tree for the Kurixalus frogs, members of which occur across parts of South, South-east and East Asia. A clade on Taiwan, represented by Kurixalus idiootocus and the Kurixalus eiffingeri species complex, is deemed to have been resident since the middle Cenozoic; its closest congeners are in southern Indochina (not in the adjacent parts of south-east China), and the split between the two is dated at 32.8 Ma. Furthermore, a sub-population of Kurixalus eiffingeri is believed to have colonized islands in the western Ryukyus c. 13.5 Ma. There is, however, a problem with this scenario: the landmass regarded as modern-day Taiwan has existed only for 4-5 million years (it results from a young and ongoing tectonic-plate collision). Assuming the Kurixalus phylogeny and the dating of its branchings are correct, then a palaeobiogeographical scenario involving an older, alternative land surface with later transfer to Taiwan, possibly involving over-water dispersal, would reconcile the biology, but testing this may be difficult/impossible. If the ages of the nodes in the proposed tree are found to be significantly overestimated, the geology and biology might more easily be accommodated.

Taxonomic reassessment of two subspecies of Chinese skink in Taiwan based on morphological and molecular investigations (Squamata, Scincidae)

The Chinese skink, Plestiodonchinensis (Gray, 1838), is widely distributed across continental China, Taiwan, the Korean Peninsula, and offshore islets, and consists of several subspecies. Here morphological and molecular methods have been used to reassess the taxonomic status and distributions of P.c.formosensis (Van Denburgh, 1912) and P.c.leucostictus (Hikida, 1988), which are endemic to Taiwan and Green Island (an islet off the east coast of Taiwan), respectively. It can be confirmed that the eastern Taiwanese populations of P.c.formosensis exhibit similar juvenile color patterning and genetic composition to the islet subspecies P.c.leucostictus, and are distinct from consubspecific populations in western Taiwan. Therefore, the eastern Taiwanese populations are assigned to P.c.leucostictus, and this subspecies is recognized as a distinct species, Plestiodonleucostictus (Hikida, 1988), based on their unique juvenile coloration and highly divergent DNA sequences. Our results also revealed that P.c.formosensis in western Taiwan is close to nominotypical subspecies from the continent, suggesting the necessity of a comprehensive taxonomic analysis in the future.

Subspecific Differentiation Events of Montane Stag Beetles (Coleoptera, Lucanidae) Endemic to Formosa Island

Taxonomic debates have been carrying on for decades over Formosan stag beetles, which consist of a high proportion of endemic species and subspecies featuring morphological variations associated with local adaptation. With the influence of periodical Pleistocene glaciations and the presence of several mountain ranges, the genetic differentiation and taxonomic recognition, within this medium-size island, of two endemic subspecies for each of four montane stag beetles, i.e. Lucanus ogakii, L. kanoi, Prismognathus davidis, and Neolucanus doro, has been an appealing issue. Based on monophyletic lineages and population structure, possible divergent scenarios have been proposed to clarify the subspecific status for each of the above mentioned stag beetles. Phylogenetic inferences based on COI+16S rDNA+28S rDNA of 240 Formosan lucanids have confirmed most species are monophyletic groups; and the intraspecific (<2%) and interspecific (>2%) genetic distances of the two mitochondrial genes could be applied concordantly for taxonomic identification. On account of Bayesian-based species delimitation, geographic distribution, population structure, and sequence divergences, the subspecific status for L. ogakii, L. kanoi, and Pri. davidis are congruent with their geographic distribution in this island; and the calibration time based on the mitochondrial genes shows the subspecific split events occurred 0.7–1 million years ago. In addition, a more complicated scenario, i.e. genetic differentiation including introgression/hybridization events, might have occurred among L. ogakii, L. kanoi, and L. maculifemoratus. The geological effects of mountain hindrance accompanied by periodical glaciations could have been vital in leading to the geographical subspecific differentiation of these montane stag beetles.

Geographic variation in advertisement calls of a Microhylid frog - testing the role of drift and ecology

Acoustic signals for mating are important traits that could drive population differentiation and speciation. Ecology may play a role in acoustic divergence through direct selection (e.g., local adaptation to abiotic environment), constraint of correlated traits (e.g., acoustic traits linked to another trait under selection), and/or interspecific competition (e.g., character displacement). However, genetic drift alone can also drive acoustic divergence. It is not always easy to differentiate the role of ecology versus drift in acoustic divergence. In this study, we tested the role of ecology and drift in shaping geographic variation in the advertisement calls of Microhyla fissipes. We examined three predictions based on ecological processes: (1) the correlation between temperature and call properties across M. fissipes populations; (2) the correlation between call properties and body size across M. fissipes populations; and (3) reproductive character displacement (RCD) in call properties between M. fissipes populations that are sympatric with and allopatric to a congener M. heymonsi. To test genetic drift, we examined correlations among call divergence, geographic distance, and genetic distance across M. fissipes populations. We recorded the advertisement calls from 11 populations of M. fissipes in Taiwan, five of which are sympatrically distributed with M. heymonsi. We found geographic variation in both temporal and spectral properties of the advertisement calls of M. fissipes. However, the call properties were not correlated with local temperature or the callers' body size. Furthermore, we did not detect RCD. By contrast, call divergence, geographic distance, and genetic distance between M. fissipes populations were all positively correlated. The comparisons between phenotypic Q_st (P_st) and F_st values did not show significant differences, suggesting a role of drift. We concluded that genetic drift, rather than ecological processes, is the more likely driver for the geographic variation in the advertisement calls of M. fissipes.