Stéphane Robin, Jean-François Cosson, Julien Chiquet, Cédric Midoux, Olivier Rué, Mahendra Mariadassou, Xavier Bailly, J. Aubert, Emilie Lejal, A. Cougoul, P. Gasqui, Thomas Pollet, K. Chalvet-Monfray, Muriel Vayssier-Taussat, Agustín Estrada-Peña, Biologie moléculaire et immunologie parasitaires et fongiques (BIPAR), École nationale vétérinaire d'Alfort (ENVA)-Laboratoire de santé animale, sites de Maisons-Alfort et de Dozulé, Agence nationale de sécurité sanitaire de l'alimentation, de l'environnement et du travail (ANSES)-Agence nationale de sécurité sanitaire de l'alimentation, de l'environnement et du travail (ANSES)-Université Paris-Est Créteil Val-de-Marne - Paris 12 (UPEC UP12)-Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE), Mathématiques et Informatique Appliquées (MIA-Paris), Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE)-AgroParisTech-Université Paris-Saclay, Faculty of Veterinary Medicine, University of Zaragoza - Universidad de Zaragoza [Zaragoza], Université Paris-Saclay, INRAE, BioinfOmics, MIGALE bioinformatics facility (MIGALE), Université Paris-Saclay-Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE), Mathématiques et Informatique Appliquées du Génome à l'Environnement [Jouy-En-Josas] (MaIAGE), Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE), Université Paris-Saclay, Procédés biotechnologiques au service de l'environnement (UR PROSE), Unité Mixte de Recherche d'Épidémiologie des maladies Animales et zoonotiques (UMR EPIA), VetAgro Sup - Institut national d'enseignement supérieur et de recherche en alimentation, santé animale, sciences agronomiques et de l'environnement (VAS)-Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE), Département Santé Animale (DEPT SA), Animal, Santé, Territoires, Risques et Ecosystèmes (UMR ASTRE), Centre de Coopération Internationale en Recherche Agronomique pour le Développement (Cirad)-Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE), INRAE, Métaprogramme Méta-omiques et Ecosystèmes Microbiens (MEM), Adaptation of Agriculture and Forests to Climate Change (AAFCC), Région Ile-de-France, ANR-17-CE32-0011,NGB,Biosurveillance Next-Gen des changements dans la structure et le fonctionnement des écosystèmes(2017), AgroParisTech-Université Paris-Saclay-Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE), Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE)-Centre de Coopération Internationale en Recherche Agronomique pour le Développement (Cirad), École nationale vétérinaire - Alfort (ENVA)-Laboratoire de santé animale, sites de Maisons-Alfort et de Dozulé, and Mathématiques et Informatique Appliquées (MIA Paris-Saclay)
BackgroundTicks transmit pathogens of medical and veterinary importance, and represent an increasing threat for human and animal health. Important steps in assessing disease risk and developing possible new future control strategies involve identifying tick-borne microbes, their temporal dynamics and interactions.MethodsUsing high throughput sequencing, we studied the microbiota dynamics ofIxodes ricinusfrom 371 nymphs collected monthly over three consecutive years in a peri-urban forest. After adjusting a Poisson Log Normal model to our data set, the implementation of a principal component analysis as well as sparse network reconstruction and differential analysis allowed us to assess inter-annual, seasonal and monthly variability ofI. ricinusmicrobial communities as well as their interactions.ResultsAround 75% of the detected sequences belonged to five genera known to be maternally inherited bacteria in arthropods and potentially circulating in ticks: CandidatusMidichloria, Rickettsia, Spiroplasma, ArsenophonusandWolbachia. The structure of theI. ricinusmicrobiota was temporally variable with interannual recurrence and seemed to be mainly driven by OTUs belonging to environmental genera. The total network analysis revealed a majority of positive (partial) correlations. We identified strong relationships between OTUs belonging toWolbachiaandArsenophonus, betraying the presence of the parasitoid waspIxodiphagus hookeriin ticks, and the well known arthropod symbiontSpiroplasma, previously documented to be involved in the defense against parasitoid wasp inDrosophila melanogaster. Other associations were observed between the tick symbiont CandidatusMidichloriaand pathogens belonging toRickettsia, probablyRickettsia helvetica. More specific network analysis finally suggested that the presence of pathogens belonging to generaBorrelia, AnaplasmaandRickettsiamight disrupt microbial interactions inI. ricinus.ConclusionsHere, we identified theI. ricinusmicrobiota and documented for the first time the existence and recurrence of marked temporal shifts in the tick microbial community dynamics. We statistically showed strong relationships between the presence of some pathogens and the structure of theI. ricinusnon-pathogenic microbes. We interestingly detected close links between some tick symbionts and the potential presence of either pathogenicRickettsiaor a parasitoid in ticks. All these new findings might be very promising for the future development of new control strategies of ticks and tick-borne diseases.