Plastome variation and phylogeny of Taxillus (Loranthaceae).

Several molecular phylogenetic studies of the mistletoe family Loranthaceae have been published such that now the general pattern of relationships among the genera and their biogeographic histories are understood. Less is known about species relationships in the larger (> 10 species) genera. This study examines the taxonomically difficult genus Taxillus composed of 35-40 Asian species. The goal was to explore the genetic diversity present in Taxillus plastomes, locate genetically variable hotspots, and test these for their utility as potential DNA barcodes. Using genome skimming, complete plastomes, as well as nuclear and mitochondrial rDNA sequences, were newly generated for eight species. The plastome sequences were used in conjunction with seven publicly available Taxillus sequences and three sequences of Scurrula, a close generic relative. The Taxillus plastomes ranged from 121 to 123 kbp and encoded 90-93 plastid genes. In addition to all of the NADH dehydrogenase complex genes, four ribosomal genes, infA and four intron-containing tRNA genes were lost or pseudogenized in all of the Taxillus and Scurrula plastomes. The topologies of the plastome, mitochondrial rDNA and nuclear rDNA trees were generally congruent, though with discordance at the position of T. chinensis. Several variable regions in the plastomes were identified that have sufficient numbers of parsimony informative sites as to recover the major clades seen in the complete plastome tree. Instead of generating complete plastome sequences, our study showed that accD alone or the concatenation of accD and rbcL can be used in future studies to facilitate identification of Taxillus samples and to generate a molecular phylogeny with robust sampling within the genus.

A five-gene molecular phylogeny reveals Parapanteles Ashmead (Hymenoptera: Braconidae) to be polyphyletic as currently composed.

Parapanteles Ashmead (Braconidae: Microgastrinae) is a medium-sized genus of microgastrine wasps that was erected over a century ago and lacks a unique synapomorphic character, and its monophyly has not been tested by any means. Parapanteles usually are parasitoids of large, unconcealed caterpillars (macrolepidoptera) and have been reared from an unusually large diversity of hosts for a relatively small microgastrine genus. We used Cytochrome Oxidase I sequences ("DNA barcodes") available for Parapanteles and other microgastrines to sample the generic diversity of described and undescribed species currently placed in Parapanteles, and then sequenced four additional genes for this subsample (wingless, elongation factor 1-alpha, ribosomal subunit 28s, and NADH dehydrogenase subunit 1). We constructed individual gene trees and concatenated Bayesian and maximum-likelihood phylogenies for this 5-gene subsample. In these phylogenies, most Parapanteles species formed a monophyletic clade within another genus, Dolichogenidea, while the remaining Parapanteles species were recovered polyphyletically within several other genera. The latter likely represent misidentified members of other morphologically similar genera. Species in the monophyletic clade containing most Parapanteles parasitized caterpillars from only five families - Erebidae (Arctiinae), Geometridae, Saturniidae, Notodontidae, and Crambidae. We do not make any formal taxonomic decisions here because we were not able to include representatives of type species for Parapanteles or other relevant genera, and because we feel such decisions should be reserved until a comprehensive morphological analysis of the boundaries of these genera is accomplished.

Rivers of Indochina as potential drivers of lineage diversification in the spotted flying lizard (Draco maculatus) species complex.

Southeast Asia hosts a rich concentration of biodiversity within multiple biodiversity hotspots. Indochina, a region with remarkably high levels of in situ diversification, possesses five major rivers (Ayeyarwady, Chiang Mai, Mekong, Red, and Salween), several of which coincide with phylogenetic breaks of terrestrial taxa. Draco maculatus possesses a range that stretches across Indochina, which widespread geographic distribution along with potential discrete variation within subspecies alludes to the possibility of this taxon constituting multiple divergent lineages. Using sequence data from three mitochondrial (12S, 16S, and ND2) and three nuclear (BDNF, CMOS, and PNN) genes, we provide the first estimated phylogeny of this hypothesized species complex and examine its phylogeographic architecture with maximum likelihood and Bayes factor delimitation (BFD) approaches. Our results support multiple divergent lineages with phylogenetic breaks coincident with rivers, indicating that river barriers may be contributing to the elevated levels of in situ diversification of Indochina.

MpsAB is important for Staphylococcus aureus virulence and growth at atmospheric CO2 levels.

The mechanisms behind carbon dioxide (CO2) dependency in non-autotrophic bacterial isolates are unclear. Here we show that the Staphylococcus aureus mpsAB operon, known to play a role in membrane potential generation, is crucial for growth at atmospheric CO2 levels. The genes mpsAB can complement an Escherichia coli carbonic anhydrase (CA) mutant, and CA from E. coli can complement the S. aureus delta-mpsABC mutant. In comparison with the wild type, S. aureus mps mutants produce less hemolytic toxin and are less virulent in animal models of infection. Homologs of mpsA and mpsB are widespread among bacteria and are often found adjacent to each other on the genome. We propose that MpsAB represents a dissolved inorganic carbon transporter, or bicarbonate concentrating system, possibly acting as a sodium bicarbonate cotransporter.

Alignment-free similarity analysis for protein sequences based on fuzzy integral.

Sequence comparison is an essential part of modern molecular biology research. In this study, we estimated the parameters of Markov chain by considering the frequencies of occurrence of the all possible amino acid pairs from each alignment-free protein sequence. These estimated Markov chain parameters were used to calculate similarity between two protein sequences based on a fuzzy integral algorithm. For validation, our result was compared with both alignment-based (ClustalW) and alignment-free methods on six benchmark datasets. The results indicate that our developed algorithm has a better clustering performance for protein sequence comparison.

LMAP_S: Lightweight Multigene Alignment and Phylogeny eStimation.

BACKGROUND: Recent advances in genome sequencing technologies and the cost drop in high-throughput sequencing continue to give rise to a deluge of data available for downstream analyses. Among others, evolutionary biologists often make use of genomic data to uncover phenotypic diversity and adaptive evolution in protein-coding genes. Therefore, multiple sequence alignments (MSA) and phylogenetic trees (PT) need to be estimated with optimal results. However, the preparation of an initial dataset of multiple sequence file(s) (MSF) and the steps involved can be challenging when considering extensive amount of data. Thus, it becomes necessary the development of a tool that removes the potential source of error and automates the time-consuming steps of a typical workflow with high-throughput and optimal MSA and PT estimations. RESULTS: We introduce LMAP_S (Lightweight Multigene Alignment and Phylogeny eStimation), a user-friendly command-line and interactive package, designed to handle an improved alignment and phylogeny estimation workflow: MSF preparation, MSA estimation, outlier detection, refinement, consensus, phylogeny estimation, comparison and editing, among which file and directory organization, execution, manipulation of information are automated, with minimal manual user intervention. LMAP_S was developed for the workstation multi-core environment and provides a unique advantage for processing multiple datasets. Our software, proved to be efficient throughout the workflow, including, the (unlimited) handling of more than 20 datasets. CONCLUSIONS: We have developed a simple and versatile LMAP_S package enabling researchers to effectively estimate multiple datasets MSAs and PTs in a high-throughput fashion. LMAP_S integrates more than 25 software providing overall more than 65 algorithm choices distributed in five stages. At minimum, one FASTA file is required within a single input directory. To our knowledge, no other software combines MSA and phylogeny estimation with as many alternatives and provides means to find optimal MSAs and phylogenies. Moreover, we used a case study comparing methodologies that highlighted the usefulness of our software. LMAP_S has been developed as an open-source package, allowing its integration into more complex open-source bioinformatics pipelines. LMAP_S package is released under GPLv3 license and is freely available at https://lmap-s.sourceforge.io/.

Uncovered variability in olive moth (Prays oleae) questions species monophyly.

The olive moth -Prays oleae Bern.- remains a significant pest of olive trees showing situation dependent changes in population densities and in severity of damages. The genetic variability of olive moth was assessed on three main olive orchards regions in Portugal by three different markers (COI, nad5 and RpS5), suggesting high species diversity albeit with no obvious relation with a regional pattern nor to an identified ecological niche. Selected COI sequences obtained in this study were combined with those available in the databases for Prays genus to generate a global dataset. The reconstruction of the Prays phylogeny based on this marker revealed the need to revise Prays oleae to confirm its status of single species: COI data suggests the co-existence of two sympatric evolutionary lineages of morphologically cryptic olive moth. We show, however, that the distinct mitochondrial subdivision observed in the partial COI gene fragment is not corroborated by the other DNA sequences. There is the need of understanding this paradigm and the extent of Prays variability, as the disclosure of lineage-specific differences in biological traits between the identified lineages is fundamental for the development of appropriate pest management practices.

A genomic perspective of the pink-headed duck Rhodonessa caryophyllacea suggests a long history of low effective population size.

The first molecular phylogenetic hypothesis for the possibly extinct pink-headed duck Rhodonessa caryophyllacea unambiguously shows that it belongs to the pochard radiation that also includes the genera Aythya and Netta. It is the sister to all modern-day pochards and belongs to a lineage that branched off from the others more than 2.8 million years ago. Rhodonessa caryophyllacea is believed to never have been common in modern time and we show this has probably been the situation for as long as 100,000 years. Our results suggest that their effective population size varied between 15,000 and 25,000 individuals during the last 150,000 years of the Pleistocene. The reasons behind this are largely unknown as very little is known about the life-history and biology of this species. Presumably it is due to factors related to feeding or to breeding, but we may never know this for sure.

Phylogenetic evaluation of Amynthas earthworms from South China reveals the initial ancestral state of spermathecae.

Our knowledge of the phylogeny of the earthworm genus Amynthas under the family Megascolecidae, which is comprised of a huge number of species, is very limited compared to the better-known and much smaller family Lumbricidae. In order to investigate the phylogenetic relationships among the species within the genus Amynthas, which is the largest genus of the Megascolecidae family, nuclear and mitochondrial DNA sequences of 77 species, including 76 in-group Amynthas species collected from South China and 1 out-group species, were analyzed. A 5402bp segments composed of whole nuclear 18S rDNA and the mitochondrial genes COI, COII, ND1, 12S, and 16S was assembled from 77 species. Maximum Likelihood and Bayesian analyses of the concatenated sequences were performed. The results revealed evolution of two geographically independent lineages, both showing the ancestral state of two pairs of spermatheca (Sp.p 7/8/9). We found the species groups described by Sims and Easton (1972) to be non-monophyletic, and the origin of the parthenogenetic species group to likely be a quadthecal ancestor. These results provide modest evidence in support of an Indochinese peninsula origin of the Chinese Amynthas species and divergence of the genus once it had spread to mainland China. The findings of this study are consistent with a divergence scenario that resulted in at least one branch spreading to the Southeast of China and another branch spreading to the Southwest of China, but further research is required to confirm this interpretation of the Amynthas phylogeny.

Never judge an iguana by its spines: Systematics of the Yucatan spiny tailed iguana, Ctenosaura defensor (Cope, 1866).

Spiny tailed iguanas are highly diverse clade of lizards in Mesoamerica, ranging from northern Mexico through Panama. Utilizing 2 regions of mitochondrial DNA (1948bp) and 4 nuclear loci (2232bp) we explored the relationships between these species and the phylogeographic history of the major clades. We discovered that the lineage endemic to the Yucatan Peninsula renders the genus Ctenosaura paraphyletic. To resolve this non-monophyly, we resurrect the taxon Cachryx Cope, 1866, and provide a new diagnosis for the genus. We also find that small body-size and highly spinose tails in the species previously referred to the subgenus Enyaliosaurus, have evolved independently 3 times. Cachryx were recovered as sister to the lineage of iguanines endemic to the Galapagos Islands, and we discuss biogeographic scenarios to explain this relationship as well as those among the primary clades of Ctenosaura in Mesoamerica.

Striking pseudogenization in avian phylogenetics: Numts are large and common in falcons.

Nuclear copies of mitochondrial genes (numts) are a well-known feature of eukaryotic genomes and a concern in systematics, as they can mislead phylogenetic inferences when inadvertently used. Studies on avian numts initially based on the chicken genome suggest that numts may be uncommon and relatively short among birds. Here we ask how common numts are in falcons, based on recently sequenced genomes of the Saker falcon (Falco cherrug) and Peregrine falcon (F. peregrinus). We identified numts by BLASTN searches and then extracted CYTB, ND2 and COI sequences from them, which were then used for phylogeny inference along with several sequences from other species in Falconiformes. Our results indicate that avian numts may be much more frequent and longer than previously thought. Phylogenetic inferences revealed multiple independent nuclear insertions throughout the history of the Falconiformes, including cases of sequences available in public databases and wrongly identified as authentic mtDNA. New sequencing technologies and ongoing efforts for whole genome sequencing will provide exciting opportunities for avian numt research in the near future.

Molecular phylogeny of Miltogramminae (Diptera: Sarcophagidae): Implications for classification, systematics and evolution of larval feeding strategies.

Miltogramminae is one of the phylogenetically most poorly studied taxa of the species-rich family Sarcophagidae (Diptera). Most species are kleptoparasites in nests of solitary aculeate wasps and bees, although parasitoids and saprophagous species are also known, and the ancestral miltogrammine life habit remains unsettled. Here, we present for the first time a comprehensive phylogenetic tree consisting of 58 representatives of Miltogramminae, reconstructed using sequence data from three mitochondrial (COI, cytB, ND4) and one nuclear (Ef-1α) genes. Our phylogenetic hypothesis suggests that: (1) Miltogramminae are sister to Paramacronychiinae, (2) Miltogramminae can be divided into the "lower miltogrammines" containing two clades of mainly saprophages and a clade of "higher miltogrammines" with mainly kleptoparasitic species, (3) only three genera turn out to be non-monophyletic: Miltogramma, Senotainia and Pterella and (4) the genus Sarcotachina, which traditionally has been considered as belonging to the Paramacronychiinae, is placed in one of the clades of "lower miltogrammines". Ancestral state reconstruction of larval feeding strategy and five larval characters reveals that the ancestor of Miltogramminae was likely a saprophage retaining plesiomorphic oral ridges and a cephaloskeleton with sclerotized dorsal bridge. Synapomorphies like large pseudocephalic sensory organs and well-developed cuticular sculpture suggest that the ancestral first instar larva actively searched for a buried food supply.

Comparison of different strategies for using fossil calibrations to generate the time prior in Bayesian molecular clock dating.

Fossil calibrations are the utmost source of information for resolving the distances between molecular sequences into estimates of absolute times and absolute rates in molecular clock dating analysis. The quality of calibrations is thus expected to have a major impact on divergence time estimates even if a huge amount of molecular data is available. In Bayesian molecular clock dating, fossil calibration information is incorporated in the analysis through the prior on divergence times (the time prior). Here, we evaluate three strategies for converting fossil calibrations (in the form of minimum- and maximum-age bounds) into the prior on times, which differ according to whether they borrow information from the maximum age of ancestral nodes and minimum age of descendent nodes to form constraints for any given node on the phylogeny. We study a simple example that is analytically tractable, and analyze two real datasets (one of 10 primate species and another of 48 seed plant species) using three Bayesian dating programs: MCMCTree, MrBayes and BEAST2. We examine how different calibration strategies, the birth-death process, and automatic truncation (to enforce the constraint that ancestral nodes are older than descendent nodes) interact to determine the time prior. In general, truncation has a great impact on calibrations so that the effective priors on the calibration node ages after the truncation can be very different from the user-specified calibration densities. The different strategies for generating the effective prior also had considerable impact, leading to very different marginal effective priors. Arbitrary parameters used to implement minimum-bound calibrations were found to have a strong impact upon the prior and posterior of the divergence times. Our results highlight the importance of inspecting the joint time prior used by the dating program before any Bayesian dating analysis.

Trans-biome diversity in Australian grass-specialist lizards (Diplodactylidae: Strophurus).

Comparisons of biodiversity patterns within lineages that occur across major climate gradients and biomes, can provide insights into the relative roles that lineage history, landscape and climatic variation, and environmental change have played in shaping regional biotas. In Australia, while there has been extensive research into the origins and patterns of diversity in the Australian Arid Zone (AAZ), how diversity is distributed across this biome and the Australian Monsoonal Tropics (AMT) to the north, has been less studied. We compared the timing and patterns of diversification across this broad aridity gradient in a clade of lizards (Strophurus: phasmid geckos) that only occur in association with a unique Australian radiation of sclerophyllous grasses (Triodia: spinifex). Our results indicate that overall genetic diversity is much higher, older and more finely geographically structured within the AMT, including distantly related clades endemic to the sandstone escarpments of the Kimberley and Arnhem Plateau. Niche modelling analyses also suggest that the distribution of taxa in the AMT is more strongly correlated with variation in topographic relief than in the AAZ. The two broad patterns that we recovered - (i) lineage endemism increases as latitude decreases, and (ii) endemism is tightly correlated to rocky regions - parallel and corroborate other recent studies of habitat generalists and specialised saxicoline lineages occurring across these same regions. Early Miocene diversification estimates also suggest that, soon after Triodia grasses colonised Australia and began to diversify in the Miocene, phasmid geckos with Gondwanan ancestry shifted into these grasses, and have subsequently remained closely associated with this unique vegetation type.

The effects of Pleistocene climate change on biotic differentiation in a montane songbird clade from Wallacea.

The role of Pleistocene Ice Age in tropical diversification is poorly understood, especially in archipelagos, in which glaciation-induced sea level fluctuations may lead to complicated changes in land distribution. To assess how Pleistocene land bridges may have facilitated gene flow in tropical archipelagos, we investigated patterns of diversification in the rarely-collected rusty-bellied fantail Rhipidura teysmanni (Passeriformes: Rhipiduridae) complex from Wallacea using a combination of bioacoustic traits and whole-genome sequencing methods (dd-RADSeq). We report a biogeographic leapfrog pattern in the vocalizations of these birds, and uncover deep genomic divergence among island populations despite the presence of intermittent land connections between some. We demonstrate how rare instances of genetic introgression have affected the evolution of this species complex, and document the presence of double introgressive mitochondrial sweeps, highlighting the dangers of using only mitochondrial DNA in evolutionary research. By applying different tree inference approaches, we demonstrate how concatenation methods can give inaccurate results when investigating divergence in closely-related taxa. Our study highlights high levels of cryptic avian diversity in poorly-explored Wallacea, elucidates complex patterns of Pleistocene climate-mediated diversification in an elusive montane songbird, and suggests that Pleistocene land bridges may have accounted for limited connectivity among montane Wallacean biota.

Sequence variation in three mitochondrial genes among Gasterophilus intestinalis isolates from two distinct regions in China.

Horse gastrointestinal myiasis caused by larvae of Gasterophilus spp. has a worldwide distribution. However, little information is available on the genetic variation of Gasterophilus spp. at the molecular level. In the present study, sequence variation was examined in three mitochondrial (mt) genes, namely mt cytochrome c oxidase subunit 1 (cox1), NADH dehydrogenase subunits 5 (nad5) and large subunit ribosomal RNA (rrnL), among G. intestinalis isolates originating from Heilongjiang Provience and Xinjiang Uygur Autonomous Region, China. The complete cox1, nad5, and rrnL gene sequences were amplified by PCR separately from adult G. intestinalis individuals and the amplicons were sequenced from both directions. The size of the sequences of cox1, nad5, and rrnL genes was 1539 bp, 1717 bp, and 1321 bp, respectively. The A + T contents of the sequences were 62.7-63.9% (cox1), 70-70.7% (nad5), and 74.5-74.9% (rrnL). The intra-specific sequence variations within G. intestinalis were 0.1-3.1% for cox1, 0.4-3.4% for nad5, and 0.1-1.5% for rrnL. Phylogenetic analyses based on the cox1 sequences using the Bayesian inference indicated that all the G. intestinalis isolates grouped together with high statistical support. These findings demonstrated clearly the usefulness of mt cox1, nad5, and rrnL sequences for the molecular identification and for studying population genetics of G. intestinalis in horses.

The Chloroplast Genome of Utricularia reniformis Sheds Light on the Evolution of the ndh Gene Complex of Terrestrial Carnivorous Plants from the Lentibulariaceae Family.

Lentibulariaceae is the richest family of carnivorous plants spanning three genera including Pinguicula, Genlisea, and Utricularia. Utricularia is globally distributed, and, unlike Pinguicula and Genlisea, has both aquatic and terrestrial forms. In this study we present the analysis of the chloroplast (cp) genome of the terrestrial Utricularia reniformis. U. reniformis has a standard cp genome of 139,725bp, encoding a gene repertoire similar to essentially all photosynthetic organisms. However, an exclusive combination of losses and pseudogenization of the plastid NAD(P)H-dehydrogenase (ndh) gene complex were observed. Comparisons among aquatic and terrestrial forms of Pinguicula, Genlisea, and Utricularia indicate that, whereas the aquatic forms retained functional copies of the eleven ndh genes, these have been lost or truncated in terrestrial forms, suggesting that the ndh function may be dispensable in terrestrial Lentibulariaceae. Phylogenetic scenarios of the ndh gene loss and recovery among Pinguicula, Genlisea, and Utricularia to the ancestral Lentibulariaceae cladeare proposed. Interestingly, RNAseq analysis evidenced that U. reniformis cp genes are transcribed, including the truncated ndh genes, suggesting that these are not completely inactivated. In addition, potential novel RNA-editing sites were identified in at least six U. reniformis cp genes, while none were identified in the truncated ndh genes. Moreover, phylogenomic analyses support that Lentibulariaceae is monophyletic, belonging to the higher core Lamiales clade, corroborating the hypothesis that the first Utricularia lineage emerged in terrestrial habitats and then evolved to epiphytic and aquatic forms. Furthermore, several truncated cp genes were found interspersed with U. reniformis mitochondrial and nuclear genome scaffolds, indicating that as observed in other smaller plant genomes, such as Arabidopsis thaliana, and the related and carnivorous Genlisea nigrocaulis and G. hispidula, the endosymbiotic gene transfer may also shape the U. reniformis genome in a similar fashion. Overall the comparative analysis of the U. reniformis cp genome provides new insight into the ndh genes and cp genome evolution of carnivorous plants from Lentibulariaceae family.

Cryptic species diversity in sub-Antarctic islands: A case study of Lepidonotothen.

The marine fauna of the Southern Ocean is well known for an impressive adaptive radiation of fishes, the notothenioids. However, when compared to other marine areas, the frigid waters of the Southern Ocean also contain a seemingly large proportion of cryptic species. The documented instances of speciation in the absence of morphological change are largely observed in invertebrate taxa, in particular around peri- and sub-Antarctic islands such as South Georgia, which has been dubbed a cryptic species hotspot. This prevalence of cryptic species raises the question of how generalizable these patterns are for Antarctic vertebrates. Here we examine aspects of genotype and phenotype in an Antarctic notothenioid fish species, Lepidonotothen nudifrons, which is distributed in near shore habitats of the Antarctic Peninsula, South Orkney Islands, South Georgia, and the South Sandwich Islands. The results of our analyses show that L. nudifrons comprises two species. We highlight that cryptic species are phenomena not restricted to invertebrate lineages, raising the possibility that the species diversity of notothenioids and other Southern Ocean fishes is under-described. In addition, our findings raise several questions about the evolutionary origin and maintenance of morphological stasis in one of the most extreme habitats on earth.

Parthenogenesis through the ice ages: A biogeographic analysis of Caucasian rock lizards (genus Darevskia).

Darevskia rock lizards include both sexual and parthenogenetic species, mostly distributed in the heterogeneous and ecologically diverse Caucasus. The parthenogenetic species originated via directional hybridogenesis, with only some of the sexual species known to serve as parentals. However, it remains unclear when and where these events happened and how many parental lineages were involved. A multilocus phylogeographic analysis was performed on the parthenogens D. unisexualis, D. bendimahiensis and D. uzzeli, and their putative maternal species D. raddei. Results show the parthenogenetic species all have relatively recent origins, approximately 200-70kyr ago, and at least three hybridization events were involved in their formation. Ecological niche models identify the region where hybridization events leading to the formation of D. unisexualis took place, namely in the northeast of the current distribution. Models also suggest that the sexual D. raddei might have undergone a habitat shift between the Last Interglacial and the Last Glacial Maximum.

Phylogeography of the Vermilion Flycatcher species complex: Multiple speciation events, shifts in migratory behavior, and an apparent extinction of a Galápagos-endemic bird species.

The Vermilion Flycatcher (Pyrocephalus rubinus) is a widespread species found in North and South America and the Galápagos. Its 12 recognized subspecies vary in degree of geographic isolation, phenotypic distinctness, and migratory status. Some authors suggest that Galápagos subspecies nanus and dubius constitute one or more separate species. Observational reports of distinct differences in song also suggest separate species status for the austral migrant subspecies rubinus. To evaluate geographical patterns of diversification and taxonomic limits within this species complex, we carried out a molecular phylogenetic analysis encompassing 10 subspecies and three outgroup taxa using mitochondrial (ND2, Cyt b) and nuclear loci (ODC introns 6 through 7, FGB intron 5). We used samples of preserved tissues from museum collections as well as toe pad samples from museum skins. Galápagos and continental clades were recovered as sister groups, with initial divergence at ∼1mya. Within the continental clade, North and South American populations were sister groups. Three geographically distinct clades were recovered within South America. We detected no genetic differences between two broadly intergrading North American subspecies, mexicanus and flammeus, suggesting they should not be recognized as separate taxa. Four western South American subspecies were also indistinguishable on the basis of loci that we sampled, but occur in a region with patchy habitat, and may represent recently isolated populations. The austral migrant subspecies, rubinus, comprised a monophyletic mitochondrial clade and had many unique nuclear DNA alleles. In combination with its distinct song, exclusive song recognition behavior, different phenology, and an isolated breeding range, our data suggests that this taxon represents a separate species from other continental populations. Mitochondrial and nuclear genetic data, morphology, and behavior suggest that Galápagos forms should be elevated to two full species corresponding to the two currently recognized subspecies, nanus and dubius. The population of dubius is presumed to be extinct, and thus would represent the first documented extinction of a Galápagos-endemic bird species. Two strongly supported mitochondrial clades divide Galápagos subspecies nanus in a geographic pattern that conflicts with previous hypotheses that were based on plumage color. Several populations of nanus have recently become extinct or are in serious decline. Urgent conservation measures should seek to preserve the deep mitochondrial DNA diversity within nanus, and further work should explore whether additional forms should be recognized within nanus. Ancestral states analysis based on our phylogeny revealed that the most recent common ancestor of extant Vermilion Flycatcher populations was migratory, and that migratory behavior was lost more often than gained within Pyrocephalus and close relatives, as has been shown to be the case within Tyrannidae as a whole.