Optimising machine learning prediction of minimum inhibitory concentrations in Klebsiella pneumoniae.

Minimum Inhibitory Concentrations (MICs) are the gold standard for quantitatively measuring antibiotic resistance. However, lab-based MIC determination can be time-consuming and suffers from low reproducibility, and interpretation as sensitive or resistant relies on guidelines which change over time. Genome sequencing and machine learning promise to allow in silico MIC prediction as an alternative approach which overcomes some of these difficulties, albeit the interpretation of MIC is still needed. Nevertheless, precisely how we should handle MIC data when dealing with predictive models remains unclear, since they are measured semi-quantitatively, with varying resolution, and are typically also left- and right-censored within varying ranges. We therefore investigated genome-based prediction of MICs in the pathogen Klebsiella pneumoniae using 4367 genomes with both simulated semi-quantitative traits and real MICs. As we were focused on clinical interpretation, we used interpretable rather than black-box machine learning models, namely, Elastic Net, Random Forests, and linear mixed models. Simulated traits were generated accounting for oligogenic, polygenic, and homoplastic genetic effects with different levels of heritability. Then we assessed how model prediction accuracy was affected when MICs were framed as regression and classification. Our results showed that treating the MICs differently depending on the number of concentration levels of antibiotic available was the most promising learning strategy. Specifically, to optimise both prediction accuracy and inference of the correct causal variants, we recommend considering the MICs as continuous and framing the learning problem as a regression when the number of observed antibiotic concentration levels is large, whereas with a smaller number of concentration levels they should be treated as a categorical variable and the learning problem should be framed as a classification. Our findings also underline how predictive models can be improved when prior biological knowledge is taken into account, due to the varying genetic architecture of each antibiotic resistance trait. Finally, we emphasise that incrementing the population database is pivotal for the future clinical implementation of these models to support routine machine-learning based diagnostics.

Host association and intracellularity evolved multiple times independently in the Rickettsiales.

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.

Tick-Borne Encephalitis, Lombardy, Italy.

Tick-borne encephalitis was limited to northeast portions of Italy. We report in Lombardy, a populous region in the northwest, a chamois displaying clinical signs of tickborne encephalitis virus that had multiple virus-positive ticks attached, as well as a symptomatic man. Further, we show serologic evidence of viral circulation in the area.

Concomitant Resistance to Cefiderocol and Ceftazidime/Avibactam in Two Carbapenemase-Producing Klebsiella pneumoniae Isolates from Two Lung Transplant Patients.

In this study, we present two cases of Klebsiella pneumoniae, one KPC-33- and one NDM-1-producing, showing resistance to cefiderocol and ceftazidime/avibactam, collected in the intensive care unit of a hospital in Northern Italy from two patients who had recently undergone lung transplantation. Whole-genome sequencing was performed to investigate the molecular features of these strains.

High-throughput screening of pathogens in Ixodes ricinus removed from hosts in Lombardy, northern Italy.

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.

Population structure and antimicrobial resistance among Klebsiella isolates sampled from human, animal, and environmental sources in Ghana: a cross-sectional genomic One Health study.

BACKGROUND: One Health approaches to address the increasing threat of antimicrobial resistance (AMR) are gaining attention. However, data on the distribution and movement of bacteria and their AMR-associated genes between clinical and non-clinical sources are scarce, especially from low-income and middle-income countries. We aimed to analyse Klebsiella isolates from various sources in Ghana and compare the prevalence of AMR with datasets from two other countries. METHODS: We conducted a cross-sectional genomic One Health study. Multiple clinical, environmental, and animal sources were sampled from 78 locations (eg, hospitals, residential areas, and farms) in and around Tamale, Ghana. Clinical samples were collected through routine screening and in cases of suspected infection between March 15 and Sept 15, 2019, and samples from the wider environment were collected during a dedicated sampling effort between the dates of Aug 19, 2018, and Sept 26, 2019. Sampling locations were approximately evenly distributed from the centre of the city and steadily outwards to capture both rural and urban locations. Samples with positive growth for Klebsiella were included. Isolates of Klebsiella were obtained from the samples using Simmons citrate agar medium and characterised by antimicrobial susceptibility testing and whole-genome sequencing. A comparative analysis with Klebsiella population surveys from Pavia, Italy, and Tromsø, Norway, was performed. AMR-associated and virulence genes were detected, and the population distribution of these genes was studied. FINDINGS: Of 957 samples collected around Tamale, Ghana, 620 were positive for Klebsiella spp. 573 Klebsiella isolates were successfully sequenced, of which 370 were Klebsiella pneumoniae. Only two hospital isolates were carbapenem-resistant. Extended-spectrum ß-lactamase (ESBL) genes were relatively common among the Ghanaian clinical isolates but rare in the environmental samples. Prevalence of ESBL genes in human-hospital disease samples was 64% (14 of 22 isolates) in Ghana and 44% (four of nine isolates) in Italy, and prevalence in human-hospital carriage samples was 7% (eight of 107) in Ghana and 13% (54 of 428) in Italy; the prevalence was higher in human-hospital disease samples than in human-hospital carriage samples in both countries, and prevalence across both samples in both countries was higher than in Norway. Ghanaian isolates showed evidence of high recombination rates (recombination events compared with point mutations [r/m] 9·455) and a considerable accessory gene overlap with isolates from Italy and Norway. INTERPRETATION: Although several AMR-associated gene classes were observed relatively frequently in non-clinical sources, ESBL, carbapenemase, and virulence genes were predominantly present only in hospital samples. These results suggest that interventions should be focused on clinical settings to have the greatest effect on the prevalence and dissemination of AMR-associated genes. FUNDING: European Research Council (742158), Academy of Finland EuroHPC grant, Trond Mohn Foundation (BATTALION grant), and Wellcome Trust.

P-DOR, an easy-to-use pipeline to reconstruct bacterial outbreaks using genomics.

SUMMARY: Bacterial Healthcare-Associated Infections (HAIs) are a major threat worldwide, which can be counteracted by establishing effective infection control measures, guided by constant surveillance and timely epidemiological investigations. Genomics is crucial in modern epidemiology but lacks standard methods and user-friendly software, accessible to users without a strong bioinformatics proficiency. To overcome these issues we developed P-DOR, a novel tool for rapid bacterial outbreak characterization. P-DOR accepts genome assemblies as input, it automatically selects a background of publicly available genomes using k-mer distances and adds it to the analysis dataset before inferring a Single-Nucleotide Polymorphism (SNP)-based phylogeny. Epidemiological clusters are identified considering the phylogenetic tree topology and SNP distances. By analyzing the SNP-distance distribution, the user can gauge the correct threshold. Patient metadata can be inputted as well, to provide a spatio-temporal representation of the outbreak. The entire pipeline is fast and scalable and can be also run on low-end computers. AVAILABILITY AND IMPLEMENTATION: P-DOR is implemented in Python3 and R and can be installed using conda environments. It is available from GitHub https://github.com/SteMIDIfactory/P-DOR under the GPL-3.0 license.

Three-dimensional images reveal the impact of the endosymbiont Midichloria mitochondrii on the host mitochondria.

The hard tick, Ixodes ricinus, a main Lyme disease vector, harbors an intracellular bacterial endosymbiont. Midichloria mitochondrii is maternally inherited and resides in the mitochondria of I. ricinus oocytes, but the consequences of this endosymbiosis are not well understood. Here, we provide 3D images of wild-type and aposymbiotic I. ricinus oocytes generated with focused ion beam-scanning electron microscopy. Quantitative image analyses of endosymbionts and oocyte mitochondria at different maturation stages show that the populations of both mitochondrion-associated bacteria and bacterium-hosting mitochondria increase upon vitellogenisation, and that mitochondria can host multiple bacteria in later stages. Three-dimensional reconstructions show symbiosis-dependent morphologies of mitochondria and demonstrate complete M. mitochondrii inclusion inside a mitochondrion. Cytoplasmic endosymbiont located close to mitochondria are not oriented towards the mitochondria, suggesting that bacterial recolonization is unlikely. We further demonstrate individual globular-shaped mitochondria in the wild type oocytes, while aposymbiotic oocytes only contain a mitochondrial network. In summary, our study suggests that M. mitochondrii modulates mitochondrial fragmentation in oogenesis possibly affecting organelle function and ensuring its presence over generations.

The evolution of intramitochondriality in Midichloria bacteria.

Midichloria spp. are intracellular bacterial symbionts of ticks. Representatives of this genus colonise mitochondria in the cells of their hosts. To shed light on this unique interaction we evaluated the presence of an intramitochondrial localization for three Midichloria in the respective tick host species and generated eight high-quality draft genomes and one closed genome, showing that this trait is non-monophyletic, either due to losses or multiple acquisitions. Comparative genomics supports the first hypothesis, as the genomes of non-mitochondrial symbionts are reduced subsets of those capable of colonising the organelles. We detect genomic signatures of mitochondrial tropism, including the differential presence of type IV secretion system and flagellum, which could allow the secretion of unique effectors and/or direct interaction with mitochondria. Other genes, including adhesion molecules, proteins involved in actin polymerisation, cell wall and outer membrane proteins, are only present in mitochondrial symbionts. The bacteria could use these to manipulate host structures, including mitochondrial membranes, to fuse with the organelles or manipulate the mitochondrial network.

Antigen discovery by bioinformatics analysis and peptide microarray for the diagnosis of cystic echinococcosis.

BACKGROUND: Cystic echinococcosis (CE), caused by Echinococcus granulosus sensu lato, is a neglected zoonosis. Its diagnosis relies on imaging, supported by serology, while only imaging is useful for staging and follow-up. Since diagnostic tools and expertise are not widely available, new accurate and easily implementable assays for the diagnosis and follow-up of CE are highly needed. METHODOLOGY/PRINCIPAL FINDINGS: We aimed to identify new E. granulosus antigens through a bioinformatics selection applied to the parasite genome, followed by peptide microarray screening and validation in ELISA, using independent panels of sera from patients with hepatic CE and clinically relevant controls. From 950 proteins selected in silico, 2,379 peptides were evaluated by microarray for IgG reactivity and eight candidates selected for validation. Reactivity to one peptide was significantly higher in the CE group (p = 0.044), but had suboptimal diagnostic accuracy. CONCLUSIONS/SIGNIFICANCE: Here we performed bioinformatics analysis and peptide microarray for antigen discovery, useful for the diagnosis of CE. Eight candidates were selected and validated. Reactivity to one peptide associated to CE but had suboptimal diagnostic accuracy. Importantly, the database developed in this study may be used to identify other antigenic candidates for CE diagnosis and follow-up.

The COVID-19 Pandemic Sparked Off a Large-Scale Outbreak of Carbapenem-Resistant Acinetobacter baumannii from the Endemic Strains at an Italian Hospital.

Acinetobacter baumannii is a nosocomial pathogen that poses a serious threat due to the rise of incidence of multidrug-resistant (MDR) strains. During the COVID-19 pandemic, MDR A. baumannii clones have caused several outbreaks worldwide. Here, we describe a detailed investigation of an MDR A. baumannii outbreak that occurred at Policlinico San Matteo (Pavia, Italy). A total of 96 A. baumannii strains, isolated between January and July 2020 from 41 inpatients (both SARS-CoV-2 positive and negative) in different wards, were characterized by phenotypic and genomic analyses combining Illumina and Nanopore sequencing. Antibiotic susceptibility testing revealed that all isolates were resistant to carbapenems, and the sequence analysis attributed this to the carbapenemase gene blaOXA-23. Virulence factor screening unveiled that all strains carried determinants for biofilm formation, while plasmid analysis revealed the presence of two plasmids, one of which was ~100 kbp long and encoded a phage sequence. A core genome-based phylogeny was inferred to integrate outbreak strain genomes with background genomes from public databases and the local surveillance program. All strains belonged to the globally disseminated sequence type 2 (ST2) clone and were mainly divided into two clades. Isolates from the outbreak clustered with surveillance isolates from 2019, suggesting that the outbreak was caused by two strains that were already circulating in the hospital before the start of the pandemic. The intensive spread of A. baumannii in the hospital was enhanced by the extreme emergency situation of the first COVID-19 pandemic wave that resulted in reduced attention to infection prevention and control practices. IMPORTANCE The COVID-19 pandemic, especially during the first wave, posed a great challenge to the hospital management and generally promoted nosocomial pathogen dissemination. MDR A. baumannii can easily spread and persist for a long time on surfaces, causing outbreaks in health care settings. Infection prevention and control practices, epidemiological surveillance, and microbiological screening are fundamental in order to control such outbreaks. Here, we sequenced the genomes of 96 isolates from an outbreak of MDR A. baumannii strains using both short- and long-read technology in order to reconstruct the outbreak events in fine detail. The sequence data demonstrated that two endemic clones of MDR A. baumannii were the source of this large hospital outbreak during the first COVID-19 pandemic wave, confirming the effect of COVID-19 emergency disrupting the protection provided by the use of the standard prevention procedures.

Characterization of the bacterial microbiome of Swedish ticks through 16S rRNA amplicon sequencing of whole ticks and of individual tick organs.

BACKGROUND: The composition of the microbial flora associated with ixodid ticks has been studied in several species, revealing the importance of geographical origin, developmental stage(s) and feeding status of the tick, as well as substantial differences between tissues and organs. Studying the microbiome in the correct context and scale is therefore necessary for understanding the interactions between tick-borne pathogens and other microorganisms as well as other aspects of tick biology. METHODS: In the present study the microbial flora of whole Ixodes ricinus, I. persulcatus and I. trianguliceps ticks were analyzed with 16S rRNA amplicon sequencing. Additionally, tick organs (midguts, Malpighian tubules, ovaries, salivary glands) from flat and engorged I. ricinus female ticks were examined with the same methodology. RESULTS: The most abundant bacteria belonged to the group of Proteobacteria (Cand. Midichloria mitochondrii and Cand. Lariskella). 16S amplicon sequencing of dissected tick organs provided more information on the diversity of I. ricinus-associated microbial flora, especially when organs were collected from engorged ticks. Bacterial genera significantly associated with tick feeding status as well as genera associated with the presence of tick-borne pathogens were identified. CONCLUSIONS: These results contribute to the knowledge of microbial flora associated with ixodid ticks in their northernmost distribution limit in Europe and opens new perspectives for other investigations on the function of these bacteria, including those using other approaches like in vitro cultivation and in vitro models.

A New Benzothiazolthiazolidine Derivative, 11726172, Is Active In Vitro, In Vivo, and against Nonreplicating Cells of Mycobacterium tuberculosis.

Tuberculosis (TB) still poses a global menace as one of the deadliest infectious diseases. A quarter of the human population is indeed latently infected with Mycobacterium tuberculosis. People with latent infection have a 5 to 10% lifetime risk of becoming ill with TB, representing a reservoir for TB active infection. This is a worrisome problem to overcome in the case of relapse; unfortunately, few drugs are effective against nonreplicating M. tuberculosis cells. Novel strategies to combat TB, including its latent form, are urgently needed. In response to the lack of new effective drugs and after screening about 500 original chemical molecules, we selected a compound, 11726172, that is endowed with potent antitubercular activity against M. tuberculosis both in vitro and in vivo and importantly also against dormant nonculturable bacilli. We also investigated the mechanism of action of 11726172 by applying a multidisciplinary approach, including transcriptomic, labeled metabolomic, biochemical, and microbiological procedures. Our results represent an important step forward in the development of a new antitubercular compound with a novel mechanism of action active against latent bacilli. IMPORTANCE The discontinuation of TB services due to COVID-19 causes concern about a future resurgence of TB, also considering that latent infection affects a high number of people worldwide. To combat this situation, the identification of antitubercular compounds targeting Mycobacterium tuberculosis through novel mechanisms of action is necessary. These compounds should be active against not only replicating bacteria cells but also nonreplicating cells to limit the reservoir of latently infected people on which the bacterium can rely to spread after reactivation.

A large-scale genomic snapshot of Klebsiella spp. isolates in Northern Italy reveals limited transmission between clinical and non-clinical settings.

The Klebsiella group, found in humans, livestock, plants, soil, water and wild animals, is genetically and ecologically diverse. Many species are opportunistic pathogens and can harbour diverse classes of antimicrobial resistance genes. Healthcare-associated Klebsiella pneumoniae clones that are non-susceptible to carbapenems can spread rapidly, representing a high public health burden. Here we report an analysis of 3,482 genome sequences representing 15 Klebsiella species sampled over a 17-month period from a wide range of clinical, community, animal and environmental settings in and around the Italian city of Pavia. Northern Italy is a hotspot for hospital-acquired carbapenem non-susceptible Klebsiella and thus a pertinent setting to examine the overlap between isolates in clinical and non-clinical settings. We found no genotypic or phenotypic evidence for non-susceptibility to carbapenems outside the clinical environment. Although we noted occasional transmission between clinical and non-clinical settings, our data point to a limited role of animal and environmental reservoirs in the human acquisition of Klebsiella spp. We also provide a detailed genus-wide view of genomic diversity and population structure, including the identification of new groups.

CodY Is a Global Transcriptional Regulator Required for Virulence in Group B Streptococcus.

Group B Streptococcus (GBS) is a Gram-positive bacterium able to switch from a harmless commensal of healthy adults to a pathogen responsible for invasive infections in neonates. The signals and regulatory mechanisms governing this transition are still largely unknown. CodY is a highly conserved global transcriptional regulator that links nutrient availability to the regulation of major metabolic and virulence pathways in low-G+C Gram-positive bacteria. In this work, we investigated the role of CodY in BM110, a GBS strain representative of a hypervirulent lineage associated with the majority of neonatal meningitis. Deletion of codY resulted in a reduced ability of the mutant strain to cause infections in neonatal and adult animal models. The observed decreased in vivo lethality was associated with an impaired ability of the mutant to persist in the blood, spread to distant organs, and cross the blood-brain barrier. Notably, the codY null mutant showed reduced adhesion to monolayers of human epithelial cells in vitro and an increased ability to form biofilms, a phenotype associated with strains able to asymptomatically colonize the host. RNA-seq analysis showed that CodY controls about 13% of the genome of GBS, acting mainly as a repressor of genes involved in amino acid transport and metabolism and encoding surface anchored proteins, including the virulence factor Srr2. CodY activity was shown to be dependent on the availability of branched-chain amino acids, which are the universal cofactors of this regulator. These results highlight a key role for CodY in the control of GBS virulence.

OCCURRENCE OF EUSTRONGYLIDES EXCISUS (NEMATODA:DIOCTOPHYMATIDAE) IN EUROPEAN PERCH (PERCA FLUVIATILIS) AND GREAT CORMORANT (PHALACROCORAX CARBO) IN LAKE ANNONE, NORTHERN ITALY.

The genus Eustrongylides includes zoonotic nematodes that infect fish species and fish-eating birds of freshwater ecosystems. This study aimed to evaluate the occurrence of Eustrongylides in the paratenic host Perca fluviatilis (European perch) and in the definitive host, Phalacrocorax carbo sinensis (great cormorant), in Lake Annone, a shallow eutrophic lake located in the pre-mountainous area of the Alps in northwest Italy where wintering cormorants coexist with new breeding colonies. A total of 114 European perch and 48 cormorants were examined for the occurrence of Eustrongylides. All parasites collected were identified with microscopic examination and molecular analysis. Overall, 11 specimens of European perch (9.6%) and 13 individuals of cormorants (27%) harbored nematodes identified as fourth-stage larvae and adults of Eustrongylides excisus. The observed prevalence of Eustrongylides spp. appears to be intermediate between the higher values in cormorant breeding areas in northern Europe and the lower prevalence observed in their wintering sites in southernmost Europe. Considering the eutrophication status of freshwater ecosystems and the increasing population of the cormorants, Eustrongylides has an increasing potential range of dispersion in Europe, including Italy; thus an extensive surveillance should be carried out, especially given the zoonotic potential of this nematode.

Isolation of a Colistin-Susceptible MDR Pantoea calida Harboring the mcr-9 Gene Suggests the Silent Spread of the Resistance Factor.

Pantoea spp. are bacteria that are often detected in the environment and as symbionts of arthropods. They sporadically cause infections in humans and recently extended spectrum beta-lactamases (ESBL)- and carbapenemase-producing strains have started to emerge. In this study, we report the isolation and the complete genome sequence of a strain of Pantoea calida encoding the colistin-resistance gene mcr-9. The strain was isolated from a preterm newborn in a neonatal pathology ward. On clinical examination, his vital signs were normal and blood culture was negative. Rectal swab screening for ESBL-producing Enterobacterales allowed to isolate the bacterium, and a complete genome was obtained using both short and long read sequencing. The mcr-9 gene was found to be encoded on a IncHI2 superplasmid, which confers resistance to six classes of antibiotics, including beta lactams (ESBL). Despite the presence of mcr-9, the isolate retains susceptibility to colistin, which could be explained by the absence of compatible regulatory genes (qseBC) from the genome. The presence of the resistance gene is undetectable with the routine clinical procedures, that is, phenotypic tests. This suggests that a silent spread might be ongoing in the ward. To our knowledge, this is the first description of an MDR P. calida and of a Pantoea spp. encoding any mobile colistin resistance gene.

Molecular identification of Contracaecum rudolphii A and B (Nematoda: Anisakidae) from cormorants collected in a freshwater ecosystem of the pre-alpine area in Northern Italy.

Contracaecum rudolphii (s.l.) is a complex of sibling species with different genetic structure and ecological preference. This study reports the presence of specimens of Contracaecum rudolphii (s.l.) from sedentary and wintering cormorants (Phalacrocorax carbo sinensis) from the pre-mountain area of the Alps in Northern Italy, an important crossroads for most of the bird migration routes. A total of 48 specimens of cormorants collected from two adjacent freshwater habitats were analysed and C. rudolphii nematodes were retrieved in 100% of the examined specimens. A subsamples of 115 C. rudolphii individuals were genetically characterized and found to belong to the sibling species C. rudolphii B (n = 90) and C. rudolphii A (n = 25). C. rudolphii B were retrieved from both locations and included adults as well as larvae, while only adults of C. rudolphii A were detected, and in just one location. As expected for a freshwater environment, C. rudolphii B constitutes the largest sibling fraction, indicating that this likely is the endemic species, while cormorants originating from the breeding brackish lagoons and marine coastal environments of central and northern Europe could have brought C. rudolphii A from their breeding sites or migration stopovers.

'Candidatus Gromoviella agglomerans', a novel intracellular Holosporaceae parasite of the ciliate Paramecium showing marked genome reduction.

Holosporales are an alphaproteobacterial lineage encompassing bacteria obligatorily associated with multiple diverse eukaryotes. For most representatives, little is known on the interactions with their hosts. In this study, we characterized a novel Holosporales symbiont of the ciliate Paramecium polycaryum. This bacterium inhabits the host cytoplasm, frequently forming quite large aggregates. Possibly due to such aggregates, host cells sometimes displayed lethal division defects. The symbiont was also able to experimentally stably infect another Paramecium polycaryum strain. The bacterium is phylogenetically related with symbionts of other ciliates and diplonemids, forming a putatively fast-evolving clade within the family Holosporaceae. Similarly to many close relatives, it presents a very small genome (<600 kbp), and, accordingly, a limited predicted metabolism, implying a heavy dependence on Paramecium, thanks also to some specialized membrane transporters. Characterized features, including the presence of specific secretion systems, are overall suggestive of a mild parasitic effect on the host. From an evolutionary perspective, a potential ancestral trend towards pronounced genome reduction and possibly linked to parasitism could be inferred, at least among fast-evolving Holosporaceae, with some lineage-specific traits. Interestingly, similar convergent features could be observed in other host-associated lineages, in particular Rickettsiales among Alphaproteobacteria.