Exon capture museomics deciphers the nine-banded armadillo species complex and identifies a new species endemic to the Guiana Shield.

The nine-banded armadillo (Dasypus novemcinctus) is the most widespread xenarthran species across the Americas. Recent studies have suggested it is composed of four morphologically and genetically distinct lineages of uncertain taxonomic status. To address this issue, we used a museomic approach to sequence 80 complete mitogenomes and capture 997 nuclear loci for 71 Dasypus individuals sampled across the entire distribution. We carefully cleaned up potential genotyping errors and cross contaminations that could blur species boundaries by mimicking gene flow. Our results unambiguously support four distinct lineages within the D. novemcinctus complex. We found cases of mito-nuclear phylogenetic discordance but only limited contemporary gene flow confined to the margins of the lineage distributions. All available evidence including the restricted gene flow, phylogenetic reconstructions based on both mitogenomes and nuclear loci, and phylogenetic delimitation methods consistently supported the four lineages within D. novemcinctus as four distinct species. Comparable genetic differentiation values to other recognized Dasypus species further reinforced their status as valid species. Considering congruent morphological results from previous studies, we provide an integrative taxonomic view to recognise four species within the D. novemcinctus complex: D. novemcinctus, D. fenestratus, D. mexicanus, and D. guianensis sp. nov., a new species endemic of the Guiana Shield that we describe here. The two available individuals of D. mazzai and D. sabanicola were consistently nested within D. novemcinctus lineage and their status remains to be assessed. The present work offers a case study illustrating the power of museomics to reveal cryptic species diversity within a widely distributed and emblematic species of mammals.

Prevalence of Genetically Complex Leishmania Strains With Hybrid and Mito-Nuclear Discordance.

Approximately 20 Leishmania species are known to cause cutaneous, mucocutaneous, and visceral disorders in humans. Identification of the causative species in infected individuals is important for appropriate treatment and a favorable prognosis because infecting species are known to be the major determinant of clinical manifestations and may affect treatments for leishmaniasis. Although Leishmania species have been conventionally identified by multilocus enzyme electrophoresis, genetic analysis targeting kinetoplast and nuclear DNA (kDNA and nDNA, respectively) is now widely used for this purpose. Recently, we conducted countrywide epidemiological studies of leishmaniasis in Ecuador and Peru to reveal prevalent species using PCR-RFLP targeting nDNA, and identified unknown hybrid parasites in these countries together with species reported previously. Furthermore, comparative analyses of kDNA and nDNA revealed the distribution of parasites with mismatches between these genes, representing the first report of mito-nuclear discordance in protozoa. The prevalence of an unexpectedly high rate (~10%) of genetically complex strains including hybrid strains, in conjunction with the observation of mito-nuclear discordance, suggests that genetic exchange may occur more frequently than previously thought in natural Leishmania populations. Hybrid Leishmania strains resulting from genetic exchanges are suggested to cause more severe clinical symptoms when compared with parental strains, and to have increased transmissibility by vectors of the parental parasite species. Therefore, it is important to clarify how such genetic exchange influences disease progression and transmissibility by sand flies in nature. In addition, our aim was to identify where and how the genetic exchange resulting in the formation of hybrid and mito-nuclear discordance occurs.

The roles of vicariance and dispersal in the differentiation of two species of the Rhinella marina species complex.

The high levels of Neotropical biodiversity are commonly associated with the intense Neogene-Quaternary geological events and climate dynamics. Here, we investigate the evolutionary history of two species of Neotropical closely related amphibians (R. horribilis and R. marina). We combine published data with new mitochondrial DNA sequences and multiple nuclear markers, including 12 microsatellites. The phylogenetic analyses showed support for grouping the samples in two main clades; R. horribilis (Central America and Mexico) and R. marina (South America east of the Andes). However, the phylogenetic inferences also show an evident mito-nuclear discordance. We use Approximate Bayesian Computation (ABC) to test the role of different events in the diversification between the two groups recovered. We found that both species were affected primarily by a recent Pleistocene divergence, which was similar to the divergence estimate revealed by the Isolation-with-Migration model, under persistent bidirectional gene flow through time. We provide the first evidence that R. horribilis is differentiated from the South American R. marina at the nuclear level supporting the taxonomic status of R. horribilis, which has been controversial for more than a century.

Inferring introgression using RADseq and DFOIL : Power and pitfalls revealed in a case study of spiny lizards (Sceloporus).

Introgression is now commonly reported in studies across the Tree of Life, aided by recent advancements in data collection and analysis. Nevertheless, researchers working with nonmodel species lacking reference genomes may be stymied by a mismatch between available resources and methodological demands. In this study, we demonstrate a fast and simple approach for inferring introgression using RADseq data, and apply it to a case study involving spiny lizards (Sceloporus) from northeastern México. First, we find evidence for recurrent mtDNA introgression between the two focal species based on patterns of mito-nuclear discordance. We then test for nuclear introgression by exhaustively applying the "five-taxon" D-statistic (DFOIL ) to all relevant individuals sampled for RADseq data. In our case, this exhaustive approach (dubbed "ExDFOIL ") entails testing up to ~250,000 unique four-taxon combinations of individuals across species. To facilitate use of this ExDFOIL approach, we provide scripts for many relevant tasks, including the selection of appropriate four-taxon combinations, execution of DFOIL tests in parallel and visualization of introgression results in phylogenetic and geographic space. Using ExDFOIL , we find evidence for ancient introgression between the focal species. Furthermore, we reveal geographic variation in patterns of introgression that is consistent with patterns of mito-nuclear discordance and with recurrent introgression. Overall, our study demonstrates that the combination of DFOIL and RADseq data can effectively detect introgression under a variety of sampling conditions (for individuals, populations and loci). Importantly, we also find evidence that batch-specific error and linkage in RADseq data may mislead inferences of introgression under certain conditions.

Hidden histories of gene flow in highland birds revealed with genomic markers.

Genomic studies are revealing that divergence and speciation are marked by gene flow, but it is not clear whether gene flow has played a prominent role during the generation of biodiversity in species-rich regions of the world where vicariance is assumed to be the principal mode by which new species form. We revisit a well-studied organismal system in the Mexican Highlands, Aphelocoma jays, to test for gene flow among Mexican sierras. Prior results from mitochondrial DNA (mtDNA) largely conformed to the standard model of allopatric divergence, although there was also evidence for more obscure histories of gene flow in a small sample of nuclear markers. We tested for these 'hidden histories' using genomic markers known as ultraconserved elements (UCEs) in concert with phylogenies, clustering algorithms and newer introgression tests specifically designed to detect ancient gene flow (e.g. ABBA/BABA tests). Results based on 4303 UCE loci and 2500 informative SNPs are consistent with varying degrees of gene flow among highland areas. In some cases, gene flow has been extensive and recent (although perhaps not ongoing today), whereas in other cases there is only a trace signature of ancient gene flow among species that diverged as long as 5 million years ago. These results show how a species complex thought to be a model for vicariance can reveal a more reticulate history when a broader portion of the genome is queried. As more organisms are studied with genomic data, we predict that speciation-with-bouts-of-gene-flow will turn out to be a common mode of speciation.

Estimating the temporal and spatial extent of gene flow among sympatric lizard populations (genus Sceloporus) in the southern Mexican highlands.

Interspecific gene flow is pervasive throughout the tree of life. Although detecting gene flow between populations has been facilitated by new analytical approaches, determining the timing and geography of hybridization has remained difficult, particularly for historical gene flow. A geographically explicit phylogenetic approach is needed to determine the overlap of ancestral populations. In this study, we performed population genetic analyses, species delimitation, simulations and a recently developed approach of species tree diffusion to infer the phylogeographic history, timing and geographic extent of gene flow in lizards of the Sceloporus spinosus group. The two species in this group, S. spinosus and S. horridus, are distributed in eastern and western portions of Mexico, respectively, but populations of these species are sympatric in the southern Mexican highlands. We generated data consisting of three mitochondrial genes and eight nuclear loci for 148 and 68 individuals, respectively. We delimited six lineages in this group, but found strong evidence of mito-nuclear discordance in sympatric populations of S. spinosus and S. horridus owing to mitochondrial introgression. We used coalescent simulations to differentiate ancestral gene flow from secondary contact, but found mixed support for these two models. Bayesian phylogeography indicated more than 60% range overlap between ancestral S. spinosus and S. horridus populations since the time of their divergence. Isolation-migration analyses, however, revealed near-zero levels of gene flow between these ancestral populations. Interpreting results from both simulations and empirical data indicate that despite a long history of sympatry among these two species, gene flow in this group has only recently occurred.