Encephalitozoon cuniculi (Microspora) genome: physical map and evidence for telomere-associated rDNA units on all chromosomes

The protist Aurantiochytrium has universal subtelomeric rDNAs and is a host for mirusviruses

Viruses are the most abundant biological entities in the world's oceans, where they play important ecological and biogeochemical roles. Metagenomics is revealing new groups of eukaryotic viruses, although disconnected from known hosts. Among these are the recently described mirusviruses, which share some similarities with herpesviruses.1 50 years ago, "herpes-type" viral particles2 were found in a thraustochytrid member of the labyrinthulomycetes, a diverse group of abundant and ecologically important marine eukaryotes,3,4 but could not be further characterized by methods then available. Long-read sequencing has allowed us to connect the biology of mirusviruses and thraustochytrids. We sequenced the genome of the genetically tractable model thraustochytrid Aurantiochytrium limacinum ATCC MYA-1381 and found that its 26 linear chromosomes have an extraordinary configuration. Subtelomeric ribosomal DNAs (rDNAs) found at all chromosome ends are interspersed with long repeated sequence elements denoted as long repeated-telomere and rDNA spacers (LORE-TEARS). We identified two genomic elements that are related to mirusvirus genomes. The first is a ∼300-kbp episome (circular element 1 [CE1]) present at a high copy number. Strikingly, the second, distinct, mirusvirus-like element is integrated between two sets of rDNAs and LORE-TEARS at the left end of chromosome 15 (LE-Chr15). Similar to metagenomically derived mirusviruses, these putative A. limacinum mirusviruses have a virion module related to that of herpesviruses along with an informational module related to nucleocytoplasmic large DNA viruses (NCLDVs). CE1 and LE-Chr15 bear striking similarities to episomal and endogenous latent forms of herpesviruses, respectively, and open new avenues of research into marine virus-host interactions.

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.

Impact of Genome Reduction in Microsporidia

Microsporidia represent an evolutionary outlier in the tree of life and occupy the extreme edge of the eukaryotic domain with some of their biological features. Many of these unicellular fungi-like organisms have reduced their genomic content to potentially the lowest limit. With some of the most compacted eukaryotic genomes, microsporidia are excellent model organisms to study reductive evolution and its functional consequences. While the growing number of sequenced microsporidian genomes have elucidated genome composition and organization, a recent increase in complementary post-genomic studies has started to shed light on the impacts of genome reduction in these unique pathogens. This chapter will discuss the biological framework enabling genome minimization and will use one of the most ancient and essential macromolecular complexes, the ribosome, to illustrate the effects of extreme genome reduction on a structural, molecular, and cellular level. We outline how reductive evolution in microsporidia has shaped DNA organization, the composition and function of the ribosome, and the complexity of the ribosome biogenesis process. Studying compacted mechanisms, processes, or macromolecular machines in microsporidia illuminates their unique lifestyle and provides valuable insights for comparative eukaryotic structural biology.

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.

Simplified Transformation of Ostreococcus tauri Using Polyethylene Glycol

Ostreococcustauri is an easily cultured representative of unicellular algae (class Mamiellophyceae) that abound in oceans worldwide. Eight complete 13-22 Mb genomes of phylogenetically divergent species within this class are available, and their DNA sequences are nearly always present in metagenomic data produced from marine samples. Here we describe a simplified and robust transformation protocol for the smallest of these algae (O. tauri). Polyethylene glycol (PEG) treatment was much more efficient than the previously described electroporation protocol. Short (2 min or less) incubation times in PEG gave >104 transformants per microgram DNA. The time of cell recovery after transformation could be reduced to a few hours, permitting the experiment to be done in a day rather than overnight as used in previous protocols. DNA was randomly inserted in the O. tauri genome. In our hands PEG was 20-40-fold more efficient than electroporation for the transformation of O. tauri, and this improvement will facilitate mutagenesis of all of the dispensable genes present in the tiny O. tauri genome.

Microsporidia: Obligate Intracellular Pathogens Within the Fungal Kingdom

Microsporidia are obligate intracellular pathogens related to Fungi. These organisms have a unique invasion organelle, the polar tube, which upon appropriate environmental stimulation rapidly discharges out of the spore, pierces a host cell's membrane, and serves as a conduit for sporoplasm passage into the host cell. Phylogenetic analysis suggests that microsporidia are related to the Fungi, being either a basal branch or sister group. Despite the description of microsporidia over 150 years ago, we still lack an understanding of the mechanism of invasion, including the role of various polar tube proteins, spore wall proteins, and host cell proteins in the formation and function of the invasion synapse. Recent advances in ultrastructural techniques are helping to better define the formation and functioning of the invasion synapse. Over the past 2 decades, proteomic approaches have helped define polar tube proteins and spore wall proteins as well as the importance of posttranslational modifications such as glycosylation in the functioning of these proteins, but the absence of genetic techniques for the manipulation of microsporidia has hampered research on the function of these various proteins. The study of the mechanism of invasion should provide fundamental insights into the biology of these ubiquitous intracellular pathogens that can be integrated into studies aimed at treating or controlling microsporidiosis.

[Nephritis in a Staffordshire Terrier puppy caused by Encephalitozoon cuniculi genotype I]

An 8-week-old puppy was presented to the clinic because of anorexia, polyuria, polydipsia and azotemia. The urinary sediment displayed large amounts of protozoan-like organisms, which could be identified by PCR as Encephalitozoon cuniculi genotype I. In the immunofluorescence antibody test (IFAT) an antibody titer against E. cuniculi of 1  :  2560 was found. The dog was treated with fenbendazole over 3 weeks. After 3 months, an antibody titer against E. cuniculi could no longer be detected. The dog recovered completely. This rare case demonstrates that E. cuniculi genotype I can cause clinical disease in dogs with renal involvement, which can be successfully treated with fenbendazole.

Integrative Approach for the Reliable Detection and Specific Identification of the Microsporidium Loma morhua in Atlantic Cod (Gadus morhua)

Microsporidia are fungal parasites that infect diverse invertebrate and vertebrate hosts. Finfish aquaculture supports epizootics due to high host density and the high biotic potential of these parasites. Reliable methods for parasite detection and identification are a necessary precursor to empirical assessment of strategies to mitigate the effects of these pathogens during aquaculture. We developed an integrative approach to detect and identify Loma morhua infecting Atlantic cod. We show that the spleen is more reliable than the commonly presumed gills as best organ for parasite detection in spite of substantial morphological plasticity in xenoma complexes. We developed rDNA primers with 100% sensitivity in detecting L. morhua and with utility in distinguishing some congeneric Loma species. ITS sequencing is necessary to distinguish L. morhua from other congeneric microsporidia due to intraspecific nucleotide variation. 64% of L. morhua ITS variants from Atlantic cod have a 9-nucleotide motif that distinguishes it from Loma spp. infecting non-Gadus hosts. The remaining 36% of ITS variants from Atlantic cod are distinguished from currently represented Loma spp., particularly those infecting Gadus hosts, based on a 14-nucleotide motif. This research approach is amenable to developing templates in support of reliable detection and identification of other microsporidian parasites in fishes.

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.

Simple method of DNA stretching on glass substrate for fluorescence imaging and spectroscopy

We demonstrate a simple method of stretching DNA to its full length, suitable for optical imaging and atomic force microscopy (AFM). Two competing forces on the DNA molecules, which are the electrostatic attraction between positively charged dye molecules (YOYO-1) intercalated into DNA and the negatively charged surface of glass substrate, and the centrifugal force of the rotating substrate, are mainly responsible for the effective stretching and the dispersion of single strands of DNA. The density of stretched DNA molecules could be controlled by the concentration of the dye-stained DNA solution. Stretching of single DNA molecules was confirmed by AFM imaging and the photoluminescence spectra of single DNA molecule stained with YOYO-1 were obtained, suggesting that our method is useful for spectroscopic analysis of DNA at the single molecule level.

Ribosomal gene polymorphism in small genomes: analysis of different 16S rRNA sequences expressed in the honeybee parasite Nosema ceranae (Microsporidia)

To date, few organisms have been shown to possess variable ribosomal RNA, otherwise considered a classic example of uniformity by concerted evolution. The polymorphism for the 16S rRNA in Nosema ceranae analysed here is striking as Microsporidia are intracellular parasites which have suffered a strong reduction in their genomes and cellular organization. Moreover, N. ceranae infects the honeybee Apis mellifera, and has been associated with the colony-loss phenomenon during the last decade. The variants of 16S rRNA include single nucleotide substitutions, one base insertion-deletion, plus a tetranucleotide indel. We show that different gene variants are expressed. The polymorphic sites tend to be located in particular regions of the rRNA molecule, and the comparison to the Escherichia coli 16S rRNA secondary structure indicates that most variations probably do not preclude ribosomal activity. The fact that the polymorphisms in such a minimal organism as N. ceranae are maintained in samples collected worldwide suggest that the existence of differently expressed rRNA may play an adaptive role in the microsporidian.

Intraspecific polymorphism of rDNA among five Nosema bombycis isolates from different geographic regions in China

The microsporidian Nosema bombycis is the causative agent of pébrine, a highly infectious disease of the silkworm Bombyx mori. Three regions of the multicopy rDNA gene were examined in order to investigate the relationships among five Nosema isolates from various regions of China. Ribosomal DNA alleles are present on each of the 18 chromosomes of N. bombycis and show a high degree of variation. In this study the small subunit (SSU) rDNA, internal transcribed spacer (ITS) and intergenic spacer (IGS) regions for up to 10 different rDNA copies from each N. bombycis isolate were cloned and sequenced. As expected we see greater polymorphism in the ITS region (88 variable sites in 179 nucleotides) and IGS (200 variable sites in 279 nucleotides) than in the SSU rDNA (24 variable sites in 1232 nucleotides). Phylogenetic analysis shows greater differences between alleles within an isolate than between the same alleles from different isolates. The data reveal two very different groups, one from the Sichuan province and the other with a broad distribution including four provinces in southeast China and Japan. The Sichuan isolate does not have any rDNA alleles with sequences identical to those in the other isolates, implying that it is a separate, non-intermixing, population or perhaps a separate species from the other isolates. In light of the polymorphic nature of the rDNA alleles in N. bombycis and their presence on every chromosome, the rDNA gene may be useful for understanding the movement and ultimately the source of pébrine infections.

Extremely reduced levels of heterozygosity in the vertebrate pathogen Encephalitozoon cuniculi

The genomes of microsporidia in the genus Encephalitozoon have been extensively studied for their minimalistic features, but they have seldom been used to investigate basic characteristics of the biology of these organisms, such as their ploidy or their mode of reproduction. In the present study, we aimed to tackle this issue by mapping Illumina sequence reads against the genomes of four strains of E. cuniculi. This approach, combined with more conventional molecular biology techniques, resulted in the identification of heterozygosity in all strains investigated, a typical signature of a diploid nuclear state. In sharp contrast with similar studies recently performed on a distant microsporidian lineage (Nematocida spp.), the level of heterozygosity that we identified across the E. cuniculi genomes was found to be extremely low. This reductive intraindividual genetic variation could result from the long-term propagation of these strains under laboratory conditions, but we propose that it could also reflect an intrinsic capacity of these vertebrate pathogens to self-reproduce.

Characterization of active ribosomal RNA harboring MITEs insertion in microsporidian Nosema bombycis genome

Microsporidia are a group of obligate intracellular parasites of medical and agricultural importance, which can infect almost all animals, including human beings. Using the genome data of Nosema bombycis, four families of miniature inverted-repeat transposable elements (MITEs) in ribosomal DNA (rDNA) were characterized in the microsporidian N. bombycis and were named LSUME1, ITSME1, SSUME1, and SSUME2, respectively. The genome-wide investigation of these MITEs shows that these MITEs families distribute randomly in N. bombycis genome. All insertion sequences have conserved characteristics of MITEs, the direct repeat sequence and terminal inverted-repeat sequence at both ends of each MITEs sequence. Additionally, using the CLC RNA Workbench Software, secondary structures of rRNA containing MITEs sequence have been predicted and were located in variable region or expansion segment. Furthermore, using two different probes, one is prepared by MITE sequence only (short probe) and the other is prepared by MITE sequence flanking partial rDNA sequence (long probe); northern blotting and dot blotting have been performed to detect the transcriptional and functional activity of the rDNA containing MITEs insertion. Fortunately, we found that the rDNA, which harbors the MITE, not only can be transcripted but also can form a complete ribosome. This is an interesting thing that one gene can keep active even when it has been inserted with another sequence. But the biological and structural significance of this observation is not readily apparent.

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.

Polymorphism and recombination for rDNA in the putatively asexual microsporidian Nosema ceranae, a pathogen of honeybees

Nosema ceranae is currently one of the major pathogens of honeybees, related to the worldwide colony losses phenomenon. The genotyping of strains based on ribosomal DNA (rDNA) can be misleading if the repeated units are not identical. The analysis of cloned rDNA fragments containing the intergenic spacer (IGS) and part of the rDNA small-subunit (SSU) gene, from N. ceranae isolates from different European and Central Asia populations, revealed a high diversity of sequences. The variability involved single-nucleotide polymorphisms and insertion/deletions, resulting in 79 different haplotypes. Two sequences from the same isolate could be as different as any pair of sequences from different samples; in contrast, identical haplotypes were also found in very different geographical origins. Consequently, haplotypes cannot be organized in a consistent phylogenetic tree, clearly indicating that rDNA is not a reliable marker for the differentiation of N. ceranae strains. The results indicate that recombination between different sequences may produce new variants, which is quite surprising in microsporidia, usually considered to have an asexual mode of reproduction. The diversity of sequences and their geographical distribution indicate that haplotypes of different lineages may occasionally be present in a same cell and undergo homologue recombination, therefore suggesting a sexual haplo-diploid cycle.

The intriguing nature of microsporidian genomes

Microsporidia are a group of highly adapted unicellular fungi that are known to infect a wide range of animals, including humans and species of great economic importance. These organisms are best known for their very simple cellular and genomic features, an adaptive consequence of their obligate intracellular parasitism. In the last decade, the acquisition of a large amount of genomic and transcriptomic data from several microsporidian species has greatly improved our understanding of the consequences of a purely intracellular lifestyle. In particular, genome sequence data from these pathogens has revealed how obligate intracellular parasitism can result in radical changes in the composition and structure of nuclear genomes and how these changes can affect cellular and evolutionary mechanisms that are otherwise well conserved among eukaryotes. This article reviews our current understanding of the genome content and structure of microsporidia, discussing their evolutionary origin and cataloguing the mechanisms that have often been involved in their extreme reduction.

Genomic survey of the non-cultivatable opportunistic human pathogen, Enterocytozoon bieneusi

Enterocytozoon bieneusi is the most common microsporidian associated with human disease, particularly in the immunocompromised population. In the setting of HIV infection, it is associated with diarrhea and wasting syndrome. Like all microsporidia, E. bieneusi is an obligate, intracellular parasite, but unlike others, it is in direct contact with the host cell cytoplasm. Studies of E. bieneusi have been greatly limited due to the absence of genomic data and lack of a robust cultivation system. Here, we present the first large-scale genomic dataset for E. bieneusi. Approximately 3.86 Mb of unique sequence was generated by paired end Sanger sequencing, representing about 64% of the estimated 6 Mb genome. A total of 3,804 genes were identified in E. bieneusi, of which 1,702 encode proteins with assigned functions. Of these, 653 are homologs of Encephalitozoon cuniculi proteins. Only one E. bieneusi protein with assigned function had no E. cuniculi homolog. The shared proteins were, in general, evenly distributed among the functional categories, with the exception of a dearth of genes encoding proteins associated with pathways for fatty acid and core carbon metabolism. Short intergenic regions, high gene density, and shortened protein-coding sequences were observed in the E. bieneusi genome, all traits consistent with genomic compaction. Our findings suggest that E. bieneusi is a likely model for extreme genome reduction and host dependence.

Enterocytozoon bieneusi is a clinically significant pathogen associated with human microsporidiosis, particularly in immunocompromised individuals. E. bieneusi is widespread in mammals, and there is no effective commercial treatment for infection. The pathogen cannot be readily cultivated, and animal models are limited. We therefore undertook a sequence survey and generated the first large-scale genomic dataset for E. bieneusi, which we used to study the organization and structure of its genome and to perform a comparative analysis with Encephalitozoon cuniculi, another microsporidian whose genome has been completely sequenced. The E. bieneusi genome showed many traits associated with genome compaction including high gene density, short intergenic regions, shortened proteins, and few introns. With one exception, all E. bieneusi proteins with assigned functions had E. cuniculi homologs. We found a paucity of genes encoding proteins associated with fatty acid and carbon metabolism. The possibility that these core functions are reduced in an intracellular parasite is intriguing, but because the genome sequence of E. bieneusi is incomplete, we cannot exclude the possibility that additional proteins associated with the various metabolic pathways would be discovered in a completed genome.

Multiple rDNA units distributed on all chromosomes of Nosema bombycis

Among Microsporidia, Nosema bombycis has a novel arrangement of LSUrRNA, SSUrRNA, ITS, IGS and 5SrRNA. To determine the distribution of rDNA among the chromosomes, we performed genome-wide screening and Southern blotting with three probes (SSU, ITS and IGS). Southern blotting revealed that ribosomal RNA genes are distributed on all chromosomes of N. bombycis, which is contrary to the previous result, which concluded that the N. bombycis rRNA genes were limited to a single chromosome. This wide distribution is similar to that of the rDNA unit of Encephalitozoon cuniculi. Screening of the N. bombycis genome detected 53 LSUrRNA elements, 43 SSUrRNA elements and 36 5SrRNA elements. However, it is still difficult to determine their loci on the chromosomes as the genomic map is unfinished.

Genome diversity in microbial eukaryotes

The genomic peculiarities among microbial eukaryotes challenge the conventional wisdom of genome evolution. Currently, many studies and textbooks explore principles of genome evolution from a limited number of eukaryotic lineages, focusing often on only a few representative species of plants, animals and fungi. Increasing emphasis on studies of genomes in microbial eukaryotes has and will continue to uncover features that are either not present in the representative species (e.g. hypervariable karyotypes or highly fragmented mitochondrial genomes) or are exaggerated in microbial groups (e.g. chromosomal processing between germline and somatic nuclei). Data for microbial eukaryotes have emerged from recent genome sequencing projects, enabling comparisons of the genomes from diverse lineages across the eukaryotic phylogenetic tree. Some of these features, including amplified rDNAs, subtelomeric rDNAs and reduced genomes, appear to have evolved multiple times within eukaryotes, whereas other features, such as absolute strand polarity, are found only within single lineages.

Novel nucleomorph genome architecture in the cryptomonad genus hemiselmis

Cryptomonads are ubiquitous aquatic unicellular eukaryotes that acquired photosynthesis through the uptake and retention of a red algal endosymbiont. The nuclear genome of the red alga persists in a highly reduced form termed a nucleomorph. The nucleomorph genome of the model cryptomonad Guillardia theta has been completely sequenced and is a mere 551 kilobases (kb) in size, spread over three chromosomes. The presence of three chromosomes appears to be a universal characteristic of nucleomorph genomes in cryptomonad algae as well as in the chlorarachniophytes, an unrelated algal lineage with a nucleomorph and plastid genome derived from a green algal endosymbiont. Another feature of nucleomorph genomes in all cryptomonads and chlorarachniophytes examined thus far is the presence of subtelomeric ribosomal DNA (rDNA) repeats at the ends of each chromosome. Here we describe the first exception to this canonical nucleomorph genome architecture in the cryptomonad Hemiselmis rufescens CCMP644. Using pulsed-field gel electrophoresis (PFGE), we estimate the size of the H. rufescens nucleomorph genome to be approximately 580 kb, slightly larger than the G. theta genome. Unlike the situation in G. theta and all other known cryptomonads, sub-telomeric repeats of the rDNA cistron appear to be absent on both ends of the second largest chromosome in H. rufescens and two other members of this genus. Southern hybridizations using a variety of nucleomorph protein gene probes against PFGE-separated H. rufescens chromosomes indicate that recombination has been a major factor in shaping the karyotype and genomic structure of cryptomonad nucleomorphs.

Stripped-down DNA repair in a highly reduced parasite

Background

Encephalitozoon cuniculi is a member of a distinctive group of single-celled parasitic eukaryotes called microsporidia, which are closely related to fungi. Some of these organisms, including E. cuniculi, also have uniquely small genomes that are within the prokaryotic range. Thus, E. cuniculi has undergone a massive genome reduction which has resulted in a loss of genes from diverse biological pathways, including those that act in DNA repair.

DNA repair is essential to any living cell. A loss of these mechanisms invariably results in accumulation of mutations and/or cell death. Six major pathways of DNA repair in eukaryotes include: non-homologous end joining (NHEJ), homologous recombination repair (HRR), mismatch repair (MMR), nucleotide excision repair (NER), base excision repair (BER) and methyltransferase repair. DNA polymerases are also critical players in DNA repair processes.

Given the close relationship between microsporidia and fungi, the repair mechanisms present in E. cuniculi were compared to those of the yeast Saccharomyces cerevisiae to ascertain how the process of genome reduction has affected the DNA repair pathways.

Results

E. cuniculi lacks 16 (plus another 6 potential absences) of the 56 DNA repair genes sought via BLASTP and PSI-BLAST searches. Six of 14 DNA polymerases or polymerase subunits are also absent in E. cuniculi. All of these genes are relatively well conserved within eukaryotes. The absence of genes is not distributed equally among the different repair pathways; some pathways lack only one protein, while there is a striking absence of many proteins that are components of both double strand break repair pathways. All specialized repair polymerases are also absent.

Conclusion

Given the large number of DNA repair genes that are absent from the double strand break repair pathways, E. cuniculi is a prime candidate for the study of double strand break repair with minimal machinery. Strikingly, all of the double strand break repair genes that have been retained by E. cuniculi participate in other biological pathways.

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

Microsporidia are fungi-related obligate intracellular parasites that infect numerous animals, including man. Encephalitozoon cuniculi harbours a very small genome (2.9 Mbp) with about 2,000 coding sequences (CDSs). Most repeated CDSs are of unknown function and are distributed in subterminal regions that mark the transitions between subtelomeric rDNA units and chromosome cores. A potential multigenic family (interB) encoding proteins within a size range of 579-641 aa was investigated by PCR and RT-PCR. Thirty members were finally assigned to the E. cuniculi interB family and a predominant interB transcript was found to originate from a newly identified gene on chromosome III. Microsporidian species from eight different genera infecting insects, fishes or mammals, were tested for a possible intra-phylum conservation of interB genes. Only representatives of the Encephalitozoon, Vittaforma and Brachiola genera, differing in host range but all able to invade humans, were positive. Molecular karyotyping of Brachiola algerae showed a complex set of chromosome bands, providing a haploid genome size estimate of 15-20 Mbp. In spite of this large difference in genome complexity, B. algerae and E. cuniculi shared some similar interB gene copies and a common location of interB genes in near-rDNA subterminal regions.

Multiple rRNA Variants in a Single Spore of the Microsporidian Nosema bombi

To understand the source of the multiple DNA sequence variants of Nosema bombi ribosomal RNA (rRNA) found in a single bumble bee host, we PCR amplified, cloned, and sequenced the partial rRNA gene from 125 clones, which were derived from four out of 46 spores individually isolated from a single host by laser microdissection. At least two rRNA variants, characterized by either (GTTT)(2) or (GTTT)(3) repeat units within the internal transcribed spacer (ITS) region, were found per spore in approximately equal proportions, variants which were also found in approximately equal proportions in 55 clones of the two DNA extracts of multiple spores from the same host. Firstly, we demonstrate for the first time that DNA sequences can be obtained from single-binucleate microsporidia. Secondly, it appears that concerted evolution has not homogenized the sequences of all rRNA copies within a single N. bombi spore or even within a single nucleus. We thereby demonstrate unequivocally that two or more rRNA sequence variants exist per N. bombi spore, and urge caution in the use of multicopy rRNA genes for population genetic and phylogenetic analysis of this and other Microsporidia unless homologous copies can be reliably typed.

Complete rRNA Gene Sequences Reveal that the MicrosporidiumNosema bombiInfects Diverse Bumblebee (Bombusspp.) Hosts and Contains Multiple Polymorphic Sites

Characterisation of microsporidian species and differentiation among genetic variants of the same species has typically relied on ribosomal RNA (rRNA) gene sequences. We characterised the entire rRNA gene of a microsporidium from 11 isolates representing eight different European bumblebee (Bombus) species. We demonstrate that the microsporidium Nosema bombi infected all hosts that originated from a wide geographic area. A total of 16 variable sites (all single nucleotid polymorphisms (SNPs)) was detected in the small subunit (SSU) rRNA gene and 42 (39 SNPs and 3 indels) in the large subunit (LSU) rRNA sequence. Direct sequencing of PCR-amplified DNA products of the internal transcribed spacer (ITS) region revealed identical sequences in all isolates. In contrast, ITS fragment length determined by PAGE and sequencing of cloned amplicons gave better resolution of sequences and revealed multiple SNPs across isolates and two fragment sizes in each isolate (six short and seven long amplicon variants). Genetic variants were not unique to individual host species. Moreover, two or more sequence variants were obtained from individual bumblebee hosts, suggesting the existence of multiple, variable copies of rRNA in the same microsporidium, and contrary to that expected for a class of multi-gene family under concerted evolution theory. Our data on within-genome rRNA variability call into question the usefulness of rRNA sequences to characterise intraspecific genetic variants in the Microsporidia and other groups of unicellular organisms.

Combing DNA on CTAB-coated surfaces

A fluorescence microscope (FM) coupled with an intensified charge-coupled device (ICCD) camera was used to investigate the combing of DNA on cetyltrimethyl ammonium bromide (CTAB)-coated glass surfaces. DNA molecules can be combed uniform and straight on CTAB-coated surfaces. Different combing characteristics at different pH values were found. At lower pH (ca. 5.5), DNA molecules were stretched 30% longer than the unextended and DNA extremities bound with CTAB-coated surfaces via hydrophobic interaction. At high pH values (e.g., 6.4 and 6.5), DNA molecules were extended about 10% longer and DNA extremities bound with CTAB-coated surfaces via electrostatic attraction. At pH 6.0, DNA molecules could be extended 30% longer on 0.2-mM CTAB-coated surfaces. CTAB cationic surfactant has both a hydrophobic motif and a positively charged group. So, CTAB-coated surfaces can bind DNA extremities via hydrophobic effect or electrostatic attraction at different pH values. It was also found that combing of DNA on CTAB-coated surfaces is reversible. The number of DNA base pairs binding to CTAB-coated surfaces was calculated.

Analyses of the Ribosomal DNA Region in Nosema bombycis NIS 001

Ribosomal DNA (rDNA) containing small subunit (SSU) rDNA and both flanking regions in the entomopathogenic microsporidian Nosema bombycis NIS 001 was amplified from genomic DNA with a primer set based on the sequence of an inverse polymerase chain reaction (PCR)-derived fragment. In this fragment, SSU rDNA was divided by a 618-bp insert at nt 599, and 5S rDNA was located downstream of the SSU rDNA, fragmented by 284-bp intergenic spacer. In addition, the 48-bp 3'-end of large subunit (LSU) rDNA was located 118 bp upstream of the fragmented SSU rDNA. In the amplicon, the region upstream of the LSU rDNA was a homologue of the C-terminal CHARLIE8 transposon-like element of human GTF2IRD2. In this organism, another fragmented SSU rDNA, which was divided by a 231-bp insert at nt 50, was also detected. Both the intact (insertless) and fragmented SSU rDNAs clustered with LSU rDNA and 5S rDNA and the intergenic sequences between SSU rDNA and 5S rDNA were divergent in an organism. Reverse transcription (RT)-PCR assay indicated that not only the intact SSU rDNA but also the fragmened SSU rDNA were transcribed in N. bombycis.

Description of an ultrathin multiwire proportional chamber-based detector and application to the characterization of theSpraguea lophii(Microspora) two-dimensional genome fingerprint

Multiwire proportional chamber is a useful technology to build detectors that supersede the lack of interactivity of autoradiography in molecular biology experiments. Some drawbacks still limited the diffusion of existing instruments in biological laboratories. The major competitors are storage phosphor imaging systems. The simplified description of a radio-chromato-imager prototype (RCI) based on an original ultrathin multiwire proportional chamber is presented. It combines the advantage of the different existing technologies to present competitive properties in terms of efficiency, spatial resolution, robustness, manipulation easiness and production cost. Application of the RCI detector to molecular biology was performed by the analysis of karyotype and restriction display two-dimensional pulsed-field gel electrophoresis (KARD 2-D PFGE) data which are used to describe small eukaryotic genome structures. The comparative analysis with autoradiography was performed with the PDQuest software on Spraguea lophii (Microspora) genome fingerprints. The spot detection procedure applied to the different images leads to a similar conclusion considering the genome structure of S. lophii which appeared to be composed of 15 chromosomes for 13 karyotypic bands (200-880 kbp).

Inter-strain variability of insertion/deletion events in the Encephalitozoon cuniculi genome: a comparative KARD-PFGE analysis

We applied a two-dimensional pulsed-field gel electrophoresis procedure to the genomes of two karyotype variants assigned to two different strains of the microsporidian Encephalitozoon cuniculi, termed D (strain III) and F (strain II). Data obtained for BssHII and MluI restriction fragment length polymorphisms in each chromosome are compiled and compared to the reference strain I variant A. Six Insertion/Deletion (InDels) are found in subterminal position, some of these being characteristic of either D or F. Like in strain 1, the terminal fragments extending between each telomere and rDNA locus are conserved in length for each chromosome. They are however smaller than in reference variant. This size reduction is estimated to be 2.5 kbp for the strain III isolate and 3.5 kbp for the strain II isolate. We hypothesize that for the three E. cuniculi strains, all chromosome extremities are prone to a constant process of sequence homogenization through mitotic recombination between conserved regions.

Sequence survey of the genome of the opportunistic microsporidian pathogen, Vittaforma corneae

The microsporidian Vittaforma corneae has been reported as a pathogen of the human stratum corneum, where it can cause keratitis, and is associated with systemic infections. In addition to this direct role as an infectious, etiologic agent of human disease, V. corneae has been used as a model organism for another microsporidian, Enterocytozoon bieneusi, a frequent and problematic pathogen of HIV-infected patients that, unlike V. corneae, is difficult to maintain and to study in vitro. Unfortunately, few molecular sequences are available for V. corneae. In this study, seventy-four genome survey sequences (GSS) were obtained from genomic DNA of spores of laboratory-cultured V. corneae. Approximately, 41 discontinuous kilobases of V. corneae were cloned and sequenced to generate these GSS. Putative identities were assigned to 44 of the V. corneae GSS based on BLASTX searches, representing 21 discrete proteins. Of these 21 deduced V. corneae proteins, only two had been reported previously from other microsporidia (until the recent report of the Encephalitozoon cuniculi genome). Two of the V. corneae proteins were of particular interest, reverse transcriptase and topoisomerase IV (parC). Since the existence of transposable elements in microsporidia is controversial, the presence of reverse transcriptase in V. corneae will contribute to resolution of this debate. The presence of topoisomerase IV was remarkable because this enzyme previously had been identified only from prokaryotes. The 74 GSS included 26.7 kilobases of unique sequences from which two statistics were generated: GC content and codon usage. The GC content of the unique GSS was 42%, lower than that of another microsporidian, E. cuniculi (48% for protein-encoding regions), and substantially higher than that predicted for a third microsporidian, Spraguea lophii (28%). A comparison using the Pearson correlation coefficient showed that codon usage in V. corneae was similar to that in the yeasts, Saccharomyces cerevisiae (r = 0.79) and Shizosaccharomyces pombe (r = 0.70), but was markedly dissimilar to E. cuniculi (r = 0.19).

Combing the genome for genomic instability

Genomic instability is one of the major features of cancer cells. The clinical phenotypes associated with several human diseases have been linked to recurrent DNA rearrangements and dysfunction of DNA replication processes that involve unstable genomic regions. Analysis of these rearrangements, which are frequently submicroscopic and can lead to loss or gain of dosage-sensitive genes or gene disruption, requires the development of sensitive, high-resolution techniques. This will lead to a better understanding of the mechanisms underlying genome instability and a greater awareness of the role of chromosomal rearrangements in disease. A new technology that involves molecular combing, a method that permits straightening and aligning molecules of genomic DNA, should make possible a detailed analysis of genomic events at the level of single DNA molecules. Such a single molecule approach could help to elucidate important properties that are masked in bulk studies.

Fine analysis of the Pneumocystis carinii f. sp. carinii genome by two-dimensional pulsed-field gel electrophoresis

Pneumocystis carinii is a general designation for a group of unusual unicellular fungal parasites responsible of pneumopathy in animal hosts. Divided into several subgroups termed the 'special forms', P. carinii is prone to an extensive karyotype variation. In previous studies, the nuclear genome of these organisms has been considered to be haploid and a set of 16 chromosomes has been assigned to P. carinii f. sp. carinii, a special form known to infect rats. We report the analysis of the genome of an isolate representative of the karyotype 1 of this special form, using two-dimensional pulsed-field gel electrophoresis procedures. The 'karyotype and restriction display' (KARD) fingerprints indicated the presence of 17 different chromosomes. The haploid genome size was estimated to be 8.4 Mbp. Some homologous chromosomes were distinguished on the basis of a single restriction fragment length polymorphism, which raises the possibility of a diploid nucleus. A restriction map of the chromosome 15, characterized by two homologues with a size difference of 7 kb, was constructed. Hybridization data indicated that insertion/deletion events may have occurred within subtelomeric regions which carry genes encoding the major surface glycoprotein (MSG) of Pneumocystis.

Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi

Microsporidia are obligate intracellular parasites infesting many animal groups. Lacking mitochondria and peroxysomes, these unicellular eukaryotes were first considered a deeply branching protist lineage that diverged before the endosymbiotic event that led to mitochondria. The discovery of a gene for a mitochondrial-type chaperone combined with molecular phylogenetic data later implied that microsporidia are atypical fungi that lost mitochondria during evolution. Here we report the DNA sequences of the 11 chromosomes of the approximately 2.9-megabase (Mb) genome of Encephalitozoon cuniculi (1,997 potential protein-coding genes). Genome compaction is reflected by reduced intergenic spacers and by the shortness of most putative proteins relative to their eukaryote orthologues. The strong host dependence is illustrated by the lack of genes for some biosynthetic pathways and for the tricarboxylic acid cycle. Phylogenetic analysis lends substantial credit to the fungal affiliation of microsporidia. Because the E. cuniculi genome contains genes related to some mitochondrial functions (for example, Fe-S cluster assembly), we hypothesize that microsporidia have retained a mitochondrion-derived organelle.

Molecular characteristics and physiology of microsporidia

A survey of the molecular features of microsporidia is presented which attempts to comment on unresolved questions concerning the physiology of these amitochondrial intracellular parasites. Various transports of host-derived molecules can be predicted and trehalose appears as a potential reserve of glucose for energy metabolism. Significant insights into membrane lipids, polyamine metabolism and sporogony-specific proteins have been gained. Some species, such as Encephalitozoon cuniculi, are heterogeneous entities and harbor a small genome. Although showing a variation in genome size of 8.5-fold, microsporidia share reduced rDNA genes. Finally, data on gene organization and a possible evolutionary relationship with fungi are considered.

Microsporidia: emerging pathogenic protists

Microsporidia are eukaryotic spore forming obligate intracellular protozoan parasites first recognized over 100 years ago. These organisms infect all of the major animal groups and are now recognized as opportunistic pathogens of humans. Microsporidian spores are common in the environment and microsporidia pathogenic to humans have been found in water supplies. The genera Nosema, Vittaforma, Brachiola, Pleistophora, Encephalitozoon, Enterocytozoon, Septata (reclassified to Encephalitozoon) and Trachipleistophora have been found in human infections. These organisms have the smallest known eukaryotic genomes. Microsporidian ribosomal RNA sequences have proven useful as diagnostic tools as well as for phylogenetic analysis. Recent phylogenetic analysis suggests that Microsporidia are related to the fungi. These organisms are defined by the presence of a unique invasion organelle consisting of a single polar tube that coils around the interior of the spore. All microsporidia exhibit the same response to stimuli, that is, the polar tube discharges from the anterior pole of the spore in an explosive reaction. If the polar tube is discharged next to a cell, it can pierce the cell and transfer its sporoplasm into the cell. A technique was developed for the purification of polar tube proteins (PTPs) using differential extraction followed by reverse phase HPLC. This method was used to purify the PTPs from Glugea americanus, Encephalitozoon cuniculi, Enc. hellem and Enc. intestinalis. These PTPs demonstrate conserved characteristics such as solubility, hydrophobicity, mass, proline content and immunologic epitopes. The major PTP gene from Enc. cuniculi and Enc. hellem has been cloned and expressed in vitro. The gene sequences support the importance of ER and in the formation of the polar tube as suggested by morphologic studies. Analysis of the cloned proteins also indicates that secondary structural characteristics are conserved. These characteristics are probably important in the function of this protein during the eversion/assembly of the polar tube and in providing elasticity and resiliency for sporoplasm passage.

Sequence and analysis of chromosome I of the amitochondriate intracellular parasite Encephalitozoon cuniculi (Microspora)

A DNA sequencing program was applied to the small (<3 Mb) genome of the microsporidian Encephalitozoon cuniculi, an amitochondriate eukaryotic parasite of mammals, and the sequence of the smallest chromosome was determined. The approximately 224-kb E. cuniculi chromosome I exhibits a dyad symmetry characterized by two identical 37-kb subtelomeric regions which are divergently oriented and extend just downstream of the inverted copies of an 8-kb duplicated cluster of six genes. Each subtelomeric region comprises a single 16S-23S rDNA transcription unit, flanked by various tandemly repeated sequences, and ends with approximately 1 kb of heterogeneous telomeric repeats. The central (or core) region of the chromosome harbors a highly compact arrangement of 132 potential protein-coding genes plus two tRNA genes (one gene per 1.14 kb). Most genes occur as single copies with no identified introns. Of these putative genes, only 53 could be assigned to known functions. A number of genes from the transcription and translation machineries as well as from other cellular processes display characteristic eukaryotic signatures or are clearly eukaryote-specific.

The microsporidian Encephalitozoon

Encephalitozoon cuniculi is an attractive model system for amitochondriate intracellular eukaryotic parasites. It is characterized by a very small genome (below 3 Mbp) and a unique invasion apparatus. Furthermore, as an infectious agent, it is important in human and veterinary medicine. The compactness of its genome involves the reduction of rDNA sequences as well as of some protein-coding genes and intergenic regions. Its highly differentiated apparatus to penetrate the host cell, an extrusome-like polar tube, is composed of novel proteins and may permit various pathways of infestation. Completion of the systematic E. cuniculi sequencing project should provide an important reference system for the comparative genomics of amitochondriate and mitochondriate parasites. Further analysis of orphan genes should help to identify factors that are responsible for its intracellular parasitic way of life.

Towards the minimal eukaryotic parasitic genome

Microsporidia are well-known to infect immunocompromised patients and are also responsible for clinical syndromes in immunocompetent individuals. In recent years, evidence has been obtained in support of a very close relationship between Microsporidia and Fungi. In some species, the compaction of the genome and genes is remarkable. Thus, a systematic sequencing project has been initiated for the 2.9 Mbp genome of Encephalitozoon cuniculi, which will be useful for future comparative genomic studies.

Occurence of subtelomeric rearrangements in the genome of the microsporidian parasite Encephalitozoon cuniculi, as revealed by a new fingerprinting procedure based on two-dimensional pulsed field gel electrophoresis

In Microsporidia, mitochondria-lacking eukaryotic intracellular parasites, genomic comparisons were so far based on molecular karyotyping. The mammal-infecting species Encephalitozoon cuniculi is characterized by a very low haploid genome size (approximately 2.8 Mbp) and rather high karyotype variability. Recently, we developed a two-dimensional pulsed field gel electrophoresis (2-D PFGE) fingerprinting technique useful for constructing a restriction map fo the genome of a mouse E. cuniculi isolate (karyotype variant A). The so-called karyotype and restriction display 2-D PFGE (KARD-PFGE) protocol involved 1-D chromosome separation, digestion with a rare cutter, Klenow radiolabeling of genomic DNA and 2-D separation of restriction fragments followed by autoradiography. In order to assess its suitability for detecting polymorphic loci in E. cuniculi, we applied KARD-PFGE with either BssHII or Mlul digestion to genome analysis of two rabbit isolates representative of two different karyotype variants (A and C). The 2-D spot pattern of the rabbit isolate variant A is identical to the reference mouse isolate but differs greatly from the rabbit isolate variant C. Chromosomal restriction fragment length polymorphisms (RFLPs) provide strong evidence for homologous chromosomes and frequent DNA rearrangements within subtelomeric regions just upstream of the dispersed rDNA units closely associated with each chromosomal end.

In-gel DNA radiolabelling and two-dimensional pulsed field gel electrophoresis procedures suitable for fingerprinting and mapping small eukaryotic genomes

A simple method for complete genome radiolabelling is described, involving long-wave UV exposure of agarose-embedded chromosomal DNA and [alpha-(32)P]dCTP incorporation mediated by the Klenow fragment. Experiments on the budding yeast genome show that the labelling procedure can be coupled with two new two-dimensional pulsed field gel electrophoresis (2D-PFGE) protocols of genome analysis: (i) the KARD (karyotype and restriction display)-PFGE which provides a complete view of the fragments resulting from a single restriction of the whole genome and (ii) the DDIC (double digestion of isolated chromosome)-PFGE which is the eukaryotic counterpart of complete/complete 2D-PFGE in bacterial genomics.