The mitochondrial genome and a 60-kb nuclear DNA segment from Naegleria fowleri, the causative agent of primary amoebic meningoencephalitis

Modes and mechanisms for the inheritance of mitochondria and plastids in pathogenic protists

Pathogenic protists are responsible for many diseases that significantly impact human and animal health across the globe. Almost all protists possess mitochondria or mitochondrion-related organelles, and many contain plastids. These endosymbiotic organelles are crucial to survival and provide well-validated and widely utilised drug targets in parasitic protists such as Plasmodium and Toxoplasma. However, mutations within the organellar genomes of mitochondria and plastids can lead to drug resistance. Such mutations ultimately challenge our ability to control and eradicate the diseases caused by these pathogenic protists. Therefore, it is important to understand how organellar genomes, and the resistance mutations encoded within them, are inherited during protist sexual reproduction and how this may impact the spread of drug resistance and future therapeutic approaches to target these organelles. In this review, we detail what is known about mitochondrial and plastid inheritance during sexual reproduction across different pathogenic protists, often turning to their better studied, nonpathogenic relatives for insight.

Immunogens in Balamuthia mandrillaris: a proteomic exploration

Balamuthia mandrillaris is the causative agent of granulomatous amoebic encephalitis, a rare and often fatal infection affecting the central nervous system. The amoeba is isolated from diverse environmental sources and can cause severe infections in both immunocompromised and immunocompetent individuals. Given the limited understanding of B. mandrillaris, our research aimed to explore its protein profile, identifying potential immunogens crucial for early granulomatous amoebic encephalitis diagnosis. Cultures of B. mandrillaris and other amoebas were grown under axenic conditions, and total amoebic extracts were obtained. Proteomic analyses, including two-dimensional electrophoresis and mass spectrometry, were performed. A 50-kDa band showed a robust recognition of antibodies from immunized BALB/c mice; peptides contained in this band were matched with elongation factor-1 alpha, which emerged as a putative key immunogen. Besides, lectin blotting revealed the presence of glycoproteins in B. mandrillaris, and confocal microscopy demonstrated the focal distribution of the 50-kDa band throughout trophozoites. Cumulatively, these observations suggest the participation of the 50-kDa band in adhesion and recognition mechanisms. Thus, these collective findings demonstrate some protein characteristics of B. mandrillaris, opening avenues for understanding its pathogenicity and developing diagnostic and therapeutic strategies.

Meteora sporadica, a protist with incredible cell architecture, is related to Hemimastigophora

"Kingdom-level" branches are being added to the tree of eukaryotes at a rate approaching one per year, with no signs of slowing down.1,2,3,4 Some are completely new discoveries, whereas others are morphologically unusual protists that were previously described but lacked molecular data. For example, Hemimastigophora are predatory protists with two rows of flagella that were known since the 19th century but proved to represent a new deep-branching eukaryote lineage when phylogenomic analyses were conducted.2Meteora sporadica5 is a protist with a unique morphology; cells glide over substrates along a long axis of anterior and posterior projections while a pair of lateral "arms" swing back and forth, a motility system without any obvious parallels. Originally, Meteora was described by light microscopy only, from a short-term enrichment of deep-sea sediment. A small subunit ribosomal RNA (SSU rRNA) sequence was reported recently, but the phylogenetic placement of Meteora remained unresolved.6 Here, we investigated two cultivated Meteora sporadica isolates in detail. Transmission electron microscopy showed that both the anterior-posterior projections and the arms are supported by microtubules originating from a cluster of subnuclear microtubule organizing centers (MTOCs). Neither have a flagellar axoneme-like structure. Sequencing the mitochondrial genome showed this to be among the most gene-rich known, outside jakobids. Remarkably, phylogenomic analyses of 254 nuclear protein-coding genes robustly support a close relationship with Hemimastigophora. Our study suggests that Meteora and Hemimastigophora together represent a morphologically diverse "supergroup" and thus are important for resolving the tree of eukaryote life and early eukaryote evolution.

Genomics and transcriptomics yields a system-level view of the biology of the pathogen Naegleria fowleri

BACKGROUND

The opportunistic pathogen Naegleria fowleri establishes infection in the human brain, killing almost invariably within 2 weeks. The amoeba performs piece-meal ingestion, or trogocytosis, of brain material causing direct tissue damage and massive inflammation. The cellular basis distinguishing N. fowleri from other Naegleria species, which are all non-pathogenic, is not known. Yet, with the geographic range of N. fowleri advancing, potentially due to climate change, understanding how this pathogen invades and kills is both important and timely.

RESULTS

Here, we report an -omics approach to understanding N. fowleri biology and infection at the system level. We sequenced two new strains of N. fowleri and performed a transcriptomic analysis of low- versus high-pathogenicity N. fowleri cultured in a mouse infection model. Comparative analysis provides an in-depth assessment of encoded protein complement between strains, finding high conservation. Molecular evolutionary analyses of multiple diverse cellular systems demonstrate that the N. fowleri genome encodes a similarly complete cellular repertoire to that found in free-living N. gruberi. From transcriptomics, neither stress responses nor traits conferred from lateral gene transfer are suggested as critical for pathogenicity. By contrast, cellular systems such as proteases, lysosomal machinery, and motility, together with metabolic reprogramming and novel N. fowleri proteins, are all implicated in facilitating pathogenicity within the host. Upregulation in mouse-passaged N. fowleri of genes associated with glutamate metabolism and ammonia transport suggests adaptation to available carbon sources in the central nervous system.

CONCLUSIONS

In-depth analysis of Naegleria genomes and transcriptomes provides a model of cellular systems involved in opportunistic pathogenicity, uncovering new angles to understanding the biology of a rare but highly fatal pathogen.

Single-cell genomics unveils a canonical origin of the diverse mitochondrial genomes of euglenozoans

Background

The supergroup Euglenozoa unites heterotrophic flagellates from three major clades, kinetoplastids, diplonemids, and euglenids, each of which exhibits extremely divergent mitochondrial characteristics. Mitochondrial genomes (mtDNAs) of euglenids comprise multiple linear chromosomes carrying single genes, whereas mitochondrial chromosomes are circular non-catenated in diplonemids, but circular and catenated in kinetoplastids. In diplonemids and kinetoplastids, mitochondrial mRNAs require extensive and diverse editing and/or trans-splicing to produce mature transcripts. All known euglenozoan mtDNAs exhibit extremely short mitochondrial small (rns) and large (rnl) subunit rRNA genes, and absence of tRNA genes. How these features evolved from an ancestral bacteria-like circular mitochondrial genome remains unanswered.

Results

We sequenced and assembled 20 euglenozoan single-cell amplified genomes (SAGs). In our phylogenetic and phylogenomic analyses, three SAGs were placed within kinetoplastids, 14 within diplonemids, one (EU2) within euglenids, and two SAGs with nearly identical small subunit rRNA gene (18S) sequences (EU17/18) branched as either a basal lineage of euglenids, or as a sister to all euglenozoans. Near-complete mitochondrial genomes were identified in EU2 and EU17/18. Surprisingly, both EU2 and EU17/18 mitochondrial contigs contained multiple genes and one tRNA gene. Furthermore, EU17/18 mtDNA possessed several features unique among euglenozoans including full-length rns and rnl genes, six mitoribosomal genes, and nad11, all likely on a single chromosome.

Conclusions

Our data strongly suggest that EU17/18 is an early-branching euglenozoan with numerous ancestral mitochondrial features. Collectively these data contribute to untangling the early evolution of euglenozoan mitochondria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01035-y.

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.

Evolutionary History of Mitochondrial Genomes in Discoba, Including the Extreme Halophile Pleurostomum flabellatum (Heterolobosea)

Data from Discoba (Heterolobosea, Euglenozoa, Tsukubamonadida, and Jakobida) are essential to understand the evolution of mitochondrial genomes (mitogenomes), because this clade includes the most primitive-looking mitogenomes known, as well some extremely divergent genome information systems. Heterolobosea encompasses more than 150 described species, many of them from extreme habitats, but only six heterolobosean mitogenomes have been fully sequenced to date. Here we complete the mitogenome of the heterolobosean Pleurostomum flabellatum, which is extremely halophilic and reportedly also lacks classical mitochondrial cristae, hinting at reduction or loss of respiratory function. The mitogenome of P. flabellatum maps as a 57,829-bp-long circular molecule, including 40 coding sequences (19 tRNA, two rRNA, and 19 orfs). The gene content and gene arrangement are similar to Naegleria gruberi and Naegleria fowleri, the closest relatives with sequenced mitogenomes. The P. flabellatum mitogenome contains genes that encode components of the electron transport chain similar to those of Naegleria mitogenomes. Homology searches against a draft nuclear genome showed that P. flabellatum has two homologs of the highly conserved Mic60 subunit of the MICOS complex, and likely lost Mic19 and Mic10. However, electron microscopy showed no cristae structures. We infer that P. flabellatum, which originates from high salinity (313‰) water where the dissolved oxygen concentration is low, possesses a mitochondrion capable of aerobic respiration, but with reduced development of cristae structure reflecting limited use of this aerobic capacity (e.g., microaerophily).

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.

Recent advances in understanding mitochondrial genome diversity

Ever since its discovery, the double-stranded DNA contained in the mitochondria of eukaryotes has fascinated researchers because of its bacterial endosymbiotic origin, crucial role in encoding subunits of the respiratory complexes, compact nature, and specific inheritance mechanisms. In the last few years, high-throughput sequencing techniques have accelerated the sequencing of mitochondrial genomes (mitogenomes) and uncovered the great diversity of organizations, gene contents, and modes of replication and transcription found in living eukaryotes. Some early divergent lineages of unicellular eukaryotes retain certain synteny and gene content resembling those observed in the genomes of alphaproteobacteria (the inferred closest living group of mitochondria), whereas others adapted to anaerobic environments have drastically reduced or even lost the mitogenome. In the three main multicellular lineages of eukaryotes, mitogenomes have pursued diverse evolutionary trajectories in which different types of molecules (circular versus linear and single versus multipartite), gene structures (with or without self-splicing introns), gene contents, gene orders, genetic codes, and transfer RNA editing mechanisms have been selected. Whereas animals have evolved a rather compact mitochondrial genome between 11 and 50 Kb in length with a highly conserved gene content in bilaterians, plants exhibit large mitochondrial genomes of 66 Kb to 11.3 Mb with large intergenic repetitions prone to recombination, and fungal mitogenomes have intermediate sizes of 12 to 236 Kb.

trans-Cinnamic Acid Conjugated Gold Nanoparticles as Potent Therapeutics against Brain-Eating Amoeba Naegleria fowleri

Primary amoebic meningoencephalitis (PAM), a deadly brain infection, is caused by brain-eating amoeba Naegleria fowleri. The current first line of treatment against PAM is a mixture of amphotericin B, rifampin, and miltefosine. Since, no single effective drug has been developed so far, the mortality rate is above 95%. Moreover, severe adverse side effects are associated with these drugs. Nanotechnology has provided several advances in biomedical applications especially in drug delivery and diagnosis. Herein, for the first time we report antiamoebic properties of cinnamic acid (CA) and gold nanoparticles conjugated with CA (CA-AuNPs) against N. fowleri. CA-AuNPs were successfully synthesized by sodium borohydride reduction of tetrachloroauric acid. Size and morphology were determined by atomic force microscopy (AFM) while the surface plasmon resonance band was analyzed by ultraviolet-visible (UV-vis) spectrophotometry for the characterization of the nanoparticles. Amoebicidal and cytopathogenicity (host cell cytotoxicity) assays revealed that both CA and CA-AuNPs displayed significant anti- N. fowleri properties ( P < 0.05), whereas nanoparticles conjugation further enhanced the anti- N. fowleri effects of CA. This study established a potential drug lead, while CA-AuNPs appear to be promising candidate for drug discovery against PAM.

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.

Horizontally-acquired genetic elements in the mitochondrial genome of a centrohelid Marophrys sp. SRT127

Mitochondrial genomes exhibit diverse features among eukaryotes in the aspect of gene content, genome structure, and the mobile genetic elements such as introns and plasmids. Although the number of published mitochondrial genomes is increasing at tremendous speed, those of several lineages remain unexplored. Here, we sequenced the complete mitochondrial genome of a unicellular heterotrophic eukaryote, Marophrys sp. SRT127 belonging to the Centroheliozoa, as the first report on this lineage. The circular-mapped mitochondrial genome, which is 113,062 bp in length, encodes 69 genes typically found in mitochondrial genomes. In addition, the Marophrys mitochondrial genome contains 19 group I introns. Of these, 11 introns have genes for homing endonuclease (HE) and phylogenetic analyses of HEs have shown that at least five Marophrys HEs are related to those in green algal plastid genomes, suggesting intron transfer between the Marophrys mitochondrion and green algal plastids. We also discovered a putative mitochondrial plasmid in linear form. Two genes encoded in the circular-mapped mitochondrial genome were found to share significant similarities to those in the linear plasmid, suggesting that the plasmid was integrated into the mitochondrial genome. These findings expand our knowledge on the diversity and evolution of the mobile genetic elements in mitochondrial genomes.

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

Background

Members of the genus Naegleria are free-living eukaryotes with the capability to transform from the amoeboid form into resting cysts or moving flagellates in response to environmental conditions. More than 40 species have been characterized, but only Naegleria fowleri (N. fowleri) is known as a human pathogen causing primary amoebic meningoencephalitis (PAM), a fast progressing and mostly fatal disease of the central nervous system. Several studies report an involvement of phospholipases and other molecular factors, but the mechanisms involved in pathogenesis are still poorly understood. To gain a better understanding of the relationships within the genus of Naegleria and to investigate pathogenicity factors of N. fowleri, we characterized the genome of its closest non-pathogenic relative N. lovaniensis.

Results

To gain insights into the taxonomy of Naegleria, we sequenced the genome of N. lovaniensis using long read sequencing technology. The assembly of the data resulted in a 30 Mb genome including the circular mitochondrial sequence. Unravelling the phylogenetic relationship using OrthoMCL protein clustering and maximum likelihood methods confirms the close relationship of N. lovaniensis and N. fowleri. To achieve an overview of the diversity of Naegleria proteins and to assess characteristics of the human pathogen N. fowleri, OrthoMCL protein clustering including data of N. fowleri, N. lovaniensis and N. gruberi was performed. GO enrichment analysis shows an association of N. fowleri specific proteins to the GO terms “Membrane” and “Protein Secretion.”

Conclusion

In this study, we characterize the hitherto unknown genome of N. lovaniensis. With the description of the 30 Mb genome, a further piece is added to reveal the complex taxonomic relationship of Naegleria. Further, the whole genome sequencing data confirms the hypothesis of the close relationship between N. fowleri and N. lovaniensis. Therefore, the genome of N. lovaniensis provides the basis for further comparative approaches on the molecular and genomic level to unravel pathogenicity factors of its closest human pathogenic relative N. fowleri and possible treatment options for the rare but mostly fatal primary meningoencephalitis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4994-1) contains supplementary material, which is available to authorized users.

Heterogeneous Antimicrobial Susceptibility Characteristics in Pseudomonas aeruginosa Isolates from Cystic Fibrosis Patients

Patients with cystic fibrosis endure “chronic focal infections” with a variety of microorganisms. One microorganism, Pseudomonas aeruginosa, adapts to the host and develops resistance to a wide range of antimicrobials. Interestingly, as the infection progresses, multiple isogenic strains of P. aeruginosa emerge and coexist within the airways of these patients. Despite a common parental origin, the multiple strains of P. aeruginosa develop vastly different susceptibility patterns to actively used antimicrobial agents—a phenomenon we define as “heterogeneous MICs.” By sequencing pairs of P. aeruginosa isolates displaying heterogeneous MICs, we observed widespread isogenic gene lesions in drug transporters, DNA mismatch repair mechanisms, and many other structural or cellular functions. Coupled with the heterogeneous MICs, these genetic lesions demonstrated a symbiotic response to host selection and suggested evolution of a multicellular syntrophic bacterial lifestyle. Current laboratory standard interpretive criteria do not address the emergence of heterogeneous growth and susceptibilities in vitro with treatment implications.

Clinical isolates of Pseudomonas aeruginosa from patients with cystic fibrosis (CF) are known to differ from those associated with non-CF hosts by colony morphology, drug susceptibility patterns, and genomic hypermutability. Pseudomonas aeruginosa isolates from CF patients have long been recognized for their overall reduced rate of antimicrobial susceptibility, but their intraclonal MIC heterogeneity has long been overlooked. Using two distinct cohorts of clinical strains (n = 224 from 56 CF patients, n = 130 from 68 non-CF patients) isolated in 2013, we demonstrated profound Etest MIC heterogeneity in CF P. aeruginosa isolates in comparison to non-CF P. aeruginosa isolates. On the basis of whole-genome sequencing of 19 CF P. aeruginosa isolates from 9 patients with heterogeneous MICs, the core genome phylogenetic tree confirmed the within-patient CF P. aeruginosa clonal lineage along with considerable coding sequence variability. No extrachromosomal DNA elements or previously characterized antibiotic resistance mutations could account for the wide divergence in antimicrobial MICs between P. aeruginosa coisolates, though many heterogeneous mutations in efflux and porin genes and their regulators were present. A unique OprD sequence was conserved among the majority of isolates of CF P. aeruginosa analyzed, suggesting a pseudomonal response to selective pressure that is common to the isolates. Genomic sequence data also suggested that CF pseudomonal hypermutability was not entirely due to mutations in mutL, mutS, and uvr. We conclude that the net effect of hundreds of adaptive mutations, both shared between clonally related isolate pairs and unshared, accounts for their highly heterogeneous MIC variances. We hypothesize that this heterogeneity is indicative of the pseudomonal syntrophic-like lifestyle under conditions of being “locked” inside a host focal airway environment for prolonged periods.

IMPORTANCE Patients with cystic fibrosis endure “chronic focal infections” with a variety of microorganisms. One microorganism, Pseudomonas aeruginosa, adapts to the host and develops resistance to a wide range of antimicrobials. Interestingly, as the infection progresses, multiple isogenic strains of P. aeruginosa emerge and coexist within the airways of these patients. Despite a common parental origin, the multiple strains of P. aeruginosa develop vastly different susceptibility patterns to actively used antimicrobial agents—a phenomenon we define as “heterogeneous MICs.” By sequencing pairs of P. aeruginosa isolates displaying heterogeneous MICs, we observed widespread isogenic gene lesions in drug transporters, DNA mismatch repair machinery, and many other structural or cellular functions. Coupled with the heterogeneous MICs, these genetic lesions demonstrated a symbiotic response to host selection and suggested evolution of a multicellular syntrophic bacterial lifestyle. Current laboratory standard interpretive criteria do not address the emergence of heterogeneous growth and susceptibilities in vitro with treatment implications.

Mitochondrial Genome Evolution and a Novel RNA Editing System in Deep-Branching Heteroloboseids

Discoba (Excavata) is an evolutionarily important group of eukaryotes that includes Jakobida, with the most bacterial-like mitochondrial genomes known, and Euglenozoa, many of which have extensively fragmented mitochondrial genomes. However, little is known about the mitochondrial genomes of Heterolobosea, the third main group of Discoba. Here, we studied two heteroloboseids—an undescribed amoeba “BB2” and Pharyngomonas kirbyi. Phylogenomic analysis revealed that they form a clade that is a sister group to all other Heterolobosea. We characterized the mitochondrial genomes of BB2 and P. kirbyi, which encoded 44 and 48 putative protein-coding genes respectively. Their gene contents were similar to that of Naegleria. In BB2, mitochondrially encoded RNAs were heavily edited, with ∼500 mononucleotide insertion events, mostly guanosines. These insertions always have the same identity as an adjacent nucleotide. Editing occurs in all ribosomal RNAs and protein-coding transcripts except one, and half of the transfer RNAs. Analysis of Illumina deep-sequencing data suggested that this RNA editing is very accurate and efficient, and most likely co-transcriptional. The dissimilarity of this editing process to other RNA editing phenomena in discobids, as well as its apparent absence in P. kirbyi, suggest that this remarkably extensive system of insertional editing evolved independently in the BB2 lineage, after its divergence from the P. kirbyi lineage.

Biology and pathogenesis of Naegleria fowleri

Naegleria fowleri is a protist pathogen that can cause lethal brain infection. Despite decades of research, the mortality rate related with primary amoebic meningoencephalitis owing to N. fowleri remains more than 90%. The amoebae pass through the nose to enter the central nervous system killing the host within days, making it one of the deadliest opportunistic parasites. Accordingly, we present an up to date review of the biology and pathogenesis of N. fowleri and discuss needs for future research against this fatal infection.

Aerobic mitochondria of parasitic protists: Diverse genomes and complex functions

In this review the main features of the mitochondria of aerobic parasitic protists are discussed. While the best characterized organelles are by far those of kinetoplastid flagellates and Plasmodium, we also consider amoebae Naegleria and Acanthamoeba, a ciliate Ichthyophthirius and related lineages. The simplistic view of the mitochondrion as just a power house of the cell has already been abandoned in multicellular organisms and available data indicate that this also does not apply for protists. We discuss in more details the following mitochondrial features: genomes, post-transcriptional processing, translation, biogenesis of iron-sulfur complexes, heme metabolism and the electron transport chain. Substantial differences in all these core mitochondrial features between lineages are compatible with the view that aerobic protists harbor organelles that are more complex and flexible than previously appreciated.

Group II Intron-Mediated Trans-Splicing in the Gene-Rich Mitochondrial Genome of an Enigmatic Eukaryote, Diphylleia rotans

Although mitochondria have evolved from a single endosymbiotic event, present day mitochondria of diverse eukaryotes display a great range of genome structures, content and features. Group I and group II introns are two features that are distributed broadly but patchily in mitochondrial genomes across branches of the tree of eukaryotes. While group I intron-mediated trans-splicing has been reported from some lineages distantly related to each other, findings of group II intron-mediated trans-splicing has been restricted to members of the Chloroplastida. In this study, we found the mitochondrial genome of the unicellular eukaryote Diphylleia rotans possesses currently the second largest gene repertoire. On the basis of a probable phylogenetic position of Diphylleia, which is located within Amorphea, current mosaic gene distribution in Amorphea must invoke parallel gene losses from mitochondrial genomes during evolution. Most notably, although the cytochrome c oxidase subunit (cox) 1 gene was split into four pieces which located at a distance to each other, we confirmed that a single mature mRNA that covered the entire coding region could be generated by group II intron-mediated trans-splicing. This is the first example of group II intron-mediated trans-splicing outside Chloroplastida. Similar trans-splicing mechanisms likely work for bipartitely split cox2 and nad3 genes to generate single mature mRNAs. We finally discuss origin and evolution of this type of trans-splicing in D. rotans as well as in eukaryotes.

Computational identification of putative miRNAs and their target genes in pathogenic amoeba Naegleria fowleri

Naegleria fowleri is a parasitic unicellular free living eukaryotic amoeba. The parasite spreads through contaminated water and causes primary amoebic meningoencephalitis (PAM). Therefore, it is of interest to understand its molecular pathogenesis. Hence, we analyzed the parasite genome for miRNAs (microRNAs) that are non-coding, single stranded RNA molecules. We identified 245 miRNAs using computational methods in N. fowleri, of which five miRNAs are conserved. The predicted miRNA targets were analyzed by using miRanda (software) and further studied the functions by subsequently annotating using AmiGo (a gene ontology web tool).

Clinical metagenomic identification of Balamuthia mandrillaris encephalitis and assembly of the draft genome: the continuing case for reference genome sequencing

Background

Primary amoebic meningoencephalitis (PAM) is a rare, often lethal, cause of encephalitis, for which early diagnosis and prompt initiation of combination antimicrobials may improve clinical outcomes.

Methods

In this study, we sequenced a full draft assembly of the Balamuthia mandrillaris genome (44.2 Mb in size) from a rare survivor of PAM, and recovered the mitochondrial genome from six additional Balamuthia strains. We also used unbiased metagenomic next-generation sequencing (NGS) and SURPI bioinformatics analysis to diagnose an ultimately fatal case of Balamuthia mandrillaris encephalitis in a 15-year-old girl.

Results and Discussion

Comparative analysis of the mitochondrial genome and high-copy number genes from six additional Balamuthia mandrillaris strains demonstrated remarkable sequence variation, and the closest Balamuthia homologs corresponded to other amoebae, hydroids, algae, slime molds, and peat moss. Real-time NGS testing of hospital day 6 CSF and brain biopsy samples detected Balamuthia on the basis of high-quality hits to 16S and 18S ribosomal RNA sequences present in the National Center for Biotechnology Information (NCBI) nt reference database. The presumptive diagnosis of PAM by visualization of amoebae on brain biopsy histopathology and NGS analysis was subsequently confirmed at the US Centers for Disease Control and Prevention (CDC) using a Balamuthia-specific PCR assay. Retrospective analysis of a day 1 CSF sample revealed that more timely identification of Balamuthia by metagenomic NGS, potentially resulting in a better clinical outcome, would have required availability of the complete genome sequence.

Conclusions

These results underscore the diverse evolutionary origins of Balamuthia mandrillaris, provide new targets for diagnostic assay development, and will facilitate further investigations of the biology and pathogenesis of this eukaryotic pathogen. The failure to identify PAM from a day 1 sample without a fully sequenced Balamuthia genome in the database highlights the critical importance of whole-genome reference sequences for microbial detection by metagenomic NGS.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-015-0235-2) contains supplementary material, which is available to authorized users.

An ensemble strategy that significantly improves de novo assembly of microbial genomes from metagenomic next-generation sequencing data

Next-generation sequencing (NGS) approaches rapidly produce millions to billions of short reads, which allow pathogen detection and discovery in human clinical, animal and environmental samples. A major limitation of sequence homology-based identification for highly divergent microorganisms is the short length of reads generated by most highly parallel sequencing technologies. Short reads require a high level of sequence similarities to annotated genes to confidently predict gene function or homology. Such recognition of highly divergent homologues can be improved by reference-free (de novo) assembly of short overlapping sequence reads into larger contigs. We describe an ensemble strategy that integrates the sequential use of various de Bruijn graph and overlap-layout-consensus assemblers with a novel partitioned sub-assembly approach. We also proposed new quality metrics that are suitable for evaluating metagenome de novo assembly. We demonstrate that this new ensemble strategy tested using in silico spike-in, clinical and environmental NGS datasets achieved significantly better contigs than current approaches.

Gene content evolution in Discobid mitochondria deduced from the phylogenetic position and complete mitochondrial genome of Tsukubamonas globosa

The unicellular eukaryotic assemblage Discoba (Excavata) comprises four lineages: the Heterolobosea, Euglenozoa, Jakobida, and Tsukubamonadida. Discoba has been considered as a key assemblage for understanding the early evolution of mitochondrial (mt) genomes, as jakobids retain the most gene-rich (i.e., primitive) genomes compared with any other eukaryotes determined to date. However, to date, mt genome sequences have been completed for only a few groups within Discoba, including jakobids, two closely related heteroloboseans, and kinetoplastid euglenozoans. The Tsukubamonadida is the least studied lineage, as the order was only recently established with the description of a sole representative species, Tsukubamonas globosa. The evolutionary relationship between T. globosa and other discobids has yet to be resolved, and no mt genome data are available for this particular organism. Here, we use a “phylogenomic” approach to resolve the relationship between T. globosa, heteroloboseans, euglenozoans, and jakobids. In addition, we have characterized the mt genome of T. globosa (48,463 bp in length), which encodes 52 putative protein-coding and 29 RNA genes. By mapping the gene repertoires of discobid mt genomes onto the well-resolved Discoba tree, we model gene loss events during the evolution of discobid mt genomes.

Detection of biomarkers of pathogenic Naegleria fowleri through mass spectrometry and proteomics

Emerging methods based on mass spectrometry (MS) can be used in the rapid identification of microorganisms. Thus far, these practical and rapidly evolving methods have mainly been applied to characterize prokaryotes. We applied matrix-assisted laser-desorption-ionization-time-of-flight mass spectrometry MALDI-TOF MS in the analysis of whole cells of 18 N. fowleri isolates belonging to three genotypes. Fourteen originated from the cerebrospinal fluid or brain tissue of primary amoebic meningoencephalitis patients and four originated from water samples of hot springs, rivers, lakes or municipal water supplies. Whole Naegleria trophozoites grown in axenic cultures were washed and mixed with MALDI matrix. Mass spectra were acquired with a 4700 TOF-TOF instrument. MALDI-TOF MS yielded consistent patterns for all isolates examined. Using a combination of novel data processing methods for visual peak comparison, statistical analysis and proteomics database searching we were able to detect several biomarkers that can differentiate all species and isolates studied, along with common biomarkers for all N. fowleri isolates. Naegleria fowleri could be easily separated from other species within the genus Naegleria. A number of peaks detected were tentatively identified. MALDI-TOF MS fingerprinting is a rapid, reproducible, high-throughput alternative method for identifying Naegleria isolates. This method has potential for studying eukaryotic agents.

Complex patterns of gene fission in the eukaryotic folate biosynthesis pathway

Shared derived genomic characters can be useful for polarizing phylogenetic relationships, for example, gene fusions have been used to identify deep-branching relationships in the eukaryotes. Here, we report the evolutionary analysis of a three-gene fusion of folB, folK, and folP, which encode enzymes that catalyze consecutive steps in de novo folate biosynthesis. The folK-folP fusion was found across the eukaryotes and a sparse collection of prokaryotes. This suggests an ancient derivation with a number of gene losses in the eukaryotes potentially as a consequence of adaptation to heterotrophic lifestyles. In contrast, the folB-folK-folP gene is specific to a mosaic collection of Amorphea taxa (a group encompassing: Amoebozoa, Apusomonadida, Breviatea, and Opisthokonta). Next, we investigated the stability of this character. We identified numerous gene losses and a total of nine gene fission events, either by break up of an open reading frame (four events identified) or loss of a component domain (five events identified). This indicates that this three gene fusion is highly labile. These data are consistent with a growing body of data indicating gene fission events occur at high relative rates. Accounting for these sources of homoplasy, our data suggest that the folB-folK-folP gene fusion was present in the last common ancestor of Amoebozoa and Opisthokonta but absent in the Metazoa including the human genome. Comparative genomic data of these genes provides an important resource for designing therapeutic strategies targeting the de novo folate biosynthesis pathway of a variety of eukaryotic pathogens such as Acanthamoeba castellanii.

Missing genes, multiple ORFs, and C-to-U type RNA editing in Acrasis kona (Heterolobosea, Excavata) mitochondrial DNA

Discoba (Excavata) is an ancient group of eukaryotes with great morphological and ecological diversity. Unlike the other major divisions of Discoba (Jakobida and Euglenozoa), little is known about the mitochondrial DNAs (mtDNAs) of Heterolobosea. We have assembled a complete mtDNA genome from the aggregating heterolobosean amoeba, Acrasis kona, which consists of a single circular highly AT-rich (83.3%) molecule of 51.5 kb. Unexpectedly, A. kona mtDNA is missing roughly 40% of the protein-coding genes and nearly half of the transfer RNAs found in the only other sequenced heterolobosean mtDNAs, those of Naegleria spp. Instead, over a quarter of A. kona mtDNA consists of novel open reading frames. Eleven of the 16 protein-coding genes missing from A. kona mtDNA were identified in its nuclear DNA and polyA RNA, and phylogenetic analyses indicate that at least 10 of these 11 putative nuclear-encoded mitochondrial (NcMt) proteins arose by direct transfer from the mitochondrion. Acrasis kona mtDNA also employs C-to-U type RNA editing, and 12 homologs of DYW-type pentatricopeptide repeat (PPR) proteins implicated in plant organellar RNA editing are found in A. kona nuclear DNA. A mapping of mitochondrial gene content onto a consensus phylogeny reveals a sporadic pattern of relative stasis and rampant gene loss in Discoba. Rampant loss occurred independently in the unique common lineage leading to Heterolobosea + Tsukubamonadida and later in the unique lineage leading to Acrasis. Meanwhile, mtDNA gene content appears to be remarkably stable in the Acrasis sister lineage leading to Naegleria and in their distant relatives Jakobida.

Genome-wide identification of pathogenicity factors of the free-living amoeba Naegleria fowleri

BACKGROUND

The free-living amoeba Naegleria fowleri is the causative agent of the rapidly progressing and typically fatal primary amoebic meningoencephalitis (PAM) in humans. Despite the devastating nature of this disease, which results in > 97% mortality, knowledge of the pathogenic mechanisms of the amoeba is incomplete. This work presents a comparative proteomic approach based on an experimental model in which the pathogenic potential of N. fowleri trophozoites is influenced by the compositions of different media.

RESULTS

As a scaffold for proteomic analysis, we sequenced the genome and transcriptome of N. fowleri. Since the sequence similarity of the recently published genome of Naegleria gruberi was far lower than the close taxonomic relationship of these species would suggest, a de novo sequencing approach was chosen. After excluding cell regulatory mechanisms originating from different media compositions, we identified 22 proteins with a potential role in the pathogenesis of PAM. Functional annotation of these proteins revealed, that the membrane is the major location where the amoeba exerts its pathogenic potential, possibly involving actin-dependent processes such as intracellular trafficking via vesicles.

CONCLUSION

This study describes for the first time the 30 Mb-genome and the transcriptome sequence of N. fowleri and provides the basis for the further definition of effective intervention strategies against the rare but highly fatal form of amoebic meningoencephalitis.

The perils of pathogen discovery: origin of a novel parvovirus-like hybrid genome traced to nucleic acid extraction spin columns

Next-generation sequencing was used for discovery and de novo assembly of a novel, highly divergent DNA virus at the interface between the Parvoviridae and Circoviridae. The virus, provisionally named parvovirus-like hybrid virus (PHV), is nearly identical by sequence to another DNA virus, NIH-CQV, previously detected in Chinese patients with seronegative (non-A-E) hepatitis. Although we initially detected PHV in a wide range of clinical samples, with all strains sharing ∼99% nucleotide and amino acid identity with each other and with NIH-CQV, the exact origin of the virus was eventually traced to contaminated silica-binding spin columns used for nucleic acid extraction. Definitive confirmation of the origin of PHV, and presumably NIH-CQV, was obtained by in-depth analyses of water eluted through contaminated spin columns. Analysis of environmental metagenome libraries detected PHV sequences in coastal marine waters of North America, suggesting that a potential association between PHV and diatoms (algae) that generate the silica matrix used in the spin columns may have resulted in inadvertent viral contamination during manufacture. The confirmation of PHV/NIH-CQV as laboratory reagent contaminants and not bona fide infectious agents of humans underscores the rigorous approach needed to establish the validity of new viral genomes discovered by next-generation sequencing.

Infections with free-living amebae

Acanthamoeba spp., Balamuthia mandrillaris, and Naegleria fowleri are mitochondria-bearing, free-living eukaryotic amebae that have been known to cause infections of the central nervous system (CNS) of humans and other animals. Several species of Acanthamoeba belonging to several different genotypes cause an insidious and chronic disease, granulomatous amebic encephalitis (GAE), principally in immunocompromised hosts including persons infected with HIV/AIDS. Acanthamoeba spp., belonging to mostly group 2, also cause infection of the human cornea, Acanthamoeba keratitis. Balamuthia mandrillaris causes GAE in both immunocompromised and immunocompetent hosts mostly in the very young or very old individuals. Both Acanthamoeba spp. and B. mandrillaris also cause a disseminated disease including the lungs, skin, kidneys, and uterus. Naegleria fowleri, on the other hand, causes an acute and fulminating, necrotizing infection of the CNS called primary amebic meningoencephalitis (PAM) in children and young adults with a history of recent exposure to warm fresh water. Additionally, another free-living ameba Sappinia pedata, previously described as S. diploidea, also has caused a single case of amebic meningoencephalitis. In this review the biology of these amebae, clinical manifestations, molecular and immunological diagnosis, and epidemiological features associated with GAE and PAM are discussed.