De Novo Transcriptome Profiling of Naegleria fowleri Trophozoites and Cysts via RNA Sequencing

Abstract

Naegleria fowleri is a pathogenic free-living amoeba, commonly found around the world in warm, fresh water and soil. N. fowleri trophozoites can infect humans by entering the brain through the nose and causing usually fatal primary amebic meningoencephalitis (PAM). Trophozoites can encyst to survive under unfavorable conditions such as cold temperature, starvation, and desiccation. Recent technological advances in genomics and bioinformatics have provided unique opportunities for the identification and pre-validation of pathogen-related and environmental resistance through improved understanding of the biology of pathogenic N. fowleri trophozoites and cysts at a molecular level. However, genomic and transcriptomic data on differential expression genes (DEGs) between trophozoites and cysts of N. fowleri are very limited. Here, we report transcriptome Illumina RNA sequencing (RNA-seq) for N. fowleri trophozoites and cysts and de novo transcriptome assembly. RNA-seq libraries were generated from RNA extracted from N. fowleri sampled from cysts, and a reference transcriptome was generated through the assembly of trophozoite data. In the database, the assembly procedure resulted in 42,220 contigs with a mean length of 11,254 nucleotides and a C+G content of 37.21%. RNA sequencing showed that 146 genes in cysts of N. fowleri indicated 2-fold upregulation in comparison with trophozoites of N. fowleri, and 163 genes were downregulated; these genes were found to participate in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. The KEGG pathway included metabolic (131 sequences) and genetic information processing (66 sequences), cellular processing (43 sequences), environmental information processing (22 sequences), and organismal system (20 sequences) pathways. On the other hand, an analysis of 11,254 sequences via the Gene Ontology database showed that their annotations contained 1069 biological processes including the cellular process (228 sequences) and metabolic process (214 sequences); 923 cellular components including cells (240 sequences) and cell parts (225 sequences); and 415 molecular functions including catalytic activities (195 sequences) and binding processes (186 sequences). Differential expression levels increased in cysts of N. fowleri compared to trophozoites of N. fowleri, which were mainly categorized as serine/threonine protease, kinase, and lipid metabolism-related proteins. These results may provide new insights into pathogen-related genes or environment-resistant genes in the pathogenesis of N. fowleri.