Phylogeography and conservation genetics of the hellbender salamander (Cryptobranchus alleganiensis)

We investigated hellbender phylogeography through phylogenetic analyses of individuals sampled from 16 locations throughout their range in the eastern United States. Analyses were conducted on concatenated cytochrome-oxidase I (COI), cytochrome-b (Cytb) and NADH dehydrogenase subunit 4 (ND4) mtDNA sequence, totaling 2160 nucleotides. Hellbender haplotypes differed by 0.1 to 5.8% maximum likelihood (ML) corrected sequence divergence. Phylogenetic analyses reveal that hellbenders are separated into 8 reciprocally monophyletic populations or clades differentiated by a minimum of 0.7 to 5.4% sequence divergence, each of which constitutes a separate Management Unit (MU). High among population divergence and reciprocal monophyly suggest that female-mediated gene flow is severely restricted or non-existent among each MU. Hellbenders are currently divided into two subspecies, Cryptobranchus alleganiensis alleganiensis and C. a. bishopi based on morphological characters. The phylogenetic analyses presented here strongly indicate that these subspecies are paraphyletic. Management priorities for the hellbender should be reconsidered in light of these new molecular data. Results from Bayesian rooting indicate the root of the hellbender mtDNA tree lies on the branch leading to hellbender haplotypes from the Current, Eleven Point and New Rivers. The rooted tree suggests that a common ancestor in the southern Ozarks and/or southern Appalachians gave rise to northern hellbender populations, consistent with a Pleistocene refuge hypothesis.     

Selection of evolutionary models for phylogenetic hypothesis testing using parametric methods

Recent molecular studies have incorporated the parametric bootstrap method to test a priori hypotheses when the results of molecular based phylogenies are in conflict with these hypotheses. The parametric bootstrap requires the specification of a particular substitutional model, the parameters of which will be used to generate simulated, replicate DNA sequence data sets. It has been both suggested that, (a) the method appears robust to changes in the model of evolution, and alternatively that, (b) as realistic model of DNA substitution as possible should be used to avoid false rejection of a null hypothesis. Here we empirically evaluate the effect of suboptimal substitution models when testing hypotheses of monophyly with the parametric bootstrap using data sets of mtDNA cytochrome oxidase I and II (COI and COII) sequences for Macaronesian Calathus beetles, and mitochondrial 16S rDNA and nuclear ITS2 sequences for European Timarcha beetles. Whether a particular hypothesis of monophyly is rejected or accepted appears to be highly dependent on whether the nucleotide substitution model being used is optimal. It appears that a parameter rich model is either equally or less likely to reject a hypothesis of monophyly where the optimal model is unknown. A comparison of the performance of the Kishino–Hasegawa (KH) test shows it is not as severely affected by the use of suboptimal models, and overall it appears to be a less conservative method with a higher rate of failure to reject null hypotheses.     

The endangered Sonoran topminnow: Examination of species and ESUs using three mtDNA genes

There has been controversy over the species status of Sonoran topminnows and debate about the presence of ESUs in the Gila topminnow. From examination of sequence variation at 2626 base pairs over three mtDNA genes, we found a 29 (1.1%) nucleotide genetic difference between Gila and Yaqui topminnows. This provides strong support that these two taxa are separate species, Poeciliopsis occidentalis (Gila topminnow) and P. sonoriensis (Yaqui topminnow) and have been separated for approximately one million years. All the Gila topminnows within Arizona have the same sequence for the three mtDNA genes, that is, there is not reciprocal monophyly for mtDNA sequence data for the two previously designated ESUs. However, evidence of the unique habitat for Monkey Spring, its long-term isolation from other Gila topminnow habitats, and the presence of unique fish and invertebrate taxa in Monkey Spring support the designation of the Monkey Spring topminnows as an ESU. Finally, theoretical considerations using molecular data and estimates of heterozygosity and genetic distance for nuclear genes between populations of the Gila topminnow show that the lack of mtDNA variation is not inconsistent with the level and pattern of nuclear genetic variation observed.     

Base Compositional Bias and Phylogenetic Analyses: A Test of the “Flying DNA” Hypothesis

Phylogenetic methods can produce biased estimates of phylogeny when base composition varies along different lineages. Pettigrew (1994,Curr. Biol.4:277–280) has suggested that base composition bias is responsible for the apparent support for the monophyly of bats (Chiroptera: megabats and microbats) from several different nuclear and mitochondrial genes. Pettigrew's “flying DNA” hypothesis makes several predictions: (1) that metabolic constraints associated with flying result in elevated levels of adenine and thymine throughout the genome of both megabats and microbats, (2) that the resulting base compositional bias in bats is sufficient to mislead phylogenetic methods and account for the support for bat monophyly from several nuclear and mitochondrial genes, and (3) that phylogenetic analysis using pairwise distances corrected for compositional bias should eliminate the support for bat monophyly. We tested these predictions by analyzing DNA sequences from two nuclear and three mitochondrial genes. The predicted base compositional bias does not appear to exist in some of the genes, and in other genes the differences in AT content are very small. Analyses under a wide diversity of criteria and models of evolution, including analyses that take base composition into account (using log-determinant distances), all strongly support bat monophyly. Moreover, simulation analyses indicate that even extreme bias toward AT-base composition in bats would be insufficient to explain the observed levels of support for bat monophyly. These analyses provide no support for the “flying DNA” hypothesis, whereas the monophyly of bats appears to be well supported by the DNA sequence data.     

Finding stories in noise: Mitochondrial portraits from RAD data

Mitochondrial DNA (mtDNA) has formed the backbone of phylogeographic research for many years; however, recent trends focus on genome‐wide analyses. One method proposed for calibrating inferences from noisy next‐generation data, such as RAD sequencing, is to compare these results with analyses of mitochondrial sequences. Most researchers using this approach appear to be unaware that many single nucleotide polymorphisms (SNPs) identified from genome‐wide sequence data are themselves mitochondrial, or assume that these are too few to bias analyses. Here, we demonstrate two methods for mining mitochondrial markers using RAD sequence data from three South African species of yellowfish, Labeobarbus. First, we use a rigorous SNP discovery pipeline using the program stacks , to identify variant sites in mtDNA, which we then combine into haplotypes. Second, we directly map sequence reads against a mitochondrial genome reference. This method allowed us to reconstruct up to 98% of the Labeobarbus mitogenome. We validated these mitogenome reconstructions through blast database searches and by comparison with cytochrome b gene sequences obtained through Sanger sequencing. Finally, we investigate the organismal consequences of these data including ancient genetic exchange and a recent translocation among populations of L. natalensis, as well as interspecific hybridization between L. aeneus and L. kimberleyensis.     

Higher‐level phylogeny of the insect order Hemiptera: is Auchenorrhyncha really paraphyletic?

The higher‐level phylogeny of the order Hemiptera remains a contentious topic in insect systematics. The controversy is chiefly centred on the unresolved question of whether or not the hemipteran suborder Auchenorrhyncha (including the extant superfamilies Fulgoroidea, Membracoidea, Cicadoidea and Cercopoidea) is a monophyletic lineage. Presented here are the results of a multilocus molecular phylogenetic investigation of relationships among the major hemipteran lineages, designed specifically to address the question of Auchenorrhyncha monophyly in the context of broad taxonomic sampling across Hemiptera. Phylogenetic analyses (maximum parsimony, maximum likelihood and Bayesian inference) were based on DNA nucleotide sequence data from seven gene regions (18S rDNA, 28S rDNA, histone H3, histone 2A, wingless, cytochrome c oxidase I and NADH dehydrogenase subunit 4) generated from 86 in‐group exemplars representing all major lineages of Hemiptera (plus seven out‐group taxa). All combined analyses of these data recover the monophyly of Auchenorrhyncha, and also support the monophyly of each of the following lineages: Hemiptera, Sternorrhyncha, Heteropterodea, Heteroptera, Fulgoroidea, Cicadomorpha, Membracoidea, Cercopoidea and Cicadoidea. Also presented is a review of the major lines of morphological and molecular evidence for and against the monophyly of Auchenorrhyncha.     

Population subdivision among desert bighorn sheep (Oviscanadensis) ewes revealed by mitochondrial DNA analysis

We used behavioural observations and mitochondrial DNA (mtDNA) sequence analysis to examine demographic and genetic structure within and among home-range groups of desert bighorn sheep (Oviscanadensis) ewes in the Peninsular Ranges of southern California, USA. We identified substantial genetic variation in the first 515 bp of the mtDNA control region and determined that seven haplotypes were distributed in a nonrandom fashion among these ewe subpopulations. Although a significant (P < 0.01) amount of mtDNA variation (33%) was partitioned among home-range groups, we did not find strong evidence for matrilineal substructuring within these groups. Based on analyses of molecular variance, and comparisons of behavioural associations and distances between centres of activity, we concluded that within a given home-range group, bighorn sheep ewes generally associate with other ewes based on their availability rather than their matrilineal relationships. Our results also supported the conclusion that multiple ewe subpopulations exist within the Peninsular Ranges, and that these subpopulations are the most basic demographic and genetic units.     

Incongruent gene trees, complex evolutionary processes, and the phylogeny of a group of North American minnows (Hybognathus Agassiz 1855)

Hybognathus is a putatively monophyletic group of North American minnows containing seven extant species. Although much is known about the taxonomy, biology, and life history of Hybognathus species, their phylogenetic relations with each other remain unclear. We address this problem with partial sequences of mitochondrial cytochrome-b, 16S rRNA, and ND4 mtDNA genes and nuclear growth hormone (GH) and S7 introns from representatives of all Hybognathus species and four outgroup taxa. Phylogenetic analyses of these data corroborated previous studies on the monophyly of seven recognized species of Hybognathus, and indicated weak support for monophyly of the genus. Topological tests, however, revealed significant (all P < 0.001) conflict among molecular data sets. Ad hoc removal of taxa from topologies and subsequent testing indicated that incongruence was localized to two specific ingroup taxa (H. hayi and H. regius) and suggested that the conflict is a function of underlying processes that have generated the observed phylogenetic patterns. Hybridization (ancestral), which is commonplace in cyprinids, may best explain topological disagreements among datasets; however, retention of ancestral polymorphism and natural selection remain as alternative hypotheses. Our results highlight methodological and topological problems associated with estimating interspecific phylogenies from multiple genes.     

Phylogeny of iguanian lizards inferred from 29 nuclear loci, and a comparison of concatenated and species-tree approaches for an ancient, rapid radiation

Iguanian lizards form a diverse clade whose members have been the focus of many comparative studies of ecology, behavior, and evolution. Despite the importance of phylogeny to such studies, interrelationships among many iguanian clades remain uncertain. Within the Old World clade Acrodonta, Agamidae is sometimes found to be paraphyletic with respect to Chamaeleonidae, and recent molecular studies have produced conflicting results for many major clades. Within the largely New World clade Pleurodonta, relationships among the 12 currently recognized major subclades (mostly ranked as families) have been largely unresolved or poorly supported in previous studies. To clarify iguanian evolutionary history, we first infer phylogenies using concatenated maximum-likelihood (ML) and Bayesian analyses of DNA sequence data from 29 nuclear protein-coding genes for 47 iguanian and 29 outgroup taxa. We then estimate a relaxed-clock Bayesian chronogram for iguanians using BEAST. All three methods produce identical topologies. Within Acrodonta, we find strong support for monophyly of Agamidae with respect to Chamaeleonidae, and for almost all relationships within agamids. Within Pleurodonta, we find strong Bayesian support for almost all relationships, and strong ML support for some interfamilial relationships and for monophyly of almost all families (excepting Polychrotidae). Our phylogenetic results suggest a non-traditional biogeographic scenario in which pleurodonts originated in the Northern Hemisphere and subsequently spread southward into South America. The pleurodont portion of the tree is characterized by several very short, deep branches, raising the possibility of deep coalescences that may confound concatenated analyses. We therefore also use 27 of these genes to implement a coalescent-based species-tree approach for pleurodonts. Although this analysis strongly supports monophyly of the pleurodont families, interfamilial relationships are generally different from those in the concatenated tree, and support is uniformly poor. However, a species-tree analysis using only the seven most variable loci yields higher support and more congruence with the concatenated tree. This suggests that low support in the 27-gene species-tree analysis may be an artifact of the many loci that are uninformative for very short branches. This may be a general problem for the application of species-tree methods to rapid radiations, even with phylogenomic data sets. Finally, we correct the non-monophyly of Polychrotidae by recognizing the pleurodont genus Anolis (sensu lato) as a separate family (Dactyloidae), and we correct the non-monophyly of the agamid genus Physignathus by resurrection of the genus Istiurus for the former Physignathus lesueurii.     

Molecular analysis off the phylogeny off 11 genera off the Corvidae

Phylogenetic analysis of 11 genera of Corvidae has been undertaken using a sequence of 925 bp of the cytochrome b gene of the mtDNA. Both maximum-parsimony and maximum-likelihood analyses were performed, and different weighting schemes have been used according to the saturation observed in the data (mostly at the third codon position). The cytochrome b data provide good evidence for the monophyly of the Corvidae, including the Piapiac Ptilostomus afer in spite of its aberrant morphology. In all analyses, the choughs are the first lineage emerging in the family; this result conflicts with the traditional chough-crow group, as well as the hypothesis of the ancestral position of the jays. The monophyly of the American jays is well supported by the molecular data, unlike the Palearctic jays, which are not closely related. The other intergeneric relationships are not strongly supported by bootstrap proportions, and are not congruent in all different weightings and methods, except for groups of species of the same or very closely related genera.