Genetic structure of desert ground squirrels over a 20-degree-latitude transect from Oregon through the Baja California peninsula

Climatic dissimilarity associated with phylogenetic breaks

Shared phylogenetic breaks often are associated with clear geographic barriers but some common phylogeographic breaks may lack obvious underlying mechanisms. A phylogenetic break involving multiple taxa was found in the Baja California Peninsula that was associated with a past sea barrier. However, geological evidence is lacking for this barrier’s past existence, and despite its current absence, the genetic breaks have persisted. This work explores the relationships between the current climatic niches for matrilineages of 11 vertebrate species as a possible explanation for the current geographic partitioning of matrilineages. We evaluated the climatic occupancy of each matrilineage through ecological niche models, background similarity, niche overlap, niche divergence, and Mantel tests. We found disparities in the climatic occupancy between north and south matrilineage of each taxon. Northern matrilineages are associated with lower temperatures and winter rains, while southern matrilineages reside in areas with higher temperatures and summer rains.

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.

The endemic insular and peninsular species Chaetodipus spinatus (Mammalia, Heteromyidae) breaks patterns for Baja California

The Baja California peninsula is the second longest, most geographically isolated peninsula on Earth. Its physiography and the presence of many surrounding islands has facilitated studies of the underlying patterns and drivers of genetic structuring for a wide spectrum of organisms. Chaetodipus spinatus is endemic to the region and occurs on 12 associated islands, including 10 in the Gulf of California and two in the Pacific Ocean. This distribution makes it a model species for evaluating natural historical barriers. We test hypotheses associated with the relationship between the range of the species, patterns in other species, and its relationship to Pleistocene-Holocene climatic changes. We analyzed sequence data from mtDNA genes encoding cytochrome b (Cytb) and cytochrome c oxidase subunits I (COI) and III (COIII) in 26 populations including all 12 islands. The matrilineal genealogy, statistical parsimony network and Bayesian skyline plot indicated an origin of C. spinatus in the southern part of the peninsula. Our analyses detected several differences from the common pattern of peninsular animals: no mid-peninsula break exists, Isla Carmen hosts the most divergent population, the population on an ancient southern Midriff island does not differ from peninsular populations, and a mtDNA peninsular discordance occurs near Loreto.

Establishing a new animal model for hepadnaviral infection: susceptibility of Chinese Marmota-species to woodchuck hepatitis virus infection

Hepatitis B virus infection (HBV) is a major medical problem in China. The lack of a suitable infection model in China is recognized as an obstacle for research on HBV in China. Chinese Marmota-species is phylogenetically closely related to Marmota monax, thus, it might be suitable to serve as an animal model for HBV infection. Therefore, we attempted to prove the claim about the existence of woodchuck hepatitis virus (WHV)-like viruses in Chinese Marmota-species and to determine the susceptibility of these species to experimental WHV infection. In the present study, 653 sera from three Chinese Marmota-species, Marmota himalayana, Marmota baibacina and Marmota bobak, were screened for WHV-like viruses by serological and molecular assays. The susceptibility to WHV of three species was investigated by experimental infection and monitored by testing of anti-WHc and WHsAg by ELISA, detection of WHV DNA by PCR, and detection of WHV replication intermediates and antigens in liver samples. No evidence for the existence of a genetically closely related virus to WHV in three Chinese Marmota-species was found by serological assays and PCR. M. himalayana was susceptible to WHV infection as inoculated animals became positive for anti-WHc, WHsAg and WHV DNA. Further, WHV replication intermediates and proteins were detected in liver samples. In contrast, M. baibacina remained negative for tested virological parameters. M. bobak species showed a limited susceptibility to WHV. Our data do not support early reports about WHV-like viruses in China. M. himalayana is suitable for the establishment of a model for hepadnaviral infection.

Fossils and phylogeny uncover the evolutionary history of a unique antipredator behaviour

Recently, two squirrel species (Spermophilus spp.) were discovered to anoint their bodies with rattlesnake scent as a means of concealing their odour from these chemosensory predators. In this study, we tested multiple species with predator scents (rattlesnake and weasel) to determine the prevalence of scent application across the squirrel phylogeny. We reconstructed the evolutionary history of the behaviour using a phylogenetic analysis and fossil records of historic predator co-occurrence. Squirrels with historical and current rattlesnake co-occurrence all applied rattlesnake scent, whereas no relationship existed between weasel scent application and either weasel or rattlesnake co-occurrence. This was surprising because experimental tests confirmed rattlesnake and weasel scent were both effective at masking prey odour from hunting rattlesnakes (the primary predator of squirrels). Ancestral reconstructions and fossil data suggest predator scent application in squirrels is ancient in origin, arising before co-occurrences with rattlesnakes or weasels in response to some other, now extinct, chemosensory predator.

Fire and ice: genetic structure of the Uinta ground squirrel (Spermophilus armatus) across the Yellowstone hotspot

The range of the Uinta ground squirrel, Spermophilus armatus, is centred over one of the most tectonically active regions today, the Yellowstone hotspot. We document the role of Quaternary tectonic and climatic history on the genetic structure of this species by screening museum and extant individuals throughout its range. Phylogeographic, divergence time, and demographic analyses of partial mitochondrial cytochrome b and control region DNA sequences yield insight into the cadence of evolution across three spatiotemporal scales: (i) a relatively deep intraspecific divergence of S. armatus into three lineages coincident with the last major volcanic eruption in the region and maintained by the Snake River Plain; (ii) demographic expansion in two lineages corresponding to the time of last deglaciation of the region; and (iii) a recent (< 50 years) local extinction of the third lineage coincident with climatic change and conversion of habitat for agricultural purposes in eastern Idaho. Beyond these inferences, our study highlights the unique value of museum material to phylogeography, and shows that small mammal recolonization of previously glaciated montane 'islands' differs from northward postglacial expansion observed in areas previously covered by continental ice sheets. Montane 'islands' may harbour high genetic diversity because of admixture and recurrent expansion/extinction.

Two waves of diversification in mammals and reptiles of Baja California revealed by hierarchical Bayesian analysis

Many species inhabiting the Peninsular Desert of Baja California demonstrate a phylogeographic break at the mid-peninsula, and previous researchers have attributed this shared pattern to a single vicariant event, a mid-peninsular seaway. However, previous studies have not explicitly considered the inherent stochasticity associated with the gene-tree coalescence for species preceding the time of the putative mid-peninsular divergence. We use a Bayesian analysis of a hierarchical model to test for simultaneous vicariance across co-distributed sister lineages sharing a genealogical break at the mid-peninsula. This Bayesian method is advantageous over traditional phylogenetic interpretations of biogeography because it considers the genetic variance associated with the coalescent and mutational processes, as well as the among-lineage demographic differences that affect gene-tree coalescent patterns. Mitochondrial DNA data from six small mammals and six squamate reptiles do not support the perception of a shared vicariant history among lineages exhibiting a north-south divergence at the mid-peninsula, and instead support two events differentially structuring genetic diversity in this region.

Phylogeographical structure in the subterranean tuco-tuco Ctenomys talarum (Rodentia: Ctenomyidae): contrasting the demographic consequences of regional and habitat-specific histories

In this work we examined the phylogeography of the South American subterranean herbivorous rodent Ctenomys talarum (Talas tuco-tuco) using mitochondrial DNA (mtDNA) control region (D-loop) sequences, and we assessed the geographical genetic structure of this species in comparison with that of subterranean Ctenomys australis, which we have shown previously to be parapatric to C. talarum and to also live in a coastal sand dune habitat. A significant apportionment of the genetic variance among regional groups indicated that putative geographical barriers, such as rivers, substantially affected the pattern of genetic structure in C. talarum. Furthermore, genetic differentiation is consistent with a simple model of isolation by distance, possibly evidencing equilibrium between gene flow and local genetic drift. In contrast, C. australis showed limited hierarchical partitioning of genetic variation and departed from an isolation-by-distance pattern. Mismatch distributions and tests of neutrality suggest contrasting histories of these two species: C. talarum appears to be characterized by demographic stability and no significant departures from neutrality, whereas C. australis has undergone a recent demographic expansion and/or departures from strict neutrality in its mtDNA.

Phylogeography of the Western Lyresnake (Trimorphodon biscutatus): testing aridland biogeographical hypotheses across the Nearctic-Neotropical transition

The Western Lyresnake (Trimorphodon biscutatus) is a widespread, polytypic taxon inhabiting arid regions from the warm deserts of the southwestern United States southward along the Pacific versant of Mexico to the tropical deciduous forests of Mesoamerica. This broadly distributed species provides a unique opportunity to evaluate a priori biogeographical hypotheses spanning two major distinct biogeographical realms (the Nearctic and Neotropical) that are usually treated separately in phylogeographical analyses. I investigated the phylogeography of T. biscutatus using maximum likelihood and Bayesian phylogenetic analysis of mitochondrial DNA (mtDNA) from across this species' range. Phylogenetic analyses recovered five well-supported clades whose boundaries are concordant with existing geographical barriers, a pattern consistent with a model of vicariant allopatric divergence. Assuming a vicariance model, divergence times between mitochondrial lineages were estimated using Bayesian relaxed molecular clock methods calibrated using geological information from putative vicariant events. Divergence time point estimates were bounded by broad confidence intervals, and thus these highly conservative estimates should be considered tentative hypotheses at best. Comparison of mtDNA lineages and taxa traditionally recognized as subspecies based on morphology suggest this taxon is comprised of multiple independent lineages at various stages of divergence, ranging from putative secondary contact and hybridization to sympatry of 'subspecies'.

Phylogeography of the montane caddisflyDrusus discolor: evidence for multiple refugia and periglacial survival

We studied the genetic population structure and phylogeography of the montane caddisfly Drusus discolor across its entire range in central and southern Europe. The species is restricted to mountain regions and exhibits an insular distribution across the major mountain ranges. Mitochondrial sequence data (COI) of 254 individuals from the entire species range is analysed to reveal population genetic structure. The data show little molecular variation within populations and regions, but distinct genetic differentiation between mountain ranges. Most populations are significantly differentiated based on F(ST) and exact tests of population differentiation and most haplotypes are unique to a single mountain range. Phylogenetic analyses reveal deep divergence between geographically isolated lineages. Combined, these results suggest that past fragmentation is the prominent process structuring the populations across Europe. We use tests of selective neutrality and mismatch distributions, to study the demographic population history of regions with haplotype overlap. The high level of genetic differentiation between mountain ranges and estimates of demographic history provide evidence for the existence of multiple glacial refugia, including several in central Europe. The study shows that these aquatic organisms reacted differently to Pleistocene cooling than many terrestrial species. They persisted in numerous refugia over multiple glacial cycles, allowing many local endemic clades to form.

Phylogeography of the black fly Simulium tani (Diptera: Simuliidae) from Thailand as inferred from mtDNA sequences

Intraspecific phylogeography has been used widely as a tool to infer population history. However, little attention has been paid to Southeast Asia despite its importance in terms of biodiversity. Here we used the cytochrome oxidase I gene of mitochondrial DNA (mtDNA) for a phylogeographic study of 147 individuals of the black fly Simulium tani from Thailand. The mtDNA revealed high genetic differentiation between the major geographical regions of north, east and central/south Thailand. Mismatch distributions indicate population expansions during the mid-Pleistocene and the late Pleistocene suggesting that current population structure and diversity may be due in part to the species' response to Pleistocene climatic fluctuations. The genealogical structure of the haplotypes, high northern diversity and maximum-likelihood inference of historical migration rates, suggest that the eastern and central/southern populations originated from northern populations in the mid-Pleistocene. Subsequently, the eastern region had had a largely independent history but the central/southern population may be largely the result of recent (c. 100,000 years ago) expansion, either from the north again, or from a relictual population in the central region. Cytological investigation revealed that populations from the south and east have two overlapping fixed chromosomal inversions. Since these populations also share ecological characteristics it suggests that inversions are involved in ecological adaptation. In conclusion both contemporary and historical ecological conditions are playing an important role in determining population genetic structure and diversity.

Studies of morphological and molecular phylogenetic divergence in spiders (Araneae: Homalonychus) from the American southwest, including divergence along the Baja California Peninsula

Comparative phylogenetic and phylogeographic analyses have revealed a pervasive midpeninsular divergence in the mitochondrial genealogies of numerous vertebrate taxa distributed on the Baja California Peninsula. In this study, we extend the investigation of regional vicariance in Baja California to an arthropod taxon by examining patterns of phylogenetic and morphological divergence in the spider genus Homalonychus (Araneae, Homalonychidae). We analyzed data from two mtDNA genes (16S rRNA and NADH dehydrogenase subunit (1) and a nuclear gene (28S rRNA) using maximum parsimony and Bayesian phylogenetic analyses, and also conducted geometric morphometric analyses employing landmark data on male and female genitalia. Genes and morphology both reveal a deep split across the Colorado River and Gulf of California, separating Homalonychus selenopoides on the east side of river from its congener Homalonychus theologus on the west side of the river, including the Baja California Peninsula. Along the north-south axis of the Baja Peninsula, an apparently more recent midpeninsular phylogenetic break is evident within H. theologus in the mitochondrial genome and in female genitalia. However, there is no measurable divergence between northern and southern populations in either nuclear DNA or male genitalia. We suggest that this discordance between datasets reflects either a difference in rates of evolution between male versus female systems, or that male-based nuclear gene flow is obscuring a phylogenetic split that is fixed in the female-based systems. Our findings provide additional support for a midpeninsular Baja divergence event, although the timing and geological evidence for such an event remain elusive.