Characterization of Rosculus vilicus sp. nov., a rhizarian amoeba interacting with Mycobacterium avium subsp. paratuberculosis

Free-living amoebae are described as potential reservoirs for pathogenic bacteria in the environment. It has been hypothesized that this might be the case for Mycobacterium avium subsp. paratuberculosis, the bacterium responsible for paratuberculosis. In a previous work, we isolated an amoeba from a water sample in the environment of infected cattle and showed that this amoeba was associated with Mycobacterium avium subsp. paratuberculosis. While a partial 18S rRNA gene has allowed us to suggest that this amoeba was Rosculus-like, at that time we were not able to sub-cultivate it. In the present study, we succeeded in cultivating this strain at 20–25°C. This amoeba is among the smallest (5–7 μm) described. The sequencing of the whole genome allowed us to extract the full 18S rRNA gene and propose this strain as a new species of the Rosculus genus, i.e., R. vilicus. Of note, the mitochondrial genome is particularly large (184,954 bp). Finally, we showed that this amoeba was able to phagocyte Mycobacterium avium subsp. paratuberculosis and that the bacterium was still observed within amoebae after at least 3 days. In conclusion, we characterized a new environmental amoeba species at the cellular and genome level that was able to interact with Mycobacterium avium subsp. paratuberculosis. As a result, R. vilicus is a potential candidate as environmental reservoir for Mycobacterium avium subsp. paratuberculosis but further experiments are needed to test this hypothesis.

1 Introduction

Free-living amoebae (FLA) are single-celled eukaryotes, omnipresent in water and soil environments (Rodríguez-Zaragoza, 1994; Samba-Louaka et al., 2019). These protists share common morphological features, notably pseudopods that allow them to move upon surfaces and to thrive by phagocytosis. Even though they are likely to play an important role in microbial ecology due to their ability to ingest bacteria as their main food sources (Jürgens and Matz, 2002), FLA have yet to be adequately described.

FLA constitutes a highly polyphyletic group of microbial eukaryotes and can be affiliated not only to Amoebozoa, the sole supergroup exclusively represented by amoeboid microorganisms, but also to Opisthokonta, Heterolobosea and Rhizaria. Among those three large eukaryotic groups, FLA consist mostly in small clades of amoebae, featuring a varied array of cell morphologies (Bass et al., 2009; Samba-Louaka et al., 2019). Among Rhizaria, the most notorious amoebae representatives are Vampyrellids, a group of predatory amoebae found in freshwater and marine environments. They can feed on bacteria and other eukaryotic preys (Hess and Suthaus, 2022). Additionally, other lesser known rhizarian amoebae were recently described. Most notably, molecular analyses highlighted the existence of a potentially very diverse clade, the Sainouroidea. The latter encompasses several genera of small-sized amoeba, such as Rosculus, which can be found in cattle, snake, fish feces, freshwater and soil (Schuler et al., 2018).

While FLA ingest bacteria by phagocytosis, some bacteria can resist phagocytosis and, consequently, might more specifically be able to resist phagocytosis by the immune cells of animals. This is due to a high level of conservation between phagocytosis mechanisms of FLA and animal cells (Escoll et al., 2013). This resistance was initially described for Legionella pneumophila, and many similar observations were made later on, involving diverse bacteria such as Mycobacterium spp., including Mycobacterium avium subsp. paratuberculosis (MAP) the bacterium responsible for paratuberculosis (Drancourt, 2014; Claeys and Robinson, 2018). In addition, amoeba-resistant bacteria might become more virulent insofar as they are adapted to the phagocyte environment and express virulence genes, thereby exemplifying the concept of coincidental evolution (Primm et al., 2004; Sun et al., 2018). Consequently, interest in FLA was revived by their capacity to become reservoirs of pathogenic bacteria in the environment.

It has been hypothesized that amoebae might be environmental reservoirs of Mycobacteria (Salah et al., 2009). To test this, we tried to isolate amoebae from water in the environment of infected cattle and look for the association between MAP, the bacterium responsible for paratuberculosis, and indigenous amoebae (Samba-Louaka et al., 2018). One amoeba was detected and identified, via 18S rRNA partial gene sequencing, as a Rosculus-like amoeba belonging to the supergroup Rhizaria. However, we were not able at that time to maintain the growth of this strain.

In this work, we ultimately identified optimal conditions for growing and isolating the strain associated with MAP. We fully sequenced the genome of this amoeba and analysis of the 18S rRNA gene sequence allowed us to confirm this amoeba among the Rosculus genus. Finally, an interaction study was performed to test the ability of this amoeba to phagocytose MAP.

2 Materials and methods

2.1 Isolation and cultivation of amoeba

This amoeba was previously isolated from an environmental water sample from the environment of infected cattle (Samba-Louaka et al., 2018). A sample (cryotube) frozen at that time was used to obtain a pure culture by critical dilution. Then, the cultivation was maintained within PAS buffer (Na₃C₆H₅O₇ 1 g/L, MgSO₄:7H₂O 0.4 M, CaCL₂:2H₂O 0.05 M, Na₂HPO₄:2H₂O 0.25 M, KH₂PO₄ 0.25 M, pH 6.5) containing E. coli (5.10⁸ bacteria/mL) as nutrient.

To determine the optimal growing temperatures, Rosculus vilicus was cultured at 20, 25, 30 or 37°C for 24 h. Cells were detached and counted after an incubation of 3, 6 and 24 h. The experiment was repeated three times.

2.2 DNA extraction and genome sequencing

From pure culture, DNA was extracted in order to sequence this amoeba genome with two sequencing technologies, short and long-read approaches. Briefly, lysis buffer was added directly on adherent R. vilicus. Then, DNA was extracted with Blood & Tissue kit (Qiagen) following the manufacturer’s recommendations and quantified using the Qubit fluorometer (Thermo Fischer Scientific).

For the long-read approach, 3 μg of DNA was used for library preparation using the Oxford Nanopore Technologies ligation kit (SQK-LSK109) according to the manufacturer’s protocol with the Long Fragment Buffer (LFB). The library was then loaded on a flow cell (R9.4.1) and sequenced with MinION Mk1C. Base calling was done using Guppy v6.5.7 in high-accuracy mode. Simultaneously, genomic DNA was sent to SeqCenter (Pittsburgh, PA, USA) for library preparation and short read sequencing on an Illumina NovaSeq 6000. Demultiplexing, quality control and adapter trimming were performed with BCL convert software at SeqCenter.

2.3 Phylogenetic and genome analysis

We strived to follow a hybrid-based approach inspired by recently published protocols (Wick et al., 2023) to assemble a high-quality draft genome. For this, initial assembly was performed using Flye v. 2.9.1 (Kolmogorov et al., 2019) with long reads only, displaying quality score (Q) > 9 after base calling. A step of self-polishing was then implemented using Medaka v 1.7.2. Illumina reads were trimmed based on Q > 30 using Cutadapt v 4.4 (Martin, 2011), and for additional polishing of contigs using Polypolish v 0.5 (Wick and Holt, 2022). The resulting, fully polished assembly was then binned using Metabat2 V 2.15 (Kang et al., 2019). Manual inspections of bins were performed to ensure proper binning of sequences. For inferring phylogeny, Barrnap v 1.2.2 was used to identify and extract gene sequences coding for ribosomal RNA.¹ Full-length 18S rRNA gene sequence was used to reconstruct phylogeny using IQ-TREE v 2.0.3 (Minh et al., 2020). The optimal substitution model was chosen based on the output of ModelFinder (Kalyaanamoorthy et al., 2017), according to the best Bayesian information criterion score. The robustness of nodes was tested by 1,000 iterations of conventional bootstraps and 1,000 iterations of the Shimodaira–Hasegawa approximate likelihood ratio test.

Mitochondrial genome was extracted based on joint analysis of graph assemblies, contig coverage and GC content. After manual and putative identification, the supposedly mitochondrial genome was further confirmed as such using multiple dedicated annotation tools, i.e., MITOS 2 (Donath et al., 2019), Mfannot² and RNAweasel³ to identify mitochondrial proteins and RNAs. The nuclear genome and mitochondrial genome data were deposited at NCBI under the BioProject ID PRJNA1029052.

2.4 Morphological characterization of amoeba strain

Trophozoites and cyst-like forms from R. vilicus were observed by microscopy (Olympus IX73 inverted microscope, 1,000x magnification) using a differential interference contrast. All sizes were measured using ImageJ (Schneider et al., 2012). The morphological comparison between trophozoite and cyst-like forms was also assessed by electron microscopy. For scanning electron microscopy (SEM), cells were cryofixed in nitrogen with VCM (Leica), sublimed and 4 nm wolfram coated in ACE 600 (Leica). Each step was performed under vacuum and samples were transferred from one station to another using VCT 500 (Leica). Samples were observed at −100°C in scanning electron microscope Teneo VolumeScope (Thermo Fisher) to 3.5 kV and 10 kV.

For transmission electron microscopy (TEM), amoebae were fixed with 2.5% glutaraldehyde and 2% osmium tetroxide. Dehydration was carried out using different concentrations of acetone. Epon resin were used for impregnation and embedding process. After cutting with ultramicrotome UC6 (Leica), samples were contrasted with uranyl acetate and lead citrate. Samples observations were proceeded on transmission electron microscope (JEOL 1010) at 80 keV. Images were acquired with Quemesa camera (Olympus).

2.5 Infection experiments

To study the permissiveness of R. vilicus to MAP infection, an infection experiment was set up using a MAP K10-GFP strain (Li et al., 2005). R. vilicus (2.10⁵ cells) were added in each well of a 6-well plate in a 2 mL PAS buffer and incubated for 1 h at 20°C. MAP clumps were disaggregated with 3 passages in 27 gauges needle, centrifuged for 5 min at 200 g and added to R. vilicus at MOI 10 during 2 h at 20°C. Amikacin (100 μg/mL) was added for 2 h to remove extracellular mycobacteria before the washing step. Finally, the infected cells were incubated in PAS buffer added with E. coli (4.10⁸ bacteria/mL) for 3 days at 20°C. Colocalization of MAP K10-GFP with R. vilicus was observed under fluorescence microscope (Olympus IX51 inverted fluorescence microscope).

3 Results

3.1 Rosculus vilicus grows preferentially between 20 and 25°C

We previously isolated a Rosculus-like amoeba that we cultured at 30°C but we ultimately lost over subcultures, suggesting at that time that growth conditions were not appropriate for this isolate. To determine the optimal temperature, we cultivated this amoeba at 20, 25, 30, and 37°C. As demonstrated in this study, the optimal growth temperatures ranged from 20°C to 25°C. Indeed, at 30°C the development was inhibited while at 37°C, the total number of cells decreased, suggesting cell death (Figure 1).

FIGURE 1

Figure 1. R. vilicus optimal growth was at 20°C and 25°C. Growth comparison of R. vilicus was tested by counting number of cells per ml at different incubation temperatures. The results represent three independent counts, and errors bars represent the standard error of the mean (±SEM). Statistical analysis was performed using Kruskal-Wallis test and uncorrected Dunn’s test (*p < 0.05).

3.2 Rosculus vilicus shares morphological features with Rosculus genus

To characterize this isolate, microscopic studies were performed. Differential interferential contrast micrographs showed a monopodial lobose amoeba with a well-defined frontal hyaline area. Trophozoite cells usually displayed rapid eruptive locomotion with some lateral eruption (Figures 2A–C). An uroidal structure could be observed on some trophozoites, the opposite of that found in the hyaline area (Figures 2A–C). Average length was estimated at 7.1 μm (range of 4.2–13.8 μm) and average width was 4.8 μm (range of 2.7–12.1 μm). During observations, some amoebae adopted a floating form, more rounded than trophozoite (6.0 μm average length and 4.9 μm average width), harboring a long stem, about 18.9 μm on average (10.2–31.4 μm) (Figures 2F,G). All amoeba cells contained a large vacuole and several granularities. Given the ability of amoebae to differentiate into cyst forms, we incubated our amoebae in an encysting buffer for 4 days. We observed rounded cells with an average width of 4.8 μm (Figures 2D,E). The observation of trophozoites and cyst-like forms under a scanning electron microscope confirmed their estimated size (Figures 3A–D). The cell surface of cyst-like forms displayed one ostiole (Figure 3D).

FIGURE 2

Figure 2. R. vilicus adopts different forms. Trophozoite forms (A–C), rounded cells (D,E) and floating forms with a stem (F,G). Bar length represents 10 μm.

FIGURE 3

Figure 3. R. vilicus external morphology viewed by scanning electron microscopy. Trophozoite forms (A,B) and cyst-like forms (C,D).

Transmission electron micrographs enabled us to observe the ultrastructure of both trophozoites and cyst-like forms. The cell membrane in both trophozoites and cyst-like forms consisted of a thin single membrane (Figures 4A,B). Regarding trophozoites, the cytoplasm contained a large nucleus bordered by a dense chromatin, large mitochondria with discoidal mitochondrial cristae and several vacuoles (Figure 4A). Interestingly, in encystment conditions, the extracellular medium contained free bacteria and numerous vesicle-like structures filled with bacteria (Figure 4C).

FIGURE 4

Figure 4. R. vilicus internal ultrastructure viewed by transmission electron microscopy; trophozoite form (A) cyst-like form (B) and extracellular medium observed for encystment condition containing vesicles and bacteria (C). N, nucleus; M, mitochondria; V, vesicle containing bacteria; B, bacteria.

3.3 Rosculus vilicus is a new species with a large mitochondrial genome

As no genome was described in the Rosculus genus and even in the Sainouroidea clade, we decided to further characterize the isolate on this aspect. Whole genome sequencing was undertaken, using a hybrid, short and long-read sequencing approach. Based on this, the assembled genome of this isolate was predicted to be 40.8 Mb long, showing 37.7% GC enrichment (Table 1). While E. coli genome (used as a food source) was recovered in the sequencing data, no other bacterial genome could be retrieved. In this hybrid assembly, particular focus was given to recovery of the mitochondrial genome. One contiguous, highly covered sequence was recovered from the assembly, consisting in a 184,954 bp long sequence, with overall 24.9% GC enrichment. While it was unusually long, even when compared with other mitochondrial genomes from other related protists from the Cercozoa group (Wideman et al., 2020) we are confident that it corresponds to a bona fide mitochondrion. Indeed, mitochondrial genome annotation revealed gene repertoires similar to mitochondria from other cercozoans, while showing even coverage and reproducible assembly features, even with short reads only (Supplementary Figure S1).

TABLE 1

Table 1. Statistics and features of R. vilicus nuclear and mitochondrial genome assemblies.

This sequencing effort brought into focus the first draft genome for the whole Sainouroidea. Because of this, however, comparative analyses of gene and genome sequences was clearly hampered, and consequently restricted to widely used marker genes. Given the resulting paucity of genomic information on this clade (and on Cercozoa in general), in order to investigate the taxonomy and phylogenetic positioning of our isolate, we focused on 18S rRNA gene sequences. The full-length 18S rRNA gene sequence (2,128 bp) used for BLASTn search showed its highest match with a sequence designated as ‘uncultured Cercozoa isolate 1’ (98.18% identity, 59% of query coverage). This sequence indeed corresponds to the previously published partial 18S rRNA gene sequence of this isolate (Samba-Louaka et al., 2018).

To confirm the belonging of this amoeba to the Rosculus genus, a phylogenetic tree of the 18S rRNA gene was constructed, using a set of rRNA sequences from all available Sainouroidea representatives as of June 2023. Phylogeny inferences clearly placed our isolate within the Rosculus genus (Figure 5). The relationship between the sequence from our isolate, compared to other known sequences of the Rosculus genus, justified our proposal of a new species, which we tentatively named Rosculus vilicus sp. nov., pertaining to the farm where the amoeba was isolated.

FIGURE 5

Figure 5. R. vilicus is likely a new species, closely related to R. vulgaris. Phylogeny of Sainouroidea rooted with Filoreta as outgroup based on 18S rRNA gene.

3.4 Rosculus vilicus interacts with MAP for at least a few days

Rosculus vilicus was previously isolated in association with a MAP (Samba-Louaka et al., 2018). Here, we evaluated the permissiveness of R. vilicus to the MAP K10 strain. An interaction experiment was carried out with the MAP K10-GFP strain to follow the fate of the bacteria by fluorescence microscopy. The result showed that R. vilicus co-localized with MAP K10-GFP at least 72 h post-infection (indicated by white arrows, Figures 6A–D). SomeMAP K10 were localized extracellularly (Figures 6A,B). Furthermore, video-microscopy experiments highlighted R. vilicus moving by ameboid movement with MAP K10 following the same course (Supplementary Figures S2, S3). These results suggest an association for a period of at least 3 days between R. vilicus and MAP.

FIGURE 6

Figure 6. R. vilicus co-localized with MAP K10-GFP. Micrographs with magnification x400 show R. vilicus in co-localization with MAP K10-GFP at 4 h (A) 24 h (B) 48 h (C) and 72 h post-infection (D). The interaction is highlighted with white arrows or framed and magnified three times.

4 Discussion

It is well-known that FLA can host pathogenic and non-pathogenic bacteria, thereby playing a potential role on their persistence in the environment. However, in the environment the interactions between FLA and MAP have been very poorly investigated (Drancourt, 2014).

In a previous study, aimed at isolating FLA associated with MAP in the environment of infected cattle, we isolated a strain, which we identified at the time as being somehow related to Rosculus, from a drinking trough (Samba-Louaka et al., 2018). This strain was grown for a short period of time, but lost along subcultures, most likely indicating that the growth conditions were not fully adequate at that time. Here, the first challenge was to identify the conditions for optimal and perennial cultivation of this FLA. Different incubation temperatures were tested, clearly showing an optimum at 20–25°C. This result is in accordance with the cultivation conditions described for other Rosculus species (Schuler et al., 2018). Current knowledge on Sainouroidea, and more particularly on Rosculus spp., suggests that they are amphizoic (Bass et al., 2016). Members of the Rosculus genus were repeatedly isolated from both soil and fecal matter (Schuler et al., 2018) and sometimes even directly from rectal sampling of the animal host (Hawes, 1955). This would indicate that Rosculus are not only able to grow in the environment, but also in the feces and even in the digestive system of their respective hosts. In the case of Rosculus vilicus, our inability to grow at 37°C may be surprising, as one may expect its ability to thrive at temperatures close to the bovine digestive system (Burfeind et al., 2010). The presence of complex bacterial consortia and oxygenation levels may play an important role in growth capabilities of this amoeba and would necessitate further characterization. Also, it is possible that Rosculus is present in feces but may be able to grow only in the environment when the temperature drops to 20–25°C and not in the gut at 37°C.

The morphology of R. vilicus, observed under light microscopy, displays features comparable to the Rosculus strain previously described in the literature. For example, its small size is a morphological trait found in the other Rosculus isolates described. In line with this, we observed that R. vilicus trophozoites were on average 7.1 μm long and 4.8 μm wide, while cyst-like forms were round and on average 4.8 μm in diameter. These observations correspond to the measurements of previously isolated Rosculus strains. Depending on the species, they are between 2.5 and 13 μm long and between 4.2 and 5.5 μm wide in their locomotive forms (Hawes, 1963; Bass et al., 2009; Tyml and Dyková, 2018). Trophozoite forms of Rosculus were shown to have a rapid eruptive locomotion with an elongate hyaline pseudopod. Endoplasm is described as crowded with vacuoles, resulting in a “sandy appearance” (Hawes, 1963; Bass et al., 2009). In our observations, similar morphological traits were noticed; rapid eruptive locomotion with one hyaline pseudopod and dense granuloplasm. No flagella were observed on R. vilicus, as described on the other Rosculus strains (Hawes, 1963; Bass et al., 2009; Tyml and Dyková, 2018). However, a new morphological trait, hitherto undescribed in Rosculus species, was a long stem on floating forms (Figures 3A,B). To our knowledge, it had not been described previously, and while the associated function is not obvious, we may hypothesize that adopting such a morphology may impact cell buoyancy, thereby facilitating dispersal along with the flow of liquid.

Cyst-like forms observed under light microscopy were round or oval-shaped, with refringent wall, and correspond to previously observed cystic form. Observation by TEM was more difficult because few cells were observed in our encystment conditions. The cells observed had a thin cell wall, contrary to what had been observed before, and which is generally found in amoeba encystment forms (Hawes, 1963; Bass et al., 2009; Tyml and Dyková, 2018). Despite the lack of thick cell wall, the observation of an ostiole with SEM suggests that it is really a cyst form.

As is the case with most protist clades, genomic information on FLA is dramatically lacking. This should stimulate a systemic attempt at producing genomic data for novel protists isolates (Sibbald and Archibald, 2017). In line with this intent, we proceeded with the genomic characterization of R. vilicus, which has yielded the first genomic data for the Sainouroidea and adds to the (poorly) known genomic diversity of rhizarian protists. The draft genome of R. vilicus, assembled through a combination of long and short read data, provided an estimate of its size, which is ca. 40.8 Mbp long. This genome size is comparable with other FLA from distinct groups, such as Acanthamoeba castellanii (43.8 Mbp; Matthey-Doret et al., 2022), Vermamoeba vermiformis (39.5 Mbp; Delafont et al., 2021) or Naegleria gruberi (41 Mbp, Fritz-Laylin et al., 2010). Further studies providing transcriptomic data on R. vilicus could provide an initial annotation of this genome, a prerequisite for in-depth studies of gene repertoire. Aside from this nuclear genome, our hybrid assembly approach enabled us identify and isolate the mitochondrial genome of R. vilicus, which is unusually long (185 kbp). This size estimate far exceeds the size of other mitochondrial genomes from the Cercozoa, for which recent estimates indicated sizes around 40 to 50 kbp (Wideman et al., 2020). While significantly larger than other known sequences from the same clade, mitochondrial gene repertoire of R. vilicus is highly comparable to that of other cercozoans and seems to once again prove that we have identified and correctly assembled this organellar genome (Supplementary Figure S1). To our knowledge, the characterization of R. vilicus mitochondrion represents the largest known mitochondrial genome ever described for a single-celled organism, thereby providing an interesting model for investigation of the dynamics of mitochondrial genome evolution. The novelty of this isolate should be seen as another argument in favor of the systematic implementation of genome sequencing for all novel protist isolates, as it represents a huge source of original findings.

In the process of assembling R. vilicus genome, we noted that no bacteria (apart from E. coli used as food source) was recovered. This suggests that R. vilicus does not live in symbiotic associations with other bacteria, as is the case for numerous other FLA (Samba-Louaka et al., 2019). However, one cannot rule out the possibility that the growth conditions we used, and more globally a switch from environmental to in vitro conditions, may have contributed to a loss of this association, on which we may have missed out in this study.

The strain isolated from water trough of infected cattle herds with paratuberculosis was a good candidate to verify our hypothesis that free-living amoebae could be MAP vectors. Since the first evidence the ability of mycobacteria species to survive in amoebae (Prasad and Gupta, 1978) different Mycobacterium species have been investigated regarding their interactions with phagocytic protozoan such as Acanthamoeba, Dictyostelium or Tetrahymena. Several Mycobacterium species were able to resist amoebal digestion, suggesting a role of amoebae in mycobacterial infections (Thomas and McDonnell, 2007). Most of these findings documented the involvement of amoebozoans and heteroloboseans representative in such interactions (Balczun and Scheid, 2017; Shi et al., 2021). However, recent work has demonstrated that numerous rhizarian FLA could also bear intracellular bacteria, thereby confirming that tight FLA-bacteria interactions indeed occur in all these taxonomic groups (Pohl et al., 2021; Solbach et al., 2021).

Few studies have related MAP association with amoebae. Among them, it was shown that incubation of MAP with Acanthamoeba polyphaga and Acanthamoeba castellanii, led to a decreasing MAP number in the first few post-internalization days, but increased after an extended period of several weeks (Mura et al., 2006; Whan et al., 2006; Phillips et al., 2020). Our results provide evidence that MAP can also persist for at least for 3 days in association with R. vilicus. However, it was a rare event, as only a few amoebae were associated with the MAP K10-GFP and it was not possible to distinguish if the bacteria were inside (phagocytized) or outside the amoebae. This suggests, added to previous observations using different FLA, that MAP might interact with a wide range of amoebae distantly related.

In conclusion, our study enabled to recover R. vilicus, a new species of Rosculus that grows preferentially between 20 and 25°C. Microscopic observations showed a small amoeba with morphological traits representative of the genus Rosculus. Whole genome sequencing added genomic data to the poorly known genomic diversity of rhizarian protists. Unexpectedly, this work highlighted the particularities mitochondrial genome size. R. vilicus can interact with MAP for at least a few days. This result suggests the possible role of amoebae in the persistence of MAP and contributes new evidence about the protists range that interact with MAP in environment.

4.1 Description of Rosculus vilicus sp. nov

Description of Rosculus vilicus [vi.li.cus. N.L. adj.] sp. nov.

Taxonomic assignment: Eukaryota; TSAR; Rhizaria; Cercozoa; Sainouroidea; Guttulinopsidae; Rosculus.

Diagnosis: Cercozoan, limax-shaped amoeba. Single nucleus. Found in freshwater. Trophozoites: 4.2–13.8 μm (average: 7.1 μm) in length and 2.7–12.1 μm (average: 4.8 μm) in width. Trophozoites can adopt a more rounded floating form which typically displays a 10.2–31.4 μm stem (average: 18.9 μm). The amoebae can form a cyst-like structure when starved (average: 4.8 μm). The amoebae feed on the bacteria Escherichia coli.

Type strain: Isolate AJ1 (available upon request). Isolated from cattle drinking trough located in Indre et Loire department, France. The amoebae were isolated on a non-nutrient agar plate seeded with E. coli and routinely cultivated on PAS buffer with E. coli (OD 0.5) at 20°C.

Etymology: vilicus from villa (“country house”) + icus; pertaining to the farm.

The genome size of this isolate is was predicted to be 40.8 Mb long, showing 37.7% GC enrichment. This is the first genome of Rosculus spp. and there is no genome available in related species shown in Figure 5. The sequence was deposited at NCBI under the BioProject ID PRJNA1029052.

Data availability statement

The data presented in this study are deposited in the NCBI database under accession numbers OR644361 and PRJNA1029052.

Author contributions

AJ: Conceptualization, Data curation, Writing – original draft, Writing – review & editing, Formal analysis, Investigation. VD: Writing – original draft, Writing – review & editing, Conceptualization, Data curation, Methodology, Formal analysis. J-LM: Writing – review & editing, Project administration, Funding acquisition. FB: Writing – review & editing, Conceptualization, Funding acquisition. AS-L: Writing – original draft, Writing – review & editing, Conceptualization, Data curation, Methodology, Formal analysis. YH: Conceptualization, Data curation, Funding acquisition, Supervision, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Agence Nationale de la Recherche (ANR-22-CE35-0010, « Amitub », 01-01-2027).

Acknowledgments

We thank Emile Bere (ImageUp, Université de Poitiers) for his assistance with microscopy experiments and Thibault Chautrand for his technical assistance. We also thank Bouziane Moumen from the Ebioinf bioinformatics platform CNRS UMR7267, EBI.

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.

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.2023.1324985/full#supplementary-material

Supplementary Figure S1 | The unusually large mitochondrion of R. vilicus shows a gene repertoire comparable to, and phylogenetic affinities with, other cercozoans. Analyses of mitochondrial contig coverage suggests that its unusual size does not result from artefactual assembly (A). Gene and transfer RNA repertoire from R. vilicus almost perfectly mirror that of other cercozoans, as identified by Wideman et al. (2020). For tRNA, a red-filled square indicates an independent codon reassignment. The ‘extras’ category indicates the unusual codons N (auu) in A (uau) in B (B). Maximum likelihood phylogeny of the COX1 protein (332 amino acid residues) shows clear affiliation of R. vilicus with other cercozoans, congruent with nuclear rRNA phylogeny (C).

Supplementary Figure S2 | R. vilicus and MAP K10-GFP co-localize on the video taken 4 hours post-infection.

Supplementary Figure S3 | R. vilicus and MAP K10-GFP co-localize on the video taken 72 hours post-infection.

Footnotes

1. ^https://github.com/tseemann/barrnap

2. ^https://megasun.bch.umontreal.ca/apps/mfannot/

3. ^https://megasun.bch.umontreal.ca/apps/rnaweasel/

References

Balczun, C., and Scheid, P. L. (2017). Free-living amoebae as hosts for and vectors of intracellular microorganisms with public health significance. Viruses 9:65. doi: 10.3390/v9040065

PubMed Abstract | CrossRef Full Text | Google Scholar

Bass, D., Chao, E. E.-Y., Nikolaev, S., Yabuki, A., Ishida, K.-I., Berney, C., et al. (2009). Phylogeny of novel naked filose and reticulose cercozoa: granofilosea cl. n. and proteomyxidea revised. Protist 160, 75–109. doi: 10.1016/j.protis.2008.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Bass, D., Silberman, J. D., Brown, M. W., Pearce, R. A., Tice, A. K., Jousset, A., et al. (2016). Coprophilic amoebae and flagellates, including Guttulinopsis, Rosculus and Helkesimastix, characterise a divergent and diverse rhizarian radiation and contribute to a large diversity of faecal-associated protists. Environ. Microbiol. 18, 1604–1619. doi: 10.1111/1462-2920.13235

CrossRef Full Text | Google Scholar

Burfeind, O., von Keyserlingk, M. A. G., Weary, D. M., Veira, D. M., and Heuwieser, W. (2010). Short communication: repeatability of measures of rectal temperature in dairy cows. J. Dairy Sci. 93, 624–627. doi: 10.3168/jds.2009-2689

CrossRef Full Text | Google Scholar

Claeys, T. A., and Robinson, R. T. (2018). The many lives of nontuberculous mycobacteria. J. Bacteriol. 200:e00739-17. doi: 10.1128/JB.00739-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Delafont, V., Gasque, M., and Hechard, Y. (2021). Simultaneous hybrid genome sequencing of Vermamoeba vermiformis and its Dependentiae endosymbiont Vermiphilus pyriformis. bioRxiv. doi: 10.1101/2021.04.23.440484

CrossRef Full Text | Google Scholar

Donath, A., Jühling, F., Al-Arab, M., Bernhart, S. H., Reinhardt, F., Stadler, P. F., et al. (2019). Improved annotation of protein-coding genes boundaries in metazoan mitochondrial genomes. Nucleic Acids Res. 47, 10543–10552. doi: 10.1093/nar/gkz833

PubMed Abstract | CrossRef Full Text | Google Scholar

Drancourt, M. (2014). Looking in amoebae as a source of mycobacteria. Microb. Pathog. 77, 119–124. doi: 10.1016/j.micpath.2014.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Escoll, P., Rolando, M., Gomez-Valero, L., and Buchrieser, C. (2013). From amoeba to macrophages: exploring the molecular mechanisms of Legionella pneumophila infection in both hosts. Curr. Top. Microbiol. Immunol. 376, 1–34. doi: 10.1007/82_2013_351

CrossRef Full Text | Google Scholar

Fritz-Laylin, L. K., Prochnik, S. E., Ginger, M. L., Dacks, J. B., Carpenter, M. L., Field, M. C., et al. (2010). The genome of Naegleria gruberi illuminates early eukaryotic versatility. Cells 140, 631–642. doi: 10.1016/j.cell.2010.01.032

CrossRef Full Text | Google Scholar

Hawes, R. S. J. (1955). A Limax-Amœba from the rectum of the grass snake, Natrix natrix, as a facultative aerobe in vitro. Nature 175, 779–780. doi: 10.1038/175779b0

CrossRef Full Text | Google Scholar

Hawes, R. S. (1963). On Rosculus ithacus gen. N., Sp. N. (protozoa, amoebina), with special reference to its mitosis and phylogenetic relations. J. Morphol. 113, 139–149. doi: 10.1002/jmor.1051130202

CrossRef Full Text | Google Scholar

Hess, S., and Suthaus, A. (2022). The vampyrellid amoebae (Vampyrellida, Rhizaria). Protist 173:125854. doi: 10.1016/j.protis.2021.125854

CrossRef Full Text | Google Scholar

Jürgens, K., and Matz, C. (2002). Predation as a shaping force for the phenotypic and genotypic composition of planktonic bacteria. Antonie Van Leeuwenhoek 81, 413–434. doi: 10.1023/A:1020505204959

CrossRef Full Text | Google Scholar

Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., von Haeseler, A., and Jermiin, L. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589. doi: 10.1038/nmeth.4285

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, D., Li, F., Kirton, E. S., Thomas, A., Egan, R. S., An, H., et al. (2019). MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ 7:e7359. doi: 10.7287/peerj.preprints.27522v1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolmogorov, M., Yuan, J., Lin, Y., and Pevzner, P. A. (2019). Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 37, 540–546. doi: 10.1038/s41587-019-0072-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Bannantine, J. P., Zhang, Q., Amonsin, A., May, B. J., Alt, D., et al. (2005). The complete genome sequence of Mycobacterium avium subspecies paratuberculosis. Proc. Natl. Acad. Sci. U. S. A. 102, 12344–12349. doi: 10.1073/pnas.0505662102

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J 17, 10–12. doi: 10.14806/ej.17.1.200

CrossRef Full Text | Google Scholar

Matthey-Doret, C., Colp, M. J., Escoll, P., Thierry, A., Moreau, P., Curtis, B., et al. (2022). Chromosome-scale assemblies of Acanthamoeba castellanii genomes provide insights into Legionella pneumophila infection-related chromatin reorganization. Genome Res. 32, 1698–1710. doi: 10.1101/gr.276375.121

PubMed Abstract | CrossRef Full Text | Google Scholar

Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., von Haeseler, A., et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37, 1530–1534. doi: 10.1093/molbev/msaa015

PubMed Abstract | CrossRef Full Text | Google Scholar

Mura, M., Bull, T. J., Evans, H., Sidi-Boumedine, K., McMinn, L., Rhodes, G., et al. (2006). Replication and long-term persistence of bovine and human strains of Mycobacterium avium subsp. paratuberculosis within Acanthamoeba polyphaga. Appl. Environ. Microbiol. 72, 854–859. doi: 10.1128/AEM.72.1.854-859.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillips, I. L., Everman, J. L., Bermudez, L. E., and Danelishvili, L. (2020). Acanthamoeba castellanii as a screening tool for Mycobacterium avium subspecies paratuberculosis virulence factors with relevance in macrophage infection. Microorganisms 8:1571. doi: 10.3390/microorganisms8101571

CrossRef Full Text | Google Scholar

Pohl, N., Solbach, M. D., and Dumack, K. (2021). The wastewater protist Rhogostoma minus (Thecofilosea, Rhizaria) is abundant, widespread, and hosts Legionellales. Water Res. 203:117566. doi: 10.1016/j.watres.2021.117566

PubMed Abstract | CrossRef Full Text | Google Scholar

Prasad, B. N. K., and Gupta, S. K. (1978). Preliminary report on the engulement and retention of mycobacteria by trophozoites of exenically grown Acanthamoeba castellanii douglas, 1930. Curr. Sci. 47, 245–247.

Google Scholar

Primm, T. P., Lucero, C. A., and Falkinham, J. O. 3rd (2004). Health impacts of environmental mycobacteria. Clin. Microbiol. Rev. 17, 98–106. doi: 10.1128/CMR.17.1.98-106.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodríguez-Zaragoza, S. (1994). Ecology of free-living amoebae. Crit. Rev. Microbiol. 20, 225–241. doi: 10.3109/10408419409114556

CrossRef Full Text | Google Scholar

Salah, I. B., Ghigo, E., and Drancourt, M. (2009). Free-living amoebae, a training field for macrophage resistance of mycobacteria. Clin. Microbiol. Infect. 15, 894–905. doi: 10.1111/j.1469-0691.2009.03011.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Samba-Louaka, A., Delafont, V., Rodier, M.-H., Cateau, E., and Héchard, Y. (2019). Free-living amoebae and squatters in the wild: ecological and molecular features. FEMS Microbiol. Rev. 43, 415–434. doi: 10.1093/femsre/fuz011

CrossRef Full Text | Google Scholar

Samba-Louaka, A., Robino, E., Cochard, T., Branger, M., Delafont, V., Aucher, W., et al. (2018). Environmental Mycobacterium avium subsp. paratuberculosis hosted by free-living amoebae. Front. Cell. Infect. Microbiol. 8:28. doi: 10.3389/fcimb.2018.00028

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, C. A., Rasband, W. S., and Eliceiri, K. W. (2012). NIH image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. doi: 10.1038/nmeth.2089

PubMed Abstract | CrossRef Full Text | Google Scholar

Schuler, G. A., Tice, A. K., Pearce, R. A., Foreman, E., Stone, J., Gammill, S., et al. (2018). Phylogeny and classification of novel diversity in Sainouroidea (Cercozoa, Rhizaria) sheds light on a highly diverse and divergent clade. Protist 169, 853–874. doi: 10.1016/j.protis.2018.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, Y., Queller, D. C., Tian, Y., Zhang, S., Yan, Q., He, Z., et al. (2021). The ecology and evolution of amoeba-bacterium interactions. Appl. Environ. Microbiol. 87:e01866-20. doi: 10.1128/AEM.01866-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Sibbald, S. J., and Archibald, J. M. (2017). More protist genomes needed. Nat Ecol Evol 1:145. doi: 10.1038/s41559-017-0145

PubMed Abstract | CrossRef Full Text | Google Scholar

Solbach, M. D., Bonkowski, M., and Dumack, K. (2021). Novel endosymbionts in rhizarian amoebae imply universal infection of unrelated free-living amoebae by Legionellales. Front. Cell. Infect. Microbiol. 11:642216. doi: 10.3389/fcimb.2021.642216

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, S., Noorian, P., and McDougald, D. (2018). Dual role of mechanisms involved in resistance to predation by protozoa and virulence to humans. Front. Microbiol. 9:1017. doi: 10.3389/fmicb.2018.01017

CrossRef Full Text | Google Scholar

Thomas, V., and McDonnell, G. (2007). Relationship between mycobacteria and amoebae: ecological and epidemiological concerns. Lett. Appl. Microbiol. 45, 349–357. doi: 10.1111/j.1472-765X.2007.02206.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tyml, T., and Dyková, I. (2018). Sappinia sp. (Amoebozoa: Thecamoebida) and Rosculus sp. (SAR: Cercozoa) isolated from king penguin guano collected in the Subantarctic (South Georgia, Salisbury plain) and their coexistence in culture. J. Eukaryot. Microbiol. 65, 544–555. doi: 10.1111/jeu.12500

CrossRef Full Text | Google Scholar

Whan, L., Grant, I. R., and Rowe, M. T. (2006). Interaction between Mycobacterium avium subsp. paratuberculosis and environmental protozoa. BMC Microbiol. 6:63. doi: 10.1186/1471-2180-6-63

PubMed Abstract | CrossRef Full Text | Google Scholar

Wick, R. R., and Holt, K. E. (2022). Polypolish: short-read polishing of long-read bacterial genome assemblies. PLoS Comput. Biol. 18:e1009802. doi: 10.1371/journal.pcbi.1009802

PubMed Abstract | CrossRef Full Text | Google Scholar

Wick, R. R., Judd, L. M., and Holt, K. E. (2023). Assembling the perfect bacterial genome using Oxford Nanopore and Illumina sequencing. PLoS Comput. Biol. 19:e1010905. doi: 10.1371/journal.pcbi.1010905

PubMed Abstract | CrossRef Full Text | Google Scholar

Wideman, J. G., Monier, A., Rodríguez-Martínez, R., Leonard, G., Cook, E., Poirier, C., et al. (2020). Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists. Nat. Microbiol. 5, 154–165. doi: 10.1038/s41564-019-0605-4

CrossRef Full Text | Google Scholar