Gain and loss of multiple functionally related, horizontally transferred genes in the reduced genomes of two microsporidian parasites

How to overcome constraints imposed by microsporidian genome features to ensure gene prediction?

Since the advent of sequencing techniques and due to their continuous evolution, it has become easier and less expensive to obtain the complete genome sequence of any organism. Nevertheless, to elucidate all biological processes governing organism development, quality annotation is essential. In genome annotation, predicting gene structure is one of the most important and captivating challenges for computational biology. This aspect of annotation requires continual optimization, particularly for genomes as unusual as those of microsporidia. Indeed, this group of fungal-related parasites exhibits specific features (highly reduced gene sizes, sequences with high rate of evolution) linked to their evolution as intracellular parasites, requiring the implementation of specific annotation approaches to consider all these features. This review aimed to outline these characteristics and to assess the increasingly efficient approaches and tools that have enhanced the accuracy of gene prediction for microsporidia, both in terms of sensitivity and specificity. Subsequently, a final part will be dedicated to postgenomic approaches aimed at reinforcing the annotation data generated by prediction software. These approaches include the characterization of other understudied genes, such as those encoding regulatory noncoding RNAs or very small proteins, which also play crucial roles in the life cycle of these microorganisms.

Genome evolution in intracellular parasites: Microsporidia and Apicomplexa

Microsporidia and Apicomplexa are eukaryotic, single-celled, intracellular parasites with huge public health and economic importance. Typically, these parasites are studied separately, emphasizing their uniqueness and diversity. In this review, we explore the huge amount of genomic data that has recently become available for the two groups. We compare and contrast their genome evolution and discuss how their transitions to intracellular life may have shaped it. In particular, we explore genome reduction and compaction, genome expansion and ploidy, gene shuffling and rearrangements, and the evolution of centromeres and telomeres.

MicroAnnot: A Dedicated Workflow for Accurate Microsporidian Genome Annotation

With nearly 1700 species, Microsporidia represent a group of obligate intracellular eukaryotes with veterinary, economic and medical impacts. To help understand the biological functions of these microorganisms, complete genome sequencing is routinely used. Nevertheless, the proper prediction of their gene catalogue is challenging due to their taxon-specific evolutionary features. As innovative genome annotation strategies are needed to obtain a representative snapshot of the overall lifestyle of these parasites, the MicroAnnot tool, a dedicated workflow for microsporidian sequence annotation using data from curated databases of accurately annotated microsporidian genes, has been developed. Furthermore, specific modules have been implemented to perform small gene (<300 bp) and transposable element identification. Finally, functional annotation was performed using the signature-based InterProScan software. MicroAnnot's accuracy has been verified by the re-annotation of four microsporidian genomes for which structural annotation had previously been validated. With its comparative approach and transcriptional signal identification method, MicroAnnot provides an accurate prediction of translation initiation sites, an efficient identification of transposable elements, as well as high specificity and sensitivity for microsporidian genes, including those under 300 bp.

Biogenesis, inheritance, and 3D ultrastructure of the microsporidian mitosome

During the reductive evolution of obligate intracellular parasites called microsporidia, a tiny remnant mitochondrion (mitosome) lost its typical cristae, organellar genome, and most canonical functions. Here, we combine electron tomography, stereology, immunofluorescence microscopy, and bioinformatics to characterise mechanisms of growth, division, and inheritance of this minimal mitochondrion in two microsporidia species (grown within a mammalian RK13 culture-cell host). Mitosomes of Encephalitozoon cuniculi (2-12/cell) and Trachipleistophora hominis (14-18/nucleus) displayed incremental/non-phasic growth and division and were closely associated with an organelle identified as equivalent to the fungal microtubule-organising centre (microsporidian spindle pole body; mSPB). The mitosome-mSPB association was resistant to treatment with microtubule-depolymerising drugs nocodazole and albendazole. Dynamin inhibitors (dynasore and Mdivi-1) arrested mitosome division but not growth, whereas bioinformatics revealed putative dynamins Drp-1 and Vps-1, of which, Vps-1 rescued mitochondrial constriction in dynamin-deficient yeast (Schizosaccharomyces pombe). Thus, microsporidian mitosomes undergo incremental growth and dynamin-mediated division and are maintained through ordered inheritance, likely mediated via binding to the microsporidian centrosome (mSPB).

Molecular and cellular characterization of four putative nucleotide transporters from the shrimp microsporidian Enterocytozoon hepatopenaei (EHP)

Microsporidia are obligate intracellular parasites that lost several enzymes required in energy production. The expansion of transporter families in these organisms enables them to hijack ATP from hosts. In this study, nucleotide transporters of the microsporidian Enterocytozoon hepatopenaei (EHP), which causes slow growth in economically valuable Penaeus shrimp, were characterized. Analysis of the EHP genome suggested the presence of four putative nucleotide transporter genes, namely EhNTT1, EhNTT2, EhNTT3, and EhNTT4. Sequence alignment revealed four charged amino acids that are conserved in previously characterized nucleotide transporters. Phylogenetic analysis suggested that EhNTT1, 3, and 4 were derived from one horizontal gene transfer event, which was independent from that of EhNTT2. Localization of EhNTT1 and EhNTT2 using immunofluorescence analysis revealed positive signals within the envelope of developing plasmodia and on mature spores. Knockdown of EhNTT2 by double administration of sequence specific double-stranded RNA resulted in a significant reduction in EHP copy numbers, suggesting that EhNTT2 is crucial for EHP replication in shrimp. Taken together, the insight into the roles of NTTs in microsporidian proliferation can provide the biological basis for the development of alternative control strategies for microsporidian infection in shrimp.

Fungal taxonomy: A puzzle with many missing pieces

Fungi are multifaceted organisms found in almost all ecosystems on Earth, where they establish various types of symbiosis with other living beings. Despite being recognized by humans since ancient times, and the high number of works delving into their biology and ecology, much is still unknown about these organisms. Some criteria classically used for their study are nowadays limited, generating confusion in categorizing them, and even more, when trying to understand their genealogical relationships. To identify species within Fungi, phenotypic characters to date are not sufficient, and to construct a broad phylogeny or a phylogeny of a particular group, there are still gaps affecting the generated trees, making them unstable and easily debated. For health professionals, fungal identification at lower levels such as genus and species, is enough to select the most appropriate therapy for their control, understand the epidemiology of clinical pictures associated, and recognize outbreaks and antimicrobial resistance. However, the taxonomic location within the kingdom, information with apparently little relevance, can allow phylogenetic relationships to be established between fungal taxa, facilitating the understanding of their biology, distribution in nature, and pathogenic potential evolution. Advances in molecular biology and computer science techniques from the last 30 years have led to crucial changes aiming to establish the criteria to define a fungal species, allowing us to reach a kind of stable phylogenetic construction. However, there is still a long way to go, and it requires the joint work of the scientific community at a global level and support for basic research.

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.

Structural and Functional Annotation of Hypothetical Proteins from the Microsporidia Species Vittaforma corneae ATCC 50505 Using in silico Approaches

Microsporidia are spore-forming eukaryotes that are related to fungi but have unique traits that set them apart. They have compact genomes as a result of evolutionary gene loss associated with their complete dependency on hosts for survival. Despite having a relatively small number of genes, a disproportionately high percentage of the genes in microsporidia genomes code for proteins whose functions remain unknown (hypothetical proteins-HPs). Computational annotation of HPs has become a more efficient and cost-effective alternative to experimental investigation. This research developed a robust bioinformatics annotation pipeline of HPs from Vittaforma corneae, a clinically important microsporidian that causes ocular infections in immunocompromised individuals. Here, we describe various steps to retrieve sequences and homologs and to carry out physicochemical characterization, protein family classification, identification of motifs and domains, protein-protein interaction network analysis, and homology modelling using a variety of online resources. Classification of protein families produced consistent findings across platforms, demonstrating the accuracy of annotation utilizing in silico methods. A total of 162 out of 2034 HPs were fully annotated, with the bulk of them categorized as binding proteins, enzymes, or regulatory proteins. The protein functions of several HPs from Vittaforma corneae were accurately inferred. This improved our understanding of microsporidian HPs despite challenges related to the obligate nature of microsporidia, the absence of fully characterized genes, and the lack of homologous genes in other systems.

Heterologous Expressed NbSWP12 from Microsporidia Nosema bombycis Can Bind with Phosphatidylinositol 3-phosphate and Affect Vesicle Genesis

Microsporidia are a big group of single-celled obligate intracellular organisms infecting most animals and some protozoans. These minimalist eukaryotes lack numerous genes in metabolism and vesicle trafficking. Here, we demonstrated that the spore wall protein NbSWP12 of microsporidium Nosema bombycis belongs to Bin/Amphiphysin/Rvs (BAR) protein family and can specifically bind with phosphatidylinositol 3-phosphate [Ptdlns(3)P]. Since Ptdlns(3)P is involved in endosomal vesicle biogenesis and trafficking, we heterologous expressed NbSWP12 in yeast Saccharomyces cerevisiae and proved that NbSWP12 can target the cell membrane and endocytic vesicles. Nbswp12 transformed into Gvp36 (a BAR protein of S. cerevisiae) deletion mutant rescued the defect phenotype of vesicular traffic. This study identified a BAR protein function in vesicle genesis and sorting and provided clues for further understanding of how microsporidia internalize nutrients and metabolites during proliferation.

Identification and subcellular localization analysis of membrane protein Ycf 1 in the microsporidian Nosema bombycis

Microsporidia are obligate intracellular parasites that can infect a wide range of vertebrates and invertebrates including humans and insects, such as silkworm and bees. The microsporidium Nosema bombycis can cause pebrine in Bombyx mori, which is the most destructive disease in the sericulture industry. Although membrane proteins are involved in a wide range of cellular functions and part of many important metabolic pathways, there are rare reports about the membrane proteins of microsporidia up to now. We screened a putative membrane protein Ycf 1 from the midgut transcriptome of the N. bombycis-infected silkworm. Gene cloning and bioinformatics analysis showed that the Ycf 1 gene contains a complete open reading frame (ORF) of 969 bp in length encoding a 322 amino acid polypeptide that has one signal peptide and one transmembrane domain. Indirect immunofluorescence results showed that Ycf 1 protein is distributed on the plasma membrane. Expression pattern analysis showed that the Ycf 1 gene expressed in all developmental stages of N. bombycis. Knockdown of the Ycf 1 gene by RNAi effectively inhibited the proliferation of N. bombycis. These results indicated that Ycf 1 is a membrane protein and plays an important role in the life cycle of N. bombycis.

Identification and Characterization of Three Spore Wall Proteins of Enterocytozoon Bieneusi

Enterocytozoon bieneusi is the most common microsporidian pathogen in farm animals and humans. Although several spore wall proteins (SWPs) of other human-pathogenic microsporidia have been identified, SWPs of E. bieneusi remain poorly characterized. In the present study, we identified the sequences of three E. bieneusi SWPs from whole genome sequence data, expressed them in Escherichia coli, generated a monoclonal antibody (mAb) against one of them (EbSWP1), and used the mAb in direct immunofluorescence detection of E. bieneusi spores in fecal samples. The amino acid sequence of EbSWP1 shares some identity to EbSWP2 with a BAR2 domain, while the sequence of EbSWP3 contains a MICSWaP domain. No cross-reactivity among the EbSWPs was demonstrated using the polyclonal antibodies generated against them. The mAb against EbSWP1 was shown to react with E. bieneusi spores in fecal samples. Using chromotrope 2R staining-based microscopy as the gold standard, the sensitivity and specificity of the direct immunofluorescence for the detection of E. bieneusi were 91.4 and 73.7%. Data generated from the study could be useful in the characterization of E. bieneusi and immunological detection of the pathogen.

A new microsporidian parasite, Ordospora pajunii sp. nov (Ordosporidae), of Daphnia longispina highlights the value of genomic data for delineating species boundaries

Speciation is a complex and continuous process that makes the delineation of species boundaries a challenging task in particular in species with little morphological differentiation, such as parasites. In this case, the use of genomic data is often necessary, such as for the intracellular Microsporidian parasites. Here we characterize the genome of a gut parasite of the cladoceran Daphnia longispina (isolate FI-F-10), which we propose as a new species within the genus Ordospora: O. pajunii sp. nov (Ordosporidae). FI-F-10 closest relative, O. colligata is only found in D. magna. Both microsporidian species share several morphological features. Although it is not possible to estimate divergence times for Microsporidia due to the lack of fossil records and accelerated evolutionary rates, we base our proposal on the phylogenomic and genomic distances between both microsporidian lineages. Phylogenomic reconstruction shows that FI-F-10 forms an early diverging branch basal to the cluster that contains all known O. colligata strains. Whole-genome comparisons show that FI-F-10 presents a greater divergence at the sequence level than observed among O. colligata strains, and its genomic Average Nucleotide Identity (ANI) values against O. colligata are beyond the intra-specific range previously established for yeast and prokaryotes. Our data confirm that the ANI metrics are useful for fine genetic divergence calibration across Microsporidia taxa. In combination with phylogenetic and ecological data, genome-based metrics provide a powerful approach to delimitate species boundaries.

Current Therapy and Therapeutic Targets for Microsporidiosis

Microsporidia are obligate intracellular, spore-forming parasitic fungi which are grouped with the Cryptomycota. They are both opportunistic pathogens in humans and emerging veterinary pathogens. In humans, they cause chronic diarrhea in immune-compromised patients and infection is associated with increased mortality. Besides their role in pébrine in sericulture, which was described in 1865, the prevalence and severity of microsporidiosis in beekeeping and aquaculture has increased markedly in recent decades. Therapy for these pathogens in medicine, veterinary, and agriculture has become a recent focus of attention. Currently, there are only a few commercially available antimicrosporidial drugs. New therapeutic agents are needed for these infections and this is an active area of investigation. In this article we provide a comprehensive summary of the current as well as several promising new agents for the treatment of microsporidiosis including: albendazole, fumagillin, nikkomycin, orlistat, synthetic polyamines, and quinolones. Therapeutic targets which could be utilized for the design of new drugs are also discussed including: tubulin, type 2 methionine aminopeptidase, polyamines, chitin synthases, topoisomerase IV, triosephosphate isomerase, and lipase. We also summarize reports on the utility of complementary and alternative medicine strategies including herbal extracts, propolis, and probiotics. This review should help facilitate drug development for combating microsporidiosis.

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.

OUP accepted manuscript

DNA repair is an important component of genome integrity and organisms with reduced repair capabilities tend to accumulate mutations at elevated rates. Microsporidia are intracellular parasites exhibiting high levels of genetic divergence postulated to originate from the lack of several proteins, including the heterotrimeric Rad9–Rad1–Hus1 DNA repair clamp. Microsporidian species from the Encephalitozoonidae have undergone severe streamlining with small genomes coding for about 2,000 proteins. The highly divergent sequences found in Microsporidia render functional inferences difficult such that roughly half of these 2,000 proteins have no known function. Using a structural homology-based annotation approach combining protein structure prediction and tridimensional similarity searches, we found that the Rad9–Rad1–Hus1 DNA clamp is present in Microsporidia, together with many other components of the DNA repair machinery previously thought to be missing from these organisms. Altogether, our results indicate that the DNA repair machinery is present and likely functional in Microsporidia.

Microsporidiosis in Humans

Microsporidia are obligate intracellular pathogens identified ∼150 years ago as the cause of pébrine, an economically important infection in silkworms. There are about 220 genera and 1,700 species of microsporidia, which are classified based on their ultrastructural features, developmental cycle, host-parasite relationship, and molecular analysis. Phylogenetic analysis suggests that microsporidia are related to the fungi, being grouped with the Cryptomycota as a basal branch or sister group to the fungi. Microsporidia can be transmitted by food and water and are likely zoonotic, as they parasitize a wide range of invertebrate and vertebrate hosts. Infection in humans occurs in both immunocompetent and immunodeficient hosts, e.g., in patients with organ transplantation, patients with advanced human immunodeficiency virus (HIV) infection, and patients receiving immune modulatory therapy such as anti-tumor necrosis factor alpha antibody. Clusters of infections due to latent infection in transplanted organs have also been demonstrated. Gastrointestinal infection is the most common manifestation; however, microsporidia can infect virtually any organ system, and infection has resulted in keratitis, myositis, cholecystitis, sinusitis, and encephalitis. Both albendazole and fumagillin have efficacy for the treatment of various species of microsporidia; however, albendazole has limited efficacy for the treatment of Enterocytozoon bieneusi. In addition, immune restoration can lead to resolution of infection. While the prevalence rate of microsporidiosis in patients with AIDS has fallen in the United States, due to the widespread use of combination antiretroviral therapy (cART), infection continues to occur throughout the world and is still seen in the United States in the setting of cART if a low CD4 count persists.

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.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

The use of flow cytometry for fungal nuclear DNA quantification

Genome size information is sparse across fungi, with information being available for less than 2000 species. So far, most records have been obtained using static, microscope-based cytometry methods or derived from genome sequencing projects. Flow cytometry is now considered the state-of-the-art method for obtaining genome size measurements, and appropriate methods and DNA standards are available, enabling the analysis of most genome size ranges in a rapid, robust and inexpensive way. The average fungal genome size is 60 Mbp, but sizes vary across phylogeny, ranging from 2.2 (Encephalitozoon romaleae) to 3706 Mbp (Jafnea semitosta). In several fungal clades, genome size expansion seems to accompany evolution either to plant mutualism or to plant parasitism (particularly biotrophy), and fungi that interact with plants seem to have larger genomes than saprobes and those that interact with animals. Whereas flow cytometry for nuclear DNA quantification is routinely employed in plant sciences for genome size and ploidy studies, its use in fungal biology is still infrequent. Appropriate standards, methods and best practices are described here, with the aim of stimulating a more generalized and widespread use of flow cytometry for fungal genome size measurement.

Secretomes of medically important fungi reflect morphological and phylogenetic diversity

Secretome represents a main target for understanding the mechanisms of fungal adaptation. In the present study, we focus on the secretomes of fungi associated with infections in humans and other mammals in order to explore relationships between the diverse morphological and phylogenetic groups. Almost all the mammalian pathogenic fungi analyzed have secretome sizes smaller than 1000 proteins and, secreted proteins comprise between 5% and 10% of the total proteome. As expected, the correlation pattern between the secretome size and the total proteome was similar to that described in previous secretome studies of fungi. With regard to the morphological groups, minimum secretome sizes of less than 250 secreted proteins and low values for the fraction of secreted proteins are shown in mammalian pathogenic fungi with reduced proteomes such as microsporidia, atypical fungi and some species of yeasts and yeast-like fungi (Malassezia). On the other hand, filamentous fungi have significantly more secreted proteins and the highest numbers are present in species of filamentous fungi that also are plant or insect pathogens (Fusarium verticilloides, Fusarium oxysporum and Basidiobolus meristosporus). With respect to phylogeny, there are also variations in secretome size across fungal subphyla: Microsporidia, Taphrinomycotina, Ustilagomycotina and Saccharomycotina contain small secretomes; whereas larger secretomes are found in Agaricomycotina, Pezizomycotina, Mucoromycotina and Entomophthoromycotina. Finally, principal component analysis (PCA) was conducted on the complete secretomes. The PCA results revealed that, in general, secretomes of fungi belonging to the same morphological group or subphyla cluster together. In conclusion, our results point out that in medically important fungi there is a relationship between the secretome and the morphological group or phylogenetic classification.

Advances in understanding the evolution of fungal genome architecture

Diversity within the fungal kingdom is evident from the wide range of morphologies fungi display as well as the various ecological roles and industrial purposes they serve. Technological advances, particularly in long-read sequencing, coupled with the increasing efficiency and decreasing costs across sequencing platforms have enabled robust characterization of fungal genomes. These sequencing efforts continue to reveal the rampant diversity in fungi at the genome level. Here, we discuss studies that have furthered our understanding of fungal genetic diversity and genomic evolution. These studies revealed the presence of both small-scale and large-scale genomic changes. In fungi, research has recently focused on many small-scale changes, such as how hypermutation and allelic transmission impact genome evolution as well as how and why a few specific genomic regions are more susceptible to rapid evolution than others. High-throughput sequencing of a diverse set of fungal genomes has also illuminated the frequency, mechanisms, and impacts of large-scale changes, which include chromosome structural variation and changes in chromosome number, such as aneuploidy, polyploidy, and the presence of supernumerary chromosomes. The studies discussed herein have provided great insight into how the architecture of the fungal genome varies within species and across the kingdom and how modern fungi may have evolved from the last common fungal ancestor and might also pave the way for understanding how genomic diversity has evolved in all domains of life.

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.

Horizontal Gene Transfer to a Defensive Symbiont with a Reduced Genome in a Multipartite Beetle Microbiome

Symbiotic mutualisms of bacteria and animals are ubiquitous in nature, running a continuum from facultative to obligate from the perspectives of both partners. The loss of functions required for living independently but not within a host gives rise to reduced genomes in many symbionts. Although the phenomenon of genome reduction can be explained by existing evolutionary models, the initiation of the process is not well understood. Here, we describe the microbiome associated with the eggs of the beetle Lagria villosa, consisting of multiple bacterial symbionts related to Burkholderia gladioli, including a reduced-genome symbiont thought to be the exclusive producer of the defensive compound lagriamide. We show that the putative lagriamide-producing symbiont is the only member of the microbiome undergoing genome reduction and that it has already lost the majority of its primary metabolism and DNA repair pathways. The key step preceding genome reduction in the symbiont was likely the horizontal acquisition of the putative lagriamide lga biosynthetic gene cluster. Unexpectedly, we uncovered evidence of additional horizontal transfers to the symbiont's genome while genome reduction was occurring and despite a current lack of genes needed for homologous recombination. These gene gains may have given the genome-reduced symbiont a selective advantage in the microbiome, especially given the maintenance of the large lga gene cluster despite ongoing genome reduction.IMPORTANCE Associations between microorganisms and an animal, plant, or fungal host can result in increased dependence over time. This process is due partly to the bacterium not needing to produce nutrients that the host provides, leading to loss of genes that it would need to live independently and to a consequent reduction in genome size. It is often thought that genome reduction is aided by genetic isolation-bacteria that live in monocultures in special host organs, or inside host cells, have less access to other bacterial species from which they can obtain genes. Here, we describe exposure of a genome-reduced beetle symbiont to a community of related bacteria with nonreduced genomes. We show that the symbiont has acquired genes from other bacteria despite going through genome reduction, suggesting that isolation has not yet played a major role in this case of genome reduction, with horizontal gene gains still offering a potential route for adaptation.

Microsporidia with Vertical Transmission Were Likely Shaped by Nonadaptive Processes

Microsporidia have the leanest genomes among eukaryotes, and their physiological and genomic simplicity has been attributed to their intracellular, obligate parasitic life-style. However, not all microsporidia genomes are small or lean, with the largest dwarfing the smallest ones by at least an order of magnitude. To better understand the evolutionary mechanisms behind this genomic diversification, we explore here two clades of microsporidia with distinct life histories, Ordospora and Hamiltosporidium, parasitizing the same host species, Daphnia magna. Based on seven newly assembled genomes, we show that mixed-mode transmission (the combination of horizontal and vertical transmission), which occurs in Hamiltosporidium, is found to be associated with larger and AT-biased genomes, more genes, and longer intergenic regions, as compared with the exclusively horizontally transmitted Ordospora. Furthermore, the Hamiltosporidium genome assemblies contain a variety of repetitive elements and long segmental duplications. We show that there is an excess of nonsynonymous substitutions in the microsporidia with mixed-mode transmission, which cannot be solely attributed to the lack of recombination, suggesting that bursts of genome size in these microsporidia result primarily from genetic drift. Overall, these findings suggest that the switch from a horizontal-only to a mixed mode of transmission likely produces population bottlenecks in Hamiltosporidium species, therefore reducing the effectiveness of natural selection, and allowing their genomic features to be largely shaped by nonadaptive processes.

Molecular Dating of the Emergence of Anaerobic Rumen Fungi and the Impact of Laterally Acquired Genes

The anaerobic gut fungi (AGF), or Neocallimastigomycota, inhabit the rumen and alimentary tract of herbivorous mammals, where they play important roles in the degradation of plant fiber. Comparative genomic and phylogenomic analyses of the AGF have long been hampered by their fastidious growth condition, as well as their large (up to 200 Mb) and AT-biased (78 to 84%) genomes. We sequenced 21 AGF transcriptomes and combined them with 5 available AGF genome sequences to explore their evolutionary relationships, time their divergence, and characterize gene gain/loss patterns associated with their evolution. We estimate that the most recent common ancestor of the AGF diverged 66 (±10) million years ago, a time frame that coincides with the evolution of grasses (Poaceae), as well as the mammalian transition from insectivory to herbivory. The concordance of independent estimations suggests that AGF have been important in shaping the success of mammalian herbivory transition by improving the efficiency of energy acquisition from recalcitrant plant materials. Comparative genomics identified multiple lineage-specific genes in the AGF, two of which were acquired from rumen gut bacteria and animal hosts via horizontal gene transfer (HGT). A third AGF domain, plant-like polysaccharide lyase, represents a novel gene in fungi that potentially aids AGF to degrade pectin. Analysis of genomic and transcriptomic sequences confirmed both the presence and expression of these lineage-specific genes in nearly all AGF clades. These genetic elements may contribute to the exceptional abilities of AGF to degrade plant biomass and enable metabolism of the rumen microbes and animal hosts.IMPORTANCE Anaerobic fungi living in the rumen of herbivorous mammals possess an extraordinary ability to degrade plant biomass. We examined the origin and genomic composition of these poorly characterized anaerobic gut fungi using both transcriptome and genomic data. Phylogenomics and molecular dating analyses found remarkable concurrence of the divergence times of the rumen fungi, the forage grasses, and the dietary shift of ancestral mammals from primarily insectivory to herbivory. Comparative genomics identified unique machinery in these fungi to utilize plant polysaccharides. The rumen fungi were also identified with the ability to code for three protein domains with putative functions in plant pectin degradation and microbial defense, which were absent from all other fungal organisms (examined over 1,000 fungal genomes). Two of these domains were likely acquired from rumen gut bacteria and animal hosts separately via horizontal gene transfer. The third one is a plant-like polysaccharide lyase, representing a unique fungal enzyme with potential pectin breakdown abilities.

Genetic and Genome Analyses Reveal Genetically Distinct Populations of the Bee Pathogen Nosema ceranae from Thailand

The recent global decline in Western honeybee (Apis mellifera) populations is of great concern for pollination and honey production worldwide. Declining honeybee populations are frequently infected by the microsporidian pathogen Nosema ceranae. This species was originally described in the Asiatic honeybee (Apis cerana), and its identification in global A. mellifera hives could result from a recent host transfer. Recent genome studies have found that global populations of this parasite are polyploid and that humans may have fueled their global expansion. To better understand N. ceranae biology, we investigated its genetic diversity within part of their native range (Thailand) and among different hosts (A. mellifera, A. cerana) using both PCR and genome-based methods. We find that Thai N. ceranae populations share many SNPs with other global populations and appear to be clonal. However, in stark contrast with previous studies, we found that these populations also carry many SNPs not found elsewhere, indicating that these populations have evolved in their current geographic location for some time. Our genome analyses also indicate the potential presence of diploidy within Thai populations of N. ceranae.

Multilocus Sequence Typing and Population Genetic Analysis of Enterocytozoon bieneusi: Host Specificity and Its Impacts on Public Health

Microsporidia comprise a large class of unicellular eukaryotic pathogens that are medically and agriculturally important, but poorly understood. There have been nearly 1,500 microsporidian species described thus far, which are variable in biology, genetics, genomics, and host specificity. Among those, Enterocytozoon bieneusi is the well-known species responsible for the most recorded cases of human microsporidian affections. The pathogen can colonize a broad range of mammals and birds and most of the animals surveyed share some genotypes with humans, posing a threat to public health. Based on DNA sequence analysis of the ribosomal internal transcribed spacer (ITS) and phylogenetic analysis, several hundreds of E. bieneusi genotypes have been defined and clustered into different genetic groups with varied levels of host specificity. However, single locus-based typing using ITS might have insufficient resolution to discriminate among E. bieneusi isolates with complex genetic or hereditary characteristics and to assess the elusive reproduction or transmission modes of the organism, highlighting the need for exploration and application of multilocus sequence typing (MLST) and population genetic tools. The present review begins with a primer on microsporidia and major microsporidian species, briefly introduces the recent advances on E. bieneusi ITS genotyping and phylogeny, summarizes recent MLST and population genetic data, analyzes the inter- and intragroup host specificity at the MLST level, and interprets the public health implications of host specificity in zoonotic or cross-species transmission of this ubiquitous fungus.

Dicer regulates Nosema ceranae proliferation in honeybees

Nosema ceranae is a microsporidian parasite that infects the honeybee midgut epithelium. The protein-coding gene Dicer is lost in most microsporidian genomes but is present in N. ceranae. By feeding infected honeybees with small interfering RNA targeting the N. ceranae gene coding Dicer (siRNA-Dicer), we found that N. ceranae spore loads were significantly reduced. In addition, over 10% of total parasite protein-coding genes showed significantly divergent expression profiles after siRNA-Dicer treatment. Parasite genes for cell proliferation, ABC transporters and hexokinase were downregulated at 3 days postinfection, a key point in the middle of parasite replication cycles. In addition, genes involved in metabolic pathways of honeybees and N. ceranae showed significant co-expression. Furthermore, the siRNA-Dicer treatment partly reversed the expression patterns of honeybee genes. The honeybee gene mucin-2-like showed significantly upregulation in the siRNA-Dicer group compared with the infection group continually at 4, 5 and 6 days postinfection, suggesting that the siRNA-Dicer feeding promoted the strength of the mucus barrier resulted from interrupted parasite proliferation. As the gene Dicer broadly regulates N. ceranae proliferation and honeybee metabolism, our data suggest the RNA interference pathway is an important infection strategy for N. ceranae.

Sense-antisense gene overlap is probably a cause for retaining the few introns in Giardia genome and the implications

Background

It is widely accepted that the last eukaryotic common ancestor and early eukaryotes were intron-rich and intron loss dominated subsequent evolution, thus the presence of only very few introns in some modern eukaryotes must be the consequence of massive loss. But it is striking that few eukaryotes were found to have completely lost introns. Despite extensive research, the causes of massive intron losses remain elusive. Actually the reverse question -- how the few introns can be retained under the evolutionary selection pressure of intron loss -- is equally significant but was rarely studied, except that it was conjectured that the essential functions of some introns prevent their loss. The situation that extremely few (eight) spliceosome-mediated cis-spliced introns present in the relatively simple genome of Giardia lamblia provides an excellent opportunity to explore this question.

Results

Our investigation found three types of distribution patterns of the few introns in the intron-containing genes: ancient intron in ancient gene, later-evolved intron in ancient gene, and later-evolved intron in later-evolved gene, which can reflect to some extent the dynamic evolution of introns in Giardia. Without finding any special features or functional importance of these introns responsible for their retention, we noticed and experimentally verified that some intron-containing genes form sense-antisense gene pairs with transcribable genes on their complementary strands, and that the introns just reside in the overlapping regions.

Conclusions

In Giardia’s evolution, despite constant evolutionary selection pressure of intron loss, intron gain can still occur in both ancient and later-evolved genes, but only a few introns are retained; at least the evolutionary retention of some of the introns might not be due to the functional constraint of the introns themselves but the causes outside of introns, such as the constraints imposed by other genomic functional elements overlapping with the introns. These findings can not only provide some clues to find new genomic functional elements -- in the areas overlapping with introns, but suggest that “functional constraint” of introns may not be necessarily directly associated with intron loss and gain, and that the real functions are probably still outside of our current knowledge.

Reviewers

This article was reviewed by Mikhail Gelfand, Michael Gray, and Igor Rogozin.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0226-5) contains supplementary material, which is available to authorized users.

Evolutionary Genomics of Metchnikovella incurvata (Metchnikovellidae): An Early Branching Microsporidium

Metchnikovellids are highly specialized hyperparasites, which infect and reproduce inside gregarines (Apicomplexa) inhabiting marine invertebrates. Their phylogenetic affiliation was under constant discussion until recently, when analysis of the first near-complete metchnikovellid genome, that of Amphiamblys sp., placed it in a basal position with respect to most Microsporidia. Microsporidia are a highly diversified lineage of extremely reduced parasites related to Rozellida (Rozellosporidia = Rozellomycota = Cryptomycota) within the Holomycota clade of Opisthokonta. By sequencing DNA from a single-isolated infected gregarine cell we obtained an almost complete genome of a second metchnikovellid species, and the first one of a taxonomically described and well-documented species, Metchnikovella incurvata. Our phylogenomic analyses show that, despite being considerably divergent from each other, M. incurvata forms a monophyletic group with Amphiamplys sp., and confirm that metchnikovellids are one of the deep branches of Microsporidia. Comparative genomic analysis demonstrates that, like most Microsporidia, metchnikovellids lack mitochondrial genes involved in energy transduction and are thus incapable of synthesizing their own ATP via mitochondrial oxidative phosphorylation. They also lack the horizontally acquired ATP transporters widespread in most Microsporidia. We hypothesize that a family of mitochondrial carrier proteins evolved to transport ATP from the host into the metchnikovellid cell. We observe the progressive reduction of genes involved in DNA repair pathways along the evolutionary path of Microsporidia, which might explain, at least partly, the extremely high evolutionary rate of the most derived species. Our data also suggest that genome reduction and acquisition of novel genes co-occurred during the adaptation of Microsporidia to their hosts.

Evolution of Dicer and Argonaute orthologs in microsporidian parasites

Microsporidia are a group of intracellular parasites which infect animal hosts. The infection can broadly influence the hosts' metabolism, growth as well as immune responses. Recently, a functional RNAi pathway was suggested from the microsporidia parasite Nosema ceranae, whereby the gene Dicer showed strong impact on spore proliferation. Based on sequenced microsporidian species, the RNAi gene orthologs have only been annotated for a few species. In order to study the selection of RNAi gene Dicer and Argonaute orthologs from microsporidian genomes, a phylogenetic analysis was performed based on single copy orthologs of 21 microsporidian parasite species. Of the 21 studied parasite species, 11 parasite species maintained Dicer and Argonaute orthologs, which were further used to build the gene trees. The gene Dicer and Argonaute orthologs were either both maintained or both lost. The topology structures between the phylogenetic Dicer, Argonaute and species trees were consistent. The results suggest that the gene Dicer and Argonaute were selected as a unit, which were selectively maintained/lost during the lineage divergence. The study provides general insights on the selection of RNAi pathway in microsporidian parasites and the evolution of parasitism.

Invertebrate host responses to microsporidia infections

Microsporidia are a group of fungi-like intracellular and unicellular parasites, which infect nearly all animals. As "master parasites", over 1400 microsporidian species have been described to date. Microsporidia infections in economical invertebrates (e.g., silkworm, shrimp) cause huge financial losses, while other microsporidia infections in daphnia, nematode, locust, honeybee and mosquito play important roles in the regulation of their population size. Research investigating invertebrate host responses following microsporidia infections has yielded numerous interesting results, especially pertaining to the innate immune response to these pathogens. In this review, we comparatively summarize the invertebrate host responses to various microsporidia infections. We discuss numerous critical events in host responses including ubiquitin-mediated resistance, production of reactive oxygen species, melanization and innate immune pathways, and the increased basic metabolism and the accumulation of juvenile hormone in infected hosts. Recent studies progressing our understanding of microsporidia infection are also highlighted. Collectively, these advances shed more light on general rules of invertebrate host immune responses and pathogenesis mechanisms of microsporidia, and concurrently offer valuable clues for further research on the crosstalk between hosts and intracellular pathogens.

De novo synthesis of thiamine (vitamin B1) is the ancestral state in Plasmodium parasites – evidence from avian haemosporidians

Parasites often have reduced genomes as their own genes become redundant when utilizing their host as a source of metabolites, thus losing their own de novo production of metabolites. Primate malaria parasites can synthesize vitamin B1 (thiamine) de novo but rodent malaria and other genome-sequenced apicomplexans cannot, as the three essential genes responsible for this pathway are absent in their genomes. The unique presence of functional thiamine synthesis genes in primate malaria parasites and their sequence similarities to bacterial orthologues, have led to speculations that this pathway was horizontally acquired from bacteria. Here we show that the genes essential for the de novo synthesis of thiamine are found also in avian Plasmodium species. Importantly, they are also present in species phylogenetically basal to all mammalian and avian Plasmodium parasites, i.e. Haemoproteus. Furthermore, we found that these genes are expressed during the blood stage of the avian malaria infection, indicating that this metabolic pathway is actively transcribed. We conclude that the ability to synthesize thiamine is widespread among haemosporidians, with a recent loss in the rodent malaria species.

Comparative genomics of microsporidian genomes reveals a minimal non-coding RNA set and new insights for transcription in minimal eukaryotic genomes

Microsporidia are ubiquitous intracellular pathogens whose opportunistic nature led to their increased recognition with the rise of the AIDS pandemic. As the RNA world was largely unexplored in this parasitic lineage, we developed a dedicated in silico methodology to carry out exhaustive identification of ncRNAs across the Encephalitozoon and Nosema genera. Thus, the previously missing U1 small nuclear RNA (snRNA) and small nucleolar RNAs (snoRNAs) targeting only the LSU rRNA were highlighted and were further validated using 5' and 3'RACE-PCR experiments. Overall, the 15 ncRNAs that were found shared between Encephalitozoon and Nosema spp. may represent the minimal core set required for parasitic life. Interestingly, the systematic presence of a CCC- or GGG-like motif in 5' of all ncRNA and mRNA gene transcripts regardless of the RNA polymerase involved suggests that the RNA polymerase machineries in microsporidia species could use common factors. Our data provide additional insights in accordance with the simplification processes observed in these reduce genomes and underline the usefulness of sequencing closely related species to help identify highly divergent ncRNAs in these parasites.

A Novel Spore Wall Protein from Antonospora locustae (Microsporidia: Nosematidae) Contributes to Sporulation

Microsporidia are obligate intracellular parasites, existing in a wide variety of animal hosts. Here, we reported AlocSWP2, a novel protein identified from the spore wall of Antonospora locustae (formerly, Nosema locustae, and synonym, Paranosema locustae), containing four cysteines that are conserved among the homologues of several Microspodian pathogens in insects and mammals. AlocSWP2 was detected in the wall of mature spores via indirect immunofluorescence assay. In addition, immunocytochemistry localization experiments showed that the protein was observed in the wall of sporoblasts, sporonts, and meronts during sporulation within the host body, also in the wall of mature spores. AlocSWP2 was not detected in the fat body of infected locust until the 9th day after inoculating spores via RT‐PCR experiments. Furthermore, the survival percentage of infected locusts injected with dsRNA of AlocSWP2 on the 15th, 16th, and 17th days after inoculation with microsporidian were significantly higher than those of infected locusts without dsRNA treatment. Conversely, the amount of spores in locusts infected with A. locustae after treated with RNAi AlocSWP2 was significantly lower than those of infected locusts without RNAi of this gene. This novel spore wall protein from A. locustae may be involved in sporulation, thus contributing to host mortality.

Molecular Testing in Emerging Infectious Diseases

It was widely believed in the late 1960s that infectious diseases had been conquered by vaccines and antibiotics and humans were no longer under threat by microbial pathogens. Yet, since that time more than 60 pathogens have been discovered that can cause serious emerging infectious diseases. Molecular methods have played critical roles in the discovery, monitoring, and clinical diagnostics of emerging pathogens. In this chapter, we present well-recognized emerging pathogens. We provide examples of the utility of molecular assays in research and clinical care of emerging infectious diseases. We also discuss some theoretical and practical limitations of molecular tests and the future prospects of expanding molecular diagnostics for emerging pathogens based on new advances of knowledge and technologies.

The Genome of Nosema sp. Isolate YNPr: A Comparative Analysis of Genome Evolution within the Nosema/Vairimorpha Clade

The microsporidian parasite designated here as Nosema sp. Isolate YNPr was isolated from the cabbage butterfly Pieris rapae collected in Honghe Prefecture, Yunnan Province, China. The genome was sequenced by Illumina sequencing and compared to those of two related members of the Nosema/Vairimorpha clade, Nosema ceranae and Nosema apis. Based upon assembly statistics, the Nosema sp. YNPr genome is 3.36 x 10⁶bp with a G+C content of 23.18% and 2,075 protein coding sequences. An “ACCCTT” motif is present approximately 50-bp upstream of the start codon, as reported from other members of the clade and from Encephalitozoon cuniculi, a sister taxon. Comparative small subunit ribosomal DNA (SSU rDNA) analysis as well as genome-wide phylogenetic analysis confirms a closer relationship between N. ceranae and Nosema sp. YNPr than between the two honeybee parasites N. ceranae and N. apis. The more closely related N. ceranae and Nosema sp. YNPr show similarities in a number of structural characteristics such as gene synteny, gene length, gene number, transposon composition and gene reduction. Based on transposable element content of the assemblies, the transposon content of Nosema sp. YNPr is 4.8%, that of N. ceranae is 3.7%, and that of N. apis is 2.5%, with large differences in the types of transposons present among these 3 species. Gene function annotation indicates that the number of genes participating in most metabolic activities is similar in all three species. However, the number of genes in the transcription, general function, and cysteine protease categories is greater in N. apis than in the other two species. Our studies further characterize the evolution of the Nosema/Vairimorpha clade of microsporidia. These organisms maintain variable but very reduced genomes. We are interested in understanding the effects of genetic drift versus natural selection on genome size in the microsporidia and in developing a testable hypothesis for further studies on the genomic ecology of this group.

Contribution of Multiple Inter-Kingdom Horizontal Gene Transfers to Evolution and Adaptation of Amphibian-Killing Chytrid, Batrachochytrium dendrobatidis

Amphibian populations are experiencing catastrophic declines driven by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Although horizontal gene transfer (HGT) facilitates the evolution and adaptation in many fungi by conferring novel function genes to the recipient fungi, inter-kingdom HGT in Bd remains largely unexplored. In this study, our investigation detects 19 bacterial genes transferred to Bd, including metallo-beta-lactamase and arsenate reductase that play important roles in the resistance to antibiotics and arsenates. Moreover, three probable HGT gene families in Bd are from plants and one gene family coding the ankyrin repeat-containing protein appears to come from oomycetes. The observed multi-copy gene families associated with HGT are probably due to the independent transfer events or gene duplications. Five HGT genes with extracellular locations may relate to infection, and some other genes may participate in a variety of metabolic pathways, and in doing so add important metabolic traits to the recipient. The evolutionary analysis indicates that all the transferred genes evolved under purifying selection, suggesting that their functions in Bd are similar to those of the donors. Collectively, our results indicate that HGT from diverse donors may be an important evolutionary driver of Bd, and improve its adaptations for infecting and colonizing host amphibians.

TE-domestication and horizontal transfer in a putative Nef-AP1mu mimic of HLA-A cytoplasmic domain re-trafficking

Genes of the major histocompatibility complex (MHC; also called HLA in human) are polymorphic elements in the genomes of sharks to humans. Class-I and class-II MHC loci appear responsible for much of the genetic linkage to myriad disease states via the capacity to bind short (~8-15 a.a.) peptides of a given pathogen's proteome, or in some cases, the altered proteomes of cancerous cells, and even (in autoimmunity) certain nominal 'self' peptides (Janeway, 2004).(1) Unfortunately, little is known about how the canonical structure of the MHC-I/-II peptide-presenting gene evolved, particularly since beyond ~500 Mya (sharks) no paralogs exist.(2,3) We previously reported that HLA-A isotype alleles with the α1-helix, R65 motif, are wide-spread in phylogeny, but that the α 2-helix, H151R motif, has apparently segregated out of most species. Surprisingly, an uncharacterized orf in T. syrichta (Loc-103275158) encoded R151, but within a truncated A-23 like gene containing 5'- and 3'- footprints of the transposon (TE), tigger-1; the extant tarsier A-23 allele is totally missing exon-3 and part-of exon-4; together, suggesting TE-mediated inactivation of an intact/ancestral A-23 allele (Murray, 2015a).(4) The unique Loc-103275158 orf encodes a putative 15-exon transcript with no apparent paralogs throughout phylogeny. However, an HLA-A11 like gene in M. leucophaeus with a shortened C-terminal domain, and an HLA-A like orf in C. atys with two linked α1/α2/α3 domains, both contain a second transmembrane segment, which is conserved in Loc-103275158. Thus, we could model the putative protein with its Nef-like tail domain docked to its MHC-I like α3 domain (i.e., on the same side of a membrane). This modeled tertiary structure is strikingly similar to the solved structure of the Nef:MHC-I CD:AP1mu transporter (Jia, 2012).(5) Nef:AP1mu binds the CD of MHC-I in trafficking MHC-I away from the trans-golgi and into the endocytic pathway in HIV-1 infected cells. The CD loop of the Loc-103275158 provisional protein conserved the nominal MHC-I CD tyrosine phosphorylation site, and it has an N-terminal SH3 domain that we docked in one conformation to its internal Nef-like domain. Here, we suggest that phosphorylation of the protein's CD-loop signals an exchange between the internal Nef-like domain and a lentiviral-Nef for binding the N-terminal SH3 domain - freeing the Nef-like domain to bind MHC-I CD. Since the 5'-tigger sequence encodes part of the pseudo α1/α2 MHC-I domain, and the 3'-tigger part of the Nef-like domain, we speculate that transposition proceeded phylogenetically disparate horizontal transfers, involving adjacent 5'- and 3'- parasitic footprints, which we also found in the Loc-103275158 orf.

Laterally Transferred Gene Recruited as a Venom in Parasitoid Wasps

Parasitoid wasps use venom to manipulate the immunity and metabolism of their host insects in a variety of ways to provide resources for their offspring. Yet, how genes are recruited and evolve to perform venom functions remain open questions. A recently recognized source of eukaryotic genome innovation is lateral gene transfer (LGT). Glycoside hydrolase family 19 (GH19) chitinases are widespread in bacteria, microsporidia, and plants where they are used in nutrient acquisition or defense, but have previously not been known in metazoans. In this study, a GH19 chitinase LGT is described from the unicellular microsporidia/Rozella clade into parasitoid wasps of the superfamily Chalcidoidea, where it has become recruited as a venom protein. The GH19 chitinase is present in 15 species of chalcidoid wasps representing four families, and phylogenetic analysis indicates that it was laterally transferred near or before the origin of Chalcidoidea (∼95 Ma). The GH19 chitinase gene is highly expressed in the venom gland of at least seven species, indicating a role in the complex host manipulations performed by parasitoid wasp venom. RNAi knockdown in the model parasitoid Nasonia vitripennis reveals that-following envenomation-the GH19 chitinase induces fly hosts to upregulate genes involved in an immune response to fungi. A second, independent LGT of GH19 chitinase from microsporidia into mosquitoes was also found, also supported by phylogenetic reconstructions. Besides these two LGT events, GH19 chitinase is not found in any other sequenced animal genome, or in any fungi outside the microsporidia/Rozella clade.

Horizontally acquired genes in early-diverging pathogenic fungi enable the use of host nucleosides and nucleotides

Horizontal gene transfer (HGT) among bacteria, archaea, and viruses is widespread, but the extent of transfers from these lineages into eukaryotic organisms is contentious. Here we systematically identify hundreds of genes that were likely acquired horizontally from a variety of sources by the early-diverging fungal phyla Microsporidia and Cryptomycota. Interestingly, the Microsporidia have acquired via HGT several genes involved in nucleic acid synthesis and salvage, such as those encoding thymidine kinase (TK), cytidylate kinase, and purine nucleotide phosphorylase. We show that these HGT-derived nucleic acid synthesis genes tend to function at the interface between the metabolic networks of the host and pathogen. Thus, these genes likely play vital roles in diversifying the useable nucleic acid components available to the intracellular parasite, often through the direct capture of resources from the host. Using an in vivo viability assay, we also demonstrate that one of these genes, TK, encodes an enzyme that is capable of activating known prodrugs to their active form, which suggests a possible treatment route for microsporidiosis. We further argue that interfacial genes with well-understood activities, especially those horizontally transferred from bacteria or viruses, could provide medical treatments for microsporidian infections.

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.

The Prediction and Validation of Small CDSs Expand the Gene Repertoire of the Smallest Known Eukaryotic Genomes

The proper prediction of the gene catalogue of an organism is essential to obtain a representative snapshot of its overall lifestyle, especially when it is not amenable to culturing. Microsporidia are obligate intracellular, sometimes hard to culture, eukaryotic parasites known to infect members of every animal phylum. To date, sequencing and annotation of microsporidian genomes have revealed a poor gene complement with highly reduced gene sizes. In the present paper, we investigated whether such gene sizes may have induced biases for the methodologies used for genome annotation, with an emphasis on small coding sequence (CDS) gene prediction. Using better delineated intergenic regions from four Encephalitozoon genomes, we predicted de novo new small CDSs with sizes ranging from 78 to 255 bp (median 168) and corroborated these predictions by RACE-PCR experiments in Encephalitozoon cuniculi. Most of the newly found genes are present in other distantly related microsporidian species, suggesting their biological relevance. The present study provides a better framework for annotating microsporidian genomes and to train and evaluate new computational methods dedicated at detecting ultra-small genes in various organisms.

Microsporidia: Eukaryotic Intracellular Parasites Shaped by Gene Loss and Horizontal Gene Transfers

Microsporidia are eukaryotic parasites of many animals that appear to have adapted to an obligate intracellular lifestyle by modifying the morphology and content of their cells. Living inside other cells, they have lost many, or all, metabolic functions, resulting in genomes that are always gene poor and often very small. The minute content of microsporidian genomes led many to assume that these parasites are biochemically static and uninteresting. However, recent studies have demonstrated that these organisms can be surprisingly complex and dynamic. In this review I detail the most significant recent advances in microsporidian genomics and discuss how these have affected our understanding of many biological aspects of these peculiar eukaryotic intracellular pathogens.

Comparative Analysis of Codon Usage Bias Patterns in Microsporidian Genomes

The sub-3 Mbp genomes from microsporidian species of the Encephalitozoon genus are the smallest known among eukaryotes and paragons of genomic reduction and compaction in parasites. However, their diminutive stature is not characteristic of all Microsporidia, whose genome sizes vary by an order of magnitude. This large variability suggests that different evolutionary forces are applied on the group as a whole. In this study, we have compared the codon usage bias (CUB) between eight taxonomically distinct microsporidian genomes: Encephalitozoon intestinalis, Encephalitozoon cuniculi, Spraguea lophii, Trachipleistophora hominis, Enterocytozoon bieneusi, Nematocida parisii, Nosema bombycis and Nosema ceranae. While the CUB was found to be weak in all eight Microsporidia, nearly all (98%) of the optimal codons in S. lophii, T. hominis, E. bieneusi, N. parisii, N. bombycis and N. ceranae are fond of A/U in third position whereas most (64.6%) optimal codons in the Encephalitozoon species E. intestinalis and E. cuniculi are biased towards G/C. Although nucleotide composition biases are likely the main factor driving the CUB in Microsporidia according to correlation analyses, directed mutational pressure also likely affects the CUB as suggested by ENc-plots, correspondence and neutrality analyses. Overall, the Encephalitozoon genomes were found to be markedly different from the other microsporidians and, despite being the first sequenced representatives of this lineage, are uncharacteristic of the group as a whole. The disparities observed cannot be attributed solely to differences in host specificity and we hypothesize that other forces are at play in the lineage leading to Encephalitozoon species.

Gene transfers shaped the evolution of de novo NAD+ biosynthesis in eukaryotes

NAD⁺ is an essential molecule for life, present in each living cell. It can function as an electron carrier or cofactor in redox biochemistry and energetics, and serves as substrate to generate the secondary messenger cyclic ADP ribose and nicotinic acid adenine dinucleotide phosphate. Although de novo NAD⁺ biosynthesis is essential, different metabolic pathways exist in different eukaryotic clades. The kynurenine pathway starting with tryptophan was most likely present in the last common ancestor of all eukaryotes, and is active in fungi and animals. The aspartate pathway, detected in most photosynthetic eukaryotes, was probably acquired from the cyanobacterial endosymbiont that gave rise to chloroplasts. An evolutionary analysis of enzymes catalyzing de novo NAD⁺ biosynthesis resulted in evolutionary trees incongruent with established organismal phylogeny, indicating numerous gene transfers. Endosymbiotic gene transfers probably introduced the aspartate pathway into eukaryotes and may have distributed it among different photosynthetic clades. In addition, several horizontal gene transfers substituted eukaryotic genes with bacterial orthologs. Although horizontal gene transfer is accepted as a key mechanism in prokaryotic evolution, it is supposed to be rare in eukaryotic evolution. The essential metabolic pathway of de novo NAD⁺ biosynthesis in eukaryotes was shaped by numerous gene transfers.