Ancient and novel small RNA pathways compensate for the loss of piRNAs in multiple independent nematode lineages.

Small RNA pathways act at the front line of defence against transposable elements across the Eukaryota. In animals, Piwi interacting small RNAs (piRNAs) are a crucial arm of this defence. However, the evolutionary relationships among piRNAs and other small RNA pathways targeting transposable elements are poorly resolved. To address this question we sequenced small RNAs from multiple, diverse nematode species, producing the first phylum-wide analysis of how small RNA pathways evolve. Surprisingly, despite their prominence in Caenorhabditis elegans and closely related nematodes, piRNAs are absent in all other nematode lineages. We found that there are at least two evolutionarily distinct mechanisms that compensate for the absence of piRNAs, both involving RNA-dependent RNA polymerases (RdRPs). Whilst one pathway is unique to nematodes, the second involves Dicer-dependent RNA-directed DNA methylation, hitherto unknown in animals, and bears striking similarity to transposon-control mechanisms in fungi and plants. Our results highlight the rapid, context-dependent evolution of small RNA pathways and suggest piRNAs in animals may have replaced an ancient eukaryotic RNA-dependent RNA polymerase pathway to control transposable elements.We thank Sylviane Moss for high-throughput sequencing support. We thank Charles Bradshaw for help with computation and IT. We thank Marie-Anne Felix and Frank Jiggins for critical comments on the manuscript. We thank Matt Berriman (Wellcome Trust Sanger Centre, Hinxton, Cambridge, UK) for allowing us to use unpublished genomic sequencing data for N. brasiliensis. We thank Einhardt Schierenberg (University of Cologne, Germany) and Werner Armonies (Alfred Wegener Institute, Sylt, Germany) for help with collection of E. brevis.This is the final version of the article, originally published in PLoS Biology, 2015, 13(2): e1002061. doi:10.1371/journal.pbio.100206

Alpha and beta subunits of pyruvate dehydrogenase E1 from the microsporidian Nosema locustae: mitochondrionderived carbon metabolism in microsporidia.

Abstract Microsporidia are highly adapted eukaryotic intracellular parasites that infect a variety of animals. Microsporidia contain no recognisable mitochondrion, but recently have been shown to have evolved from fungi and to possess heat shock protein genes derived from mitochondria. These findings make it clear that microsporidian ancestors were mitochondrial, yet it remains unknown whether they still contain the organelle, and if so what its role in microsporidian metabolism might be. Here we have characterised genes encoding the alpha and beta subunits of pyruvate dehydrogenase complex E1 (PDH, EC 1.2.4.1) from the microsporidian Nosema locustae. All other amitochondriate eukaryotes studied to date have lost the PDH complex and replaced it with pyruvate:ferredoxin oxidoreductase (PFOR). Nevertheless, molecular phylogeny shows that these Nosema enzymes are most closely related to mitochondrial PDH from other eukaryotes, demonstrating that elements of mitochondrial metabolism have been retained in microsporidia, and that PDH has not been wholly lost. However, there is still no evidence for a mitochondrion in microsporidia, and neither PDH subunit is predicted to encode an amino terminal leader sequence that could function as a mitochondrion-targeting transit peptide, raising questions as to whether these proteins function in a relic organelle or in the cytosol. Moreover, it is also unclear whether these proteins remain part of the PDH complex, or whether they have been retained for another purpose. We propose that microsporidia may utilise a unique pyruvate decarboxylation pathway involving PDH, demonstrating once again the diversity of core metabolism in amitochondriate eukaryotes

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.Botany, Department ofScience, Faculty ofReviewedFacult

Evolutionary relationship between translation initiation factor eIF-2gamma and selenocysteine-specific elongation factor SELB: change of function in translation factors

Abstract. Eubacterial and eukaryotic translation initiation systems have very little in common, and therefore the evolutionary events that gave rise to these two disparate systems are difficult to ascertain. One common feature is the presence of initiation, elongation, and release factors belonging to a large GTPase superfamily. One of these initiation factors, the ␥ subunit of initiation factor 2 (eIF-2␥), is found only in eukaryotes and archaebacteria. We have sequenced eIF-2␥ gene fragments from representative diplomonads, parabasalia, and microsporidia and used these new sequences together with new archaebacterial homologues to examine the phylogenetic position of eIF-2␥ within the GTPase superfamily. The archaebacterial and eukaryotic eIF-2␥ proteins are found to be very closely related, and are in turn related to SELB, the selenocysteine-specific elongation factor from eubacteria. The overall topology of the GTPase tree further suggests that the eIF-2␥/SELB group may represent an ancient subfamily of GTPases that diverged prior to the last common ancestor of extant life

ESTs from the microsporidian Edhazardia aedis

Background: 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 Encephalitozoon cuniculi and the locust parasite Antonospora locustae). We have conducted an EST survey of the mosquito parasite Edhazardia aedis, which has an estimated genome size several times that of more well-studied species. The only other microsporidian EST project is from A. locustae, and serves as a basis for comparison with E. aedis. Results The spore transcriptomes of A. locustae and E. aedis 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, E. aedis seems to lack the multi-gene transcripts present in A. locustae and E. cuniculi. We also present the first documented case of transcription of a transposable element in microsporidia. Conclusion Although E. aedis and A. locustae 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 E. aedis that have not been found in other microsporidia suggests that extreme genome reduction and compaction is lineage specific and not typical of all microsporidia.Botany, Department ofScience, Faculty ofNon UBCReviewedFacult

Transcriptome analysis of the parasite Encephalitozoon cuniculi: an in-depth examination of pre-mRNA splicing in a reduced eukaryote

Background: The microsporidian Encephalitozoon cuniculi possesses one of the most reduced and compacted eukaryotic genomes. Reduction in this intracellular parasite has affected major cellular machinery, including the loss of over fifty core spliceosomal components compared to S. cerevisiae. To identify expression changes throughout the parasite’s life cycle and also to assess splicing in the context of this reduced system, we examined the transcriptome of E. cuniculi using Illumina RNA-seq. Results We observed that nearly all genes are expressed at three post-infection time-points examined. A large fraction of genes are differentially expressed between the first and second (37.7%) and first and third (43.8%) time-points, while only four genes are differentially expressed between the latter two. Levels of intron splicing are very low, with 81% of junctions spliced at levels below 50%. This is dramatically lower than splicing levels found in two other fungal species examined. We also describe the first case of alternative splicing in a microsporidian, an unexpected complexity given the reduction in spliceosomal components. Conclusions Low levels of splicing observed are likely the result of an inefficient spliceosome; however, at least in one case, splicing appears to be playing a functional role. Although several RNA decay genes are encoded in E. cuniculi, the lack of a few key players could be reducing decay levels and therefore increasing the proportion of unspliced transcripts. Significant proportions of genes are differentially expressed in the first forty-eight hours but not after, indicative of genetic changes that precede the intracellular to infective stage transition.Botany, Department ofNon UBCScience, Faculty ofReviewedFacult

Bacterial Catalase in the Microsporidian Nosema locustae: Implications for Microsporidian Metabolism and Genome Evolution

Microsporidia constitute a group of extremely specialized intracellular parasites that infect virtually all animals. They are highly derived, reduced fungi that lack several features typical of other eukaryotes, including canonical mitochondria, flagella, and peroxisomes. Consistent with the absence of peroxisomes in microsporidia, the recently completed genome of the microsporidian Encephalitozoon cuniculi lacks a gene for catalase, the major enzymatic marker for the organelle. We show, however, that the genome of the microsporidian Nosema locustae, in contrast to that of E. cuniculi, encodes a group II large-subunit catalase. Surprisingly, phylogenetic analyses indicate that the N. locustae catalase is not specifically related to fungal homologs, as one would expect, but is instead closely related to proteobacterial sequences. This finding indicates that the N. locustae catalase is derived by lateral gene transfer from a bacterium. The catalase gene is adjacent to a large region of the genome that appears to be far less compact than is typical of microsporidian genomes, a characteristic which may make this region more amenable to the insertion of foreign genes. The N. locustae catalase gene is expressed in spores, and the protein is detectable by Western blotting. This type of catalase is a particularly robust enzyme that has been shown to function in dormant cells, indicating that the N. locustae catalase may play some functional role in the spore. There is no evidence that the N. locustae catalase functions in a cryptic peroxisome