<title>Background</title> <p>Microsporidia are single-celled eukaryotic intracellular parasites that are related to fungi. Currently, over 1200 species have been identified, infecting animals from nearly every phylum, including commercially important species such as honeybees and fish, as well as humans [##UREF##0##1##]. Inside host cells, microsporidia proliferate as vegetative stages (meronts, schizonts) which eventually produce spores that are released when the host cell lyses. Spores possess a unique host cell invasion apparatus called the polar filament, which is forcefully everted upon germination to form a tube and can pierce a nearby host cell [##UREF##0##1##]. The tube then acts as a conduit allowing the contents of the spore to be injected into the host cell's cytoplasm, where the parasite undergoes vegetative replication.</p> <p>Microsporidia are a diverse group of organisms, and vary greatly in the complexity of their life cycles. For instance, <italic>Encephalitozoon cuniculi </italic>and <italic>Antonospora locustae </italic>produce only one type of spore (uninucleate in the former and binucleate in the latter), and complete their entire life cycles inside one host individual, while <italic>Amblyospora californica </italic>requires two host groups (mosquitoes and microcrustacea) and produces three morphologically and functionally discrete spore types [##UREF##0##1##].</p> <p>Microsporida possess some of the smallest primary nuclear genomes known (as tiny as 2.3 Mbp). The only microsporidian whose genome has been completely sequenced is the human parasite, <italic>E. cuniculi</italic>. At a meager 2.9 Mbp, <italic>E. cuniculi</italic>'s genome is extremely compact, with only 2000 genes [##REF##11719806##2##]. A small genome sequence survey (GSS) project has been conducted on <italic>A. locustae</italic>, a locust parasite that has been approved as a biological control agent in the United States [##REF##15186746##3##]. <italic>A. locustae</italic>'s genome is roughly 5.4 Mbp in size [##UREF##1##4##], or about twice the size of <italic>E. cuniculi</italic>'s genome. Despite the genome size difference, both genomes appear to be structured in much the same way. Genes are closely packed (nearly one gene per kilobase), are small in size compared to homologues in animals and fungi, and are intron-poor. There is also a much greater degree of synteny between these two organisms than would be expected given their phylogenetic relationship, which implies that although microsporidian genes are fast-evolving, genomic rearrangements occur only rarely [##REF##15186746##3##] (See Fig. ##FIG##0##1##).</p> <p>However, we have very little information on microsporidian genomes of larger sizes. <italic>Edhazardia aedis </italic>is a microsporidian that infects <italic>Aedes aegypti</italic>, the mosquito vector of the Dengue hemorrhagic and yellow fever viruses. <italic>E. aedis </italic>has been intensively studied as a viable biological control agent for <italic>A. aegypti </italic>[##UREF##2##5##] and has a genome estimated to be to be many times larger than that of <italic>E. cuniculi</italic>. There are several possible explanations for this difference: <italic>E. aedis </italic>may have more genes that control its complex life cycle. Genes may also be longer, more widely spaced, and contain more introns than <italic>E. cuniculi </italic>[##REF##11719806##2##].</p> <p>Morphological studies conducted on <italic>E. aedis </italic>have revealed at least four different types of spores – two uninucleate and two binucleate [##REF##2724178##6##,##UREF##3##7##]. The two types of uninucleate spore types differ morphologically but possess similar pyriform shapes. However, the cell division events from which they arise differ. Spores produced via mitosis are roughly 8.5 μm in length, whereas spores produced via meiosis (meiospores) are about 7.5 μm. Small binucleate spores (~6.5 μm in length) that have short polar filaments are formed first, followed by the production of larger binucleate spores (~9 μm in length) that are ovoid in shape. Meiospore formation is usually abortive and rarely produces normal spores [##REF##2724178##6##].</p> <p><italic>E. aedis' </italic>life cycle is moderately complex and involves two generations of the mosquito host. It begins when a uninucleate spore is ingested by a mosquito larva from the environment. Once in the gut, the spore germinates and begins to multiply in the host tissue. Within 48 hours, small binucleate spores are formed that are responsible for spread to other tissues. Orally infected larvae generally exhibit reduced growth, and may die before reaching maturity if the parasite load is high, thus releasing more spores into the environment. However, if the infection load is sufficiently small, the larva will mature into an adult mosquito and survive to reproduce [##REF##2724178##6##]. If the adult mosquito is female, large binucleate spores will develop in her ovaries and will infect oocytes, thus passing the infection on to the next generation where the majority of mortality occurs in larvae. Little is known about the factors that modulate the transition from one phase in the life cycle to the next, or about the changes in gene expression that occur during these transitions.</p> <p>It is also possible that the difference in genome size between <italic>E. aedis </italic>and <italic>E. cuniculi </italic>or <italic>A. locustae </italic>may have less to do with the number of genes, and more to do with genome architecture. <italic>E. aedis </italic>genes could be longer, more widely spaced, and contain more introns than <italic>E. cuniculi </italic>[##REF##11719806##2##]. In an effort to learn more about <italic>E. aedis</italic>' genome, a GSS of &gt;200 kbp was conducted [##UREF##4##8##]. This study concluded that <italic>E. aedis</italic>' genome structure is very different from those of <italic>E. cuniculi </italic>and <italic>A. locustae</italic>. A large portion of the genome is occupied by non-coding DNA and genes are not closely packed together, although the existence of local areas of compaction could not be ruled out.</p> <p>Previous examinations of ESTs from microsporidia have only been conducted on microsporidia with small genomes. These transcripts possessed unusual features that are atypical in eukaryotes. Examinations of ESTs from <italic>A. locustae </italic>[##REF##16037215##9##] and transcripts from <italic>E. cuniculi </italic>revealed numerous multi-gene transcripts. These transcripts are different from prokaryotic operons, as the proteins encoded by the transcript do not have related functions and are often not encoded on the same DNA strand. Many transcripts encode only a portion of one gene, while the other is present in its entirety [##REF##16037215##9##,##UREF##5##10##]. The reason for this phenomenon is not known, but it has been suggested that transcriptional control elements have been lost (or moved into adjacent genes) during the process of genome compaction [##UREF##5##10##].</p> <p>As <italic>E. aedis</italic>' genome and life cycle are very different from <italic>E. cuniculi </italic>and <italic>A. locustae</italic>, it is reasonable to assume that the transcript structure and number of genes present may differ as well. In this study, we describe the first survey of ESTs from a microsporidian with a much larger genome size and complex life cycle. In sequencing over 1300 transcripts, we have elucidated more of <italic>E. aedis</italic>' genome content, and have gained a profile of its transcript structure and composition. Surprisingly, the <italic>E. aedis </italic>uninucleate spore transcriptome is remarkably similar to that of <italic>A. locustae</italic>.</p>
<title>Methods</title> <p>Uninucleate <italic>E. aedis </italic>spores were grown and harvested from <italic>A. aegypti </italic>larvae as described previously [##REF##7650718##27##].</p> <p><italic>E. aedis </italic>spores were lysed in Ambion's plant RNA isolation aid and lysis/binding solution from an Ambion RNAqueous kit using a bead beater operating at 2500 rpm for 6 minutes with glass beads. RNA was extracted from the resulting supernatant using the RNAqueous kit. A microquantity cDNA library was constructed by Marligen, using the pExpress-1 vector. 1307 clones with an average insert size of 1.5 kb were uni-directionally sequenced using an automated capillary sequencer. Sequences were manually edited and analyzed using Sequencher 4.2 software. Proteins encoded by the transcripts were identified via BLASTX [##REF##9254694##28##] searches performed on the NCBI website (Genbank). Transcripts were identified as encoding a particular protein when BLASTX hits to Genbank proteins had e-values of 10<sup>-4 </sup>or lower. Transcripts were scored as "present in other microsporidia" when the best BLASTX hit was a gene present in other microsporidia or when the best hit was a gene that has a microsporidian homologue, and the homologue was identified in other microsporidia by BLASTing the <italic>E. aedis </italic>transcript against available microsporidian data. Putative <italic>E. aedis</italic>-specific genes are transcripts that contain open reading frames at least 100 base pairs in size and do not have any BLASTX hits with e-values lower than 10<sup>-3</sup>. In order to facilitate access to the EST sequences, they were uploaded and annotated by the dbEST website [##REF##17202165##17##].</p>
<title>Results</title> <title>Overview</title> <p>Sequences were deposited into the Genbank EST database and have the accession numbers <ext-link ext-link-type="gen" xlink:href="FG063843">FG063843</ext-link> to <ext-link ext-link-type="gen" xlink:href="FG065106">FG065106</ext-link>. From the 1307 clones sequenced, 133 unique genes were found; 55 were represented by a single transcript, while the remaining 78 were represented by two or more. 97 of the 133 unique genes are present in other microsporidia (See Table ##TAB##0##1##), while 10 are present in other (non-microsporidian) organisms (See Table ##TAB##1##2##), 18 are putatively <italic>E. aedis</italic>-specific and 8 have no apparent open reading frames. Coding sequences contained 43% G+C while 5' and 3' untranslated regions possessed 27% and 26%, respectively.</p> <p>Approximately a quarter of the transcripts analyzed coded for Hsp70. Almost all of the Hsp70 sequences were most similar to the "heat shock related 70 kDa protein" found in <italic>E. cuniculi </italic>(<ext-link ext-link-type="gen" xlink:href="NP_597563">NP_597563</ext-link>). Single nucleotide variation exists between sequences, usually as 3<sup>rd </sup>position synonymous substitutions. Where non-synonymous substitutions exist, they are always a single nucleotide and there are no indels between sequences. Mitochondrial-type and DNAK-like Hsp70s were also represented.</p> <p>Genes were assigned to COG categories to allow for comparison with <italic>A. locustae</italic>. Figure ##FIG##1##2## illustrates the percentages of total <italic>E. aedis </italic>transcripts that are dedicated to each COG category. Total <italic>A. locustae </italic>transcripts are provided for comparative purposes. As the randomness of the library is uncertain, it is possible that some transcripts are artificially overrepresented. It is therefore more informative to examine unique transcripts (ie. counting multiple transcripts for the same gene only once) rather than total transcripts. Figure ##FIG##2##3## displays the percentages of unique <italic>E. aedis </italic>and <italic>A. locustae </italic>transcripts dedicated to each category. Surprisingly, the values are similar and sometimes identical (maximum difference between <italic>E. aedis </italic>and <italic>A. locustae </italic>categories is 5%).</p> <p>Notable transcripts include a retrotransposon that is similar to LTR retrotransposons present in <italic>Sorghum bicolor </italic>(<ext-link ext-link-type="gen" xlink:href="AAD27571">AAD27571</ext-link>) and <italic>Nosema bombycis </italic>(<ext-link ext-link-type="gen" xlink:href="ABE26655">ABE26655</ext-link>). All belong to the Ty3/Gypsy family of retrotransposons. <italic>E. aedis </italic>also possesses a methionine aminopeptidase 2 gene (MetAP-2), which is present in <italic>E. cuniculi</italic>. There were several transcripts present that appear homologous to proteins found in various eukaryotes, but are absent in other microsporidia examined to date. These include hypothetical or unknown proteins found in <italic>Oryza</italic>, <italic>Danio </italic>and <italic>Plasmodium</italic>, as well as genes encoding proteins with identified functions, such as an adenosine kinase, a lysine-tRNA ligase and an L-asparaginase (See Table ##TAB##1##2##). In addition, <italic>E. aedis </italic>encodes a putative hydrolase-like protein that is present in <italic>A. locustae</italic>, but absent in <italic>E. cuniculi</italic>.</p> <p><italic>E. cuniculi </italic>and <italic>A. locustae </italic>both contain a small number of introns in their genomes and consequently, they have retained a minimal set of splicing machinery. These two organisms are not closely related (See Fig. ##FIG##0##1##), but they do share a few conserved introns [##UREF##6##11##]. Therefore, there is reason to suspect that some of these introns may also be present in <italic>E. aedis</italic>. Fortunately, seven transcripts of the gene encoding ribosomal protein L5 (which contains an intron in <italic>E. cuniculi</italic>) were recovered from the <italic>E. aedis </italic>library. These sequences were used to design primers to amplify the L5 gene from genomic DNA. It was found that the <italic>E. aedis </italic>L5 gene does not contain an intron.</p> <title>Transcript structure</title> <p>As <italic>E. aedis </italic>is an intracellular parasite and therefore cannot be easily cultured, RNA was limited and the library could not be constructed in a 5' cap-dependent manner. Therefore, nearly all of the inserts encoding the same gene were of different lengths, and most were 5' truncated. However, some of <italic>E. aedis</italic>' transcripts appear to have very long 5' untranslated regions (UTRs) of several hundred base pairs. To further assess transcript structure, cap-dependent 5' RACE (rapid amplification of cDNA ends) was conducted on transcripts from a moderately represented gene, glucosamine fructose-6-phosphate aminotransferase. 5' RACE confirmed that transcript lengths for this gene do vary, with 5' UTRs ranging from 255 to 348 bp (See Fig. ##FIG##3##4##).</p> <p>Contrary to the variable start sites of the transcripts, nearly all appear to have identical end sites. The notable exceptions are the heat shock related 70 kDa protein transcripts, which have somewhat variable 3' polyadenylation sites.</p> <p>There were frequently single nucleotide differences between sequences in contigs, but these differences were usually restricted to silent third position substitutions. In instances where the substitutions are not silent, they are conservative amino acid substitutions. These differences could represent different copies of the same gene or different alleles within the population (UTRs were not available in most cases to determine which).</p>
<title>Discussion</title> <title>Comparing microsporidian transcriptomes</title> <p>This is the second microsporidian EST project to be conducted and the first from a microsporidian possessing a large genome, allowing for a meaningful comparison of microsporidian spore transcriptomes. Despite the vast differences in genome size and life cycle complexity between <italic>E. aedis </italic>and <italic>A. locustae</italic>, their transcriptomes are highly similar in their compositions. The proportions of unique transcripts encoding proteins devoted to the "protein destination" COG category in both <italic>E. aedis </italic>and <italic>A. locustae </italic>are relatively large (19% and 16%, respectively) (See Fig. ##FIG##2##3##). It is interesting to note that proteomic work correlates with these results, as the number of proteins in <italic>E. cuniculi </italic>devoted to the "protein destination" COG category form a large percentage of the total proteins present (~28%) that have known functions [##REF##16691553##12##].</p> <p>When the total number of unique genes found in <italic>E. aedis </italic>and <italic>A. locustae </italic>are compared based on COG category classification, the percentages in each category are close to identical (See Fig. ##FIG##2##3##). The largest differences lie in the categories of cellular organization and biogenesis, cellular communication and signal transduction and cell rescue, defense, cell death and aging. One notable difference between the two spore transcriptomes is that no transposable elements were recovered in the <italic>A. locustae </italic>ESTs, whereas <italic>E. aedis </italic>transcribes a retrotransposon of the Ty3/gypsy family. Transposable elements have been previously reported to exist in the genomes of <italic>Nosema bombycis </italic>[##REF##16797019##13##], <italic>Spraguea lophii </italic>[##REF##9390396##14##], <italic>Brachiola algerae </italic>and <italic>E. aedis </italic>[##UREF##4##8##] (See below). To the best of our knowledge, this is the first instance of documented transposable element transcription in microsporidia, and could indicate active transposition.</p> <p>Nearly 8% of the unique transcripts from <italic>E. aedis </italic>encode genes that are present in various eukaryotes, but are absent from other microsporidia. The existence of these genes has several possible explanations. Sequence data from microsporidia is scarce, and the only completely sequenced genome is that of <italic>E. cuniculi</italic>. Therefore, it is currently impossible to assert that these genes are absent in any microsporidia other than <italic>E. cuniculi</italic>. The possibility exists that they were present in the genome of the microsporidian ancestor, and were lost during genome reduction/compaction events in <italic>E. cuniculi</italic>. These genes could also have arisen from lateral transfer events or they could have come to resemble genes in other organisms by chance or by convergence. Parsimoniously, the first explanation seems most likely, therefore, these data seem to suggest that the ancestor of microsporidia was not, indeed, compact to the extent of <italic>E. cuniculi</italic>.</p> <p>The MetAP-2 protein is a target for drug therapy in <italic>E. cuniculi </italic>[##UREF##7##15##]. The <italic>E. aedis </italic>copy of the MetAP-2 gene is very similar to that present in <italic>E. cuniculi</italic>, and contains the amino acid residues that bind the drug fumagillin as well as those believed to coordinate metals. Like <italic>E. cuniculi</italic>, <italic>E. aedis </italic>lacks a polylysine tract at the N-terminus of the MetAP-2 protein that is present in animals, other fungi and plants. This tract plays a role in hindering the phosphorylation of eukaryotic initiation factor 2α (eIF2α), and its absence indicates that the microsporidian proteins likely lack this function [##UREF##7##15##].</p> <p>Although our work indicates that the <italic>E. aedis </italic>L5 gene does not contain an intron like its <italic>E. cuniculi </italic>homologue (see Results, above), there is reason to believe that there are introns elsewhere in the genome. There are several transcripts encoding proteins that act in pre-mRNA splicing: an arginine/serine rich pre-mRNA splicing factor (<ext-link ext-link-type="gen" xlink:href="NP_597487">NP_597487</ext-link> in <italic>E. cuniculi</italic>), a pre-mRNA splicing factor (<ext-link ext-link-type="gen" xlink:href="NP_586183">NP_586183</ext-link> in <italic>E. cuniculi</italic>) and a U5 associated snRNP (<ext-link ext-link-type="gen" xlink:href="NP_586393">NP_586393</ext-link> in <italic>E. cuniculi</italic>). These genes comprise 2.2% of the total unique genes found.</p> <title>Hsp70</title> <p>Roughly 28% of total <italic>E. aedis </italic>transcripts encoded some form of Hsp70, a heat shock protein that assists in the folding of other proteins. Hsp70 helps prevent proteins from becoming insoluble and also plays a role in various other intracellular processes, such as apoptosis [##REF##15770419##16##]. The action of Hsp70 allows mutant proteins to continue functioning by being refolded instead of being degraded, which necessitates the costly synthesis of more protein. The number of Hsp70 transcripts in the <italic>E. aedis </italic>ESTs is an order of magnitude higher than was found in <italic>A. locustae </italic>(2%) [##REF##17202165##17##]. We are cautious in this interpretation as we have not quantitatively assessed the transcription level of Hsp70 in <italic>E. aedis</italic>, and it is likely that transcripts of this protein are somewhat overrepresented in the library.</p> <p>Although no <italic>E. cuniculi </italic>ESTs have been published, Brosson et al. [##REF##16691553##12##] investigated the proteins present in spores. Hsp70 constitutes a moderate amount of all protein present. Brosson and his colleagues classified all proteins based on their COG categories, and found that all "protein destination" proteins together comprise 21% of <italic>E. cuniculi</italic>'s proteome. Intriguingly, Brosson et al.'s [##REF##16691553##12##] experiments indicate that of the four copies of Hsp70 in <italic>E. cuniculi</italic>, the predominately expressed copy of Hsp70 in <italic>E. cuniculi </italic>is homologous to the highly represented transcript in <italic>E. aedis</italic>. In <italic>A. locustae</italic>, the most highly transcribed copy was most similar to the abundantly transcribed copy in <italic>E. aedis </italic>as well [##REF##17202165##17##]. Therefore, it is likely that microsporidia employ similar primary mechanisms to ensure proper folding of proteins.</p> <p>In other parasites and endosymbionts, such as <italic>Buchnera aphidicola</italic>, Hsp70 is also highly expressed [##REF##12791133##18##] and may constitute up to 10% of the protein contained in the cell at any one time. In species that lead parasitic or endosymbiotic lifestyles, genetic drift and relaxed selection pressure frequently lead to an increased mutation rate. The need for Hsp70 in order for proteins to fold correctly seems to increase with both the size and number of mutations in the protein [##REF##15770419##16##]. Although microsporidian genomes appear to have had little rearrangement, the nucleotide mutation rate seems to be high in this group of organisms [##REF##15599510##19##,##REF##10767522##20##]. Microsporidia could, therefore, contain elevated levels of Hsp70 in order to allow folding of mutant proteins.</p> <title>Transposable elements</title> <p>One of the <italic>E. aedis </italic>ESTs closely matches the integrase domain of the Ty3/gypsy family of retrotransposons. Several of these elements were identified in a GSS of <italic>E. aedis </italic>[##UREF##4##8##] and a few other microsporidian species, but to the best of our knowledge, this is the first instance in which transcripts of any microsporidian retrotransposon have been found. Transcripts could be indicative of active transposition occurring in <italic>E. aedis</italic>' genome.</p> <p>Ty3/gypsy retrotransposons exist in many organisms ranging from the microsporidia <italic>Spraguea lophii </italic>[##REF##9390396##14##], <italic>Brachiola algerae </italic>[##UREF##4##8##], and <italic>Nosema bombycis </italic>[##REF##16797019##13##] to <italic>Saccharomyces</italic>, <italic>Drosophila </italic>and <italic>Sorghum</italic>. Ty3 elements have been well characterized in budding yeast, and exist in 1–4 copies per genome, where they are transcribed by RNA polymerase III. Transcription typically occurs only in haploid cells in the presence of mating pheromones [##REF##7705653##21##]. The <italic>N. bombycis </italic>genome contains at least 8 different retrotransposons in the Ty3/gypsy family, but unlike yeast, they are not exclusively located upstream of tRNAs [##REF##16797019##13##]. Nearly all <italic>N. bombycis </italic>retrotransposons encode a polyprotein containing 5 domains, which exist in a defined order: Gag, protease, reverse transcriptase, RnaseH and integrase. As many of the sequences in the <italic>E. aedis </italic>library appear to be 5' truncated, it is possible that the other domains upstream of the integrase in the polyprotein are also present in genomic DNA. Indeed, the GSS project revealed sequences matching the reverse transcriptase domain [##UREF##4##8##]. Although the microsporidian <italic>Vittaforma corneae </italic>is also known to possess at least one transposable element [##REF##12425527##22##], it belongs to a different family than those present in <italic>E. aedis </italic>- the L1 family present in humans.</p> <p>The only completely sequenced microsporidian genome, that of <italic>E. cuniculi </italic>[##REF##11719806##2##], is completely devoid of transposable elements. The existence of similar transposable elements (of the Ty3/gypsy family) in the distantly related <italic>S. lophii, N. bombycis</italic>, <italic>B. algerae </italic>and <italic>E. aedis </italic>(See Fig. ##FIG##0##1##) implies that this element may have been present in the genome of the ancestor of microsporidia. Therefore, the process of genome compaction that gave rise to the <italic>E. cuniculi </italic>genome likely involved purging transposable elements.</p> <p>It has been suggested that transposable elements may act to reorganize genes within the genome. Xu et al. [##REF##16797019##13##] compared regions of synteny between <italic>N. bombycis </italic>and <italic>E. cuniculi </italic>chromosomes, as selection appears to be acting to retain gene synteny among microsporidia, even if they are only distantly related [##REF##15186746##3##]. In <italic>N. bombycis</italic>, transposable elements flank these syntenic regions [##REF##16797019##13##]. If <italic>E. aedis' </italic>large genome is partially a product of transposable element proliferation, one would expect much less synteny between this species and other microsporidia. Perhaps future research will elucidate other roles that transposable elements have played in shaping microsporidian genomes, especially since the minute genome of <italic>E. cuniculi </italic>seems to lack them, while they are present in larger genomes.</p> <p>The functions that these transposable elements perform in a given genome are cryptic at best, but evidence is emerging that they may be more than just simply parasitic DNA. Peaston et al. [##REF##15469847##23##] recently discovered that a class of mouse retrotransposons appears to regulate gene expression in embryos.</p> <title>Transcript structure</title> <p>Transcripts in <italic>A. locustae </italic>typically contain more than one gene. These transcripts do not necessarily contain complete open reading frames for all genes and the genes are frequently in opposite orientations [##REF##16037215##9##]. It is not known how many proteins are made from each transcript or whether this situation is typical for microsporidia, but recent work by Corradi et al. [##UREF##5##10##] suggests that <italic>E. cuniculi </italic>also possesses multi-gene transcripts.</p> <p>Unlike <italic>A. locustae </italic>and <italic>E. cuniculi</italic>, <italic>E. aedis </italic>appears to transcribe very few multi-gene transcripts, if any at all. This is not unexpected, given that <italic>E. aedis </italic>genes appear to be separated by large intergenic spaces [##UREF##4##8##]. The <italic>E. aedis </italic>GSS could not rule out the possibility that local areas of compacted genes might exist [##UREF##4##8##]. Given the lack of multi-gene transcripts identified, this seems increasingly unlikely.</p> <p>Also contrary to what is found in <italic>A. locustae</italic>, nearly all of <italic>E. aedis</italic>' transcripts encode proteins in a positive frame (&lt;1% are in a negative frame, compared to 17% in <italic>A. locustae</italic>) [##REF##16037215##9##]. Although antisense transcripts are used in many organisms (possibly also <italic>A. locustae</italic>) to suppress translation, it appears unlikely that this type of regulation occurs in <italic>E. aedis</italic>. Conversely, the large number of antisense transcripts in <italic>A. locustae </italic>may be due to a lack of transcriptional regulation resulting from genome compaction.</p> <p><italic>E. aedis</italic>' transcripts seem to start at multiple locations upstream of the start codon (5' UTR length is 180 bp on average) but terminate at the same position with a relatively short 3' UTR (51 bp on average) (See, for example, Figure ##FIG##3##4##). This is more in line with transcription in <italic>E. cuniculi </italic>and contrasts with the situation in <italic>A. locustae</italic>, where transcripts start directly upstream of the translation initiation site, but often terminate much farther downstream in the adjacent gene [##UREF##5##10##]. For comparison, the yeast <italic>S. cerevisiae </italic>contains much shorter 5' UTRs than 3' UTRs (15–75 and ~144 bp, respectively [##REF##15905473##24##,##REF##9889288##25##]), a common trend seen in other fungi, plants and animals. The reason for this reversal is unknown, since 3' UTRs are ubiquitously used as translation regulators. It is likely that <italic>E. aedis </italic>lacks some of the translational control mechanisms present in other fungi, plants and animals [##REF##12575997##26##].</p>
<title>Conclusion</title> <p>This is the first examination of ESTs from a microsporidian containing a large genome. The extent of genome compaction in the microsporidian ancestor is not known, but the presence of genes in <italic>E. aedis </italic>that have not been found in other microsporidia suggests that extreme reduction and compaction occurred only in specific lineages. Surprisingly, <italic>E. aedis </italic>has a predicted uninucleate spore transcriptome that is highly similar to that of the distantly related <italic>A. locustae</italic>, although the two species have diverse life cycles and genome sizes.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Microsporidia are a group of parasites related to fungi that infect a wide variety of animals and have gained recognition from the medical community in the past 20 years due to their ability to infect immuno-compromised humans. Microsporidian genomes range in size from 2.3 to 19.5 Mbp, but almost all of our knowledge comes from species that have small genomes (primarily from the human parasite <italic>Encephalitozoon cuniculi </italic>and the locust parasite <italic>Antonospora locustae</italic>). We have conducted an EST survey of the mosquito parasite <italic>Edhazardia aedis</italic>, which has an estimated genome size several times that of more well-studied species. The only other microsporidian EST project is from <italic>A. locustae</italic>, and serves as a basis for comparison with <italic>E. aedis</italic>.</p> <title>Results</title> <p>The spore transcriptomes of <italic>A. locustae </italic>and <italic>E. aedis </italic>were compared and the numbers of unique transcripts that belong to each COG (Clusters of Orthologous Groups of proteins) category differ by at most 5%. The transcripts themselves have widely varying start sites and encode a number of proteins that have not been found in other microsporidia examined to date. However, <italic>E. aedis </italic>seems to lack the multi-gene transcripts present in <italic>A. locustae </italic>and <italic>E. cuniculi</italic>. We also present the first documented case of transcription of a transposable element in microsporidia.</p> <title>Conclusion</title> <p>Although <italic>E. aedis </italic>and <italic>A. locustae </italic>are distantly related, have very disparate life cycles and contain genomes estimated to be vastly different sizes, their patterns of transcription are similar. The architecture of the ancestral microsporidian genome is unknown, but the presence of genes in <italic>E. aedis </italic>that have not been found in other microsporidia suggests that extreme genome reduction and compaction is lineage specific and not typical of all microsporidia.</p>
<title>Authors' contributions</title> <p>EEG extracted RNA from the <italic>E. aedis </italic>spores, performed and interpreted the sequence analyses and drafted the manuscript. JJB cultivated insect larvae and harvested <italic>E. aedis </italic>spores. NMF conceived of this study, contributed to the interpretation of the results and helped draft the manuscript.</p>
<title>Acknowledgements</title> <p>EEG's work is supported by an NSERC Postgraduate Scholarship and research in the Fast lab is funded by an NSERC Discovery Grant to NMF.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>The phylogenetic relationships between several microsporidia.</bold> Species that house transposable elements belonging to the Ty3/gypsy family are highlighted in blue, while species containing LTR transposons are highlighted in yellow. Genome sizes are indicated to the right of each species. (Adapted from Slamovits et al., 2004.)</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Total <italic>E. aedis </italic>transcripts represented by COG (Clusters of Orthologous Groups of proteins) category with and without Hsp70.</bold> Total <italic>A. locustae </italic>transcripts are provided for comparison. (<italic>A. locustae </italic>data adapted from Williams et al., 2005.)</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Unique <italic>E. aedis </italic>transcripts represented by COG category.</bold> Unique <italic>A. locustae </italic>transcripts are provided for comparison. (<italic>A. locustae </italic>data adapted from Williams et al., 2005.)</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>5' RACE conducted on a moderately represented transcript in <italic>E. aedis </italic>reveals multiple transcription start sites.</bold> ESTs are depicted in green and RACE products in purple. The predicted translational start codon is indicated by the orange arrow. As indicated, the <italic>E. cuniculi </italic>homologue of this gene contains 128 amino acids at N-terminus that appear to be absent in <italic>E. aedis</italic>.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Gene Name</bold></td><td align="left"><bold>Species</bold></td><td align="left"><bold>Genbank Accession Number</bold></td></tr></thead><tbody><tr><td align="left">16S rRNA GENE</td><td align="left"><italic>Brachiola algerae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AM422905">AM422905</ext-link></td></tr><tr><td align="left">1-ACYL-SN-GLYCEROL-3-PHOSPHATE ACYLTRANSFERASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586146">NP_586146</ext-link></td></tr><tr><td align="left">26S PROTEASOME REGULATORY SUBUNIT 4</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586091">NP_586091</ext-link></td></tr><tr><td align="left">26S PROTEASOME REGULATORY SUBUNIT 6</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586128">NP_586128</ext-link></td></tr><tr><td align="left">26S PROTEASOME REGULATORY SUBUNIT 8</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_955738">XP_955738</ext-link></td></tr><tr><td align="left">40S RIBOSOMAL PROTEIN S2</td><td align="left"><italic>Leishmania infantum</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_001466537">XP_001466537</ext-link></td></tr><tr><td align="left">40S RIBOSOMAL PROTEIN S3</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_955676">XP_955676</ext-link></td></tr><tr><td align="left">40S RIBOSOMAL PROTEIN S4</td><td align="left"><italic>Mycetophagus quadripustulatus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="CAJ17168">CAJ17168</ext-link></td></tr><tr><td align="left">40S RIBOSOMAL PROTEIN SA or P40</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_584728">NP_584728</ext-link></td></tr><tr><td align="left">60S RIBOSOMAL PROTEIN L3</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597630">NP_597630</ext-link></td></tr><tr><td align="left">60S RIBOSOMAL PROTEIN L4</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597213">NP_597213</ext-link></td></tr><tr><td align="left">60S RIBOSOMAL PROTEIN L5</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_585846">NP_585846</ext-link></td></tr><tr><td align="left">6-PHOSPHOFRUCTOKINASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597579">NP_597579</ext-link></td></tr><tr><td align="left">ABC TRANSPORTER (MITOCHONDRIAL TYPE) <bold>#1</bold></td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586426">NP_586426</ext-link></td></tr><tr><td align="left">ABC TRANSPORTER (MITOCHONDRIAL TYPE) <bold>#2</bold></td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586426">NP_586426</ext-link></td></tr><tr><td align="left">ACTIN</td><td align="left"><italic>Blakeslea trispora</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AAW32475">AAW32475</ext-link></td></tr><tr><td align="left">ARGININE/SERINE RICH PRE-mRNA SPLICING FACTOR</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597487">NP_597487</ext-link></td></tr><tr><td align="left">ASSOCIATED WITH RAN (NUCLEAR IMPORT/EXPORT) FUNCTION FAMILY MEMBER</td><td align="left"><italic>Caenorhabditis elegans</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_499369">NP_499369</ext-link></td></tr><tr><td align="left">ATP SYNTHASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_955732">XP_955732</ext-link></td></tr><tr><td align="left">BELONGS TO THE ABC TRANSPORTER SUPERFAMILY</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597462">NP_597462</ext-link></td></tr><tr><td align="left">cAMP-DEPENDENT PROTEIN KINASE TYPE 1 REGULATORY CHAIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597223">NP_597223</ext-link></td></tr><tr><td align="left">CASEIN KINASE 1 HOMOLOG (INVOLVED IN DNA REPAIR</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597600">NP_597600</ext-link></td></tr><tr><td align="left">CATION-TRANSPORTING ATPase</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586078">NP_586078</ext-link></td></tr><tr><td align="left">CHOLINE PHOSPHATE CYTIDYLYLTRANSFERASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586276">NP_586276</ext-link></td></tr><tr><td align="left">DNA REPLICATION LICENSING FACTOR MCM2</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_584768">NP_584768</ext-link></td></tr><tr><td align="left">DNA REPLICATION LICENSING FACTOR OF THE MCM FAMILY MCM6</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597420">NP_597420</ext-link></td></tr><tr><td align="left">DNA REPLICATION LICENSING FACTOR OF THE MCM FAMILY MCM7</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_585977">NP_585977</ext-link></td></tr><tr><td align="left">DNAJ PROTEIN HOMOLOG 2</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586004">NP_586004</ext-link></td></tr><tr><td align="left">DNAK-LIKE PROTEIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586489">NP_586489</ext-link></td></tr><tr><td align="left">EUKARYOTIC TRANSLATION INITIATION FACTOR 4A</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_955671">XP_955671</ext-link></td></tr><tr><td align="left">FIBRILLARIN (34kDa NUCLEOLAR PROTEIN)</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586197">NP_586197</ext-link></td></tr><tr><td align="left">GENERAL TRANSCRIPTION FACTOR</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597292">NP_597292</ext-link></td></tr><tr><td align="left">GLUCOSAMINE FRUCTOSE-6-PHOSPHATE AMINOTRANSFERASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586057">NP_586057</ext-link></td></tr><tr><td align="left">GLYCERALDEHYDE-3-PHOSPHATE DEHYDROGENASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586008">NP_586008</ext-link></td></tr><tr><td align="left">GUANINE NUCLEOTIDE BINDING PROTEIN BETA SUBUNIT</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597241">NP_597241</ext-link></td></tr><tr><td align="left">HEAT SHOCK RELATED 70 kDa PROTEIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597563">NP_597563</ext-link></td></tr><tr><td align="left">HEAT-SHOCK PROTEIN HSP90 HOMOLOG</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_584635">NP_584635</ext-link></td></tr><tr><td align="left">HISTIDYL tRNA SYNTHETASE</td><td align="left"><italic>Antonospora locustae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AAT12372">AAT12372</ext-link></td></tr><tr><td align="left">HISTONE ACETYLTRANSFERASE TYPE B SUBUNIT 2</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586003">NP_586003</ext-link></td></tr><tr><td align="left">HISTONE DEACETYLASE 1</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597645">NP_597645</ext-link></td></tr><tr><td align="left">HISTONE DEACETYLASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_955621">XP_955621</ext-link></td></tr><tr><td align="left">HISTONE H3</td><td align="left"><italic>Mus musculus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="JQ1983">JQ1983</ext-link></td></tr><tr><td align="left">HSP 101 RELATED PROTEIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586448">NP_586448</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN ECU02_0840</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_584609">NP_584609</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN ECU02_0950</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_584620">NP_584620</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN ECU06_0450</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_585801">NP_585801</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN ECU06_1280</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_585884">NP_585884</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN ECU07_0530</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_585981">NP_585981</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN ECU08_1500</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597278">NP_597278</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN ECU09_0740</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_955628">XP_955628</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN ECU09_1700</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_955723">XP_955723</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN ECU11_1720</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586478">NP_586478</ext-link></td></tr><tr><td align="left">LIM DOMAIN-CONTAINING PROTEIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586340">NP_586340</ext-link></td></tr><tr><td align="left">LONG CHAIN FATTY ACID CoA LIGASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586206">NP_586206</ext-link></td></tr><tr><td align="left">METHIONINE AMINOPEPTIDASE TYPE 2</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586190">NP_586190</ext-link></td></tr><tr><td align="left">METHIONINE PERMEASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_585905">NP_585905</ext-link></td></tr><tr><td align="left">NIFS-LIKE PROTEIN (CYSTEINE DESULFURASE) INVOLVED IN IRON-SULFUR CLUSTER SYNTHESIS</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586483">NP_586483</ext-link></td></tr><tr><td align="left">NUCLEAR SER/THR PROTEIN PHOSPHATASE PP1-1 GAMMA CATALYTIC SUBUNIT</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597385">NP_597385</ext-link></td></tr><tr><td align="left">P68-LIKE PROTEIN (DEAD BOX FAMILY OF RNA HELICASES)</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597238">NP_597238</ext-link></td></tr><tr><td align="left">PEPTIDE CHAIN RELEASE FACTOR SUBUNIT 1</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597376">NP_597376</ext-link></td></tr><tr><td align="left">PEPTIDE ELONGATION FACTOR 2</td><td align="left"><italic>Glugea plecoglossi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="BAA11470">BAA11470</ext-link></td></tr><tr><td align="left">PHOSPHATIDYLINOSITOL TRANSFER PROTEIN, ALPHA</td><td align="left"><italic>Danio rerio</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_957229">NP_957229</ext-link></td></tr><tr><td align="left">PHOSPHOMANNOMUTASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597365">NP_597365</ext-link></td></tr><tr><td align="left">POLYADENYLATE-BINDING PROTEIN 2</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586226">NP_586226</ext-link></td></tr><tr><td align="left">POLYPROTEIN</td><td align="left"><italic>Sorghum bicolor</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AAD27571">AAD27571</ext-link></td></tr><tr><td align="left">PRE-mRNA SPLICING FACTOR</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586183">NP_586183</ext-link></td></tr><tr><td align="left">PROTEIN KINASE B-LIKE PROTEIN</td><td align="left"><italic>Plasmodium falciparum</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AAT06260">AAT06260</ext-link></td></tr><tr><td align="left">PROTEIN TRANSPORT PROTEIN SEC23 HOMOLOG (COPII COAT)</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586385">NP_586385</ext-link></td></tr><tr><td align="left">PUTATIVE HYDROLASE-LIKE PROTEIN</td><td align="left"><italic>Antonospora locustae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AAU11090">AAU11090</ext-link></td></tr><tr><td align="left">PUTATIVE ZINC FINGER PROTEIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597297">NP_597297</ext-link></td></tr><tr><td align="left">SER/THR PROTEIN PHOSPHATASE 2-A</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_584753">NP_584753</ext-link></td></tr><tr><td align="left">SER/THR PROTEIN PHOSPHATASE PP2-A REGULATORY SUBUNIT B</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597423">NP_597423</ext-link></td></tr><tr><td align="left">SERINE/THREONINE PROTEIN KINASE (REQUIRED FOR ACTIN RING AND SEPTATION)</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_965898">XP_965898</ext-link></td></tr><tr><td align="left">SIMILAR TO DNAJ-LIKE PROTEIN</td><td align="left"><italic>Nasonia vitripennis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_001602403">XP_001602403</ext-link></td></tr><tr><td align="left">SIMILARITY TO 14-3-3 PROTEIN 1</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597610">NP_597610</ext-link></td></tr><tr><td align="left">SIMILARITY TO ADP/ATP CARRIER PROTEIN</td><td align="left"><italic>Paranosema grylli</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="CAI30461">CAI30461</ext-link></td></tr><tr><td align="left">SIMILARITY TO CDC20 (WD-REPEAT PROTEIN)</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597660">NP_597660</ext-link></td></tr><tr><td align="left">SIMILARITY TO Hsp70-RELATED PROTEIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_584537">NP_584537</ext-link></td></tr><tr><td align="left">SIMILARITY TO HYPOTHETICAL INTEGRAL MEMBRANE PROTEIN YQ55_CAEEL</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597662">NP_597662</ext-link></td></tr><tr><td align="left">SIMILARITY TO HYPOTHETICAL PROTEIN YAAT_BACSU</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597532">NP_597532</ext-link></td></tr><tr><td align="left">SIMILARITY TO HYPOTHETICAL PROTEIN YB36_METJA</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597239">NP_597239</ext-link></td></tr><tr><td align="left">SIMILARITY TO PUTATIVE AMINOACID TRANSPORTER YEU9_yeast</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_584803">NP_584803</ext-link></td></tr><tr><td align="left">SIMILARITY TO SKT5 PROTEIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586349">NP_586349</ext-link></td></tr><tr><td align="left">SIMILARITY TO TRANSCRIPTION INITIATION FACTOR TFIIA</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597616">NP_597616</ext-link></td></tr><tr><td align="left">STE12 TRANSCRIPTION FACTOR</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586509">NP_586509</ext-link></td></tr><tr><td align="left">STRUCTURE-SPECIFIC RECOGNITION PROTEIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586030">NP_586030</ext-link></td></tr><tr><td align="left">T COMPLEX PROTEIN 1 SUBUNIT BETA</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_955601">XP_955601</ext-link></td></tr><tr><td align="left">THREONYL tRNA SYNTHETASE <bold>#1</bold></td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586084">NP_586084</ext-link></td></tr><tr><td align="left">THREONYL tRNA SYNTHETASE <bold>#2</bold></td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586084">NP_586084</ext-link></td></tr><tr><td align="left">TRANSLATION ELONGATION FACTOR 1 ALPHA</td><td align="left"><italic>Glugea plecoglossi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="BAA12288">BAA12288</ext-link></td></tr><tr><td align="left">TRIOSE PHOSPHATE ISOMERASE</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586329">NP_586329</ext-link></td></tr><tr><td align="left">TUBULIN BETA CHAIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_597591">NP_597591</ext-link></td></tr><tr><td align="left">U5 ASSOCIATED snRNP</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586393">NP_586393</ext-link></td></tr><tr><td align="left">UNNAMED PROTEIN PRODUCT (Hsp70)</td><td align="left"><italic>Candida glabrata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_445544">XP_445544</ext-link></td></tr><tr><td align="left">VACUOLAR ATP SYNTHASE CATALYTIC SUBUNIT A</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586434">NP_586434</ext-link></td></tr><tr><td align="left">VACUOLAR ATP SYNTHASE SUBUNIT B</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_586219">NP_586219</ext-link></td></tr><tr><td align="left">ZINC FINGER PROTEIN</td><td align="left"><italic>Encephalitozoon cuniculi</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_584833">NP_584833</ext-link></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Gene Name</bold></td><td align="left"><bold>Species Name</bold></td><td align="left"><bold>Genbank Accession Number</bold></td></tr></thead><tbody><tr><td align="left">60S RIBOSOMAL PROTEIN L2</td><td align="left"><italic>Babesia bovis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_001612300">XP_001612300</ext-link></td></tr><tr><td align="left">ADENOSINE KINASE</td><td align="left"><italic>Homo sapiens</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AAA97893">AAA97893</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN</td><td align="left"><italic>Candida albicans</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_717148">XP_717148</ext-link></td></tr><tr><td align="left">HYPOTHETICAL PROTEIN PY5484</td><td align="left"><italic>Plasmodium yoelii yoelii</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_725949">XP_725949</ext-link></td></tr><tr><td align="left">L-ASPARIGINASE</td><td align="left"><italic>Dirofilaria immitis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="Q9U518">Q9U518</ext-link></td></tr><tr><td align="left">LYSINE tRNA LIGASE</td><td align="left"><italic>Saccharomyces cerevisiae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="CAA39699">CAA39699</ext-link></td></tr><tr><td align="left">PUTATIVE VESICULAR TRANSPORT FACTOR USO1P</td><td align="left"><italic>Candida albicans</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_710120">XP_710120</ext-link></td></tr><tr><td align="left">SEC63 DOMAIN CONTAINING PROTEIN</td><td align="left"><italic>Trichomonas vaginalis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_001580151">XP_001580151</ext-link></td></tr><tr><td align="left">PROTEIN PHOSPHATASE 2B</td><td align="left"><italic>Cryptosporidium hominis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="XP_666159">XP_666159</ext-link></td></tr><tr><td align="left">WD-40 REPEAT FAMILY PROTEIN</td><td align="left"><italic>Arabidopsis thaliana</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NP_201533">NP_201533</ext-link></td></tr></tbody></table></table-wrap>
<table-wrap-foot><p>The 97 unique <italic>E. aedis </italic>transcripts which are homologous to genes present in other microsporidia are listed above. Species names and Genbank accession numbers of top BLASTX hits are indicated. Bold text in the "Gene Name" column indicates instances where two different transcripts both had the same top BLASTX hit. Underlining indicates a copy of Hsp70 that is most similar to a protein that remains unnamed in Genbank.</p></table-wrap-foot> <table-wrap-foot><p>The 10 genes which are present in <italic>E. aedis </italic>but absent from other microsporidia are listed above. The species names and Genbank accession numbers of the top BLASTX hit for each gene are also listed.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-296-1"/> <graphic xlink:href="1471-2164-9-296-2"/> <graphic xlink:href="1471-2164-9-296-3"/> <graphic xlink:href="1471-2164-9-296-4"/>
[{"surname": ["Wittner", "Weiss"], "given-names": ["M", "LM"], "source": ["The Microsporidia and Microsporidosis"], "year": ["1999"], "publisher-name": ["Washington, D. C.: ASM Press"]}, {"surname": ["Streett"], "given-names": ["DA"], "article-title": ["Analysis of "], "italic": ["Nosema locustae "], "source": ["J Invertebr Pathol"], "year": ["1994"], "volume": ["63"], "fpage": ["301"], "lpage": ["303"], "pub-id": ["10.1006/jipa.1994.1055"]}, {"surname": ["Becnel"], "given-names": ["JJ"], "italic": ["Edhazardia aedis ", "Aedes aegypti "], "source": ["Proc Vth Int Colloq Invertebr Pathol Microb Control: Adelaide, Australia"], "year": ["1990"], "fpage": ["56"], "lpage": ["60"]}, {"surname": ["Johnson", "Becnel", "Undeen"], "given-names": ["MA", "JJ", "AH"], "article-title": ["A new sporulation sequence in "], "italic": ["Edhazardia aedis ", "Aedes aegypti "], "source": ["J Invertebr Pathol"], "year": ["1997"], "volume": ["70"], "fpage": ["69"], "lpage": ["75"], "pub-id": ["10.1006/jipa.1997.4678"]}, {"surname": ["Williams", "Lee", "Becnel", "Weiss", "Fast", "Keeling"], "given-names": ["BAP", "RCH", "JJ", "LM", "NM", "PJ"], "article-title": ["Genome sequence surveys of "], "italic": ["Brachiola algerae ", "Edhazardia aedis "], "source": ["BMC Genomics"]}, {"surname": ["Corradi", "Gangaeva", "Keeling"], "given-names": ["N", "A", "PJ"], "article-title": ["Comparative profiling of overlapping transcription in the compacted genomes of microsporidia "], "italic": ["Antonospora locustae ", "Encephalitozoon cuniculi"], "source": ["Genomics"]}, {"surname": ["Limpright", "Fast"], "given-names": ["VO", "NM"], "article-title": ["personal communication"]}, {"surname": ["Pandrea", "Mittleider", "Brindley", "Didier", "Robertson"], "given-names": ["I", "D", "PJ", "ES", "DL"], "article-title": ["Phylogenetic relationships of methionine aminopeptidase 2 among "], "italic": ["Encephalitozoon "], "source": ["Mol Biochem Parasit"], "year": ["2005"], "volume": ["140"], "fpage": ["141"], "lpage": ["152"], "pub-id": ["10.1016/j.molbiopara.2004.12.006"]}]
BMC Genomics. 2008 Jun 20; 9:296