The genome of Naegleria lovaniensis, the basis for a comparative approach to unravel pathogenicity factors of the human pathogenic amoeba N. fowleri

Naegleria genus pangenome reveals new structural and functional insights into the versatility of these free-living amoebae

Introduction

Free-living amoebae of the Naegleria genus belong to the major protist clade Heterolobosea and are ubiquitously distributed in soil and freshwater habitats. Of the 47 Naegleria species described, N. fowleri is the only one being pathogenic to humans, causing a rare but fulminant primary amoebic meningoencephalitis. Some Naegleria genome sequences are publicly available, but the genetic basis for Naegleria diversity and ability to thrive in diverse environments (including human brain) remains unclear.

Methods

Herein, we constructed a high-quality Naegleria genus pangenome to obtain a comprehensive catalog of genes encoded by these amoebae. For this, we first sequenced, assembled, and annotated six new Naegleria genomes.

Results and Discussion

Genome architecture analyses revealed that Naegleria may use genome plasticity features such as ploidy/aneuploidy to modulate their behavior in different environments. When comparing 14 near-to-complete genome sequences, our results estimated the theoretical Naegleria pangenome as a closed genome, with 13,943 genes, including 3,563 core and 10,380 accessory genes. The functional annotations revealed that a large fraction of Naegleria genes show significant sequence similarity with those already described in other kingdoms, namely Animalia and Plantae. Comparative analyses highlighted a remarkable genomic heterogeneity, even for closely related strains and demonstrate that Naegleria harbors extensive genome variability, reflected in different metabolic repertoires. If Naegleria core genome was enriched in conserved genes essential for metabolic, regulatory and survival processes, the accessory genome revealed the presence of genes involved in stress response, macromolecule modifications, cell signaling and immune response. Commonly reported N. fowleri virulence-associated genes were present in both core and accessory genomes, suggesting that N. fowleri's ability to infect human brain could be related to its unique species-specific genes (mostly of unknown function) and/or to differential gene expression. The construction of Naegleria first pangenome allowed us to move away from a single reference genome (that does not necessarily represent each species as a whole) and to identify essential and dispensable genes in Naegleria evolution, diversity and biology, paving the way for further genomic and post-genomic studies.

Distribution and Current State of Molecular Genetic Characterization in Pathogenic Free-Living Amoebae

Free-living amoebae (FLA) are protozoa widely distributed in the environment, found in a great diversity of terrestrial biomes. Some genera of FLA are linked to human infections. The genus Acanthamoeba is currently classified into 23 genotypes (T1-T23), and of these some (T1, T2, T4, T5, T10, T12, and T18) are known to be capable of causing granulomatous amoebic encephalitis (GAE) mainly in immunocompromised patients while other genotypes (T2, T3, T4, T5, T6, T10, T11, T12, and T15) cause Acanthamoeba keratitis mainly in otherwise healthy patients. Meanwhile, Naegleria fowleri is the causative agent of an acute infection called primary amoebic meningoencephalitis (PAM), while Balamuthia mandrillaris, like some Acanthamoeba genotypes, causes GAE, differing from the latter in the description of numerous cases in patients immunocompetent. Finally, other FLA related to the pathologies mentioned above have been reported; Sappinia sp. is responsible for one case of amoebic encephalitis; Vermamoeba vermiformis has been found in cases of ocular damage, and its extraordinary capacity as endocytobiont for microorganisms of public health importance such as Legionella pneumophila, Bacillus anthracis, and Pseudomonas aeruginosa, among others. This review addressed issues related to epidemiology, updating their geographic distribution and cases reported in recent years for pathogenic FLA.

The transcriptome of Balamuthia mandrillaris trophozoites for structure-guided drug design

Balamuthia mandrillaris, a pathogenic free-living amoeba, causes cutaneous skin lesions as well as granulomatous amoebic encephalitis, a 'brain-eating' disease. As with the other known pathogenic free-living amoebas (Naegleria fowleri and Acanthamoeba species), drug discovery efforts to combat Balamuthia infections of the central nervous system are sparse; few targets have been validated or characterized at the molecular level, and little is known about the biochemical pathways necessary for parasite survival. Current treatments of encephalitis due to B. mandrillaris lack efficacy, leading to case fatality rates above 90%. Using our recently published methodology to discover potential drugs against pathogenic amoebas, we screened a collection of 85 compounds with known antiparasitic activity and identified 59 compounds that impacted the growth of Balamuthia trophozoites at concentrations below 220 µM. Since there is no fully annotated genome or proteome of B. mandrillaris, we sequenced and assembled its transcriptome from a high-throughput RNA-sequencing (RNA-Seq) experiment and located the coding sequences of the genes potentially targeted by the growth inhibitors from our compound screens. We determined the sequence of 17 of these target genes and obtained expression clones for 15 that we validated by direct sequencing. These will be used in the future in combination with the identified hits in structure guided drug discovery campaigns to develop new approaches for the treatment of Balamuthia infections.

Primary Amoebic Meningoencephalitis by Naegleria fowleri: Pathogenesis and Treatments

Naegleria fowleri is a free-living amoeba (FLA) that is commonly known as the "brain-eating amoeba." This parasite can invade the central nervous system (CNS), causing an acute and fulminating infection known as primary amoebic meningoencephalitis (PAM). Even though PAM is characterized by low morbidity, it has shown a mortality rate of 98%, usually causing death in less than two weeks after the initial exposure. This review summarizes the most recent information about N. fowleri, its pathogenic molecular mechanisms, and the neuropathological processes implicated. Additionally, this review includes the main therapeutic strategies described in case reports and preclinical studies, including the possible use of immunomodulatory agents to decrease neurological damage.

Comparative Genomic and Transcriptomic Analysis of Naegleria fowleri Clinical and Environmental Isolates

Out of over 40 species of Naegleria, which are free-living thermophilic amebae found in freshwater and soil worldwide, only Naegleria fowleri infects humans, causing primary amebic meningoencephalitis (PAM), a typically fatal brain disease. To understand the population structure of Naegleria species and the genetic relationships between N. fowleri isolates and to detect pathogenic factors, we characterized 52 novel clinical and environmental N. fowleri genomes and a single Naegleria lovaniensis strain, along with transcriptomic data for a subset of 37 N. fowleri isolates. Whole-genome analysis of 56 isolates from three Naegleria species (N. fowleri, N. lovaniensis, and Naegleria gruberi) identified several genes unique to N. fowleri that have previously been linked to the pathogenicity of N. fowleri, while other unique genes could be associated with novel pathogenicity factors in this highly fatal pathogen. Population structure analysis estimated the presence of 10 populations within the three Naegleria species, of which 7 populations were within N. fowleri. The whole-nuclear-genome (WNG) phylogenetic analysis showed an overall geographical clustering of N. fowleri isolates, with few exceptions, and provided higher resolution in identifying potential clusters of isolates beyond that of the traditional locus typing. There were only 34 genes that showed significant differences in gene expression between the clinical and environmental isolates. Genomic data generated in this study can be used for developing rapid molecular assays and to conduct future population-based global genomic analysis and will also be a valuable addition to genomic reference databases, where shotgun metagenomics data from routine water samples could be searched for the presence of N. fowleri strains. IMPORTANCE N. fowleri, the only known Naegleria species to infect humans, causes fatal brain disease. PAM cases from 1965 to 2016 showed <20 cases per year globally. Out of approximately 150 cases in North America since 1962, only four PAM survivors are known, yielding a >97% case fatality rate, which is critically high. Although the pathogenesis of N. fowleri has been studied for the last 50 years, pathogenetic factors that lead to human infection and breaching the blood-brain barrier remain unknown. In addition, little is known regarding the genomic diversity both within N. fowleri isolates and among Naegleria species. In this study, we generated novel genome sequences and performed comparative genomic and transcriptomic analysis of a set of 52 N. fowleri draft genome sequences from clinical and environmental isolates derived from all over the world in the last 53 years, which will help shape future genome-wide studies and develop sensitive assays for routine surveillance.

Current report on the prevalence of free-living amoebae (FLA) in natural hot springs: a systematic review

The occurrence of potentially pathogenic free-living amoebae (FLA) in natural hot springs is considered a public health concern. FLAs are known to cause serious health outcomes to a wide spectrum of mammalian hosts. The present study aimed to provide the distribution of isolated cases of FLAs in hot springs through a systematic review process of available published articles online. Relevant studies are published between January 2010 and January 2020 involving the isolation of Naegleria spp., Acanthamoeba spp., Balamuthia spp., Sappinia spp., and Vermamoeba spp. in natural hot springs in the United States, South America, North America, Europe, Asia, and Africa. Articles were identified through a search of PubMed and Google Scholar databases. Out of 94 articles screened, a total of 20 articles are included in the study with consideration of established inclusion and exclusion criteria. The most common FLAs isolated in hot springs are Acanthamoeba spp. (134; 48.5%) and Naegleria spp. (127; 46.0%). Other FLAs isolated in hot springs include Balamuthia spp. (2; 0.7%) and Vermamoeba spp. (13; 4.7%). FLA in hot springs used for recreational and medical purposes is a potential source of infection. It is recommended that strict surveillance and maintenance of hot springs be implemented to prevent potential future infection.

Application of the omics sciences to the study of Naegleria fowleri, Acanthamoeba spp., and Balamuthia mandrillaris: current status and future projections

In this review, we focus on the sequenced genomes of the pathogens Naegleria fowleri, Acanthamoeba spp. and Balamuthia mandrillaris, and the remarkable discoveries regarding the pathogenicity and genetic information of these organisms, using techniques related to the various omics branches like genomics, transcriptomics, and proteomics. Currently, novel data produced through comparative genomics analyses and both differential gene and protein expression in these free-living amoebas have allowed for breakthroughs to identify genes unique to N. fowleri, genes with active transcriptional activity, and their differential expression in conditions of modified virulence. Furthermore, orthologous genes of the various nuclear genomes within the Naegleria and Acanthamoeba genera have been clustered. The proteome of B. mandrillaris has been reconstructed through transcriptome data, and its mitochondrial genome structure has been thoroughly described with a unique characteristic that has come to light: a type I intron with the capacity of interrupting genes through its self-splicing ribozymes activity. With the integration of data derived from the diverse omic sciences, there is a potential approximation that reflects the molecular complexity required for the identification of virulence factors, as well as crucial information regarding the comprehension of the molecular mechanisms with which these interact. Altogether, these breakthroughs could contribute to radical advances in both the fields of therapy design and medical diagnosis in the foreseeable future.

The extrachromosomal elements of the Naegleria genus: How little we know

There are currently 47 characterized species in the Naegleria genus of free-living amoebae. Each amoeba has thousands of extrachromosomal elements that are closed circular structures comprised of a single ribosomal DNA (rDNA) copy and a large non-rDNA sequence. Despite the presence of putative open reading frames and introns, ribosomal RNA is the only established transcript. A single origin of DNA replication (ori) has been mapped within the non-rDNA sequence for one species (N. gruberi), a finding that strongly indicates that these episomes replicate independently of the cell's chromosomal DNA component. This article reviews that which has been published about these interesting DNA elements and by analyzing available sequence data, discusses the possibility that different phylogenetically related clusters of Naegleria species individually conserve ori structures and suggests where the rRNA promoter and termination sites may be located.

Complete Genome Sequence of the Naegleria fowleri (Strain LEE) Closed Circular Extrachromosomal Ribosomal DNA Element

The circular extrachromosomal ribosomal DNA (rDNA) element of Naegleria fowleri strain LEE was molecularly cloned and fully sequenced. The element comprises 15,786 bp and encodes a single copy of the organism’s rDNA cistron. The nonribosomal sequence contains five potential open reading frames, two large direct repeat sequences, and numerous smaller repeated-sequence regions.

The circular extrachromosomal ribosomal DNA (rDNA) element of Naegleria fowleri strain LEE was molecularly cloned and fully sequenced. The element comprises 15,786 bp and contains a single copy of the organism’s rDNA cistron. The nonribosomal sequence contains five potential open reading frames, two large direct repeat sequences, and numerous smaller repeated-sequence regions.

Free-living amoebae and other neglected protistan pathogens: Health emergency signals?

Protistan parasites have an undisputed global health impact. However, outside of a few key exceptions, e.g. the agent of malaria, most of these infectious agents are neglected as important health threats. The Symposium entitled "Free-living amoebae and neglected pathogenic protozoa: health emergency signals?" held at the European Congress of Protistology in Rome, July 2019, brought together researchers addressing scientific and clinical questions about some of these fascinating organisms. Topics presented included the molecular basis of pathogenicity in Acanthamoeba; genomics of Naegleria fowleri; and epidemiology of poorly diagnosed enteric protistan species, including Giardia, Cryptosporidium, Blastocystis, Dientamoeba. The Symposium aim was to excite the audience about the opportunities and challenges of research in these underexplored organisms and to underline the public health implications of currently under-appreciated protistan infections. The major take home message is that any knowledge that we gain about these organisms will allow us to better address them, in terms of monitoring and treatment, as sources of future health emergencies.

Proteomics and Lipidomics Investigations to Decipher the Behavior of Willaertia magna C2c Maky According to Different Culture Modes

Willaertia magna C2c Maky is a free-living amoeba that has demonstrated its ability to inhibit the intracellular multiplication of some Legionella pneumophila strains, which are pathogenic bacteria inhabiting the aquatic environment. The Amoeba, an industry involved in the treatment of microbiological risk in the water and plant protection sectors, has developed a natural biocide based on the property of W. magna to manage the proliferation of the pathogen in cooling towers. In axenic liquid medium, amoebas are usually cultivated in adhesion on culture flask. However, we implemented a liquid culture in suspension using bioreactors in order to produce large quantities of W. magna. In order to investigate the culture condition effects on W. magna, we conducted a study based on microscopic, proteomics and lipidomics analyzes. According to the culture condition, amoeba exhibited two different phenotypes. The differential proteomics study showed that amoebas seemed to promote the lipid metabolism pathway in suspension culture, whereas we observed an upregulation of the carbohydrate pathway in adherent culture. Furthermore, we observed an over-regulation of proteins related to the cytoskeleton for W. magna cells grown in adhesion. Regarding the lipid analysis, suspension and adhesion cell growth showed comparable lipid class compositions. However, the differential lipid analysis revealed differences that confirmed cell phenotype differences observed by microscopy and predicted by proteomics. Overall, this study provides us with a better insight into the biology and molecular processes of W. magna in different culture lifestyles.

Ultrastructural, Cytochemical, and Comparative Genomic Evidence of Peroxisomes in Three Genera of Pathogenic Free-Living Amoebae, Including the First Morphological Data for the Presence of This Organelle in Heteroloboseans

Peroxisomes perform various metabolic processes that are primarily related to the elimination of reactive oxygen species and oxidative lipid metabolism. These organelles are present in all major eukaryotic lineages, nevertheless, information regarding the presence of peroxisomes in opportunistic parasitic protozoa is scarce and in many cases it is still unknown whether these organisms have peroxisomes at all. Here, we performed ultrastructural, cytochemical, and bioinformatic studies to investigate the presence of peroxisomes in three genera of free-living amoebae from two different taxonomic groups that are known to cause fatal infections in humans. By transmission electron microscopy, round structures with a granular content limited by a single membrane were observed in Acanthamoeba castellanii, Acanthamoeba griffini, Acanthamoeba polyphaga, Acanthamoeba royreba, Balamuthia mandrillaris (Amoebozoa), and Naegleria fowleri (Heterolobosea). Further confirmation for the presence of peroxisomes was obtained by treating trophozoites in situ with diaminobenzidine and hydrogen peroxide, which showed positive reaction products for the presence of catalase. We then performed comparative genomic analyses to identify predicted peroxin homologues in these organisms. Our results demonstrate that a complete set of peroxins-which are essential for peroxisome biogenesis, proliferation, and protein import-are present in all of these amoebae. Likewise, our in silico analyses allowed us to identify a complete set of peroxins in Naegleria lovaniensis and three novel peroxin homologues in Naegleria gruberi. Thus, our results indicate that peroxisomes are present in these three genera of free-living amoebae and that they have a similar peroxin complement despite belonging to different evolutionary lineages.

Insight into the Lifestyle of Amoeba Willaertia magna during Bioreactor Growth Using Transcriptomics and Proteomics

Willaertia magna C2c maky is a thermophilic free-living amoeba strain that showed ability to eliminate Legionella pneumophila, a pathogenic bacterium living in the aquatic environment. The amoeba industry has proposed the use of Willaertia magna as a natural biocide to control L. pneumophila proliferation in cooling towers. Here, transcriptomic and proteomic studies were carried out in order to expand knowledge on W. magna produced in a bioreactor. Illumina RNA-seq generated 217 million raw reads. A total of 8790 transcripts were identified, of which 6179 and 5341 were assigned a function through comparisons with National Center of Biotechnology Information (NCBI) reference sequence and the Clusters of Orthologous Groups of proteins (COG) databases, respectively. To corroborate these transcriptomic data, we analyzed the W. magna proteome using LC–MS/MS. A total of 3561 proteins were identified. The results of transcriptome and proteome analyses were highly congruent. Metabolism study showed that W. magna preferentially consumed carbohydrates and fatty acids to grow. Finally, an in-depth analysis has shown that W. magna produces several enzymes that are probably involved in the metabolism of secondary metabolites. Overall, our multi-omic study of W. magna opens the way to a better understanding of the genetics and biology of this amoeba.

Identification of differential protein recognition pattern between Naegleria fowleri and Naegleria lovaniensis

Many pathogenicity factors are involved in the development of primary amoebic meningoencephalitis (PAM) caused by N fowleri. However, most of them are not exclusive for N fowleri and they have not even been described in other nonpathogenic Naegleria species. Therefore, the objective of this work was to identify differential proteins and protein pattern recognition between Naegleria fowleri and Naegleria lovaniensis using antibodies anti-N fowleri as strategy to find vaccine candidates against meningoencephalitis. Electrophoresis and Western blots conventional and 2-DE were performed for the identification of antigenic proteins, and these were analysed by the mass spectrometry technique. The results obtained in 2-DE gels and Western blot showed very notable differences in spot intensity between these two species, specifically those with relative molecular weight of 100, 75, 50 and 19 kDa. Some spots corresponding to these molecular weights were identified as actin fragment, myosin II, heat shock protein, membrane protein Mp2CL5 among others, with differences in theoretical post-translational modifications. In this work, we found differences in antigenic proteins between both species, proteins that could be used for a further development of vaccines against N fowleri infection.

Description of Virulent Factors and Horizontal Gene Transfers of Keratitis-Associated Amoeba Acanthamoeba Triangularis by Genome Analysis

Acanthamoeba triangularis strain SH 621 is a free-living amoeba belonging to Acanthamoeba ribo-genotype T4. This ubiquitous protist is among the free-living amoebas responsible for Acanthamoeba keratitis, a severe infection of human cornea. Genome sequencing and genomic comparison were carried out to explore the biological functions and to better understand the virulence mechanism related to the pathogenicity of Acanthamoeba keratitis. The genome assembly harbored a length of 66.43 Mb encompassing 13,849 scaffolds. The analysis of predicted proteins reported the presence of 37,062 ORFs. A complete annotation revealed 33,168 and 16,605 genes that matched with NCBI non-redundant protein sequence (nr) and Cluster of Orthologous Group of proteins (COG) databases, respectively. The Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) annotation reported a great number of genes related to carbohydrate, amino acid and lipid metabolic pathways. The pangenome performed with 8 available amoeba genomes belonging to genus Acanthamoeba revealed a core genome containing 843 clusters of orthologous genes with a ratio core genome/pangenome of less than 0.02. We detected 48 genes related to virulent factors of Acanthamoeba keratitis. Best hit analyses in nr database identified 99 homologous genes shared with amoeba-resisting microorganisms. This study allows the deciphering the genome of a free-living amoeba with medical interest and provides genomic data to better understand virulence-related Acanthamoeba keratitis.

Investigation of potential pathogenicity of Willaertia magna by investigating the transfer of bacteria pathogenicity genes into its genome

Willaertia magna c2c maky is a thermophilic amoeba closely related to the genus Naegleria. This free-living amoeba has the ability to eliminate Legionella pneumophila, which is an amoeba-resisting bacterium living in an aquatic environment. To prevent the proliferation of L. pneumophila in cooling towers, the use of W. magna as natural biocide has been proposed. To provide a better understanding of the W. magna genome, whole-genome sequencing was performed through the study of virulence factors and lateral gene transfers. This amoeba harbors a genome of 36.5 megabases with 18,519 predicted genes. BLASTp analyses reported protein homology between 136 W. magna sequences and amoeba-resistant microorganisms. Horizontal gene transfers were observed based on the basis of the phylogenetic reconstruction hypothesis. We detected 15 homologs of N. fowleri genes related to virulence, although these latter were also found in the genome of N. gruberi, which is a non-pathogenic amoeba. Furthermore, the cytotoxicity test performed on human cells supports the hypothesis that the strain c2c maky is a non-pathogenic amoeba. This work explores the genomic repertory for the first draft genome of genus Willaertia and provides genomic data for further comparative studies on virulence of related pathogenic amoeba, N. fowleri.

Nanopore sequencing improves the draft genome of the human pathogenic amoeba Naegleria fowleri

Naegleria fowleri is an environmental protist found in soil and warm freshwater sources worldwide and is known for its ability to infect humans and causing a rapid and mostly fatal primary amoebic meningoencephalitis. When contaminated water enters the nose, the facultative parasite follows the olfactory nerve and enters the brain by crossing the cribriform plate where it causes tissue damage and haemorrhagic necrosis. Although N. fowleri has been studied for several years, the mechanisms of pathogenicity are still poorly understood. Furthermore, there is a lack of knowledge on the genomic level and the current reference assembly is limited in contiguity. To improve the draft genome and to investigate pathogenicity factors, we sequenced the genome of N. fowleri using Oxford Nanopore Technology (ONT). Assembly and polishing of the long reads resulted in a high-quality draft genome whose N50 is 18 times higher than the previously published genome. The prediction of potentially secreted proteins revealed a large proportion of enzymes with a hydrolysing function, which could play an important role during the pathogenesis and account for the destructive nature of primary amoebic meningoencephalitis. The improved genome provides the basis for further investigation unravelling the biology and the pathogenic potential of N. fowleri.

Free-living amoebae and squatters in the wild: ecological and molecular features

Free-living amoebae are protists frequently found in water and soils. They feed on other microorganisms, mainly bacteria, and digest them through phagocytosis. It is accepted that these amoebae play an important role in the microbial ecology of these environments. There is a renewed interest for the free-living amoebae since the discovery of pathogenic bacteria that can resist phagocytosis and of giant viruses, underlying that amoebae might play a role in the evolution of other microorganisms, including several human pathogens. Recent advances, using molecular methods, allow to bring together new information about free-living amoebae. This review aims to provide a comprehensive overview of the newly gathered insights into (1) the free-living amoeba diversity, assessed with molecular tools, (2) the gene functions described to decipher the biology of the amoebae and (3) their interactions with other microorganisms in the environment.