Jeffrey W. Streicher

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake.

Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct “Seq-to-SSR” approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus, and Clarks Nutcracker, Nucifraga columbiana), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as

How Should Genes and Taxa be Sampled for Phylogenomic Analyses with Missing Data? An Empirical Study in Iguanian Lizards

10 per sample – a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.

Diversification and asymmetrical gene flow across time and space: lineage sorting and hybridization in polytypic barking frogs

Targeted sequence capture is becoming a widespread tool for generating large phylogenomic data sets to address difficult phylogenetic problems. However, this methodology often generates data sets in which increasing the number of taxa and loci increases amounts of missing data. Thus, a fundamental (but still unresolved) question is whether sampling should be designed to maximize sampling of taxa or genes, or to minimize the inclusion of missing data cells. Here, we explore this question for an ancient, rapid radiation of lizards, the pleurodont iguanians. Pleurodonts include many well-known clades (e.g., anoles, basilisks, iguanas, and spiny lizards) but relationships among families have proven difficult to resolve strongly and consistently using traditional sequencing approaches. We generated up to 4921 ultraconserved elements with sampling strategies including 16, 29, and 44 taxa, from 1179 to approximately 2.4 million characters per matrix and approximately 30% to 60% total missing data. We then compared mean branch support for interfamilial relationships under these 15 different sampling strategies for both concatenated (maximum likelihood) and species tree (NJst) approaches (after showing that mean branch support appears to be related to accuracy). We found that both approaches had the highest support when including loci with up to 50% missing taxa (matrices with ~40-55% missing data overall). Thus, our results show that simply excluding all missing data may be highly problematic as the primary guiding principle for the inclusion or exclusion of taxa and genes. The optimal strategy was somewhat different for each approach, a pattern that has not been shown previously. For concatenated analyses, branch support was maximized when including many taxa (44) but fewer characters (1.1 million). For species-tree analyses, branch support was maximized with minimal taxon sampling (16) but many loci (4789 of 4921). We also show that the choice of these sampling strategies can be critically important for phylogenomic analyses, since some strategies lead to demonstrably incorrect inferences (using the same method) that have strong statistical support. Our preferred estimate provides strong support for most interfamilial relationships in this important but phylogenetically challenging group.

Molecular phylogeny of microhylid frogs (Anura: Microhylidae) with emphasis on relationships among New World genera

Young species complexes that are widespread across ecologically disparate regions offer important insights into the process of speciation because of their relevance to how local adaptation and gene flow influence diversification. We used mitochondrial DNA and up to 28 152 genomewide single nucleotide polymorphisms from polytypic barking frogs (Craugastor augusti complex) to infer phylogenetic relationships and test for the signature of introgressive hybridization among diverging lineages. Our phylogenetic reconstructions suggest (i) a rapid Pliocene–Pleistocene radiation that produced at least nine distinct lineages and (ii) that geographic features of the arid Central Mexican Plateau contributed to two independent northward expansions. Despite clear lineage differentiation (many private alleles and high between‐lineage FST scores), D‐statistic tests, which differentiate introgression from ancestral polymorphism, allowed us to identify two putative instances of reticulate gene flow. Partitioned D‐statistics provided evidence that these events occurred in the same direction between clades but at different points in time. After correcting for geographic distance, we found that lineages involved in hybrid gene flow interactions had higher levels of genetic variation than independently evolving lineages. These findings suggest that the nature of hybrid compatibility can be conserved overlong periods of evolutionary time and that hybridization between diverging lineages may contribute to standing levels of genetic variation.

Multilocus molecular phylogenetic analysis of the montane Craugastor podiciferus species complex (Anura: Craugastoridae) in Isthmian Central America

BackgroundOver the last ten years we have seen great efforts focused on revising amphibian systematics. Phylogenetic reconstructions derived from DNA sequence data have played a central role in these revisionary studies but have typically under-sampled the diverse frog family Microhylidae. Here, we present a detailed phylogenetic study focused on expanding previous hypotheses of relationships within this cosmopolitan family. Specifically, we placed an emphasis on assessing relationships among New World genera and those taxa with uncertain phylogenetic affinities (i.e., incertae sedis).ResultsOne mitochondrial and three nuclear genes (about 2.8 kb) were sequenced to assess phylogenetic relationships. We utilized an unprecedented sampling of 200 microhylid taxa representing 91% of currently recognized subfamilies and 95% of New World genera. Our analyses do not fully resolve relationships among subfamilies supporting previous studies that have suggested a rapid early diversification of this clade. We observed a close relationship between Synapturanus and Otophryne of the subfamily Otophryninae. Within the subfamily Gastrophryninae relationships between genera were well resolved.ConclusionOtophryninae is distantly related to all other New World microhylids that were recovered as a monophyletic group, Gastrophryninae. Within Gastrophryninae, five genera were recovered as non-monophyletic; we propose taxonomic re-arrangements to render all genera monophyletic. This hypothesis of relationships and updated classification for New World microhylids may serve as a guide to better understand the evolutionary history of this group that is apparently subject to convergent morphological evolution and chromosome reduction. Based on a divergence analysis calibrated with hypotheses from previous studies and fossil data, it appears that microhylid genera inhabiting the New World originated during a period of gradual cooling from the late Oligocene to mid Miocene.

Phylogenomic analyses reveal novel relationships among snake families

The Craugastor podiciferus complex is a group of phenotypically polymorphic direct-developing frogs that inhabit the Talamancan highlands of Costa Rica and Panama. The montane distribution of this group creates natural allopatry among members and offers an excellent opportunity to explore geographic models of speciation. Using a multilocus approach, we obtained data from one nuclear (c-myc) and three mitochondrial (12S, 16S, and COI) gene regions from 40 individuals within the C. podiciferus complex. Molecular phylogenetic analyses revealed a basal split that placed samples from western Panama as sister to Costa Rican (CR) samples, corroborating a previous suggestion that the former lineage may represent an undescribed species. Within the CR clades we found six distinct haplogroups whose distributions largely corresponded to geographic features and included instances of sympatry. Divergence estimates were used to develop a preliminary evolutionary timeframe for the diversification of the C. podiciferus complex. Based on collective evidence, we hypothesize that movement of the CR haplogroups has occurred between currently isolated areas of suitable habitat via second order climatic fluctuations during the Pleistocene. The levels of genetic differentiation within the C. podiciferus complex are remarkable given the relatively small geographic area (ca. 8000 km(2)) of occurrence. This diversity emphasizes the need for further study and taxonomic revision to aid in conservation planning for this complex which, like many amphibians, has experienced recent population declines.

Thousands of microsatellite loci from the venomous coralsnake Micrurus fulvius and variability of select loci across populations and related species

Snakes are a diverse and important group of vertebrates. However, relationships among the major groups of snakes have remained highly uncertain, with recent studies hypothesizing very different (and typically weakly supported) relationships. Here, we address family-level snake relationships with new phylogenomic data from 3776 nuclear loci from ultraconserved elements (1.40million aligned base pairs, 52% missing data overall) sampled from 29 snake species that together represent almost all families, a dataset ∼100 times larger than used in previous studies. We found relatively strong support from species-tree analyses (NJst) for most relationships, including three largely novel clades: (1) a clade uniting the boas, pythons and their relatives, (2) a clade placing cylindrophiids and uropeltids with this clade, and (3) a clade uniting bolyeriids (Round Island boas) with pythonids and their relatives (xenopeltids and loxocemids). Relationships among families of advanced snakes (caenophidians) were also strongly supported. The results show the potential for phylogenomic analyses to resolve difficult groups, but also show a surprising sensitivity of the analyses to the inclusion or exclusion of outgroups.

Evolutionary relationships amongst polymorphic direct-developing frogs in the Craugastor rhodopis Species Group (Anura: Craugastoridae)

Studies of population genetics increasingly use next‐generation DNA sequencing to identify microsatellite loci in nonmodel organisms. There are, however, relatively few studies that validate the feasibility of transitioning from marker development to experimental application across populations and species. North American coralsnakes of the Micrurus fulvius species complex occur in the United States and Mexico, and little is known about their population structure and phylogenetic relationships. This absence of information and population genetics markers is particularly concerning because they are highly venomous and have important implications on human health. To alleviate this problem in coralsnakes, we investigated the feasibility of using 454 shotgun sequences for microsatellite marker development. First, a genomic shotgun library from a single individual was sequenced (approximately 7.74 megabases; 26 831 reads) to identify potentially amplifiable microsatellite loci (PALs). We then hierarchically sampled 76 individuals from throughout the geographic distribution of the species complex and examined whether PALs were amplifiable and polymorphic. Approximately half of the loci tested were readily amplifiable from all individuals, and 80% of the loci tested for variation were variable and thus informative as population genetic markers. To evaluate the repetitive landscape characteristics across multiple snakes, we also compared microsatellite content between the coralsnake and two other previously sampled snakes, the venomous copperhead (Agkistrodon contortrix) and Burmese python (Python molurus).

Rapid range expansion in the Great Plains narrow-mouthed toad (Gastrophryne olivacea) and a revised taxonomy for North American microhylids

The Craugastor rhodopis Species Group includes two leaf-litter frog species (C. loki and C. rhodopis). These direct-developing frogs inhabit tropical regions of Mexico and northern Central America. Characterizing diversity within the group has been difficult due to high levels of phenotypic polymorphism within and between species. Because of these polymorphisms, each taxon has junior synonyms. Using a fragment of mitochondrial DNA (mtDNA), we investigated genetic diversity in the C. rhodopis Species Group. We then examined type specimens (including types of junior synonyms) to match nomenclature to geographically circumscribed genetic clusters. Our molecular analyses revealed four major lineages within the C. rhodopis Species Group: (1) a widely distributed clade in western Mexico, (2) a highland clade in eastern Mexico, (3) a widely distributed lowland clade occurring in eastern Mexico, Guatemala and El Salvador, and (4) a haplotype from Volcán San Martín in Veracruz, Mexico. We identified the first clade as C. occidentalis, a taxon currently placed in the ecologically similar but phylogenetically distant C. mexicanus Species Series. In light of this we place C. occidentalis in the C. rhodopis Species Group and designate a lectotype and paralectotype for the species. The second and third clades inhabiting eastern Mexico and northern Central America correspond to C. rhodopis and C. loki, respectively. Additionally, we examined the taxonomic distribution of certain colour pattern traits and compensatory mutations in Domain III of the mtDNA 12S ribosomal RNA gene. Our recovery of the divergent Veracruz haplotype and extensive mtDNA structure within species indicates that additional taxonomic revision will be necessary.

Genetic surfing, not allopatric divergence, explains spatial sorting of mitochondrial haplotypes in venomous coralsnakes.

We investigated genetic variation within the Great Plains narrow-mouthed toad, Gastrophryne olivacea, across its geographic range in the United States and Mexico. An analysis of mitochondrial DNA (mtDNA) from 105 frogs revealed remarkably low levels of genetic diversity in individuals inhabiting the central United States and northern Mexico. We found that this widespread matrilineal lineage is divergent (ca. 2% in mtDNA) from haplotypes that originate from the western United States and western coast of Mexico. Using a dataset that included all five species of Gastrophryne and both species of the closely related genus Hypopachus, we investigated the phylogenetic placement of G. olivacea. This analysis recovered strong support that G. olivacea, the tropically distributed G. elegans, and the temperately distributed G. carolinensis constitute a monophyletic assemblage. However, the placement of G. pictiventris and G. usta render Gastrophryne paraphyletic with respect to Hypopachus. To complement our mitochondrial analysis, we examined a small fragment of nuclear DNA and recovered consistent patterns. In light of our findings we recommend (1) the resurrection of the nomen G. mazatlanensisTaylor (1943) for the disjunct western clade of G. olivacea and (2) the tentative placement of G. pictiventris and G. usta in Hypopachus. To explore possible scenarios leading to low levels of genetic diversity in G. olivacea, we used mismatch distributions and Bayesian Skyline plots to examine historic population expansion and demography. Collectively these analyses suggest that G. olivacea rapidly expanded in effective population size and geographic range during the late Pleistocene or early Holocene. This hypothesis is consistent with fossil data from northern localities and contemporary observations that suggest ongoing northern expansion. Given our findings, we suspect that the rapid range expansion of G. olivacea may have been facilitated by ecological associations with open habitats and seasonal water bodies.