Effects of Female-Specific Selection for Reproductive Investment on Male Fertility Traits.

Despite sharing an autosomal genome, the often divergent reproductive strategies of males and females cause selection to act in a sex-specific manner. Selection acting on one sex can have negative, positive, or neutral fitness consequences on the opposite sex. Here we test how female-limited selection on reproductive investment in Japanese quail (Coturnix japonica) affects male fertility-related traits. Despite there being no difference in the size of males' testes from lines selected for high female reproductive investment (H-line) or low female reproductive investment (L-line), in both lines, the left testis had a greater volume of sperm-producing tissue. Since H-line females have a larger left-side restricted oviduct, this suggests a positive genetic correlation between male and female gonad function, and that internal testis structure is a target of sexual selection. However, despite H-line males having previously been found to have greater fertilisation success in a competitive scenario, we found little evidence of a difference between the lines in sperm number, motility, velocity, length, or the number of sperm that reached the ova. Pre-copulatory cues and/or the role of seminal fluid in sperm motility may thus be more likely to contribute to the H-line male fertilisation advantage in this species.

The great tit HapMap project: A continental-scale analysis of genomic variation in a songbird.

A major aim of evolutionary biology is to understand why patterns of genomic diversity vary within taxa and space. Large-scale genomic studies of widespread species are useful for studying how environment and demography shape patterns of genomic divergence. Here, we describe one of the most geographically comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 29 populations, spanning ~30 degrees of latitude and 40 degrees of longitude - almost the entire geographical range of the European subspecies. Genome-wide variation was consistent with a recent colonisation across Europe from a South-East European refugium, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear 'islands of differentiation', even among populations with very low levels of genome-wide differentiation. Low local recombination rates were a strong predictor of high local genomic differentiation (FST), especially in island and peripheral mainland populations, suggesting that the interplay between genetic drift and recombination causes highly heterogeneous differentiation landscapes. We also detected genomic outlier regions that were confined to one or more peripheral great tit populations, probably as a result of recent directional selection at the species' range edges. Haplotype-based measures of selection were related to recombination rate, albeit less strongly, and highlighted population-specific sweeps that likely resulted from positive selection. Our study highlights how comprehensive screens of genomic variation in wild organisms can provide unique insights into spatio-temporal evolutionary dynamics.

The association between age and telomere length is age-dependent: Evidence for a threshold model of telomere length maintenance.

Telomere length and dynamics are commonly used biomarkers of somatic state, yet the role of telomeres underlying the aging process is still debated. Indeed, to date, empirical evidence for an association between age and telomere length is mixed. Here, we test if the age-dependency of the association between age and telomere length can provide a potential explanation for the reported inconsistencies across studies. To this end, we quantified telomere length by telomere restriction fragment analysis in two groups of Japanese quail (Coturnix japonica) that differed in their age distribution. One group consisted of young adults only, whereas the second group consisted of adults across a wide range of ages. In the young adults group, there was a highly significant negative association between telomere length and age, whereas no association between age and telomere length was found in the all-ages adults group. This difference between groups was not due to telomere length-dependent selective disappearance. Our results shows that the association between telomere length and age is age-dependent and suggest that the costs and benefits associated with telomere maintenance are dynamic across an individual's life course.

Rapid decline of prenatal maternal effects with age is independent of postnatal environment in a precocial bird.

Maternal effects are an important source of phenotypic variation with potentially large fitness consequences, but how their importance varies with the quality of the environment across an individual's ontogeny is poorly understood. We bred Japanese quail (Coturnix japonica) of known pedigree and experimentally manipulated the quality of the offspring diet, to estimate the importance of prenatal maternal effects in shaping variation in body mass from hatching to adulthood. Maternal genetic effects on body mass at hatching were strong, and largely caused by variation in egg mass, but their importance rapidly declined with age. Whereas there was a large effect of diet on growth, this did not affect the decline of maternal effects variance. The importance of additive genetic and residual variance increased with age, with the latter being considerably larger in the poor diet treatment. Hence, we found no evidence for prenatal maternal effect by postnatal environment interactions and that prenatal maternal effects are rapidly replaced by direct additive genetic and residual effects when offspring start to develop outside the egg. Thereby these results shed new light on the dynamics of the role of maternal versus offspring genes across ontogeny and environments.

Increased male-induced harm in response to female-limited selection: interactive effects between intra- and interlocus sexual conflict?

Interlocus sexual conflict (IRSC) occurs because of shared interactions that have opposite effects on male and female fitness. Typically, it is assumed that loci involved in IRSC have sex-limited expression and are thus not directly affected by selective pressures acting on the other sex. However, if loci involved in IRSC have pleiotropic effects in the other sex, intersexual selection can shape the evolutionary dynamics of conflict escalation and resolution, as well as the evolution of reproductive traits linked to IRSC loci, and vice versa. Here we used an artificial selection approach in Japanese quail (Coturnix japonica) to test if female-limited selection on reproductive investment affects the amount of harm caused by males during mating. We found that males originating from lines selected for high female reproductive investment caused more oxidative damage in the female reproductive tract than males originating from lines selected for low female reproductive investment. This male-induced damage was specific to the oviduct and not found in other female tissues, suggesting that it was ejaculate-mediated. Our results suggest that intersexual selection shapes the evolution of IRSC and that male-induced harm may contribute to the maintenance of variation in female reproductive investment.

Ixodes ricinus density, Borrelia prevalence and the density of infected nymphs along an urban-rural gradient in southern England.

Ticks are found across a range of habitats, with woodland being particularly important for high densities and prevalence of Borrelia infection. Assessments of risk in urban woodland can be difficult if there are low densities and small sample sizes for Borrelia prevalence estimates. This study targeted six urban woodlands with established tick populations, as well as six woodlands in peri-urban zones and six woodlands in rural zones in and around the cities of Bath and Southampton, in the South of England. Nymph densities were estimated, and 100 nymphs were tested from each of the 18 woodlands studied. Ixodes ricinus ticks were found in all woodlands surveyed, and overall density of nymphs (DON) per 100 m2 was 18.17 in urban woodlands, 26.0 in peri-urban woodlands and 17.67 in rural woodlands. Out of 600 nymphs tested across urban woodlands, 10.3% were infected with Borrelia. The same proportion of nymphs collected in rural woodlands were positive for Borrelia. In peri-urban woodlands, 10.8% of nymphs tested positive. Across both cities combined, density of infected nymphs (DIN) was 2.73 per 100 m2 in peri-urban woodland, 1.87 per 100 m2 in urban woodland and 1.82 per 100 m2 in rural woodland. Overall, DON, Borrelia prevalence and DIN did not differ significantly along an urban-rural gradient. This suggests the risk of Lyme borreliosis transmission could be similar, or perhaps even elevated in urban woodland if there is higher public footfall, subsequent contact with ticks and less awareness of the risks. This is particularly important from a public health perspective, as Borrelia garinii dominated across the gradient and this genospecies is linked to neuroborreliosis.

Impact of green space connectivity on urban tick presence, density and Borrelia infected ticks in different habitats and seasons in three cities in southern England.

Understanding the effects of local habitat and wider landscape connectivity factors on tick presence, nymph density and Borrelia species (spp.) prevalence in the tick population is important for identifying the public health risk from Lyme borreliosis. This multi-city study collected data in three southern England cities (Bath, Bristol, and Southampton) during spring, summer, and autumn in 2017. Focusing specifically on urban green space used for recreation which were clearly in urbanised areas, 72 locations were sampled. Additionally, geospatial datasets on urban green space coverage within 250 m and 1 km of sampling points, as well as distance to woodland were incorporated into statistical models. Distance to woodland was negatively associated with tick presence and nymph density, particularly during spring and summer. Furthermore, we observed an interaction effect between habitat and season for tick presence and nymph density, with woodland habitat having greater tick presence and nymph density during spring. Borrelia spp. infected Ixodes ricinus were found in woodland, woodland edge and under canopy habitats in Bath and Southampton. Overall Borrelia spp. prevalence in nymphs was 2.8%, similar to wider UK studies assessing prevalence in Ixodes ricinus in rural areas. Bird-related Borrelia genospecies dominated across sites, suggesting bird reservoir hosts may be important in urban green space settings for feeding and infecting ticks. Whilst overall density of infected nymphs across the three cities was low (0.03 per 100 m2), risk should be further investigated by incorporating data on tick bites acquired in urban settings, and subsequent Lyme borreliosis transmission.

Tick range expansion to higher elevations: does Borrelia burgdorferi sensu lato facilitate the colonisation of marginal habitats?

BACKGROUND: Parasites can alter host and vector phenotype and thereby affect ecological processes in natural populations. Laboratory studies have suggested that Borrelia burgdorferi sensu lato, the causative agent of human Lyme borreliosis, may induce physiological and behavioural alterations in its main tick vector in Europe, Ixodes ricinus, which increase the tick's mobility and survival under challenging conditions. These phenotypic alterations may allow I. ricinus to colonise marginal habitats ('facilitation hypothesis'), thereby fuelling the ongoing range expansion of I. ricinus towards higher elevations and latitudes induced by climate change. To explore the potential for such an effect under natural conditions, we studied the prevalence of B. burgdorferi s.l. in questing I. ricinus and its variation with elevation in the Swiss Alps. RESULTS: We screened for B. burgdorferi s.l. infection in questing nymphs of I. ricinus (N = 411) from 15 sites between 528 and 1774 m.a.s.l to test if B. burgdorferi s.l. prevalence is higher at high elevations (i.e. in marginal habitats). Opposite of what is predicted under the facilitation hypothesis, we found that B. burgdorferi s.l. prevalence in I. ricinus nymphs decreased with increasing elevation and that Borrelia prevalence was 12.6% lower in I. ricinus nymphs collected at the range margin compared to nymphs in the core range. But there was no association between Borrelia prevalence and elevation within the core range of I. ricinus. Therefore the observed pattern was more consistent with a sudden decrease in Borrelia prevalence above a certain elevation, rather than a gradual decline with increasing elevation across the entire tick range. CONCLUSIONS: In conclusion, we found no evidence that B. burgdorferi s.l.-induced alterations of I. ricinus phenotype observed in laboratory studies facilitate the colonisation of marginal habitats in the wild. Rather, ticks in marginal habitats are substantially less likely to harbour the pathogen. These findings have implications for a better understanding of eco-evolutionary processes in natural host-parasite systems, as well as the assessment of Lyme borreliosis risk in regions where I. ricinus is newly emerging.

Reproductive Strategies Affect Telomere Dynamics across the Life Course.

AbstractBecause parental care has a heritable basis, the benefits of receiving increased parental provisioning early in life are genetically linked to the costs of providing increased parental provisioning at adulthood. Reproductive strategies thus result in distinct cost-benefit syndromes across the life course that may shape individual health and aging trajectories. Here we used an artificial selection approach in Japanese quail (Coturnix japonica) to test how reproductive strategies affect telomere length, a biomarker of somatic state, at different life stages. We show that males but not females from lines selected for low maternal investment (i.e., developing in a relatively small egg) had shorter telomeres at birth. These patterns were still weakly present at the end of the juvenile growth period. In contrast, significantly shorter telomeres were found in reproductively active adult birds from the high-investment lines, suggesting that telomere attrition was accelerated in these individuals once they had become reproductively active. Our study shows that reproductive strategies differentially affect telomere dynamics across the life course, highlighting the role of cross-generational constraints in shaping individual aging trajectories.

MHC class II genotype-by-pathogen genotype interaction for infection prevalence in a natural rodent-Borrelia system.

MHC genes are extraordinarily polymorphic in most taxa. Host-pathogen coevolution driven by negative frequency-dependent selection (NFDS) is one of the main hypotheses for the maintenance of such immunogenetic variation. Here, we test a critical but rarely tested assumption of this hypothesis-that MHC alleles affect resistance/susceptibility to a pathogen in a strain-specific way, that is, there is a host genotype-by-pathogen genotype interaction. In a field study of bank voles naturally infected with the tick-transmitted bacterium Borrelia afzelii, we tested for MHC class II (DQB) genotype-by-B. afzelii strain interactions for infection prevalence between 10 DQB alleles and seven strains. One allele (DQB*37) showed an interaction, such that voles carrying DQB*37 had higher prevalence of two strains and lower prevalence of one strain than individuals without the allele. These findings were corroborated by analyses of strain composition of infections, which revealed an effect of DQB*37 in the form of lower ß diversity among infections in voles carrying the allele. Taken together, these results provide rare support at the molecular genetic level for a key assumption of models of antagonistic coevolution through NFDS.

Maternally transferred thyroid hormones and life-history variation in birds.

In vertebrates, thyroid hormones (THs) play an important role in the regulation of growth, development, metabolism, photoperiodic responses and migration. Maternally transferred THs are important for normal early phase embryonic development when embryos are not able to produce endogenous THs. Previous studies have shown that variation in maternal THs within the physiological range can influence offspring phenotype. Given the essential functions of maternal THs in development and metabolism, THs may be a mediator of life-history variation across species. We tested the hypothesis that differences in life histories are associated with differences in maternal TH transfer across species. Using birds as a model, we specifically tested whether maternally transferred yolk THs covary with migratory status, developmental mode and traits related to pace-of-life (e.g. basal metabolic rate, maximum life span). We collected un-incubated eggs (n = 1-21 eggs per species, median = 7) from 34 wild and captive bird species across 17 families and six orders to measure yolk THs [both triiodothyronine (T3) and thyroxine (T4)], compiled life-history trait data from the literature and used Bayesian phylogenetic mixed models to test our hypotheses. Our models indicated that both concentrations and total amounts of the two main forms of THs (T3 and T4) were higher in the eggs of migratory species compared to resident species, and total amounts were higher in the eggs of precocial species, which have longer prenatal developmental periods, than in those of altricial species. However, maternal yolk THs did not show clear associations with pace-of-life-related traits, such as fecundity, basal metabolic rate or maximum life span. We quantified interspecific variation in maternal yolk THs in birds, and our findings suggest higher maternal TH transfer is associated with the precocial mode of development and migratory status. Whether maternal THs represent a part of the mechanism underlying the evolution of precocial development and migration or a consequence of such life histories is currently unclear. We therefore encourage further studies to explore the physiological mechanisms and evolutionary processes underlying these patterns.

Intralocus conflicts associated with a supergene.

Chromosomal inversions frequently underlie major phenotypic variation maintained by divergent selection within and between sexes. Here we examine whether and how intralocus conflicts contribute to balancing selection stabilizing an autosomal inversion polymorphism in the ruff Calidris pugnax. In this lekking shorebird, three male mating morphs (Independents, Satellites and Faeders) are controlled by an inversion-based supergene. We show that in a captive population, Faeder females, who are smaller and whose inversion haplotype has not undergone recombination, have lower average reproductive success in terms of laying rate, egg size, and offspring survival than Independent females, who lack the inversion. Satellite females, who carry a recombined inversion haplotype and have intermediate body size, more closely resemble Independent than Faeder females in reproductive performance. We inferred that the lower reproductive output of Faeder females is most likely balanced by higher than average reproductive success of individual Faeder males. These findings suggest that intralocus conflicts may play a major role in the evolution and maintenance of supergene variants.

Questing Ixodes ricinus ticks and Borrelia spp. in urban green space across Europe: A review.

For more than three decades, it has been recognized that Ixodes ricinus ticks occur in urban green space in Europe and that they harbour multiple pathogens linked to both human and animal diseases. Urban green space use for health and well-being, climate mitigation or biodiversity goals is promoted, often without consideration for the potential impact on tick encounters or tick-borne disease outcomes. This review synthesizes the results of over 100 publications on questing I. ricinus and Borrelia spp. infections in ticks in urban green space in 24 European countries. It presents data on several risk indicators for Lyme borreliosis and highlights key research gaps and recommendations for future studies. Across Europe, mean density of I. ricinus in urban green space was 6.9 (range; 0.1-28.8) per 100 m2 and mean Borrelia prevalence was 17.3% (range; 3.1%-38.1%). Similar density estimates were obtained for nymphs, which had a Borrelia prevalence of 14.2% (range; 0.5%-86.7%). Few studies provided data on both questing nymph density and Borrelia prevalence, but those that did found an average of 1.7 (range; 0-5.6) Borrelia-infected nymphs per 100 m2 of urban green space. Although a wide range of genospecies were reported, Borrelia afzelii was the most common in most parts of Europe, except for England where B. garinii was more common. The emerging pathogen Borrelia miyamotoi was also found in several countries, but with a much lower prevalence (1.5%). Our review highlights that I. ricinus and tick-borne Borrelia pathogens are found in a wide range of urban green space habitats and across several seasons. The impact of human exposure to I. ricinus and subsequent Lyme borreliosis incidence in urban green space has not been quantified. There is also a need to standardize sampling protocols to generate better baseline data for the density of ticks and Borrelia prevalence in urban areas.

Intergenerational Costs of Oxidative Stress: Reduced Fitness in Daughters of Mothers That Experienced High Levels of Oxidative Damage during Reproduction.

AbstractParental condition transfer effects occur when the parents' physiological state during reproduction affects offspring performance. Oxidative damage may mediate such effects, yet evidence that oxidative damage experienced by parents during reproduction negatively affects offspring fitness is scarce and limited to early life stages. We show in Japanese quail (Coturnix japonica) that maternal levels of oxidative damage, measured during reproduction, negatively predict the number of offspring produced by daughters. This maternal effect on daughters' reproductive success was mediated by an effect on hatching success rather than on the number of eggs laid by daughters. We also observed a negative association between fathers' oxidative damage levels and the number of eggs laid by daughters but a positive association between fathers' oxidative damage levels and the hatching success of those eggs. These opposing paternal effects canceled each other out, resulting in no overall effect on the number of offspring produced by daughters. No significant association between a female's own level of oxidative damage during reproduction and her reproductive success was observed. Our results suggest that oxidative damage experienced by parents is a better predictor of an individual's reproductive performance than oxidative damage experienced by the individual itself. Although the mechanisms underlying these parental condition transfer effects are currently unknown, changes in egg composition or (epi)genetic alterations of gametes may play a role. These findings highlight the importance of an intergenerational perspective when quantifying costs of physiological stress.

Elevational Changes in Bacterial Microbiota Structure and Diversity in an Arthropod-Disease Vector.

Environmental conditions change rapidly along elevational gradients and have been found to affect community composition in macroscopic taxa, with lower diversity typically observed at higher elevations. In contrast, microbial community responses to elevation are still poorly understood. Specifically, the effects of elevation on vector-associated microbiota have not been studied to date, even though the within-vector microbial community is known to influence vector competence for a range of zoonotic pathogens. Here we characterize the structure and diversity of the bacterial microbiota in an important zoonotic disease vector, the sheep tick Ixodes ricinus, along replicated elevational gradient (630-1673 m) in the Swiss Alps. 16S rRNA sequencing of the whole within-tick bacterial microbiota of questing nymphs and adults revealed a decrease in Faith's phylogenetic microbial alpha diversity with increasing elevation, while beta diversity analyses revealed a lower variation in microbial community composition at higher elevations. We also found a higher microbial diversity later in the season and significant differences in microbial diversity among tick life stages and sexes, with lowest microbial alpha diversity observed in adult females. No associations between tick genetic diversity and bacterial diversity were observed. Our study demonstrates systematic changes in tick bacterial microbiota diversity along elevational gradients. The observed patterns mirror diversity changes along elevational gradients typically observed in macroscopic taxa, and they highlight the key role of environmental factors in shaping within-host microbial communities in ectotherms.

Sex-specific effects of experimental ectoparasite infestation on telomere length in great tit nestlings.

Telomere length is a biomarker of biological ageing and lifespan in various vertebrate taxa. Evidence is accumulating that telomeres shorten more rapidly when an individual is exposed to environmental stressors. Parasites are potent selective agents that can cause physiological stress directly or indirectly through the activation of the host's immune system. Yet to date, empirical evidence for a role of parasites in telomere dynamics in natural populations is limited. Here, we show experimentally that exposure to ectoparasitic hen fleas (Ceratophyllus gallinae) during growth results in shorter telomeres in female, but not male, great tit (Parus major) nestlings. Females had longer telomeres than males when growing up in experimentally deparasitized nests but, likely because of the sex-specific effects of ectoparasitism on telomere length, this sexual dimorphism was absent in birds growing up in experimentally infested nests. Our results provide the first experimental evidence for a role of ectoparasitism in telomere dynamics in a natural vertebrate population, and suggest that the costs of infection manifest in sex-specific ways.

Selection for Divergent Reproductive Investment Affects Neuron Size and Foliation in the Cerebellum.

The cerebellum has a highly conserved internal circuitry, but varies greatly in size and morphology within and across species. Despite this variation, the underlying volumetric changes among the layers of the cerebellar cortex or their association with Purkinje cell numbers and sizes is poorly understood. Here, we examine intraspecific scaling relationships and variation in the quantitative neuroanatomy of the cerebellum in Japanese quail (Coturnix japonica) selected for high or low reproductive investment. As predicted by the circuitry of the cerebellum, the volumes of the constituent layers of the cerebellar cortex were strongly and positively correlated with one another and with total cerebellar volume. The number of Purkinje cells also significantly and positively co-varied with total cerebellar volume and the molecular layer, but not the granule cell layer or white matter volumes. Purkinje cell size and cerebellar foliation did not significantly covary with any cerebellar measures, but differed significantly between the selection lines. Males and females from the high-investment lines had smaller Purkinje cells than males and females from the low-investment lines and males from the high-investment lines had less folded cerebella than quail from the low-investment lines. These results suggest that within species, the layers of the cerebellum increase in a coordinated fashion, but Purkinje cell size and cerebellar foliation do not increase proportionally with overall cerebellum size. In contrast, selection for differential reproductive investment affects Purkinje cell size and cerebellar foliation, but not other quantitative measures of cerebellar anatomy.

Artificial selection reveals the role of transcriptional constraints in the maintenance of life history variation.

The trade-off between reproduction and self-maintenance is a cornerstone of life history theory, yet its proximate underpinnings are elusive. Here, we used an artificial selection approach to create replicated lines of Japanese quail (Coturnix japonica) that differ genetically in their reproductive investment. Whole transcriptome sequencing revealed that females from lines selected for high reproductive output show a consistent upregulation of genes associated with reproduction but a simultaneous downregulation of immune genes. Concordant phenotypic differences in immune function (i.e., specific antibody response against keyhole limpet hemocyanin) were observed between the selection lines, even in males who do not provide parental care. Our findings demonstrate the key role of obligate transcriptional constraints in the maintenance of life history variation. These constraints set fundamental limits to productivity and health in natural and domestic animal populations.

Combining genome-wide association study and FST -based approaches to identify targets of Borrelia-mediated selection in natural rodent hosts.

Recent advances in high-throughput sequencing technologies provide opportunities to gain novel insights into the genetic basis of phenotypic trait variation. Yet to date, progress in our understanding of genotype-phenotype associations in nonmodel organisms in general and natural vertebrate populations in particular has been hampered by small sample sizes typically available for wildlife populations and a resulting lack of statistical power, as well as a limited ability to control for false-positive signals. Here we propose to combine a genome-wide association study (GWAS) and FST -based approach with population-level replication to partly overcome these limitations. We present a case study in which we used this approach in combination with genotyping-by-sequencing (GBS) single nucleotide polymorphism (SNP) data to identify genomic regions associated with Borrelia afzelii resistance or susceptibility in the natural rodent host of this Lyme disease-causing spirochete, the bank vole (Myodes glareolus). Using this combined approach we identified four consensus SNPs located in exonic regions of the genes Slc26a4, Tns3, Wscd1 and Espnl, which were significantly associated with the voles' Borrelia infectious status within and across populations. Functional links between host responses to bacterial infections and most of these genes have previously been demonstrated in other rodent systems, making them promising new candidates for the study of evolutionary host responses to Borrelia emergence. Our approach is applicable to other systems and may facilitate the identification of genetic variants underlying disease resistance or susceptibility, as well as other ecologically relevant traits, in wildlife populations.

Bacterial microbiota composition of Ixodes ricinus ticks: the role of environmental variation, tick characteristics and microbial interactions.

Ecological factors, host characteristics and/or interactions among microbes may all shape the occurrence of microbes and the structure of microbial communities within organisms. In the past, disentangling these factors and determining their relative importance in shaping within-host microbiota communities has been hampered by analytical limitations to account for (dis)similar environmental preferences ('environmental filtering'). Here we used a joint species distribution modelling (JSDM) approach to characterize the bacterial microbiota of one of the most important disease vectors in Europe, the sheep tick Ixodes ricinus, along ecological gradients in the Swiss Alps. Although our study captured extensive environmental variation along elevational clines, the explanatory power of such large-scale ecological factors was comparably weak, suggesting that tick-specific traits and behaviours, microhabitat and -climate experienced by ticks, and interactions among microbes play an important role in shaping tick microbial communities. Indeed, when accounting for shared environmental preferences, evidence for significant patterns of positive or negative co-occurrence among microbes was found, which is indicative of competition or facilitation processes. Signals of facilitation were observed primarily among human pathogens, leading to co-infection within ticks, whereas signals of competition were observed between the tick endosymbiont Spiroplasma and human pathogens. These findings highlight the important role of small-scale ecological variation and microbe-microbe interactions in shaping tick microbial communities and the dynamics of tick-borne disease.