The order Rickettsiales (Alphaproteobacteria) encompasses multiple diverse lineages of host-associated bacteria, including pathogens, reproductive manipulators, and mutualists. Here, in order to understand how intracellularity and host association originated in this order, and whether they are ancestral or convergently evolved characteristics, we built a large and phylogenetically-balanced dataset that includes de novo sequenced genomes and a selection of published genomic and metagenomic assemblies. We perform detailed functional reconstructions that clearly indicates “late” and parallel evolution of obligate host-association in different Rickettsiales lineages. According to the depicted scenario, multiple independent horizontal acquisitions of transporters led to the progressive loss of biosynthesis of nucleotides, amino acids and other metabolites, producing distinct conditions of host-dependence. Each clade experienced a different pattern of evolution of the ancestral arsenal of interaction apparatuses, including development of specialised effectors involved in the lineage-specific mechanisms of host cell adhesion and/or invasion.

Ticks are important vectors of many pathogens in Europe, where the most impactful species is Ixodes ricinus. Recently, the geographical distribution of this tick species has been expanding, resulting in an increased risk of human exposure to tick bites. With the present study, we aimed to screen 350 I. ricinus specimens collected from humans and wild animals (mainly ungulates), to have a broader understanding of the tick-borne pathogens circulating in the Lombardy region, in northern Italy. To do so, we took advantage of a high-throughput real-time microfluidic PCR approach to screen ticks in a cost-effective and time-saving manner. Molecular analysis of the dataset revealed the presence of four genera of bacteria and two genera of protozoa: in ungulates, 77 % of collected ticks carried Anaplasma phagocytophilum, while the most common pathogen species in ticks removed from humans were those belonging to Borrelia burgdorferi sensu lato group (7.6 %). We also detected other pathogenic microorganisms, such as Rickettisa monacensis, Rickettsia helvetica, Neoehrlichia mikurensis, Babesia venatorum, and Hepatozoon martis. Besides, we also reported the presence of the pathogenic agent Borrelia miyamotoi in the area (1.4 % overall). The most common dual co-infection detected in the same tick individual involved A. phagocytophilum and Rickettsia spp. Our study provided evidence of the circulation of different tick-borne pathogens in a densely populated region in Italy.

Summary: Members of the bacterial order Rickettsiales are obligatorily associated with a wide range of eukaryotic hosts. Their evolutionary trajectories, in particular concerning the origin of shared or differential traits among distant sub‐lineages, are still poorly understood. Here, we characterized a novel Rickettsiales bacterium associated with the ciliate Paramecium tredecaurelia and phylogenetically related to the Rickettsia genus. Its genome encodes significant lineage‐specific features, chiefly the mevalonate pathway gene repertoire, involved in isoprenoid precursor biosynthesis. Not only this pathway has never been described in Rickettsiales , it also is very rare among bacteria, though typical in eukaryotes, thus likely representing a horizontally acquired trait. The presence of these genes could enable an efficient exploitation of host‐derived intermediates for isoprenoid synthesis. Moreover, we hypothesize the reversed reactions could have replaced canonical pathways for producing acetyl‐CoA, essential for phospholipid biosynthesis. Additionally, we detected phylogenetically unrelated mevalonate pathway genes in metagenome‐derived Rickettsiales sequences, likely indicating evolutionary convergent effects of independent horizontal gene transfer events. Accordingly, convergence, involving both gene acquisitions and losses, is highlighted as a relevant evolutionary phenomenon in Rickettsiales , possibly favoured by plasticity and comparable lifestyles, representing a potentially hidden origin of other more nuanced similarities among sub‐lineages.

Ticks require bacterial symbionts for the provision of necessary compounds that are absent in their hematophagous diet. Such symbionts are frequently vertically transmitted and, most commonly, belong to the Coxiella genus, which also includes the human pathogen Coxiella burnetii. This genus can be divided in four main clades, presenting partial but incomplete cocladogenesis with the tick hosts. Here, we report the genome sequence of a novel Coxiella, endosymbiont of the African tick Amblyomma nuttalli, and the ensuing comparative analyses. Its size (∼1 Mb) is intermediate between symbionts of Rhipicephalus species and other Amblyomma species. Phylogenetic analyses show that the novel sequence is the first genome of the B clade, the only one for which no genomes were previously available. Accordingly, it allows to draw an enhanced scenario of the evolution of the genus, one of parallel genome reduction of different endosymbiont lineages, which are now at different stages of reduction from a more versatile ancestor. Gene content comparison allows to infer that the ancestor could be reminiscent of C. burnetii. Interestingly, the convergent loss of mismatch repair could have been a major driver of such reductive evolution. Predicted metabolic profiles are rather homogenous among Coxiella endosymbionts, in particular vitamin biosynthesis, consistently with a host-supportive role. Concurrently, similarities among Coxiella endosymbionts according to host genus and despite phylogenetic unrelatedness hint at possible host-dependent effects.

The zoonotic parasiteCryptosporidium parvumis a global cause of gastrointestinal disease in humans and ruminants. Sequence analysis of the highly polymorphicgp60gene enabled the classification ofC. parvumisolates into multiple groups (e.g., IIa, IIc, Id) and a large number of subtypes. In Europe, subtype IIaA15G2R1 is largely predominant and has been associated with many water- and food-borne outbreaks. In this study, we generated new whole-genome sequence (WGS) data from 123 human- and ruminant-derived isolates collected in 13 European countries and included other available WGS data from Europe, Egypt, China, and the United States (n = 72) in the largest comparative genomics study to date. We applied rigorous filters to exclude mixed infections and analyzed a data set from 141 isolates from the zoonotic groups IIa (n = 119) and IId (n = 22). Based on 28,047 high-quality, biallelic genomic SNPs, we identified three distinct and strongly supported populations: Isolates from China (IId) and Egypt (IIa and IId) formed population 1; a minority of European isolates (IIa and IId) formed population 2; and the majority of European (IIa, including all IIaA15G2R1 isolates) and all isolates from the United States (IIa) clustered in population 3. Based on analyses of the population structure, population genetics, and recombination, we show that population 3 has recently emerged and expanded throughout Europe to then, possibly from the United Kingdom, reach the United States, where it also expanded. The reason(s) for the successful spread of population 3 remain elusive, although genes under selective pressure uniquely in this population were identified.