For several decades, the Midwestern USA has been impacted by blacklegged tick (Ixodes scapularis) range expansion, which, as the main vector of the Lyme disease-causing bacterium Borrelia burgdorferi, is linked to a regional increase in Lyme disease incidence. Earlier studies of genetic differentiation of blacklegged ticks have not tested detailed hypotheses about range expansion in the Midwest, despite the importance of this topic to public health. We addressed this gap by investigating the origin and environmental factors that influenced blacklegged tick establishment and spread in the Midwestern region. By analyzing fine-scale spatial population genomic data, we find low genetic differentiation consistent with the known recent range expansion. However, within Wisconsin, blacklegged ticks have unique genetic ancestries that differ from other Midwestern regions, suggesting multiple origins. Our data provide evidence for blacklegged tick sources in northern Wisconsin contributing to the recent expansion. In addition, we find a distinctive mixture of ancestry along the Mississippi River in southwestern Wisconsin and in Indiana, which was previously not identified. The most recently invaded populations in Michigan exhibit sharp genetic divergence from Wisconsin and Indiana samples despite their proximity, warranting further examination of their genetic origin and expansion processes. Lastly, landscape factors contribute to significant reductions in gene flow, potentially limiting genetic exchange and disease transmission within Midwestern states. This new knowledge of blacklegged tick range expansion processes can improve vector surveillance, pest management, and public health related to tick-borne disease risks.

Tick vectors and tick-borne disease are increasingly impacting human populations globally. An important challenge is to understand tick movement patterns, as this information can be used to improve management and predictive modeling of tick population dynamics. Evolutionary analysis of genetic divergence, gene flow, and local adaptation provides insight on movement patterns at large spatiotemporal scales. We develop low coverage, whole genome resequencing data for 92 samples representing range-wide variation in blacklegged ticks, Ixodes scapularis, across the U.S. Through analysis of population genomic data, we find that tick populations are structured geographically, with gradual isolation by distance separating three population clusters in the northern U.S., southeastern U.S., and a unique cluster represented by a sample from Tennessee. Populations in the northern U.S. underwent population contractions during the last glacial period and diverged from southern populations at least 50 thousand years ago. Genome scans of selection provide strong evidence of local adaptation at genes responding to host defenses, blood-feeding, and environmental variation. In addition, we explore the potential of low coverage genome sequencing of whole-tick samples for documenting the diversity of microbial pathogens. Metagenomic analyses recover important tick-borne pathogens and their geographical variation, including the higher prevalence of Borrelia burgdorferi in northern populations, geographic strain variation in Rickettsia species, and the rare occurrence of B. miyamotoi and Anaplasma phagocytophilum. The combination of restricted pathogen distribution, isolation by distance, and local adaptation in blacklegged ticks demonstrates that gene flow, including recent expansion, is limited to geographical scales of a few hundred kilometers.