Multiomics Reveals Symbionts, Pathogens, and Tissue-Specific Microbiome of Blacklegged Ticks (Ixodes scapularis) from a Lyme Disease Hot Spot in Southeastern Ontario, Canada.

Ticks in the family Ixodidae are important vectors of zoonoses, including Lyme disease (LD), which is caused by spirochete bacteria from the Borreliella (Borrelia) burgdorferi sensu lato complex. The blacklegged tick (Ixodes scapularis) continues to expand across Canada, creating hot spots of elevated LD risk at the leading edge of its expanding range. Current efforts to understand the risk of pathogen transmission associated with I. scapularis in Canada focus primarily on targeted screens, while natural variation in the tick microbiome remains poorly understood. Using multiomics consisting of 16S metabarcoding and ribosome-depleted, whole-shotgun RNA transcriptome sequencing, we examined the microbial communities associated with adult I. scapularis (n = 32), sampled from four tissue types (whole tick, salivary glands, midgut, and viscera) and three geographical locations within a LD hot spot near Kingston, Ontario, Canada. The communities consisted of both endosymbiotic and known or potentially pathogenic microbes, including RNA viruses, bacteria, and a Babesia sp. intracellular parasite. We show that ß-diversity is significantly higher between the bacterial communities of individual tick salivary glands and midguts than that of whole ticks. Linear discriminant analysis effect size (LEfSe) determined that the three potentially pathogenic bacteria detected by V4 16S rRNA sequencing also differed among dissected tissues only, including a Borrelia strain from the B. burgdorferi sensu lato complex, Borrelia miyamotoi, and Anaplasma phagocytophilum. Importantly, we find coinfection of I. scapularis by multiple microbes, in contrast to diagnostic protocols for LD, which typically focus on infection from a single pathogen of interest (B. burgdorferi sensu stricto). IMPORTANCE As a vector of human health concern, blacklegged ticks (Ixodes scapularis) transmit pathogens that cause tick-borne diseases (TBDs), including Lyme disease (LD). Several hot spots of elevated LD risk have emerged across Canada as I. scapularis expands its range. Focusing on a hot spot in southeastern Ontario, we used high-throughput sequencing to characterize the microbiome of whole ticks and dissected salivary glands and midguts. Compared with whole ticks, salivary glands and midguts were more diverse and associated with distinct bacterial communities that are less dominated by Rickettsia endosymbiont bacteria and are enriched for pathogenic bacteria, including a B. burgdorferi sensu lato-associated Borrelia sp., Borrelia miyamotoi, and Anaplasma phagocytophilum. We also found evidence of coinfection of I. scapularis by multiple pathogens. Overall, our study highlights the challenges and opportunities associated with the surveillance of the microbiome of I. scapularis for pathogen detection using metabarcoding and metatranscriptome approaches.

Small-molecule metabolome identifies potential therapeutic targets against COVID-19.

Respiratory viruses are transmitted and acquired via the nasal mucosa, and thereby may influence the nasal metabolome composed of biochemical products produced by both host cells and microbes. Studies of the nasal metabolome demonstrate virus-specific changes that sometimes correlate with viral load and disease severity. Here, we evaluate the nasopharyngeal metabolome of COVID-19 infected individuals and report several small molecules that may be used as potential therapeutic targets. Specimens were tested by qRT-PCR with target primers for three viruses: Influenza A (INFA), respiratory syncytial virus (RSV), and SARS-CoV-2, along with unaffected controls. The nasopharyngeal metabolome was characterized using an LC-MS/MS-based screening kit capable of quantifying 141 analytes. A machine learning model identified 28 discriminating analytes and correctly categorized patients with a viral infection with an accuracy of 96% (R2 = 0.771, Q2 = 0.72). A second model identified 5 analytes to differentiate COVID19-infected patients from those with INFA or RSV with an accuracy of 85% (R2 = 0.442, Q2 = 0.301). Specifically, Lysophosphatidylcholines-a-C18:2 (LysoPCaC18:2) concentration was significantly increased in COVID19 patients (P < 0.0001), whereas beta-hydroxybutyric acid, Methionine sulfoxide, succinic acid, and carnosine concentrations were significantly decreased (P < 0.0001). This study demonstrates that COVID19 infection results in a unique nasopharyngeal metabolomic signature with carnosine and LysoPCaC18:2 as potential therapeutic targets.

Evidence for continent-wide convergent evolution and stasis throughout 150 y of a biological invasion.

The extent to which evolution can rescue a species from extinction, or facilitate range expansion, depends critically on the rate, duration, and geographical extent of the evolutionary response to natural selection. Adaptive evolution can occur quickly, but the duration and geographical extent of contemporary evolution in natural systems remain poorly studied. This is particularly true for species with large geographical ranges and for timescales that lie between "long-term" field experiments and the fossil record. Here, we introduce the Virtual Common Garden (VCG) to investigate phenotypic evolution in natural history collections while controlling for phenotypic plasticity in response to local growing conditions. Reconstructing 150 y of evolution in Lythrum salicaria (purple loosestrife) as it invaded North America, we analyze phenology measurements of 3,429 herbarium records, reconstruct growing conditions from more than 12 million local temperature records, and validate predictions across three common gardens spanning 10° of latitude. We find that phenological clines have evolved repeatedly throughout the range, during the first century of evolution. Thereafter, the rate of microevolution stalls, recapitulating macroevolutionary stasis observed in the fossil record. Our study demonstrates that preserved specimens are a critical resource for investigating limits to evolution in natural populations. Our results show how natural selection and trade-offs measured in field studies predict adaptive divergence observable in herbarium specimens over 15 decades at a continental scale.

Effects of species interactions on the potential for evolution at species' range limits.

Species' ranges are limited by both ecological and evolutionary constraints. While there is a growing appreciation that ecological constraints include interactions among species, like competition, we know relatively little about how interactions contribute to evolutionary constraints at species' niche and range limits. Building on concepts from community ecology and evolutionary biology, we review how biotic interactions can influence adaptation at range limits by impeding the demographic conditions that facilitate evolution (which we term a 'demographic pathway to adaptation'), and/or by imposing evolutionary trade-offs with the abiotic environment (a 'trade-offs pathway'). While theory for the former is well-developed, theory for the trade-offs pathway is not, and empirical evidence is scarce for both. Therefore, we develop a model to illustrate how fitness trade-offs along biotic and abiotic gradients could affect the potential for range expansion and niche evolution following ecological release. The model shows that which genotypes are favoured at species' range edges can depend strongly on the biotic context and the nature of fitness trade-offs. Experiments that characterize trade-offs and properly account for biotic context are needed to predict which species will expand their niche or range in response to environmental change. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Genome report: a draft genome of Alliaria petiolata (garlic mustard) as a model system for invasion genetics.

The emerging field of invasion genetics examines the genetic causes and consequences of biological invasions, but few study systems are available that integrate deep ecological knowledge with genomic tools. Here, we report on the de novo assembly and annotation of a genome for the biennial herb Alliaria petiolata (M. Bieb.) Cavara and Grande (Brassicaceae), which is widespread in Eurasia and invasive across much of temperate North America. Our goal was to sequence and annotate a genome to complement resources available from hundreds of published ecological studies, a global field survey, and hundreds of genetic lines maintained in Germany and Canada. We sequenced a genotype (EFCC3-3-20) collected from the native range near Venice, Italy, and sequenced paired-end and mate pair libraries at ∼70 × coverage. A de novo assembly resulted in a highly continuous draft genome (N50 = 121 Mb; L50 = 2) with 99.7% of the 1.1 Gb genome mapping to scaffolds of at least 50 Kb in length. A total of 64,770 predicted genes in the annotated genome include 99% of plant BUSCO genes and 98% of transcriptome reads. Consistent with previous reports of (auto)hexaploidy in western Europe, we found that almost one-third of BUSCO genes (390/1440) mapped to two or more scaffolds despite <2% genome-wide average heterozygosity. The continuity and gene space quality of our draft assembly will enable molecular and functional genomic studies of A. petiolata to address questions relevant to invasion genetics and conservation strategies.

Temporal Dynamics and Evolution of SARS-CoV-2 Demonstrate the Necessity of Ongoing Viral Genome Sequencing in Ontario, Canada.

Genome-wide variation in SARS-CoV-2 reveals evolution and transmission dynamics which are critical considerations for disease control and prevention decisions. Here, we review estimates of the genome-wide viral mutation rates, summarize current COVID-19 case load in the province of Ontario, Canada (5 January 2021), and analyze published SARS-CoV-2 genomes from Ontario (collected prior to 24 November 2020) to test for more infectious genetic variants or lineages. The reported mutation rate (∼10-6 nucleotide [nt]-1 cycle-1) for SARS-CoV-2 is typical for coronaviruses. Analysis of published SARS-CoV-2 genomes revealed that the G614 spike protein mutation has dominated infections in Ontario and that SARS-CoV-2 lineages present in Ontario have not differed significantly in their rate of spread. These results suggest that the SARS-CoV-2 population circulating in Ontario has not changed significantly to date. However, ongoing genome monitoring is essential for identification of new variants and lineages that may contribute to increased viral transmission.

Convergence Research for Emerging Zoonoses.

Global emergence of vector-borne and zoonotic diseases presents a rapidly growing 'wicked' problem. We outline the need for a transdisciplinary research program that is grounded in ecological and evolutionary theory but integrates fundamentally with research perspectives spanning the health, social, and natural sciences.

Phylogenomics reveals viral sources, transmission, and potential superinfection in early-stage COVID-19 patients in Ontario, Canada.

The emergence and rapid global spread of SARS-CoV-2 demonstrates the importance of infectious disease surveillance, particularly during the early stages. Viral genomes can provide key insights into transmission chains and pathogenicity. Nasopharyngeal swabs were obtained from thirty-two of the first SARS-CoV-2 positive cases (March 18-30) in Kingston Ontario, Canada. Viral genomes were sequenced using Ion Torrent (n = 24) and MinION (n = 27) sequencing platforms. SARS-CoV-2 genomes carried forty-six polymorphic sites including two missense and three synonymous variants in the spike protein gene. The D614G point mutation was the predominate viral strain in our cohort (92.6%). A heterozygous variant (C9994A) was detected by both sequencing platforms but filtered by the ARTIC network bioinformatic pipeline suggesting that heterozygous variants may be underreported in the SARS-CoV-2 literature. Phylogenetic analysis with 87,738 genomes in the GISAID database identified global origins and transmission events including multiple, international introductions as well as community spread. Reported travel history validated viral introduction and transmission inferred by phylogenetic analysis. Molecular epidemiology and evolutionary phylogenetics may complement contact tracing and help reconstruct transmission chains of emerging diseases. Earlier detection and screening in this way could improve the effectiveness of regional public health interventions to limit future pandemics.

Integrating morphological characters, molecular markers, and distribution patterns to assess the identity of Blepharis species from Jordan.

BACKGROUND: Blepharis constitutes an important part of the vegetation of the Jordanian arid and semi-arid regions, yet whether one or more species of this genus occurs in the Jordanian area is uncertain. We addressed this question by assessing morphological characters and testing Inter-Simple Sequence Repeat (ISSR) markers from three populations of Blepharis: two northern (lower slopes of Kufranjah valley and the Dead Sea region) and one southern (Wadi al Yutm). RESULTS: Shoots from randomly chosen Blepharis plants were harvested from each of the three populations for morphological and molecular analyses. In the northern populations, spikes were lax and bract width was significantly shorter than length of the longest lateral spine compared to the southern population. A multivariate linear discriminant analysis distinguished the northern populations from the southern one by internode length, bract width, longest lateral spine length, and bract width to spine length ratio. The ISSR analysis revealed that 44 markers across eight primers were polymorphic with major allele frequency of 83.6% and an average of 5.5 polymorphic markers per primer. The genetic resemblance among individuals ranged from 0.27 to 0.96. The three Blepharis populations were accordingly clustered into two distinct groups, similar to the analysis of morphological differences and corresponding with the "northern" and "southern" population designations. CONCLUSIONS: Our results strongly indicate the occurrence of two discrete Blepharis species in Jordan and reject the hypothesis that the genus is represented by only one species. We propose that the Blepharis species in Jordan are B. attenutata Napper (represented by the northern populations) and B. ciliaris (L.) B. L. Burtt (represented by the southern population). These findings are important for informing and revising floristic work within the region and an updated key has been included in our findings.

Understanding and monitoring the consequences of human impacts on intraspecific variation.

Intraspecific variation is a major component of biodiversity, yet it has received relatively little attention from governmental and nongovernmental organizations, especially with regard to conservation plans and the management of wild species. This omission is ill-advised because phenotypic and genetic variations within and among populations can have dramatic effects on ecological and evolutionary processes, including responses to environmental change, the maintenance of species diversity, and ecological stability and resilience. At the same time, environmental changes associated with many human activities, such as land use and climate change, have dramatic and often negative impacts on intraspecific variation. We argue for the need for local, regional, and global programs to monitor intraspecific genetic variation. We suggest that such monitoring should include two main strategies: (i) intensive monitoring of multiple types of genetic variation in selected species and (ii) broad-brush modeling for representative species for predicting changes in variation as a function of changes in population size and range extent. Overall, we call for collaborative efforts to initiate the urgently needed monitoring of intraspecific variation.

Invasions and extinctions through the looking glass of evolutionary ecology.

Invasive and endangered species reflect opposite ends of a spectrum of ecological success, yet they experience many similar eco-evolutionary challenges including demographic bottlenecks, hybridization and novel environments. Despite these similarities, important differences exist. Demographic bottlenecks are more transient in invasive species, which (i) maintains ecologically relevant genetic variation, (ii) reduces mutation load, and (iii) increases the efficiency of natural selection relative to genetic drift. Endangered species are less likely to benefit from admixture, which offsets mutation load but also reduces fitness when populations are locally adapted. Invading species generally experience more benign environments with fewer natural enemies, which increases fitness directly and also indirectly by masking inbreeding depression. Adaptive phenotypic plasticity can maintain fitness in novel environments but is more likely to evolve in invasive species encountering variable habitats and to be compromised by demographic factors in endangered species. Placed in an eco-evolutionary context, these differences affect the breadth of the ecological niche, which arises as an emergent property of antagonistic selection and genetic constraints. Comparative studies of invasions and extinctions that apply an eco-evolutionary perspective could provide new insights into the environmental and genetic basis of ecological success in novel environments and improve efforts to preserve global biodiversity.This article is part of the themed issue 'Human influences on evolution, and the ecological and societal consequences'.

Phenological shifts of native and invasive species under climate change: insights from the Boechera-Lythrum model.

Warmer and drier climates have shifted phenologies of many species. However, the magnitude and direction of phenological shifts vary widely among taxa, and it is often unclear when shifts are adaptive or how they affect long-term viability. Here, we model evolution of flowering phenology based on our long-term research of two species exhibiting opposite shifts in floral phenology: Lythrum salicaria, which is invasive in North America, and the sparse Rocky Mountain native Boechera stricta Genetic constraints are similar in both species, but differences in the timing of environmental conditions that favour growth lead to opposite phenological shifts under climate change. As temperatures increase, selection is predicted to favour earlier flowering in native B. stricta while reducing population viability, even if populations adapt rapidly to changing environmental conditions. By contrast, warming is predicted to favour delayed flowering in both native and introduced L. salicaria populations while increasing long-term viability. Relaxed selection from natural enemies in invasive L. salicaria is predicted to have little effect on flowering time but a large effect on reproductive fitness. Our approach highlights the importance of understanding ecological and genetic constraints to predict the ecological consequences of evolutionary responses to climate change on contemporary timescales.This article is part of the themed issue 'Human influences on evolution, and the ecological and societal consequences'.

Contemporary evolution during invasion: evidence for differentiation, natural selection, and local adaptation.

Biological invasions are 'natural' experiments that can improve our understanding of contemporary evolution. We evaluate evidence for population differentiation, natural selection and adaptive evolution of invading plants and animals at two nested spatial scales: (i) among introduced populations (ii) between native and introduced genotypes. Evolution during invasion is frequently inferred, but rarely confirmed as adaptive. In common garden studies, quantitative trait differentiation is only marginally lower (~3.5%) among introduced relative to native populations, despite genetic bottlenecks and shorter timescales (i.e. millennia vs. decades). However, differentiation between genotypes from the native vs. introduced range is less clear and confounded by nonrandom geographic sampling; simulations suggest this causes a high false-positive discovery rate (>50%) in geographically structured populations. Selection differentials (¦s¦) are stronger in introduced than in native species, although selection gradients (¦ß¦) are not, consistent with introduced species experiencing weaker genetic constraints. This could facilitate rapid adaptation, but evidence is limited. For example, rapid phenotypic evolution often manifests as geographical clines, but simulations demonstrate that nonadaptive trait clines can evolve frequently during colonization (~two-thirds of simulations). Additionally, QST-FST studies may often misrepresent the strength and form of natural selection acting during invasion. Instead, classic approaches in evolutionary ecology (e.g. selection analysis, reciprocal transplant, artificial selection) are necessary to determine the frequency of adaptive evolution during invasion and its influence on establishment, spread and impact of invasive species. These studies are rare but crucial for managing biological invasions in the context of global change.

Rapid adaptation to climate facilitates range expansion of an invasive plant.

Adaptation to climate, evolving over contemporary time scales, could facilitate rapid range expansion across environmental gradients. Here, we examine local adaptation along a climatic gradient in the North American invasive plant Lythrum salicaria. We show that the evolution of earlier flowering is adaptive at the northern invasion front where it increases fitness as much as, or more than, the effects of enemy release and the evolution of increased competitive ability. However, early flowering decreases investment in vegetative growth, which reduces fitness by a factor of 3 in southern environments where the North American invasion commenced. Our results demonstrate that local adaptation can evolve quickly during range expansion, overcoming environmental constraints on propagule production.

Genetic trade-offs and conditional neutrality contribute to local adaptation.

Divergent natural selection promotes local adaptation and can lead to reproductive isolation of populations in contrasting environments; however, the genetic basis of local adaptation remains largely unresolved in natural populations. Local adaptation might result from antagonistic pleiotropy, where alternate alleles are favoured in distinct habitats, and polymorphism is maintained by selection. Alternatively, under conditional neutrality some alleles may be favoured in one environment but neutral at other locations. Antagonistic pleiotropy maintains genetic variation across the landscape; however, there is a systematic bias against discovery of antagonistic pleiotropy because the fitness benefits of local alleles need to be significant in at least two environments. Here, we develop a generally applicable method to investigate polygenic local adaptation and identify loci that are the targets of selection. This approach evaluates allele frequency changes after selection at loci across the genome to distinguish antagonistic pleiotropy from conditional neutrality and deleterious variation. We investigate local adaptation at the qualitative trait loci (QTL) level in field experiments, in which we expose 177 F(6) recombinant inbred lines and parental lines of Boechera stricta (Brassicaceae) to their parental environments over two seasons. We demonstrate polygenic selection for native alleles in both environments, with 2.8% of the genome exhibiting antagonistic pleiotropy and 8% displaying conditional neutrality. Our study strongly supports antagonistic pleiotropy at one large-effect flowering phenology QTL (nFT): native homozygotes had significantly greater probabilities of flowering than foreign homozygotes in both parental environments. Such large-scale field studies are essential to elucidate the genetic basis of adaptation in natural populations.

Phenotypic plasticity and adaptive evolution contribute to advancing flowering phenology in response to climate change.

Anthropogenic climate change has already altered the timing of major life-history transitions, such as the initiation of reproduction. Both phenotypic plasticity and adaptive evolution can underlie rapid phenological shifts in response to climate change, but their relative contributions are poorly understood. Here, we combine a continuous 38 year field survey with quantitative genetic field experiments to assess adaptation in the context of climate change. We focused on Boechera stricta (Brassicaeae), a mustard native to the US Rocky Mountains. Flowering phenology advanced significantly from 1973 to 2011, and was strongly associated with warmer temperatures and earlier snowmelt dates. Strong directional selection favoured earlier flowering in contemporary environments (2010-2011). Climate change could drive this directional selection, and promote even earlier flowering as temperatures continue to increase. Our quantitative genetic analyses predict a response to selection of 0.2 to 0.5 days acceleration in flowering per generation, which could account for more than 20 per cent of the phenological change observed in the long-term dataset. However, the strength of directional selection and the predicted evolutionary response are likely much greater now than even 30 years ago because of rapidly changing climatic conditions. We predict that adaptation will likely be necessary for long-term in situ persistence in the context of climate change.

Origin, fate, and architecture of ecologically relevant genetic variation.

Recent advances in molecular genetics combined with field manipulations are yielding new insight into the origin, evolutionary fate, and genetic architecture of phenotypic variation in natural plant populations, with two surprising implications for the evolution of plant genomes. First, genetic loci exhibiting antagonistic pleiotropy across natural environments appear rare relative to loci that are adaptive in one or more environments and neutral elsewhere. These 'conditionally neutral' alleles should sweep to fixation when they arise, yet genome comparisons find little evidence for such selective sweeps. Second, genes under biotic selection tend to be of larger effect than genes under abiotic selection. Recent theory suggests this may be a consequence of high gene flow among populations under selection for local adaptation.

Population divergence along lines of genetic variance and covariance in the invasive plant Lythrum salicaria in eastern North America.

Evolution during biological invasion may occur over contemporary timescales, but the rate of evolutionary change may be inhibited by a lack of standing genetic variation for ecologically relevant traits and by fitness trade-offs among them. The extent to which these genetic constraints limit the evolution of local adaptation during biological invasion has rarely been examined. To investigate genetic constraints on life-history traits, we measured standing genetic variance and covariance in 20 populations of the invasive plant purple loosestrife (Lythrum salicaria) sampled along a latitudinal climatic gradient in eastern North America and grown under uniform conditions in a glasshouse. Genetic variances within and among populations were significant for all traits; however, strong intercorrelations among measurements of seedling growth rate, time to reproductive maturity and adult size suggested that fitness trade-offs have constrained population divergence. Evidence to support this hypothesis was obtained from the genetic variance-covariance matrix (G) and the matrix of (co)variance among population means (D), which were 79.8% (95% C.I. 77.7-82.9%) similar. These results suggest that population divergence during invasive spread of L. salicaria in eastern North America has been constrained by strong genetic correlations among life-history traits, despite large amounts of standing genetic variation for individual traits.

Evolutionary constraints on adaptive evolution during range expansion in an invasive plant.

Biological invasions may expose populations to strong selection for local adaptation along geographical gradients in climate. However, evolution during contemporary timescales can be constrained by low standing genetic variation and genetic correlations among life-history traits. We examined limits to local adaptation associated with northern migration of the invasive wetland plant purple loosestrife (Lythrum salicaria) using a selection model incorporating a trade-off between flowering time and size at reproduction, and common garden experiments of populations sampled along a latitudinal transect of approximately 1200 km in eastern North America. A strong trade-off between flowering time and size at reproduction caused early-flowering plants to be smaller with reduced seed production in northern populations. Northward spread was associated with a decline in genetic variance within populations and an increase in genetic skew for flowering time and size, with limited genetic variation for small, early-flowering genotypes. These patterns were predicted by our selection model of local adaptation to shorter growing seasons and were not consistent with expectations from non-adaptive processes. Reduced fecundity may limit population growth and rates of spread in northern populations. Identifying genetic constraints on key life-history traits can provide novel insights into invasion dynamics and the causes of range limits in introduced species.

Common garden comparisons of native and introduced plant populations: latitudinal clines can obscure evolutionary inferences.

Common garden studies are increasingly used to identify differences in phenotypic traits between native and introduced genotypes, often ignoring sources of among-population variation within each range. We re-analyzed data from 32 common garden studies of 28 plant species that tested for rapid evolution associated with biological invasion. Our goals were: (i) to identify patterns of phenotypic trait variation among populations within native and introduced ranges, and (ii) to explore the consequences of this variation for how differences between the ranges are interpreted. We combined life history and physiologic traits into a single principal component (PCALL) and also compared subsets of traits related to size, reproduction, and defense (PCSIZE, PCREP, and PCDEF, respectively). On average, introduced populations exhibited increased growth and reproduction compared to native conspecifics when latitude was not included in statistical models. However, significant correlations between PC-scores and latitude were detected in both the native and introduced ranges, indicating population differentiation along latitudinal gradients. When latitude was explicitly incorporated into statistical models as a covariate, it reduced the magnitude and reversed the direction of the effect for PCALL and PCSIZE. These results indicate that unrecognized geographic clines in phenotypic traits can confound inferences about the causes of evolutionary change in invasive plants.