InterB multigenic family, a gene repertoire associated with subterminal chromosome regions of Encephalitozoon cuniculi and conserved in several human-infecting microsporidian species

Genomic and phenotypic evolution of nematode-infecting microsporidia

Microsporidia are a large phylum of intracellular parasites that can infect most types of animals. Species in the Nematocida genus can infect nematodes including Caenorhabditis elegans, which has become an important model to study mechanisms of microsporidia infection. To understand the genomic properties and evolution of nematode-infecting microsporidia, we sequenced the genomes of nine species of microsporidia, including two genera, Enteropsectra and Pancytospora, without any previously sequenced genomes. Core cellular processes, including metabolic pathways, are mostly conserved across genera of nematode-infecting microsporidia. Each species encodes unique proteins belonging to large gene families that are likely used to interact with host cells. Most strikingly, we observed one such family, NemLGF1, is present in both Nematocida and Pancytospora species, but not any other microsporidia. To understand how Nematocida phenotypic traits evolved, we measured the host range, tissue specificity, spore size, and polar tube length of several species in the genus. Our phylogenetic analysis shows that Nematocida is composed of two groups of species with distinct traits and that species with longer polar tubes infect multiple tissues. Together, our work details both genomic and trait evolution between related microsporidia species and provides a useful resource for further understanding microsporidia evolution and infection mechanisms.

Telomere-to-Telomere genome assemblies of human-infecting Encephalitozoon species

BACKGROUND

Microsporidia are diverse spore forming, fungal-related obligate intracellular pathogens infecting a wide range of hosts. This diversity is reflected at the genome level with sizes varying by an order of magnitude, ranging from less than 3 Mb in Encephalitozoon species (the smallest known in eukaryotes) to more than 50 Mb in Edhazardia spp. As a paradigm of genome reduction in eukaryotes, the small Encephalitozoon genomes have attracted much attention with investigations revealing gene dense, repeat- and intron-poor genomes characterized by a thorough pruning of molecular functions no longer relevant to their obligate intracellular lifestyle. However, because no Encephalitozoon genome has been sequenced from telomere-to-telomere and since no methylation data is available for these species, our understanding of their overall genetic and epigenetic architectures is incomplete.

METHODS

In this study, we sequenced the complete genomes from telomere-to-telomere of three human-infecting Encephalitozoon spp. -E. intestinalis ATCC 50506, E. hellem ATCC 50604 and E. cuniculi ATCC 50602- using short and long read platforms and leveraged the data generated as part of the sequencing process to investigate the presence of epigenetic markers in these genomes. We also used a mixture of sequence- and structure-based computational approaches, including protein structure prediction, to help identify which Encephalitozoon proteins are involved in telomere maintenance, epigenetic regulation, and heterochromatin formation.

RESULTS

The Encephalitozoon chromosomes were found capped by TTAGG 5-mer telomeric repeats followed by telomere associated repeat elements (TAREs) flanking hypermethylated ribosomal RNA (rRNA) gene loci featuring 5-methylcytosines (5mC) and 5-hemimethylcytosines (5hmC), themselves followed by lesser methylated subtelomeres and hypomethylated chromosome cores. Strong nucleotide biases were identified between the telomeres/subtelomeres and chromosome cores with significant changes in GC/AT, GT/AC and GA/CT contents. The presence of several genes coding for proteins essential to telomere maintenance, epigenetic regulation, and heterochromatin formation was further confirmed in the Encephalitozoon genomes.

CONCLUSION

Altogether, our results strongly support the subtelomeres as sites of heterochromatin formation in Encephalitozoon genomes and further suggest that these species might shutdown their energy-consuming ribosomal machinery while dormant as spores by silencing of the rRNA genes using both 5mC/5hmC methylation and facultative heterochromatin formation at these loci.

Comparative Genomics of Microsporidia

The microsporidia are a phylum of intracellular parasites that represent the eukaryotic cell in a state of extreme reduction, with genomes and metabolic capabilities embodying eukaryotic cells in arguably their most streamlined state. Over the past 20 years, microsporidian genomics has become a rapidly expanding field starting with sequencing of the genome of Encephalitozoon cuniculi, one of the first ever sequenced eukaryotes, to the current situation where we have access to the data from over 30 genomes across 20+ genera. Reaching back further in evolutionary history, to the point where microsporidia diverged from other eukaryotic lineages, we now also have genomic data for some of the closest known relatives of the microsporidia such as Rozella allomycis, Metchnikovella spp. and Amphiamblys sp. Data for these organisms allow us to better understand the genomic processes that shaped the emergence of the microsporidia as a group. These intensive genomic efforts have revealed some of the processes that have shaped microsporidian cells and genomes including patterns of genome expansions and contractions through gene gain and loss, whole genome duplication, differential patterns of invasion and purging of transposable elements. All these processes have been shown to occur across short and longer time scales to give rise to a phylum of parasites with dynamic genomes with a diversity of sizes and organisations.

Enterocytozoon bieneusi of animals-With an 'Australian twist'

Enterocytozoon bieneusi is a microsporidian microorganism that causes intestinal disease in animals including humans. E. bieneusi is an obligate intracellular pathogen, typically causing severe or chronic diarrhoea, malabsorption and/or wasting. Currently, E. bieneusi is recognised as a fungus, although its exact classification remains contentious. The transmission of E. bieneusi can occur from person to person and/or animals to people. Transmission is usually via the faecal-oral route through E. bieneusi spore-contaminated water, environment or food, or direct contact with infected individuals. Enterocytozoon bieneusi genotypes are usually identified and classified by PCR-based sequencing of the internal transcribed spacer region (ITS) of nuclear ribosomal DNA. To date, ~600 distinct genotypes of E. bieneusi have been recorded in ~170 species of animals, including various orders of mammals and reptiles as well as insects in >40 countries. Moreover, E. bieneusi has also been found in recreational water, irrigation water, and treated raw- and waste-waters. Although many studies have been conducted on the epidemiology of E. bieneusi, prevalence surveys of animals and humans are scant in some countries, such as Australia, and transmission routes of individual genotypes and related risk factors are poorly understood. This article/chapter reviews aspects of the taxonomy, biology and epidemiology of E. bieneusi; the diagnosis, treatment and prevention of microsporidiosis; critically appraises the naming system for E. bieneusi genotypes as well as the phylogenetic relationships of these genotypes; provides new insights into the prevalence and genetic composition of E. bieneusi populations in animals in parts of Australia using molecular epidemiological tools; and proposes some areas for future research in the E. bieneusi/microsporidiosis field.

Ricin-B-lectin enhances microsporidia Nosema bombycis infection in BmN cells from silkworm Bombyx mori

Nosema bombycis is an obligate intracellular parasitic fungus that utilizes a distinctive mechanism to infect Bombyx mori Spore germination can be used for host cell invasion; however, the detailed mechanism remains to be elucidated. The ricin-B-lectin (RBL) gene is significantly differentially regulated after N. bombycis spore germination, and NbRBL might play roles in spore germination and infection. In this study, the biological function of NbRBL was examined. Protein sequence analysis showed that NbRBL is a secreted protein that attaches to carbohydrates. The relative expression level of the NbRBL gene was low during the first 30 h post-infection (hpi) in BmN cells, and high expression was detected from 42 hpi. Gene cloning, prokaryotic expression, and antibody preparation for NbRBL were performed. NbRBL was detected in total and secreted proteins using western blot analysis. Subcellular localization analysis showed that NbRBL is an intracellular protein. Spore adherence and infection assays showed that NbRBL could enhance spore adhesion to BmN cells; the proliferative activities of BmN cells incubated with anti-NbRBL were higher than those in negative control groups after N. bombycis infection; and the treatment groups showed less damage from spore invasion. We therefore, propose that NbRBL is released during spore germination, enhances spore adhesion to BmN cells, and contributes to spore invasion.

Subtelomere organization in the genome of the microsporidian Encephalitozoon cuniculi: patterns of repeated sequences and physicochemical signatures

Background

The microsporidian Encephalitozoon cuniculi is an obligate intracellular eukaryotic pathogen with a small nuclear genome (2.9 Mbp) consisting of 11 chromosomes. Although each chromosome end is known to contain a single rDNA unit, the incomplete assembly of subtelomeric regions following sequencing of the genome identified only 3 of the 22 expected rDNA units. While chromosome end assembly remains a difficult process in most eukaryotic genomes, it is of significant importance for pathogens because these regions encode factors important for virulence and host evasion.

Results

Here we report the first complete assembly of E. cuniculi chromosome ends, and describe a novel mosaic structure of segmental duplications (EXT repeats) in these regions. EXT repeats range in size between 3.5 and 23.8 kbp and contain four multigene families encoding membrane associated proteins. Twenty-one recombination sites were identified in the sub-terminal region of E. cuniculi chromosomes. Our analysis suggests that these sites contribute to the diversity of chromosome ends organization through Double Strand Break repair mechanisms. The region containing EXT repeats at chromosome extremities can be differentiated based on gene composition, GC content, recombination sites density and chromosome landscape.

Conclusion

Together this study provides the complete structure of the chromosome ends of E. cuniculi GB-M1, and identifies important factors, which could play a major role in parasite diversity and host-parasite interactions. Comparison with other eukaryotic genomes suggests that terminal regions could be distinguished precisely based on gene content, genetic instability and base composition biais. The diversity of processes assciated with chromosome extremities and their biological consequences, as they are presented in the present study, emphasize the fact that great effort will be necessary in the future to characterize more carefully these regions during whole genome sequencing efforts.

Electronic supplementary material

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

Contrasting host-pathogen interactions and genome evolution in two generalist and specialist microsporidian pathogens of mosquitoes

Obligate intracellular pathogens depend on their host for growth yet must also evade detection by host defenses. Here we investigate host adaptation in two Microsporidia, the specialist Edhazardia aedis and the generalist Vavraia culicis, pathogens of disease vector mosquitoes. Genomic analysis and deep RNA-Seq across infection time courses reveal fundamental differences between these pathogens. E. aedis retains enhanced cell surface modification and signalling capacity, upregulating protein trafficking and secretion dynamically during infection. V. culicis is less dependent on its host for basic metabolites and retains a subset of spliceosomal components, with a transcriptome broadly focused on growth and replication. Transcriptional profiling of mosquito immune responses reveals that response to infection by E. aedis differs dramatically depending on the mode of infection, and that antimicrobial defensins may play a general role in mosquito defense against Microsporidia. This analysis illuminates fundamentally different evolutionary paths and host interplay of specialist and generalist pathogens.

Reduction and expansion in microsporidian genome evolution: new insights from comparative genomics

Microsporidia are an abundant group of obligate intracellular parasites of other eukaryotes, including immunocompromised humans, but the molecular basis of their intracellular lifestyle and pathobiology are poorly understood. New genomes from a taxonomically broad range of microsporidians, complemented by published expression data, provide an opportunity for comparative analyses to identify conserved and lineage-specific patterns of microsporidian genome evolution that have underpinned this success. In this study, we infer that a dramatic bottleneck in the last common microsporidian ancestor (LCMA) left a small conserved core of genes that was subsequently embellished by gene family expansion driven by gene acquisition in different lineages. Novel expressed protein families represent a substantial fraction of sequenced microsporidian genomes and are significantly enriched for signals consistent with secretion or membrane location. Further evidence of selection is inferred from the gain and reciprocal loss of functional domains between paralogous genes, for example, affecting transport proteins. Gene expansions among transporter families preferentially affect those that are located on the plasma membrane of model organisms, consistent with recruitment to plug conserved gaps in microsporidian biosynthesis and metabolism. Core microsporidian genes shared with other eukaryotes are enriched in orthologs that, in yeast, are highly expressed, highly connected, and often essential, consistent with strong negative selection against further reduction of the conserved gene set since the LCMA. Our study reveals that microsporidian genome evolution is a highly dynamic process that has balanced constraint, reductive evolution, and genome expansion during adaptation to an extraordinarily successful obligate intracellular lifestyle.

The genome of Spraguea lophii and the basis of host-microsporidian interactions

Microsporidia are obligate intracellular parasites with the smallest known eukaryotic genomes. Although they are increasingly recognized as economically and medically important parasites, the molecular basis of microsporidian pathogenicity is almost completely unknown and no genetic manipulation system is currently available. The fish-infecting microsporidian Spraguea lophii shows one of the most striking host cell manipulations known for these parasites, converting host nervous tissue into swollen spore factories known as xenomas. In order to investigate the basis of these interactions between microsporidian and host, we sequenced and analyzed the S. lophii genome. Although, like other microsporidia, S. lophii has lost many of the protein families typical of model eukaryotes, we identified a number of gene family expansions including a family of leucine-rich repeat proteins that may represent pathogenicity factors. Building on our comparative genomic analyses, we exploited the large numbers of spores that can be obtained from xenomas to identify potential effector proteins experimentally. We used complex-mix proteomics to identify proteins released by the parasite upon germination, resulting in the first experimental isolation of putative secreted effector proteins in a microsporidian. Many of these proteins are not related to characterized pathogenicity factors or indeed any other sequences from outside the Microsporidia. However, two of the secreted proteins are members of a family of RICIN B-lectin-like proteins broadly conserved across the phylum. These proteins form syntenic clusters arising from tandem duplications in several microsporidian genomes and may represent a novel family of conserved effector proteins. These computational and experimental analyses establish S. lophii as an attractive model system for understanding the evolution of host-parasite interactions in microsporidia and suggest an important role for lineage-specific innovations and fast evolving proteins in the evolution of the parasitic microsporidian lifecycle.

Microsporidia are unusual intracellular parasites that infect a broad range of animal cells. In comparison to their fungal relatives, microsporidian genomes have shrunk during evolution, encoding as few as 2000 proteins. This minimal molecular repertoire makes them a reduced model system for understanding host-parasite interactions. A number of microsporidian genomes have now been sequenced, but the lack of a system for genetic manipulation makes it difficult to translate these data into a better understanding of microsporidian biology. Here we present a deep sequencing project of Spraguea lophii, a fish-infecting microsporidian that is abundantly available from environmental samples. We use our sequence data combined with germination protocols and complex-mix proteomics to identify proteins released by the cell at the earliest stage of germination, representing potential pathogenicity factors. We profile the RNA expression pattern of germinating cells and identify a set of highly transcribed hypothetical genes. Our study provides new insight into the importance of uncharacterized, lineage-specific and/or fast evolving proteins in microsporidia and provides new leads for the investigation of virulence factors in these enigmatic parasites.

Microsporidia and 'the art of living together'

Parasitism, aptly defined as one of the 'living-together' strategies (Trager, 1986), presents a dynamic system in which the parasite and its host are under evolutionary pressure to evolve new and specific adaptations, thus enabling the coexistence of the two closely interacting partners. Microsporidia are very frequently encountered obligatory intracellular protistan parasites that can infect both animals and some protists and are a consummate example of various aspects of the 'living-together' strategy. Microsporidia, relatives of fungi in the superkingdom Opisthokonta, belong to the relatively small group of parasites for which the host cell cytoplasm is the site of both reproduction and maturation. The structural and physiological reduction of their vegetative stage, together with the manipulation of host cell physiology, enables microsporidia to live in the cytosolic environment for most of their life cycle in a way resembling endocytobionts. The ability to form structurally complex spores and the invention and assembly of a unique injection mechanism enable microsporidia to disperse within host tissues and between host organisms, resulting in long-lasting infections. Microsporidia have adapted their genomes to the intracellular way of life, evolved strategies how to obtain nutrients directly from the host and how to manipulate not only the infected cells, but also the hosts themselves. The enormous variability of host organisms and their tissues provide microsporidian parasites a virtually limitless terrain for diversification and ecological expansion. This review attempts to present a general overview of microsporidia, emphasising some less known and/or more recently discovered facets of their biology.

Spraguea lophii (Microsporidia) parasite of the teleost fish, Lophius piscatorius from Tunisian coasts: evidence for an extensive chromosome length polymorphism

A microsporidian of the genus Spraguea was found parasitizing the nervous tissues of Lophius piscatorius collected from various localities in the Mediterranean coastal areas of Tunisia. The tissue localization, the infection focus aspect and sporal dimorphism are characteristics of Spraguea lophii species. Molecular data based on partial sequence of SSUrRNA encoding gene shows few nucleotide polymorphisms, compared to all described Spraguea isolates. Molecular karyotype obtained on pulsed field gel electrophoresis (1D-PFGE) shows a profile with 14 stained bands in the range of 230-880 kbp and a genome size estimated to 6.700 kbp. The rare cutter endonuclease MluI KARD 2-D-PFGE fingerprint shows an extensive chromosome length polymorphism, but the number of chromosome is unchanged and consists of 15 different molecules. The extensive chromosome length polymorphism is associated to a reduced number of genetic events.

Extreme reduction and compaction of microsporidian genomes

Microsporidia are fungi-related obligate intracellular parasites with a highly reduced and compact genome, as for Encephalitozoon species which harbor a genome smaller than 3 Mbp. Genome compaction is reflected by high gene density and, for larger microsporidian genomes, size variation is due to repeat elements that do not drastically affect gene density. Furthermore, these pathogens present strong host dependency illustrated by extensive gene loss. Such adaptations associated with genome compaction induced gene size reduction but also simplification of cellular processes such as transcription. Thus, microsporidia are excellent models for eukaryotic genome evolution and gene expression in the context of host-pathogen relationships.

The complete sequence of the smallest known nuclear genome from the microsporidian Encephalitozoon intestinalis

The genome of the microsporidia Encephalitozoon cuniculi is widely recognized as a model for extreme reduction and compaction. At only 2.9 Mbp, the genome encodes approximately 2,000 densely packed genes and little else. However, the nuclear genome of its sister, Encephalitozoon intestinalis, is even more reduced; at 2.3 Mbp, it represents a 20% reduction from an already severely compacted genome, raising the question, what else can be lost? In this paper, we describe the complete sequence of the E. intestinalis genome and its comparison with that of E. cuniculi. The two species share a conserved gene content, order and density over most of their genomes. The exceptions are the subtelomeric regions, where E. intestinalis chromosomes are missing large gene blocks of sequence found in E. cuniculi. In the remaining gene-dense chromosome 'cores', the diminutive intergenic sequences and introns are actually more highly conserved than the genes themselves, suggesting that they have reached the limits of reduction for a fully functional genome.

A comparison of related genomes provides valuable information about how they evolve. Here, the complete sequence of the smallest known nuclear genome from the microsporidia E. intestinalis is described and compared with its larger sister E. cuniculi, revealing what parts are indispensable in even the most reduced genomes.

Unique physiology of host-parasite interactions in microsporidia infections

Microsporidia are intracellular parasites of all major animal lineages and have a described diversity of over 1200 species and an actual diversity that is estimated to be much higher. They are important pathogens of mammals, and are now one of the most common infections among immunocompromised humans. Although related to fungi, microsporidia are atypical in genomic biology, cell structure and infection mechanism. Host cell infection involves the rapid expulsion of a polar tube from a dormant spore to pierce the host cell membrane and allow the direct transfer of the spore contents into the host cell cytoplasm. This intimate relationship between parasite and host is unique. It allows the microsporidia to be highly exploitative of the host cell environment and cause such diverse effects as the induction of hypertrophied cells to harbour prolific spore development, host sex ratio distortion and host cell organelle and microtubule reorganization. Genome sequencing has revealed that microsporidia have achieved this high level of parasite sophistication with radically reduced proteomes and with many typical eukaryotic pathways pared-down to what appear to be minimal functional units. These traits make microsporidia intriguing model systems for understanding the extremes of reductive parasite evolution and host cell manipulation.

Genomic analyses of the microsporidian Nosema ceranae, an emergent pathogen of honey bees

Recent steep declines in honey bee health have severely impacted the beekeeping industry, presenting new risks for agricultural commodities that depend on insect pollination. Honey bee declines could reflect increased pressures from parasites and pathogens. The incidence of the microsporidian pathogen Nosema ceranae has increased significantly in the past decade. Here we present a draft assembly (7.86 MB) of the N. ceranae genome derived from pyrosequence data, including initial gene models and genomic comparisons with other members of this highly derived fungal lineage. N. ceranae has a strongly AT-biased genome (74% A+T) and a diversity of repetitive elements, complicating the assembly. Of 2,614 predicted protein-coding sequences, we conservatively estimate that 1,366 have homologs in the microsporidian Encephalitozoon cuniculi, the most closely related published genome sequence. We identify genes conserved among microsporidia that lack clear homology outside this group, which are of special interest as potential virulence factors in this group of obligate parasites. A substantial fraction of the diminutive N. ceranae proteome consists of novel and transposable-element proteins. For a majority of well-supported gene models, a conserved sense-strand motif can be found within 15 bases upstream of the start codon; a previously uncharacterized version of this motif is also present in E. cuniculi. These comparisons provide insight into the architecture, regulation, and evolution of microsporidian genomes, and will drive investigations into honey bee–Nosema interactions.

Honey bee colonies are in decline in many parts of the world, in part due to pressures from a diverse assemblage of parasites and pathogens. The range and prevalence of the microsporidian pathogen Nosema ceranae has increased significantly in the past decade. Here we describe the N. ceranae genome, presenting genome traits, gene models and regulatory motifs. N. ceranae has an extremely reduced and AT-biased genome, yet one with substantial numbers of repetitive elements. We identify novel genes that appear to be conserved among microsporidia but undetected outside this phylum, which are of special interest as potential virulence factors for these obligate pathogens. A previously unrecognized motif is found upstream of many start codons and likely plays a role in gene regulation across the microsporidia. These and other comparisons provide insight into the architecture, regulation, and evolution of microsporidian genomes, and provide the first genetic tools for understanding how this pathogen interacts with honey bee hosts.

Genome sequence surveys of Brachiola algerae and Edhazardia aedis reveal microsporidia with low gene densities

BACKGROUND

Microsporidia are well known models of extreme nuclear genome reduction and compaction. The smallest microsporidian genomes have received the most attention, but genomes of different species range in size from 2.3 Mb to 19.5 Mb and the nature of the larger genomes remains unknown.

RESULTS

Here we have undertaken genome sequence surveys of two diverse microsporidia, Brachiola algerae and Edhazardia aedis. In both species we find very large intergenic regions, many transposable elements, and a low gene-density, all in contrast to the small, model microsporidian genomes. We also find no recognizable genes that are not also found in other surveyed or sequenced microsporidian genomes.

CONCLUSION

Our results demonstrate that microsporidian genome architecture varies greatly between microsporidia. Much of the genome size difference could be accounted for by non-coding material, such as intergenic spaces and retrotransposons, and this suggests that the forces dictating genome size may vary across the phylum.