Ancient hybridization and mitochondrial capture between two species of chipmunks

Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe

BACKGROUND

Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry.

METHODS AND RESULTS

MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank.

CONCLUSIONS

Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

Power of Bayesian and Heuristic Tests to Detect Cross-Species Introgression with Reference to Gene Flow in the Tamias quadrivittatus Group of North American Chipmunks

In the past two decades, genomic data have been widely used to detect historical gene flow between species in a variety of plants and animals. The Tamias quadrivittatus group of North America chipmunks, which originated through a series of rapid speciation events, are known to undergo massive amounts of mitochondrial introgression. Yet in a recent analysis of targeted nuclear loci from the group, no evidence for cross-species introgression was detected, indicating widespread cytonuclear discordance. The study used the heuristic method HYDE to detect gene flow, which may suffer from low power. Here we use the Bayesian method implemented in the program BPP to re-analyze these data. We develop a Bayesian test of introgression, calculating the Bayes factor via the Savage-Dickey density ratio using the Markov chain Monte Carlo (MCMC) sample under the model of introgression. We take a stepwise approach to constructing an introgression model by adding introgression events onto a well-supported binary species tree. The analysis detected robust evidence for multiple ancient introgression events affecting the nuclear genome, with introgression probabilities reaching 63%. We estimate population parameters and highlight the fact that species divergence times may be seriously underestimated if ancient cross-species gene flow is ignored in the analysis. We examine the assumptions and performance of HYDE and demonstrate that it lacks power if gene flow occurs between sister lineages or if the mode of gene flow does not match the assumed hybrid-speciation model with symmetrical population sizes. Our analyses highlight the power of likelihood-based inference of cross-species gene flow using genomic sequence data. [Bayesian test; BPP; chipmunks; introgression; MSci; multispecies coalescent; Savage-Dickey density ratio.].

Unexpected genetic integrity boosts hope for the conservation of the red-legged partridge (Alectoris rufa, Galliformes) in Italy

The red-legged partridge (Alectoris rufa) is a medium-sized galliform endemic to southwestern Europe. In the easternmost part of the species' range, the population inhabiting Elba Island (Tuscan Archipelago National Park, Italy) is of undisputed conservation value. While we found nuclear and maternal DNA introgression with the exotic chukar partridge (A. chukar) in previous studies based on microsatellite DNA (n = 25) and two mitochondrial markers (n = 103), respectively, we disclosed a limited or null admixture in a few Elban partridges (n = 4) in a recent genomic investigation relying on 168,675 Single Nucleotide Polymorphisms (SNPs). We herein carried out an extended microsatellite DNA survey including additional 65 samples (total, 90) and six loci (total, 11) to determine both spatial structure and genetic integrity of local A. rufa. A sharp divergence between the subpopulations inhabiting the two sides of the island was disclosed, and the microsatellites indicated that all Elban partridges were not admixed with the chukar, thus fully reflecting the picture inferred using SNPs. We hypothesized that the spreading of chukar genes was constrained by negative selection, with the persistence of only the maternal lineage being indicative of thermal adaptation. The two subpopulations should be treated as distinct Management Units, and an envisaged plan to secure a stock onto nearby Pianosa Island could not only warrant endurance of the Elban population but also establish a source of valuable founders for the ex-situ management of the species in Italy. Our study exemplifies how a suitable samples/loci combination is the key to solve wildlife issues dealing with introgression.

Characterising Mitochondrial Capture in an Iberian Shrew

Mitochondrial introgression raises questions of biogeography and of the extent of reproductive isolation and natural selection. Previous phylogenetic work on the Sorex araneus complex revealed apparent mitonuclear discordance in Iberian shrews, indicating past hybridisation of Sorex granarius and the Carlit chromosomal race of S. araneus, enabling introgression of the S. araneus mitochondrial genome into S. granarius. To further study this, we genetically typed 61 Sorex araneus/coronatus/granarius from localities in Portugal, Spain, France, and Andorra at mitochondrial, autosomal, and sex-linked loci and combined our data with the previously published sequences. Our data are consistent with earlier data indicating that S. coronatus and S. granarius are the most closely related of the three species, confirming that S. granarius from the Central System mountain range in Spain captured the mitochondrial genome from a population of S. araneus. This mitochondrial capture event can be explained by invoking a biogeographical scenario whereby S. araneus was in contact with S. granarius during the Younger Dryas in central Iberia, despite the two species currently having disjunct distributions. We discuss whether selection favoured S. granarius with an introgressed mitochondrial genome. Our data also suggest recent hybridisation and introgression between S. coronatus and S. granarius, as well as between S. araneus and S. coronatus.

Genomic resolution of cryptic species diversity in chipmunks

Discovery of cryptic species is essential to understand the process of speciation and assessing the impacts of anthropogenic stressors. Here, we used genomic data to test for cryptic species diversity within an ecologically well-known radiation of North American rodents, western chipmunks (Tamias). We assembled a de novo reference genome for a single species (Tamias minimus) combined with new and published targeted sequence-capture data for 21,551 autosomal and 493 X-linked loci sampled from 121 individuals spanning 22 species. We identified at least two cryptic lineages corresponding with an isolated subspecies of least chipmunk (T. minimus grisescens) and with a restricted subspecies of the yellow-pine chipmunk (Tamias amoenus cratericus) known only from around the extensive Craters of the Moon lava flow. Additional population-level sequence data revealed that the so-called Crater chipmunk is a distinct species that is abundant throughout the coniferous forests of southern Idaho. This cryptic lineage does not appear to be most closely related to the ecologically and phenotypically similar yellow-pine chipmunk but does show evidence for recurrent hybridization with this and other species.

UCE sequencing-derived mitogenomes reveal the timing of mitochondrial replacement in Malagasy shrew tenrecs (Afrosoricida, Tenrecidae, Microgale)

Malagasy shrew tenrecs (Microgale) have increasingly been used to study speciation genetics over the last years. A previous study recently uncovered gene flow between the Shrew-toothed shrew tenrec (M. soricoides) and sympatric southern population of the Pale shrew tenrec (M. fotsifotsy). This gene flow has been suggested to be accompanied by complete mitochondrial replacement in M. fotsifotsy. To explore the temporal framework of this replacement, we assembled mitogenomes from publicly available sequencing data of ultra-conserved elements. We were able to assemble complete and partial mitogenomes for 19 specimens from five species of shrew tenrecs, which represents a multifold increase in mitogenomic resources available for all tenrecs. Phylogenetic inferences and sequence simulations support the close relationship between the mitochondrial lineages of M. soricoides and the southern population of M. fotsifotsy. Based on the nuclear divergence of northern and southern populations of M. fotsifotsy and the mitochondrial divergence between the latter and M. soricoides, there was a mean time window for replacement of ~ 350,000 years. This timeframe implies that the effective size of the ancestral M. fotsifotsy southern population was less 70,000.

Mitochondrial capture and subsequent genetic divergence generates a novel haplogroup: evidence from ancient and ongoing hybridization in mule and white-tailed deer

Odocoileus virginianus (white-tailed deer) and O. hemionus (mule deer) are sympatric across much of North America. Molecular evidence suggests that up to 24% of individuals in some populations are a product of hybrid ancestry. Several studies have alluded to ancient and recent introgression between Odocoileus spp.; however, no divergence dates were proposed. Herein, phylogenetic analyses of DNA sequences obtained from the mitochondrial Cytochrome b gene in 690 individuals identified three clades corresponding to black-tailed deer, white-tailed deer, or a unique combination of both white-tailed deer and mule deer. White-tailed deer and mule deer diverged from a common ancestor of approximately 3.13 mya followed by an ancient hybridization event of approximately 1.32 mya, in which the white-tailed deer mitochondrial genome was “captured” by mule deer. This hybridization event produced a novel haplogroup for white-tailed deer and mule deer located west of the Appalachian Mountains and east of the Cascade Range, south to Veracruz, Mexico, and north to the Yukon Territory, Canada. The ancestral mule deer-like mitochondrial genome appears to be restricted to black-tailed deer distributed along the western portion of the Cascade and Sierra Nevada Ranges of the United States and Canada, whereas the ancestral white-tailed deer-like mitochondrial genome is restricted to the eastern United States and portions of Latin America and Caribbean regions. The “captured mitochondrial genome” has continued on an independent evolutionary trajectory and represents a unique and broadly distributed haplogroup that is 7.25% and 2.84% different from the ancestral mule deer and ancestral white-tailed deer haplogroups, respectively.

Gene flow in phylogenomics: Sequence capture resolves species limits and biogeography of Afromontane forest endemic frogs from the Cameroon Highlands

Puddle frogs of the Phrynobatrachus steindachneri species complex are a useful group for investigating speciation and phylogeography in Afromontane forests of the Cameroon Volcanic Line, western Central Africa. The species complex is represented by six morphologically relatively cryptic mitochondrial DNA lineages, only two of which are distinguished at the species level - southern P. jimzimkusi and Lake Oku endemic P. njiomock, leaving the remaining four lineages identified as 'P. steindachneri'. In this study, the six mtDNA lineages are subjected to genomic sequence capture analyses and morphological examination to delimit species and to study biogeography. The nuclear DNA data (387 loci; 571,936 aligned base pairs) distinguished all six mtDNA lineages, but the topological pattern and divergence depths supported only four main clades: P. jimzimkusi, P. njiomock, and only two divergent evolutionary lineages within the four 'P. steindachneri' mtDNA lineages. One of the two lineages is herein described as a new species, P. amieti sp. nov. Reticulate evolution (hybridization) was detected within the species complex with morphologically intermediate hybrid individuals placed between the parental species in phylogenomic analyses, forming a ladder-like phylogenetic pattern. The presence of hybrids is undesirable in standard phylogenetic analyses but is essential and beneficial in the network multispecies coalescent. This latter approach provided insight into the reticulate evolutionary history of these endemic frogs. Introgressions likely occurred during the Middle and Late Pleistocene climatic oscillations, due to the cyclic connections (likely dominating during cold glacials) and separations (during warm interglacials) of montane forests. The genomic phylogeographic pattern supports the separation of the southern (Mt. Manengouba to Mt. Oku) and northern mountains at the onset of the Pleistocene. Further subdivisions occurred in the Early Pleistocene, separating populations from the northernmost (Tchabal Mbabo, Gotel Mts.) and middle mountains (Mt. Mbam, Mt. Oku, Mambilla Plateau), as well as the microendemic lineage restricted to Lake Oku (Mt. Oku). This unique model system is highly threatened as all the species within the complex have exhibited severe population declines in the past decade, placing them on the brink of extinction. In addition, Mount Oku is identified to be of particular conservation importance because it harbors three species of this complex. We, therefore, urge for conservation actions in the Cameroon Highlands to preserve their diversity before it is too late.

A Study of Hybridization Between Marmota baibacina and M. sibirica in Their Secondary Contact Zone in Mongolian Altai

The role of hybridization as one of the factors of speciation in mammals has been underestimated for a long time, but now there is a lot of data on its impact in mammalian evolution. Hybridization of species often occurs in their secondary contact zones, which is a natural model for testing factors that ensure species integrity. Studies of hybrid zones are increasingly revealing the essential role of ecological and behavioral features both in initiating crossbreeding and in maintaining interspecific barriers. We studied the hybridization of two species of marmots Marmota baibacina and M. sibirica in the zone of sympatry in Mongolian Altai Mountains. We used a bioacoustic approach to determine the localization of individuals of different species and their cohabitation sites. Genetic typing with two diploid nuclear markers and one marker each of paternal and maternal lines was used to identify hybrids. Habitat preferences of marmots were studied to understand the conditions for the formation of heterospecific pairs. We found a high proportion of hybrid individuals in boulder screes where conditions for the formation of heterospecific pairs probably exist. Our data indicate the viability and fertility of F1 hybrids and their descendants. We hypothesize that the environmental preferences and behavioral features of both species of marmots are important factors that both create conditions for hybridization and limit hybrid dispersal.

Cytonuclear discordance in the crowned-sparrows, Zonotrichia atricapilla and Zonotrichia leucophrys

The golden-crowned (Zonotrichia atricapilla) and white-crowned (Z. leucophrys) sparrows have been presented as a compelling case for rapid speciation. They display divergence in song and plumage with overlap in their breeding ranges implying reproductive isolation, but have almost identical mitochondrial genomes. Previous research proposed hybridization and subsequent mitochondrial introgression as an alternate explanation, but lacked robust nuclear gene trees to distinguish between introgression and incomplete lineage sorting. We test for signatures of these processes between Z. atricapilla and Z. leucophrys, and investigate the relationships among Z. leucophrys subspecies, using mitochondrial sequencing and a reduced representation nuclear genomic dataset. Contrary to the paraphyly evident in mitochondrial gene trees, we confirmed the reciprocal monophyly of Z. atricapilla and Z. leucophrys using large panels of single nucleotide polymorphisms (SNPs). The pattern of cytonuclear discordance is consistent with limited, historical hybridization and mitochondrial introgression, rather than a recent origin and incomplete lineage sorting between recent sister species. We found evidence of nuclear phylogeographic structure within Z. leucophrys with two distinct clades. Altogether, our results indicate deeper divergences between Z. atricapilla and Z. leucophrys than inferred using mitochondrial markers. Our results demonstrate the limitations of relying solely on mitochondrial DNA for taxonomy, and raise questions about the possibility of selection on the mitochondrial genome during temperature oscillations (e.g. during the Pleistocene). Historical mitochondrial introgression facilitated by past environmental changes could cause erroneous dating of lineage splitting in other taxa when based on mitochondrial DNA alone.

Complex histories of gene flow and a mitochondrial capture event in a nonsister pair of birds

Hybridization, introgression, and reciprocal gene flow during speciation, specifically the generation of mitonuclear discordance, are increasingly observed as parts of the speciation process. Genomic approaches provide insight into where, when, and how adaptation operates during and after speciation and can measure historical and modern introgression. Whether adaptive or neutral in origin, hybridization can cause mitonuclear discordance by placing the mitochondrial genome of one species (or population) in the nuclear background of another species. The latter, introgressed species may eventually have its own mtDNA replaced or “captured” by other species across its entire geographical range. Intermediate stages in the capture process should be observable. Two nonsister species of Australasian monarch‐flycatchers, Spectacled Monarch (Symposiachrus trivirgatus) mostly of Australia and Indonesia and Spot‐winged Monarch (S. guttula) of New Guinea, present an opportunity to observe this process. We analysed thousands of single nucleotide polymorphisms (SNPs) derived from ultraconserved elements of all subspecies of both species. Mitochondrial DNA sequences of Australian populations of S. trivirgatus form two paraphyletic clades, one being sister to and presumably introgressed by S. guttula despite little nuclear signal of introgression. Population genetic analyses (e.g., tests for modern and historical gene flow and selection) support at least one historical gene flow event between S. guttula and Australian S. trivirgatus. We also uncovered introgression from the Maluku Islands subspecies of S. trivirgatus into an island population of S. guttula, resulting in apparent nuclear paraphyly. We find that neutral demographic processes, not adaptive introgression, are the most likely cause of these complex population histories. We suggest that a Pleistocene extinction of S. guttula from mainland Australia resulted from range expansion by S. trivirgatus.

Estimating hybridization as a consequence of climatic fluctuations: A case study of two geometridae species

In some cases, the phylogenetic analysis based on the mitochondrial DNA (mtDNA) and the nuclear DNA (ncDNA) are discordant. There are three major causes of the discordance within insects, including hybridization, incomplete lineage sorting (ILS) and infection by Wolbachia. In this study, we used a combination of multilocus and coalescent analyses to explore these processes occurred during the evolutionary history of Limbatochlamys rosthorni Rothschild, 1894 and Limbatochlamys pararosthorni Han and Xue, 2005. The ncDNA phylogenetic tree supported two reciprocally monophyletic species, whereas the mtDNA results failed to reveal such a structure and revealed an extensive level of admixture between two species. Because of very low Wolbachia infection rates (<20%), we firstly excluded this reason for the discordance. The fixed nucleotide differences and large genetic distances (1.5-2.5%) at the ncDNA genes suggested that the lineage sorting process between these two species is nearly complete and two species have experienced a prolonged period of independent evolution. Thus, we secondly excluded ILS. Sharing haplotypes, mtDNA gene flow occurring and the transitional samples with morphological features supported hybridization. The distribution contraction during glaciations and postglacial distribution expansion might have facilitated hybridization. Taken together, our study indicates that the current genetic structure of L. rosthorni and L. pararosthorni is the results of contraction and fragmentation into separated refugia during glaciations, followed by postglacial expansion and admixture.

Diversification, Introgression, and Rampant Cytonuclear Discordance in Rocky Mountains Chipmunks (Sciuridae: Tamias)

Evidence from natural systems suggests that hybridization between animal species is more common than traditionally thought, but the overall contribution of introgression to standing genetic variation within species remains unclear for most animal systems. Here, we use targeted exon-capture to sequence thousands of nuclear loci and complete mitochondrial genomes from closely related chipmunk species in the Tamias quadrivittatus group that are distributed across the Great Basin and the central and southern Rocky Mountains of North America. This recent radiation includes six overlapping, ecologically distinct species (T. canipes, T. cinereicollis, T. dorsalis, T. quadrivittatus, T. rufus, and T. umbrinus) that show evidence for widespread introgression across species boundaries. Such evidence has historically been derived from a handful of markers, typically focused on mitochondrial loci, to describe patterns of introgression; consequently, the extent of introgression of nuclear genes is less well characterized. We conducted a series of phylogenomic and species-tree analyses to resolve the phylogeny of six species in this group. In addition, we performed several population genomic analyses to characterize nuclear genomes and infer coancestry among individuals. Furthermore, we used emerging quartets-based approaches to simultaneously infer the species tree (SVDquartets) and identify introgression (HyDe). We found that, in spite of rampant introgression of mitochondrial genomes between some species pairs (and sometimes involving up to three species), there appears to be little to no evidence for nuclear introgression. These findings mirror other genomic results where complete mitochondrial capture has occurred between chipmunk species in the absence of appreciable nuclear gene flow. The underlying causes of recurrent massive cytonuclear discordance remain unresolved in this group but mitochondrial DNA appears highly misleading of population histories as a whole. Collectively, it appears that chipmunk species boundaries are largely impermeable to nuclear gene flow and that hybridization, while pervasive with respect to mtDNA, has likely played a relatively minor role in the evolutionary history of this group.

Disentangling lousy relationships: Comparative phylogenomics of two sucking louse lineages parasitizing chipmunks

The evolution of obligate parasites is often interpreted in light of their hosts' evolutionary history. An expanded approach is to examine the histories of multiple lineages of parasites that inhabit similar environments on a particular host lineage. Western North American chipmunks (genus Tamias) have a broad distribution, a history of divergence with gene flow, and host two species of sucking lice (Anoplura), Hoplopleura arboricola and Neohaematopinus pacificus. From total genomic sequencing, we obtained sequences of over 1100 loci sampled across the genomes of these lice to compare their evolutionary histories and examine the roles of host association in structuring louse relationships. Within each louse species, clades are largely associated with closely related chipmunk host species. Exceptions to this pattern appear to have a biogeographic component, but differ between the two louse species. Phylogenetic relationships among these major louse clades, in both species, are not congruent with chipmunk relationships. In the context of host associations, each louse lineage has a different evolutionary history, supporting the hypothesis that host-parasite assemblages vary both across the landscape and with the taxa under investigation. In addition, the louse Hoplopleura erratica (parasitizing the eastern Tamias striatus) is embedded within H. arboricola, rendering it paraphyletic. This phylogenetic result, together with comparable divergences within H. arboricola, indicate a need for taxonomic revision. Both host divergence and biogeographic components shape parasite diversification as demonstrated by the distinctive diversification patterns of these two independently evolving lineages that parasitize the same hosts.

No evidence for phylosymbiosis in western chipmunk species

Phylosymbiosis refers to a congruent pattern between the similarity of microbiomes of different species and the branching pattern of the host phylogeny. Phylosymbiosis has been detected in a variety of vertebrate and invertebrate hosts, but has only been assessed in geographically isolated populations. We tested for phylosymbiosis in eight (sub)species of western chipmunks with overlapping ranges and ecological niches; we used a nuclear (Acrosin) and a mitochondrial (CYTB) phylogenetic marker because there are many instances of mitochondrial introgression in chipmunks. We predicted that similarity among microbiomes increases with: (1) increasing host mitochondrial relatedness, (2) increasing host nuclear genome relatedness and (3) decreasing geographic distance among hosts. We did not find statistical evidence supporting phylosymbiosis in western chipmunks. Furthermore, in contrast to studies of other mammalian microbiomes, similarity of chipmunk microbiomes is not predominantly determined by host species. Sampling site explained most variation in microbiome composition, indicating an important role of local environment in shaping microbiomes. Fecal microbiomes of chipmunks were dominated by Bacteroidetes (72.2%), followed by Firmicutes (24.5%), which is one of the highest abundances of Bacteroidetes detected in wild mammals. Future work will need to elucidate the effects of habitat, ecology and host genomics on chipmunk microbiomes.

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.

Beyond DNA barcoding: The unrealized potential of genome skim data in sample identification

Genetic tools are increasingly used to identify and discriminate between species. One key transition in this process was the recognition of the potential of the ca 658bp fragment of the organelle cytochrome c oxidase I (COI) as a barcode region, which revolutionized animal bioidentification and lead, among others, to the instigation of the Barcode of Life Database (BOLD), containing currently barcodes from >7.9 million specimens. Following this discovery, suggestions for other organellar regions and markers, and the primers with which to amplify them, have been continuously proposed. Most recently, the field has taken the leap from PCR-based generation of DNA references into shotgun sequencing-based "genome skimming" alternatives, with the ultimate goal of assembling organellar reference genomes. Unfortunately, in genome skimming approaches, much of the nuclear genome (as much as 99% of the sequence data) is discarded, which is not only wasteful, but can also limit the power of discrimination at, or below, the species level. Here, we advocate that the full shotgun sequence data can be used to assign an identity (that we term for convenience its "DNA-mark") for both voucher and query samples, without requiring any computationally intensive pretreatment (e.g. assembly) of reads. We argue that if reference databases are populated with such "DNA-marks," it will enable future DNA-based taxonomic identification to complement, or even replace PCR of barcodes with genome skimming, and we discuss how such methodology ultimately could enable identification to population, or even individual, level.

One Hundred Mitochondrial Genomes of Cicadas

Mitochondrial genomes can provide valuable information on the biology and evolutionary histories of their host organisms. Here, we present and characterize the complete coding regions of 107 mitochondrial genomes (mitogenomes) of cicadas (Insecta: Hemiptera: Auchenorrhyncha: Cicadoidea), representing 31 genera, 61 species, and 83 populations. We show that all cicada mitogenomes retain the organization and gene contents thought to be ancestral in insects, with some variability among cicada clades in the length of a region between the genes nad2 and cox1, which encodes 3 tRNAs. Phylogenetic analyses using these mitogenomes recapitulate a recent 5-gene classification of cicadas into families and subfamilies, but also identify a species that falls outside of the established taxonomic framework. While protein-coding genes are under strong purifying selection, tests of relative evolutionary rates reveal significant variation in evolutionary rates across taxa, highlighting the dynamic nature of mitochondrial genome evolution in cicadas. These data will serve as a useful reference for future research into the systematics, ecology, and evolution of the superfamily Cicadoidea.

The efficiency of universal mitochondrial DNA barcodes for species discrimination of Pomacea canaliculata and Pomacea maculata

Invasive apple snails, Pomacea canaliculata and P. maculata, have a widespread distribution globally and are regarded as devastating pests of agricultural wetlands. The two species are morphologically similar, which hinders species identification via morphological approaches and species-specific management efforts. Advances in molecular genetics may contribute effective diagnostic tools to potentially resolve morphological ambiguity. DNA barcoding has revolutionized the field of taxonomy by providing an alternative, simple approach for species discrimination, where short sections of DNA, the cytochrome c oxidase subunit I (COI) gene in particular, are used as ‘barcodes’ to delineate species boundaries. In our study, we aimed to assess the effectiveness of two mitochondrial markers, the COI and 16S ribosomal deoxyribonucleic acid (16S rDNA) markers for DNA barcoding of P. canaliculata and P. maculata. The COI and 16S rDNA sequences of 40 Pomacea specimens collected from six localities in Peninsular Malaysia were analyzed to assess their barcoding performance using phylogenetic methods and distance-based assessments. The results confirmed both markers were suitable for barcoding P. canaliculata and P. maculata. The phylogenies of the COI and 16S rDNA markers demonstrated species-specific monophyly and were largely congruent with the exception of one individual. The COI marker exhibited a larger barcoding gap (6.06–6.58%) than the 16S rDNA marker (1.54%); however, the magnitude of barcoding gap generated within the barcoding region of the 16S rDNA marker (12-fold) was bigger than the COI counterpart (approximately 9-fold). Both markers were generally successful in identifying P. canaliculata and P. maculata in the similarity-based DNA identifications. The COI + 16S rDNA concatenated dataset successfully recovered monophylies of P. canaliculata and P. maculata but concatenation did not improve individual datasets in distance-based analyses. Overall, although both markers were successful for the identification of apple snails, the COI molecular marker is a better barcoding marker and could be utilized in various population genetic studies of P. canaliculata and P. maculata.

Genome-wide SNPs clarify lineage diversity confused by coloration in coralsnakes of the Micrurus diastema species complex (Serpentes: Elapidae)

New world coralsnakes of the genus Micrurus are a diverse radiation of highly venomous and brightly colored snakes that range from North Carolina to Argentina. Species in this group have played central roles in developing and testing hypotheses about the evolution of mimicry and aposematism. Despite their diversity and prominence as model systems, surprisingly little is known about species boundaries and phylogenetic relationships within Micrurus, which has substantially hindered meaningful analyses of their evolutionary history. Here we use mitochondrial genes together with thousands of nuclear genomic loci obtained via ddRADseq to study the phylogenetic relationships and population genomics of a subclade of the genus Micrurus: The M. diastema species complex. Our results indicate that prior species and species-group inferences based on morphology and color pattern have grossly misguided taxonomy, and that the M. diastema complex is not monophyletic. Based on our analyses of molecular data, we infer the phylogenetic relationships among species and populations, and provide a revised taxonomy for the group. Two non-sister species-complexes with similar color patterns are recognized, the M. distans and the M. diastema complexes, the first being basal to the monadal Micrurus and the second encompassing most North American monadal taxa. We examined all 13 species, and their respective subspecies, for a total of 24 recognized taxa in the M. diastema species complex. Our analyses suggest a reduction to 10 species, with no subspecific designations warranted, to be a more likely estimate of species diversity, namely, M. apiatus, M. browni, M. diastema, M. distans, M. ephippifer, M. fulvius, M. michoacanensis, M. oliveri, M. tener, and one undescribed species.

Temporal genomic contrasts reveal rapid evolutionary responses in an alpine mammal during recent climate change

Many species have experienced dramatic changes in their abundance and distribution during recent climate change, but it is often unclear whether such ecological responses are accompanied by evolutionary change. We used targeted exon sequencing of 294 museum specimens (160 historic, 134 modern) to generate independent temporal genomic contrasts spanning a century of climate change (1911-2012) for two co-distributed chipmunk species: an endemic alpine specialist (Tamias alpinus) undergoing severe range contraction and a stable mid-elevation species (T. speciosus). Using a novel analytical approach, we reconstructed the demographic histories of these populations and tested for evidence of recent positive directional selection. Only the retracting species showed substantial population genetic fragmentation through time and this was coupled with positive selection and substantial shifts in allele frequencies at a gene, Alox15, involved in regulation of inflammation and response to hypoxia. However, these rapid population and gene-level responses were not detected in an analogous temporal contrast from another area where T. alpinus has also undergone severe range contraction. Collectively, these results highlight that evolutionary responses may be variable and context dependent across populations, even when they show seemingly synchronous ecological shifts. Our results demonstrate that temporal genomic contrasts can be used to detect very recent evolutionary responses within and among contemporary populations, even in the face of complex demographic changes. Given the wealth of specimens archived in natural history museums, comparative analyses of temporal population genomic data have the potential to improve our understanding of recent and ongoing evolutionary responses to rapidly changing environments.

DNA Barcoding and Taxonomic Challenges in Describing New Putative Species: Examples from Sootywing and Cloudywing Butterflies (Lepidoptera: Hesperiidae)

DNA barcoding has resulted in the ‘discovery’ of a vast number of new species and subspecies. Assigning formal scientific names to these taxa remains a major challenge. Names sometimes are newly designated. Alternatively, available valid names can be resurrected from synonymy, based on barcode analyses together with classical taxonomic characters. For the most part, however, new putative species revealed by barcoding studies go undescribed. This situation is most often attributed to insufficient taxonomic expertise with the authors conducting the study, together with a critical lack of formally trained taxonomists. However, even with formal training, and additional supportive data from morphological, ecological or life history characters, other factors can arise that impede new species descriptions. In the present paper, several specific taxonomic challenges that have arisen from barcode analyses in two groups of skipper butterflies (Lepidoptera: Hesperiidae), the Sootywings (Pholisora catullus and P. mejicanus) and the Coyote Cloudywing (Achalarus toxeus) are highlighted and discussed. Both P. catullus and A. toxeus show relatively large intraspecific genetic divergences of barcodes (2–3%) which suggests the possibility of previously unrecognized cryptic speciation within each group. Some of the challenges to providing formal names and clarifying taxonomic status of these cryptic taxa could be largely overcome by (1) barcoding type specimens, (2) clarifying imprecise and often vague or suspect type localities, and (3) by conducting in-depth comparative studies on genitalic morphology.

The genomic impact of historical hybridization with massive mitochondrial DNA introgression

BACKGROUND

The extent to which selection determines interspecific patterns of genetic exchange enlightens the role of adaptation in evolution and speciation. Often reported extensive interspecific introgression could be selection-driven, but also result from demographic processes, especially in cases of invasive species replacements, which can promote introgression at their invasion front. Because invasion and selective sweeps similarly mold variation, population genetics evidence for selection can only be gathered in an explicit demographic framework. The Iberian hare, Lepus granatensis, displays in its northern range extensive mitochondrial DNA introgression from L. timidus, an arctic/boreal species that it replaced locally after the last glacial maximum. We use whole-genome sequencing to infer geographic and genomic patterns of nuclear introgression and fit a neutral model of species replacement with hybridization, allowing us to evaluate how selection influenced introgression genome-wide, including for mtDNA.

RESULTS

Although the average nuclear and mtDNA introgression patterns contrast strongly, they fit a single demographic model of post-glacial invasive replacement of timidus by granatensis. Outliers of elevated introgression include several genes related to immunity, spermatogenesis, and mitochondrial metabolism. Introgression is reduced on the X chromosome and in low recombining regions.

CONCLUSIONS

General nuclear and mtDNA patterns of introgression can be explained by purely demographic processes. Hybrid incompatibilities and interplay between selection and recombination locally modulate levels of nuclear introgression. Selection promoted introgression of some genes involved in conflicts, either interspecific (parasites) or possibly cytonuclear. In the latter case, nuclear introgression could mitigate the potential negative effects of alien mtDNA on mitochondrial metabolism and male-specific traits.

Mito-nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi-permeable. Using mitochondrial cytochrome b sequences and genome-wide loci developed via ddRAD-seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5-54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.

Our love-hate relationship with DNA barcodes, the Y2K problem, and the search for next generation barcodes

DNA barcodes are very useful for species identification especially when identification by traditional morphological characters is difficult. However, the short mitochondrial and chloroplast barcodes currently in use often fail to distinguish between closely related species, are prone to lateral transfer, and provide inadequate phylogenetic resolution, particularly at deeper nodes. The deficiencies of short barcode identifiers are similar to the deficiencies of the short year identifiers that caused the Y2K problem in computer science. The resolution of the Y2K problem was to increase the size of the year identifiers. The performance of conventional mitochondrial COI barcodes for phylogenetics was compared with the performance of complete mitochondrial genomes and nuclear ribosomal RNA repeats obtained by genome skimming for a set of caddisfly taxa (Insect Order Trichoptera). The analysis focused on Trichoptera Family Hydropsychidae, the net-spinning caddisflies, which demonstrates many of the frustrating limitations of current barcodes. To conduct phylogenetic comparisons, complete mitochondrial genomes (15 kb each) and nuclear ribosomal repeats (9 kb each) from six caddisfly species were sequenced, assembled, and are reported for the first time. These sequences were analyzed in comparison with eight previously published trichopteran mitochondrial genomes and two triochopteran rRNA repeats, plus outgroup sequences from sister clade Lepidoptera (butterflies and moths). COI trees were not well-resolved, had low bootstrap support, and differed in topology from prior phylogenetic analyses of the Trichoptera. Phylogenetic trees based on mitochondrial genomes or rRNA repeats were well-resolved with high bootstrap support and were largely congruent with each other. Because they are easily sequenced by genome skimming, provide robust phylogenetic resolution at various phylogenetic depths, can better distinguish between closely related species, and (in the case of mitochondrial genomes), are backwards compatible with existing mitochondrial barcodes, it is proposed that mitochondrial genomes and rRNA repeats be used as next generation DNA barcodes.

Comparative Phylogenomic Assessment of Mitochondrial Introgression among Several Species of Chipmunks (Tamias)

Many species are not completely reproductively isolated, resulting in hybridization and genetic introgression. Organellar genomes, such as those derived from mitochondria (mtDNA) and chloroplasts, introgress frequently in natural systems; however, the forces shaping patterns of introgression are not always clear. Here, we investigate extensive mtDNA introgression in western chipmunks, focusing on species in the Tamias quadrivittatus group from the central and southern Rocky Mountains. Specifically, we investigate the role of selection in driving patterns of introgression. We sequenced 51 mtDNA genomes from six species and combine these sequences with other published genomic data to yield annotated mitochondrial reference genomes for nine species of chipmunks. Genomic characterization was performed using a series of molecular evolutionary and phylogenetic analyses to test protein-coding genes for positive selection. We fit a series of maximum likelihood models using a model-averaging approach, assessed deviations from neutral expectations, and performed additional tests to search for codons under the influence of selection. We found no evidence for positive selection among these genomes, suggesting that selection has not been the driving force of introgression in these species. Thus, extensive mtDNA introgression among several species of chipmunks likely reflects genetic drift of introgressed alleles in historically fluctuating populations.

The genomic ancestry, landscape genetics and invasion history of introduced mice in New Zealand

The house mouse (Mus musculus) provides a fascinating system for studying both the genomic basis of reproductive isolation, and the patterns of human-mediated dispersal. New Zealand has a complex history of mouse invasions, and the living descendants of these invaders have genetic ancestry from all three subspecies, although most are primarily descended from M. m. domesticus. We used the GigaMUGA genotyping array (approximately 135 000 loci) to describe the genomic ancestry of 161 mice, sampled from 34 locations from across New Zealand (and one Australian city-Sydney). Of these, two populations, one in the south of the South Island, and one on Chatham Island, showed complete mitochondrial lineage capture, featuring two different lineages of M. m. castaneus mitochondrial DNA but with only M. m. domesticus nuclear ancestry detectable. Mice in the northern and southern parts of the North Island had small traces (approx. 2-3%) of M. m. castaneus nuclear ancestry, and mice in the upper South Island had approximately 7-8% M. m. musculus nuclear ancestry including some Y-chromosomal ancestry-though no detectable M. m. musculus mitochondrial ancestry. This is the most thorough genomic study of introduced populations of house mice yet conducted, and will have relevance to studies of the isolation mechanisms separating subspecies of mice.

Ancient Mitochondrial Capture as Factor Promoting Mitonuclear Discordance in Freshwater Fishes: A Case Study in the Genus Squalius (Actinopterygii, Cyprinidae) in Greece

Hybridization and incomplete lineage sorting are common confounding factors in phylogeny and speciation resulting in mitonuclear disparity. Mitochondrial introgression, a particular case of hybridization, may, in extreme cases, lead to replacement of the mitochondrial genome of one species with that of another (mitochondrial capture). We investigated mitochondrial introgression involving two species of the cyprinid genus Squalius in the western Peloponnese region of Greece using molecular and morphological data. We found evidence of complete mitochondrial introgression of Squalius keadicus into two populations recognized as Squalius peloponensis from the Miras and Pamissos River basins and a divergence of mitochondrial genomes of S. keadicus from the Evrotas basin from that of the introgressed populations dating from the Pleistocene. Secondary contact among basins is a possible factor in connection of the species and the introgression event. Morphological analyses support the hypothesis of mitochondrial introgression, as S. keadicus was different from the other three populations recognized as S. peloponensis, although significant differences were found among the four populations. Isolation by geographical barriers arose during Pleistocene in the western Peloponnese were the source of the evolution of the two reciprocally monophyletic subclades found in the S. keadicus mitochondrial clade, and the morphological differences found among the four populations. Along with the lack of structure in the nuclear genome in the three populations ascribed to S. peloponensis, this suggests an incipient speciation process occurring in these Squalius species in the western Peloponnese.

Expected pairwise congruence among gene trees under the coalescent model

Although it is widely appreciated that gene trees may differ from the overall species tree and from one another due to various evolutionary processes (e.g., incomplete lineage sorting (ILS), horizontal gene transfer, etc.), the extent of this incongruence is rarely quantified and discussed. Here we consider the expected amount of incongruence arising from ILS, as modeled by the coalescent process. In particular, we compute the probability that two gene trees randomly sampled from the same species tree agree with one another as well as the distribution of the Robinson-Foulds distance between them, for species trees with three to eight taxa. We demonstrate that, as expected under the coalescent model, the amount of discordance is affected by species tree-specific factors such as speciation times and effective population sizes for the species under consideration. Our results highlight the fact that substantial discordance may occur, even when the number of species is very small, which has implications both for larger taxon samples and for any method that uses estimated gene trees as the basis for further statistical inference. The amount of incongruence is substantial enough that such methods may need to be modified to account for variability in the underlying gene trees.

Temporal and spatial mosaics: deep host association and shallow geographic drivers shape genetic structure in a widespread pinworm, Rauschtineria eutamii

Climate and host demographic cycling often shape both parasite genetic diversity and host distributions, processes that transcend a history of strict host-parasite association. We explored host associations and histories based on an evaluation of mitochondrial and nuclear sequences to reveal the underlying history and genetic structure of a pinworm, Rauschtineria eutamii, infecting 10 species of western North American chipmunks (Rodentia:Tamias, subgenus Neotamias). Rauschtineria eutamii contains divergent lineages influenced by the diversity of hosts and variation across the complex topography of western North America. We recovered six reciprocally monophyletic R. eutamii mitochondrial clades, largely supported by nuclear gene trees, exhibiting divergence levels comparable to intraspecific variation reported for other nematodes. Phylogenetic relationships among pinworm clades suggest that R. eutamii colonized an ancestral lineage of western chipmunks and lineages persisted during historical isolation in diverging Neotamias species or species groups. Pinworm diversification, however, is incongruent and asynchronous relative to host diversification. Secondarily, patterns of shallow divergence were shaped by geography through events of episodic colonization reflecting an interaction of taxon pulses and ecological fitting among assemblages in recurrent sympatry. Pinworms occasionally infect geographically proximal host species; however, host switching may be unstable or ephemeral, as there is no signal of host switching in the deeper history of R. eutamii.

Evolutionary history of endemic Sulawesi squirrels constructed from UCEs and mitogenomes sequenced from museum specimens

Background

The Indonesian island of Sulawesi has a complex geological history. It is composed of several landmasses that have arrived at a near modern configuration only in the past few million years. It is the largest island in the biodiversity hotspot of Wallacea—an area demarcated by the biogeographic breaks between Wallace’s and Lydekker’s lines. The mammal fauna of Sulawesi is transitional between Asian and Australian faunas. Sulawesi’s three genera of squirrels, all endemic (subfamily Nannosciurinae: Hyosciurus, Rubrisciurus and Prosciurillus), are of Asian origin and have evolved a variety of phenotypes that allow a range of ecological niche specializations. Here we present a molecular phylogeny of this radiation using data from museum specimens. High throughput sequencing technology was used to generate whole mitochondrial genomes and a panel of nuclear ultraconserved elements providing a large genome-wide dataset for inferring phylogenetic relationships.

Results

Our analysis confirmed monophyly of the Sulawesi taxa with deep divergences between the three endemic genera, which predate the amalgamation of the current island of Sulawesi. This suggests lineages may have evolved in allopatry after crossing Wallace’s line. Nuclear and mitochondrial analyses were largely congruent and well supported, except for the placement of Prosciurillus murinus. Mitochondrial analysis revealed paraphyly for Prosciurillus, with P. murinus between or outside of Hyosciurus and Rubrisciurus, separate from other species of Prosciurillus. A deep but monophyletic history for the four included species of Prosciurillus was recovered with the nuclear data.

Conclusions

The divergence of the Sulawesi squirrels from their closest relatives dated to ~9.7–12.5 million years ago (MYA), pushing back the age estimate of this ancient adaptive radiation prior to the formation of the current conformation of Sulawesi. Generic level diversification took place around 9.7 MYA, opening the possibility that the genera represent allopatric lineages that evolved in isolation in an ancient proto-Sulawesian archipelago. We propose that incongruence between phylogenies based on nuclear and mitochondrial sequences may have resulted from biogeographic discordance, when two allopatric lineages come into secondary contact, with complete replacement of the mitochondria in one species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0650-z) contains supplementary material, which is available to authorized users.

Nuclear introgression without mitochondrial introgression in two turtle species exhibiting sex-specific trophic differentiation

Despite the presence of reproductive barriers between species, interspecific gene introgression has been documented in a range of natural systems. Comparing patterns of genetic introgression in biparental versus matrilineal markers can potentially reveal sex-specific barriers to interspecific gene flow. Hybridization has been documented in the freshwater turtles Graptemys geographica and G. pseudogeographica, whose ranges are largely sympatric. Morphological differentiation between the species is restricted to females, with female G. geographica possessing large heads and jaws compared to the narrow heads of G. pseudogeographica females. If hybrid females are morphologically intermediate, they may be less successful at exploiting parental feeding niches, thereby limiting the introgression of maternally inherited, but not biparental, molecular markers. We paired sequence data with stable isotope analysis and examined sex-specific genetic introgression and trophic differentiation in sympatric populations of G. geographica and G. pseudogeographica. We observed introgression from G. pseudogeographica into G. geographica at three nuclear loci, but not at the mitochondrial locus. Analysis of ∂(15)N and ∂(13)C was consistent with species differences in trophic positioning in females, but not males. These results suggest that ecological divergence in females may reduce the opportunity for gene flow in this system.

Revision of widespread red squirrels (genus: Tamiasciurus) highlights the complexity of speciation within North American forests

Integration of molecular methods, ecological modeling, and statistical hypothesis testing are increasing our understanding of differentiation within species and phylogenetic relationships among species by revealing environmental connections to evolutionary processes. Within mammals, novel diversity is being discovered and characterized as more complete geographic sampling is coupled with newer multi-disciplinary approaches. North American red squirrels exemplify a forest obligate genus whose species are monitored as indicators of forest ecosystem condition, yet phylogenetic relationships reflecting evolutionary history within this genus remain tentative. Through testing of competing systematic and niche-based divergence hypotheses, we recognize three species, Tamiasciurus douglasii, T. hudsonicus, and T. fremonti. Our data provide evidence of regional differences in evolutionary dynamics and continental gradients of complexity that are important both for future management and for investigating multiple pathways that can lead to the formation of new species.

Towards a uniform nomenclature for ground squirrels: the status of the Holarctic chipmunks

The chipmunks are a Holarctic group of ground squirrels currently allocated to the genus

Phylogenomic evidence for ancient hybridization in the genomes of living cats (Felidae)

Inter-species hybridization has been recently recognized as potentially common in wild animals, but the extent to which it shapes modern genomes is still poorly understood. Distinguishing historical hybridization events from other processes leading to phylogenetic discordance among different markers requires a well-resolved species tree that considers all modes of inheritance and overcomes systematic problems due to rapid lineage diversification by sampling large genomic character sets. Here, we assessed genome-wide phylogenetic variation across a diverse mammalian family, Felidae (cats). We combined genotypes from a genome-wide SNP array with additional autosomal, X- and Y-linked variants to sample ∼150 kb of nuclear sequence, in addition to complete mitochondrial genomes generated using light-coverage Illumina sequencing. We present the first robust felid time tree that accounts for unique maternal, paternal, and biparental evolutionary histories. Signatures of phylogenetic discordance were abundant in the genomes of modern cats, in many cases indicating hybridization as the most likely cause. Comparison of big cat whole-genome sequences revealed a substantial reduction of X-linked divergence times across several large recombination cold spots, which were highly enriched for signatures of selection-driven post-divergence hybridization between the ancestors of the snow leopard and lion lineages. These results highlight the mosaic origin of modern felid genomes and the influence of sex chromosomes and sex-biased dispersal in post-speciation gene flow. A complete resolution of the tree of life will require comprehensive genomic sampling of biparental and sex-limited genetic variation to identify and control for phylogenetic conflict caused by ancient admixture and sex-biased differences in genomic transmission.

Barcoding gap, but no support for cryptic speciation in the earthworm Aporrectodea longa (Clitellata: Lumbricidae)

DNA-barcoding, using the mitochondrial marker COI, has been found successful for the identification of specimens in many animal groups, but may not be suited for species discovery and delimitation if used alone. In this study, we investigate whether two observed COI haplogroups in the earthworm Aporrectodea longa correspond to two cryptic species or if the variation is intraspecific. This is done by complementing COI with two nuclear markers, ITS2 and Histone 3. The variation is studied using distance methods, parsimony networks and Bayesian coalescent trees, and the statistical distinctness of the groups is tested on gene trees using the genealogical sorting index, Rosenberg's P_AB and Rodrigo et al.'s P_(RD). We also applied multilocus species delimitation based on the multispecies coalescence model. The two haplogroups were found in COI, and all tests except P_(RD) found them to be significantly distinct. However, in ITS2, the same groups were not recovered in any analyses or tests. H3 was invariable in A. longa, and was, therefore, included only in the multilocus analysis, which preferred a model treating A. longa as one species over a model splitting it into two. We also compared two measurements of size, body length, and no. of segments between the groups. No difference in body length was found, and although a significant difference in no. of segments was noted the haplogroup with the lower mean showed both the highest and the lowest value. When combined, these results led us to the conclusion that there is no support for the separation of A. longa into two cryptic species. This study again highlights the importance of complementing mitochondrial barcodes with more data when establishing species boundaries.

A protocol for targeted enrichment of intron-containing sequence markers for recent radiations: A phylogenomic example from Heuchera (Saxifragaceae)

PREMISE OF THE STUDY

Phylogenetic inference is moving to large multilocus data sets, yet there remains uncertainty in the choice of marker and sequencing method at low taxonomic levels. To address this gap, we present a method for enriching long loci spanning intron-exon boundaries in the genus Heuchera.

METHODS

Two hundred seventy-eight loci were designed using a splice-site prediction method combining transcriptomic and genomic data. Biotinylated probes were designed for enrichment of these loci. Reference-based assembly was performed using genomic references; additionally, chloroplast and mitochondrial genomes were used as references for off-target reads. The data were aligned and subjected to coalescent and concatenated phylogenetic analyses to demonstrate support for major relationships.

RESULTS

Complete or nearly complete (>99%) sequences were assembled from essentially all loci from all taxa. Aligned introns showed a fourfold increase in divergence as opposed to exons. Concatenated analysis gave decisive support to all nodes, and support was also high and relationships mostly similar in the coalescent analysis. Organellar phylogenies were also well-supported and conflicted with the nuclear signal.

DISCUSSION

Our approach shows promise for resolving a recent radiation. Enrichment for introns is highly successful with little or no sequencing dropout at low taxonomic levels despite higher substitution and indel frequencies, and should be exploited in studies of species complexes.

Complex history of isolation and gene flow in hoary, Olympic, and endangered Vancouver Island marmots

Climate change resulting in a reduction of alpine habitat is believed to pose a considerable risk to alpine-dependent species, including many marmots. Hoary marmots (Marmota caligata) range throughout much of the mountainous Pacific Northwest (PNW) and Rocky Mountains while the closely related Olympic and Vancouver Island marmots (M. olympus and M. vancouverensis, respectively) are restricted to small isolated regions of the PNW. The endemic Vancouver Island marmot is currently classified as Critically Endangered and the Olympic marmot has recently experienced dramatic population declines. Previous phylogenetic studies of PNW marmot species have had limited power as they focused on resolving interspecific relationships, implicitly assumed an absence of gene flow among currently recognized species, included relatively few individuals, and relied heavily or entirely on mitochondrial DNA. We sequenced 2 mitochondrial and 4 nuclear markers from 167 hoary, 4 Vancouver Island, and 5 Olympic marmots in order to investigate phylogenetic relationships and historic gene flow among these species. We recovered 2 monophyletic (and predominantly allopatric) mitochondrial clades of hoary marmots that are not sister groups. Instead, Vancouver Island marmots formed a monophyletic mitochondrial sister clade to 1 of the hoary marmot clades. Nuclear loci did not recover the 2 mitochondrial clades of hoary marmots and suggest that Vancouver Island marmots may have experienced mitochondrial introgression from coastal mainland hoary marmots. Additionally, our nuclear results suggest possible gene flow between hoary and Olympic marmots despite different chromosomal formulas. Rather than resolving what has previously been considered a straightforward 3-taxon phylogenetic question, our findings suggest a complicated history of rapid divergence of the 3 species followed by intermittent and possibly ongoing gene flow between hoary marmots and both Olympic and Vancouver Island marmots. These results therefore have significant implications for the conservation of the latter 2 species, both of which are conservation concerns.