Phylogeographic Structure of the White-Footed Mouse and the Deer Mouse, Two Lyme Disease Reservoir Hosts in Québec

Thermoregulatory trade-offs underlie the effects of warming summer temperatures on deer mice

Climate warming could challenge the ability of endotherms to thermoregulate and maintain normal body temperature (Tb), but the effects of warming summer temperatures on activity and thermoregulatory physiology in many small mammals remain poorly understood. We examined this issue in deer mice (Peromyscus maniculatus), an active nocturnal species. Mice were exposed in the lab to simulated seasonal warming, in which an environmentally realistic diel cycle of ambient temperature (Ta) was gradually warmed from spring conditions to summer conditions (controls were maintained in spring conditions). Activity (voluntary wheel running) and Tb (implanted bio-loggers) were measured throughout, and indices of thermoregulatory physiology (thermoneutral zone, thermogenic capacity) were assessed after exposure. In control mice, activity was almost entirely restricted to the night-time, and Tb fluctuated ∼1.7°C between daytime lows and night-time highs. Activity, body mass and food consumption were reduced and water consumption was increased in later stages of summer warming. This was accompanied by strong Tb dysregulation that culminated in a complete reversal of the diel pattern of Tb variation, with Tb reaching extreme highs (∼40°C) during daytime heat but extreme lows (∼34°C) at cooler night-time temperatures. Summer warming was also associated with reduced ability to generate body heat, as reflected by decreased thermogenic capacity and decreased mass and uncoupling protein (UCP1) content of brown adipose tissue. Our findings suggest that thermoregulatory trade-offs associated with daytime heat exposure can affect Tb and activity at cooler night-time temperatures, impacting the ability of nocturnal mammals to perform behaviours important for fitness in the wild.

An ecological approach to predict areas with established populations of Ixodes scapularis in Quebec, Canada

Public health management of Lyme disease (LD) is a dynamic challenge in Canada. Climate warming is driving the northward expansion of suitable habitat for the tick vector, Ixodes scapularis. Information about tick population establishment is used to inform the risk of LD but is challenged by sampling biases from surveillance data. Misclassifying areas as having no established tick population underestimates the LD risk classification. We used a logistic regression model at the municipal level to predict the probability of I. scapularis population establishment based on passive tick surveillance data during the period of 2010-2017 in southern Quebec. We tested for the effect of abiotic and biotic factors hypothesized to influence tick biology and ecology. Additional variables controlled for sampling biases in the passive surveillance data. In our final selected model, tick population establishment was positively associated with annual cumulative degree-days > 0°C, precipitation and deer density, and negatively associated with coniferous and mixed forest types. Sampling biases from passive tick surveillance were controlled for using municipal population size and public health instructions on tick submissions. The model performed well as indicated by an area under the curve (AUC) of 0.92, sensitivity of 86% and specificity of 81%. Our model enables prediction of I. scapularis population establishment in areas which lack data from passive tick surveillance and may improve the sensitivity of LD risk categorization in these areas. A more sensitive system of LD risk classification is important for increasing awareness and use of protective measures employed against ticks, and decreasing the morbidity associated with LD.

Emerging Tick-Borne Pathogens in Central Canada: Recent Detections of Babesia odocoilei and Rickettsia rickettsii

Background: The spread of emerging tick-borne pathogens has steadily increased in Canada with the widespread establishment of tick vectors and vertebrate hosts. At present, Borrelia burgdorferi, the bacterium causing Lyme disease, is the most common tick-borne pathogen in Canada and primarily transmitted by Ixodes scapularis. A low prevalence of other emerging tick-borne pathogens, such as Anaplasma phagocytophilum, Babesia species, Borrelia miyamotoi, and Francisella tularensis have also been detected through surveillance efforts in Canada. Although Rickettsia rickettsii has been historically detected in Haemaphysalis leporispalustris in Canada, the current prevalence and geographic extent of this pathogen is unknown. Material and Methods: In this study, we assessed the presence and prevalence of several emerging tick-borne pathogens in ticks and hosts collected through tick dragging and small mammal trapping in Central Canada. Results: Nested PCR testing detected three pathogen species in ticks, with Babesia odocoilei and B. burgdorferi in I. scapularis in addition to R. rickettsii in H. leporispalustris. Three pathogen species were detected in small mammals by nested PCR including B. odocoilei in Blarina brevicauda, Babesia microti in Peromyscus leucopus, and a Hepatozoon species in P. leucopus and Peromyscus maniculatus. B. burgdorferi and Babesia species were the pathogens most often detected in our samples, suggesting they are widely distributed across Central Canada. We also detected B. odocoilei and R. rickettsii beyond their known geographic distribution. Conclusions: Our results provide evidence that emerging tick-borne pathogens may be present outside defined risk areas identified by current surveillance efforts in Canada. As a result, emerging tick-borne pathogens introduced by the dispersal of infected ticks by migratory birds or maintained by hosts and vectors through cryptic transmission cycles may go undetected. More comprehensive testing including all tick life stages and additional tick-borne pathogens will help detect the spread and potential risk of emerging or re-emerging tick-borne pathogens for human and wildlife populations throughout Canada.

Effectiveness of immunization with multi-component bacterial immunomodulator in foals at 35th day of life

The aim of the study was to investigate the mechanisms leading to immunization through the use of a multicomponent bacterial immunomodulator and to find out the relationship between the TLR 4 receptor with selected parameters of innate immunity and to acquire immunity. The study was conducted on 18 Polish Pony Horses foals divided into two study groups: control (n = 9) and experimental (n = 9). Foals from the experimental group received intramuscular duplicate injection of 5 ml of multi-component bacterial immunomodular at 35 and 40 days of age. RNA isolated from venous blood was used to evaluate the expression of TLR4 genes using RT-PCR. Concentration of Il-6, IL-10, IgM and IgG2 was determined by the ELISA method in blood plasma. Immunostimulation had a significant impact on the level of genes expression for TLR4 expression and IL-6 concentration. No effect of stimulation on IgM and IgG2 concentrations was found. The expression of TLR4 genes as well as the levels of interleukins could be modulated by stimulation with a pharmacological agent multi-component bacterial immunomodulator. The experiment demonstrated a strong positive correlation between TLR4 gene expression and IL-6 concentration and between TLR4 gene expression and IgM concentration.

Borrelia burgdorferi and Borrelia miyamotoi in Atlantic Canadian wildlife

Borrelia burgdorferi and Borrelia miyamotoi are tick-vectored zoonotic pathogens maintained in wildlife species. Tick populations are establishing in new areas globally in response to climate change and other factors. New Brunswick is a Canadian maritime province at the advancing front of tick population establishment and has seen increasing numbers of ticks carrying B. burgdorferi, and more recently B. miyamotoi. Further, it is part of a region of Atlantic Canada with wildlife species composition differing from much of continental North America and little information exists as to the presence and frequency of infection of Borrelia spp. in wildlife in this region. We used a citizen science approach to collect a wide range of animals including migratory birds, medium-sized mammals, and small mammals. In total we tested 339 animals representing 20 species for the presence of B. burgdorferi and B. miyamotoi. We have developed new nested PCR primers and a protocol with excellent specificity for detecting both of these Borrelia species, both single and double infections, in tissues and organs of various wildlife species. The positive animals were primarily small non-migratory mammals, approximately twice as many were infected with B. burgdorferi than B. miyamotoi and one animal was found infected with both. In addition to established reservoir species, the jumping mouse (Napaeozapus insignis) was found frequently infected; this species had the highest infection prevalence for both B. burgdorferi and B. miyamotoi and has not previously been identified as an important carrier for either Borrelia species. Comprehensive testing of tissues found that all instances of B. burgdorferi infection were limited to one tissue within the host, whereas two of the five B. miyamotoi infections were diffuse and found in multiple systems. In the one coinfected specimen, two fetuses were also recovered and found infected with B. miyamotoi. This presumptive transplacental transmission suggests that vertical transmission in mammals is possible. This finding implies that B. miyamotoi could rapidly spread into wildlife populations, as well as having potential human health implications.

Ixodes scapularis (Ixodida: Ixodidae) Parasitizing an Unlikely Host: Big Brown Bats, Eptesicus fuscus (Chiroptera: Vespertilionidae), in New York State, USA

Ixodes scapularis Say is a three-host tick that has been recorded feeding on over 150 different species of terrestrial vertebrates (mammals, birds, and reptiles). This tick is found throughout the northeastern, coastal southeastern, and upper midwestern United States and is considered the most significant vector of tick-borne pathogens to humans in North America. Despite its ubiquity and broad host range, I. scapularis previously has not been reported feeding on bats (Chiroptera). However, during 2019 and 2020, larvae and nymphs of I. scapularis were recovered from big brown bats, Eptesicus fuscus (Palisot de Beauvois), at four locations in rural New York State, USA. All Ixodes infested bats were injured and found on the ground; therefore, parasitism by I. scapularis was likely opportunistic. Nonetheless, the large number of pathogens known to be associated with bats and the frequency with which I. scapularis bites people suggest that this host-tick relationship is of at least potential epidemiological significance.

Modelling Spatiotemporal Patterns of Lyme Disease Emergence in Québec

Lyme disease is a growing public health problem in Québec. Its emergence over the last decade is caused by environmental and anthropological factors that favour the survival of Ixodes scapularis, the vector of Lyme disease transmission. The objective of this study was to estimate the speed and direction of human Lyme disease emergence in Québec and to identify spatiotemporal risk patterns. A surface trend analysis was conducted to estimate the speed and direction of its emergence based upon the first detected case of Lyme disease in each municipality in Québec since 2004. A cluster analysis was also conducted to identify at-risk regions across space and time. These analyses were reproduced for the date of disease onset and date of notification for each case of Lyme disease. It was estimated that Lyme disease is spreading northward in Québec at a speed varying between 18 and 32 km/year according to the date of notification and the date of disease onset, respectively. A significantly high risk of disease was found in seven clusters identified in the south-west of Québec in the sociosanitary regions of Montérégie and Estrie. The results obtained in this study improve our understanding of the spatiotemporal patterns of Lyme disease in Québec, which can be used for proactive, targeted interventions by public and clinical health authorities.

Surveillance of Ticks and Tick-Borne Pathogens in Suburban Natural Habitats of Central Maryland

Lyme and other tick-borne diseases are increasing in the eastern United States and there is a lack of research on integrated strategies to control tick vectors. Here we present results of a study on tick-borne pathogens detected from tick vectors and rodent reservoirs from an ongoing 5-yr tick suppression study in the Lyme disease-endemic state of Maryland, where human-biting tick species, including Ixodes scapularis Say (Acari: Ixodidae) (the primary vector of Lyme disease spirochetes), are abundant. During the 2017 tick season, we collected 207 questing ticks and 602 ticks recovered from 327 mice (Peromyscus spp. (Rodentia: Cricetidae)), together with blood and ear tissue from the mice, at seven suburban parks in Howard County. Ticks were selectively tested for the presence of the causative agents of Lyme disease (Borrelia burgdorferi sensu lato [s.l.]), anaplasmosis (Anaplasma phagocytophilum), babesiosis (Babesia microti), ehrlichiosis (Ehrlichia ewingii, Ehrlichia chaffeensis, and 'Panola Mountain' Ehrlichia) and spotted fever group rickettsiosis (Rickettsia spp.). Peromyscus ear tissue and blood samples were tested for Bo. burgdorferi sensu stricto (s.s), A. phagocytophilum, Ba. microti, and Borrelia miyamotoi. We found 13.6% (15/110) of questing I. scapularis nymphs to be Bo. burgdorferi s.l. positive and 1.8% (2/110) were A. phagocytophilum positive among all sites. Borrelia burgdorferi s.s. was found in 71.1% (54/76) of I. scapularis nymphs removed from mice and 58.8% (194/330) of captured mice. Results from study on tick abundance and pathogen infection status in questing ticks, rodent reservoirs, and ticks feeding on Peromyscus spp. will aid efficacy evaluation of the integrated tick management measures being implemented.

Long-term stress level in a small mammal species undergoing range expansion

Rapid climate change is currently altering species distribution ranges. Evaluating the long-term stress level in wild species undergoing range expansion may help better understanding how species cope with the changing environment. Here, we focused on the white-footed mouse (Peromyscus leucopus), a widespread small mammal species in North-America whose distribution range is rapidly shifting northward. We evaluated long-term stress level in several populations of P. leucopus in Quebec (Canada), from the northern edge of the species distribution to more core populations in Southern Quebec. We first tested the hypothesis that populations at the range margin are under higher stress than more established populations in the southern region of our study area. We then compared four measures of long-term stress level to evaluate the congruence between these commonly used methods. We did not detect any significant geographical trend in stress level across our study populations of P. leucopus. Most notably, we found no clear congruence between the four measures of stress level we used, and conclude that these four commonly used methods are not equivalent, thereby not comparable across studies.

Genomic Signatures of Selection along a Climatic Gradient in the Northern Range Margin of the White-Footed Mouse (Peromyscus leucopus)

Identifying genetic variation involved in thermal adaptation is likely to yield insights into how species adapt to different climates. Physiological and behavioral responses associated with overwintering (e.g., torpor) are thought to serve important functions in climate adaptation. In this study, we use 2 isolated Peromyscus leucopus lineages on the northern margin of the species range to identify single nucleotide polymorphisms (SNPs) showing a strong environmental association and test for evidence of parallel evolution. We found signatures of clinal selection in each lineage, but evidence of parallelism was limited, with only 2 SNPs showing parallel allele frequencies across transects. These parallel SNPs map to a gene involved in protection against iron-dependent oxidative stress (Fxn) and to a gene with unknown function but containing a forkhead-associated domain (Fhad1). Furthermore, within transects, we find significant clinal patterns in genes enriched for functions associated with glycogen homeostasis, synaptic function, intracellular Ca2+ balance, H3 histone modification, as well as the G2/M transition of cell division. Our results are consistent with recent literature on the cellular and molecular basis of climate adaptation in small mammals and provide candidate genomic regions for further study.

A new mouse of the Peromyscus maniculatus species complex (Cricetidae) from the highlands of central Mexico

The Peromyscus maniculatus species complex is one of the most widespread group of small mammals in North America. However, the taxonomy and phylogenetic relationships among its constituent taxa remain unclear. As part of a revision of Peromyscus specimens from the highlands of the Trans-Mexican Volcanic Belt in central Mexico, we identified five individuals collected in 1968 that differed externally from other Peromyscus specimens, although morphologically similar to P. labecula and P. melanotis, both latter in the P. maniculatus species complex. Based on cranial measurements and mitochondrial DNA sequences, we aimed to more accurately determine the phylogenetic relationships and the taxonomic status of these individuals. Molecular phylogenetic analyses showed that the specimens formed a monophyletic clade sister to the P. maniculatus species complex. Pairwise genetic distances between those specimens and other species within the P. maniculatus species complex were greater than 7.91%. In addition, morphological analyses clearly distinguished the test specimens from P. melanotis and P. labecula. Based on the results of our molecular and morphological analyses, we conclude that these specimens represent an undescribed species of the P. maniculatus species complex, which we describe herein.

Conservation of the genome-wide recombination rate in white-footed mice

Despite being linked to the fundamental processes of chromosome segregation and offspring diversification, meiotic recombination rates vary within and between species. Recent years have seen progress in quantifying recombination rate evolution across multiple temporal and genomic scales. Nevertheless, the level of variation in recombination rate within wild populations-a key determinant of evolution in this trait-remains poorly documented on the genomic scale. To address this notable gap, we used immunofluorescent cytology to quantify genome-wide recombination rates in males from a wild population of the white-footed mouse, Peromyscus leucopus. For comparison, we measured recombination rates in a second population of male P. leucopus raised in the laboratory and in male deer mice from the subspecies Peromyscus maniculatus bairdii. Although we found differences between individuals in the genome-wide recombination rate, levels of variation were low-within populations, between populations, and between species. Quantification of synaptonemal complex length and crossover positions along chromosome 1 using a novel automated approach also revealed conservation in broad-scale crossover patterning, including strong crossover interference. We propose stabilizing selection targeting recombination or correlated processes as the explanation for these patterns.

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.

Infection resistance and tolerance in Peromyscus spp., natural reservoirs of microbes that are virulent for humans

The widely-distributed North American species Peromyscus leucopus and P. maniculatus of cricetine rodents are, between them, important natural reservoirs for several zoonotic diseases of humans: Lyme disease, human granulocytic anaplasmosis, babesiosis, erhlichiosis, hard tickborne relapsing fever, Powassan virus encephalitis, hantavirus pulmonary syndrome, and plague. While these infections are frequently disabling and sometimes fatal for humans, the peromyscines display little pathology and apparently suffer few consequences, even when prevalence of persistent infection in a population is high. While these Peromyscus spp. are unable to clear some of the infections, they appear to have partial resistance, which limits the burden of the pathogen. In addition, they display traits of infection tolerance, which reduces the damage of the infection. Research on these complementary resistance and tolerance phenomena in Peromyscus has relevance both for disease control measures targeting natural reservoirs and for understanding the mechanisms of the comparatively greater sickness of many humans with these and other infections.