Delimiting species without nuclear monophyly in Madagascar's mouse lemurs

Taking the discovery approach in integrative taxonomy: decrypting a complex of narrow-endemic Alpine harvestmen (Opiliones: Phalangiidae: Megabunus)

Species delimitation is fundamental for biological studies, yet precise delimitation is not an easy task, and every involved approach has an inherent failure rate. Integrative taxonomy, a method that merges multiple lines of evidence, can profoundly contribute to reliable alpha-taxonomy and shed light on the processes behind speciation. In this study, we explored and validated species limits in a group of closely related Megabunus harvestmen (Eupnoi, Phalangiidae) endemic to the European Alps. Without a priori species hypotheses, we used multiple sources of inference, including mitochondrial and multilocus nuclear DNA, morphometrics and chemistry. The results of these discovery approaches revealed morphological crypsis and multiple new species within two of the five hitherto known species. Based on our analyses, we discussed the most plausible evolutionary scenarios, invoked the most reasonable species hypotheses and validated the new species limits. Building upon the achieved rigour, three new species, Megabunus cryptobergomas Muster and Wachter sp. nov., Megabunus coelodonta Muster and Steiner sp. nov., and Megabunus lentipes Muster and Komposch sp. nov., are formally described. In addition, we provide a dichotomous morphological key to the Megabunus species of the Alps. Our work demonstrates the suitability of integrative, discovery-based approaches in combination with validation approaches to precisely characterize species and enabled us to implement nomenclatural consequences for this genus.

High-stakes species delimitation in eyeless cave spiders (Cicurina, Dictynidae, Araneae) from central Texas

A remarkable radiation of completely eyeless, cave-obligate spider species (Cicurina) has been described from limestone caves of Texas. This radiation includes over 50 described species, with a large number of hypothesized single-cave endemics, and four species listed as US Federally Endangered. Because of this conservation importance, species delimitation in the group is 'high-stakes'- it is imperative that species hypotheses are data rich, objective, and robust. This study focuses on a complex of four cave-dwelling Cicurina distributed on the northwestern edge of Austin, Texas. Several of the existing species hypotheses in this complex are weak, based on morphological comparisons of small samples of adult female specimens; one species description (for C. wartoni) is based on a single adult specimen. Species limits in this group were newly assessed using morphological, mitochondrial and nuclear DNA sequence data evidence, analysed using a variety of approaches. All data support a clear lineage separation between C. buwata versus the C. travisae complex (including C. travisae, C. wartoni and C. reddelli). Observed congruence across multiple analyses indicate that the C. travisae complex represents a single species, and the formal species synonymy presented here has important conservation implications. The integrative framework utilized in this study serves as a potential model for other Texas cave Cicurina, including US Federally Endangered species. More generally, this study illustrates how and why taxon-focused conservation efforts must prioritize modern species delimitation research (if the existing taxonomy is weak), before devoting precious downstream resources to conservation efforts. The study also highlights the issue of taxonomic type II error that diversity biologists increasingly face as species delimitation moves into the genomics era.

Molecular evolutionary characterization of a V1R subfamily unique to strepsirrhine primates

Vomeronasal receptor genes have frequently been invoked as integral to the establishment and maintenance of species boundaries among mammals due to the elaborate one-to-one correspondence between semiochemical signals and neuronal sensory inputs. Here, we report the most extensive sample of vomeronasal receptor class 1 (V1R) sequences ever generated for a diverse yet phylogenetically coherent group of mammals, the tooth-combed primates (suborder Strepsirrhini). Phylogenetic analysis confirms our intensive sampling from a single V1R subfamily, apparently unique to the strepsirrhine primates. We designate this subfamily as V1Rstrep. The subfamily retains extensive repertoires of gene copies that descend from an ancestral gene duplication that appears to have occurred prior to the diversification of all lemuriform primates excluding the basal genus Daubentonia (the aye-aye). We refer to the descendent clades as V1Rstrep-α and V1Rstrep-β. Comparison of the two clades reveals different amino acid compositions corresponding to the predicted ligand-binding site and thus potentially to altered functional profiles between the two. In agreement with previous studies of the mouse lemur (genus, Microcebus), the majority of V1Rstrep gene copies appear to be intact and under strong positive selection, particularly within transmembrane regions. Finally, despite the surprisingly high number of gene copies identified in this study, it is nonetheless probable that V1R diversity remains underestimated in these nonmodel primates and that complete characterization will be limited until high-coverage assembled genomes are available.

How do you solve a problem like Letharia? A new look at cryptic species in lichen-forming fungi using Bayesian clustering and SNPs from multilocus sequence data

The inclusion of molecular data is increasingly an integral part of studies assessing species boundaries. Analyses based on predefined groups may obscure patterns of differentiation, and population assignment tests provide an alternative for identifying population structure and barriers to gene flow. In this study, we apply population assignment tests implemented in the programs STRUCTURE and BAPS to single nucleotide polymorphisms from DNA sequence data generated for three previous studies of the lichenized fungal genus Letharia. Previous molecular work employing a gene genealogical approach circumscribed six species-level lineages within the genus, four putative lineages within the nominal taxon L. columbiana (Nutt.) J.W. Thomson and two sorediate lineages. We show that Bayesian clustering implemented in the program STRUCTURE was generally able to recover the same six putative Letharia lineages. Population assignments were largely consistent across a range of scenarios, including: extensive amounts of missing data, the exclusion of SNPs from variable markers, and inferences based on SNPs from as few as three gene regions. While our study provided additional evidence corroborating the six candidate Letharia species, the equivalence of these genetic clusters with species-level lineages is uncertain due, in part, to limited phylogenetic signal. Furthermore, both the BAPS analysis and the ad hoc ΔK statistic from results of the STRUCTURE analysis suggest that population structure can possibly be captured with fewer genetic groups. Our findings also suggest that uneven sampling across taxa may be responsible for the contrasting inferences of population substructure. Our results consistently supported two distinct sorediate groups, ‘L. lupina’ and L. vulpina, and subtle morphological differences support this distinction. Similarly, the putative apotheciate species ‘L. lucida’ was also consistently supported as a distinct genetic cluster. However, additional studies will be required to elucidate the relationships of other L. columbiana s.l. populations with the two sorediate genetic clusters.

Multilocus coalescent analyses reveal the demographic history and speciation patterns of mouse lemur sister species

Background

Debate continues as to whether allopatric speciation or peripatric speciation through a founder effect is the predominant force driving evolution in vertebrates. The mouse lemurs of Madagascar are a system in which evolution has generated a large number of species over a relatively recent time frame. Here, we examine speciation patterns in a pair of sister species of mouse lemur, Microcebus murinus and M. griseorufus. These two species have ranges that are disparately proportioned in size, with M. murinus showing a much more extensive range that marginally overlaps that of M. griseorufus. Given that these two species are sister taxa, the asymmetric but overlapping geographic ranges are consistent with a model of peripatric speciation. To test this hypothesis, we analyze DNA sequence data from four molecular markers using coalescent methods. If the peripatric speciation model is supported, we predict substantially greater genetic diversity in M. murinus, relative to M. griseorufus. Further, we expect a larger effective population size in M. murinus and in the common ancestor of the two species than in M. griseorufus, with a concomitant decrease in gene tree/species tree incongruence in the latter and weak signs of demographic expansion in M. murinus.

Results

Our results reject a model of peripatric divergence. Coalescent effective population size estimates were similar for both extant species and larger than that estimated for their most recent common ancestor. Gene tree results show similar levels of incomplete lineage sorting within species with respect to the species tree, and locus-specific estimates of genetic diversity are concordant for both species. Multilocus demographic analyses suggest range expansions for M. murinus, with this species also experiencing more recent population declines over the past 160 thousand years.

Conclusions

Results suggest that speciation occurred in allopatry from a common ancestor narrowly distributed throughout southwest Madagascar, with subsequent range expansion for M. murinus. Population decline in M. murinus is likely related to patterns of climate change in Madagascar throughout the Pleistocene, potentially exacerbated by continual anthropogenic perturbation. Genome-level data are needed to quantify the role of niche specialization and adaptation in shaping the current ranges of these species.

A phylogeographic evaluation of the Amolops mantzorum species group: cryptic species and plateau uplift

Phylogeographic analysis brings organism phylogeny, regional geological/climatic history, and population demography together, and provides insights into species history and speciation. Using DNA sequence data from a mitochondrial marker (cytochrome b) and a nuclear marker (cmyc intron 2), we examined phylogeography of the Amolops mantzorum species group. We first constructed Bayesian and maximum-likelihood gene trees and medium-joint networks for the recovered haplotypes, and estimated divergence time for each major lineage. Species boundaries were then delineated using the general mixed Yule-coalescent model (GMYC) and a Bayesian species-delimitation method (BP&P). Potential gene flow between putative species was also estimated using the isolation-with-migration model. Furthermore, species-distribution modeling was used to probe linkage between climatic changes and speciation. Lastly, extended Bayesian skyline plotting was employed to reveal historical demography for each putative species. Our analyses clearly delimit nine species in the group, including four well-recognized species and five putative species, of which three are often listed as synonyms of A. mantzorum. The nominal species A. mantzorum may in fact include two cryptic species. Divergence-time estimates align the speciation events with the recent intense uplift of the Tibetan Plateau in the last 3.6 million years. The species-distribution modeling also reveals different habitat preferences among species that are potentially linked to climatic changes associated with the uplift. Furthermore, three species sustained continuous population growth through the last glaciation, while others expanded only after the last glaciation. The eastern escarpment of Tibet is a biodiversity hotspot; its rich species diversity, unique topography, and well-established geological history provide an excellent system for evolutionary studies.

Pleistocene speciation in North American lichenized fungi and the impact of alternative species circumscriptions and rates of molecular evolution on divergence estimates

Pleistocene climatic fluctuations influenced patterns of genetic variation and promoted speciation across a wide range of species groups. Lichens are commonly found in habitats that were directly impacted by glacial cycles; however, the role of Pleistocene climate in driving speciation in most lichen symbionts remains unclear. This uncertainty is due in part to limitations in our ability to accurately recognize independently evolving lichen-forming fungal lineages and a lack of relevant fossil calibrations. Using a coalescent-based species tree approach, we estimated divergence times for two sister clades in the genus Xanthoparmelia (Parmeliaceae) restricted to western North America. We assessed the influence of two different species circumscription scenarios and various locus-specific rates of molecular evolution on divergence estimates. Species circumscriptions were validated using the program BP&P. although speciation was generally supported in both scenarios, divergence times differed between traditional species circumscriptions and those based on genetic data, with more recent estimates resulting from the former. Similarly, rates of evolution for different loci resulted in variable divergence time estimates. However, our results unambiguously indicate that diversification in the sampled Xanthoparmelia clades occurred during the Pleistocene. Our study highlights the potential impact of ambiguous species circumscriptions and uncertain rates of molecular evolution on estimating divergence times within a multilocus species tree framework.

Phylogeographic analysis of the true lemurs (genus Eulemur) underlines the role of river catchments for the evolution of micro-endemism in Madagascar

INTRODUCTION

Due to its remarkable species diversity and micro-endemism, Madagascar has recently been suggested to serve as a biogeographic model region. However, hypothesis-based tests of various diversification mechanisms that have been proposed for the evolution of the island's micro-endemic lineages are still limited. Here, we test the fit of several diversification hypotheses with new data on the broadly distributed genus Eulemur using coalescent-based phylogeographic analyses.

RESULTS

Time-calibrated species tree analyses and population genetic clustering resolved the previously polytomic species relationships among eulemurs. The most recent common ancestor of eulemurs was estimated to have lived about 4.45 million years ago (mya). Divergence date estimates furthermore suggested a very recent diversification among the members of the "brown lemur complex", i.e. former subspecies of E. fulvus, during the Pleistocene (0.33-1.43 mya). Phylogeographic model comparisons of past migration rates showed significant levels of gene flow between lineages of neighboring river catchments as well as between eastern and western populations of the redfronted lemur (E. rufifrons).

CONCLUSIONS

Together, our results are concordant with the centers of endemism hypothesis (Wilmé et al. 2006, Science 312:1063-1065), highlight the importance of river catchments for the evolution of Madagascar's micro-endemic biota, and they underline the usefulness of testing diversification mechanisms using coalescent-based phylogeographic methods.

Cryptic genetic diversity is paramount in small-bodied amphibians of the genus Euparkerella (Anura: Craugastoridae) endemic to the Brazilian Atlantic forest

Morphological similarity associated to restricted distributions and low dispersal abilities make the direct developing “Terrarana” frogs of the genus Euparkerella a good model for examining diversification processes. We here infer phylogenetic relationships within the genus Euparkerella, using DNA sequence data from one mitochondrial and four nuclear genes coupled with traditional Bayesian phylogenetic reconstruction approaches and more recent coalescent methods of species tree inference. We also used Bayesian clustering analysis and a recent Bayesian coalescent-based approach specifically to infer species delimitation. The analysis of 39 individuals from the four known Euparkerella species uncovered high levels of genetic diversity, especially within the two previously morphologically-defined E. cochranae and E. brasiliensis. Within these species, the gene trees at five independent loci and trees from combined data (concatenated dataset and the species tree) uncovered six deeply diverged and geographically coherent evolutionary units, which may have diverged between the Miocene and the Pleistocene. These six units were also uncovered in the Bayesian clustering analysis, and supported by the Bayesian coalescent-based species delimitation (BPP), and Genealogical Sorting Index (GSI), providing thus strong evidence for underestimation of the current levels of diversity within Euparkerella. The cryptic diversity now uncovered opens new opportunities to examine the origins and maintenance of microendemism in the context of spatial heterogeneity and/or human induced fragmentation of the highly threatened Brazilian Atlantic forest hotspot.

True lemurs…true species - species delimitation using multiple data sources in the brown lemur complex

Background

Species are the fundamental units in evolutionary biology. However, defining them as evolutionary independent lineages requires integration of several independent sources of information in order to develop robust hypotheses for taxonomic classification. Here, we exemplarily propose an integrative framework for species delimitation in the “brown lemur complex” (BLC) of Madagascar, which consists of seven allopatric populations of the genus Eulemur (Primates: Lemuridae), which were sampled extensively across northern, eastern and western Madagascar to collect fecal samples for DNA extraction as well as recordings of vocalizations. Our data base was extended by including museum specimens with reliable identification and locality information for skull shape and pelage color analysis.

Results

Between-group analyses of principal components revealed significant heterogeneity in skull shape, pelage color variation and loud calls across all seven populations. Furthermore, post-hoc statistical tests between pairs of populations revealed considerable discordance among different data sets for different dyads. Despite a high degree of incomplete lineage sorting among nuclear loci, significant exclusive ancestry was found for all populations, except for E. cinereiceps, based on one mitochondrial and three nuclear genetic loci.

Conclusions

Using several independent lines of evidence, our results confirm the species status of the members of the BLC under the general lineage concept of species. More generally, the present analyses demonstrate the importance and value of integrating different kinds of data in delimiting recently evolved radiations.

How to fail at species delimitation

Species delimitation is the act of identifying species-level biological diversity. In recent years, the field has witnessed a dramatic increase in the number of methods available for delimiting species. However, most recent investigations only utilize a handful (i.e. 2-3) of the available methods, often for unstated reasons. Because the parameter space that is potentially relevant to species delimitation far exceeds the parameterization of any existing method, a given method necessarily makes a number of simplifying assumptions, any one of which could be violated in a particular system. We suggest that researchers should apply a wide range of species delimitation analyses to their data and place their trust in delimitations that are congruent across methods. Incongruence across the results from different methods is evidence of either a difference in the power to detect cryptic lineages across one or more of the approaches used to delimit species and could indicate that assumptions of one or more of the methods have been violated. In either case, the inferences drawn from species delimitation studies should be conservative, for in most contexts it is better to fail to delimit species than it is to falsely delimit entities that do not represent actual evolutionary lineages.

New metrics for comparison of taxonomies reveal striking discrepancies among species delimitation methods in Madascincus lizards

Delimiting and describing species is fundamental to numerous biological disciplines such as evolution, macroecology, and conservation. Delimiting species as independent evolutionary lineages may and often does yield different outcomes depending on the species criteria applied, but methods should be chosen that minimize the inference of objectively erroneous species limits. Several protocols exploit single-gene or multi-gene coalescence statistics, assignment tests or other rationales related to nuclear DNA (nDNA) allele sharing to automatically delimit species. We apply seven different species delimitation protocols to a taxonomically confusing group of Malagasy lizards (Madascincus), and compare the resulting taxonomies with two newly developed metrics: the Taxonomic index of congruence C tax which quantifies the congruence between two taxonomies, and the Relative taxonomic resolving power index R tax which quantifies the potential of an approach to capture a high number of species boundaries. The protocols differed in the total number of species proposed, between 9 and 34, and were also highly incongruent in placing species boundaries. The Generalized Mixed Yule-Coalescent approach captured the highest number of potential species boundaries but many of these were clearly contradicted by extensive nDNA admixture between sympatric mitochondrial DNA (mtDNA) haplotype lineages. Delimiting species as phenotypically diagnosable mtDNA clades failed to detect two cryptic species that are unambiguous due to a lack of nDNA gene flow despite sympatry. We also consider the high number of species boundaries and their placement by multi-gene Bayesian species delimitation as poorly reliable whereas the Bayesian assignment test approach provided a species delimitation highly congruent with integrative taxonomic practice. The present study illustrates the trade-off in taxonomy between reliability (favored by conservative approaches) and resolving power (favored by inflationist approaches). Quantifying excessive splitting is more difficult than quantifying excessive lumping, suggesting a priority for conservative taxonomies in which errors are more liable to be detected and corrected by subsequent studies.

Multilocus species delimitation in a complex of morphologically conserved trapdoor spiders (mygalomorphae, antrodiaetidae, aliatypus)

Species are a fundamental unit for biological studies, yet no uniform guidelines exist for determining species limits in an objective manner. Given the large number of species concepts available, defining species can be both highly subjective and biased. Although morphology has been commonly used to determine species boundaries, the availability and prevalence of genetic data has allowed researchers to use such data to make inferences regarding species limits. Genetic data also have been used in the detection of cryptic species, where other lines of evidence (morphology in particular) may underestimate species diversity. In this study, we investigate species limits in a complex of morphologically conserved trapdoor spiders (Mygalomorphae, Antrodiaetidae, Aliatypus) from California. Multiple approaches were used to determine species boundaries in this highly genetically fragmented group, including both multilocus discovery and validation approaches (plus a chimeric approach). Additionally, we introduce a novel tree-based discovery approach using species trees. Results suggest that this complex includes multiple cryptic species, with two groupings consistently recovered across analyses. Due to incongruence across analyses for the remaining samples, we take a conservative approach and recognize a three species complex, and formally describe two new species (Aliatypus roxxiae, sp. nov. and Aliatypus starretti, sp. nov.). This study helps to clarify species limits in a genetically fragmented group and provides a framework for identifying and defining the cryptic lineage diversity that prevails in many organismal groups.

Genetic diversity, phylogeny and evolution of alkaloid sequestering in Cuban miniaturized frogs of the Eleutherodactylus limbatus group

The miniaturized frogs of the Eleutherodactylus limbatus group, endemic to Cuba, have recently been shown to contain alkaloids in their skin, representing the fifth independent evolutionary origin of alkaloid sequestering in anurans. Based on a comprehensive sampling of all known species of the group we here assess their range-wide phylogeography using DNA sequences of the mitochondrial cytochrome b and the nuclear Rag-1 gene. We found E. etheridgei, E. cubanus, and E. orientalis to be genetically highly distinct, with uncorrected pair-wise distances >12% in the cytochrome b gene and >5% in the 16S rRNA gene, and without haplotype sharing in the Rag-1 gene. On the contrary, E. iberia, E. jaumei, E. limbatus, and one unnamed lineage from Gran Piedra mountains had lower mitochondrial divergences to each other (<6% in cytochrome b; <3% in 16S) and extensive Rag-1 haplotype sharing, suggesting that their species status requires careful revision. All species and lineages of the group occur in the mountain areas of eastern Cuba and only one of them, E. limbatus, colonized central and western Cuba. Populations of this latter lineage are only weakly genetically differentiated, and its range expansion to western Cuba is probably Pleistocenic, demonstrating that also miniaturized frog species are able of considerable range expansions (over hundreds of kilometers) in short time spans. A phylogeny based on about 5 kb of mitochondrial and nuclear genes places, with maximum support, E. etheridgei basalmost in the group, followed by E. cubanus, E. orientalis, and the E. iberia-jaumei-limbatus clade. The screening of skin alkaloids revealed the presence of alkaloids in all lineages of the group, but alkaloids known to be sequestered by other anurans were not found in the basalmost E. etheridgei. Furthermore, individuals of the derived lineages of the E. iberia-jaumei-limbatus clade had a much larger variety of different alkaloid compounds (3-18) than the other three species (0-4). Traces of alkaloids were found in other species only remotely related to the E. limbatus group, i.e., E. caspari, E. goini, E. ricordii, and E. tetajulia. According to these findings, it can be hypothesized that the capacity to sequester dietary alkaloids into the skin evolved first in an ancestral, litter-dwelling and small-sized Cuban Eleutherodactylus, followed by evolution of miniaturization and diurnality in the early representatives of the E. limbatus group, and eventually by increased brightness of color with probable aposematic function and a highly efficient mechanism of alkaloid sequestering.

An approach to identify putative hybrids in the 'coalescent stochasticity zone', as exemplified in the African plant genus Streptocarpus (Gesneriaceae)

The inference of phylogenetic relationships is often complicated by differing evolutionary histories of independently-inherited markers. The causes of the resulting gene tree incongruence can be challenging to identify, often relying on coalescent simulations dependent on unverifiable assumptions. We investigated alternative techniques using the South African rosulate species of Streptocarpus as a study group. Two independent gene trees - from the nuclear ITS region and from three concatenated plastid regions (trnL-F, rpl20-rps12 and trnC-D) - displayed widespread, strongly supported incongruence. We investigated the causes by detecting genetic exchange across morphological borders using morphological optimizations and genetic exchange across species boundaries using the genealogical sorting index. Incongruence between gene trees was associated with ancestral shifts in growth form (in four species) but not in pollination syndrome, suggesting introgression limited by reproductive barriers. Genealogical sorting index calculations showed polyphyly of two additional species, while individuals of all others were significantly associated. In one case the association was stronger according to the internal transcribed spacer data than according to the plastid data, which, given the smaller effective population size of the plastid, may also indicate introgression. These approaches offer alternative ways to identify potential hybridization events where incomplete lineage sorting cannot be rejected using simulations.

Macroevolutionary dynamics and historical biogeography of primate diversification inferred from a species supermatrix

Phylogenetic relationships, divergence times, and patterns of biogeographic descent among primate species are both complex and contentious. Here, we generate a robust molecular phylogeny for 70 primate genera and 367 primate species based on a concatenation of 69 nuclear gene segments and ten mitochondrial gene sequences, most of which were extracted from GenBank. Relaxed clock analyses of divergence times with 14 fossil-calibrated nodes suggest that living Primates last shared a common ancestor 71-63 Ma, and that divergences within both Strepsirrhini and Haplorhini are entirely post-Cretaceous. These results are consistent with the hypothesis that the Cretaceous-Paleogene mass extinction of non-avian dinosaurs played an important role in the diversification of placental mammals. Previous queries into primate historical biogeography have suggested Africa, Asia, Europe, or North America as the ancestral area of crown primates, but were based on methods that were coopted from phylogeny reconstruction. By contrast, we analyzed our molecular phylogeny with two methods that were developed explicitly for ancestral area reconstruction, and find support for the hypothesis that the most recent common ancestor of living Primates resided in Asia. Analyses of primate macroevolutionary dynamics provide support for a diversification rate increase in the late Miocene, possibly in response to elevated global mean temperatures, and are consistent with the fossil record. By contrast, diversification analyses failed to detect evidence for rate-shift changes near the Eocene-Oligocene boundary even though the fossil record provides clear evidence for a major turnover event ("Grande Coupure") at this time. Our results highlight the power and limitations of inferring diversification dynamics from molecular phylogenies, as well as the sensitivity of diversification analyses to different species concepts.

Parallel tagged amplicon sequencing reveals major lineages and phylogenetic structure in the North American tiger salamander (Ambystoma tigrinum) species complex

Modern analytical methods for population genetics and phylogenetics are expected to provide more accurate results when data from multiple genome-wide loci are analysed. We present the results of an initial application of parallel tagged sequencing (PTS) on a next-generation platform to sequence thousands of barcoded PCR amplicons generated from 95 nuclear loci and 93 individuals sampled across the range of the tiger salamander (Ambystoma tigrinum) species complex. To manage the bioinformatic processing of this large data set (344 330 reads), we developed a pipeline that sorts PTS data by barcode and locus, identifies high-quality variable nucleotides and yields phased haplotype sequences for each individual at each locus. Our sequencing and bioinformatic strategy resulted in a genome-wide data set with relatively low levels of missing data and a wide range of nucleotide variation. structure analyses of these data in a genotypic format resulted in strongly supported assignments for the majority of individuals into nine geographically defined genetic clusters. Species tree analyses of the most variable loci using a multi-species coalescent model resulted in strong support for most branches in the species tree; however, analyses including more than 50 loci produced parameter sampling trends that indicated a lack of convergence on the posterior distribution. Overall, these results demonstrate the potential for amplicon-based PTS to rapidly generate large-scale data for population genetic and phylogenetic-based research.

Hiding deep in the trees: discovery of divergent mitochondrial lineages in Malagasy chameleons of the Calumma nasutum group

We conducted a comprehensive molecular phylogenetic study for a group of chameleons from Madagascar (Chamaeleonidae: Calumma nasutum group, comprising seven nominal species) to examine the genetic and species diversity in this widespread genus. Based on DNA sequences of the mitochondrial gene (ND2) from 215 specimens, we reconstructed the phylogeny using a Bayesian approach. Our results show deep divergences among several unnamed mitochondrial lineages that are difficult to identify morphologically. We evaluated lineage diversification using a number of statistical phylogenetic methods (general mixed Yule-coalescent model; SpeciesIdentifier; net p-distances) to objectively delimit lineages that we here consider as operational taxonomic units (OTUs), and for which the taxonomic status remains largely unknown. In addition, we compared molecular and morphological differentiation in detail for one particularly diverse clade (the C. boettgeri complex) from northern Madagascar. To assess the species boundaries within this group we used an integrative taxonomic approach, combining evidence from two independent molecular markers (ND2 and CMOS), together with genital and other external morphological characters, and conclude that some of the newly discovered OTUs are separate species (confirmed candidate species, CCS), while others should best be considered as deep conspecific lineages (DCLs). Our analysis supports a total of 33 OTUs, of which seven correspond to described species, suggesting that the taxonomy of the C. nasutum group is in need of revision.

The lemur revolution starts now: the genomic coming of age for a non-model organism

Morris Goodman was a revolutionary. Together with a mere handful of like-minded scientists, Morris established himself as a leader in the molecular phylogenetic revolution of the 1960s. The effects of this revolution are most evident in this journal, which he founded in 1992. Happily for lemur biologists, one of Morris Goodman's primary interests was in reconstructing the phylogeny of the primates, including the tooth-combed Lorisifomes of Africa and Asia, and the Lemuriformes of Madagascar (collectively referred to as the suborder Strepsirrhini). This paper traces the development of molecular phylogenetic and evolutionary genetic trends and methods over the 50-year expanse of Morris Goodman's career, particularly as they apply to our understanding of lemuriform phylogeny, biogeography, and biology. Notably, this perspective reveals that the lemuriform genome is sufficiently rich in phylogenetic signal such that the very earliest molecular phylogenetic studies - many of which were conducted by Goodman himself - have been validated by contemporary studies that have exploited advanced computational methods applied to phylogenomic scale data; studies that were beyond imagining in the earliest days of phylogeny reconstruction. Nonetheless, the frontier still beckons. New technologies for gathering and analyzing genomic data will allow investigators to build upon what can now be considered a nearly-known phylogeny of the Lemuriformes in order to ask innovative questions about the evolutionary mechanisms that generate and maintain the extraordinary breadth and depth of biological diversity within this remarkable clade of primates.

PhyloMarker–-A Tool for Mining Phylogenetic Markers through Genome Comparison: Application of the Mouse Lemur (Genus Microcebus) Phylogeny

Molecular phylogeny is a fundamental tool to understanding the evolution of all life forms. One common issue faced by molecular phylogeny is the lack of sufficient molecular markers. Here, we present PhyloMarker, a phylogenomic tool designed to find nuclear gene markers for the inference of phylogeny through multiple genome comparison. Around 800 candidate markers were identified by PhyloMarker through comparison of partial genomes of Microcebus and Otolemur. In experimental tests of 20 randomly selected markers, nine markers were successfully amplified by PCR and directly sequenced in all 17 nominal Microcebus species. Phylogenetic analyses of the sequence data obtained for 17 taxa and nine markers confirmed the distinct lineage inferred from previous mtDNA data. PhyloMarker has also been used by other projects including the herons (Ardeidae, Aves) phylogeny and the Wood mice (Muridae, Mammalia) phylogeny. All source code and sample data are made available at http://bioinfo-srv1.awh.unomaha.edu/phylomarker /.

Pervasive and ongoing positive selection in the vomeronasal-1 receptor (V1R) repertoire of mouse lemurs

Chemosensory genes are frequently the target of positive selection and are often present in large gene families, but little is known about heterogeneity of selection in these cases and its relation to function. Here, we use the vomeronasal-1 receptor (V1R) repertoire of mouse lemurs (Microcebus spp.) as a model system to study patterns of selection of chemosensory genes at several different levels. Mouse lemurs are small nocturnal strepsirrhine primates and have a large (~200 loci) repertoire of V1R loci that are likely important for intraspecific pheromonal communication and interspecific interactions, for example, recognition of predator cues. We investigated signals and patterns of positive selection among the 105 identified full length V1R loci in the gray mouse lemur and within 7 V1R loci amplified across multiple mouse lemur species. Phylogenetic reconstructions of published sequences revealed at least nine monophyletic clusters of V1Rs in gray mouse lemurs that have diversified since the split between lemurs and lorisoid primates. A large majority of clusters evolved under significant positive selection. Similar results were found in V1Rs of closely related greater galagos. Comparison with function of related V1R clusters in mice suggested a potential relationship between receptor function and strength of selection. Interestingly, most codons identified as being under positive selection are located in the extracellular domains of the receptors and hence likely indicate the position of residues involved in ligand binding. Positive selection was also detected within five V1R loci (=71% of analyzed loci) sequenced from 6 to 10 mouse lemur species, indicating ongoing selection within the genus, which may be related to sexual selection and, potentially, speciation processes. Variation in strength of positive selection on V1Rs showed no simple relationship to cluster size. The diversity of V1R loci in mouse lemurs reflects their adaptive evolution and is most likely related to the fundamental relevance of olfactory communication and predator recognition in these primates. Overall, adaptive evolution is the predominant mode of evolution of V1R loci at all levels, and the substantial heterogeneity in the strength of selection may be related to receptor function.

Integration of molecular, ecological, morphological and endosymbiont data for species delimitation within the Pnigalio soemius complex (Hymenoptera: Eulophidae)

Integrative taxonomy is a recently developed approach that uses multiple lines of evidence such as molecular, morphological, ecological and geographical data to test species limits, and it stands as one of the most promising approaches to species delimitation in taxonomically difficult groups. The Pnigalio soemius complex (Hymenoptera: Eulophidae) represents an interesting taxonomical and ecological study case, as it is characterized by a lack of informative morphological characters, deep mitochondrial divergence, and is susceptible to infection by parthenogenesis-inducing Rickettsia. We tested the effectiveness of an integrative taxonomy approach in delimiting species within the P. soemius complex. We analysed two molecular markers (COI and ITS2) using different methods, performed multivariate analysis on morphometric data and exploited ecological data such as host-plant system associations, geographical separation, and the prevalence, type and effects of endosymbiont infection. The challenge of resolving different levels of resolution in the data was met by setting up a formal procedure of data integration within and between conflicting independent lines of evidence. An iterative corroboration process of multiple sources of data eventually indicated the existence of several cryptic species that can be treated as stable taxonomic hypotheses. Furthermore, the integrative approach confirmed a trend towards host specificity within the presumed polyphagous P. soemius and suggested that Rickettsia could have played a major role in the reproductive isolation and genetic diversification of at least two species.

Concatenation and concordance in the reconstruction of mouse lemur phylogeny: an empirical demonstration of the effect of allele sampling in phylogenetics

The systematics and speciation literature is rich with discussion relating to the potential for gene tree/species tree discordance. Numerous mechanisms have been proposed to generate discordance, including differential selection, long-branch attraction, gene duplication, genetic introgression, and/or incomplete lineage sorting. For speciose clades in which divergence has occurred recently and rapidly, recovering the true species tree can be particularly problematic due to incomplete lineage sorting. Unfortunately, the availability of multilocus or "phylogenomic" data sets does not simply solve the problem, particularly when the data are analyzed with standard concatenation techniques. In our study, we conduct a phylogenetic study for a nearly complete species sample of the dwarf and mouse lemur clade, Cheirogaleidae. Mouse lemurs (genus, Microcebus) have been intensively studied over the past decade for reasons relating to their high level of cryptic species diversity, and although there has been emerging consensus regarding the evolutionary diversity contained within the genus, there is no agreement as to the inter-specific relationships within the group. We attempt to resolve cheirogaleid phylogeny, focusing especially on the mouse lemurs, by employing a large multilocus data set. We compare the results of Bayesian concordance methods with those of standard gene concatenation, finding that though concatenation yields the strongest results as measured by statistical support, these results are found to be highly misleading. By employing an approach where individual alleles are treated as operational taxonomic units, we show that phylogenetic results are substantially influenced by the selection of alleles in the concatenation process.

Delimiting species using multilocus data: diagnosing cryptic diversity in the southern cavefish, Typhlichthys subterraneus (Teleostei: Amblyopsidae)

A major challenge facing biodiversity conservation and management is that a significant portion of species diversity remains undiscovered or undescribed. This is particularly evident in subterranean animals in which species delimitation based on morphology is difficult because differentiation is often obscured by phenotypic convergence. Multilocus genetic data constitute a valuable source of information for species delimitation in such organisms, but until recently, few methods were available to objectively test species delimitation hypotheses using genetic data. Here, we use recently developed methods for discovering and testing species boundaries and relationships using a multilocus dataset in a widely distributed subterranean teleost fish, Typhlichthys subterraneus, endemic to Eastern North America. We provide evidence that species diversity in T. subterraneus is currently underestimated and that the picture of a single, widely distributed species is not supported. Rather, several morphologically cryptic lineages comprise the diversity in this clade, including support for the recognition of T. eigenmanni. The high number of cryptic species in Typhlichthys highlights the utility of multilocus genetic data in delimiting species, particularly in lineages that exhibit slight morphological disparity, such as subterranean organisms. However, results depend on sampling of individuals and loci; this issue needs further study.

Hybridization of mouse lemurs: different patterns under different ecological conditions

Background

Several mechanistic models aim to explain the diversification of the multitude of endemic species on Madagascar. The island's biogeographic history probably offered numerous opportunities for secondary contact and subsequent hybridization. Existing diversification models do not consider a possible role of these processes. One key question for a better understanding of their potential importance is how they are influenced by different environmental settings. Here, we characterized a contact zone between two species of mouse lemurs, Microcebus griseorufus and M. murinus, in dry spiny bush and mesic gallery forest that border each other sharply without intermediate habitats between them. We performed population genetic analyses based on mtDNA sequences and nine nuclear microsatellites and compared the results to a known hybrid zone of the same species in a nearby wide gradient from dry spiny bush over transitional forest to humid littoral forest.

Results

In the spiny-gallery system, Microcebus griseorufus is restricted to the spiny bush; Microcebus murinus occurs in gallery forest and locally invades the dryer habitat of its congener. We found evidence for bidirectional introgressive hybridization, which is closely linked to increased spatial overlap within the spiny bush. Within 159 individuals, we observed 18 hybrids with mitochondrial haplotypes of both species. Analyses of simulated microsatellite data indicate that we identified hybrids with great accuracy and that we probably underestimated their true number. We discuss short-term climatic fluctuations as potential trigger for the dynamic of invasion and subsequent hybridization. In the gradient hybrid zone in turn, long-term aridification could have favored unidirectional nuclear introgression from Microcebus griseorufus into M. murinus in transitional forest.

Conclusions

Madagascar's southeastern transitional zone harbors two very different hybrid zones of mouse lemurs in different environmental settings. This sheds light on the multitude of opportunities for the formation of hybrid zones and indicates an important influence of environmental factors on secondary contact and hybridization. Our findings suggest that hybridization could enhance the adaptability of mouse lemurs without necessarily leading to a loss of distinctiveness. They point to a potential role of hybridization in Madagascar's diversification history that requires further investigation.

Isolation-driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae)

Lineage, or true 'species', trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene-tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single-locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re-evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent-based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.

Multilocus phylogeny and species delimitation within the Natterer's bat species complex in the Western Palearctic

Delimiting species is a crucial issue for many biological disciplines and is of primary importance for designing effective conservation plans. Traditional taxonomy based on morphological characters can be misled by the presence of phenotypic plesiomorphism or adaptative convergence. The use of multiple locus genetic data appears thus as a powerful tool for recognizing species boundaries. In this study, we used six nuclear introns and two mitochondrial markers to conduct a phylogenetic study of the Myotis nattereri species complex in the Western Palearctic. We combined tree-based and non-tree-based analyses, and also used concatenated phylogenetic methods of the separated nuclear and mitochondrial dataset as well as a recent coalescence-based multilocus approach. The strong concordance between the results of the analyses conducted confirms that M. nattereri is a paraphyletic group that is composed of four well-differentiated lineages in the study area. In the framework of the unified species concept, these four clades can be confidently considered as four valid species. This recognition of new cryptic species in the Western Mediterranean region shows that the biodiversity of this well-studied area is still not fully understood.

On species delimitation: yet another lemur species or just genetic variation?

Background

Although most taxonomists agree that species are independently evolving metapopulation lineages that should be delimited with several kinds of data, the taxonomic practice in Malagasy primates (Lemuriformes) looks quite different. Several recently described lemur species are based solely on evidence of genetic distance and diagnostic characters of mitochondrial DNA sequences sampled from a few individuals per location. Here we explore the validity of this procedure for species delimitation in lemurs using published sequence data.

Results

We show that genetic distance estimates and Population Aggregation Analysis (PAA) are inappropriate for species delimitation in this group of primates. Intra- and interspecific genetic distances overlapped in 14 of 17 cases independent of the genetic marker used. A simulation of a fictive taxonomic study indicated that for the mitochondrial D-loop the minimum required number of individuals sampled per location is 10 in order to avoid false positives via PAA.

Conclusions

Genetic distances estimates and PAA alone should not be used for species delimitation in lemurs. Instead, several nuclear and sex-specific loci should be considered and combined with other data sets from morphology, ecology or behavior. Independent of the data source, sampling should be done in a way to ensure a quantitative comparison of intra- and interspecific variation of the taxa in question. The results of our study also indicate that several of the recently described lemur species should be reevaluated with additional data and that the number of good species among the currently known taxa is probably lower than currently assumed.