Phylogeny, divergence times and species limits of spiny lizards (Sceloporus magister species group) in western North American deserts and Baja California

Gene Flow and Isolation in the Arid Nearctic Revealed by Genomic Analyses of Desert Spiny Lizards

The opposing forces of gene flow and isolation are two major processes shaping genetic diversity. Understanding how these vary across space and time is necessary to identify the environmental features that promote diversification. The detection of considerable geographic structure in taxa from the arid Nearctic has prompted research into the drivers of isolation in the region. Several geographic features have been proposed as barriers to gene flow, including the Colorado River, Western Continental Divide (WCD), and a hypothetical Mid-Peninsular Seaway in Baja California. However, recent studies suggest that the role of barriers in genetic differentiation may have been overestimated when compared to other mechanisms of divergence. In this study, we infer historical and spatial patterns of connectivity and isolation in Desert Spiny Lizards (Sceloporus magister) and Baja Spiny Lizards (Sceloporus zosteromus), which together form a species complex composed of parapatric lineages with wide distributions in arid western North America. Our analyses incorporate mitochondrial sequences, genomic-scale data, and past and present climatic data to evaluate the nature and strength of barriers to gene flow in the region. Our approach relies on estimates of migration under the multispecies coalescent to understand the history of lineage divergence in the face of gene flow. Results show that the S. magister complex is geographically structured, but we also detect instances of gene flow. The WCD is a strong barrier to gene flow, while the Colorado River is more permeable. Analyses yield conflicting results for the catalyst of differentiation of peninsular lineages in S. zosteromus. Our study shows how large-scale genomic data for thoroughly sampled species can shed new light on biogeography. Furthermore, our approach highlights the need for the combined analysis of multiple sources of evidence to adequately characterize the drivers of divergence.

Molecular and morphometric analyses reveal host-specific cryptic speciation in a mite species, Tetranychus neocaledonicus (Andre, 1933) (Acari: Tetranychidae)

Host- and habitat-induced morphological shape and size variations are common in phytophagous and parasitic taxa. Several integrated morphological and molecular techniques have been commonly used to understand host-induced morpho-cryptic species forms. Compared to other arthropods, cryptic speciation was more common in Acari. This study focused on the host-specific morphological cryptic shape and size variations of Tetranychus neocaledonicus, collected from moringa and cassava hosts. We used geometric morphometric analysis to uncover the shape and size of inter-and intra-spider mite populations, and discovered that host-specific shape and size variations existed in spider mites regardless of sex. Interestingly, there was no phylogenetic signal in spider mites, implying that the morpho-cryptic speciation of T. neocaledonicus is solely based on the host-induced selection. The molecular clock hypothesis was accepted in our CO1 and 18s rRNA phylogeny analyses, and spider mites collected from both hosts were genetically less diverse. We conclude that T. neocaledonicus exhibited morphologically detectable cryptic population diversity in each host but that these populations are evolutionarily young form. Apart from these host-induced variations, we also monitored the impact of the clearing agent (lactic acid) on the shape and size of T. neocaledonicus; from this study, we proved that the clearing agent significantly alters the taxonomically important morphological traits of spider mites irrespective of the mites' sex, as confirmed by multivariate statistical analysis. This is the first study report to investigated the host-induced morphological variations of spider mites and the impact of a clearing agent.

Diversity and conservation of terrestrial vertebrates (birds, mammals, and reptiles) of Sierra Cucapá, Mexicali, Baja California, Mexico

Knowledge about the biodiversity of Baja California has been obtained mainly from natural protected areas (NPAs), while some unprotected natural areas have been poorly studied. The Sierra Cucapá in the northeast of the peninsula is one example. The objectives of this study are 1) to integrate existing knowledge of bird, mammal, and reptile diversity in Cucapá from public databases, citizen science platforms, and information generated from fieldwork, 2) to identify the spatial distribution of records in the study area, 3) to compare the composition of vertebrate species of Cucapá with that of NPAs of northern part of the peninsula, and 4) to assess the biological conservation value of Cucapá. We obtained records of 150 species of native vertebrates (102 birds, 34 mammals, and 14 reptiles) of which 10 species of birds, four mammals, and seven reptiles are included in a risk extinction category. The different sources of information contributed in a complementary way to the species inventories. Large areas in western and northern Cucapá lack records. The total difference in species composition between Cucapá and nearby NPAs ranged between 58 and 69% for birds, 61 and 79% for mammals, and 69 and 87% for reptiles. The species richness of Cucapá, its particular species composition, the presence of species in risk extinction categories, and the number and size of unexplored areas indicate that this area represents an opportunity for biological conservation in the northern part of the Peninsula. This work provides compelling data for the protection of Cucapá.

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.

Allopatric divergence and secondary contact with gene flow: a recurring theme in rattlesnake speciation

The study of recently diverged lineages whose geographical ranges come into contact can provide insight into the early stages of speciation and the potential roles of reproductive isolation in generating and maintaining species. Such insight can also be important for understanding the strategies and challenges for delimiting species within recently diverged species complexes. Here, we use mitochondrial and nuclear genetic data to study population structure, gene flow and demographic history across a geographically widespread rattlesnake clade, the western rattlesnake species complex (Crotalus cerberus, Crotalus viridis, Crotalus oreganus and relatives), which contains multiple lineages with ranges that overlap geographically or contact one another. We find evidence that the evolutionary history of this group does not conform to a bifurcating tree model and that pervasive gene flow has broadly influenced patterns of present-day genetic diversity. Our results suggest that lineage diversity has been shaped largely by drift and divergent selection in isolation, followed by secondary contact, in which reproductive isolating mechanisms appear weak and insufficient to prevent introgression, even between anciently diverged lineages. The complexity of divergence and secondary contact with gene flow among lineages also provides new context for why delimiting species within this complex has been difficult and contentious historically.

Genomic divergence in allopatric Northern Cardinals of the North American warm deserts is linked to behavioral differentiation

Biogeographic barriers are considered important in initiating speciation through geographic isolation, but they rarely indiscriminately and completely reduce gene flow across entire communities. Explicitly demonstrating which factors are associated with gene-flow levels across barriers would help elucidate how speciation is initiated and isolation maintained. Here, we investigated the association of behavioral isolation on population differentiation in Northern Cardinals (Cardinalis cardinalis) distributed across the Cochise Filter Barrier, a region of transitional habitat which separates the Sonoran and Chihuahuan deserts of North America. Using genomewide markers, we modeled demographic history by fitting the data to isolation and isolation-with-migration models. The best-fit model indicated that desert populations diverged in the Pleistocene with low, historic, and asymmetric gene flow across the barrier. We then tested behavioral isolation using reciprocal call-broadcast experiments to compare song recognition between deserts, controlling for song dialect changes within deserts. We found that male Northern Cardinals in both deserts were most aggressive to local songs and failed to recognize across-barrier songs. A correlation of genomic differentiation and strong song discrimination is consistent with a model where speciation is initiated across a barrier and maintained by behavioral isolation.

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.

Phylogeography of endemic Xantus' hummingbird (Hylocharis xantusii) shows a different history of vicariance in the Baja California Peninsula

Studies of phylogeographic patterns provide insight into the processes driving lineage divergence in a particular region. To identify the processes that caused phylogeographic breaks, it is necessary to use historical information and a set of appropriate molecular data to explain current patterns. To understand the influence of geological or ecological processes on the phylogeography of the only species of hummingbird endemic to the Baja California Peninsula, Hylocharis xantusii, mitochondrial DNA sequences of three concatenated genes (Cyt-b, COI and ND2; 2297bp in total) in 100 individuals were analyzed. The spatial analyses of genetic variation showed phylogeographic structure consisting of a north, central and south regions. According to estimated divergence times, two vicariant events are supported, permanent separation of the peninsula and formation of the Gulf of California at 5mya and temporary isolation of the southern region at the Isthmus of La Paz at 3mya. The temporal frame of genetic differentiation of intraspecific haplotypes indicates that 90% of haplotypes diverged within the last 500,000years, with a population expansion 80,000years ago. Only four haplotypes diverged ∼2.2 my and occurred in the south (Hxan_36, 38 and 45), and north (Hxan_45 and 56) regions; only haplotype 45 is shared between south and north populations. These regions also have the most recent haplotypes from 12,500 to 16,200years ago, and together with high levels of genetic diversity, we suggest two refuge areas, the Northern and Southern regions. Our results indicate that the phylogeographic pattern first results from vicariance processes, then is followed by historical and recent climate fluctuations that influenced conditions on the peninsula, and it is also related to oases distribution. This study presents the first investigation of phylogeography of the peninsular' endemic Xantus' hummingbird.

Phylogenomics of a rapid radiation: is chromosomal evolution linked to increased diversification in north american spiny lizards (Genus Sceloporus)?

Background

Resolving the short phylogenetic branches that result from rapid evolutionary diversification often requires large numbers of loci. We collected targeted sequence capture data from 585 nuclear loci (541 ultraconserved elements and 44 protein-coding genes) to estimate the phylogenetic relationships among iguanian lizards in the North American genus Sceloporus. We tested for diversification rate shifts to determine if rapid radiation in the genus is correlated with chromosomal evolution.

Results

The phylogenomic trees that we obtained for Sceloporus using concatenation and coalescent-based species tree inference provide strong support for the monophyly and interrelationships among nearly all major groups. The diversification analysis supported one rate shift on the Sceloporus phylogeny approximately 20–25 million years ago that is associated with the doubling of the speciation rate from 0.06 species/million years (Ma) to 0.15 species/Ma. The posterior probability for this rate shift occurring on the branch leading to the Sceloporus species groups exhibiting increased chromosomal diversity is high (posterior probability = 0.997).

Conclusions

Despite high levels of gene tree discordance, we were able to estimate a phylogenomic tree for Sceloporus that solves some of the taxonomic problems caused by previous analyses of fewer loci. The taxonomic changes that we propose using this new phylogenomic tree help clarify the number and composition of the major species groups in the genus. Our study provides new evidence for a putative link between chromosomal evolution and the rapid divergence and radiation of Sceloporus across North America.

Electronic supplementary material

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

Phylogeography of Pogonomyrmex barbatus and P. rugosus harvester ants with genetic and environmental caste determination

We present a phylogeographic study of at least six reproductively isolated lineages of new world harvester ants within the Pogonomyrmex barbatus and P. rugosus species group. The genetic and geographic relationships within this clade are complex: Four of the identified lineages show genetic caste determination (GCD) and are divided into two pairs. Each pair has evolved under a mutualistic system that necessitates sympatry. These paired lineages are dependent upon one another because their GCD requires interlineage matings for the production of F1 hybrid workers, and intralineage matings are required to produce queens. This GCD system maintains genetic isolation among these interdependent lineages, while simultaneously requiring co-expansion and emigration as their distributions have changed over time. It has also been demonstrated that three of these four GCD lineages have undergone historical hybridization, but the narrower sampling range of previous studies has left questions on the hybrid parentage, breadth, and age of these groups. Thus, reconstructing the phylogenetic and geographic history of this group allows us to evaluate past insights and hypotheses and to plan future inquiries in a more complete historical biogeographic context. Using mitochondrial DNA sequences sampled across most of the morphospecies' ranges in the U.S.A. and Mexico, we conducted a detailed phylogeographic study. Remarkably, our results indicate that one of the GCD lineage pairs has experienced a dramatic range expansion, despite the genetic load and fitness costs of the GCD system. Our analyses also reveal a complex pattern of vicariance and dispersal in Pogonomyrmex harvester ants that is largely concordant with models of late Miocene, Pliocene, and Pleistocene range shifts among various arid-adapted taxa in North America.

Testing Taxon Tenacity of Tortoises: evidence for a geographical selection gradient at a secondary contact zone

We examined a secondary contact zone between two species of desert tortoise, Gopherus agassizii and G. morafkai. The taxa were isolated from a common ancestor during the formation of the Colorado River (4–8 mya) and are a classic example of allopatric speciation. However, an anomalous population of G. agassizii comes into secondary contact with G. morafkai east of the Colorado River in the Black Mountains of Arizona and provides an opportunity to examine reinforcement of species' boundaries under natural conditions. We sampled 234 tortoises representing G. agassizii in California (n - 103), G. morafkai in Arizona (n - 78), and 53 individuals of undetermined assignment in the contact zone including and surrounding the Black Mountains. We genotyped individuals for 25 STR loci and determined maternal lineage using mtDNA sequence data. We performed multilocus genetic clustering analyses and used multiple statistical methods to detect levels of hybridization. We tested hypotheses about habitat use between G. agassizii and G. morafkai in the region where they co-occur using habitat suitability models. Gopherus agassizii and G. morafkai maintain independent taxonomic identities likely due to ecological niche partitioning, and the maintenance of the hybrid zone is best described by a geographical selection gradient model.

Comparative species divergence across eight triplets of spiny lizards (Sceloporus) using genomic sequence data

Species divergence is typically thought to occur in the absence of gene flow, but many empirical studies are discovering that gene flow may be more pervasive during species formation. Although many examples of divergence with gene flow have been identified, few clades have been investigated in a comparative manner, and fewer have been studied using genome-wide sequence data. We contrast species divergence genetic histories across eight triplets of North American Sceloporus lizards using a maximum likelihood implementation of the isolation–migration (IM) model. Gene flow at the time of species divergence is modeled indirectly as variation in species divergence time across the genome or explicitly using a migration rate parameter. Likelihood ratio tests (LRTs) are used to test the null model of no gene flow at speciation against these two alternative gene flow models. We also use the Akaike information criterion to rank the models. Hundreds of loci are needed for the LRTs to have statistical power, and we use genome sequencing of reduced representation libraries to obtain DNA sequence alignments at many loci (between 340 and 3,478; mean = 1,678) for each triplet. We find that current species distributions are a poor predictor of whether a species pair diverged with gene flow. Interrogating the genome using the triplet method expedites the comparative study of species divergence history and the estimation of genetic parameters associated with speciation.

Cryptic diversity in host-associated populations of Tetra pinnatifidae (Acari: Eriophyoidea): what do morphometric, mitochondrial and nuclear data reveal and conceal?

Traditional morphology-based taxonomy of eriophyoid mites (Acari: Eriophyoidea) has been challenged by molecular-based technologies in the detection of cryptic species. However, the implications of such cryptic diversity appear to differ when methods based on different types of data are used. Here, samples of a host-associated eriophyoid mite species, Tetra pinnatifidae, collected from different host plants and localities are evaluated. The congruence of results based on morphometric (32 characters), mitochondrial (16S), and nuclear (28S) data were evaluated and showed a host-associated cryptic diversity dividing this morphospecies into several groups/clades that were morphometrically indistinguishable. In comparison, the 16S data confirmed cryptic speciation and intra-clade host-associated diversity, while 28S did not. In contrast, 28S data revealed potential gene flow between host-associated populations. High mitochondrial divergence, as well as low nuclear and morphological divergence indicated very recent stage of cryptic diversity of this eriophyoid mite.

Species delimitation and digit number in a North African skink

Delimitation of species is an important and controversial area within evolutionary biology. Many species boundaries have been defined using morphological data. New genetic approaches now offer more objective evaluation and assessment of the reliability of morphological variation as an indicator that speciation has occurred. We examined geographic variation in morphology of the continuously distributed skink Chalcides mionecton from Morocco and used Bayesian analyses of nuclear and mitochondrial DNA (mtDNA) loci to examine: (i) their concordance with morphological patterns, (ii) support for species delimitation, (iii) timing of speciation, and (iv) levels of gene flow between species. Four digit individuals were found at sites between Cap Rhir (in the south) and the northern extreme of the range, whereas five-digit individuals were found in two disjunct areas: (i) south of Cap Rhir and (ii) the north of the range where they were often syntopic with four-digit individuals. The pattern of variation in generalized body dimensions was largely concordant with that in digit number, suggesting two general morphotypes. Bayesian analyses of population structure showed that individuals from sites south of Cap Rhir formed one genetic cluster, but that northern four- and five-digit individuals clustered together. Statistical support for delimitation of these genetic clusters into two species was provided by a recent Bayesian method. Phylogenetic-coalescent dating with external time calibrations indicates that speciation was relatively recent, with a 95% posterior interval of 0.46-2.66 mya. This postdates equivalent phylogenetic dating estimates of sequence divergence by approximately 1 Ma. Statistical analyses of a small number of independent loci provide important insights into the history of the speciation process in C. mionecton and support delimitation of populations into two species with distributions that are spatially discordant with patterns of morphological variation.

Pleistocene speciation with and without gene flow in Euphaea damselflies of subtropical and tropical East Asian islands

Climatic oscillations during the Pleistocene period could have had a profound impact on the origin of tropical species by the alternation of allopatric isolation and interpopulation gene flow cycles. However, whether tropical speciation involves strictly allopatric isolation, or proceeds in the face of homogenizing gene flow, is relatively unclear. Here, we investigated geographical modes of speciation in four closely related Euphaea damselfly species endemic to the subtropical and tropical East Asian islands using coalescent analyses of a multilocus data set. The reconstructed phylogenies demonstrated distinct species status for each of the four species and the existence of two sister species pairs, Euphaea formosa/E. yayeyamana and E. decorata/E. ornata. The species divergence time of the sibling Euphaea damselflies dates back to within the last one Mya of the Middle to Lower Pleistocene. The speciation between the populous E. formosa of Taiwan and the less numerous E. yayeyamana of the Yaeyama islands occurred despite significant bidirectional, asymmetric gene flow, which is strongly inconsistent with a strictly allopatric model. In contrast, speciation of the approximately equal-sized populations of E. decorata of the southeast Asian mainland and E. ornata of Hainan is inferred to have involved allopatric divergence without gene flow. Our findings suggest that differential selection of natural or sexual environments is a prominent driver of species divergence in subtropical E. formosa and E. yayeyamana; whereas for tropical E. decorata and E. ornata at lower latitudes, allopatric isolation may well be a pivotal promoter of species formation.

Multi-locus estimates of population structure and migration in a fence lizard hybrid zone

A hybrid zone between two species of lizards in the genus Sceloporus (S. cowlesi and S. tristichus) on the Mogollon Rim in Arizona provides a unique opportunity to study the processes of lineage divergence and merging. This hybrid zone involves complex interactions between 2 morphologically and ecologically divergent subspecies, 3 chromosomal groups, and 4 mitochondrial DNA (mtDNA) clades. The spatial patterns of divergence between morphology, chromosomes and mtDNA are discordant, and determining which of these character types (if any) reflects the underlying population-level lineages that are of interest has remained impeded by character conflict. The focus of this study is to estimate the number of populations interacting in the hybrid zone using multi-locus nuclear data, and to then estimate the migration rates and divergence time between the inferred populations. Multi-locus estimates of population structure and gene flow were obtained from 12 anonymous nuclear loci sequenced for 93 specimens of Sceloporus. Population structure estimates support two populations, and this result is robust to changes to the prior probability distribution used in the Bayesian analysis and the use of spatially-explicit or non-spatial models. A coalescent analysis of population divergence suggests that gene flow is high between the two populations, and that the timing of divergence is restricted to the Pleistocene. The hybrid zone is more accurately described as involving two populations belonging to S. tristichus, and the presence of S. cowlesi mtDNA haplotypes in the hybrid zone is an anomaly resulting from mitochondrial introgression.

Phylogeography of supralittoral rocky intertidal Ligia isopods in the pacific region from central California to central Mexico

Background

Ligia isopods are widely distributed in the Pacific rocky intertidal shores from central California to central Mexico, including the Gulf of California. Yet, their biological characteristics restrict them to complete their life cycles in a very narrow range of the rocky intertidal supralittoral. Herein, we examine phylogeographic patterns of Ligia isopods from 122 localities between central California and central Mexico. We expect to find high levels of allopatric diversity. In addition, we expect the phylogeographic patterns to show signatures of past vicariant events that occurred in this geologically dynamic region.

Methodology/Principal Findings

We sequenced two mitochondrial genes (Cytochrome Oxidase I and 16S ribosomal DNA). We conducted Maximum Likelihood and Bayesian phylogenetic analyses. We found many divergent clades that, in general, group according to geography. Some of the most striking features of the Ligia phylogeographic pattern include: (1) deep mid-peninsular phylogeographic breaks on the Pacific and Gulf sides of Baja peninsula; (2) within the Gulf lineages, the northern peninsula is most closely related to the northern mainland, while the southern peninsula is most closely related to the central-southern mainland; and, (3) the southernmost portion of the peninsula (Cape Region) is most closely related to the southernmost portion of mainland.

Conclusions/Significance

Our results shed light on the phylogenetic relationships of Ligia populations in the study area. This study probably represents the finest-scale phylogeographic examination for any organism to date in this region. Presence of highly divergent lineages suggests multiple Ligia species exist in this region. The phylogeographic patterns of Ligia in the Gulf of California and Baja peninsula are incongruent with a widely accepted vicariant scenario among phylogeographers, but consistent with aspects of alternative geological hypotheses and phylo- and biogeographic patterns of several other taxa. Our findings contribute to the ongoing debate regarding the geological origin of this important biogeographic region.

Cryptic failure of partitioned Bayesian phylogenetic analyses: lost in the land of long trees

Partitioned Bayesian phylogenetic analyses of routine genetic data sets, constructed using MrBayes (Ronquist and Huelsenbeck 2003), can become trapped in regions of parameter space characterized by unrealistically long trees and distorted partition rate multipliers. Such analyses commonly fail to reach stationarity during hundreds of millions of generations of sampling-many times longer than most published analyses. Some data sets are so prone to this problem that paired MrBayes runs begun from different starting trees repeatedly find the same incorrect long-tree solutions and consequently pass the most commonly employed tests of stationarity, including the average standard deviation of split frequencies (ASDSF) and the potential scale reduction factor (PSRF) statistics offered by MrBayes (Gelman and Rubin 1992). In these situations, failure to reach stationarity is recognizable only in light of prior knowledge of model parameters, such as the expectation that third-codon-position sites usually evolve fastest in protein-coding genes. The conditions that lead to the long-tree problem are frequently encountered in phylogenetic studies today, and I present 6 demonstration examples from the literature. Although the effects on tree length (TL) are often dramatic, effects on topology appear to be subtle. Susceptibility to the problem is sometimes predicted by the difference between the true TL and the starting TL. In some cases, the problems described here can be avoided or reduced by manipulation of the starting TL and/or by adjustments to the prior on branch lengths. In more difficult situations, accurate branch length estimation may not be possible with Bayesian methods because of dependence of the solution on the branch length prior.

Delimiting species boundaries within the Neotropical bamboo Otatea (Poaceae: Bambusoideae) using molecular, morphological and ecological data

Species delimitation is a task that has engaged taxonomists for more than two centuries. Recently, it has been demonstrated that molecular data and ecological niche modeling are useful in species delimitation. In this paper multiple data sets (molecular, morphological, ecological) were utilized to set limits for the species belonging to the Neotropical bamboo Otatea, because there is disagreement about species circumscriptions and also because the genus has an interesting distribution, with most of its populations in Mexico and a single disjunct population in Colombia. Molecular and morphological phylogenetic analyses recovered trees with conflicting topologies. Tree-based morphological and character-based analyses recognized the same entities. Ecological niche models and PCA/MANOVAS agreed with the recognition of the same entities that resulted from the morphological analyses. Morphological analyses retrieved clades supported by diagnostic characters and coherent geographical distributions. Based on these results seven entities should be recognized in Otatea, instead of the three previously described species.

Mustela or Vison? Evidence for the taxonomic status of the American mink and a distinct biogeographic radiation of American weasels

The American mink's relationship to the weasels in Mustela has been uncertain. Karyological, morphological, and phylogenetic comparisons to Eurasian Mustela support placing the mink outside the genus as Neovison vison. However, genetic comparisons that incorporate other endemic American Mustela suggest the interpretation of N. vison's position to Mustela has been handicapped by biased geographic sampling. Here, we analyzed mitochondrial cytochrome-b from all weasels endemic to the Americas, including two poorly known South American species (M. felipei, M. africana), weasels native to North America (M. vison, M. frenata, M. nigripes), Mustela migrant to North America (M. erminea, M. nivalis), palearctic Mustela, and other American members of Mustelidae. Bayesian and likelihood inference methods were used to construct a phylogeny of Mustela, and relaxed Bayesian phylogenetic techniques estimated ages of divergence within the genus using priors calibrated by fossil ages. Our analyses show that the American mink and the smaller Mustela endemic to the Americas represent a distinct phylogenetic heritage apart from their Eurasian cousins, and biogeographic barriers like the Bering and Panamanian land bridges have influenced the evolutionary history of Mustela in the Americas.

Multigene analysis of phylogenetic relationships and divergence times of primate sucking lice (Phthiraptera: Anoplura)

Cospeciation between hosts and parasites offers a unique opportunity to use information from parasites to infer events in host evolutionary history. Although lice (Insecta: Phthiraptera) are known to cospeciate with their hosts and have frequently served as important markers to infer host evolutionary history, most molecular studies are based on only one or two markers. Resulting phylogenies may, therefore, represent gene histories (rather than species histories), and analyses of multiple molecular markers are needed to increase confidence in the results of phylogenetic analyses. Herein, we phylogenetically examine nine molecular markers in primate sucking lice (Phthiraptera: Anoplura) and we use these markers to estimate divergence times among louse lineages. Individual and combined analyses of these nine markers are, for the most part, congruent, supporting relationships hypothesized in previous studies. Only one marker, the nuclear protein-coding gene Histone 3, has a significantly different tree topology compared to the other markers. The disparate evolutionary history of this marker, however, has no significant effect on topology or nodal support in the combined phylogenetic analyses. Therefore, phylogenetic results from the combined data set likely represent a solid hypothesis of species relationships. Additionally, we find that simultaneous use of multiple markers and calibration points provides the most reliable estimates of louse divergence times, in agreement with previous studies estimating divergences among species. Estimates of phylogenies and divergence times also allow us to verify the results of [Reed, D.L., Light, J.E., Allen, J.M., Kirchman, J.J., 2007. Pair of lice lost or parasites regained: the evolutionary history of anthropoid primate lice. BMC Biol. 5, 7.]; there was probable contact between gorilla and archaic hominids roughly 3 Ma resulting in a host switch of Pthirus lice from gorillas to archaic hominids. Thus, these results provide further evidence that data from cospeciating organisms can yield important information about the evolutionary history of their hosts.

Birth of a biome: insights into the assembly and maintenance of the Australian arid zone biota

The integration of phylogenetics, phylogeography and palaeoenvironmental studies is providing major insights into the historical forces that have shaped the Earth's biomes. Yet our present view is biased towards arctic and temperate/tropical forest regions, with very little focus on the extensive arid regions of the planet. The Australian arid zone is one of the largest desert landform systems in the world, with a unique, diverse and relatively well-studied biota. With foci on palaeoenvironmental and molecular data, we here review what is known about the assembly and maintenance of this biome in the context of its physical history, and in comparison with other mesic biomes. Aridification of Australia began in the Mid-Miocene, around 15 million years, but fully arid landforms in central Australia appeared much later, around 1-4 million years. Dated molecular phylogenies of diverse taxa show the deepest divergences of arid-adapted taxa from the Mid-Miocene, consistent with the onset of desiccation. There is evidence of arid-adapted taxa evolving from mesic-adapted ancestors, and also of speciation within the arid zone. There is no evidence for an increase in speciation rate during the Pleistocene, and most arid-zone species lineages date to the Pliocene or earlier. The last 0.8 million years have seen major fluctuations of the arid zone, with large areas covered by mobile sand dunes during glacial maxima. Some large, vagile taxa show patterns of recent expansion and migration throughout the arid zone, in parallel with the ice sheet-imposed range shifts in Northern Hemisphere taxa. Yet other taxa show high lineage diversity and strong phylogeographical structure, indicating persistence in multiple localised refugia over several glacial maxima. Similar to the Northern Hemisphere, Pleistocene range shifts have produced suture zones, creating the opportunity for diversification and speciation through hybridisation, polyploidy and parthenogenesis. This review highlights the opportunities that development of arid conditions provides for rapid and diverse evolutionary radiations, and re-enforces the emerging view that Pleistocene environmental change can have diverse impacts on genetic structure and diversity in different biomes. There is a clear need for more detailed and targeted phylogeographical studies of Australia's arid biota and we suggest a framework and a set of a priori hypotheses by which to proceed.