Pleistocene climatic cycling drives intra-specific diversification in the intermediate horseshoe bat (Rhinolophus affinis) in Southern China

Phylogeographic Pattern of the Assassin Bug Sycanus bifidus Inferred from Mitochondrial Genomes and Nuclear Genes

The assassin bug Sycanus bifidus has a wide distribution across southern China. This study explored its distribution and evolution by analyzing mitochondrial and nuclear ribosomal RNA genes, revealing how Pleistocene climate and geological changes shaped its phylogeography. We identified two main clades, A and B, that diverged in the Middle Pleistocene. Hainan Island's populations form a unique group within Clade A, suggesting that the Qiongzhou Strait served as a dispersal corridor during glaciation. Rising sea levels likely separated the Hainan population afterward. Ecological niche modeling showed that both populations have been viable since the last interglacial period, with demographic analyses indicating possible expansions during the Middle and Late Pleistocene, driven by favorable climates. This study highlights the significant effects of Pleistocene sea-level and climatic changes on the distribution and evolution of S. bifidus in China.

Taxonomic status of the Lesser Woolly Horseshoe bats (Chiroptera, Rhinolophidae, Rhinolophus beddomei) in peninsular India and Sri Lanka

Insularity provides ample opportunities for species diversification. Sri Lanka is home to a large diversity of species, many of which are endemic but morphologically similar to species found in southern India, due to recent speciation events, suggesting a complex evolutionary history. However, in some taxa although morphological diversity has been noted, the genetic level variations are minimal. Among the wide-ranging horseshoe bats such a phenomenon is noted. In this study, we used bioacoustics, morphometric and molecular data to evaluate the relationships between the taxa of lesser woolly horseshoe bats in the India and Sri Lanka. Our study reveals that the two taxa-Rhinolophus beddomei Andersen, 1905 and here we have validated the existing subspecies from peninsular India and R. sobrinus Andersen, 1918 from Sri Lanka are genetically very close to R. perniger Hodgson, 1843. Currently the taxa-beddomei and sobrinus are recognized as subspecies of Rhinolophus beddomei Andersen, 1905. We provide a detailed description of the taxa beddomei and sobrinus as the original descriptions are limited in nature.

Global patterns of phylogenetic diversity and transmission of bat coronavirus

Bats are reservoirs for multiple coronaviruses (CoVs). However, the phylogenetic diversity and transmission of global bat-borne CoVs remain poorly understood. Here, we performed a Bayesian phylogeographic analysis based on 3,594 bat CoV RdRp gene sequences to study the phylogenetic diversity and transmission of bat-borne CoVs and the underlying driving factors. We found that host-switching events occurred more frequently for α-CoVs than for β-CoVs, and the latter was highly constrained by bat phylogeny. Bat species in the families Molossidae, Rhinolophidae, Miniopteridae, and Vespertilionidae had larger contributions to the cross-species transmission of bat CoVs. Regions of eastern and southern Africa, southern South America, Western Europe, and Southeast Asia were more frequently involved in cross-region transmission events of bat CoVs than other regions. Phylogenetic and geographic distances were the most important factors limiting CoV transmission. Bat taxa and global geographic hotspots associated with bat CoV phylogenetic diversity were identified, and bat species richness, mean annual temperature, global agricultural cropland, and human population density were strongly correlated with the phylogenetic diversity of bat CoVs. These findings provide insight into bat CoV evolution and ecological transmission among bat taxa. The identified hotspots of bat CoV evolution and transmission will guide early warnings of bat-borne CoV zoonotic diseases.

The evolutionary history and ancestral biogeographic range estimation of old-world Rhinolophidae and Hipposideridae (Chiroptera)

BACKGROUND

Family Rhinolophidae (horseshoe bats), Hipposideridae (leaf-nosed bats) and Rhinonycteridae (trident bats) are exclusively distributed in the Old-World, and their biogeography reflects the complex historic geological events throughout the Cenozoic. Here we investigated the origin of these families and unravel the conflicting family origin theories using a high resolution tree covering taxa from each zoogeographic realm from Africa to Australia. Ancestral range estimations were performed using a probabilistic approach implemented in BioGeoBEARS with subset analysis per biogeographic range [Old-World as whole, Australia-Oriental-Oceania (AOO) and Afrotropical-Madagascar-Palearctic (AMP)].

RESULT

Our result supports an Oriental origin for Rhinolophidae, whereas Hipposideridae originated from the Oriental and African regions in concordance with fossil evidence of both families. The fossil evidence indicates that Hipposideridae has diversified across Eurasia and the Afro-Arabian region since the Middle Eocene. Meanwhile, Rhinonycteridae (the sister family of Hipposideridae) appears to have originated from the Africa region splitting from the common ancestor with Hipposideridae in Africa. Indomalaya is the center of origin of Rhinolophidae AOO lineages, and Indomalayan + Philippines appears to be center of origin of Hipposideridae AOO lineage indicating allopatric speciation and may have involved jump-dispersal (founder-event) speciation within AOO lineage. Wallacea and the Philippines may have been used as stepping stones for dispersal towards Oceania and Australia from the Oriental region. Multiple colonization events via different routes may have occurred in the Philippines (i.e., Palawan and Wallacea) since the Late Miocene. The colonization of Rhinolophidae towards Africa from Asia coincided with the estimated time of Tethys Ocean closure around the Oligocene to Miocene (around 27 Ma), allowing species to disperse via the Arabian Peninsula. Additionally, the number of potential cryptic species in Rhinolophidae in Southeast Asia may have increased since Plio-Pleistocene and late Miocene.

CONCLUSION

Overall, we conclude an Oriental origin for Rhinolophidae, and Oriental + African for Hipposideridae. The result demonstrates that complex historical events, in addition to species specific ecomorphology and specialization of ecological niches may shape current distributions.

How to Accurately Delineate Morphologically Conserved Taxa and Diagnose Their Phenotypic Disparities: Species Delimitation in Cryptic Rhinolophidae (Chiroptera)

Systematics and taxonomy are the backbone of all components of biology and ecology, yet cryptic species present a major challenge for accurate species identification. This is especially problematic as they represent a substantial portion of undiscovered biodiversity, and have implications for not only species conservation, but even assaying potential risk of zoonotic spillover. Here, we use integrative approaches to delineate potential cryptic species in horseshoe bats (Rhinolophidae), evaluate the phenotypic disparities between cryptic species, and identify key traits for their identification. We tested the use of multispecies coalescent models (MSC) using Bayesian Phylogenetic and Phylogeography (BPP) and found that BPP was useful in delineating potential cryptic species, and consistent with acoustic traits. Our results show that around 40% of Asian rhinolophid species are potentially cryptic and have not been formally described. In order to avoid potential misidentification and allow species to be accurately identified, we identified quantitative noseleaf sella and acoustic characters as the most informative traits in delineating between potential cryptic species in Rhinolophidae. This highlights the physical differences between cryptic species that are apparent in noseleaf traits which often only qualitatively described but rarely measured. Each part of the noseleaf including the sella, lateral lappets, and lancet furrows, play roles in focusing acoustic beams and thus, provide useful characteristics to identify cryptic Rhinolophus species. Finally, species delimitation for cryptic species cannot rely on genetic data alone, but such data should be complemented by other evidence, including phenotypic, acoustic data, and geographic distributions to ensure accurate species identification and delineation.

Exploring the Natural Origins of SARS-CoV-2 in the Light of Recombination

The lack of an identifiable intermediate host species for the proximal animal ancestor of SARS-CoV-2, and the large geographical distance between Wuhan and where the closest evolutionary related coronaviruses circulating in horseshoe bats (members of the Sarbecovirus subgenus) have been identified, is fueling speculation on the natural origins of SARS-CoV-2. We performed a comprehensive phylogenetic study on SARS-CoV-2 and all the related bat and pangolin sarbecoviruses sampled so far. Determining the likely recombination events reveals a highly reticulate evolutionary history within this group of coronaviruses. Distribution of the inferred recombination events is nonrandom with evidence that Spike, the main target for humoral immunity, is beside a recombination hotspot likely driving antigenic shift events in the ancestry of bat sarbecoviruses. Coupled with the geographic ranges of their hosts and the sampling locations, across southern China, and into Southeast Asia, we confirm that horseshoe bats, Rhinolophus, are the likely reservoir species for the SARS-CoV-2 progenitor. By tracing the recombinant sequence patterns, we conclude that there has been relatively recent geographic movement and cocirculation of these viruses' ancestors, extending across their bat host ranges in China and Southeast Asia over the last 100 years. We confirm that a direct proximal ancestor to SARS-CoV-2 has not yet been sampled, since the closest known relatives collected in Yunnan shared a common ancestor with SARS-CoV-2 approximately 40 years ago. Our analysis highlights the need for dramatically more wildlife sampling to: 1) pinpoint the exact origins of SARS-CoV-2's animal progenitor, 2) the intermediate species that facilitated transmission from bats to humans (if there is one), and 3) survey the extent of the diversity in the related sarbecoviruses' phylogeny that present high risk for future spillovers.

Mitonuclear mismatch alters nuclear gene expression in naturally introgressed Rhinolophus bats

Background

Mitochondrial function involves the interplay between mitochondrial and nuclear genomes. Such mitonuclear interactions can be disrupted by the introgression of mitochondrial DNA between taxa or divergent populations. Previous studies of several model systems (e.g. Drosophila) indicate that the disruption of mitonuclear interactions, termed mitonuclear mismatch, can alter nuclear gene expression, yet few studies have focused on natural populations.

Results

Here we study a naturally introgressed population in the secondary contact zone of two subspecies of the intermediate horseshoe bat (Rhinolophus affinis), in which individuals possess either mitonuclear matched or mismatched genotypes. We generated transcriptome data for six tissue types from five mitonuclear matched and five mismatched individuals. Our results revealed strong tissue-specific effects of mitonuclear mismatch on nuclear gene expression with the largest effect seen in pectoral muscle. Moreover, consistent with the hypothesis that genes associated with the response to oxidative stress may be upregulated in mitonuclear mismatched individuals, we identified several such gene candidates, including DNASE1L3, GPx3 and HSPB6 in muscle, and ISG15 and IFI6 in heart.

Conclusion

Our study reveals how mitonuclear mismatch arising from introgression in natural populations is likely to have fitness consequences. Underlying the processes that maintain mitonuclear discordance is a step forward to understand the role of mitonuclear interactions in population divergence and speciation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00424-x.

Extensive alternative splicing triggered by mitonuclear mismatch in naturally introgressed Rhinolophus bats

Mitochondrial function needs strong interactions of mitochondrial and nuclear (mitonuclear) genomes, which can be disrupted by mitonuclear mismatch due to mitochondrial DNA (mtDNA) introgression between two formerly isolated populations or taxa. This mitonuclear disruption may cause severe cellular stress in mismatched individuals. Gene expression changes and alternative splicing (AS) are two important transcriptional regulations to respond to environmental or cellular stresses. We previously identified a naturally introgressed population in the intermediate horseshoe bat (Rhinolophus affinis). Individuals from this population belong to R. a. himalayanus and share almost identical nuclear genetic background; however, some of them had mtDNA from another subspecies (R. a. macrurus). With this unique natural system, we examined gene expression changes in six tissues between five mitonuclear mismatched and five matched individuals. A small number of differentially expressed genes (DEGs) were identified, and functional enrichment analysis revealed that most DEGs were related to immune response although some may be involved in response to oxidative stress. In contrast, we identified extensive AS events and alternatively spliced genes (ASGs) between mismatched and matched individuals. Functional enrichment analysis revealed that multiple ASGs were directly or indirectly associated with energy production in mitochondria which is vital for survival of organism. To our knowledge, this is the first study to examine the role of AS in responding to cellular stress caused by mitonuclear mismatch in natural populations. Our results suggest that AS may play a more important role than gene expression regulation in responding to severe environmental or cellular stresses.

The complete mitochondrial genome of Rhinolophus affinis himalayanus

Here we generated the complete mitochondrial genome of one subspecies of R. affinis (R. affinis himalayanus) using next generation sequencing and Sanger sequencing. The length of the complete mitochondrial genome was 16,886 bp, containing 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a non-coding control region. A maximum-likelihood tree based on the 13 concatenated mitochondrial protein-coding genes of 16 Rhinolophus taxon and one outgroup Hipposideros armiger indicates that R. affinis shows a closer relationship with R. sinicus complex than with other taxa.

Echolocation call frequency variation in horseshoe bats: molecular basis revealed by comparative transcriptomics

Recently diverged taxa with contrasting phenotypes offer opportunities for unravelling the genetic basis of phenotypic variation in nature. Horseshoe bats are a speciose group that exhibit a derived form of high-duty cycle echolocation in which the inner ear is finely tuned to echoes of the narrowband call frequency. Here, by focusing on three recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis) that display divergent echolocation call frequencies, we aim to identify candidate loci putatively involved in hearing frequency variation. We used de novo transcriptome sequencing of two mainland taxa (himalayanus and macrurus) and one island taxon (hainanus) to compare expression profiles of thousands of genes. By comparing taxa with divergent call frequencies (around 15 kHz difference), we identified 252 differentially expressed genes, of which six have been shown to be involved in hearing or deafness in human/mouse. To obtain further validation of these results, we applied quantitative reverse transcription-PCR to the candidate gene FBXL15 and found a broad association between the level of expression and call frequency across taxa. The genes identified here represent strong candidate loci associated with hearing frequency variation in bats.

Origin and cross-species transmission of bat coronaviruses in China

Bats are presumed reservoirs of diverse coronaviruses (CoVs) including progenitors of Severe Acute Respiratory Syndrome (SARS)-CoV and SARS-CoV-2, the causative agent of COVID-19. However, the evolution and diversification of these coronaviruses remains poorly understood. Here we use a Bayesian statistical framework and a large sequence data set from bat-CoVs (including 630 novel CoV sequences) in China to study their macroevolution, cross-species transmission and dispersal. We find that host-switching occurs more frequently and across more distantly related host taxa in alpha- than beta-CoVs, and is more highly constrained by phylogenetic distance for beta-CoVs. We show that inter-family and -genus switching is most common in Rhinolophidae and the genus Rhinolophus. Our analyses identify the host taxa and geographic regions that define hotspots of CoV evolutionary diversity in China that could help target bat-CoV discovery for proactive zoonotic disease surveillance. Finally, we present a phylogenetic analysis suggesting a likely origin for SARS-CoV-2 in Rhinolophus spp. bats.

Genome-wide data reveal discordant mitonuclear introgression in the intermediate horseshoe bat (Rhinolophus affinis)

Closely related taxa often exhibit mitonuclear discordance attributed to introgression of mitochondrial DNA (mtDNA), yet few studies have considered the underlying causes of mtDNA introgression. Here we test for demographic versus adaptive processes as explanations for mtDNA introgression in three subspecies of the intermediate horseshoe bat (Rhinolophus affinis). We generated sequences of 1692 nuclear genes and 13 mitochondrial protein-coding genes for 48 individuals. Phylogenetic reconstructions based on 320 exon sequences and 2217 single nucleotide polymorphisms (SNPs) both revealed conflicts between the species tree and mtDNA tree. These results, together with geographic patterns of mitonuclear discordance, and shared identical or near-identical mtDNA sequences, suggest extensive introgression of mtDNA between the two parapatric mainland subspecies. Under demographic hypotheses, we would also expect to uncover traces of ncDNA introgression, however, population structure and gene flow analyses revealed little nuclear admixture. Furthermore, we found inconsistent estimates of the timing of population expansion and that of the most recent common ancestor for the clade containing introgressed haplotypes. Without a clear demographic explanation, we also examined whether introgression likely arises from adaptation. We found that five mtDNA genes contained fixed amino acid differences between introgressed and non-introgressed individuals, including putative positive selection found in one codon, although this did not show introgression. While our evidence for rejecting demographic hypotheses is arguably stronger than that for rejecting adaptation, we find no definitive support for either explanation. Future efforts will focus on larger-scale resequencing to decipher the underlying causes of discordant mitonuclear introgression in this system.

Origin and cross-species transmission of bat coronaviruses in China

Bats are presumed reservoirs of diverse coronaviruses (CoVs) including progenitors of Severe Acute Respiratory Syndrome (SARS)-CoV and SARS-CoV-2, the causative agent of COVID-19. However, the evolution and diversification of these coronaviruses remains poorly understood. We used a Bayesian statistical framework and sequence data from all known bat-CoVs (including 630 novel CoV sequences) to study their macroevolution, cross-species transmission, and dispersal in China. We find that host-switching was more frequent and across more distantly related host taxa in alpha-than beta-CoVs, and more highly constrained by phylogenetic distance for beta-CoVs. We show that inter-family and -genus switching is most common in Rhinolophidae and the genus Rhinolophus . Our analyses identify the host taxa and geographic regions that define hotspots of CoV evolutionary diversity in China that could help target bat-CoV discovery for proactive zoonotic disease surveillance. Finally, we present a phylogenetic analysis suggesting a likely origin for SARS-CoV-2 in Rhinolophus spp. bats.

From mainland to islands: colonization history in the tree frog Kurixalus (Anura: Rhacophoridae)

The origin and colonization history of Kurixalus, a genus of small arboreal tree frogs breeding exclusively in shallow swamps, is under disputed. On the basis of comprehensive sampling program, the evolutionary history of Kurixalus is investigated based on 3 mitochondrial genes. Our results indicate that the genus Kurixalus originated in the Asian mainland and subsequently arrived at its current distribution in Borneo, Taiwan, Ryukyu, and Hainan islands by a series of dispersal events. Moreover, the colonization of Taiwan from mainland Asia has occurred 2 times. The initial colonization of Taiwan occurred at 3.46–8.68 Mya (95% highest posterior density), which rejects the hypothesis that Kurixalus probably originated from Taiwan during the early Oligocene and favors the model of Neogene-origin rather than the model of Quaternary-origin for Taiwanese Kurixalus. Kurixalus eiffingeri has dispersed from Taiwan to the Ryukyus once or 2 times pending more data. Both transoceanic dispersal and landbridge dispersal have played a role in the colonization process; the former resulted in the colonization of Taiwan and the Ryukyus and the latter led to the colonization of Borneo and Hainan.

Integrated approaches to identifying cryptic bat species in areas of high endemism: The case of Rhinolophus andamanensis in the Andaman Islands

The diversity of bats worldwide includes large numbers of cryptic species, partly because divergence in acoustic traits such as echolocation calls are under stronger selection than differences in visual appearance in these nocturnal mammals. Island faunas often contain disproportionate numbers of endemic species, and hence we might expect cryptic, endemic species to be discovered relatively frequently in bats inhabiting islands. Species are best defined when multiple lines of evidence supports their diagnosis. Here we use morphometric, acoustic, and molecular phylogenetic data to show that a horseshoe bat in the Andaman Islands is distinct in all three aspects, supporting its status as a distinct species. We recommend investigation into possible new and endemic bat species on islands by using integrated approaches that provide independent lines of evidence for taxonomic distinctiveness. We provide a formal redescription of the taxon newly raised to species level, Rhinolophus andamanensis Dobson, 1872.

Vicariance and Its Impact on the Molecular Ecology of a Chinese Ranid Frog Species-Complex (Odorrana schmackeri, Ranidae)

Paleogeological events and Pleistocene climatic fluctuations have had profound influences on the genetic patterns and phylogeographic structure of species in southern China. In this study, we investigated the population genetic structure and Phylogeography of the Odorrana schmackeri species complex, mountain stream-dwelling odorous frogs, endemic to southern China. We obtained mitochondrial sequences (1,151bp) of the complete ND2 gene and two flanking tRNAs of 511 individuals from 25 sites for phylogeographic analyses. Phylogenetic reconstruction revealed seven divergent evolutionary lineages, with mean pairwise (K2P) sequence distances from 7.8% to 21.1%, except for a closer ND2 distance (3.4%). The complex geological history of southern China drove matrilineal divergence in the O. schmackeri species complex into highly structured geographical units. The first divergence between lineage A+B and other lineages (C-G) had likely been influenced by the uplift of coastal mountains of Southeast China during the Mio-Pliocene period. The subsequent divergences between the lineages C-G may have followed the formation of the Three Gorges and the intensification of the East Asian summer monsoon during the late Pliocene and early Pleistocene. Demographic analyses indicated that major lineages A and C have been experienced recent population expansion (c. 0.045–0.245 Ma) from multiple refugia prior to the Last Glacial Maximum (LGM). Molecular analysis suggest that these seven lineages may represent seven different species, three described species and four cryptic species and should at least be separated into seven management units corresponding to these seven geographic lineages for conservation.

Molecular phylogeography and ecological niche modelling of a widespread herbaceous climber, Tetrastigma hemsleyanum (Vitaceae): insights into Plio-Pleistocene range dynamics of evergreen forest in subtropical China

Warm-temperate evergreen (WTE) forest represents the typical vegetation type of subtropical China, but how its component species responded to past environmental change remains largely unknown. Here, we reconstruct the evolutionary history of Tetrastigma hemsleyanum, an herbaceous climber restricted to the WTE forest. Twenty populations were genotyped using chloroplast DNA sequences and nuclear microsatellite loci to assess population structure and diversity, supplemented by phylogenetic dating, ancestral area reconstructions and ecological niche modeling (ENM) of the species distributions during the Last Glacial Maximum (LGM) and at present. Lineages in Southwest vs Central-South-East China diverged through climate/tectonic-induced vicariance of an ancestral southern range during the early Pliocene. Long-term stability in the Southwest contrasts with latitudinal range shifts in the Central-South-East region during the early-to-mid-Pleistocene. Genetic and ENM data strongly suggest refugial persistence in situ at the LGM. Pre-Quaternary environmental changes appear to have had a persistent influence on the population genetic structure of this subtropical WTE forest species. Our findings suggest relative demographic stability of this biome in China over the last glacial-interglacial cycle, in contrast with palaeobiome reconstructions showing that this forest biome retreated to areas of today's tropical South China during the LGM.

Taxonomic implications of geographical variation in Rhinolophus affinis (Chiroptera: Rhinolophidae) in mainland Southeast Asia

BACKGROUND

Rhinolophusaffinis sensu lato isdistributed throughout Southeast Asia. The taxonomic status of forms attributed to the species is unclear due to the limited sample size with incomplete datasets and the taxa have high variation in morphology and echolocation call frequency. The aim of the study was to evaluate the distribution and taxonomic status of the subspecific forms of R. affinis in mainland Southeast Asia using large sample size with multiple datasets, including morphological, acoustic, and genetic data, both to elucidate taxonomic relationships and to test for congruence between these datasets.

RESULTS

Three morphological forms were confirmed within the region; two concur with previously recognized taxa, namely R.affinis macrurus andR.affinis superans,and are strongly supported by morphological and genetic data. The third form is morphologically distinct, but its taxonomic status remains unclear. It is probable that this third form represents a distinct taxonomic entity; however, more data are required to confirm this. R. a. macrurus is known from the north of peninsular Thailand, Cambodia, Myanmar, Laos, and Vietnam (Indochinese subregion); R. a. superans is found throughout the Thai-Malay Peninsula (Sundaic subregion); whilst the third form is presently known from east central Myanmar (Shan state) and lower northern Vietnam (Nghe An Province).

CONCLUSIONS

Our results suggest that at least three morphological forms occur in mainland Southeast Asia including one form which appears to be new to science. Echolocation call data for R. affinis are not a robust taxonomic tool as it shows a significant degree of variation which is not explained or supported by genetic and morphological findings. This study highlights significant levels of morphological variation in mainland Southeast Asia and provides an essential basis for further studies aiming to understand the population genetics, phylogeography, and taxonomy of the species.

Differential introgression among loci across a hybrid zone of the intermediate horseshoe bat (Rhinolophus affinis)

Background

Hybrid zones formed by the secondary contact of divergent lineages represent natural laboratories for studying the genetic basis of speciation. Here we tested for patterns of differential introgression among three X-linked and 11 autosomal regions to identify candidate loci related to either reproductive isolation or adaptive introgression across a hybrid zone between two Chinese mainland subspecies of the intermediate horseshoe bat Rhinolophus affinis: R. a. himalayanus and R. a. macrurus.

Results

Our results support the previous suggestion that macrurus formed when a third subspecies (R. a. hainanus) recolonized the mainland from Hainan Island, and that himalayanus is the ancestral taxon. However, this overall evolutionary history was not reflected in all loci examined, with considerable locus-wise heterogeneity seen in gene tree topologies, levels of polymorphism, genetic differentiation and rates of introgression. Coalescent simulations suggested levels of lineage mixing seen at some nuclear loci might result from incomplete lineage sorting. Isolation with migration models supported evidence of gene flow across the hybrid zone at one intronic marker of the hearing gene Prestin.

Conclusions

We suggest that phylogenetic discordance with respect to the species tree seen here is likely to arise via a combination of incomplete lineage sorting and a low incidence of introgression although we cannot rule out other explanations such as selection and recombination. Two X-linked loci and one autosomal locus were identified as candidate regions related to reproductive isolation across the hybrid zone. Our work highlights the importance of including multiple genomic regions in characterizing patterns of divergence and gene flow across a hybrid zone.

Sensory trait variation in an echolocating bat suggests roles for both selection and plasticity

Background

Across heterogeneous environments selection and gene flow interact to influence the rate and extent of adaptive trait evolution. This complex relationship is further influenced by the rarely considered role of phenotypic plasticity in the evolution of adaptive population variation. Plasticity can be adaptive if it promotes colonization and survival in novel environments and in doing so may increase the potential for future population differentiation via selection. Gene flow between selectively divergent environments may favour the evolution of phenotypic plasticity or conversely, plasticity itself may promote gene flow, leading to a pattern of trait differentiation in the presence of gene flow. Variation in sensory traits is particularly informative in testing the role of environment in trait and population differentiation. Here we test the hypothesis of ‘adaptive differentiation with minimal gene flow’ in resting echolocation frequencies (RF) of Cape horseshoe bats (Rhinolophus capensis) across a gradient of increasingly cluttered habitats.

Results

Our analysis reveals a geographically structured pattern of increasing RF from open to highly cluttered habitats in R. capensis; however genetic drift appears to be a minor player in the processes influencing this pattern. Although Bayesian analysis of population structure uncovered a number of spatially defined mitochondrial groups and coalescent methods revealed regional-scale gene flow, phylogenetic analysis of mitochondrial sequences did not correlate with RF differentiation. Instead, habitat discontinuities between biomes, and not genetic and geographic distances, best explained echolocation variation in this species. We argue that both selection for increased detection distance in relatively less cluttered habitats and adaptive phenotypic plasticity may have influenced the evolution of matched echolocation frequencies and habitats across different populations.

Conclusions

Our study reveals significant sensory trait differentiation in the presence of historical gene flow and suggests roles for both selection and plasticity in the evolution of echolocation variation in R. capensis. These results highlight the importance of population level analyses to i) illuminate the subtle interplay between selection, plasticity and gene flow in the evolution of adaptive traits and ii) demonstrate that evolutionary processes may act simultaneously and that their relative influence may vary across different environments.

Phylogeography of the Rickett's big-footed bat, Myotis pilosus (Chiroptera: Vespertilionidae): a novel pattern of genetic structure of bats in China

Background

China is characterized by complex topographic structure and dramatic palaeoclimatic changes, making species biogeography studies particularly interesting. Previous researchers have also demonstrated multiple species experienced complex population histories, meanwhile multiple shelters existed in Chinese mainland. Despite this, species phylogeography is still largely unexplored. In the present study, we used a combination of microsatellites and mitochondrial DNA (mtDNA) to investigate the phylogeography of the east Asian fish-eating bat (Myotis pilosus).

Results

Phylogenetic analyses showed that M. pilosus comprised three main lineages: A, B and C, which corresponded to distinct geographic populations of the Yangtze Plain (YTP), Sichuan Basin (SCB) and North and South of China (NSC), respectively. The most recent common ancestor of M. pilosus was dated as 0.25 million years before present (BP). Population expansion events were inferred for populations of Clade C, North China Plain region, Clade B and YunGui Plateau region at 38,700, 15,900, 4,520 and 4,520 years BP, respectively. Conflicting results were obtained from mtDNA and microsatellite analyses; strong population genetic structure was obtained from mtDNA data but not microsatellite data. The microsatellite data indicated that genetic subdivision fits an isolation-by-distance (IBD) model, but the mtDNA data failed to support this model.

Conclusions

Our results suggested that Pleistocene climatic oscillations might have had a profound influence on the demographic history of M. pilosus. Spatial genetic structures of maternal lineages that are different from those observed in other sympatric bats species may be as a result of interactions among special population history and local environmental factors. There are at least three possible refugia for M. pilosus during glacial episodes. Apparently contradictory genetic structure patterns of mtDNA and microsatellite could be explained by male-mediated gene flow among populations. This study also provides insights on the necessity of conservation of M. pilosus populations to conserve this genetic biodiversity, especially in the areas of YTP, SCB and NSC regions.

Diversification and Demography of the Oriental Garden Lizard (Calotes versicolor) on Hainan Island and the Adjacent Mainland

The Oriental garden lizard (Calotes versicolor) is one of the few non-gekkonid lizards that are geographically widespread in the tropics. We investigated its population dynamics on Hainan Island and the adjacent mainland of China and Vietnam, focusing on the impact of cyclic upheaval and submergence of land bridges during the Pleistocene. Our Bayesian phylogenetic analysis reveals two mitochondrial lineages, A and B, which are estimated to have coalesced about 0.26 million years ago (95% credibility interval: 0.05–0.61 million years ago). Lineage A contains individuals mainly from central and southern Wuzhi Mountain on Hainan Island, whereas lineage B mainly comprises individuals from other sites on the island plus the adjacent mainland. The estimated coalescence times within lineages A (0.05 million years ago) and B (0.13 million years ago) fall within a period of cyclical land-bridge formation and disappearance in the Pleistocene. A spatial analysis of molecular variance identified two distinct population groupings: I, primarily containing lineage A, and II, mainly consisting of lineage B. However, haplotypes from lineages A and B occur sympatrically, suggesting that gene flow is ongoing. Neither Wuzhi Mountain nor Qiongzhou Strait and Gulf of Tonkin act as barriers to gene flow among C. versicolor populations. Analyses of the data using mismatch distributions and extended Bayesian skyline plots provide evidence of a relatively stable population size through time for Group I, and moderate population expansions and contractions during the end of the Pleistocene for Group II. We conclude that the phylogeographical patterns of C. versicolor are the combined product of Pleistocene sea-level oscillations and nonphysical barriers to gene flow.

Lineage divergence and historical gene flow in the Chinese horseshoe bat (Rhinolophus sinicus)

Closely related taxa living in sympatry provide good opportunities to investigate the origin of barriers to gene flow as well as the extent of reproductive isolation. The only two recognized subspecies of the Chinese rufous horseshoe bat Rhinolophus sinicus are characterized by unusual relative distributions in which R. s. septentrionalis is restricted to a small area within the much wider range of its sister taxon R. s. sinicus. To determine the history of lineage divergence and gene flow between these taxa, we applied phylogenetic, demographic and coalescent analyses to multi-locus datasets. MtDNA gene genealogies and microsatellite-based clustering together revealed three divergent lineages of sinicus, corresponding to Central China, East China and the offshore Hainan Island. However, the central lineage of sinicus showed a closer relationship with septentrionalis than with other lineages of R. s. sinicus, in contrary to morphological data. Paraphyly of sinicus could result from either past asymmetric mtDNA introgression between these two taxa, or could suggest septentrionalis evolved in situ from its more widespread sister subspecies. To test between these hypotheses, we applied coalescent-based phylogenetic reconstruction and Approximate Bayesian Computation (ABC). We found that septentrionalis is likely to be the ancestral taxon and therefore a recent origin of this subspecies can be ruled out. On the other hand, we found a clear signature of asymmetric mtDNA gene flow from septentrionalis into central populations of sinicus yet no nuclear gene flow, thus strongly pointing to historical mtDNA introgression. We suggest that the observed deeply divergent lineages within R. sinicus probably evolved in isolation in separate Pleistocene refugia, although their close phylogeographic correspondence with distinct eco-environmental zones suggests that divergent selection might also have promoted broad patterns of population genetic structure.

Effect of Pleistocene climatic oscillations on the phylogeography and demography of red knobby newt (Tylototriton shanjing) from southwestern China

Factors that determine the genetic structure of species in southwestern China remain largely unknown. In this study, phylogeography and demography of Tylototriton shanjing was investigated from a mitochondrial perspective to address the role of the Quaternary ice ages in shaping phylogeographic history and genetic diversity of Yunnan. A total of 146 individuals from 19 populations across the entire range of the species were collected. We detected four maternal phylogenetic lineages corresponding to four population groups, and found that major glaciation events during the Pleistocene have triggered the intra-specific divergence. Coalescent simulations indicated that the populations retreated to different refugia located in southern Yunnan, northwestern Yunnan, the border region of western Yunnan with Myanmar, and middle-western Yunnan, respectively, during previous glacial periods in the Pleistocene, and these four refugia were not retained during the Last Glacial Maximum. Population expansions occurred during the last inter-glaciation, during which ice core and pollen data indicated that the temperature and precipitation gradually increased, and declines of population sizes started after the beginning of the Last Glacial Maximum when the climate became cooler and dryer. The paleo-drainage system had no contribution to the current genetic structure and the rivers were not dispersal barriers for this salamander.

Genealogy and Demographic History of a Widespread Amphibian throughout Indochina

Relatively little is known about spatial patterns of cryptic diversity in tropical species and the processes that generate them. Few studies examine the geographic distribution of genetic lineages in Southeast Asia, an area hypothesized to harbor substantial cryptic diversity. We investigated the evolutionary history of Asian tree frogs of the Polypedates leucomystax complex (n = 172) based on 1800 bp of the mtDNA genes ND1 and cytochrome b and tested hypotheses pertaining to climate, geology, and dispersal patterns. Analyses revealed substantial genetic diversity and lineage divergence throughout the region with evidence for widespread sympatric lineages and a general north versus south clustering. Relaxed molecular clock analysis and tests for demographic expansion identified an initial cladogenesis during the Miocene with subsequent Plio-Pleistocene diversification, with the former corresponding to periods of increased aridity and the onset of monsoonal weather systems. Rates of diversification were relatively constant until the Early Pleistocene when rates increased exponentially. We found equivocal evidence for isolation-by-distance and a potential role of some landscape features as partial barriers to dispersal. Finally, our analyses showed that divergence between insular and mainland populations occurred before Homo sapiens colonized Southeast Asia, suggesting that historical human-mediated dispersal did not drive insular diversification. Our results suggested that demographic expansion in the Late Pleistocene resulted in widespread sympatric lineages in the P. leucomystax complex throughout southern China and Indochina and further clarified the evolutionary history of lineages within P. leucomystax.

Four new bat species (Rhinolophus hildebrandtii complex) reflect Plio-Pleistocene divergence of dwarfs and giants across an Afromontane archipelago

Gigantism and dwarfism evolve in vertebrates restricted to islands. We describe four new species in the Rhinolophus hildebrandtii species-complex of horseshoe bats, whose evolution has entailed adaptive shifts in body size. We postulate that vicissitudes of palaeoenvironments resulted in gigantism and dwarfism in habitat islands fragmented across eastern and southern Africa. Mitochondrial and nuclear DNA sequences recovered two clades of R. hildebrandtii senso lato which are paraphyletic with respect to a third lineage (R. eloquens). Lineages differ by 7.7 to 9.0% in cytochrome b sequences. Clade 1 includes R. hildebrandtii sensu stricto from the east African highlands and three additional vicariants that speciated across an Afromontane archipelago through the Plio-Pleistocene, extending from the Kenyan Highlands through the Eastern Arc, northern Mozambique and the Zambezi Escarpment to the eastern Great Escarpment of South Africa. Clade 2 comprises one species confined to lowland savanna habitats (Mozambique and Zimbabwe). A third clade comprises R. eloquens from East Africa. Speciation within Clade 1 is associated with fixed differences in echolocation call frequency, and cranial shape and size in populations isolated since the late Pliocene (ca 3.74 Mya). Relative to the intermediate-sized savanna population (Clade 2), these island-populations within Clade 1 are characterised by either gigantism (South African eastern Great Escarpment and Mts Mabu and Inago in Mozambique) or dwarfism (Lutope-Ngolangola Gorge, Zimbabwe and Soutpansberg Mountains, South Africa). Sympatry between divergent clades (Clade 1 and Clade 2) at Lutope-Ngolangola Gorge (NW Zimbabwe) is attributed to recent range expansions. We propose an "Allometric Speciation Hypothesis", which attributes the evolution of this species complex of bats to divergence in constant frequency (CF) sonar calls. The origin of species-specific peak frequencies (overall range = 32 to 46 kHz) represents the allometric effect of adaptive divergence in skull size, represented in the evolution of gigantism and dwarfism in habitat islands.

Phylogeography and demographic history of Babina pleuraden (Anura, Ranidae) in southwestern China

Factors that determine genetic structure of species in southwestern China remain largely unknown. In this study, sequences of two mitochondrial genes (COI and cyt b) were determined to investigate the phylogeography and demography of Babina pleuraden, a pond frog endemic to southwestern China. A total of 262 individuals from 22 populations across the entire range of the species were collected. Our results indicate that B. pleuraden comprises five well-supported mitochondrial lineages roughly corresponding to five geographical areas. The phylogeographic structure of B. pleuraden has been shaped primarily by the unique regional responses of the Yunnan Plateau to the rapid uplift of the Qinghai-Tibetan Plateau occurred c. 2.5 Mya (B phrase of Qingzang Movement) and climatic oscillation during middle Pleistocene (c. 0.64–0.36 Mya), rather than by the paleo-drainage systems. The present wide distribution of the species has resulted from recent population expansion (c. 0.053–0.025 Mya) from multiple refugia prior to the Last Glacial Maximum, corresponding to the scenario of “refugia within refugia”.

Correlated genetic and ecological diversification in a widespread southern African horseshoe bat

The analysis of molecular data within a historical biogeographical framework, coupled with ecological characteristics can provide insight into the processes driving diversification. Here we assess the genetic and ecological diversity within a widespread horseshoe bat Rhinolophus clivosus sensu lato with specific emphasis on the southern African representatives which, although not currently recognized, were previously described as a separate species R. geoffroyi comprising four subspecies. Sequence divergence estimates of the mtDNA control region show that the southern African representatives of R. clivosus s.l. are as distinct from samples further north in Africa than they are from R. ferrumequinum, the sister-species to R. clivosus. Within South Africa, five genetically supported geographic groups exist and these groups are corroborated by echolocation and wing morphology data. The groups loosely correspond to the distributions of the previously defined subspecies and Maxent modelling shows a strong correlation between the detected groups and ecoregions. Based on molecular clock calibrations, it is evident that climatic cycling and related vegetation changes during the Quaternary may have facilitated diversification both genetically and ecologically.

Species delimitation and phylogeography of Aphonopelma hentzi (Araneae, Mygalomorphae, Theraphosidae): cryptic diversity in North American tarantulas

BACKGROUND

The primary objective of this study is to reconstruct the phylogeny of the hentzi species group and sister species in the North American tarantula genus, Aphonopelma, using a set of mitochondrial DNA markers that include the animal "barcoding gene". An mtDNA genealogy is used to consider questions regarding species boundary delimitation and to evaluate timing of divergence to infer historical biogeographic events that played a role in shaping the present-day diversity and distribution. We aimed to identify potential refugial locations, directionality of range expansion, and test whether A. hentzi post-glacial expansion fit a predicted time frame.

METHODS AND FINDINGS

A Bayesian phylogenetic approach was used to analyze a 2051 base pair (bp) mtDNA data matrix comprising aligned fragments of the gene regions CO1 (1165 bp) and ND1-16S (886 bp). Multiple species delimitation techniques (DNA tree-based methods, a "barcode gap" using percent of pairwise sequence divergence (uncorrected p-distances), and the GMYC method) consistently recognized a number of divergent and genealogically exclusive groups.

CONCLUSIONS

The use of numerous species delimitation methods, in concert, provide an effective approach to dissecting species boundaries in this spider group; as well they seem to provide strong evidence for a number of nominal, previously undiscovered, and cryptic species. Our data also indicate that Pleistocene habitat fragmentation and subsequent range expansion events may have shaped contemporary phylogeographic patterns of Aphonopelma diversity in the southwestern United States, particularly for the A. hentzi species group. These findings indicate that future species delimitation approaches need to be analyzed in context of a number of factors, such as the sampling distribution, loci used, biogeographic history, breadth of morphological variation, ecological factors, and behavioral data, to make truly integrative decisions about what constitutes an evolutionary lineage recognized as a "species".