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.

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.

Genome-Wide Analysis of Antigen 43 (Ag43) Variants: New Insights in Their Diversity, Distribution and Prevalence in Bacteria.

Antigen 43 (Ag43) expression induces aggregation and biofilm formation that has consequences for bacterial colonisation and infection. Ag43 is secreted through the Type 5 subtype "a" secretion system (T5aSS) and is a prototypical member of the family of self-associating autotransporters (SAATs). As a T5aSS protein, Ag43 has a modular architecture comprised of (i) a signal peptide, (ii) a passenger domain that can be subdivided into three subdomains (SL, EJ, and BL), (iii) an autochaperone (AC) domain, and (iv) an outer membrane translocator. The cell-surface SL subdomain is directly involved in the "Velcro-handshake" mechanism resulting in bacterial autoaggregation. Ag43 is considered to have a ubiquitous distribution in E. coli genomes and many strains harbour multiple agn43 genes. However, recent phylogenetic analyses indicated the existence of four distinct Ag43 classes exhibiting different propensities for autoaggregation and interactions. Given the knowledge of the diversity and distribution of Ag43 in E. coli genomes is incomplete, we have performed a thorough in silico investigation across bacterial genomes. Our comprehensive analyses indicate that Ag43 passenger domains cluster in six phylogenetic classes associated with different SL subdomains. The diversity of Ag43 passenger domains is a result of the association of the SL subtypes with two different EJ-BL-AC modules. We reveal that agn43 is almost exclusively present among bacterial species of the Enterobacteriaceae family and essentially in the Escherichia genus (99.6%) but that it is not ubiquitous in E. coli. The gene is typically present as a single copy but up to five copies of agn43 with different combinations of classes can be observed. The presence of agn43 as well as its different classes appeared to differ between Escherichia phylogroups. Strikingly, agn43 is present in 90% of E. coli from E phylogroup. Our results shed light on Ag43 diversity and provide a rational framework for investigating its role in E. coli ecophysiology and physiopathology.

Characterization of an easy-to-use method for the routine analysis of the central metabolism using an affordable low-resolution GC-MS system: application to Arthrospira platensis.

We developed an easy-to-use method for the routine analysis of the central metabolism using an affordable low-resolution GC-MS system run in SIM mode. The profiling approach was optimized for the derivatization protocol of some 60 targeted metabolites. The performance of two silylation reagents (MSTFA and BSTFA) that allowed the comprehensive derivatization of 42 key intermediary metabolites of the 60 initially targeted (organic acids, phosphate derivatives, monosaccharides and amino acids) was measured. The experimental results unequivocally showed that the MSTFA reagent met mandatory criteria including ease of handling (a very simple one-step protocol was developed), comprehensiveness of derivatization (the 42 compounds covered the extended metabolic pathways of the central carbon metabolism, with a coverage percentage ranging from 17% for the worst to 90% for the best result), optimized response coefficient of the whole derivatives (median value greater than the others by one order of magnitude) and repeatability of the protocol (RSD value below 25% for the whole procedure). When tested in real conditions (cyanobacteria polar extract), the experimental results showed that the profiling methodology was adequately repeatable (RSD = 35%) to ensure quantification results comparable with much more sensitive analytical techniques (capillary electrophoresis/mass spectrometry and liquid chromatography/triple quadrupole mass spectrometry system), while needing only about twice the quantity of biomass. Graphical abstract Schematic overview of an easy-to-use profiling method for the routine analysis of the central metabolism using a low-resolution GC-MS system.

Sequence capture by hybridization to explore modern and ancient genomic diversity in model and nonmodel organisms.

The recent expansion of next-generation sequencing has significantly improved biological research. Nevertheless, deep exploration of genomes or metagenomic samples remains difficult because of the sequencing depth and the associated costs required. Therefore, different partitioning strategies have been developed to sequence informative subsets of studied genomes. Among these strategies, hybridization capture has proven to be an innovative and efficient tool for targeting and enriching specific biomarkers in complex DNA mixtures. It has been successfully applied in numerous areas of biology, such as exome resequencing for the identification of mutations underlying Mendelian or complex diseases and cancers, and its usefulness has been demonstrated in the agronomic field through the linking of genetic variants to agricultural phenotypic traits of interest. Moreover, hybridization capture has provided access to underexplored, but relevant fractions of genomes through its ability to enrich defined targets and their flanking regions. Finally, on the basis of restricted genomic information, this method has also allowed the expansion of knowledge of nonreference species and ancient genomes and provided a better understanding of metagenomic samples. In this review, we present the major advances and discoveries permitted by hybridization capture and highlight the potency of this approach in all areas of biology.

Comparative genomics highlights the unique biology of Methanomassiliicoccales, a Thermoplasmatales-related seventh order of methanogenic archaea that encodes pyrrolysine.

BACKGROUND: A seventh order of methanogens, the Methanomassiliicoccales, has been identified in diverse anaerobic environments including the gastrointestinal tracts (GIT) of humans and other animals and may contribute significantly to methane emission and global warming. Methanomassiliicoccales are phylogenetically distant from all other orders of methanogens and belong to a large evolutionary branch composed by lineages of non-methanogenic archaea such as Thermoplasmatales, the Deep Hydrothermal Vent Euryarchaeota-2 (DHVE-2, Aciduliprofundum boonei) and the Marine Group-II (MG-II). To better understand this new order and its relationship to other archaea, we manually curated and extensively compared the genome sequences of three Methanomassiliicoccales representatives derived from human GIT microbiota, "Candidatus Methanomethylophilus alvus", "Candidatus Methanomassiliicoccus intestinalis" and Methanomassiliicoccus luminyensis. RESULTS: Comparative analyses revealed atypical features, such as the scattering of the ribosomal RNA genes in the genome and the absence of eukaryotic-like histone gene otherwise present in most of Euryarchaeota genomes. Previously identified in Thermoplasmatales genomes, these features are presently extended to several completely sequenced genomes of this large evolutionary branch, including MG-II and DHVE2. The three Methanomassiliicoccales genomes share a unique composition of genes involved in energy conservation suggesting an original combination of two main energy conservation processes previously described in other methanogens. They also display substantial differences with each other, such as their codon usage, the nature and origin of their CRISPRs systems and the genes possibly involved in particular environmental adaptations. The genome of M. luminyensis encodes several features to thrive in soil and sediment conditions suggesting its larger environmental distribution than GIT. Conversely, "Ca. M. alvus" and "Ca. M. intestinalis" do not present these features and could be more restricted and specialized on GIT. Prediction of the amber codon usage, either as a termination signal of translation or coding for pyrrolysine revealed contrasted patterns among the three genomes and suggests a different handling of the Pyl-encoding capacity. CONCLUSIONS: This study represents the first insights into the genomic organization and metabolic traits of the seventh order of methanogens. It suggests contrasted evolutionary history among the three analyzed Methanomassiliicoccales representatives and provides information on conserved characteristics among the overall methanogens and among Thermoplasmata.

Large scale explorative oligonucleotide probe selection for thousands of genetic groups on a computing grid: application to phylogenetic probe design using a curated small subunit ribosomal RNA gene database.

Phylogenetic Oligonucleotide Arrays (POAs) were recently adapted for studying the huge microbial communities in a flexible and easy-to-use way. POA coupled with the use of explorative probes to detect the unknown part is now one of the most powerful approaches for a better understanding of microbial community functioning. However, the selection of probes remains a very difficult task. The rapid growth of environmental databases has led to an exponential increase of data to be managed for an efficient design. Consequently, the use of high performance computing facilities is mandatory. In this paper, we present an efficient parallelization method to select known and explorative oligonucleotide probes at large scale using computing grids. We implemented a software that generates and monitors thousands of jobs over the European Computing Grid Infrastructure (EGI). We also developed a new algorithm for the construction of a high-quality curated phylogenetic database to avoid erroneous design due to bad sequence affiliation. We present here the performance and statistics of our method on real biological datasets based on a phylogenetic prokaryotic database at the genus level and a complete design of about 20,000 probes for 2,069 genera of prokaryotes.

KASpOD--a web service for highly specific and explorative oligonucleotide design.

SUMMARY: KASpOD is a web service dedicated to the design of signature sequences using a k-mer-based algorithm. Such highly specific and explorative oligonucleotides are then suitable for various goals, including Phylogenetic Oligonucleotide Arrays. AVAILABILITY: http://g2im.u-clermont1.fr/kaspod. CONTACT: eric.peyretaillade@udamail.fr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The intracellular parasite Anncaliia algerae induces a massive miRNA down-regulation in human cells.

Anncaliia algerae belongs to microsporidia, a group of obligate intracellular parasites related to fungi. These parasites are largely spread in water and food-webs and can infect a wide variety of hosts ranging from invertebrates to vertebrates including humans. In humans, microsporidian infections are mainly opportunistic as immunocompetent hosts can clear parasites naturally. Recent studies however have reported persistent microsporidian infections and have highlighted them as a risk factor in colon cancer. This may be a direct result of cell infection or may be an indirect effect of the infectious microenvironment and the host's response. In both cases, this raises the question of the effects of microsporidian infection at the host and host-cell levels. We aimed to address the question of human host intracellular response to microsporidian infection through a transcriptomic kinetic study of human foreskin fibroblasts (HFF) infected with A.algerae, a human infecting microsporidia with an exceptionally wide host range. We focused solely on host response studying both coding and small non-coding miRNA expression. Our study revealed a generalized down-regulation of cell miRNAs throughout infection with up to 547 different miRNAs downregulated at some timepoints and also transcriptomic dysregulations that could facilitate parasite development with immune and lipid metabolism genes modulation. We also hypothesize possible small nucleic acid expropriation explaining the miRNA downregulation. This work contributes to a better understanding of the dialogue that can occur between an intracellular parasite and its host at the cellular level, and can guide future studies on microsporidian infection biology to unravel the mode of action of these minimalist parasites at the tissue or host levels.We have also generated a kinetic and comprehensive transcriptomic data set of an infectious process that can help support comparative studies in the broader field of parasitology. Lastly, these results may warrant for caution regarding microsporidian exposure and persistent infections.

Identification and localization of polar tube proteins in the extruded polar tube of the microsporidian Anncaliia algerae.

Microsporidia are obligate intracellular parasites able to infect a wide range of hosts from invertebrates to vertebrates. The success of their invasion process is based on an original organelle, the polar tube, which is suddenly extruded from the spore to inoculate the sporoplasm into the host cytoplasm. The polar tube is mainly composed of proteins named polar tube proteins (PTPs). A comparative analysis allowed us to identify genes coding for 5 PTPs (PTP1 to PTP5) in the genome of the microsporidian Anncaliia algerae. While PTP1 and PTP2 are found on the whole polar tube, PTP3 is present in a large part of the extruded polar tube except at its end-terminal part. On the contrary, PTP4 is specifically detected at the end-terminal part of the polar tube. To complete PTPs repertoire, sequential sporal protein extractions were done with high concentration of reducing agents. In addition, a method to purify polar tubes was developed. Mass spectrometry analysis conducted on both samples led to the identification of a PTP3-like protein (PTP3b), and a new PTP (PTP7) only found at the extremity of the polar tube. The specific localization of PTPs asks the question of their roles in cell invasion processes used by A. algerae.

Genome sequence of "Candidatus Methanomethylophilus alvus" Mx1201, a methanogenic archaeon from the human gut belonging to a seventh order of methanogens.

We report the draft genome sequence of "Candidatus Methanomethylophilus alvus" Mx1201, a methanogen present in the human gut. It was enriched from human feces under anaerobic conditions with methanol as the substrate. Its circular genome, of around 1.7 Mb, contains genes needed for methylotrophic methanogenesis from methanol and tri-, di-, and monomethylamine.

Detecting unknown sequences with DNA microarrays: explorative probe design strategies.

Designing environmental DNA microarrays that can be used to survey the extreme diversity of microorganisms existing in nature, represents a stimulating challenge in the field of molecular ecology. Indeed, recent efforts in metagenomics have produced a substantial amount of sequence information from various ecosystems, and will continue to accumulate large amounts of sequence data given the qualitative and quantitative improvements in the next-generation sequencing methods. It is now possible to take advantage of these data to develop comprehensive microarrays by using explorative probe design strategies. Such strategies anticipate genetic variations and thus are able to detect known and unknown sequences in environmental samples. In this review, we provide a detailed overview of the probe design strategies currently available to construct both phylogenetic and functional DNA microarrays, with emphasis on those permitting the selection of such explorative probes. Furthermore, exploration of complex environments requires particular attention on probe sensitivity and specificity criteria. Finally, these innovative probe design approaches require exploiting newly available high-density microarray formats.

HiSpOD: probe design for functional DNA microarrays.

MOTIVATION: The use of DNA microarrays allows the monitoring of the extreme microbial diversity encountered in complex samples like environmental ones as well as that of their functional capacities. However, no probe design software currently available is adapted to easily design efficient and explorative probes for functional gene arrays. RESULTS: We present a new efficient functional microarray probe design algorithm called HiSpOD (High Specific Oligo Design). This uses individual nucleic sequences or consensus sequences produced by multiple alignments to design highly specific probes. Indeed, to bypass crucial problem of cross-hybridizations, probe specificity is assessed by similarity search against a large formatted database dedicated to microbial communities containing about 10 million coding sequences (CDS). For experimental validation, a microarray targeting genes encoding enzymes involved in chlorinated solvent biodegradation was built. The results obtained from a contaminated environmental sample proved the specificity and the sensitivity of probes designed with the HiSpOD program. AVAILABILITY: http://fc.isima.fr/~g2im/hispod/.

Bacterial community composition of biological degreasing systems and health risk assessment for workers.

Biological degreasing system is a new technology based on the degradation capabilities of microorganisms to remove oil, grease, or lubricants from metal parts. No data is available about the potential biological health hazards in such system. Thus, a health risk assessment linked to the bacterial populations present in this new degreasing technology is, therefore, necessary for workers. We performed both cultural and molecular approaches in several biological degreasing systems for various industrial contexts to investigate the composition and dynamics of bacterial populations. These biological degreasing systems did not work with the original bacterial populations. Indeed, they were colonized by a defined and restricted group of bacteria. This group replaced the indigenous bacterial populations known for degrading complex substrates. Klebsiella pneumoniae, Klebsiella oxytoca, Pseudomonas aeruginosa, and Pantoea agglomerans were important members of the microflora found in most of the biological degreasing systems. These bacteria might represent a potential health hazard for workers.

Detecting variants with Metabolic Design, a new software tool to design probes for explorative functional DNA microarray development.

BACKGROUND: Microorganisms display vast diversity, and each one has its own set of genes, cell components and metabolic reactions. To assess their huge unexploited metabolic potential in different ecosystems, we need high throughput tools, such as functional microarrays, that allow the simultaneous analysis of thousands of genes. However, most classical functional microarrays use specific probes that monitor only known sequences, and so fail to cover the full microbial gene diversity present in complex environments. We have thus developed an algorithm, implemented in the user-friendly program Metabolic Design, to design efficient explorative probes. RESULTS: First we have validated our approach by studying eight enzymes involved in the degradation of polycyclic aromatic hydrocarbons from the model strain Sphingomonas paucimobilis sp. EPA505 using a designed microarray of 8,048 probes. As expected, microarray assays identified the targeted set of genes induced during biodegradation kinetics experiments with various pollutants. We have then confirmed the identity of these new genes by sequencing, and corroborated the quantitative discrimination of our microarray by quantitative real-time PCR. Finally, we have assessed metabolic capacities of microbial communities in soil contaminated with aromatic hydrocarbons. Results show that our probe design (sensitivity and explorative quality) can be used to study a complex environment efficiently. CONCLUSIONS: We successfully use our microarray to detect gene expression encoding enzymes involved in polycyclic aromatic hydrocarbon degradation for the model strain. In addition, DNA microarray experiments performed on soil polluted by organic pollutants without prior sequence assumptions demonstrate high specificity and sensitivity for gene detection. Metabolic Design is thus a powerful, efficient tool that can be used to design explorative probes and monitor metabolic pathways in complex environments, and it may also be used to study any group of genes. The Metabolic Design software is freely available from the authors and can be downloaded and modified under general public license.

Identification of transcriptional signals in Encephalitozoon cuniculi widespread among Microsporidia phylum: support for accurate structural genome annotation.

BACKGROUND: Microsporidia are obligate intracellular eukaryotic parasites with genomes ranging in size from 2.3 Mbp to more than 20 Mbp. The extremely small (2.9 Mbp) and highly compact (approximately 1 gene/kb) genome of the human parasite Encephalitozoon cuniculi has been fully sequenced. The aim of this study was to characterize noncoding motifs that could be involved in regulation of gene expression in E. cuniculi and to show whether these motifs are conserved among the phylum Microsporidia. RESULTS: To identify such signals, 5' and 3'RACE-PCR experiments were performed on different E. cuniculi mRNAs. This analysis confirmed that transcription overrun occurs in E. cuniculi and may result from stochastic recognition of the AAUAAA polyadenylation signal. Such experiments also showed highly reduced 5'UTR's (<7 nts). Most of the E. cuniculi genes presented a CCC-like motif immediately upstream from the coding start. To characterize other signals involved in differential transcriptional regulation, we then focused our attention on the gene family coding for ribosomal proteins. An AAATTT-like signal was identified upstream from the CCC-like motif. In rare cases the cytosine triplet was shown to be substituted by a GGG-like motif. Comparative genomic studies confirmed that these different signals are also located upstream from genes encoding ribosomal proteins in other microsporidian species including Antonospora locustae, Enterocytozoon bieneusi, Anncaliia algerae (syn. Brachiola algerae) and Nosema ceranae. Based on these results a systematic analysis of the approximately 2000 E. cuniculi coding DNA sequences was then performed and brings to highlight that 364 translation initiation codons (18.29% of total CDSs) had been badly predicted. CONCLUSION: We identified various signals involved in the maturation of E. cuniculi mRNAs. Presence of such signals, in phylogenetically distant microsporidian species, suggests that a common regulatory mechanism exists among the microsporidia. Furthermore, 5'UTRs being strongly reduced, these signals can be used to ensure the accurate prediction of translation initiation codons for microsporidian genes and to improve microsporidian genome annotation.

Draft Genome Sequence of Tubulinosema ratisbonensis, a Microsporidian Species Infecting the Model Organism Drosophila melanogaster.

We present the draft genome sequence of Tubulinosema ratisbonensis, a microsporidium species infecting Drosophila melanogaster A total of 3,013 protein-encoding genes and an array of transposable elements were identified. This work represents a necessary step to develop a novel model of host-parasite relationships using the highly tractable genetic model D. melanogaster.

PhylArray: phylogenetic probe design algorithm for microarray.

MOTIVATION: Microbial diversity is still largely unknown in most environments, such as soils. In order to get access to this microbial 'black-box', the development of powerful tools such as microarrays are necessary. However, the reliability of this approach relies on probe efficiency, in particular sensitivity, specificity and explorative power, in order to obtain an image of the microbial communities that is close to reality. RESULTS: We propose a new probe design algorithm that is able to select microarray probes targeting SSU rRNA at any phylogenetic level. This original approach, implemented in a program called 'PhylArray', designs a combination of degenerate and non-degenerate probes for each target taxon. Comparative experimental evaluations indicate that probes designed with PhylArray yield a higher sensitivity and specificity than those designed by conventional approaches. Applying the combined PhyArray/GoArrays strategy helps to optimize the hybridization performance of short probes. Finally, hybridizations with environmental targets have shown that the use of the PhylArray strategy can draw attention to even previously unknown bacteria.

In vitro propagation of the microsporidian pathogen Brachiola algerae and studies of its chromosome and ribosomal DNA organization in the context of the complete genome sequencing project.

Brachiola algerae has a broad host spectrum from human to mosquitoes. The successful infection of two mosquito cell lines (Mos55: embryonic cells and Sua 4.0: hemocyte-like cells) and a human cell line (HFF) highlights the efficient adaptive capacity of this microsporidian pathogen. The molecular karyotype of this microsporidian species was determined in the context of the B. algerae genome sequencing project, showing that its haploid genome consists of 30 chromosomal-sized DNAs ranging from 160 to 2240 kbp giving an estimated genome size of 23 Mbp. A contig of 12,269 bp including the DNA sequence of the B. algerae ribosomal transcription unit has been built from initial genomic sequences and the secondary structure of the large subunit rRNA constructed. The data obtained indicate that B. algerae should be an excellent parasitic model to understand genome evolution in relation to infectious capacity.

ASaiM: a Galaxy-based framework to analyze microbiota data.

Background: New generations of sequencing platforms coupled to numerous bioinformatics tools have led to rapid technological progress in metagenomics and metatranscriptomics to investigate complex microorganism communities. Nevertheless, a combination of different bioinformatic tools remains necessary to draw conclusions out of microbiota studies. Modular and user-friendly tools would greatly improve such studies. Findings: We therefore developed ASaiM, an Open-Source Galaxy-based framework dedicated to microbiota data analyses. ASaiM provides an extensive collection of tools to assemble, extract, explore, and visualize microbiota information from raw metataxonomic, metagenomic, or metatranscriptomic sequences. To guide the analyses, several customizable workflows are included and are supported by tutorials and Galaxy interactive tours, which guide users through the analyses step by step. ASaiM is implemented as a Galaxy Docker flavour. It is scalable to thousands of datasets but also can be used on a normal PC. The associated source code is available under Apache 2 license at https://github.com/ASaiM/framework and documentation can be found online (http://asaim.readthedocs.io). Conclusions: Based on the Galaxy framework, ASaiM offers a sophisticated environment with a variety of tools, workflows, documentation, and training to scientists working on complex microorganism communities. It makes analysis and exploration analyses of microbiota data easy, quick, transparent, reproducible, and shareable.