ABSTRACT: Tick-borne diseases (TBDs) continue to emerge and re-emerge in several regions of the world, highlighting the need for novel and effective control strategies. The development of effective strategies requires a better understanding of TBDs ecology, and given the complexity of these systems, interdisciplinary approaches are required. In recent years, the microbiome of vectors has received much attention, mainly because associations between native microbes and pathogens may provide a new promising path towards the disruption of pathogen transmission. However, we still do not fully understand how host genetics and environmental factors interact to shape the microbiome of organisms, or how pathogenic microorganisms affect the microbiome and vice versa. The integration of different lines of evidence may be the key to improve our understanding of TBDs ecology. In that context, we generated microbiome and pathogen presence data for Dermacentor variabilis, and integrated those data sets with population genetic data, and metadata for the same individual tick specimens. Clustering and multivariate statistical methods were used to combine, analyze, and visualize data sets. Interpretation of the results is challenging, likely due to the low levels of genetic diversity and the high abundance of a few taxa in the microbiome. Francisella was dominant in almost all ticks, regardless of geography or sex. Nevertheless, our results showed that, overall, ticks from different geographic regions differ in their microbiome composition. Additionally, DNA of Rickettsia rhipicephali, R. montanensis, R. bellii, and Anaplasma spp., was detected in D. variabilis specimens. This is the first study that successfully generated microbiome, population genetics, and pathogen presence data from the same individual ticks, and that attempted to combine the different lines of evidence. The approaches and pre-processing steps used can be applied to a variety of taxa, and help better understand ecological processes in biological systems.
