Multilocus phylogeny of alligator lizards ( Elgaria , Anguidae): Testing mtDNA introgression as the source of discordant molecular phylogenetic hypotheses

An updated list of the Mexican herpetofauna: with a summary of historical and contemporary studies

The growth in our knowledge of the diversity of the herpetofauna of Mexico has occurred over the period of approximately 445 years from the work of Francisco Hernández to that of a broad multinational array of present-day herpetologists. The work of this huge group of people has established Mexico as one of the most significant centers of herpetofaunal biodiversity in the world. This status is the result of a complex orography, in addition to diverse habitats and environments and the biogeographic history of Mexico. The current herpetofauna consists of 1,421 native and introduced species, allocated to 220 genera, and 61 families. This figure is comprised of 1,405 native species and 16 non-native species (as of April 2023). The non-native species include two anurans, 13 squamates, and one turtle. The level of endemism is very high, presently lying at 63%, with this level expected to increase with time. Species richness varies among the 32 federal entities in the country, from a low of 50 in Tlaxcala to a high of 492 in Oaxaca. Amphibian species richness by state-level can be envisioned as comprising three levels of low, medium, and high, with the lowest levels occurring in the Peninsula of Baja California, a group of seven states in north-central and central Mexico, and a group of three states in the Yucatan Peninsula, with the highest levels occupying the southern states of Guerrero, Puebla, Veracruz, Oaxaca, and Chiapas, and the medium level in the remaining states of the country. Reptile species richness also can be allocated to three categories, with the lowest level occupying Baja California Sur, a group of central states, and the states of the Yucatan Peninsula, and the highest level found in a cluster of the states of Veracruz, Guerrero, Oaxaca, and Chiapas. Knowledge of the Mexican herpetofauna will continue to grow with additional studies on systematics, conservation, and the construction of checklists at various levels.

Phylogeography and population genetics of a widespread cold-adapted ant, Prenolepis imparis

As species arise, evolve, and diverge, they are shaped by forces that unfold across short and long time scales and at both local and vast geographic scales. It is rare, however, to be able document this history across broad sweeps of time and space in a single species. Here, we report the results of a continental-scale phylogenomic analysis across the entire range of a widespread species. We analyzed sequences of 1,402 orthologous Ultraconserved Element (UCE) loci from 75 individuals to identify population genetic structure and historical demographic patterns across the continent-wide range of a cold-adapted ant, the winter ant, Prenolepis imparis. We recovered five well-supported, genetically isolated clades representing lineages that diverged from 8.2-2.2 million years ago. These include: 1) an early diverging lineage located in Florida, 2) a lineage that spans the southern United States, 3) populations that extend across the midwestern and northeastern United States, 4) populations from the western United States, and 5) populations in southwestern Arizona and Mexico. Population genetic analyses revealed little or no gene flow among these lineages, but patterns consistent with more recent gene flow among populations within lineages, and localized structure with migration in the western United States. High support for five major geographic lineages and lack of evidence of contemporary gene flow indicate in situ diversification across the species' range, producing relatively ancient lineages that persisted through subsequent climate change and glaciation during the Quaternary.

Phylogenomics of alligator lizards elucidate diversification patterns across the Mexican Transition Zone and support the recognition of a new genus

Genomic data continue to advance our understanding of species limits and biogeographic patterns. However, there is still no consensus regarding appropriate methods of phylogenomic analysis that make the best use of these heterogeneous data sets. In this study, we used thousands of ultraconserved element (UCE) loci from alligator lizards in the genus Gerrhonotus to compare and contrast species trees inferred using multiple contemporary methods and provide a time frame for biological diversification across the Mexican Transition Zone (MTZ). Concatenated maximum likelihood (ML) and Bayesian analyses provided highly congruent results, with differences limited to poorly supported nodes. Similar topologies were inferred from coalescent analyses in Bayesian Phylogenetics and Phylogeography and SVDquartets, albeit with lower support for some nodes. All divergence times fell within the Miocene, linking speciation to local Neogene vicariance and/or global cooling trends following the mid-Miocene Climatic Optimum. We detected a high level of genomic divergence for a morphologically distinct species restricted to the arid mountains of north-eastern Mexico, and erected a new genus to better reflect evolutionary history. In summary, our results further advocate leveraging the strengths and weaknesses of concatenation and coalescent methods, provide evidence for old divergences for alligator lizards, and indicate that the MTZ continues to harbour substantial unrecognized diversity.

A new extinct species of alligator lizard (Squamata: Elgaria) and an expanded perspective on the osteology and phylogeny of Gerrhonotinae

BACKGROUND

Alligator lizards (Gerrhonotinae) are a well-known group of extant North American lizard. Although many fossils were previously referred to Gerrhonotinae, most of those fossils are isolated and fragmentary cranial elements that could not be placed in a precise phylogenetic context, and only a handful of known fossils are articulated skulls. The fossil record has provided limited information on the biogeography and phylogeny of Gerrhonotinae.

RESULTS

We redescribe a nearly complete articulated fossil skull from the Pliocene sediments of the Anza-Borrego Desert in southern California, and refer the specimen to the alligator lizard genus Elgaria. The fossil is a representative of a newly described species, Elgaria peludoverde. We created a morphological matrix to assess the phylogeny of alligator lizards and facilitate identifications of fossil gerrhonotines. The matrix contains a considerably expanded taxonomic sample relative to previous morphological studies of gerrhonotines, and we sampled two specimens for many species to partially account for intraspecific variation. Specimen-based phylogenetic analyses of our dataset using Bayesian inference and parsimony inferred that Elgaria peludoverde is part of crown Elgaria. The new species is potentially related to the extant species Elgaria kingii and Elgaria paucicarinata, but that relationship was not strongly supported, probably because of extensive variation among Elgaria. We explored several alternative biogeographic scenarios implied by the geographic and temporal occurrence of the new species and its potential phylogenetic placements.

CONCLUSIONS

Elgaria peludoverde is the first described extinct species of Elgaria and provides new information on the biogeographic history and diversification of Elgaria. Our research expands the understanding of phylogenetic relationships and biogeography of alligator lizards and strengthens the foundation of future investigations. The osteological data and phylogenetic matrix that we provided will be critical for future efforts to place fossil gerrhonotines. Despite limited intraspecific sampled sizes, we encountered substantial variation among gerrhonotines, demonstrating the value of exploring patterns of variation for morphological phylogenetics and for the phylogenetic placement of fossils. Future osteological investigations on the species we examined and on species we did not examine will continue to augment our knowledge of patterns of variation in alligator lizards and aid in phylogenetics and fossil placement.

Variation in the skulls of Elgaria and Gerrhonotus (Anguidae, Gerrhonotinae) and implications for phylogenetics and fossil identification

BACKGROUND

There are limited data on intra- and interspecific osteological variation for many squamate clades. Those data are relevant for phylogenetic analyses that use osteological characters and for apomorphic identifications of fossils. We investigate whether morphological features in the skulls of extant gerrhonotine lizards can be used to distinguish taxa at the species- and genus-level and assess whether newly discovered intra- and interspecific osteological variation alters the utility of previously reported apomorphic features. We examined skulls of species belonging to the gerrhonotine genera Elgaria and Gerrhonotus. These genera contain 17 extant species, but the cranial osteology of only a few species was previously examined. As a result, intra- and interspecific osteological variation of these gerrhonotines is poorly understood.

METHODS

We employed high-resolution x-ray computed tomography (CT) to scan 25 alcohol-preserved specimens. We provide data on the skulls of all eight species of Elgaria, four for the first time, and five species of Gerrhonotus, three for the first time. We examined 3-D reconstructed skulls of the scanned specimens as well as dry, traditionally prepared skeletons (when they were available).

RESULTS

We found that the purported diagnostic utility of many previously described morphological features is impacted because of substantial morphological variation between and within species. We present an assessment of osteological differences that may be useful to differentiate species of Elgaria and Gerrhonotus, many of which are present on isolated cranial elements commonly recovered as fossils, including the premaxilla, maxilla, parietal, pterygoid, prootic, dentary, and surangular. We demonstrate the importance of documenting patterns of osteological variation using large sample sizes, and the utility of examining disarticulated cranial elements of the squamate skull to identify diagnostic morphology. This study adds to a growing body of literature suggesting that extensive documentation of morphological variation is needed to further our understanding of the phylogenetic and diagnostic utility of morphological features across vertebrate clades. Efforts in that direction likely will benefit from examination of disarticulated skeletal elements.

Mitochondrial genome of Diploderma micangshanense and its implications for phylogeny of the genus Diploderma

The lizard Diploderma micangshanense, which belongs to the family Agamidae is endemic to China. Here, we determined the complete mitogenome of D. micangshanense using an Illumina Hiseq X Ten sequencer. This mitogenome’s structure is a typical circular molecule of 16,467 bp in length, consisting of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a control region. The overall base composition of D. micangshanensis is 34.1% A, 23.64% T, 13.62% C, and 28.64% G with a slight AT bias of 57.74%. Most mitochondrial genes except ND6 and seven tRNAs were encoded on the heavy strand. Notably, the trnP gene was encoded on the heavy strand instead of its typical light strand position, providing an example of gene inversion in vertebrate mitogenomes. Phylogenetic analysis indicated that D. micangshanensis had a close relationship with D. zhaoermii.

A new species of Cynomops (Chiroptera: Molossidae) from the northwestern slope of the Andes

The systematics and taxonomy of the broadly distributed bats of the genus Cynomops has changed considerably in the last few years. Among the major changes, Cynomops abrasus was split into two species of large-bodied forms (Cynomops mastivus and C. abrasus) distributed east of the Andes. However, large Colombian specimens identified as C. abrasus from the western side of the Andes had yet to be included in any revisionary work. Phylogenetic analysis performed in this study, using mtDNA sequences (Cytochrome-b), revealed that these Colombian individuals are more closely related to Cynomops greenhalli. Morphological and molecular data allowed us to recognize populations from western Colombia, western Ecuador and northwestern Peru, as members of a new species of Cynomops. Characters that allow for its differentiation from C. greenhalli include a larger forearm, paler but more uniform ventral pelage, more globular braincase, and well-developed zygomatic processes of the maxilla (almost reaching the postorbital constriction). This study serves as another example of the importance of including multiple lines of evidence in the recognition of a new species. Given its rarity and the advanced transformation of its habitat, this new species is particularly important from a conservation perspective.

DNA barcoding reveals cryptic diversity in the underestimated genus Triplophysa (Cypriniformes: Cobitidae, Nemacheilinae) from the northeastern Qinghai-Tibet Plateau

Background

The northeastern part of the Qinghai-Tibet Plateau (QTP) presents a high number of plateau loach species. As one of the three major groups of fishes distributed on the QTP, plateau loach has high ecological value. However, the taxonomy and systematics of these fish are still controversial, and a large number of new species have been reported. The reason for this phenomenon is that the degree of morphological variation is low, the phylogenetic information provided by morphological and anatomical features used for species identification is relatively poor, and many cryptic species are observed. Based on the high-density sampling points from the biodiversity hotspots surveyed, this study aims to evaluate the biodiversity of plateau loach in the northeastern part of the QTP and reveal the hidden diversity by comparing morphological species with molecular operational taxonomic units (MOTUs).

Results

After careful identification and comparison of the morphology and DNA barcoding of 1630 specimens, 22 species were identified, with 20 considered valid local species and two identified as new species that had not been previously described. Based on the combination of morphological and molecular methods, a total of 24 native species were found, two of which were cryptic species: Triplophysa robusta sp1 and Triplophysa minxianensis sp1. Fourteen of the 24 species form clusters of barcodes that allow them to be reliably identified. The remaining cases involved 10 closely related species, including rapidly differentiated species and species that seemed to have experienced incomplete lineage sorting or showed introgressions.

Conclusions

The results highlight the need to combine traditional taxonomies with molecular methods to correctly identify species, especially closely related species, such as the plateau loach. This study provides a basis for protecting the biodiversity of plateau loach.

Developmental morphology and evolution of extraembryonic membranes of lizards and snakes (Reptilia, Squamata)

Amniote embryos are supported and nourished by a suite of tissues, the extraembryonic membranes, that provide vascular connections to the egg contents. Oviparous reptiles share a basic pattern of development inherited from a common ancestor; a vascular chorioallantoic membrane, functioning as a respiratory organ, contacts the eggshell and a vascular yolk sac membrane conveys nutrients to the embryo. Squamates (lizards, snakes) have evolved a novel variation in morphogenesis of the yolk sac that results in a unique structure, the yolk cleft/isolated yolk mass complex. This structure is a source of phylogenetic variation in architecture of the extraembryonic membranes among oviparous squamates. The yolk cleft/isolated yolk mass complex is retained in viviparous species and influences placental architecture. The aim of this paper is to review extraembryonic membrane development and morphology in oviparous and related viviparous squamates to explore patterns of variation. The survey includes all oviparous species for which data are available (11 species; 4 families). Comparisons with viviparous species encompass six independent origins of viviparity. The comparisons reveal that both phylogeny and reproductive mode influence variation in extraembryonic membrane development and that phylogenetic variation influences placental evolution. Models of the evolution of squamate placentation have relied primarily on comparisons between independently derived viviparous species. The inclusion of oviparous species in comparative analyses largely supports these models, yet exposes convergent patterns of evolution that become apparent when phylogenetic variation is recognized.

Phylogeography of the iconic Australian red-tailed black-cockatoo (Calyptorhynchus banksii) and implications for its conservation

Advances in sequencing technologies have revolutionized wildlife conservation genetics. Analysis of genomic data sets can provide high-resolution estimates of genetic structure, genetic diversity, gene flow, and evolutionary history. These data can be used to characterize conservation units and to effectively manage the genetic health of species in a broad evolutionary context. Here we utilize thousands of genome-wide single-nucleotide polymorphisms (SNPs) and mitochondrial DNA to provide the first genetic assessment of the Australian red-tailed black-cockatoo (Calyptorhynchus banksii), a widespread bird species comprising populations of varying conservation concern. We identified five evolutionarily significant units, which are estimated to have diverged during the Pleistocene. These units are only partially congruent with the existing morphology-based subspecies taxonomy. Genetic clusters inferred from mitochondrial DNA differed from those based on SNPs and were less resolved. Our study has a range of conservation and taxonomic implications for this species. In particular, we provide advice on the potential genetic rescue of the Endangered and restricted-range subspecies C. b. graptogyne, and propose that the western C. b. samueli population is diagnosable as a separate subspecies. The results of our study highlight the utility of considering the phylogeographic relationships inferred from genome-wide SNPs when characterizing conservation units and management priorities, which is particularly relevant as genomic data sets become increasingly accessible.

Introgression obscures lineage boundaries and phylogeographic history in the western banded gecko, Coleonyx variegatus (Squamata: Eublepharidae)

The geomorphological formation of the Baja California peninsula and the Gulf of California is a principal driver of diversification for the reptiles of North America’s warm deserts. The western banded gecko, Coleonyx variegatus, is distributed throughout the Mojave, Sonoran and Peninsular deserts. In this study we use multilocus sequence data to address deep phylogeographic structure within C. variegatus. Analyses of mtDNA data recover six divergent clades throughout the range of C. variegatus. Topology of the mtDNA gene tree suggests separate origins of peninsular populations with an older lineage in the south and a younger one in the north. In contrast, analyses of multilocus nuclear data provide support for four lineages, corresponding to the subspecies C. v. abbotti, C. v. peninsularis, C. v. sonoriensis and C. v. variegatus. Phylogenetic analyses of the nuclear data recover C. v. abbotti and C. v. peninsularis as a clade, indicating a single origin of the peninsular populations. Discordance between the nuclear and mtDNA data is largely the result of repeated episodes of mtDNA introgression that have obscured both lineage boundaries and biogeographic history. Dating analyses of the combined nuclear and mtDNA data suggest that the peninsular clade diverged from the continental group in the Late Miocene.

Out of the Hengduan Mountains: Molecular phylogeny and historical biogeography of the Asian water snake genus Trimerodytes (Squamata: Colubridae)

The Asian water snake genus Trimerodytes (formerly Sinonatrix) is endemic to East and Southeast Asia. Although several species have been included in various phylogenetic studies previously, the evolution and relationships among members of this genus as a whole remain unexplored. In this study, we report the sequencing two protein-coding mitochondrial gene fragments (MTCYB and ND2) and three nuclear genes (c-mos, NT3, and Rag1), reconstruct interspecific phylogeny, and explore biogeography for the genus Trimerodytes. Both Bayesian inference and maximum likelihood analyses consistently recover the monophyly of Trimerodytes with strong support, with T. yapingi the sister-group to the remaining species. The divergence date and ancestral area estimation suggest that Trimerodytes likely originated in Hengduan Mountains (eastern Tibetan Plateau) in western China at 23.93 Ma (95% HPD: 17.09-31.30), and intraspecific divergence began at about 4.23 Ma (95% HPD: 2.74-6.10). Analyses support the validity of T. yunnanensis.

Multiple convergent events created a nominal widespread species: Triplophysa stoliczkae (Steindachner, 1866) (Cobitoidea: Nemacheilidae)

BACKGROUND

Triplophysa stoliczkae is the most widespread species in the genus Triplophysa and may have originated from morphological convergence. To understand the evolutionary history of T. stoliczkae, we employed a multilocus approach to investigate the phylogenetics and the morphological evolution of T. stoliczkae on the Qinghai-Tibetan Plateau.

RESULTS

All phylogenetic analyses (two mitochondrial and five nuclear loci), a genealogical sorting index and species tree inferences suggested that T. stoliczkae consists of distinct lineages that were not closest relatives. The time estimation indicated that the divergence events between "T. stoliczkae" and other Triplophysa species occurred from approximately 0.10 to 4.51 Ma. The ancestral state analyses supported the independent evolution of T. stoliczkae morphology in distinct lineages. The morphometric analysis and convergence estimates demonstrated significant phenotypic convergence among "T. stoliczkae" lineages.

CONCLUSIONS

Triplophysa stoliczkae includes 4 different lineages with similar morphologies. The increasingly harsh environments that have occurred since the Pliocene have driven the occurrences of scrape-feeding fish in the genus Triplophysa. Morphological adaptations associated with scrape-feeding behavior resulted in convergences and the artificial lumping of four different species in the nominal taxon T. stoliczkae. A taxonomic revision for T. stoliczkae is needed.

Introgressive hybridization in a Spiny-Tailed Iguana, Ctenosaura pectinata, and its implications for taxonomy and conservation

Introgression, the transmission of genetic material of one taxon into another through hybridization, can have various evolutionary outcomes. Previous studies have detected signs of introgression between western populations of the Mexican endemic and threatened spiny-tailed iguana, Ctenosaura pectinata. However, the extent of this phenomenon along the geographic distribution of the species is unknown. Here, we use multilocus data together with detailed geographic sampling to (1) define genotypic clusters within C. pectinata; (2) evaluate geographic concordance between maternally and biparentally inherited markers; (3) examine levels of introgression between genotypic clusters, and (4) suggest taxonomic modifications in light of this information. Applying clustering methods to genotypes of 341 individuals from 49 localities of C. pectinata and the closely related C. acanthura, we inferred the existence of five genotypic clusters. Contact zones between genotypic clusters with signatures of interbreeding were detected, showing different levels of geographic discordance with mtDNA lineages. In northern localities, mtDNA and microsatellites exhibit concordant distributions, supporting the resurrection of C. brachylopha. Similar concordance is observed along the distribution of C. acanthura, confirming its unique taxonomic identity. Genetic and geographic concordance is also observed for populations within southwestern Mexico, where the recognition of a new species awaits in depth taxonomic revision. In contrast, in western localities a striking pattern of discordance was detected where up to six mtDNA lineages co-occur with only two genotypic clusters. Given that the type specimen originated from this area, we suggest that individuals from western Mexico keep the name C. pectinata. Our results have profound implications for conservation, management, and forensics of Mexican iguanas.

Phylogeography and niche modelling: reciprocal enlightenment

Phylogeography examines the spatial genetic structure of species. Environmental niche modelling (or ecological niche modelling; ENM) examines the environmental limits of a species’ ecological niche. These two fields have great potential to be used together. ENM can shed light on how phylogeographical patterns develop and help identify possible drivers of spatial structure that need to be further investigated. Specifically, ENM can be used to test for niche differentiation among clades, identify factors limiting individual clades and identify barriers and contact zones. It can also be used to test hypotheses regarding the effects of historical and future climate change on spatial genetic patterns by projecting niches using palaeoclimate or future climate data. Conversely, phylogeographical information can populate ENM with within-species genetic diversity. Where adaptive variation exists among clades within a species, modelling their niches separately can improve predictions of historical distribution patterns and future responses to climate change. Awareness of patterns of genetic diversity in niche modelling can also alert conservationists to the potential loss of genetically diverse areas in a species’ range. Here, we provide a simplistic overview of both fields, and focus on their potential for integration, encouraging researchers on both sides to take advantage of the opportunities available.

Molecular systematics of the Triplophysa robusta (Cobitoidea) complex: Extensive gene flow in a depauperate lineage

Gene flow between populations assumed to be isolated frequently leads to incorrect inferences of evolutionary history. Understanding gene flow and its causes has long been a key topic in evolutionary biology. In this study, we explored the evolutionary history of the Triplophysa robusta complex, using a combination of multilocus analyses and coalescent simulation. Our multilocus approach detected conspicuous mitonuclear discordances in the T. robusta complex. Mitochondrial results showed reticular clades, whereas the nuclear results corresponded with the morphological data. Coalescent simulation indicated that gene flow was the source of these discordances. Molecular clock analysis combined with geological processes suggest that intense geological upheavals have shaped a complicated evolutionary history for the T. robusta complex since the late Miocene, causing extensive gene flow which has distorted the molecular systematics of the T. robusta complex. We suggest that frequent gene flow may restrict speciation in the T. robusta complex, leading to such a depauperate lineage. Based on this comprehensive understanding, we provide our proposals for taxonomic revision of the T. robusta complex.

The skull of the gerrhonotine lizard Elgaria panamintina (Squamata: Anguidae)

We provide the first description of the skull, osteoderms, and hyoid apparatus of the poorly known alligator lizard Elgaria panamintina, and compare the cranial osteology of that species to the widespread and well-studied taxon Elgaria multicarinata. Patterns of morphological variation resulting from ontogenetic transformations and pathology are discussed. We employed x-ray computed tomography (CT) scans to examine two adult specimens of Elgaria panamintina and two adult specimens of Elgaria multicarinata, in addition to examining multiple traditionally prepared skeletal specimens of the latter species. CT scans provide simultaneous study of both articulated and disarticulated elements, allowing us to describe and document the morphology of the skull with exceptional precision and detail. The description of the skull of Elgaria panamintina serves as a generalization for all Elgaria; here we provide the first complete description of the skull of this genus for future uses in morphological and phylogenetic studies of both extant species and fossils.

Supervised machine learning reveals introgressed loci in the genomes of Drosophila simulans and D. sechellia

Hybridization and gene flow between species appears to be common. Even though it is clear that hybridization is widespread across all surveyed taxonomic groups, the magnitude and consequences of introgression are still largely unknown. Thus it is crucial to develop the statistical machinery required to uncover which genomic regions have recently acquired haplotypes via introgression from a sister population. We developed a novel machine learning framework, called FILET (Finding Introgressed Loci via Extra-Trees) capable of revealing genomic introgression with far greater power than competing methods. FILET works by combining information from a number of population genetic summary statistics, including several new statistics that we introduce, that capture patterns of variation across two populations. We show that FILET is able to identify loci that have experienced gene flow between related species with high accuracy, and in most situations can correctly infer which population was the donor and which was the recipient. Here we describe a data set of outbred diploid Drosophila sechellia genomes, and combine them with data from D. simulans to examine recent introgression between these species using FILET. Although we find that these populations may have split more recently than previously appreciated, FILET confirms that there has indeed been appreciable recent introgression (some of which might have been adaptive) between these species, and reveals that this gene flow is primarily in the direction of D. simulans to D. sechellia.

Understanding the extent to which species or diverged populations hybridize in nature is crucially important if we are to understand the speciation process. Accordingly numerous research groups have developed methodology for finding the genetic evidence of such introgression. In this report we develop a supervised machine learning approach for uncovering loci which have introgressed across species boundaries. We show that our method, FILET, has greater accuracy and power than competing methods in discovering introgression, and in addition can detect the directionality associated with the gene flow between species. Using whole genome sequences from Drosophila simulans and Drosophila sechellia we show that FILET discovers quite extensive introgression between these species that has occurred mostly from D. simulans to D. sechellia. Our work highlights the complex process of speciation even within a well-studied system and points to the growing importance of supervised machine learning in population genetics.

Conflicting Evolutionary Histories of the Mitochondrial and Nuclear Genomes in New World Myotis Bats

The rapid diversification of Myotis bats into more than 100 species is one of the most extensive mammalian radiations available for study. Efforts to understand relationships within Myotis have primarily utilized mitochondrial markers and trees inferred from nuclear markers lacked resolution. Our current understanding of relationships within Myotis is therefore biased towards a set of phylogenetic markers that may not reflect the history of the nuclear genome. To resolve this, we sequenced the full mitochondrial genomes of 37 representative Myotis, primarily from the New World, in conjunction with targeted sequencing of 3648 ultraconserved elements (UCEs). We inferred the phylogeny and explored the effects of concatenation and summary phylogenetic methods, as well as combinations of markers based on informativeness or levels of missing data, on our results. Of the 294 phylogenies generated from the nuclear UCE data, all are significantly different from phylogenies inferred using mitochondrial genomes. Even within the nuclear data, quartet frequencies indicate that around half of all UCE loci conflict with the estimated species tree. Several factors can drive such conflict, including incomplete lineage sorting, introgressive hybridization, or even phylogenetic error. Despite the degree of discordance between nuclear UCE loci and the mitochondrial genome and among UCE loci themselves, the most common nuclear topology is recovered in one quarter of all analyses with strong nodal support. Based on these results, we re-examine the evolutionary history of Myotis to better understand the phenomena driving their unique nuclear, mitochondrial, and biogeographic histories.