Nuclear gene phylogeography using PHASE: dealing with unresolved genotypes, lost alleles, and systematic bias in parameter estimation

The role of Pleistocene dispersal in shaping species richness of sky island wintergreens from the Himalaya-Hengduan Mountains

In addition to topography and climate, biogeographic dispersal has been considered to influence plant diversity in the Himalaya-Hengduan Mountains (HHM), yet, the mode and tempo of sky island dispersal and its influence on species richness has been little explored. Through phylogenetic analysis of Gaultheria ser. Trichophyllae, a sky island alpine clade within the HHM, we test the hypothesis that dispersal has affected current local species richness. We inferred the dynamics of biogeographic dispersal with correlation tests on direction, distance, occurrence time, and regional species richness. We found that G. ser. Trichophyllae originated at the end of the Miocene and mostly dispersed toward higher longitudes (eastward). In particular, shorter intra-regional eastward dispersals and longer inter-regional westward dispersals were most frequently observed. We detected a prevalence of eastward intra-region dispersals in both glacial periods and interglacials. These dispersals may have been facilitated by the reorganization of paleo-drainages and monsoon intensification through time. We suggest that the timing of dispersal corresponding to glacial periods and the prevalence of intra-region dispersal, rather than dispersal frequency, most influenced the pattern of species richness of G. ser. Trichophyllae. This study facilitates a more comprehensive understanding of biodiversity in the sky islands within the HHM.

Genetic signatures of lineage fusion closely resemble population decline

Accurate interpretation of the genetic signatures of past demographic events is crucial for reconstructing evolutionary history. Lineage fusion (complete merging, resulting in a single panmictic population) is a special case of secondary contact that is seldom considered. Here, the circumstances under which lineage fusion can be distinguished from population size constancy, growth, bottleneck, and decline were investigated. Multi-locus haplotype data were simulated under models of lineage fusion with different divergence versus sampling lag times (D:L ratios). These pseudo-observed datasets also differed in their allocation of a fixed amount of sequencing resources (number of sampled alleles, haplotype length, number of loci). Distinguishability of lineage fusion versus each of 10 untrue non-fusion scenarios was quantified based on six summary statistics (neutrality tests). Some datasets were also analyzed using extended Bayesian skyline plots. Results showed that signatures of lineage fusion very closely resemble those of decline-high distinguishability was generally limited to the most favorable scenario (D:L = 9), using the most sensitive summary statistics (F S and Z nS), coupled with the optimal sequencing resource allocation (maximizing number of loci). Also, extended Bayesian skyline plots often erroneously inferred population decline. Awareness of the potential for lineage fusion to carry the hallmarks of population decline is critical.

Temporal and spatial diversification along the Amazonia-Cerrado transition in Neotropical treefrogs of the Boana albopunctata species group

Despite extensive research on biodiversity in Neotropical forests, biodiversity in seasonally dry, open biomes in South America has been underestimated until recently. We leverage a widespread group, Boana albopunctata, to uncover cryptic lineages and investigate the timing of diversification in Neotropical anurans with a focus on dry diagonal biomes (Cerrado, Caatinga and Chaco) and the ecotone between Amazonia and the Cerrado. We inferred a multilocus phylogeny of the B. albopunctata species group that includes 15 of 18 described species, recovered two cryptic species, and reconstructed the timing of diversification among species distributed across multiple South American biomes. One new potential species (B. aff. steinbachi), sampled in the Amazonian state of Acre, clustered within the B. calcara-fasciata species complex and is close to B. steinbachi. A second putative new species (B. aff. multifasciata), sampled in the Amazonia-Cerrado ecotone, is closely related to B. multifasciata. Lastly, we place a recently identified Cerrado lineage (B. aff. albopuncata) into the B. albopunctata species group phylogeny for the first time. Our ancestral range reconstruction showed that species in the B. albopuctata group likely dispersed from Amazonia-Cerrado into the dry-diagonal and Atlantic Forest. Intraspecies demography showed, for both B. raniceps and B. albopunctata, signs of rapid expansion across the dry diagonal. Similarly, for one clade of B. multifasciata, our analyses support an invasion of the Cerrado from Amazonia, followed by a rapid expansion across the open diagonal biomes. Thus, our study recovers several recent divergences along the Amazonia-Cerrado ecotone in northern Brazil. Tectonic uplift and erosion in the late Miocene and climate oscillations in the Pleistocene corresponded with estimated divergence times in the dry diagonal and Amazonia-Cerrado ecotone. Our study highlights the importance of these threatened open formations in the generation of biodiversity in the Neotropics.

SNP-based Phylogenomic Inference in Holarctic Ground Squirrels (Urocitellus)

Resolution of rapid evolutionary radiations requires harvesting maximal signal from phylogenomic datasets. However, studies of non-model clades often target conserved loci that are characterized by reduced information content, which can negatively affect gene tree precision and species tree accuracy. Single nucleotide polymorphism (SNP)-based methods are an underutilized but potentially valuable tool for estimating phylogeny and divergence times because they do not rely on resolved gene trees, allowing information from many or all variant loci to be leveraged in species tree reconstruction. We evaluated the utility of SNP-based methods in resolving phylogeny of Holarctic ground squirrels (Urocitellus), a radiation that has been difficult to disentangle, even in prior phylogenomic studies. We inferred phylogeny from a dataset of >3,000 ultraconserved element loci (UCEs) using two methods (SNAPP, SVDquartets) and compared our results with a new mitogenome phylogeny. We also systematically evaluated how phasing of UCEs improves per-locus information content, and inference of topology and other parameters within each of these SNP-based methods. Phasing improved topological resolution and branch length estimation at shallow levels (within species complexes), but less so at deeper levels, likely reflecting true uncertainty due to ancestral polymorphisms segregating in these rapidly diverging lineages. We resolved several key clades in Urocitellus and present targeted opportunities for future phylogenomic inquiry. Our results extend the roadmap for use of SNPs to address vertebrate radiations and support comparative analyses at multiple temporal scales.

Genome-wide assessment elucidates connectivity and the evolutionary history of the highly dispersive marine invertebrate Littoraria flava (Littorinidae: Gastropoda)

An important goal of marine population genetics is to understand how spatial connectivity patterns are influenced by historical and evolutionary factors. In this study, we evaluate the demographic history and population structure of Littoraria flava, a highly dispersive marine gastropod in the Brazilian intertidal zone. To test the hypotheses that the species has (1) historically high levels of gene flow on a macrogeographical spatial scale and (2) a distribution in rocky shores that consists of subpopulations, we collected specimens along the Brazilian coastline and combined different sets of genetic markers (mitochondrial DNA, ITS-2 and single nucleotide polymorphisms) with niche-based modelling to predict its palaeodistribution. Low genetic structure was observed, as well as high gene flow over long distances. The demographic analyses suggest that L. flava has had periods of population bottlenecks followed by expansion. According to both palaeodistribution and coalescent simulations, these expansion events occurred during the Pleistocene interglacial cycles (21 kya) and the associated climatic changes were the probable drivers of the distribution of the species. This is the first phylogeographical study of a marine gastropod on the South American coast based on genomic markers associated with niche modelling.

Evolutionary history of Sundaland shrews (Eulipotyphla: Soricidae: Crocidura) with a focus on Borneo

The hyperdiverse shrew genus Crocidura is one of few small mammal genera distributed across Sundaland and all of its boundaries. This represents a rare opportunity to study the geological history of this region through the evolutionary history of these shrews. We generate a phylogeny of all recognized species of Sundaland Crocidura and show that most speciation events took place during the Pleistocene, prior to the inundation of the Sunda Shelf around 400 000 years ago. We find east–west differentiation within two separate lineages on Borneo, and that the current taxonomy of its two endemic species does not reflect evolutionary history, but ecophenotypic variation of plastic traits related to elevation. Sulawesi shrews are monophyletic, with a single notable exception: the black-footed shrew (C. nigripes). We show that the black-footed shrew diverged from its relatives on Borneo recently, suggesting a human-assisted breach of Wallace’s line. Overall, the number of Crocidura species, especially on Borneo, probably remains an underestimate.

Genetic, morphological and acoustic differentiation of African trident bats (Rhinonycteridae: Triaenops)

Rhinonycteridae (trident bats) are a small Palaeotropical family of insectivorous bats allied to Hipposideridae. Their taxonomy has been in a state of flux. Here, we use mitochondrial and nuclear sequences to evaluate species relationships, confirming the monophyly of both Triaenops and Paratriaenops. Although most Triaenops afer specimens are recovered as a group, mitochondrial analyses strongly support some Kenyan individuals as members of Triaenops persicus. Analyses of four nuclear introns (ACOX2, COPS7A, RODGI and STAT5A) strongly support the mitochondrial topology. Morphometric analysis of the skull, external morphology and echolocation calls confirm that the Triaenops from the Rift Valley in Kenya (Nakuru, Baringo and Pokot counties) are distinct from typical T. afer in coastal (Kilifi and Kwale counties) or interior (Laikipia and Makueni counties) colonies. We interpret these analyses to indicate that two species of Triaenops occur in East Africa: T. afer in coastal regions along the Indian Ocean and in the highlands of central Kenya and Ethiopia, and T. persicus in the Rift Valley of Kenya. Although they appear widely disjunct from Middle Eastern populations, Kenyan T. persicus might be more widely distributed in the Rift Valley; they are somewhat differentiated from Middle Eastern populations in terms of both cranial morphology and vocalizations.

Habitat fragmentation differentially shapes neutral and immune gene variation in a tropical bird species

Habitat fragmentation is a major cause of biodiversity loss, responsible for an alteration of intraspecific patterns of neutral genetic diversity and structure. Although neutral genetic variation can be informative for demographic inferences, it may be a poor predictor of adaptive genetic diversity and thus of the consequences of habitat fragmentation on selective evolutionary processes. In this context, we contrasted patterns of genetic diversity and structure of neutral loci (microsatellites) and immune genes (i.e., toll-like receptors) in an understorey bird species, the wedge-billed woodcreeper Glyphorynchus spirurus. The objectives were (1) to investigate forest fragmentation effects on population genetic diversity, (2) to disentangle the relative role of demography (genetic drift and migration) and selection, and (3) to assess whether immunogenetic patterns could be associated with variation of ectoparasite (i.e., ticks) pressures. Our results revealed an erosion of neutral genetic diversity and a substantial genetic differentiation among fragmented populations, resulting from a decrease in landscape connectivity and leading to the divergence of distinct genetic pools at a small spatial scale. Patterns of genetic diversity observed for TLR4 and TLR5 were concordant with neutral genetic patterns, whereas those observed for TLR3 and TLR21 were discordant. This result underlines that the dominant evolutionary force shaping immunogenetic diversity (genetic drift vs. selection) may be different depending on loci considered. Finally, tick prevalence was higher in fragmented environments. We discussed the hypothesis that pathogen selective pressures may contribute to maintain adaptive genetic diversity despite the negative demographic effect of habitat fragmentation on neutral genetic diversity.

Molecular systematics and phylogeography of a widespread Neotropical avian lineage: evidence for cryptic speciation with protracted gene flow throughout the Late Quaternary

Here we use an integrative approach, including coalescent-based methods, isolation–migration and species distribution models, to infer population structure, divergence times and diversification in the two species of the genus Cymbilaimus (Aves, Thamnophilidae). Our results support a recent and rapid diversification with both incomplete lineage sorting and gene flow shaping the evolutionary history of Cymbilaimus. The spatio-temporal pattern of cladogenesis suggests that Cymbilaimus originated in the north/western portion of cis-Andean South America and then diversified into the Brazilian Shield and Central America after consolidation of the modern Amazonian drainage and the Andean range. This evolutionary scenario is explained by cycles of range expansion and dispersal, followed by isolation, and recurrent gene flow, during the last 1.2 Myr. Our results agree with those recently reported for other closely related suboscine lineages, whereby the window of introgression between closely related taxa remains open for up to a few million years after their original split. In Cymbilaimus, introgression was recurrent between C. lineatus and C. sanctaemariae, even after they acquired vocal and ecological differentiation, supporting the claim that at least in Neotropical suboscines, full reproductive compatibility may take millions of years to evolve and cannot be interpreted as synonymous with a lack of speciation.

Impact of the Mid-Pleistocene Revolution and Anthropogenic Factors on the Dispersion of Asian Black-Spined Toads (Duttaphrynus melanostictus)

Divergence-time estimation critically improves the understanding of biogeography processes underlying the distribution of species, especially when fossil data is not available. We hypothesise that the Asian black-spined toad, Duttaphrynus melanostictus, expanded into the Eastern Indomalaya following the Quaternary glaciations with the subsequent colonisation of new landscapes during the Last Glacial Maximum. Divergence dating inferred from 364 sequences of mitochondrial tRNAGly ND3 supported the emergence of a common ancestor to the three D. melanostictus clades around 1.85 (±0.77) Ma, matching with the Lower to Mid-Pleistocene transition. Duttaphrynus melanostictus then dispersed into Southeast Asia from the central Indo-Pacific and became isolated in the Southern Sundaic and Wallacea regions 1.43 (±0.10) Ma through vicariance as a result of sea level oscillations. The clade on the Southeast Asian mainland then colonised the peninsula from Myanmar to Vietnam and expanded towards Southeastern China at the end of the Mid-Pleistocene Revolution 0.84 (±0.32) Ma. Population dynamics further highlight an expansion of the Southeast Asian mainland population towards Taiwan, the Northeastern edge of the species' range after the last interglacial, and during the emergence of the Holocene human settlements around 7000 BP. Thus, the current divergence of D. melanostictus into three segregated clades was mostly shaped by Quaternary glaciations, followed by natural dispersion events over land bridges and accelerated by anthropogenic activities.

Molecular phylogenetics of slit-faced bats (Chiroptera: Nycteridae) reveal deeply divergent African lineages

The bat family Nycteridae contains only the genus Nycteris, which comprises 13 currently recognized species from Africa and the Arabian Peninsula, one species from Madagascar, and two species restricted to Malaysia and Indonesia in South-East Asia. We investigated genetic variation, clade membership, and phylogenetic relationships in Nycteridae with broad sampling across Africa for most clades. We sequenced mitochondrial cytochrome b (cytb) and four independent nuclear introns (2,166 bp) from 253 individuals. Although our samples did not include all recognized species, we recovered at least 16 deeply divergent monophyletic lineages using independent mitochondrial and multilocus nuclear datasets in both gene tree and species tree analyses. Mean pairwise uncorrected genetic distances among species-ranked Nycteris clades (17% for cytb and 4% for concatenated introns) suggest high levels of phylogenetic diversity in Nycteridae. We found a large number of designated clades whose members are distributed wholly or partly in East Africa (10 of 16 clades), indicating that Nycteris diversity has been historically underestimated and raising the possibility that additional unsampled and/or undescribed Nycteris species occur in more poorly sampled Central and West Africa. Well-resolved mitochondrial, concatenated nuclear, and species trees strongly supported African ancestry for SE Asian species. Species tree analyses strongly support two deeply diverged subclades that have not previously been recognized, and these clades may warrant recognition as subgenera. Our analyses also strongly support four traditionally recognized species groups of Nycteris. Mitonuclear discordance regarding geographic population structure in Nycteris thebaica appears to result from male-biased dispersal in this species. Our analyses, almost wholly based on museum voucher specimens, serve to identify species-rank clades that can be tested with independent datasets, such as morphology, vocalizations, distributions, and ectoparasites. Our analyses highlight the need for a comprehensive revision of Nycteridae.

Genetic variation and relationships among Afrotropical species of Myotis (Chiroptera: Vespertilionidae)

The genus Myotis is nearly cosmopolitan and the second-most speciose genus of mammals, but its Afrotropical members are few and poorly known. We analyzed phylogenetic and phylogeographic relationships of six of the eight known Afrotropical species using Cytb and sequences from four nuclear introns. Using Bayesian and maximum-likelihood approaches to generate single-locus, concatenated, and species trees, we confirmed prior evidence that the clade containing Afrotropical Myotis also contains both Palearctic and Indomalayan members. Additionally, we demonstrate that M. bocagii is sister to the Indian Ocean species M. anjouanensis, that this group is sister to M. tricolor and the Palearctic M. emarginatus, and find evidence suggesting that M. welwitschii is the earliest-diverging Afrotropical species and sister to the remainder. Although M. tricolor and M. welwitschii are both currently regarded as monotypic, both mitochondrial and nuclear data sets document significant, largely concordant geographic structure in each. Evidence for the distinction of two lineages within M. tricolor is particularly strong. On the other hand, geographic structure is lacking in M. bocagii, despite the current recognition of two subspecies in that species. Additional geographic sampling (especially at or near type localities), finer-scale sampling (especially in zones of sympatry), and integrative taxonomic assessments will be needed to better document this radiation and refine its nomenclature.

Oceanographic barriers, divergence, and admixture: Phylogeography and taxonomy of two putative subspecies of short-finned pilot whale

Genomic phylogeography plays an important role in describing evolutionary processes and their geographic, ecological, or cultural drivers. These drivers are often poorly understood in marine environments, which have fewer obvious barriers to mixing than terrestrial environments. Taxonomic uncertainty of some taxa (e.g., cetaceans), due to the difficulty in obtaining morphological data, can hamper our understanding of these processes. One such taxon, the short-finned pilot whale, is recognized as a single global species but includes at least two distinct morphological forms described from stranding and drive hunting in Japan, the "Naisa" and "Shiho" forms. Using samples (n = 735) collected throughout their global range, we examine phylogeographic patterns of divergence by comparing mitogenomes and nuclear SNP loci. Our results suggest three types within the species: an Atlantic Ocean type, a western/central Pacific and Indian Ocean (Naisa) type, and an eastern Pacific Ocean and northern Japan (Shiho) type. mtDNA control region differentiation indicates these three types form two subspecies, separated by the East Pacific Barrier: Shiho short-finned pilot whale, in the eastern Pacific Ocean and northern Japan, and Naisa short-finned pilot whale, throughout the remainder of the species' distribution. Our data further indicate two diverging populations within the Naisa subspecies, in the Atlantic Ocean and western/central Pacific and Indian Oceans, separated by the Benguela Barrier off South Africa. This study reveals a process of divergence and speciation within a globally-distributed, mobile marine predator, and indicates the importance of the East Pacific Barrier to this evolutionary process.

Amplicon sequencing of 42 nuclear loci supports directional gene flow between South Pacific populations of a hydrothermal vent limpet

In the past few decades, population genetics and phylogeographic studies have improved our knowledge of connectivity and population demography in marine environments. Studies of deep-sea hydrothermal vent populations have identified barriers to gene flow, hybrid zones, and demographic events, such as historical population expansions and contractions. These deep-sea studies, however, used few loci, which limit the amount of information they provided for coalescent analysis and thus our ability to confidently test complex population dynamics scenarios. In this study, we investigated population structure, demographic history, and gene flow directionality among four Western Pacific hydrothermal vent populations of the vent limpet Lepetodrilus aff. schrolli. These vent sites are located in the Manus and Lau back-arc basins, currently of great interest for deep-sea mineral extraction. A total of 42 loci were sequenced from each individual using high-throughput amplicon sequencing. Amplicon sequences were analyzed using both genetic variant clustering methods and evolutionary coalescent approaches. Like most previously investigated vent species in the South Pacific, L. aff. schrolli showed no genetic structure within basins but significant differentiation between basins. We inferred significant directional gene flow from Manus Basin to Lau Basin, with low to no gene flow in the opposite direction. This study is one of the very few marine population studies using >10 loci for coalescent analysis and serves as a guide for future marine population studies.

Alternated selection mechanisms maintain adaptive diversity in different demographic scenarios of a large carnivore

Background

Different population trajectories are expected to impact the signature of neutral and adaptive processes at multiple levels, challenging the assessment of the relative roles of different microevolutionary forces. Here, we integrate adaptive and neutral variability patterns to disentangle how adaptive diversity is driven under different demographic scenarios within the Iberian wolf (Canis lupus) range. We studied the persistent, the expanding and a small, isolated group within the Iberian wolf population, using 3 MHC class II genes (DRB1, DQA1, and DQB1), which diversity was compared with 39 microsatellite loci.

Results

Both the persistent and the expanding groups show evidence of balancing selection, revealed by a significant departure from neutrality at MHC loci, significant higher observed and expected heterozygosity and lower differentiation at MHC than at neutral loci, and signs of positive selection. However, despite exhibiting a significantly higher genetic diversity than the isolated group, the persistent group did not show significant excess of MHC heterozygotes. The expanding group, while showing a similar level of genetic diversity than the persistent group, displays by contrast a significant excess of MHC heterozygotes, which is compatible with the heterozygote advantage mechanism. Results are not clear regarding the role of drift and selection in the isolated group due to the small size of this population. Although diversity indices of MHC loci correspond to neutral expectations in the isolated group, accelerated MHC divergence, revealed by a higher differentiation at MHC than neutral loci, may indicate diversifying selection.

Conclusion

Different selective pressures were observed in the three different demographic scenarios, which are possibly driven by different selection mechanisms to maintain adaptive diversity.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1420-5) contains supplementary material, which is available to authorized users.

Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.

Diversity and relationships of Ampedini Gistel, 1848 (Coleoptera : Elateridae) in Switzerland and Europe

The click-beetle family Elateridae is an abundant and ecologically important component of the saproxylic fauna, and many species are used as key indicators of old-growth forests, a rare and endangered habitat in Europe. Among Elateridae, the tribe Ampedini includes several species with unclear taxonomic status. Here we examined the species status and the relationships between 26 of the 29 Ampedini species recorded in Switzerland using both mitochondrial (CO1) and nuclear (ITS2) markers, and compared our results with the CO1 data publicly available for Europe. Our tree-based species-delimitation analysis was largely congruent with traditional, morphology-based species, with the ITS2 and the combined (CO1+ITS2) datasets corresponding more precisely to morphology than the CO1 dataset. As expected, species with unclear or debated taxonomic status generally corresponded to closely related DNA-based species pairs. Our results support the taxonomic status of some of these ambiguous species, while others require further investigations, including a more exhaustive sampling and new morphological examinations. The Ampedini species tree provided here, which is the first attempt of a DNA-based Ampedini phylogeny, did not support the genera Brachygonus and Ampedus as monophyletic, but further investigations are necessary to confirm this result.

Allele Phasing Greatly Improves the Phylogenetic Utility of Ultraconserved Elements

Advances in high-throughput sequencing techniques now allow relatively easy and affordable sequencing of large portions of the genome, even for nonmodel organisms. Many phylogenetic studies reduce costs by focusing their sequencing efforts on a selected set of targeted loci, commonly enriched using sequence capture. The advantage of this approach is that it recovers a consistent set of loci, each with high sequencing depth, which leads to more confidence in the assembly of target sequences. High sequencing depth can also be used to identify phylogenetically informative allelic variation within sequenced individuals, but allele sequences are infrequently assembled in phylogenetic studies. Instead, many scientists perform their phylogenetic analyses using contig sequences which result from the de novo assembly of sequencing reads into contigs containing only canonical nucleobases, and this may reduce both statistical power and phylogenetic accuracy. Here, we develop an easy-to-use pipeline to recover allele sequences from sequence capture data, and we use simulated and empirical data to demonstrate the utility of integrating these allele sequences to analyses performed under the multispecies coalescent model. Our empirical analyses of ultraconserved element locus data collected from the South American hummingbird genus Topaza demonstrate that phased allele sequences carry sufficient phylogenetic information to infer the genetic structure, lineage divergence, and biogeographic history of a genus that diversified during the last 3 myr. The phylogenetic results support the recognition of two species and suggest a high rate of gene flow across large distances of rainforest habitats but rare admixture across the Amazon River. Our simulations provide evidence that analyzing allele sequences leads to more accurate estimates of tree topology and divergence times than the more common approach of using contig sequences.

Contrasting phylogeographic pattern among Eudyptes penguins around the Southern Ocean

Since at least the middle-Miocene, the Antarctic Polar Front (APF) and the Subtropical Front (STF) appear to have been the main drivers of diversification of marine biota in the Southern Ocean. However, highly migratory marine birds and mammals challenge this paradigm and the importance of oceanographic barriers. Eudyptes penguins range from the Antarctic Peninsula to subantarctic islands and some of the southernmost subtropical islands. Because of recent diversification, the number of species remains uncertain. Here we analyze two mtDNA (HVRI, COI) and two nuclear (ODC, AK1) markers from 13 locations of five putative Eudyptes species: rockhopper (E. filholi, E. chrysocome, and E. moseleyi), macaroni (E. chrysolophus) and royal penguins (E. schlegeli). Our results show a strong phylogeographic structure among rockhopper penguins from South America, subantarctic and subtropical islands supporting the recognition of three separated species of rockhopper penguins. Although genetic divergence was neither observed among macaroni penguins from the Antarctic Peninsula and sub-Antarctic islands nor between macaroni and royal penguins, population genetic analyses revealed population genetic structure in both cases. We suggest that the APF and STF can act as barriers for these species. While the geographic distance between colonies might play a role, their impact/incidence on gene flow may vary between species and colonies.

First steps towards assessing the evolutionary history and phylogeography of a widely distributed Neotropical grassland bird (Motacillidae: Anthus correndera)

Grasslands in southern South America are extensive ecosystems which harbor a unique biodiversity; however, studies on the evolution of their taxa are scarce. Here we studied the phylogeography and population history of the Correndera Pipit (Anthus correndera), a grassland specialist bird with a large breeding distribution in southern South America, with the goals of investigating its phylogeographic history and relate it to the historical development of South American grasslands. The mitochondrial NADH dehydrogenase subunit II gene (ND2) was sequenced in 66 individuals from 19 localities and the intron 9 of the sex-linked gene for aconitase (ACOI9) was sequenced from a subset of those individuals, including all five subspecies of A. correndera, as well as the closely related A. antarcticus. Phylogenetic analysis revealed two distinct lineages within the complex: the first (A) corresponding to Andean subspecies A. c. calcaratus and A. c. catamarcae and the second (B) including birds traditionally assigned to A. c. correndera, A. c. chilensis, A. c. grayi and some individuals of A. c. catamarcae. A. antarcticus is nested within this second lineage. These results were also supported by evidence of niche divergence for variables associated with precipitation. The oldest split between clade A and B was estimated at c. 0.37 Mya, during the middle Pleistocene. Species distribution models for the present and the Last Glacial Maximum (LGM) suggest that grassland areas in southern South America remained relatively stable, in contrast to the general view of a reduction in grassland cover in South America since the LGM. Recent divergences and low phylogeographic structure (for lowland vs. highland geographic groups, intra-population genetic variance was greater than inter-groups; e.g., for ACOI9: 95.47% and ND2: 51.51% respectively), suggest widespread gene flow between lowland populations.

Asymmetric hybridization in Cordulegaster (Odonata: Cordulegastridae): Secondary postglacial contact and the possible role of mechanical constraints

Two Cordulegaster dragonflies present in Italy, the Palaearctic and northern distributed Cordulegaster boltonii and the endemic to the south of the peninsula Cordulegaster trinacriae, meet in central Italy and give rise to individuals of intermediate morphology. By means of mitochondrial and nuclear markers and of Geometric Morphometrics applied to sexual appendages, we defined i) the geographical boundaries between the two species in Italy and ii) we determined the presence, the extent, and the genetic characteristics of the hybridization. Genetic data evidenced asymmetric hybridization with the males of C. trinacriae able to mate both interspecifically and intraspecifically. The results contrast with expectations under neutral gene introgression and sexual selection. This data, along with the morphological evidence of significant differences in size and shape of sexual appendages between the males of the two species, seem indicative of the role of mechanical constraints in intraspecific matings. The origin of the two species is dated about to 1.32 Mya and the hybridization resulted related to range expansion of the two species after Last Glacial Maximum and this led to the secondary contact between the two taxa in central Italy. At last, our results indicate that the range of C. trinacriae, a threatened and protected species, has been moving northward probably driven by climate changes. As a result, the latter species is currently intruding into the range of C. boltonii. The hybrid area is quite extended and the hybrids seem well adapted to the environment. From a conservation point of view, even if C. trinacriae has a strong genetic identity, the discovery of hybridization between the two species should be considered in a future species management.

Evolutionary insights from comparative transcriptome and transcriptome-wide coalescence analyses in Tetrastigma hemsleyanum

BACKGROUND

Tetrastigma hemsleyanum is of great medicinal importance and used as a model system to address the evolutionary history of warm-temperate evergreen (WTE) forest biomes in East Asia over Neogene time scales. However, further studies on the neutral and adaptive divergence processes of T. hemsleyanum are currently impeded by a lack of genomic resources. In this study, we de novo assembled and annotated a reference transcriptome for two cpDNA lineages (Central-South-East vs. Southwest) of T. hemsleyanum. We further used comparative genomic and multilocus coalescent approaches to investigate the tempo and mode of lineage diversification in T. hemsleyanum.

RESULTS

A total of 52,838 and 65,197 unigenes with an N50 of 1,667 and 1,841 bp for Central-South-East (CSE) and Southwest (SW) lineages, respectively, were recovered, and 6,692 putative orthologs were identified between the two lineages. Estimation of K_a/K_s ratios for these orthologs revealed that ten genes had K_a/K_s values significantly greater than 0.5 (P < 0.05), whereas 2,099 (K_a/K_s < 0.5, P < 0.05) were inferred to be under purifying selection. Based on three bioinformatic strategies, we identified a total of 1,018 single-copy nuclear genes (SCNGs) from the orthologs. We successfully designed eight nuclear gene primer pairs with high intraspecific variation (e.g. h_T = 0.923, π_T = 1.68×10^-3), when surveyed across a subset of T. hemsleyanum individuals. Concordant with the previous cpDNA data, the haplotype networks constructed for most nuclear gene loci clearly identified the two lineages. A multilocus coalescence analysis suggested that the separation between the two lineages appears to have occurred during the mid-Pliocene. Despite their ancient divergence, both lineages experienced expansion at rather localized scales and have continued to exchange genes at a low rate.

CONCLUSIONS

This study demonstrated the utility of transcriptome sequencing as a basis for SCNG development in non-model species and the advantages of integrating multiple nuclear loci for phylogeographic and phylogenetic studies.

DNA barcoding mosquitoes: advice for potential prospectors

Mosquitoes’ importance as vectors of pathogens that drive disease underscores the importance of precise and comparable methods of taxa identification among their species. While several molecular targets have been used to study mosquitoes since the initiation of PCR in the 1980s, its application to mosquito identification took off in the early 1990s. This review follows the research's recent journey into the use of mitochondrial DNA (mtDNA) cytochrome oxidase 1 (COI or COX1) as a DNA barcode target for mosquito species identification – a target whose utility for discriminating mosquitoes is now escalating. The pros and cons of using a mitochondrial genome target are discussed with a broad sweep of the mosquito literature suggesting that nuclear introgressions of mtDNA sequences appear to be uncommon and that the COI works well for distantly related taxa and shows encouraging utility in discriminating more closely related species such as cryptic/sibling species groups. However, the utility of COI in discriminating some closely related groups can be problematic and investigators are advised to proceed with caution as problems with incomplete lineage sorting and introgression events can result in indistinguishable COI sequences appearing in reproductively independent populations. In these – if not all – cases, it is advisable to run a nuclear marker alongside the mtDNA and thus the utility of the ribosomal DNA – and in particular the internal transcribed spacer 2 – is also briefly discussed as a useful counterpoint to the COI.

Comparative phylogeography of Aedes mosquitoes and the role of past climatic change for evolution within Africa

The study of demographic processes involved in species diversification and evolution ultimately provides explanations for the complex distribution of biodiversity on earth, indicates regions important for the maintenance and generation of biodiversity, and identifies biological units important for conservation or medical consequence. African and forest biota have both received relatively little attention with regard to understanding their diversification, although one possible mechanism is that this has been driven by historical climate change. To investigate this, we implemented a standard population genetics approach along with Approximate Bayesian Computation, using sequence data from two exon-primed intron-crossing (EPIC) nuclear loci and mitochondrial cytochrome oxidase subunit I, to investigate the evolutionary history of five medically important and inherently forest dependent mosquito species of the genus Aedes. By testing different demographic hypotheses, we show that Aedes bromeliae and Aedes lilii fit the same model of lineage diversification, admixture, expansion, and recent population structure previously inferred for Aedes aegypti. In addition, analyses of population structure show that Aedes africanus has undergone lineage diversification and expansion while Aedes hansfordi has been impacted by population expansion within Uganda. This congruence in evolutionary history is likely to relate to historical climate-driven habitat change within Africa during the late Pleistocene and Holocene epoch. We find differences in the population structure of mosquitoes from Tanzania and Uganda compared to Benin and Uganda which could relate to differences in the historical connectivity of forests across the continent. Our findings emphasize the importance of recent climate change in the evolution of African forest biota.

Pleistocene climatic oscillations in Neotropical open areas: Refuge isolation in the rodent Oxymycterus nasutus endemic to grasslands

Pleistocene climatic oscillations favoured the expansion of grassland ecosystems and open vegetation landscapes throughout the Neotropics, and influenced the evolutionary history of species adapted to such environments. In this study, we sampled populations of the rodent Oxymycterus nasutus endemic to open areas in the Pampas and Atlantic Forest biomes to assess the tempo and mode of population divergence using an integrative approach, including coalescence theory, ecological niche models, and morphometry. Our results indicated that these O. nasutus populations exhibited high levels of genetic structure. Six major mtDNA clades were found, structuring these biomes into distinct groups. Estimates of their divergence times was indicated to be 0.571 myr. The high degree of genetic structure is reflected in the analyses of geometric morphometric; skull differences between lineages in the two ecoregions were detected. During the last glacial maximum, there was a strong increase in suitable abiotic conditions for O. nasutus. Distinct molecular markers revealed a population expansion over time, with a possible demographic retraction during the post-glacial period. Considering that all clades coalesce with the last interglacial maximum, our results indicated that reduction in suitable conditions during this period may have resulted in a possible vicariance associated with refuge isolation.

Phylogeography of Schizopygopsis stoliczkai (Cyprinidae) in Northwest Tibetan Plateau area

Schizopygopsis stoliczkai (Cyprinidae, subfamily Schizothoracinae) is one of the major freshwater fishes endemic to the northwestern margin of the Tibetan Plateau. In the current study, we used mitochondrial DNA markers cytochrome b (Cyt b) and 16S rRNA (16S), as well as the nuclear marker, the second intron of the nuclear beta-actin gene (Act2), to uncover the phylogeography of S. stoliczkai. In total, we obtained 74 haplotypes from 403 mitochondrial concatenated sequences. The mtDNA markers depict the phylogenetic structures of S. stoliczkai, which consist of clade North and clade South. The split time of the two clades is dated back to 4.27 Mya (95% HPD = 1.96-8.20 Mya). The estimated split time is earlier than the beginning of the ice age of Pleistocene (2.60 Mya), suggesting that the northwestern area of the Tibetan Plateau probably contain at least two glacial refugia for S. stoliczkai. SAMOVA supports the formation of four groups: (i) the Karakash River group; (ii) The Lake Pangong group; (iii) the Shiquan River group; (iv) the Southern Basin group. Clade North included Karakash River, Lake Pangong, and Shiquan River groups, while seven populations of clade South share the haplotypes. Genetic diversity, star-like network, BSP analysis, as well as negative neutrality tests indicate recent expansions events of S. stoliczkai. Conclusively, our results illustrate the phylogeography of S. stoliczkai, implying the Shiquan River is presumably the main refuge for S. stoliczkai.

Mitochondrial and nuclear DNA reveals reticulate evolution in hares (Lepus spp., Lagomorpha, Mammalia) from Ethiopia

For hares (Lepus spp., Leporidae, Lagomorpha, Mammalia) from Ethiopia no conclusive molecular phylogenetic data are available. To provide a first molecular phylogenetic model for the Abyssinian Hare (Lepus habessinicus), the Ethiopian Hare (L. fagani), and the Ethiopian Highland Hare (L. starcki) and their evolutionary relationships to hares from Africa, Eurasia, and North America, we phylogenetically analysed mitochondrial ATPase subunit 6 (ATP6; n = 153 / 416bp) and nuclear transferrin (TF; n = 155 / 434bp) sequences of phenotypically determined individuals. For the hares from Ethiopia, genotype composition at twelve microsatellite loci (n = 107) was used to explore both interspecific gene pool separation and levels of current hybridization, as has been observed in some other Lepus species. For phylogenetic analyses ATP6 and TF sequences of Lepus species from South and North Africa (L. capensis, L. saxatilis), the Anatolian peninsula and Europe (L. europaeus, L. timidus) were also produced and additional TF sequences of 18 Lepus species retrieved from GenBank were included as well. Median joining networks, neighbour joining, maximum likelihood analyses, as well as Bayesian inference resulted in similar models of evolution of the three species from Ethiopia for the ATP6 and TF sequences, respectively. The Ethiopian species are, however, not monophyletic, with signatures of contemporary uni- and bidirectional mitochondrial introgression and/ or shared ancestral polymorphism. Lepus habessinicus carries mtDNA distinct from South African L. capensis and North African L. capensis sensu lato; that finding is not in line with earlier suggestions of its conspecificity with L. capensis. Lepus starcki has mtDNA distinct from L. capensis and L. europaeus, which is not in line with earlier suggestions to include it either in L. capensis or L. europaeus. Lepus fagani shares mitochondrial haplotypes with the other two species from Ethiopia, despite its distinct phenotypic and microsatellite differences; moreover, it is not represented by a species-specific mitochondrial haplogroup, suggesting considerable mitochondrial capture by the other species from Ethiopia or species from other parts of Africa. Both mitochondrial and nuclear sequences indicate close phylogenetic relationships among all three Lepus species from Ethiopia, with L. fagani being surprisingly tightly connected to L. habessinicus. TF sequences suggest close evolutionary relationships between the three Ethiopian species and Cape hares from South and North Africa; they further suggest that hares from Ethiopia hold a position ancestral to many Eurasian and North American species.

The demographic response of a deciduous shrub (the Indigofera bungeana complex, Fabaceae) to the Pleistocene climate changes in East Asia

East Asia harbors the highest level of floristic diversity among the world’s temperate regions. Despite the increase in phylogeographic studies of temperate plants in East Asia, far less attention has been paid to widely distributed deciduous shrubs that widespread across several floral regions. We sequenced two chloroplast DNA (cpDNA) fragments (ndhJ-trnF and trnD-trnT) and one nuclear DNA (Pgk1) of 472 individuals from 51 populations of such a group, the Indigofera bungeana complex. We used population genetic data as well as ecological niche modelling to examine the evolutionary history and glacial refugia during the Last Glacial Maximum (LGM) of this group. We recovered 133 cpDNA and 68 nuclear haplotypes. The star-phylogeny of the recovered cpDNA and nuclear haplotypes and demographic analyses suggested distinct range expansion of I. bungeana complex have occurred during the early and middle Pleistocene. The climate change of the LGM might have affected little on the distribution of this complex based on the niche modelling. However, these climate changes and geographic isolation probably resulted in fixtures of the private haplotypes and genetic differentiations between regions. Our results suggested that this arid-tolerant species complex may have different responses to the Quaternary climate changes with those climate-sensitive species.

Phylogeography of the sergeants Abudefduf sexfasciatus and A. vaigiensis reveals complex introgression patterns between two widespread and sympatric Indo-West Pacific reef fishes

On evolutionary timescales, sea level oscillations lead to recurrent spatio-temporal variation in species distribution and population connectivity. In this situation, applying classical concepts of biogeography is challenging yet necessary to understand the mechanisms underlying biodiversity in highly diverse marine ecosystems such as coral reefs. We aimed at studying the outcomes of such complex biogeographic dynamics on reproductive isolation by sampling populations across a wide spatial range of a species-rich fish genus: the sergeants (Pomacentridae: Abudefduf). We generated a mutlilocus data set that included ten morpho-species from 32 Indo-West Pacific localities. We observed a pattern of mito-nuclear discordance in two common and widely distributed species: Abudefduf sexfasciatus and Abudefduf vaigiensis. The results showed three regional sublineages (Indian Ocean, Coral Triangle region, western Pacific) in A. sexfasciatus (0.6-1.5% divergence at cytb). The other species, A. vaigiensis, is polyphyletic and consists of three distinct genetic lineages (A, B and C) (9% divergence at cytb) whose geographic ranges overlap. Although A. vaigiensis A and A. sexfasciatus were found to be distinct based on nuclear information, A. vaigiensis A was found to be nested within A. sexfasciatus in the mitochondrial gene tree. A. sexfasciatus from the Coral Triangle region and A. vaigiensis A were not differentiated from each other at the mitochondrial locus. We then used coalescent-based simulation to characterize a spatially widespread but weak gene flow between the two species. We showed that these fishes are good candidates to investigate the evolutionary complexity of the discrepancies between phenotypic and genetic similarity in closely related species.

A large historical refugium explains spatial patterns of genetic diversity in a Neotropical savanna tree species

BACKGROUND AND AIMS

The relative role of Pleistocene climate changes in driving the geographic distribution and genetic diversity of South American species is not well known, especially from open biomes such as the Cerrado, the most diverse tropical savanna, encompassing high levels of endemism. Here the effects of Quaternary climatic changes on demographic history, distribution dynamics and genetic diversity of Dimorphandra mollis, an endemic tree species widely distributed in the Cerrado, were investigated.

METHODS

A total of 38 populations covering most of the distribution of D. mollis were analysed using internal transcribed spacer (ITS) sequences and nuclear microsatellite variation [simple sequence repeats (SSRs)]. The framework incorporated statistical phylogeography, coalescent analyses and ecological niche modelling (ENM).

KEY RESULTS

Different signatures of Quaternary climatic changes were found for ITS sequences and SSRs corresponding to different time slices. Coalescent analyses revealed large and constant effective population sizes, with high historical connectivity among the populations for ITS sequences and low effective population sizes and gene flow with recent population retraction for SSRs. ENMs indicated a slight geographical range retraction during the Last Glacial Maximum. A large historical refugium across central Brazil was predicted. Spatially explicit analyses showed a spatial cline pattern in genetic diversity related to the paleodistribution of D. mollis and to the centre of its historical refugium.

CONCLUSIONS

The complex genetic patterns found in D. mollis are the result of a slight geographical range retraction during the Last Glacial Maximum followed by population expansion to the east and south from a large refugium in the central part of the Cerrado. This historical refugium is coincident with an area predicted to be climatically stable under future climate scenarios. The identified refugium should be given high priority in conservation polices to safeguard the evolutionary potential of the species under predicted future climatic changes.

Evolution in African tropical trees displaying ploidy-habitat association: The genus Afzelia (Leguminosae)

Polyploidy has rarely been documented in rain forest trees but it has recently been found in African species of the genus Afzelia (Leguminosae), which is composed of four tetraploid rain forest species and two diploid dry forest species. The genus Afzelia thus provides an opportunity to examine how and when polyploidy and habitat shift occurred in Africa, and whether they are associated. In this study, we combined three plastid markers (psbA, trnL, ndhF), two nuclear markers (ribosomal ITS and the single-copy PEPC E7 gene), plastomes (obtained by High Throughput Sequencing) and morphological traits, with an extensive taxonomic and geographic sampling to explore the evolutionary history of Afzelia. Both nuclear DNA and morphological vegetative characters separated diploid from tetraploid lineages. Although the two African diploid species were well differentiated genetically and morphologically, the relationships among the tetraploid species were not resolved. In contrast to the nuclear markers, plastid markers revealed that one of the diploid species forms a well-supported clade with the tetraploids, suggesting historical hybridisation, possibly in relation with genome duplication (polyploidization) and habitat shift from dry to rain forests. Molecular dating based on fossil-anchored gene phylogenies indicates that extant Afzelia started diverging c. 14.5 or 20Ma while extant tetraploid species started diverging c. 7.0 or 9.4Ma according to plastid and nuclear DNA, respectively. Additional studies of tropical polyploid plants are needed to assess whether the ploidy-habitat association observed in African Afzelia would reflect a role of polyploidization in niche divergence in the tropics.

Multi-locus sequence data illuminate demographic drivers of Pleistocene speciation in semi-arid southern Australian birds (Cinclosoma spp.)

BACKGROUND

During the Pleistocene, shifts of species distributions and their isolation in disjunct refugia led to varied outcomes in how taxa diversified. Some species diverged, others did not. Here, we begin to address another facet of the role of the Pleistocene in generating today's diversity. We ask which processes contributed to divergence in semi-arid southern Australian birds. We isolated 11 autosomal nuclear loci and one mitochondrial locus from a total of 29 specimens of the sister species pair, Chestnut Quail-thrush Cinclosoma castanotum and Copperback Quail-thrush C. clarum.

RESULTS

A population clustering analysis confirmed the location of the current species boundary as a well-known biogeographical barrier in southern Australia, the Eyrean Barrier. Coalescent-based analyses placed the time of species divergence to the Middle Pleistocene. Gene flow between the species since divergence has been low. The analyses suggest the effective population size of the ancestor was 54 to 178 times smaller than populations since divergence. This contrasts with recent multi-locus studies in some other Australian birds (butcherbirds, ducks) where a lack of phenotypic divergence was accompanied by larger historical population sizes. Post-divergence population size histories of C. clarum and C. castanotum were inferred using the extended Bayesian skyline model. The population size of C. clarum increased substantially during the late Pleistocene and continued to increase through the Last Glacial Maximum and Holocene. The timing of this expansion across its vast range is broadly concordant with that documented in several other Australian birds. In contrast, effective population size of C. castanotum was much more constrained and may reflect its smaller range and more restricted habitat east of the Eyrean Barrier compared with that available to C. clarum to the west.

CONCLUSIONS

Our results contribute to awareness of increased population sizes, following significant contractions, as having been important in shaping diversity in Australian arid and semi-arid zones. Further, we improve knowledge of the role of Pleistocene climatic shifts in areas of the planet that were not glaciated at that time but which still experienced that period's cyclical climatic fluctuations.

Colonization and diversification of aquatic insects on three Macaronesian archipelagos using 59 nuclear loci derived from a draft genome

The study of processes driving diversification requires a fully sampled and well resolved phylogeny, although a lack of phylogenetic markers remains a limitation for many non-model groups. Multilocus approaches to the study of recent diversification provide a powerful means to study the evolutionary process, but their application remains restricted because multiple unlinked loci with suitable variation for phylogenetic or coalescent analysis are not available for most non-model taxa. Here we identify novel, putative single-copy nuclear DNA (nDNA) phylogenetic markers to study the colonization and diversification of an aquatic insect species complex, Cloeon dipterum L. 1761 (Ephemeroptera: Baetidae), in Macaronesia. Whole-genome sequencing data from one member of the species complex were used to identify 59 nDNA loci (32,213 base pairs), followed by Sanger sequencing of 29 individuals sampled from 13 islands of three Macaronesian archipelagos. Multispecies coalescent analyses established six putative species. Three island species formed a monophyletic clade, with one species occurring on the Azores, Europe and North America. Ancestral state reconstruction indicated at least two colonization events from the mainland (to the Canaries, respectively Azores) and one within the archipelago (between Madeira and the Canaries). Random subsets of the 59 loci showed a positive linear relationship between number of loci and node support. In contrast, node support in the multispecies coalescent tree was negatively correlated with mean number of phylogenetically informative sites per locus, suggesting a complex relationship between tree resolution and marker variability. Our approach highlights the value of combining genomics, coalescent-based phylogeography, species delimitation, and phylogenetic reconstruction to resolve recent diversification events in an archipelago species complex.

Molecular phylogeny and timing of diversification in Alpine Rhithrogena (Ephemeroptera: Heptageniidae)

BACKGROUND

Larvae of the Holarctic mayfly genus Rhithrogena Eaton, 1881 (Ephemeroptera, Heptageniidae) are a diverse and abundant member of stream and river communities and are routinely used as bio-indicators of water quality. Rhithrogena is well diversified in the European Alps, with a number of locally endemic species, and several cryptic species have been recently detected. While several informal species groups are morphologically well defined, a lack of reliable characters for species identification considerably hampers their study. Their relationships, origin, timing of speciation and mechanisms promoting their diversification in the Alps are unknown.

RESULTS

Here we present a species-level phylogeny of Rhithrogena in Europe using two mitochondrial and three nuclear gene regions. To improve sampling in a genus with many cryptic species, individuals were selected for analysis according to a recent DNA-based taxonomy rather than traditional nomenclature. A coalescent-based species tree and a reconstruction based on a supermatrix approach supported five of the species groups as monophyletic. A molecular clock, mapped on the most resolved phylogeny and calibrated using published mitochondrial evolution rates for insects, suggested an origin of Alpine Rhithrogena in the Oligocene/Miocene boundary. A diversification analysis that included simulation of missing species indicated a constant speciation rate over time, rather than any pronounced periods of rapid speciation. Ancestral state reconstructions provided evidence for downstream diversification in at least two species groups.

CONCLUSIONS

Our species-level analyses of five gene regions provide clearer definitions of species groups within European Rhithrogena. A constant speciation rate over time suggests that the paleoclimatic fluctuations, including the Pleistocene glaciations, did not significantly influence the tempo of diversification of Alpine species. A downstream diversification trend in the hybrida and alpestris species groups supports a previously proposed headwater origin hypothesis for aquatic insects.

discomark: nuclear marker discovery from orthologous sequences using draft genome data

High-throughput sequencing has laid the foundation for fast and cost-effective development of phylogenetic markers. Here we present the program discomark, which streamlines the development of nuclear DNA (nDNA) markers from whole-genome (or whole-transcriptome) sequencing data, combining local alignment, alignment trimming, reference mapping and primer design based on multiple sequence alignments to design primer pairs from input orthologous sequences. To demonstrate the suitability of discomark, we designed markers for two groups of species, one consisting of closely related species and one group of distantly related species. For the closely related members of the species complex of Cloeon dipterum s.l. (Insecta, Ephemeroptera), the program discovered a total of 78 markers. Among these, we selected eight markers for amplification and Sanger sequencing. The exon sequence alignments (2526 base pairs) were used to reconstruct a well-supported phylogeny and to infer clearly structured haplotype networks. For the distantly related species, we designed primers for the insect order Ephemeroptera, using available genomic data from four sequenced species. We developed primer pairs for 23 markers that are designed to amplify across several families. The discomark program will enhance the development of new nDNA markers by providing a streamlined, automated approach to perform genome-scale scans for phylogenetic markers. The program is written in Python, released under a public licence (GNU GPL version 2), and together with a manual and example data set available at: https://github.com/hdetering/discomark.

Local endemism and within-island diversification of shrews illustrate the importance of speciation in building Sundaland mammal diversity

Island systems are important models for evolutionary biology because they provide convenient, discrete biogeographic units of study. Continental islands with a history of intermittent dry land connections confound the discrete definitions of islands and have led zoologists to predict (1) little differentiation of terrestrial organisms among continental shelf islands and (2) extinction, rather than speciation, to be the main cause of differences in community composition among islands. However, few continental island systems have been subjected to well-sampled phylogeographic studies, leaving these biogeographic assumptions of connectivity largely untested. We analyzed nine unlinked loci from shrews of the genus Crocidura from seven mountains and two lowland localities on the Sundaic continental shelf islands of Sumatra and Java. Coalescent species delimitation strongly supported all currently recognized Crocidura species from Sumatra (six species) and Java (five species), as well as one undescribed species endemic to each island. We find that nearly all species of Crocidura in the region are endemic to a single island and several of these have their closest relative(s) on the same island. Intra-island genetic divergence among allopatric, conspecific populations is often substantial, perhaps indicating species-level diversity remains underestimated. One recent (Pleistocene) speciation event generated two morphologically distinct, syntopic species on Java, further highlighting the prevalence of within-island diversification. Our results suggest that both between- and within-island speciation processes generated local endemism in Sundaland, supplementing the traditional view that the region's fauna is relictual and primarily governed by extinction. This article is protected by copyright. All rights reserved.

Variation of partial transferrin sequences and phylogenetic relationships among hares (Lepus capensis, Lagomorpha) from Tunisia

North African hares are currently included in cape hares, Lepus capensis sensu lato, a taxon that may be considered a superspecies or a complex of closely related species. The existing molecular data, however, are not unequivocal, with mtDNA control region sequences suggesting a separate species status and nuclear loci (allozymes, microsatellites) revealing conspecificity of L. capensis and L. europaeus. Here, we study sequence variation in the intron 6 (468 bp) of the transferrin nuclear gene, of 105 hares with different coat colour from different regions in Tunisia with respect to genetic diversity and differentiation, as well as their phylogenetic status. Forty-six haplotypes (alleles) were revealed and compared phylogenetically to all available TF haplotypes of various Lepus species retrieved from GenBank. Maximum Likelihood, neighbor joining and median joining network analyses concordantly grouped all currently obtained haplotypes together with haplotypes belonging to six different Chinese hare species and the African scrub hare L. saxatilis. Moreover, two Tunisian haploypes were shared with L. capensis, L timidus, L. sinensis, L. yarkandensis, and L. hainanus from China. These results indicated the evolutionary complexity of the genus Lepus with the mixing of nuclear gene haplotypes resulting from introgressive hybridization or/and shared ancestral polymorphism. We report the presence of shared ancestral polymorphism between North African and Chinese hares. This has not been detected earlier in the mtDNA sequences of the same individuals. Genetic diversity of the TF sequences from the Tunisian populations was relatively high compared to other hare populations. However, genetic differentiation and gene flow analyses (AMOVA, F_ST, Nm) indicated little divergence with the absence of geographically meaningful phylogroups and lack of clustering with coat colour types. These results confirm the presence of a single hare species in Tunisia, but a sound inference on its phylogenetic position would require additional nuclear markers and numerous geographically meaningful samples from Africa and Eurasia.