THE INTERFACE BETWEEN PHYLOGENETICS AND POPULATION GENETICS: INVESTIGATING GENE TREES, SPECIES TREES, AND POPULATION DYNAMICS IN THE PHYLLOPHAGA FRATERNA SPECIES GROUP

Multilocus Phylogeography and Species Delimitation in the Cumberland Plateau Salamander, Plethodon kentucki: Incongruence among Data Sets and Methods

Species are a fundamental unit of biodiversity, yet can be challenging to delimit objectively. This is particularly true of species complexes characterized by high levels of population genetic structure, hybridization between genetic groups, isolation by distance, and limited phenotypic variation. Previous work on the Cumberland Plateau Salamander, Plethodon kentucki, suggested that it might constitute a species complex despite occupying a relatively small geographic range. To examine this hypothesis, we sampled 135 individuals from 43 populations, and used four mitochondrial loci and five nuclear loci (5693 base pairs) to quantify phylogeographic structure and probe for cryptic species diversity. Rates of evolution for each locus were inferred using the multidistribute package, and time calibrated gene trees and species trees were inferred using BEAST 2 and *BEAST 2, respectively. Because the parameter space relevant for species delimitation is large and complex, and all methods make simplifying assumptions that may lead them to fail, we conducted an array of analyses. Our assumption was that strongly supported species would be congruent across methods. Putative species were first delimited using a Bayesian implementation of the GMYC model (bGMYC), Geneland, and Brownie. We then validated these species using the genealogical sorting index and BPP. We found substantial phylogeographic diversity using mtDNA, including four divergent clades and an inferred common ancestor at 14.9 myr (95% HPD: 10.8-19.7 myr). By contrast, this diversity was not corroborated by nuclear sequence data, which exhibited low levels of variation and weak phylogeographic structure. Species trees estimated a far younger root than did the mtDNA data, closer to 1.0 myr old. Mutually exclusive putative species were identified by the different approaches. Possible causes of data set discordance, and the problem of species delimitation in complexes with high levels of population structure and introgressive hybridization, are discussed.

Misleading phylogenetic inferences based on single-exemplar sampling in the turtle genus Pseudemys

Reconstructing species trees for clades containing weakly delimited or incorrectly identified taxa is one of the most serious challenges facing systematists because building phylogenetic trees is generally predicated on correctly identifying species membership for the terminals in an analysis. A common practice, particularly in large-scale phylogenetic analyses, is to use single-exemplar sampling under the implicit assumption that the resulting phylogenetic trees will be poorly supported if the sampled taxa are not good species. We examine this fundamental assumption in the North American turtle genus Pseudemys, a group of common, widely distributed freshwater turtles whose species boundaries and phylogenetic relationships have challenged systematists for over half a century. We sequenced 10 nuclear and three mitochondrial genes from the nine currently recognized species and subspecies of Pseudemys using geographically-widespread sampling of each taxon, and analyzed the resulting 86-individual data set using population-genetic and phylogenetic methods. We found little or no evidence supporting the division of Pseudemys into its currently recognized species/subspecies. Rather, our data strongly suggest that the group has been oversplit and contains fewer species than currently recognized. Even so, when we conducted 100 replicated, single-exemplar phylogenetic analyses of these same nine taxa, most Bayesian trees were well resolved, had high posterior probabilities, and yet returned completely conflicting topologies. These analyses suggest that phylogenetic analyses based on single-exemplar sampling may recover trees that depend on the individuals that are sampled, rather than the underlying species tree that systematists assume they are estimating. Our results clearly indicate that final resolution of Pseudemys will require an integrated analysis of morphology and historical biogeographic data coupled with extensive geographic sampling and large amounts of molecular data, and we do not recommend taxonomic changes based on our analyses. If our 100-tree resampling experiments generalize to other taxa, they suggest that single-exemplar phylogenies should be interpreted with caution, particularly for groups where species are shallowly diverged or inadequately delimited.

An approach to identify putative hybrids in the 'coalescent stochasticity zone', as exemplified in the African plant genus Streptocarpus (Gesneriaceae)

The inference of phylogenetic relationships is often complicated by differing evolutionary histories of independently-inherited markers. The causes of the resulting gene tree incongruence can be challenging to identify, often relying on coalescent simulations dependent on unverifiable assumptions. We investigated alternative techniques using the South African rosulate species of Streptocarpus as a study group. Two independent gene trees - from the nuclear ITS region and from three concatenated plastid regions (trnL-F, rpl20-rps12 and trnC-D) - displayed widespread, strongly supported incongruence. We investigated the causes by detecting genetic exchange across morphological borders using morphological optimizations and genetic exchange across species boundaries using the genealogical sorting index. Incongruence between gene trees was associated with ancestral shifts in growth form (in four species) but not in pollination syndrome, suggesting introgression limited by reproductive barriers. Genealogical sorting index calculations showed polyphyly of two additional species, while individuals of all others were significantly associated. In one case the association was stronger according to the internal transcribed spacer data than according to the plastid data, which, given the smaller effective population size of the plastid, may also indicate introgression. These approaches offer alternative ways to identify potential hybridization events where incomplete lineage sorting cannot be rejected using simulations.

Detecting cryptic speciation in the widespread and morphologically conservative carpet chameleon (Furcifer lateralis) of Madagascar

Species delimitation within recently evolved groups can be challenging because species may be difficult to distinguish morphologically. Following the General Lineage Concept, we apply a multiple evidence approach to assess species limits within the carpet chameleon Furcifer lateralis, which is endemic to Madagascar and exported in large numbers for the pet trade. Cryptic speciation within F. lateralis was considered likely because this species (1) has a vast distribution, (2) occupies exceptionally diverse habitats and (3) exhibits subtle regional differences in morphology. Phylogenetic trees reconstructed using nuclear and mitochondrial genes recovered three well-supported clades corresponding with geography. Morphological results based on canonical variates analysis show that these clades exhibit subtle differences in head casque morphology. Ecological niche modelling results found that these phylogenetic groups also occupy unique environmental space and exhibit patterns of regional endemism typical of other endemic reptiles. Combined, our findings provide diverse yet consistent evidence for the existence of three species. Consequently, we elevate the subspecies F. lateralis major to species rank and name a new species distributed in northern and western Madagascar. Initial ecological divergence, associated with speciation of F. lateralis in humid eastern habitat, fits the Ecographic Constraint model for species diversification in Madagascar. By contrast, the second speciation event provides some support for the Riverine Barrier model, with the Mangoky River possibly causing initial isolation between species. These findings thus support two contrasting models of speciation within closely related species and demonstrate the utility of applying a combined-evidence approach for detecting cryptic speciation.

Cryptic diversity within the major trypanosomiasis vector Glossina fuscipes revealed by molecular markers

BACKGROUND

The tsetse fly Glossina fuscipes s.l. is responsible for the transmission of approximately 90% of cases of human African trypanosomiasis (HAT) or sleeping sickness. Three G. fuscipes subspecies have been described, primarily based upon subtle differences in the morphology of their genitalia. Here we describe a study conducted across the range of this important vector to determine whether molecular evidence generated from nuclear DNA (microsatellites and gene sequence information), mitochondrial DNA and symbiont DNA support the existence of these taxa as discrete taxonomic units.

PRINCIPAL FINDINGS

The nuclear ribosomal Internal transcribed spacer 1 (ITS1) provided support for the three subspecies. However nuclear and mitochondrial sequence data did not support the monophyly of the morphological subspecies G. f. fuscipes or G. f. quanzensis. Instead, the most strongly supported monophyletic group was comprised of flies sampled from Ethiopia. Maternally inherited loci (mtDNA and symbiont) also suggested monophyly of a group from Lake Victoria basin and Tanzania, but this group was not supported by nuclear loci, suggesting different histories of these markers. Microsatellite data confirmed strong structuring across the range of G. fuscipes s.l., and was useful for deriving the interrelationship of closely related populations.

CONCLUSION/SIGNIFICANCE

We propose that the morphological classification alone is not used to classify populations of G. fuscipes for control purposes. The Ethiopian population, which is scheduled to be the target of a sterile insect release (SIT) programme, was notably discrete. From a programmatic perspective this may be both positive, given that it may reflect limited migration into the area or negative if the high levels of differentiation are also reflected in reproductive isolation between this population and the flies to be used in the release programme.