# Babesia bovis infection alters the composition and assembly of Rhipicephalus microplus midgut microbiota

Iván Corona-Guerrero, Apolline Maitre, Lianet Abuin-Denis, Rodrigo Morales-García, Consuelo Almazán, Dasiel Obregón, Alejandro Cabezas-Cruz, Juan Mosqueda

<figure><img src="image_1.png" /></figure>

## To cite this version:

Iván Corona-Guerrero, Apolline Maitre, Lianet Abuin-Denis, Rodrigo Morales-García, Consuelo Almazán, et al.. Babesia bovis infection alters the composition and assembly of Rhipicephalus microplus midgut microbiota. Frontiers in Microbiology, 2025, 16 (1), pp.96. 10.3389/fmicb.2025.1608409. hal-05315931

HAL Id: hal-05315931

https://hal.science/hal-05315931v1

Submitted on 15 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN ACCESS

EDITED BY
Shahid Karim,
The University of Southern Mississippi,
United States

REVIEWED BY
Vipin Rana,
University of Maryland, College Park,
United States
Surendra Raj Sharma,
The University of North Carolina at Chapel
Hill, United States

*CORRESPONDENCE
Juan Mosqueda

✉ joel.mosqueda@uaq.mx

Alejandro Cabezas-Cruz

✉ alejandro.cabezas@vet-alfort.fr

RECEIVED 09 April 2025

ACCEPTED 31 July 2025

PUBLISHED 02 September 2025

CITATION

Corona-Guerrero I, Maitre A, Abuin-Denis L, Morales-García R, Almazán C, Obregón D, Cabezas-Cruz A and Mosqueda J (2025) *Babesia bovis* infection alters the composition and assembly of *Rhipicephalus microplus* midgut microbiota. *Front. Microbiol.* 16:1608409. doi: 10.3389/fmicb.2025.1608409

COPYRIGHT

© 2025 Corona-Guerrero, Maitre, Abuin-Denis, Morales-García, Almazán, Obregón, Cabezas-Cruz and Mosqueda. This is an open-access article distributed under the terms of the [Creative Commons Attribution License (CC BY)](https://creativecommons.org/licenses/by/4.0/). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

# *Babesia bovis* infection alters the composition and assembly of *Rhipicephalus microplus* midgut microbiota

Iván Corona-Guerrero¹,², Apolline Maitre³,⁴,⁵,
Lianet Abuin-Denis³,⁶, Rodrigo Morales-García¹,
Consuelo Almazán¹, Dasiel Obregón⁷,
Alejandro Cabezas-Cruz³* and Juan Mosqueda¹,⁸*

¹Immunology and Vaccines Laboratory, C. A. Natural Sciences School, Autonomous University of Querétaro, Querétaro, Mexico. ²Doctorado en Ciencias Biológicas, Natural Sciences School, Autonomous University of Querétaro, Querétaro, Mexico. ³ANSES, INRAE, École Nationale Vétérinaire d'Alfort, UMR BIPAR, Laboratoire de Santé Animale, Maisons-Alfort, France. ⁴INRAE, UR 0045 Laboratoire de Recherches sur le Développement de l'Élevage (SELMET-LRDE), Corte, France. ⁵EA 7310, Laboratoire de Virologie, Université de Corse, Corte, France. ⁶Animal Biotechnology Department, Center for Genetic Engineering and Biotechnology, Havana, Cuba. ⁷School of Environmental Sciences, University of Guelph, Guelph, ON, Canada. ⁸C.A. Salud Animal y Microbiología Ambiental. Facultad de Ciencias Naturales, Universidad Autónoma de Querétaro, Querétaro, Mexico

**Introduction:** *Babesia bovis* is one of the main causative agents of bovine babesiosis. Livestock farmers are constantly struggling to control the population of the tick vector and reduce babesiosis outbreaks. For this reason, the development of new control strategies is necessary. Tick microbiota consists of a diverse group of symbiotic, commensal, and pathogenic microorganisms. It has been shown that altering the microbiota population prevents the transmission of apicomplexan pathogens. This work represents a primary exploratory approach to determine the changes *B. bovis* infection causes in the microbiota of *R. microplus*.

**Methods:** Two calves were infested with *R. microplus* larvae; next, one of the calves was splenectomized and infected with *Babesia bovis*. Fifteen days after the infestation, engorged females were collected from each calf. Collected ticks were separated into two groups: 0 h and 72 h. Ticks from the 0 h group were dissected to extract their midgut the same day they were collected, while midgut dissection of the other group was done after 72 h of incubation. Thus, samples were separated into 4 experimental groups depending on their infection status and the time of the dissection. Total DNA was purified and the V4 region of the bacterial 16S rRNA gene was sequenced using Illumina MiSeq technology.

**Results:** Data analysis showed fewer complex networks with reduced connectivity in infected ticks compared to the uninfected group. In both groups, the tick microbiota networks showed reduced node density at 72 h post-repletion. Different keystone taxa were found in all groups, indicating that midgut microbiota assembly is influenced by both tick developmental stage and the infection with *B. bovis*.
**Discussion:** Results of this work aim to serve as a steppingstone in the development of anti-tick microbiota vaccines capable of impairing both the life cycle of *R. microplus* and *B. bovis* transmission.

KEYWORDS

Rhipicephalus microplus, *Babesia bovis*, microbiota, network analysis, keystone bacteria

# 1 Introduction

Ticks are obligate hematophagous ectoparasites that inhabit temperate, tropical, and subtropical regions worldwide (Sonenshine and Roe, 2014). These parasites are extremely prolific, capable of remaining in the environment and infesting multiple hosts, making them efficient vectors for pathogens (Estrada-Peña, 2015). While typically displaying tropism toward their main host, ticks can occasionally infest incidental hosts in severe starvation cases (Anderson and Magnarelli, 2008).

The ixodid tick, *Rhipicephalus microplus* inhabits tropical and subtropical regions around the world and has a host tropism toward cattle (Esteve-Gasent et al., 2020). Particularly, Mexico is a country where most of its territory has a suitable climate for the development of *R. microplus* and, thus, the tick is present in a large fraction of the national territory (Diario Oficial de la Federación, 1995). According to the status of the control of *R. microplus*, Mexican authorities have divided the country into three zones: the free zone (there is no natural presence of the tick), the eradication zone (systematic treatments and constant epidemiologic vigilance are carried out to eradicate the tick) and the control zone (epidemiologic surveillance and treatments are performed to lower the impact of the tick) (Diario Oficial de la Federación, 1995). *Rhipicephalus microplus* is present in the totality of the control and eradication zones, which consist of approximately 69.40% of the country's territory (SENASICA, 2023).

The infestations caused by *R. microplus* can be extremely severe for livestock, causing economic losses due to weight loss, milk production decrease and damaged leather (Grisi et al., 2014). Ticks are also efficient vectors for pathogens due to their biology, ecology, and life cycle characteristics (Estrada-Peña, 2015). Tick prevalence prevents livestock farmers from introducing high-production breeds such as Holstein, Shorthorn or Hereford-Cebu due to their susceptibility to tick infestations in comparison to low production-high resistance breeds like Brahman and Brahman-crossbreeds (Jonsson et al., 2000; O'Kelly and Spiers, 1976). Effective tick control strategies are necessary to lower the population of the vector and reduce the incidence of tick-borne diseases such as bovine babesiosis.

*Babesia bovis* and *Babesia bigemina* are the most common causative agents of bovine babesiosis in Mexico (Alvarez et al., 2004). These apicomplexan parasites cause clinical manifestations such as acute fever, hemolytic anemia, hemoglobinuria, jaundice, anorexia, abortions, and weight loss and can result in the death of the animal (Bock et al., 2004). *Babesia bovis* causes aggregation of infected erythrocytes in lung and brain capillary blood vessels resulting in a higher mortality rate in susceptible adults (Brown and Palmer, 1999).

There are multiple strategies to control the tick population and *Babesia* infection. Acaricides have been universally used by livestock farmers for being cost-effective solutions to tick infestations, however, acaricides have important drawbacks such as safety concerns for non-human animals, soil and water accumulation, persistence in meat and selection of multi-resistant tick strains (Graf et al., 2004). These constraints made necessary the development of new integral strategies to control the tick population. An effective vaccine-centered strategy was developed by multiple groups, in which Bm86 of *R. microplus* (de la Fuente et al., 1999). Although the application of Bm86 vaccines was successful in reducing the use of acaricides to control tick infestations, its use has been limited due to the variability of Bm86 between tick strains (Merino et al., 2013).

Immune prophylaxis approaches have been employed to reduce the effect of *Babesia* on cattle. The first *Babesia* vaccines were made using live, attenuated parasites (Timms et al., 1983). Although live vaccines confer long-lasting protective immunity, they have short shelf life and the risk of reverting to a virulent phenotype (Timms et al., 1990). Recombinant vaccines, though safer, offer reduced long-term immune responses, carry a risk of reverting to *Babesia* pathogenicity, requiring cattle to be re-immunized to maintain protection against the pathogen (Jaramillo Ortiz et al., 2019; Jerzak et al., 2023).

A novel strategy employing vaccines that target the tick microbiota has been proposed. This approach is supported by evidence demonstrating that shifts in microbiota composition can impair vector feeding performance, offspring viability, and pathogen transmission capacity (Aželytė et al., 2022; Li et al., 2018; Mateos-Hernández et al., 2020; Narasimhan et al., 2017; Wu-Chuang et al., 2023). The objective of anti-tick microbiota vaccines is modulating the microbial community by immunizing the vertebrate host with a keystone taxon from microbiome to produce antibodies that, when ingested in the tick's bloodmeal, will target specific taxon, disturbing the tick microbiota assembly and functionality (Mateos-Hernández et al., 2021).

Tick microbiota is comprised of diverse bacteria, fungi, archaea, protozoa, and viruses (Xia et al., 2015; X.-L. Zhang et al., 2022). Obligate symbionts like *Wolbachia spp.* are considered core microbiota, while microorganisms such as *Babesia* or *Anaplasma*, obtained through interactions with the environment are considered part of the facultative microbiota (Bonnet and Pollet, 2021; Nikoh et al., 2014). The microorganisms of the tick microbiota perform multiple functions such as synthesis of vitamin D and other cofactors to compensate for the lack of these nutrients in the blood meal (Duron et al., 2018; Moran and Baumann, 2000).

To pinpoint the keystone taxa in ticks, it is essential to comprehensively characterize their tissue-specific microbiota,
including analyses of community composition, assembly
processes, and functional interactions. During each life stage,
multiple factors—including vertically transmitted symbionts,
environmental exposures during questing, host skin microbiota,
and acquired pathogens—dynamically reshape the microbial
community (Garcia Guizzo et al., 2023). Notably, certain pathogens
can induce dysbiosis within the vector, thereby enhancing their
infectivity and survival (Tolker-Nielsen, 2015).

To date, no studies have compared changes in microbiota
composition in response to *Babesia* infection in *R. microplus* ticks.
Understanding the microbiota of *R. microplus* is essential for the
development of anti-tick microbiota-based vaccines. The present
work represents a first exploratory study aiming to characterize
the microbiota composition of *R. microplus* at two stages of the
*Babesia bovis* infection cycle and to assess how the parasite alters the
tick's microbial network. Data obtained from this study will serve as
reference for the design of anti-tick microbiota vaccines.

**2 Materials and methods**

**2.1 Rhipicephalus microplus infestation and B. bovis infection**

Two European breed (Bos taurus) calves at least 6 months old and weighing less than 150 kg were selected for the experiment. The calves used for this experiment were obtained from the experimental ranch Valle del Guadiana in Durango, Mexico (part of the *R. microplus* control zone). Absence of anti-*Babesia* antibodies was confirmed by immunofluorescence antibody test (IFAT, not shown). Selected bovines were transferred to the Autonomous University of Queretaro (UAQ) Animal Infectomics Building and partially immobilized using a cattle chute to prevent the bovines from removing the ticks by grooming while allowing them to stretch, feed and lay down in non-stressing positions. Bovines had access to food and water *ad libitum*.

Because *Rhipicephalus microplus* is a one-host tick species, infestations in cattle must be initiated from the larval stage. Maintaining all the ticks in the same host removes the variable of facultative microbiota obtained from the interaction with different hosts and establishes each individual tick as a biological replica. Each calf was infested with 0.25 g of Media Joya strain *R. microplus* larvae (approx. 5,000 larvae). The larvae fed on the calves for 21-25 days, until they developed into fully engorged adult females. At 7 days post-infestation, one of the calves underwent a splenectomy via a left lateral approach at the Large Animal Veterinary Hospital of UAQ. Eight days after the splenectomy, it was infected with *Babesia bovis* via both intravenous and intramuscular routes, using two 2-ml vials containing 1 × 10⁸ infected erythrocytes (Hernández-Arvizu et al., 2024).

Clinical signs were monitored daily during the infection period (temperature, parasitemia, hematocrit, and other clinical signs related to *B. bovis* infection). Once the engorged ticks were collected from the infected calf, imidocarb dipropionate was administered at a dose of 3 mg/kg as therapeutic treatment for babesiosis. The Bioethics Committee of the Faculty of Natural Sciences at the Autonomous University of Querétaro (protocol number 053FCN2023) approved the non-human animal experimentation.

Starting at 21 days after the initial infestations, engorged female ticks were collected daily during a 5-day period. Collected ticks were transferred to the Immunology and Vaccines Research Laboratory (LINVAS) of the UAQ, Airport Campus to be processed. Ticks were washed by submersion in 10% benzalkonium chlorine for 5 min, rinsed with distilled water and dried with disposable paper towels. Half of the disinfected ticks were dissected the same day they were collected while the other half were incubated at 28°C and 80% relative humidity to be dissected after 72 h. At the end of the collection, ticks were separated into 4 groups: non-infected ticks dissected the same day of collection (0HN), non-infected ticks dissected 72 h after collection (72HN), infected ticks dissected the same day of collection (0HI), and infected ticks dissected 72 h after collection (72HI).

**2.2 Tick dissection and DNA extraction**

Ticks were fixed dorsally to petri dishes using double-sided tape. Ticks were opened by the posterior region of the body with a longitudinal cut to extract the midgut. Extracted tissues were washed using sterile phosphate-buffered saline (PBS, pH = 7.2) and stored individually in sterile PBS at -80°C. As a negative (blank) control for potential reagent or handling contamination, 500 µl of the same sterile PBS used for tissue rinsing was aliquoted, frozen, and carried through the entire DNA-extraction and sequencing workflow alongside the tick samples.

Tick organs were thawed in ice bath and centrifuged at
1000 × g for 30 s and the PBS was removed. Tick tissue was
frozen using liquid nitrogen, then 180 µl of ATL buffer and
20 µl of proteinase K solution from the DNeasy Blood and Tissue
Kit (QIAGEN, USA) were added. Immediately after, tissue was
macerated using a sterile micro pestle. Following that, extraction
was performed following instructions from the manufacturer. DNA
extraction was also performed in negative control PBS to check for
contamination. DNA integrity was verified by running an agarose
gel electrophoresis. A NanoDrop 2000 was used to determine DNA
concentration and purity.

**2.3 Amplicon sequencing of bacterial 16S rRNA gene**

DNA samples extracted from individual organs, with a 260/280 absorbance ratio ≥2.8 and concentration ≥10 ng/µl were sent to Novogene (Novogene Corporation Inc., Beijing, China). The V4 region of the 16S rRNA gene of bacteria and archaea was sequenced using the universal primers 505F/806R. All the samples were performed in the same batch. Analyses were performed on 7, 10, 8 and 9 samples from de 0HN, 72HN, 0HI and 72HI groups, respectively. These samples were sequenced using the pair-end approach in an Illumina MiSeq platform (Illumina, California, USA) service.

**2.4 Bioinformatic analysis**

The raw sequences were processed using the QIIME2 (version
2022.11) pipelines (Bolyen et al., 2019). Denoising, quality
trimming, and merging processes were made using DADA2 method (Callahan et al., 2016). Taxonomic classifications of ASVs were done using a classify-sklearn naïve Bayes taxonomic classifier, which was built upon the SILVA database (release 138) (Quast et al., 2012). The taxonomic data table was examined at species level to corroborate that *Babesia* mitochondrial 16s rRNA genes were identified only in infected samples. Next, the taxonomic data table was aggregated at the genus level and used for the assessment of microbiome composition and network analysis.

## 2.5 Statistical analysis

Microbial diversity metrics were calculated and ASVs level. Alpha-diversity analyses were based on microbial richness (observed ASVs) and evenness (Pielou evenness) for the two-time points under both conditions (uninfected vs. infected). Differences in alpha-diversity metrics between groups were assessed with the pairwise Kruskal-Wallis test (p < 0.05) within QIIME2 (Bolyen et al., 2019). Beta diversity measures were based on Bray-Curtis dissimilarity index (Bray and Curtis, 1957). Comparison between groups were performed using PERMANOVA test (p < 0.05), while beta dispersion (intragroup dispersion) was calculated and compared using ANOVA test (p < 0.05), using the betadisper function and the Vegan script (Oksanen et al., 2022) implemented in R (R Core Team, 2023). Bray-Curtis's dissimilarity was plotted using Principal Coordinate Analysis (PCoA; Gower, 2015). Then, we used Tukey's Honest Significant Difference (HSD) Test to make pair-wised comparisons of beta-dispersion between groups (Tukey, 1949). To assess differences in taxonomic composition between groups, the taxonomic composition table at the genera level and filtered to exclude the rare taxa (<10 sequences, and prevalence <20%). The ALDEx2 R package was used to transform the raw number of reads into center-log ratio values. This allows us to quantify and compare each taxa abundance within and between samples while avoiding an overall bacterial load bias (Fernandes et al., 2013). Only taxa with significant differences (p < 0.05) were selected for representation on a heatmaps of the differential abundance taxa, created using the "Heatplus" package in R (Ploner, 2024).

## 2.6 Co-occurrence networks analyses

Co-occurrence networks were built with the read number values from the taxonomic table using the SparCC algorithm from the SpiecEasi package (Kurtz et al., 2015) in R (version 4.3.1) (R Core Team, 2023), using a SparCC cutoff of 0.75. Networks were visualized and processed using Gephi (version 0.10) (Bastian et al., 2009). Next, we used the NetCoMi approach to determine significant differences in network metrics among groups using the Jaccard Index test, built into the NetCoMi R package to compare degree, betweenness centrality (BC), closeness centrality, eigenvector centrality and hub taxa of the most central nodes of each network (Peschel et al., 2020).

### 2.6.1 Network robustness

Robustness was evaluated by node removal to compare the network stability between different conditions. Network resistance

to random and direct node removal was tested using the swan_combinatory function of the NetSwan package (Lhomme, 2015). This involved calculating the fraction of nodes that needed to be removed to cause a connectivity loss of 80%, following both directed (e.g., based on betweenness, cascading effects and degree centrality) and random removal strategies. Directed node removal was performed by removing the highest connected node (degree method), removing the node with higher BC (betweenness method), and removing the node with the higher BC value, then recalculating BC before removing the next (cascading method). The standard error for loss of connectivity was calculated, considering variability, using a threshold of 0.975, using bootstrapping methods to calculate confidence intervals.

### 2.6.2 Determination of keystone taxa

To determine keystone taxa three criteria were employed: Firstly, these taxa had to be ubiquitous, signifying that were present across all samples within a particular group. Secondly, their eigenvector centrality had to be higher than 0.75 (Mateos-Hernández et al., 2021; Wu-Chuang et al., 2023). Eigenvector centrality is a metric for measuring the influence of a node within a network. Lastly, the mean relative abundance of these taxa needed to exceed that of the mean relative abundance of all taxa in the experimental group. To conduct this analysis, eigenvector centrality values were extracted from Gephi 0.10 software. All taxa were plotted on a two-axis graph using Graphpad Prism, indicating the cutoff points for CLR and Eigenvector centrality.

# 3 Results

## 3.1 Clinical signs of experimental calves

Starting the day of the infection, hematocrit, body temperature and parasitemia were checked daily. No abnormal clinical signs were observed on the uninfected calf. On the infected individual, no change on the clinical signs was observed from day 0 to day 5 post infection. On day 6, the first *Babesia bovis* trophozoite was detected in Giemsa-stained blood smears by light microscopy. Parasitemia was below 0.01% from day 0 to 7, increasing to 0.01% on day 8, to 0.2% on day 9 and to 4.60% on day 10. Body temperature remained at 38°C until day 7 post infection, increasing to 38.5°C on day 8, 40.6°C on day 9 and to 41.6°C on day 10. Hematocrit remained at 34% up until day 8, lowering to 27% on days 9 and 10. Imidocarb propionate was given to the infected calf on day 10 in response to acute clinical babesiosis.

## 3.2 Alpha and beta diversity metrics of R. microplus microbiota

Quality control results from Novogene showed that all samples passed the pre-amplification test, hence having the DNA concentration and purity necessary to carry out the Illumina sequencing process. Negative control samples showed no amplification, indicating no contamination was present.

To provide molecular confirmation of tick infection, we re-
quested the 16S-amplicon dataset for apicomplexan sequences
<figure><img src="image_4.png" /><figcaption>FIGURE 1</figcaption></figure>

Analysis of alpha and beta diversity of midgut microbiota. (A) Comparison of taxa richness showed a significant difference between the 0HI and 72HI groups (Kruskall-Wallis, *p* ≤ 0.05), the rest of the groups showed no significant difference (Kruskall-Wallis, *p* ≥ 0.05). (B) Pairwise comparison of taxa evenness showed no significant difference between samples (Kruskall-Wallis, *p* ≥ 0.05). (C) PCoA plot of multivariate dispersion analysis showed that all groups shared a similar microbial community composition (Tukey's HSD, *p* > 0.05). (D) Comparison of the relative abundance of taxa with significant differences between groups. Relative abundance data was transformed using CLR and expressed as a color gradient (Kruskall-Wallis, *p* < 0.05).

and detected *Babesia bovis* mitochondrial 16S rRNA fragments (≥93.1% identity to GenBank U06105.1) exclusively in the infected groups (0HI and 72HI; 116-146594 reads per sample) and absent from all uninfected samples. These parasite reads, removed prior to the bacterial network analysis as stipulated by our genus-level pipeline, are provided in Supplementary Table 1 together with their FASTA sequences.

Pairwise comparison of alpha diversity metrics showed significant difference in ASVs richness between 0HI and 72HI groups (Kruskall-Wallis, *p* ≤ 0.05) (Figure 1A). No significant difference was observed between the rest of the groups (Kruskall-Wallis, *p* > 0.05). Pairwise comparison of taxa evenness showed no significant difference between groups (Kruskall-Wallis, *p* > 0.05) (Figure 1B). The microbial community composition showed no significant difference between groups according to Permanova (*F* = 99, *p* > 0.5). The pairwise comparison also showed no statistical difference between groups (Tukey's HSD, *p* > 0.05) (Figure 1C).

## 3.3 Differential taxonomic composition

By examining the eukaryotic mitochondrial 16S rRNA genes identified during the taxonomic classification step, we detected

*Babesia* mitochondrial reads. These reads were present in samples from the 0HI and 72HI groups but were absent in the 0HN and 72HN groups. These results confirm the *B. bovis* infection in the infected groups.

A total of 228 unique bacterial genera were detected among all groups. Of these, 149 were found in the 0HN group, 216 in 72HN, 82 in 0HI, and 129 in 72HI (Supplementary Table 2). Among all samples, 46 taxa were shared, including: *Bacteroides*, *Borrelia*, *Coxiella* and *Escherichia-Shigella* (Supplementary Table 2). When comparing uninfected samples, we found 142 taxa shared between the 0 and 72 h groups, while 8 taxa were unique to the 0HN and 75 were unique to 72HN. In the *B. bovis*-infected samples, 58 taxa were shared, with 25 taxa unique to the 0HI group and 71 unique to the 72HI group. We found 6 taxa with a significantly different relative abundances between groups (*Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium*, *Comamonas*, *Acinetobacter*, *Pseudomonas*, *Rheinheimera* and *Flavobacterium*) (Kruskall-Wallis, *p* < 0.05) (Figure 1D).

## 3.4 Microbial co-occurrence network

The 0HN, 72HN, 0HI, and 72HI networks had 128/753, 58/90, 45/64, and 44/51 nodes/edges, respectively (Figure 2 and
<figure><img src="image_5.png" /><figcaption>FIGURE 2</figcaption></figure>

Co-occurrence networks of *R. microplus* midgut microbiota. Microbial co-occurrence networks of each group of samples (cutoff = 0.75). Nodes are colored according to the network module they belong to, and their size is proportional to its eigenvector centrality. Edges are colored according to their type of correlation (green for positive correlations and red to negative correlations). (A) Network of the 0HN group, composed of 128 nodes with 753 edges (100% positive correlations). (B) Network of the 72HN group, composed of 58 nodes with 90 edges (94.4% positive correlations). (C) Network of the 0HI group, composed of 45 nodes with 64 edges (100% positive correlations). (D) Network of the 72HI group, composed of 44 nodes with 51 edges (100% positive correlations).

**Table 1).** Both 0h networks and 72HI network had 100% positive correlations, while the 72HN network had 5 negative correlations (5.6%). In the four networks, a total of 220 unique taxa were found (Supplementary Table 3). From those, 57, 11, 13 and 60 were unique to the 0HN, 72HN, 0HI and 72HI groups, respectively. At the same time, only 2 taxa were shared among all 4 groups (Figure 3).

The 0HN exhibited the highest taxa richness and connectivity. Reduced richness and connectivity in 72HN, 0HI, and 72HI aligned with increased modularity. Modularity values were 0.402 (0HN), 0.659 (72HN), 0.809 (0HI), and 0.771 (72HI), indicating increasing modular partitioning. Additionally, networks had a clustering coefficient (CC) higher than 0.6 (0.719, 0.617, 0.759, 0.624 for 0HN, 72HN, 0HI and 72HI, respectively). Lower CCs in 72HN and 72HI suggest reduced interconnectivity over time post-repletion (Table 1).

## 3.5 Network topology comparison

The NetCoMi R package was used to compare the composition and topology of networks. Topology comparison between 0HN and 72HN showed a large taxa hub present in the first network, which remains after 72 h (Figure 4A). The remaining modules of both networks have different composition and topology. Comparing the 0HI with the 72HI we observed several taxa hubs, however, these hubs were composed of different taxa in each network (Figure 4B). Comparing the networks from uninfected and infected groups at 0 h post-repletion we can see that the large taxa hub in the 0HN network is not present in the 0HI group, additionally, the rest of the taxa clusters in the 0HN network are bigger than the ones in 0HI and have different components (Figure 4C). Each of the 72HN and 72HI networks have isolated taxa clusters with different compositions, indicating
<table><caption>TABLE 1 Topological features of the co-occurrence networks. A comparative analysis of the topological features is presented comparing metrics for each of the networks from the 0HN, 72HN, 0HI, and 72HI groups.</caption><thead><tr><th>Topological features</th><th>0HN</th><th>72HN</th><th>0HI</th><th>72HI</th></tr></thead><tbody><tr><td>Nodes</td><td>128</td><td>58</td><td>45</td><td>44</td></tr><tr><td>Edges</td><td>753</td><td>90</td><td>64</td><td>51</td></tr><tr><td>Positives</td><td>753 (100%)</td><td>85 (94.4%)</td><td>64 (100%)</td><td>51 (100%)</td></tr><tr><td>Negatives</td><td>0 (0%)</td><td>5 (5.6%)</td><td>0 (0%)</td><td>0 (0%)</td></tr><tr><td>Modularity</td><td>0.402</td><td>0.659</td><td>0.809</td><td>0.771</td></tr><tr><td>Clustering coefficient</td><td>0.719</td><td>0.617</td><td>0.759</td><td>0.624</td></tr><tr><td>Network diameter</td><td>10</td><td>4</td><td>3</td><td>4</td></tr></tbody></table>

a different microbial community composition (Figure 4D). Results of the Jaccard Index tests showed a significant difference in the composition of the most central nodes of both uninfected (Table 2) and infected (Table 3) networks at 0 versus 72 post-repletion hours and that this difference was a product of the composition *per se* ($P(< = Jacc)$) and not randomness ($P(> = Jacc)$).

## 3.6 Changes of network robustness

Robustness was evaluated by calculating the percentage of connectivity loss after removing nodes. Randomly removing nodes from the network showed no difference in the connectivity loss in any group (Figure 5A). By removing higher degree nodes first, we observed that 72HN group lost the most connectivity with the least fraction of nodes removed, followed by the 72HI group. 0HN and 0HI showed no noticeable difference when removing nodes by higher degree (Figure 5B). Removing nodes with the highest BC resulted in the 72HN losing connectivity with less fraction of node removal than the rest of groups. 0HN, 0HI and 72HI groups showed no difference among each other (Figure 5C). The cascading method showed that 72HN group was the least robust, followed by 0HN and 72HI and lastly 0HI (Figure 5D).

## 3.7 Identification of keystone taxa

The 0HN group had 6 keystone taxa: Xanthobacteraceae (uncultured) (family), *Lactobacillus*, *Moraxella*, Acidobacteria subgroup 2 (phylum), *Xanthomonas* and Elsterales (order) (Figure 6A and Table 4). The 6 keystone taxa found in the 72HN group were: *Cetobacterium*, *Dolosigranulum*, *Lactobacillus*, *Moraxella*, *Muribaculaceae* and *Turicibacter* (Figure 6B and Table 4). The keystone taxa found in the 0HI group were the following: *Acidibacter*, *Anaerococcus*, *Candidatus Nomurabacteria*, *Escherichia-Shigella*, *Haliangium*, and *Ralstonia* (Figure 6C and Table 4). Lastly, keystone taxa found in the 72HI were: *Planococcaceae* (family), *Bacteroides*, *Koukoulia*, Burkholderiales SC-I-84 (order) (Figure 6D and Table 4).

# 4 Discussion

It has previously been reported that colonization by pathogenic microorganisms can alter the tick microbiota composition during the infection process (Maitre et al., 2022; Narasimhan et al., 2017). Our findings support the hypothesis that microbiota from *R. microplus* differs between uninfected ticks and ticks infected with *B. bovis*. Additionally, the microbiota composition in both uninfected and uninfected ticks changed significantly between 0 and 72 h after engorgement. Understanding the shifts in microbiota composition between these two time points is crucial for elucidating its role in *Babesia* infections. *B. bovis* invasion of the tick midgut begins when merozoites are ingested by ticks at 0 h post-repletion. Inside the tick, these merozoites develop into sexual stages (gametocytes), which subsequently transform into infective forms that penetrate the midgut epithelium. By 72 h post-repletion, these forms differentiate into kinetes, a stage coinciding with cellular effects such as the upregulation and re localization of the voltage-dependent anion channel (BmVDAC) in *R. microplus* midgut cells (Rodríguez-Hernández et al., 2012, 2015).

Tick susceptibility to pathogen infection is mediated by multiple mechanisms. First, humoral components of the tick immune system like anti-microbial peptides are capable of binding to the parasite surface and inhibit their proliferation (Tsuji et al., 2007; Yao et al., 2021). Additionally, it has been demonstrated that microbiota composition can alter multiple arthropod susceptibility to pathogen infection in ixodid models transmitting pathogens like *Rickettsia helvetica* or (Gall et al., 2016; Maitre et al., 2022), *Anaplasma phagocytophilum* (Abraham et al., 2017) and *Babesia microti* (Wei et al., 2021).

## 4.1 Decrease in microbiota diversity

During the feeding process of ixodid ticks, the midgut lumen fills with the blood meal and functions as a reservoir, separated from the midgut epithelium by the peritrophic matrix (PM; Sonenshine and Roe, 2014). Nutrients and cellular debris can selectively cross the PM and bind to membrane receptors of epithelial cells to initiate receptor-mediated endocytosis. Once internalized, the nutrients are digested intracellularly and transported to the basal lamina. This process is taxing for the epithelial tissue, which must constantly replace damaged cells, ultimately leading to tissue degradation at the onset of oviposition (Sonenshine and Roe, 2014).

Our analytical pipeline—centered log-ratio normalization with ALDEx2 followed by differential-abundance testing— was intentionally chosen to fulfill the study's primary aim: to capture the emergent, community-level effects of *Babesia bovis* infection on the *Rhipicephalus microplus* mid-gut microbiota, rather than to track quantitative shifts of individual bacterial taxa. These methods provide robust relative-abundance data that underpin the network-level reconfiguration reported here. At the same time, we recognize the complementary value of absolute-quantification approaches such as taxon-specific qPCR or droplet-digital PCR. Future investigations focusing on the
<table><caption>TABLE 2 Comparing 0HN and 72HN networks topological features using the Jaccard index test. Comparison between network metrics of networks from uninfected midgut samples at 0- and 72-h post-repletion.</caption><thead><tr><th>Uninfected 0 h vs. 72 h</th><th>Jaccard index</th><th>P (≤Jacc)</th><th>P (≥Jacc)</th></tr></thead><tbody><tr><td>Degree</td><td>0.247</td><td>0.065426</td><td>0.961308</td></tr><tr><td>Betweenness centrality</td><td>0.214</td><td>0.036709*</td><td>0.982061</td></tr><tr><td>Closeness centrality</td><td>0.188</td><td>0.001148**</td><td>0.999501</td></tr><tr><td>Eigenvector centrality</td><td>0.163</td><td>0.000127***</td><td>0.999953</td></tr><tr><td>Hub taxa</td><td>0.163</td><td>0.000127***</td><td>0.999953</td></tr></tbody></table>

<table><caption>TABLE 3 Comparing 0HI and 72HI networks topological features using the Jaccard index test. Comparison between network metrics of networks from infected midgut samples at 0- and 72-h post-repletion.</caption><thead><tr><th>Infected 0 h vs. 72 h</th><th>Jaccard index</th><th>P (≤Jacc)</th><th>P (≥Jacc)</th></tr></thead><tbody><tr><td>Degree</td><td>0.247</td><td>0.065426</td><td>0.961308</td></tr><tr><td>Betweenness centrality</td><td>0.214</td><td>0.036709*</td><td>0.982061</td></tr><tr><td>Closeness centrality</td><td>0.188</td><td>0.001148**</td><td>0.999501</td></tr><tr><td>Eigenvector centrality</td><td>0.163</td><td>0.000127***</td><td>0.999953</td></tr><tr><td>Hub taxa</td><td>0.163</td><td>0.000127***</td><td>0.999953</td></tr></tbody></table>

mechanistic roles of the differentially abundant or keystone genera identified in this study could employ these techniques to validate and extend the hypotheses generated by our community-scale analysis.

Alpha diversity metrics showed no significant differences in richness or evenness of taxa between samples, except for a difference in richness between the 72HI and 72HN groups. A PCoA analysis indicated a similar taxonomic composition among groups. These results suggest a comparable number of individual reads across samples, with a uniform proportion of taxa. However, co-occurrence network analysis revealed a significant decrease in both the number and connectivity of nodes in the

0HI group compared to the 0HN group. When comparing the 72-h and 0-h groups, we observed a decrease in the number of nodes in the microbial networks of uninfected individuals; however, this decrease in diversity was not observed in the infected groups. These findings indicate that *B. bovis* infection has a direct impact on the microbiota diversity at 0 h-post repletion. Similar results have been previously reported in *R. microplus* where the infection with *Theileria sp.* correlates with a loss of bacterial diversity and richness on the vector's tissues (Adegoke et al., 2020). Another tick model showing a similar response is *Ixodes scapularis*, where infection with *Anaplasma phagocytophilum* alters the microbiota composition (Abraham et al., 2017). Microbiota alterations caused by pathogens have also been observed in other vectors such as *Culex quinquefasciatus* where *Plasmodium relictum* infection decreases the diversity of the midgut microbial community (Aželytė et al., 2022). Different symbionts and pathogenic bacteria colonize specific tissues depending on their roles in the microbial community and the tick's life stage (Wiesinger et al., 2023; X.-Y. Zhang et al., 2022). These community dynamics might explain the decrease in diversity observed at 72 h post-repletion rather than being solely attributable to *B. bovis* infection.

## 4.2 Comparison of co-occurrence networks

The 0HN network had the highest node count of all the experimental groups. When comparing this group to 72HN, we can see that the number of nodes and edges decrease. This loss of connectivity can be attributed to the midgut degradation (Sonenshine and Roe, 2014). Additionally, a similar decrease in these topology features can be observed when comparing the 0HN to the 0HI group. Microbial networks from these two groups have a similar number of nodes and edges, and the shared taxa are highly dissimilar (sharing only 12.66% of their total taxa). This phenomenon can be attributed to damage to the PM by the *Babesia* arrowhead structure-contained enzymes (Jalovecka et al., 2018) or as collateral damage caused by non-specific

<figure><img src="image_6.png" /><figcaption>FIGURE 3</figcaption></figure>

Venn diagram indicating the percentage of shared taxa between the 0HN (yellow), 72HN (green), 0HI (red) and 72HI (purple) networks. (A) Out of 146 total taxa, 0HN and 72HN share 40, 88 are unique to 0HN and 18 to 72HN. (B) Out of 146 total taxa, 0HN and 0HI share 27, 101 are unique to 0HN and 18 to 0HI. (C) The 72HN and 72HI groups have 94 taxa between both, sharing 8, 50 being unique to 72HN and 36 to 72HI. (D) The 0HI and 72HI groups have a total of 79 taxa, sharing 10 of them, while 35 are unique to 0HI and 34 to 72HI. (E) UpSet plot showing taxa shared among the 4 groups.

<footnote>*p < 0.05, **p < 0.005, ***p < 0.005.</footnote>

<footnote>*p < 0.05, **p < 0.005, ***p < 0.005.</footnote>
<figure><img src="image_7.png" /><figcaption>FIGURE 4</figcaption></figure>

Network topology comparison using NetCoMi. Network topology comparison using the NetCoMi R package. In each comparison each dot represents a taxon, and taxa are located at the same position. (A) Comparison between OHN and 72HN networks. (B) Comparison between 0HI and 72HI networks. (C) Comparison between OHN and 0HI networks. (D) Comparison between 72HN and 72HI networks.
<figure><img src="image_8.png" /><figcaption>FIGURE 5</figcaption></figure>

Robustness test by node removal. A network vulnerability test by node removal was conducted for the 0HN (orange), 72HN (green), 0HI (red) and 72HI (purple) networks using four different methods. Nodes were removed using at random (A), by higher degree (B), by higher betweenness centrality (C) and using a cascading method (D). En each graph, the X-axis represents the fraction of nodes that were removed from the network and the Y-axis represents the fraction of connectivity that was lost in response to that removal. The solid-colored line represents the relationship between the fraction of removed nodes of each sample with the connectivity loss and the shaded area is the standard error. The horizontal black line shows 0.8 of the total connectivity and vertical-colored lines shows the fraction of nodes removed at which the network loses 0.8 of the connectivity, color code is presented in the bottom right.

innate immune response from the tick (Sonenshine and Hynes, 2008).

Seeing the significant differences in network composition in all groups, it became of interest to determine the robustness of the networks. The node removal robustness test showed that 0HN network retained more connectivity when removing a higher number of nodes compared with other networks. This could be explained by its more diverse community with complex interactions. It has been mentioned that microbial networks with these characteristics tend to be more resistant to connectivity loss in response to perturbations on the network composition (Konopka et al., 2015). This resistance to change referrers to compositional stability, not functional. We still need to determine if the difference in composition also implies a change in the metabolic function of the network.

Although we complemented SparCC-derived co-occurrence networks with robustness testing (NetSwan) and node-composition statistics (NetCoMi), we recognize that network topology can vary with the choice of inference algorithm, correlation threshold and compositional correction. A systematic sensitivity analysis—e.g., recalculating networks across a range of SparCC cut-offs and comparing results with alternative approaches such as SPIEC-EASI, CoNet, CCLasso or FlashWeave—was beyond the scope of the present work but is essential to gauge the stability of keystone-taxon identification and global metrics (modularity, clustering coefficient, eigenvector centrality). Cross-method validation and threshold testing are therefore critical next steps for consolidating community-level patterns and verifying that taxa highlighted as potential targets for anti-microbiota interventions are consistently recovered across analytical frameworks.
<figure><img src="image_9.png" /><figcaption>FIGURE 6</figcaption></figure>

Keystone taxa determination based on relative abundance and connectivity. Graphs representing the keystone taxa of the 0HN (A), 72HN (B), 0HI (C) and 72HI (D). Each graph represents the relative abundance as CLR-transformed values in the X-axis. The connectivity of each taxon is represented in the Y-axis with the eigenvector centrality value. Each taxon is represented as a dot. The vertical dotted line represents the mean of the CLR-transformed abundance data, and the horizontal dotted line represents the empiric 0.75 cutoff of the eigenvector centrality value (A–C) or the mean of the eigenvector centrality value (D). A dot was considered as keystone taxa if it was above both cutoffs. The previously reported *Escherichia-Shigella* cluster (Mateos-Hernández et al., 2021) was colored in orange in the 0HI graph (C). Dots corresponding to the known tick endosymbionts, *Borrelia* and *Coxiella*, were colored yellow in all graphs.

## 4.3 Comparison of keystone taxa

A single host can have different stable microbiota compositions. The assembly of each composition is influenced by environmental and dietary factors (Shetty et al., 2017). For each stable composition, keystone taxa can exert important influence in the microbiome function (Banerjee et al., 2018). Although previously reported symbionts of *R. microplus* such as *Coxiella* and *Borrelia* (Garcia-Guizzo et al., 2022; LaRocca et al., 2009) are present in all samples' microbial communities, their lack of connectivity in the network assembly rules them out of the keystone taxa status. Uninfected groups share *Lactobacillus* and *Moraxella* as keystone taxa at both 0 and 72-h post repletion. The exact role that *Lactobacillus spp.* plays in the microbiome of ticks has not been determined, however, it has been deemed as a probiotic for the mulberry silkworm, *Bombyx mori* (Yeruva et al., 2020). Additionally, it has been shown to be highly abundant on the midgut microbiota of the soft tick *Ornithodoros moubata*, which made it an effective anti-microbiota vaccine target (Cano-Argüelles et al., 2024). *Moraxella spp.* has been associated with lipid metabolism in insects (Li et al., 2021). Lipid metabolism is necessary for vitellogenesis in insects and ticks

during ovary maturation (Amaral Xavier et al., 2018; Tufail and Takeda, 2008), which makes this keystone taxa a potential vaccine candidate.

Infected groups don't share any keystone taxa between the 0 and 72-post repletion hours. However, the *Escherichia-Shigella* cluster (ESC) was present in the 0HI group as a keystone taxon. Although works on *B. mori* suggest the ESC could be involved in the degradation of carbohydrates, the exact role of this taxa hasn't been determined (Zheng et al., 2020). However, it has been previously used as an anti-microbiota vaccine antigen capable of altering the microbial community, potentially reducing metabolic pathways such as lysine degradation (Mateos-Hernández et al., 2021). The use of ESC as an anti-tick microbiota antigen has efficiently reduced the load of the pathogenic bacterium *Borrelia afzelii* in the vector (Wu-Chuang et al., 2023).

Other keystone taxa found in different groups might have different functions. For example, *Xanthomonas* was determined as a keystone taxon in the 0HN group. Bovine used for this experiment were obtained from field condition and plant-associated microbiota was most likely obtained from its interaction with the bovine host skin. This is what would be expected, since components
<table><caption>TABLE 4 Keystone taxa found for each co-occurrence network.</caption><thead><tr><th>Keystone taxa (0HN)</th><th>CLR</th><th>EVC</th></tr></thead><tbody><tr><td>Xanthobacteraceae (uncultured) (f)</td><td>0.842</td><td>0.828</td></tr><tr><td>Lactobacillus</td><td>2.016</td><td>0.761</td></tr><tr><td>Moraxella</td><td>0.943</td><td>0.953</td></tr><tr><td>Acidobacteria (p) subgroup 2</td><td>0.929</td><td>1.000</td></tr><tr><td>Xanthomonas</td><td>0.758</td><td>0.869</td></tr><tr><td>Elsterales (o)</td><td>0.864</td><td>0.880</td></tr><tr><th>Keystone Taxa (72HN)</th><th>CLR</th><th>EVC</th></tr><tr><td>Cetobacterium</td><td>0.965</td><td>0.790</td></tr><tr><td>Dolosigranulum</td><td>1.476</td><td>0.976</td></tr><tr><td>Lactobacillus</td><td>1.882</td><td>1.000</td></tr><tr><td>Moraxella</td><td>1.646</td><td>0.971</td></tr><tr><td>Muribaculaceae</td><td>1.515</td><td>0.961</td></tr><tr><td>Turicibacter</td><td>0.910</td><td>0.974</td></tr><tr><th>Keystone taxa (0HI)</th><th>CLR</th><th>EVC</th></tr><tr><td>Acidibacter</td><td>1.149</td><td>1.000</td></tr><tr><td>Anaerococcus</td><td>0.761</td><td>0.869</td></tr><tr><td>Candidatus Nomurabacteria</td><td>0.928</td><td>0.869</td></tr><tr><td>Escherichia-Shigella*</td><td>0.936</td><td>0.903</td></tr><tr><td>Haliangium</td><td>0.845</td><td>0.802</td></tr><tr><td>Ralstonia</td><td>0.834</td><td>0.899</td></tr><tr><th>Keystone taxa (7HI)</th><th>CLR</th><th>EVC</th></tr><tr><td>Planococcaceae (f)</td><td>0.753</td><td>0.080</td></tr><tr><td>Bacteroides</td><td>0.792</td><td>0.080</td></tr><tr><td>Koukoulia</td><td>0.994</td><td>0.104</td></tr><tr><td>Burkholderiales (o) SC-I-84</td><td>1.330</td><td>0.361</td></tr></tbody></table>

from the host microbiota are often found in the vector midgut (Choubdar et al., 2021). We found *Muribaculaceae* among the keystone taxa of the 72HN group. This taxon was previously found as keystone taxa of the salivary glands' microbiota of soft ticks, although its role on the microbiota it's not yet known (Piloto-Sardiñas et al., 2023). *Turicibacter* is another keystone taxa found in the 72HI samples, its role in the tick's midgut microbial community hasn't been explored. However, in the stag beetle *Dorcus hopei* model, *Turicibacter* is highly abundant in the gut microbiota of the larval stages and plays a role in energy metabolism and its positively related to higher growth (Lu et al., 2024).

We found *Anaerococcus* as a keystone taxon in the 0HI group. This bacterium is associated with short chain fatty acid production and promotion of anti-inflammatory state in mucous tissue of *Aedes albopictus* (Wei et al., 2023), although its role within tick microbiota is not known. *Acidibacter*, another keystone taxa found in the 0HI group, has been previously reported as a high abundance bacterium in *Haemaphysalis longicornis* (Zhang et al., 2023). Other taxa in the infected groups like *Haliangium*, *Ralstonia*, *Bacteroides* and *Koukoulia* have previously been reported in the microbiota of ticks and other arthropods (Keskin et al., 2017; Wei et al., 2013),

they haven't been reported as being highly prevalent or having a key role in their host's metabolism (Keskin et al., 2017; Wei et al., 2013).

It has been previously reported that different bacteria species of the tick microbiota are correlated to their host's susceptibility or resistance to tick borne pathogen infections (Gall et al., 2016). This effect can be achieved by manipulating their host gene expression. For example: *Anaplasma phagocytophilum* can induce the expression of the *Ixodes scapularis* anti-freeze glycoprotein (*iafgp*) gene, which thins the PM enabling the invasion of the pathogen (Abraham et al., 2017). It is important to take in count the roles of the keystone taxa when selecting the bacteria that will be used as an antigen in anti-tick microbiota vaccines. We must make sure the antigen candidates are positively correlated to tick fecundity and pathogen susceptibility. Additionally, we must determine that no other taxa in the network have redundant functions to effectively disrupt the microbiota composition (Konopka et al., 2015).

## 4.4 Implications in vaccine design

As mentioned previously, the removal of most central nodes from a network result in a quicker connection loss. Due to keystone taxa being the more central and abundant members of the microbial community, targeting those with an anti-microbiota vaccine is a specific and efficient way to remove these bacteria from the midgut microbiota (Mateos-Hernández et al., 2021). Previous studies showed that removal of one keystone species using vaccination strategies can induce major changes in the overall assembly of the microbial community (Aželytė et al., 2022). For these reasons it is important to study the composition of the tick microbiota and determine their keystone taxa to adequately design vaccines. Immunization against keystone taxa is necessary to cause a major disruption in the assembly of the tick midgut microbial community.

The present work also took into consideration the microbial community composition shift caused by the infection with *B. bovis*. The observed change in overall network composition and differential keystone taxa enables us to choose a keystone bacterium that, being used as a vaccine target, could both impair the *B. bovis* infection on the tick to stop its propagation and cause a dysbiosis state in the tick that hiders its biological fitness.

## 4.5 Most viable vaccine candidates

From the keystone taxa we identified, there are several we'd like to highlight as the most viable vaccine candidates. *Lactobacillus sp.* is a promising candidate which has been previously tested as a vaccine antigen due to being associated with bovine lactic acidosis (Shu et al., 1999). Immunized bovines produced good antibody titters and developed no adverse reactions in their rumen. However, our approach didn't specify the *Lactobacillus* species present in the midgut microbiota, so it would be necessary to isolate and identify it. Isolation of *Lactobacillus* spp. requires the use of special culture media and anaerobic conditions, which further complicates the isolation and growth of this genera.

<footnote>*Previously reported as Keystone taxa by Mateos-Hernández et al. (2021)</footnote>
*Moraxella* spp. is another viable vaccine candidate. Vaccination against *Moraxella bovis* has been proven effective to prevent bovine keratoconjunctivitis-causing infections (Angelos *et al.*, 2021). Same as with *Lactobacillus* spp., it is necessary to determine the species of *Moraxella* present in the tick midgut. *Moraxella* spp. can grow under standard conditions (Baumann *et al.*, 1968), which could facilitate the migration to large scale vaccine production if necessary.

We can also consider the ESC as a potential vaccine candidate. It's possible to use only one of the cluster's bacteria to make a vaccine. *E. coli* can grow under permissive conditions and there are multiple commercially available strains. The use of *E. coli* in experimental vaccines have been proven effective to disrupt microbiota assembly and pathogen colonization in ticks (Mateos-Hernández *et al.*, 2020; Wu-Chuang *et al.*, 2023). In this case, the challenge of identifying the exact strain present in the tick midgut remains.

## 5 Limitations

In this exploratory work, we were able to determine composition and keystone taxa of the *R. microplus* midgut microbiota. Nonetheless, our approach presents limitations that open opportunity areas for further projects. One of such limitations is that this study focused on two milestones of the blood-digestion cycle—immediately after engorgement (0 h) and the late digestive phase (72 h)—because they bracket the window in which *B. bovis* parasites are first ingested and the time when they leave the midgut cells to infect other tissues. Nevertheless, microbiota composition can fluctuate markedly during the intervening hours as the PM forms, hemoglobin catabolism progresses, and epithelial turnover accelerates. By only sampling the start and end of this interval we may have missed short-lived taxa, transient interactions, or tipping points that precede the network re-configuration reported here. Future work that couples denser time-series sampling (e.g., 6, 12, 24, 48, 72 h) with absolute-quantification methods will be essential to reconstruct these intermediate transitions, clarify the chronology of keystone-taxon emergence, and determine whether the community-level patterns identified at 72 h persist, re-equilibrate, or diverge during the oviposition period.

By knowing the composition of the midgut microbiota, it's possible to make a prediction of the metabolic profile of the bacterial assembly. By comparing the taxa present in the microbial networks and their relative abundance from data obtained through Illumina sequencing and comparing it to databases we can determine the metabolic capabilities of the microbiota. Metabolic prediction was outside the scope of this project, additionally, only the 16s rRNA marker was used, which could limit the metabolic prediction. A more comprehensive metabolic profile prediction would require the use of multiple molecular markers to include fungi, viruses and protozoa.

## 6 Conclusion

From an epidemiological standpoint, *B. bovis* is a difficult pathogen to control. It is prevalent in areas with high cattle density,

young calves are resistant to infection and serve as reservoirs, and vectors such as *R. microplus* can persist in the field for prolonged periods (Jerzak *et al.*, 2023; Spickler, 2018). Understanding the effects that *B. bovis* induces in the tick microbiota is a necessary step toward elucidating the vector-microbiota-pathogen interactions that underpin *R. microplus* vectorial capacity. In the present study, we investigated differences in midgut microbiota composition in engorged *R. microplus* females at various post-repletion time points and under different *B. bovis* infection statuses. The results indicate that *B. bovis* infection alters the structure and reduces the bacterial diversity of the *R. microplus* midgut microbiota at 0 h post-repletion. Notably, taxa associated with energy metabolism, such as *Moraxella* and *Turicibacter*, were absent in infected samples. Among the keystone taxa identified in different groups, we found bacteria previously reported as keystone species and as targets for anti-microbiota vaccines in other tick species, including *Lactobacillus*, *Moraxella*, *Muribaculaceae*, and *Escherichia-Shigella*. This information can be leveraged to advance the development of anti-tick microbiota vaccines. We hypothesize that using keystone taxa from both infected and uninfected midgut microbiota as antigens represents an optimal vaccine-based tick control strategy, potentially reducing *R. microplus* fitness and preventing the transmission of *B. bovis*.

## Data availability statement

The status of the bioproject page was changed to released and a new public URL was provided: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1216386.

## Ethics statement

The animal study was approved by the Bioethics Committee, Faculty of Natural Sciences, Autonomous University of Querétaro. The study was conducted in accordance with the local legislation and institutional requirements.

## Author contributions

IC-G: Formal Analysis, Investigation, Methodology, Writing – original draft. AM: Investigation, Methodology, Writing – review & editing. LA-D: Writing – review & editing. RM-G: Investigation, Writing – review & editing. CA: Investigation, Writing – review & editing. DO: Methodology, Software, Writing – review & editing. AC-C: Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing. JM: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – review & editing.

## Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The present work was
founded by CONAHCyT project number 321162 in collaboration with ANUIES- ECOS NORD.

## Acknowledgments

We thank the personnel of the Animal Infectomics Building for their work in ensuring the wellbeing of the experimental animals during the duration of the present work, specially to Emmanuel Castañeda Villafranca for his invaluable help with overseeing the experimental bovines.

## Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling editor SK declared a past co-authorship/collaboration [Fourteen anti-tick vaccine targets are variably conserved in cattle fever ticks -10.1186/S13071-025-06683-5 (2025)] with the author JM.

## Generative AI statement

The authors declare that no Generative AI was used in the creation of this manuscript.

## Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1608409/full#supplementary-material

## References

Abraham, N. M., Liu, L., Jutras, B. L., Yadav, A. K., Narasimhan, S., Gopalakrishnan, V., et al. (2017). Pathogen-mediated manipulation of arthropod microbiota to promote infection. *Proc. Natl. Acad. Sci. U. S. A.* 114, E781-E790. doi: 10.1073/pnas.1613422114

Adegoke, A., Kumar, D., Bobo, C., Rashid, M. I., Durrani, A. Z., Sajid, M. S., et al. (2020). Tick-Borne pathogens shape the native microbiome within tick vectors. *Microorganisms* 8:1299. doi: 10.3390/microorganisms8091299

Alvarez, J. A., Ramos, J. A., Rojas, E. E., Mosqueda, J. J., Vega, C. A., Olvera, A. M., et al. (2004). Field challenge of cattle vaccinated with a combined *Babesia bovis* and *Babesia bigemina* frozen immunogen. *Ann. N. Y. Acad. Sci.* 1026, 277–283. doi: 10.1196/annals.1307.043

Amaral Xavier, M., Tirloni, L., Pinto, A. F. M., Diedrich, J. K., Yates, J. R., Mulenga, A., et al. (2018). A proteomic insight into vitellogenesis during tick ovary maturation. *Sci. Rep.* 8:4698. doi: 10.1038/s41598-018-23090-2

Anderson, J. F., and Magnarelli, L. A. (2008). Biology of ticks. *Infect. Dis. Clin. N. Am.* 22, 195–215. doi: 10.1016/j.idc.2007.12.006

Angelos, J. A., Elizalde, P., and Griebel, P. (2021). Bovine immune responses to *Moraxella bovis* and *Moraxella bovoculi* following vaccination and natural or experimental infections. *Vet. Clin. North Am. Food Anim. Pract.* 37, 253–266. doi: 10.1016/j.cvfa.2021.03.002

Aželytė, J., Wu-Chuang, A., Žiegytė, R., Platonova, E., Mateos-Hernandez, L., Maye, J., et al. (2022). Anti-Microbiota vaccine reduces avian malaria infection within mosquito vectors. *Front. Immunol.* 13:841835. doi: 10.3389/fimmu.2022.841835

Banerjee, S., Schlaeppi, K., and Van Der Heijden, M. G. A. (2018). Keystone taxa as drivers of microbiome structure and functioning. *Nat. Rev. Microbiol.* 16, 567–576. doi: 10.1038/s41579-018-0024-1

Bastian, M., Heymann, S., and Jacomy, M. (2009). Gephi: An open source software for exploring and manipulating networks. *Proc. Int. AAAI Conf. Web Soc. Media* 3, 361–362. doi: 10.1609/icwsm.v3i1.13937

Baumann, P., Doudoroff, M., and Stanier, R. Y. (1968). Study of the moraxella group I. genus moraxella and the *Neisseria catarrhalis* group. *J. Bacteriol.* 95, 58–73. doi: 10.1128/jb.95.1.58-73.1968

Bock, R., Jackson, L., De Vos, A., and Jorgensen, W. (2004). Babesiosis of cattle. *Parasitology* 129, S247–S269. doi: 10.1017/S0031182004005190

Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. *Nat. Biotechnol.* 37, 852–857. doi: 10.1038/s41587-019-0209-9

Bonnet, S. I., and Pollet, T. (2021). Update on the intricate tango between tick microbiomes and tick-borne pathogens. *Parasite Immunol.* 43:e12813. doi: 10.1111/pim.12813

Bray, J. R., and Curtis, J. T. (1957). An ordination of the upland forest communities of Southern Wisconsin. *Ecol. Monographs* 27, 325–349. doi: 10.2307/1942268

Brown, W. C., and Palmer, G. H. (1999). Designing blood-stage vaccines against *Babesia bovis* and *B. bigemina*. *Parasitol. Today* 15, 275–281. doi: 10.1016/S0169-4758(99)01471-4

Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data. *Nat. Methods* 13, 581–583. doi: 10.1038/nmeth.3869

Cano-Argüelles, A. L., Piloto-Sardiñas, E., Maitre, A., Mateos-Hernández, L., Maye, J., Wu-Chuang, A., et al. (2024). Microbiota-driven vaccination in soft ticks: Implications for survival, fitness and reproductive capabilities in *Ornithodoros moubata*. *Mol. Ecol.* 33:e17506. doi: 10.1111/mec.17506

Choubdar, N., Karimian, F., Koosha, M., and Oshaghi, M. A. (2021). An integrated overview of the bacterial flora composition of *Hyalomma anatolicum*, the main vector of CCHF. *PLoS Negl. Trop. Dis.* 15:e0009480. doi: 10.1371/journal.pntd.009480

de la Fuente, J., Rodríguez, M., Montero, C., Redondo, M., García-García, J. C., Méndez, L., et al. (1999). Vaccination against ticks (*Boophilus spp.): The experience with the Bm86-based vaccine Gavac™. *Genet. Anal. Biomol. Eng.* 15, 143–148. doi: 10.1016/S1050-3862(99)00018-2

Diario Oficial de la Federación (1995). *Norma oficial Mexicana nom-019-200-1994, campaña nacional contra la garrapata Boophilus spp.* Mexico: Diario Oficial de la Federación.

Duron, O., Morel, O., Noël, V., Buysse, M., Binetruy, F., Lancelot, R., et al. (2018). Tick-Bacteria mutualism depends on B vitamin synthesis pathways. *Curr. Biol.* 28, 1896-1902.e5. doi: 10.1016/j.cub.2018.04.038

Esteve-Gasent, M. D., Rodríguez-Vivas, R. I., Medina, R. F., Ellis, D., Schwartz, A., Cortés García, B., et al. (2020). Research on integrated management for cattle fever ticks and bovine babesiosis in the United States and Mexico: Current status and opportunities for binational coordination. *Pathogens* 9:871. doi: 10.3390/pathogens9110871

Estrada-Peña, A. (2015). Ticks as vectors: Taxonomy, biology and ecology. *Rev. Sci. Tech. l'OIE* 34, 53–65. doi: 10.20506/rst.34.1.2345

Fernandes, A. D., Macklaim, J. M., Linn, T. G., Reid, G., and Gloor, G. B. (2013). ANOVA-Like differential expression (ALDEX) analysis for mixed population RNA-Seq. *PLoS One* 8:e67019. doi: 10.1371/journal.pone.0067019

Gall, C. A., Reif, K. E., Scoles, G. A., Mason, K. L., Mousel, M., Noh, S. M., et al. (2016). The bacterial microbiome of *Dermacentor andersoni* ticks influences pathogen susceptibility. *ISME J.* 10, 1846–1855. doi: 10.1038/ismej.2015.266
Garcia Guizzo, M., Hatalová, T., Frantová, H., Zurek, L., Kopáček, P., and Perner, J. (2023). Ixodes ricinus ticks have a functional association with *Midichloria mitochondrii*. *Front. Cell. Infect. Microbiol.* 12:1081666. doi: 10.3389/fcimb.2022.1081666

Garcia-Guizzo, M., Tirloni, L., Gonzalez, S. A., Farber, M. D., Braz, G., Parizi, L. F., et al. (2022). Coxiella endosymbiont of *Rhipicephalus microplus* modulates tick physiology with a major impact in blood feeding capacity. *Front. Microbiol.* 13:868575. doi: 10.3389/fmicb.2022.868575

Gower, J. C. (2015). *Principal coordinates analysis*. Wiley StatsRef: Statistics Reference Online. doi: 10.1002/9781118445112.stat05670.pub2

Graf, J.-F., Gogolewski, R., Leach-Bing, N., Sabatini, G. A., Molento, M. B., Bordin, E. L., et al. (2004). Tick control: An industry point of view. *Parasitology* 129, S427-S442. doi: 10.1017/S0031182004006079

Grisi, L., Leite, R. C., Martins, J. R. D. S., Barros, A. T. M. D., Andreotti, R., Cançado, P. H. D., et al. (2014). Reassessment of the potential economic impact of cattle parasites in Brazil. *Rev. Bras. Parasitol*. Vet. 23, 150–156. doi: 10.1590/S1984-29612014042

Hernández-Arvizu, E. E., Asada, M., Kawazu, S.-I., Vega, C. A., Rodríguez-Torres, A., Morales-García, R., et al. (2024). Antiparasitic evaluation of aquiluscidin, a cathelicidin obtained from *Crotalus aquilus*, and the Vcn-23 derivative peptide against *Babesia bovis*, *B. bigemina* and *B. ovata*. *Pathogens* 13:496. doi: 10.3390/pathogens13060496

Jalovecka, M., Hajdusek, O., Sojka, D., Kopacek, P., and Malandrin, L. (2018). The complexity of piroplasms life cycles. *Front. Cell. Infect. Microbiol.* 8:248. doi: 10.3389/fcimb.2018.00248

Jaramillo Ortiz, J. M., Paoletta, M. S., Gravisaco, M. J., López Arias, L. S., Montenegro, V. N., de la Fournière, S. A. M., et al. (2019). Immunisation of cattle against *Babesia bovis* combining a multi-epitope modified vaccinia Ankara virus and a recombinant protein induce strong Th1 cell responses but fails to trigger neutralising antibodies required for protection. *Ticks Tick-Borne Dis.* 10:101270. doi: 10.1016/j.ttbdis.2019.101270

Jerzak, M., Gandurski, A., Tokaj, M., Stachera, W., Szuba, M., and Dybicz, M. (2023). Advances in *Babesia* vaccine development: An overview. *Pathogens* 12:300. doi: 10.3390/pathogens12020300

Jonsson, N. N., Matschoss, A. L., Pepper, P., Green, P. E., and Ansell, J. (2000). Resistance of holstein-friesian cows to infestation by the cattle tick (*Boophilus microplus*). *Vet. Parasitol.* 89, 297–305. doi: 10.1016/S0304-4017(00)00213-2

Keskin, A., Bursali, A., Snow, D. E., Dowd, S. E., and Tekin, S. (2017). Assessment of bacterial diversity in *Hyalomma aegyptium*, *H. marginatum* and *H. excavatum* ticks through tag-encoded pyrosequencing. *Exp. Appl. Acarol.* 73, 461–475. doi: 10.1007/s10493-017-0186-y

Konopka, A., Lindemann, S., and Fredrickson, J. (2015). Dynamics in microbial communities: Unraveling mechanisms to identify principles. *ISME J.* 9, 1488–1495. doi: 10.1038/ismej.2014.251

Kurtz, Z. D., Müller, C. L., Miraldi, E. R., Littman, D. R., Blaser, M. J., and Bonneau, R. A. (2015). Sparse and compositionally robust inference of microbial ecological networks. *PLoS Comp. Biol.* 11:e1004226. doi: 10.1371/journal.pcbi.1004226

LaRocca, T. J., Holthausen, D. J., Hsieh, C., Renken, C., Mannaella, C. A., and Benach, J. L. (2009). The bactericidal effect of a complement-independent antibody is osmolytic and specific to *Borrelia*. *Proc. Natl. Acad. Sci. U. S. A.* 106, 10752–10757. doi: 10.1073/pnas.0901858106

Lhomme, S. (2015). NetSwan: Network strengths and weaknesses analysis (Version R package version 0.1) [Computer software]. Available online at: https://CRAN.R-project.org/package=NetSwan. (accessed September 25, 2024).

Li, H., Zhao, C., Yang, Y., Zhou, Z., Qi, J., and Li, C. (2021). The influence of gut microbiota on the fecundity of *Henosepilachna vigintioctopunctata* (Coleoptera: Coccinellidae). *J. Insect Sci.* 21:15. doi: 10.1093/jisesa/ieab061

Li, L.-H., Zhang, Y., and Zhu, D. (2018). Effects of antibiotic treatment on the fecundity of *Rhipicephalus haemaphysaloides* ticks. *Parasites Vectors* 11:242. doi: 10.1186/s13071-018-2807-7

Lu, Y., Chu, S., Shi, Z., You, R., and Chen, H. (2024). Marked variations in diversity and functions of gut microbiota between wild and domestic stag beetle *Dorcus hopei* *hopei*. *BMC Microbiol.* 24:24. doi: 10.1186/s12866-023-03177-1

Maitre, A., Wu-Chuang, A., Mateos-Hernández, L., Foucault-Simonin, A., Moutailler, S., Paoli, J.-C., et al. (2022). *Rickettsia helvetica* infection is associated with microbiome modulation in *Ixodes ricinus* collected from humans in Serbia. *Sci. Rep.* 12:11464. doi: 10.1038/s41598-022-15681-x

Mateos-Hernández, L., Obregón, D., Maye, J., Borneres, J., Versille, N., de la Fuente, J., et al. (2020). Anti-Tick microbiota vaccine impacts *Ixodes ricinus* performance during feeding. *Vaccines* 8:702. doi: 10.3390/vaccines8040702

Mateos-Hernández, L., Obregón, D., Wu-Chuang, A., Maye, J., Bornères, J., Versillé, N., et al. (2021). Anti-Microbiota vaccines modulate the tick microbiome in a taxon-specific manner. *Front. Immunol.* 12:704621. doi: 10.3389/fimmu.2021.704621

Merino, O., Alberdi, P., Pérez De La Lastra, J. M., and De La Fuente, J. (2013). Tick vaccines and the control of tick-borne pathogens. *Front. Cell. Infect. Microbiol.* 3:30. doi: 10.3389/fcimb.2013.00030

Moran, N. A., and Baumann, P. (2000). Bacterial endosymbionts in animals. *Curr. Opin. Microbiol.* 3, 270–275. doi: 10.1016/S1369-5274(00)00088-6

Narasimhan, S., Schuij, T. J., Abraham, N. M., Rajeevan, N., Coumou, J., Graham, M., et al. (2017). Modulation of the tick gut milieu by a secreted tick protein favors *Borrelia burgdorferi* colonization. *Nat. Commun.* 8:184. doi: 10.1038/s41467-017-00208-0

Nikoh, N., Hosokawa, T., Moriyama, M., Oshima, K., Hattori, M., and Fukatsu, T. (2014). Evolutionary origin of insect–Wolbachia nutritional mutualism. *Proc. Natl. Acad. Sci. U. S. A.* 111, 10257–10262. doi: 10.1073/pnas.1409284110

O’Kelly, J. C., and Spiers, W. G. (1976). Resistance to *Boophilus microplus* (Canestrini) in genetically different types of calves in early life. *J. Parasitol.* 62:312. doi: 10.2307/3279295

Oksanen, J., Simpson, G. L., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., et al. (2022). *vegan: Community ecology package*. Available online at: http://CRAN.R-project.org/package=vegan (accessed October 29, 2024).

Peschel, S., Müller, C. L., Boulesteix, A.-L., and Depner, M. (2020). NetCoMi: Network construction and comparison for microbiome data in R. *Brief. Bioinform.* 22:bbaa290. doi: 10.1093/bib/bbaa290

Piloto-Sardiñas, E., Cano-Argüelles, A. L., Maitre, A., Wu-Chuang, A., Mateos-Hernández, L., Corduneanu, A., et al. (2023). Comparison of salivary gland and midgut microbiome in the soft ticks *Ornithodoros erraticus* and *Ornithodoros moubata*. *Front. Microbiol.* 14:1173609. doi: 10.3389/fmicb.2023.1173609

Ploner, A. (2024). Heatplus: Heatmaps with row and/or column covariates and colored clusters (Version 3.1.20). Bioconductor Project (Fred Hutchinson Cancer Research Center, Seattle): Seattle, WA, USA.

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. *Nucleic Acids Res.* 41, D590–D596. doi: 10.1093/nar/gks1219

R Core Team (2023). *R: A language and environment for statistical computing* (Version 4.3.1). Vienna: R Foundation for Statistical Computing.

Rodríguez-Hernández, E., Mosqueda, J., Alvarez-Sánchez, M. E., Neri, A. F., Mendoza-Hernández, G., and Camacho-Nuez, M. (2012). The identification of a VDAC-like protein involved in the interaction of *Babesia bigemina* sexual stages with *Rhipicephalus microplus* midgut cells. *Vet. Parasitol.* 187, 538–541. doi: 10.1016/j.vetpar.2012.01.028

Rodríguez-Hernández, E., Mosqueda, J., León-Ávila, G., Castañeda-Ortiz, E. J., Álvarez-Sánchez, M. E., Camacho, A. D., et al. (2015). BmVDAC upregulation in the midgut of *Rhipicephalus microplus*, during infection with *Babesia bigemina*. *Vet. Parasitol.* 212, 368–374. doi: 10.1016/j.vetpar.2015.06.016

SENASICA (2023). *Situación actual Campaña Nacional para el control de la garrapata Boophilus spp.* México: Servicio Nacional de Sanidad, Inocuidad y Calidad Agroalimentaria.

Shetty, S. A., Hugenholtz, F., Lahti, L., Smidt, H., and De Vos, W. M. (2017). Intestinal microbiome landscaping: Insight in community assemblage and implications for microbial modulation strategies. *FEMS Microbiol. Rev.* 41, 182–199. doi: 10.1093/femsre/fuw045

Shu, Q., Gill, H. S., Hennessy, D. W., Leng, R. A., Bird, S. H., and Rowe, J. B. (1999). Immunisation against lactic acidosis in cattle. *Res. Vet. Sci.* 67, 65–71. doi: 10.1053/rvsc.1998.0284

Sonenshine, D. E., and Hynes, W. L. (2008). Molecular characterization and related aspects of the innate immune response in ticks. *Front. Biosci.* 13, 7046–7063. doi: 10.2741/3209

Sonenshine, D. E., and Roe, R. M. (2014). *Biology of ticks*, 2nd Edn. Oxford: Oxford University Press.

Spickler, A. R. (2018). *Bovine babesiosis*. Available online at: www.cfsp.iastate.edu/Factsheets/pdfs/bovine_babesiosis.pdf (accessed September 25, 2024).

Timms, P., Dalgliesh, R. J., Barry, D. N., Dimmock, C. K., and Rodwell, B. J. (1983). *Babesia bovis*: Comparison of culture-derived parasites, non-living antigen and conventional vaccine in the protection of cattle against heterologous challenge. *Aus. Vet. J.* 60, 75–77. doi: 10.1111/j.1751-0813.1983.tb05874.x

Timms, P., Stewart, N. P., and De Vos, A. J. (1990). Study of virulence and vector transmission of *Babesia bovis* by use of cloned parasite lines. *Infect. Immun.* 58, 2171–2176. doi: 10.1128/iai.58.7.2171-2176.1990

Tolker-Nielsen, T. (2015). Biofilm development. *Microbiol. Spect.* 3:MB-0001-2014. doi: 10.1128/microbiolspec.MB-0001-2014

Tsuji, N., Battsetseg, B., Boldbaatar, D., Miyoshi, T., Xuan, X., Oliver, J. H., et al. (2007). Babesial vector tick defensin against *Babesia sp.* parasites. *Infect. Immun.* 75, 3633–3640. doi: 10.1128/IAI.00256-07

Tufail, M., and Takeda, M. (2008). Molecular characteristics of insect vitellogenins. *J. Insect Physiol.* 54, 1447–1458. doi: 10.1016/j.jinsphys.2008.08.007

Tukey, J. W. (1949). Comparing individual means in the analysis of variance. *Biometrics* 5:99. doi: 10.2307/3001913
Wei, N., Cao, J., Zhang, H., Zhou, Y., and Zhou, I. (2021). The tick microbiota dysbiosis promote tick-borne pathogen transstadial transmission in a *Babesia microti*-infected mouse model. *Front. Cell. Infect. Microbiol.* 11:713466. doi: 10.3389/fcimb. 2021.713466

Wei, T., Hu, J., Miyanaga, K., and Tanji, Y. (2013). Comparative analysis of bacterial community and antibiotic-resistant strains in different developmental stages of the housefly (Musca domestica). *Appl. Microbiol. Biotechnol.* 97, 1775-1783. doi: 10.1007/ s00253-012-4024-1

Wei, X., Peng, H., Li, Y., Meng, B., Wang, S., Bi, S., et al. (2023). Pyrethroids exposure alters the community and function of the internal microbiota in *Aedes albopictus*. *Ecotoxicol. Environ. Safety* 252:114579. doi: 10.1016/j.ecoenv.2023. 114579

Wiesinger, A., Wenderlein, J., Ulrich, S., Hiereth, S., Chitimia-Dobler, L., and Straubinger, R. K. (2023). Revealing the tick microbiome: Insights into midgut and salivary gland microbiota of female *Ixodes ricinus* ticks. *Int. J. Mol. Sci.* 24:1100. doi: 10.3390/ijms24021100

Wu-Chuang, A., Mateos-Hernandez, L., Maitre, A., Rego, R. O. M., Šima, R., Porcelli, S., et al. (2023). Microbiota perturbation by anti-microbiota vaccine reduces the colonization of *Borrelia afzelii* in *Ixodes ricinus*. *Microbiome* 11:151. doi: 10.1186/ s40168-023-01599-7

Xia, H., Hu, C., Zhang, D., Tang, S., Zhang, Z., Kou, Z., et al. (2015). Metagenomic profile of the viral communities in *Rhipicephalus spp.* ticks from Yunnan, China. *PLoS One* 10:e0121609. doi: 10.1371/journal.pone.0121609

Yao, J., Xu, Z., Sun, Z., Zhou, K., Lu, J., Mi, R., et al. (2021). Response of leucine-rich repeat domain-containing protein in *Haemaphysalis longicornis* to *Babesia microti* infection and its ligand identification. *Infect. Immun.* 89:e00268-20. doi: 10.1128/IAI. 00268-20

Yeruva, T., Vankadara, S., Ramasamy, S., and Lingaiah, K. (2020). Identification of potential probiotics in the midgut of mulberry silkworm, *Bombyx mori* through metagenomic approach. *Probiot. Antimicrob. Proteins* 12, 635-640. doi: 10.1007/ s12602-019-09580-3

Zhang, X.-L., Deng, Y.-P., Yang, T., Li, L.-Y., Cheng, T.-Y., Liu, G.-H., et al. (2022). Metagenomics of the midgut microbiome of *Rhipicephalus microplus* from China. *Parasites Vectors* 15:48. doi: 10.1186/s13071-022-05161-6

Zhang, X.-Y., Li, S.-S., Chen, K.-L., Yang, C., Zhou, X.-J., Liu, J.-Z., et al. (2022). Growth dynamics and tissue localization of a Coxiella-like endosymbiont in the tick *Haemaphysalis longicornis*. *Ticks Tick-Borne Dis.* 13:102005. doi: 10.1016/j.ttbdis.2022. 102005

Zhang, Y.-K., Li, S.-S., Yang, C., Zhang, Y.-F., and Liu, J.-Z. (2023). Mechanism of the toxic effects of tetracycline on blood meal digestion in *Haemaphysalis longicornis*. *Exp. Appl. Acarol.* 91, 681-695. doi: 10.1007/s10493-023-00 858-4

Zheng, Y., Xiao, G., Zhou, W., Gao, Y., Li, Z., Du, G., et al. (2020). Midgut microbiota diversity of potato tuber moth associated with potato tissue consumed. *BMC Microbiol.* 20:58. doi: 10.1186/s12866-020-01 740-8