Phylogeographic structure of the dunes sagebrush lizard, an endemic habitat specialist.

Phylogeographic divergence and population genetic diversity within species reflect the impacts of habitat connectivity, demographics, and landscape level processes in both the recent and distant past. Characterizing patterns of differentiation across the geographic range of a species provides insight on the roles of organismal and environmental traits in evolutionary divergence and future population persistence. This is particularly true of habitat specialists where habitat availability and resource dependence may result in pronounced genetic structure as well as increased population vulnerability. We use DNA sequence data as well as microsatellite genotypes to estimate range-wide phylogeographic divergence, historical population connectivity, and historical demographics in an endemic habitat specialist, the dunes sagebrush lizard (Sceloporus arenicolus). This species is found exclusively in dune blowouts and patches of open sand within the shinnery oak-sand dune ecosystem of southeastern New Mexico and adjacent Texas. We find evidence of phylogeographic structure consistent with breaks and constrictions in suitable habitat at the range-wide scale. In addition, we find support for a dynamic and variable evolutionary history across the range of S. arenicolus. Populations in the Monahans Sandhills have deeply divergent lineages consistent with long-term demographic stability. In contrast, populations in the Mescalero Sands are not highly differentiated, though we do find evidence of demographic expansion in some regions and relative demographic stability in others. Phylogeographic history and population genetic differentiation in this species has been shaped by the configuration of habitat patches within a geologically complex and historically dynamic landscape. Our findings identify regions as genetically distinctive conservation units as well as underscore the genetic and demographic history of different lineages of S. arenicolus.

Fire Does Not Strongly Affect Genetic Diversity or Structure of a Common Treefrog in the Endangered Florida Scrub.

Fire regimes influence natural populations of organisms in diverse ways, via direct effects on population dynamics as well as indirect effects on habitat and ecosystem processes. Although many amphibian species have evolved to persist in fire-dependent ecosystems, the effects of fire on the genetic diversity of amphibian populations remain relatively unexplored. We examined how different aspects of fire history relate to population genetic diversity and structure of an abundant anuran, Hyla femoralis, in a large, intact area of Florida scrub containing hundreds of seasonally inundated ponds. Specifically, we assessed the overall population genetic structure and examined whether variation in time since fire, fire intensity, or historical fire frequency at breeding sites explained spatial variation in genetic diversity. Based on our sampling of 17 breeding aggregations within the 2,100-ha study area, neither recent nor frequent fire reduce genetic diversity or restrict connectivity among ponds for H. femoralis. Overall, mean effective population sizes were large (average range = 68-572). We detected a positive trend between effective population size (Ne) and average intensity of the most-recent fire, with this factor explaining 42% of the variation in Ne. Our results contrast with previous studies that consistently demonstrate strong relationships between fire history and population genetic structure of scrub-associated lizard species, suggesting that H. femoralis is resilient to a wide range of fire regimes. More generally, our study contributes to understanding the roles of life-history characteristics and environmental unpredictability in shaping organisms' responses to fire.

Phylogeography of the arid-adapted Malagasy bullfrog, Laliostoma labrosum, influenced by past connectivity and habitat stability.

The rainforest biome of eastern Madagascar is renowned for its extraordinary biodiversity and restricted distribution ranges of many species, whereas the arid western region of the island is relatively species poor. We provide insight into the biogeography of western Madagascar by analyzing a multilocus phylogeographic dataset assembled for an amphibian, the widespread Malagasy bullfrog, Laliostoma labrosum. We find no cryptic species in L. labrosum (maximum 1.1% pairwise genetic distance between individuals in the 16S rRNA gene) attributable to considerable gene flow at the regional level as shown by genetic admixture in both mtDNA and three nuclear loci, especially in central Madagascar. Low breeding site fidelity, viewed as an adaptation to the unreliability of standing pools of freshwater in dry and seasonal environments, and a ubiquitous distribution within its range may underlie overall low genetic differentiation. Moreover, reductions in population size associated with periods of high aridity in western Madagascar may have purged DNA variation in this species. The mtDNA gene tree revealed seven major phylogroups within this species, five of which show mostly non-overlapping distributions. The nested positions of the northern and central mtDNA phylogroups imply a southwestern origin for all extant mtDNA lineages in L. labrosum. The current phylogeography of this species and paleo-distributions of major mtDNA lineages suggest five potential refugia in northern, western and southwestern Madagascar, likely the result of Pleistocene range fragmentation during drier and cooler climates. Lineage sorting in mtDNA and nuclear loci highlighted a main phylogeographic break between populations north and south of the Sambirano region, suggesting a role of the coastal Sambirano rainforest as a barrier to gene flow. Paleo-species distribution models and dispersal networks suggest that the persistence of some refugial populations was mainly determined by high population connectivity through space and time.

Molecular evolutionary characterization of a V1R subfamily unique to strepsirrhine primates.

Vomeronasal receptor genes have frequently been invoked as integral to the establishment and maintenance of species boundaries among mammals due to the elaborate one-to-one correspondence between semiochemical signals and neuronal sensory inputs. Here, we report the most extensive sample of vomeronasal receptor class 1 (V1R) sequences ever generated for a diverse yet phylogenetically coherent group of mammals, the tooth-combed primates (suborder Strepsirrhini). Phylogenetic analysis confirms our intensive sampling from a single V1R subfamily, apparently unique to the strepsirrhine primates. We designate this subfamily as V1Rstrep. The subfamily retains extensive repertoires of gene copies that descend from an ancestral gene duplication that appears to have occurred prior to the diversification of all lemuriform primates excluding the basal genus Daubentonia (the aye-aye). We refer to the descendent clades as V1Rstrep-α and V1Rstrep-ß. Comparison of the two clades reveals different amino acid compositions corresponding to the predicted ligand-binding site and thus potentially to altered functional profiles between the two. In agreement with previous studies of the mouse lemur (genus, Microcebus), the majority of V1Rstrep gene copies appear to be intact and under strong positive selection, particularly within transmembrane regions. Finally, despite the surprisingly high number of gene copies identified in this study, it is nonetheless probable that V1R diversity remains underestimated in these nonmodel primates and that complete characterization will be limited until high-coverage assembled genomes are available.

Review of the systematic status of Sceloporus arenicolus Degenhardt and Jones, 1972 with an estimate of divergence time.

The sagebrush lizards (Sceloporus graciosus group) consist of four taxa (S. graciosus graciosus, S. graciosus gracilis, S. graciosus vandenburgianus, and S. arenicolus) distributed in western North America. Of these, S. arenicolus is morphologically, behaviorally, and ecologically distinct as well as geographically disjunct from the other taxa, occurring only in the Mescalero-Monahans Sandhills of southeastern New Mexico and adjacent Texas. Sceloporus arenicolus is a taxon of concern because of its small range and habitat alteration due to land use practices. Understanding evolutionary relationships among members of the S. graciosus group, and especially S. arenicolus, has important implications for conservation. We examine the phylogenetic relationship of S. arenicolus relative to the three recognized subspecies of S. graciosus at mitochondrial and nuclear loci for populations sampled throughout the ranges of these taxa. Additionally, we estimate the divergence time and clade age of S. arenicolus. We find that the S. graciosus group is in need of major taxonomic revision, and also confirm that S. arenicolus is a genetically distinct and divergent lineage. These results bear important consequences for conservation and management.

Defining spatial and temporal patterns of phylogeographic structure in Madagascar's iguanid lizards (genus Oplurus).

Understanding the remarkably high species diversity and levels of endemism found among Madagascar's flora and fauna has been the focus of many studies. One hypothesis that has received much attention proposes that Quaternary climate fluctuations spurred diversification. While spatial patterns of distribution and phylogenetic relationships can provide support for biogeographic predictions, temporal estimates of divergence are required to determine the fit of these geospatial patterns to climatic or biogeographic mechanisms. We use multilocus DNA sequence data to test whether divergence times among Malagasy iguanid lizards of the subfamily Oplurinae are compatible with a hypotheses of Pliocene-Pleistocene diversification. We estimate the oplurine species tree and associated divergence times under a relaxed-clock model. In addition, we examine the phylogeographic structure and population divergence times within two sister species of Oplurus primarily distributed in the north-west and south-west of Madagascar (Oplurus cuvieri and Oplurus cyclurus, respectively). We find that divergence events among oplurine lineages occurred in the Oligocene and Miocene and are thus far older and incompatible with the hypothesis that recent climate fluctuations are related to current species diversity. However, the timing of intraspecific divergences and spatial patterns of population genetic structure within O. cuvieri and O. cyclurus suggest a role for both intrinsic barriers and recent climate fluctuations at population-level divergences. Integrating information across spatial and temporal scales allows us to identify and better understand the mechanisms generating patterns diversity.

Concatenation and concordance in the reconstruction of mouse lemur phylogeny: an empirical demonstration of the effect of allele sampling in phylogenetics.

The systematics and speciation literature is rich with discussion relating to the potential for gene tree/species tree discordance. Numerous mechanisms have been proposed to generate discordance, including differential selection, long-branch attraction, gene duplication, genetic introgression, and/or incomplete lineage sorting. For speciose clades in which divergence has occurred recently and rapidly, recovering the true species tree can be particularly problematic due to incomplete lineage sorting. Unfortunately, the availability of multilocus or "phylogenomic" data sets does not simply solve the problem, particularly when the data are analyzed with standard concatenation techniques. In our study, we conduct a phylogenetic study for a nearly complete species sample of the dwarf and mouse lemur clade, Cheirogaleidae. Mouse lemurs (genus, Microcebus) have been intensively studied over the past decade for reasons relating to their high level of cryptic species diversity, and although there has been emerging consensus regarding the evolutionary diversity contained within the genus, there is no agreement as to the inter-specific relationships within the group. We attempt to resolve cheirogaleid phylogeny, focusing especially on the mouse lemurs, by employing a large multilocus data set. We compare the results of Bayesian concordance methods with those of standard gene concatenation, finding that though concatenation yields the strongest results as measured by statistical support, these results are found to be highly misleading. By employing an approach where individual alleles are treated as operational taxonomic units, we show that phylogenetic results are substantially influenced by the selection of alleles in the concatenation process.

Increased population sampling confirms low genetic divergence among Pteropus (Chiroptera: Pteropodidae) fruit bats of Madagascar and other western Indian Ocean islands.

Fruit bats of the genus Pteropus occur throughout the Austral-Asian region west to islands off the eastern coast of Africa. Recent phylogenetic analyses of Pteropus from the western Indian Ocean found low sequence divergence and poor phylogenetic resolution among several morphologically defined species. We reexamine the phylogenetic relationships of these taxa by using multiple individuals per species. In addition, we estimate population genetic structure in two well-sampled taxa occurring on Madagascar and the Comoro Islands (P. rufus and P. seychellensis comorensis). Despite finding a similar pattern of low sequence divergence among species, increased sampling provides insight into the phylogeographic history of western Indian Ocean Pteropus, uncovering high levels of gene flow within species.

Integrating statistical genetic and geospatial methods brings new power to phylogeography.

The field of phylogeography continues to grow in terms of power and accessibility. Initially uniting population genetics and phylogenetics, it now spans disciplines as diverse as geology, statistics, climatology, ecology, physiology, and bioinformatics to name a few. One major and recent integration driving the field forward is between "statistical phylogeography" and Geographic Information Systems (GIS) (Knowles, 2009). Merging genetic and geospatial data, and their associated methodological toolkits, is helping to bring explicit hypothesis testing to the field of phylogeography. Hypotheses derived from one approach can be reciprocally tested with data derived from the other field and the synthesis of these data can help place demographic events in an historical and spatial context, guide genetic sampling, and point to areas for further investigation. Here, we present three practical examples of empirical analysis that integrate statistical genetic and GIS tools to construct and test phylogeographic hypotheses. Insights into the evolutionary mechanisms underlying recent divergences can benefit from simultaneously considering diverse types of information to iteratively test and reformulate hypotheses. Our goal is to provide the reader with an introduction to the variety of available tools and their potential application to typical questions in phylogeography with the hope that integrative methods will be more broadly and commonly applied to other biological systems and data sets.

Phylogeny and biogeography of western Indian Ocean Rousettus (Chiroptera: Pteropodidae).

We examined patterns of genetic variation in Rousettus madagascariensis from Madagascar and R. obliviosus from the Comoros (Grande Comore, Anjouan, and Mohéli). Genetic distances among individuals on the basis of 1,130 base pairs of the mitochondrial cytochrome b (Cytb) locus were estimated from specimens collected from 17 sites on Madagascar, 3 sites on Grande Comore, 3 sites on Anjouan, and 2 sites on Mohéli. We observed little variation in Madagascar and nearshore island samples (maximum 1.1%) and interisland Comoros samples (maximum 1.8%). In contrast, pairwise distances between different sampled sites on Madagascar and the Comoros varied from 8.5% to 13.2%. For 131 Malagasy animals, 69 unique haplotypes were recovered with 86 variable sites, and for 44 Comorian individuals, 17 unique haplotypes were found with 30 variable sites. No haplotype was shared between Madagascar and the Comoros, adding to previous morphological evidence that these 2 populations should be considered separate species. Cytb data showed that Rousettus populations of Madagascar (including nearshore islands) and the Comoros are respectively monophyletic and display no geographic structure in haplotype diversity, and that R. madagascariensis and R. obliviosus are strongly supported as sister to each other relative to other Rousettus species. Genotypic data from 6 microsatellite loci confirm lack of geographic structure in either of the 2 species. In pairwise tests of population differentiation, the only significant values were between samples from the Comoro Islands and Madagascar (including nearshore islands). Estimates of current and historical demographic parameters support population expansion in both the Comoros and Madagascar. These data suggest a more recent and rapid demographic expansion in Madagascar in comparison with greater population stability on the Comoros. On the basis of available evidence, open-water crossings approaching 300 km seem rarely traversed by Rousettus, and, if successful, can result in genetic isolation and subsequent differentiation.

Population differentiation of temperate amphibians in unpredictable environments.

Amphibians are a globally distributed and diverse lineage, but much of our current understanding of their population genetic structure comes from studies in mesic temperate habitats. We characterize the population genetic structure of two sympatric explosive breeding amphibians in the southwestern deserts of the United States: the Great Plains toad (Anaxyrus cognatus) and Couch's spadefoot toad (Scaphiopus couchii). For both species, we find limited genetic differentiation even between populations in adjacent valleys separated by dispersal barriers such as mountainous habitats. To understand how population genetic patterns in these two arid-adapted species compare to taxa in more mesic environments, we computed a standardized measure of population differentiation for A. cognatus, S. couchii, and for pond-breeding amphibians that inhabit mesic temperate environments. Our results indicate that the arid-adapted species have lower population genetic structure at fine and moderate scales than most other amphibian species we surveyed. We hypothesize that stochasticity in the availability of appropriate breeding sites as well as landscape homogeneity may result in increased population connectivity in desert-adapted frogs. Future work examining fine-scale population structure in amphibians from a diversity of habitats will test the generality of our findings. Intraspecific comparisons among localities with varied seasonality and habitats will be particularly useful for investigating the interaction between species-typical population dynamics and environmental characteristics as determinants of population connectivity in pond-breeding amphibians.

Nuclear and cpDNA sequences combined provide strong inference of higher phylogenetic relationships in the phlox family (Polemoniaceae).

Members of the phlox family (Polemoniaceae) serve as useful models for studying various evolutionary and biological processes. Despite its biological importance, no family-wide phylogenetic estimate based on multiple DNA regions with complete generic sampling is available. Here, we analyze one nuclear and five chloroplast DNA sequence regions (nuclear ITS, chloroplast matK, trnL intron plus trnL-trnF intergeneric spacer, and the trnS-trnG, trnD-trnT, and psbM-trnD intergenic spacers) using parsimony and Bayesian methods, as well as assessments of congruence and long branch attraction, to explore phylogenetic relationships among 84 ingroup species representing all currently recognized Polemoniaceae genera. Relationships inferred from the ITS and concatenated chloroplast regions are similar overall. A combined analysis provides strong support for the monophyly of Polemoniaceae and subfamilies Acanthogilioideae, Cobaeoideae, and Polemonioideae. Relationships among subfamilies, and thus for the precise root of Polemoniaceae, remain poorly supported. Within the largest subfamily, Polemonioideae, four clades corresponding to tribes Polemonieae, Phlocideae, Gilieae, and Loeselieae receive strong support. The monogeneric Polemonieae appears sister to Phlocideae. Relationships within Polemonieae, Phlocideae, and Gilieae are mostly consistent between analyses and data permutations. Many relationships within Loeselieae remain uncertain. Overall, inferred phylogenetic relationships support a higher-level classification for Polemoniaceae proposed in 2000.