Congruent morphological and genetic differentiation as a signature of range expansion in a fragmented landscape

Genetic footprints of a rapid and large-scale range expansion: the case of cyclic common vole in Spain

In the Anthropocene, many species are rapidly shifting their ranges in response to human-driven habitat modifications. Studying patterns and genetic signatures of range shifts helps to understand how species cope with environmental disturbances and predict future shifts in the face of global environmental change. We investigated the genetic signature of a contemporary wide-range expansion observed in the Iberian common vole Microtus arvalis asturianus shortly after a colonization event. We used mtDNA and microsatellite data to investigate patterns of genetic diversity, structure, demography, and gene flow across 57 localities covering the historical range of the species and the newly colonized area. The results showed a genetic footprint more compatible with a true range expansion (i.e. the colonization of previously unoccupied areas), than with a model of "colonization from within" (i.e. local expansions from small, unnoticed populations). Genetic diversity measures indicated that the source population was likely located at the NE of the historical range, with a declining gradient of genetic diversity towards the more recently invaded areas. At the expansion front, we observed the greatest gene flow and smallest pairwise differences between nearby localities. Both natural landscape features (rivers) and recent anthropogenic barriers (roads, railways) explained a large proportion of genetic variance among populations and had a significant impact on the colonization pathways used by voles.

Phenotypic plasticity determines differences between the skulls of tigers from mainland Asia

Tiger subspecific taxonomy is controversial because of morphological and genetic variation found between now fragmented populations, yet the extent to which phenotypic plasticity or genetic variation affects phenotypes of putative tiger subspecies has not been explicitly addressed. In order to assess the role of phenotypic plasticity in determining skull variation, we compared skull morphology among continental tigers from zoos and the wild. In turn, we examine continental tiger skulls from across their wild range, to evaluate how the different environmental conditions experienced by individuals in the wild can influence morphological variation. Fifty-seven measurements from 172 specimens were used to analyse size and shape differences among wild and captive continental tiger skulls. Captive specimens have broader skulls, and shorter rostral depths and mandible heights than wild specimens. In addition, sagittal crest size is larger in wild Amur tigers compared with those from captivity, and it is larger in wild Amur tigers compared with other wild continental tigers. The degree of phenotypic plasticity shown by the sagittal crest, skull width and rostral height suggests that the distinctive shape of Amur tiger skulls compared with that of other continental tigers is mostly a phenotypically plastic response to differences in their environments.

Morphological variation in bony structures of Hylaeamys seuanezi (Rodentia, Cricetidae, Sigmodontinae) in Southern Bahia, Brazil

Rodents of the genus Hylaeamys, are a group of cryptic species previously included in the Oryzomys capito complex. In Brazil are represented by six species, distributed in different biomes. However, the limits of their geographic distributions and taxonomy have not yet been well defined. In particular, the taxonomy for Hylaeamys seuanezi is unstable. Based on the analysis of bony structures and supported with geometric morphometric techniques, we characterize and compare the average shape from populations in four localities from Southern Bahia, Brazil. We review 145 individuals and we create morphological landmarks in skulls, mandibles, scapulae, and pelvis. In all the structures there were statistically significant differences between populations, in which the average shape from the Igrapiúna population was the most differentiated. Our results also showed differences between the pelvis of males and females, reporting for the first time sexual dimorphism for H. seuanezi. Finally, we provide a morphological diagnosis between the populations and postulate that such differences may be correlated with environmental and climatic factors that could be exerting negative pressures on H. seuanezi; as has been evidenced with other species of rodents and other mammals.

Disentangling the genetic effects of refugial isolation and range expansion in a trans-continentally distributed species

In wide-ranging taxa with historically dynamic ranges, past allopatric isolation and range expansion can both influence the current structure of genetic diversity. Considering alternate historical scenarios involving expansion from either a single refugium or from multiple refugia can be useful in differentiating the effects of isolation and expansion. Here, we examined patterns of genetic variability in the trans-continentally distributed painted turtle (Chrysemys picta). We utilized an existing phylogeographic dataset for the mitochondrial control region and generated additional data from nine populations for the mitochondrial control region (n = 302) and for eleven nuclear microsatellite loci (n = 247). We created a present-day ecological niche model (ENM) for C. picta and hindcast this model to three reconstructions of historical climate to define three potential scenarios with one, two, or three refugia. Finally, we employed spatially-explicit coalescent simulations and an approximate Bayesian computation (ABC) framework to test which scenario best fit the observed genetic data. Simulations indicated that phylogeographic and multilocus population-level sampling both could differentiate among refugial scenarios, although inferences made using mitochondrial data were less accurate when a longer coalescence time was assumed. Furthermore, all empirical genetic datasets were most consistent with expansion from a single refugium based on ABC. Our results indicate a stronger role for post-glacial range expansion, rather than isolation in allopatric refugia followed by range expansion, in structuring diversity in this species. To distinguish among complex historical scenarios, we recommend explicitly modeling the effects of range expansion and evaluating alternate refugial scenarios for wide-ranging taxa.

Admixture on the northern front: population genomics of range expansion in the white-footed mouse (Peromyscus leucopus) and secondary contact with the deer mouse (Peromyscus maniculatus)

Range expansion has genetic consequences expected to result in differentiated wave-front populations with low genetic variation and potentially introgression from a local species. The northern expansion of Peromyscus leucopus in southern Quebec provides an opportunity to test these predictions using population genomic tools. Our results show evidence of recent and post-glacial expansion. Genome-wide variation in P. leucopus indicates two post-glacial lineages are separated by the St. Lawrence River, with a more recent divergence of populations isolated by the Richelieu River. In two of three transects we documented northern populations with low diversity in at least one genetic measure, although most relationships were not significant. Consistent with bottlenecks and allele surfing during northward expansion, we document a northern-most population with low nucleotide diversity, divergent allele frequencies and the most private alleles, and observed heterozygosity indicates outcrossing. Ancestry proportions revealed putative hybrids of P. leucopus and P. maniculatus. A formal test for gene flow confirmed secondary contact, showing that a reticulate population phylogeny between P. maniculatus and P. leucopus was a better fit to the data than a bifurcating model without gene flow. Thus, we provide the first genomic evidence of gene flow between this pair of species in natural populations. Understanding the evolutionary consequences of secondary contact is an important conservation concern as climate-induced range expansions are expected to result in new hybrid zones between closely related species.

Liver microbiome of Peromyscus leucopus, a key reservoir host species for emerging infectious diseases in North America

Microbiome studies generally focus on the gut microbiome, which is composed of a large proportion of commensal bacteria. Here we propose a first analysis of the liver microbiome using next generation sequencing as a tool to detect potentially pathogenic strains. We used Peromyscus leucopus, the main reservoir host species of Lyme disease in eastern North America, as a model and sequenced V5-V6 regions of the 16S gene from 18 populations in southern Quebec (Canada). The Lactobacillus genus was found to dominate the liver microbiome. We also detected a large proportion of individuals infected by Bartonella vinsonii arupensis, a human pathogenic bacteria responsible for endocarditis, as well as Borrelia burgdorferi, the pathogen responsible for Lyme disease in North America. We then compared the microbiomes among two P. leucopus genetic clusters occurring on either side of the St. Lawrence River, and did not detect any effect of the host genotype on their liver microbiome assemblage. Finally, we report, for the first time, the presence of B. burgdorferi in a small mammal host from the northern side of the St. Lawrence River, in support of models that have predicted the northern spread of Lyme disease in Canada.

Stress level, parasite load, and movement pattern in a small-mammal reservoir host for Lyme disease

Occurrence of Lyme disease has increased rapidly in Canada in the past 5 years. The emergence of Lyme disease coincides with the range expansion of the primary host, the white-footed mouse (Peromyscus leucopus (Rafinesque, 1818)), in the region. We evaluated the effects of stress level, parasite load, and forest-patch characteristics on P. leucopus movement pattern. We found negative relations between on the one hand the adrenal gland size, a proxy for stress level, and population density, and on the other hand, home-range area and movement rate of mouse individuals, suggesting that stressed mice cannot maintain a large home range. Population density was also related with excursion (outside the forest patch) and exploration (outside the home range) rates, either directly or through its effect on home-range area and movement rate. Finally, movement rate and excursion rate were lower in individuals infested with more black-legged ticks (Ixodes scapularis Say, 1821). Our results have implication for the mechanism of Lyme disease emergence in the region: individual hosts that carry more ticks and are thus more likely to be spreading the bacterium responsible for Lyme disease are dispersing less than tick-free individuals. Monitoring of Lyme disease should thus consider how the characteristics of host communities modulate the spread of the disease across the landscape.

Climate change and habitat fragmentation drive the occurrence of Borrelia burgdorferi, the agent of Lyme disease, at the northeastern limit of its distribution

Lyme borreliosis is rapidly emerging in Canada, and climate change is likely a key driver of the northern spread of the disease in North America. We used field and modeling approaches to predict the risk of occurrence of Borrelia burgdorferi, the bacteria causing Lyme disease in North America. We combined climatic and landscape variables to model the current and future (2050) potential distribution of the black-legged tick and the white-footed mouse at the northeastern range limit of Lyme disease and estimated a risk index for B. burgdorferi from these distributions. The risk index was mostly constrained by the distribution of the white-footed mouse, driven by winter climatic conditions. The next factor contributing to the risk index was the distribution of the black-legged tick, estimated from the temperature. Landscape variables such as forest habitat and connectivity contributed little to the risk index. We predict a further northern expansion of B. burgdorferi of approximately 250–500 km by 2050 – a rate of 3.5–11 km per year – and identify areas of rapid rise in the risk of occurrence of B. burgdorferi. Our results will improve understanding of the spread of Lyme disease and inform management strategies at the most northern limit of its distribution.