Evolution of Dicer and Argonaute orthologs in microsporidian parasites

Non-coding RNAs identification and regulatory networks in pathogen-host interaction in the microsporidia congenital infection

BACKGROUND

The interaction networks between coding and non-coding RNAs (ncRNAs) including long non-coding RNA (lncRNA), covalently closed circular RNA (circRNA) and miRNA are significant to elucidate molecular processes of biological activities and interactions between host and pathogen. Congenital infection caused by vertical transmission of microsporidia N. bombycis can result in severe economic losses in the silkworm-feeding industry. However, little is known about ncRNAs that take place in the microsporidia congenital infection. Here we conducted whole-transcriptome RNA-Seq analyses to identify ncRNAs and regulatory networks for both N. bombycis and host including silkworm embryos and larvae during the microsporidia congenital infection.

RESULTS

A total of 4,171 mRNAs, 403 lncRNA, 62 circRNAs, and 284 miRNAs encoded by N. bombycis were identified, among which some differentially expressed genes formed cross-talk and are involved in N. bombycis proliferation and infection. For instance, a lncRNA/circRNA competing endogenous RNA (ceRNA) network including 18 lncRNAs, one circRNA, and 20 miRNAs was constructed to describe 14 key parasites genes regulation, such as polar tube protein 3 (PTP3), ricin-B-lectin, spore wall protein 4 (SWP4), and heat shock protein 90 (HSP90). Regarding host silkworm upon N. bombycis congenital infection, a total of 14,889 mRNAs, 3,038 lncRNAs, 19,039 circRNAs, and 3,413 miRNAs were predicted based on silkworm genome with many differentially expressed coding and non-coding genes during distinct developmental stages. Different species of RNAs form interacting network to modulate silkworm biological processes, such as growth, metamorphosis and immune responses. Furthermore, a lncRNA/circRNA ceRNA network consisting of 140 lncRNAs, five circRNA, and seven miRNAs are constructed hypothetically to describe eight key host genes regulation, such as Toll-6, Serpin-6, inducible nitric oxide synthase (iNOS) and Caspase-8. Notably, cross-species analyses indicate that parasite and host miRNAs play a vital role in pathogen-host interaction in the microsporidia congenital infection.

CONCLUSION

This is the first comprehensive pan-transcriptome study inclusive of both N. bombycis and its host silkworm with a specific focus on the microsporidia congenital infection, and show that ncRNA-mediated regulation plays a vital role in the microsporidia congenital infection, which provides a new insight into understanding the basic biology of microsporidia and pathogen-host interaction.

Genomic and phenotypic evolution of nematode-infecting microsporidia

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

Contrasting outcomes of genome reduction in mikrocytids and microsporidians

BACKGROUND

Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution.

RESULTS

At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths.

CONCLUSIONS

Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.

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.

Ricin B lectin-like proteins of the microsporidian Encephalitozoon cuniculi and Anncaliia algerae are involved in host-cell invasion

Microsporidia are obligate intracellular pathogens capable of infecting a wide variety of hosts ranging from invertebrates to vertebrates. The infection process requires a step of prior adherence of Microsporidia to the surface of host cells. A few studies demonstrated the involvement of proteins containing a ricin-B lectin (RBL) domain in parasite infection. In this study Anncalia algerae and Encephalitozoon cuniculi genomes were screened by bioinformatic analysis to identify proteins with an extracellular prediction and possessing RBL-type carbohydrate-binding domains, being both potentially relevant factors contributing to host cell adherence. Three proteins named AaRBLL-1 and AaRBLL-2 from A. algerae and EcRBLL-1 from E. cuniculi, were selected and comparative analysis of sequences suggested their belonging to a multigenic family, with a conserved structural RBL domain despite a significant amino acid sequence divergence. The production of recombinant proteins and antibodies against the three proteins allowed their subcellular localization on the spore wall and/or the polar tube. Adherence inhibition assays based on pre-treatments with recombinant proteins or antibodies highlighted the significant decrease of the proliferation of both E. cuniculi and A. algerae, strongly suggesting that these proteins are involved in the infection process.

Nosema bombycis microRNA-like RNA 8 (Nb-milR8) Increases Fungal Pathogenicity by Modulating BmPEX16 Gene Expression in Its Host, Bombyx mori

The fungus Nosema bombycis causes significant economic losses via parasitism of an economically important insect. MicroRNAs (miRNAs) play important roles in regulating host and parasite gene expression via mRNA degradation or by inhibiting protein translation. To investigate whether microRNA-like RNAs (milRNAs) regulate N. bombycis pathogenesis and to better understand the regulatory mechanisms underlying infection, we constructed small RNA libraries from N. bombycis hyphae during the schizont proliferation period. Eleven novel milRNAs were determined by RNA sequencing and stem-loop reverse transcriptase PCR (RT-PCR) assays. Moreover, a virulence-associated milRNA, Nb-milR8, was identified as critical for N. bombycis proliferation by binding and downregulating expression of its target gene, BmPEX16, in the host during infection. Silencing of Nb-milR8 or overexpression of the target BmPEX16 gene resulted in increased susceptibility of Bombyx mori to N. bombycis infection. Taken together, these results suggest that Nb-milR8 is an important virulence factor that acts as an effector to suppress host peroxidase metabolism, thereby facilitating N. bombycis proliferation. These results provide important novel insights into interactions between pathogenic fungi and their hosts.

IMPORTANCE A thorough understanding of fungal pathogen adaptations is essential for treating fungal infections. Recent studies have suggested that the role of small RNAs expressed in fungal microsporidia genomes are important for elucidating the mechanisms of fungal infections. Here, we report 11 novel microRNA-like RNAs (milRNAs) from the fungal microsporidium Nosema bombycis and identified NB-milRNAs that adaptively regulate N. bombycis proliferation. In addition, we demonstrate that N. bombycis modulates small RNA (sRNA)-mediated infection by encoding an Nb-miR8 that downregulates the expression of the host peroxidase metabolism protein BmPEX16, which is essential for peroxisome membrane biogenesis and peroxisome assembly. These results significantly contribute to our understanding of the pathogenic mechanisms of fungi, and especially microsporidia, while providing important targets for genetical engineering-based treatment of microsporidia.

Genome and Evolutionary Analysis of Nosema ceranae: A Microsporidian Parasite of Honey Bees

Microsporidia comprise a phylum of single cell, intracellular parasites and represent the earliest diverging branch in the fungal kingdom. The microsporidian parasite Nosema ceranae primarily infects honey bee gut epithelial cells, leading to impaired memory, suppressed host immune responses and colony collapse under certain circumstances. As the genome of N. ceranae is challenging to assembly due to very high genetic diversity and repetitive region, the genome was re-sequenced using long reads. We present a robust 8.8 Mbp genome assembly of 2,280 protein coding genes, including a high number of genes involved in transporting nutrients and energy, as well as drug resistance when compared with sister species Nosema apis. We also describe the loss of the critical protein Dicer in approximately half of the microsporidian species, giving new insights into the availability of RNA interference pathway in this group. Our results provided new insights into the pathogenesis of N. ceranae and a blueprint for treatment strategies that target this parasite without harming honey bees. The unique infectious apparatus polar filament and transportation pathway members can help to identify treatments to control this parasite.

Identification of novel miRNAs from the microsporidian parasite Nosema ceranae

Previously, six miRNAs were identified from the microsporidian parasite Nosema ceranae. By taking advantage of the recently updated N. ceranae and honey bee genome assemblies, we re-analyzed the deep sequencing datasets. Three novel miRNAs were identified, which were further validated by plasmid cloning and sequencing. The miRNAs correlated with significantly higher number of genes from the parasite than the host. Our data suggest the parasitic miRNAs are involved in self-regulation during the proliferation.

Transposable element abundance correlates with mode of transmission in microsporidian parasites

The extreme genome reduction and physiological simplicity of some microsporidia has been attributed to their intracellular, obligate parasitic lifestyle. Although not all microsporidian genomes are small (size range from about 2 to 50 MB), it is suggested that the size of their genomes has been streamlined by natural selection. We explore the hypothesis that vertical transmission in microsporidia produces population bottlenecks, and thus reduces the effectiveness of natural selection. Here we compare the transposable element (TE) content of 47 microsporidian genomes, and show that genome size is positively correlated with the amount of TEs, and that species that experience vertical transmission have larger genomes with higher proportion of TEs. Our findings are consistent with earlier studies inferring that nonadaptive processes play an important role in microsporidian evolution.