RNA Sequencing Reveals Widespread Transcription of Natural Antisense RNAs in Entamoeba Species.

Entamoeba is a genus of Amoebozoa that includes the intestine-colonizing pathogenic species Entamoeba histolytica. To understand the basis of gene regulation in E. histolytica from an evolutionary perspective, we have profiled the transcriptomes of its closely related species E. dispar, E. moshkovskii and E. invadens. Genome-wide identification of transcription start sites (TSS) and polyadenylation sites (PAS) revealed the similarities and differences of their gene regulatory sequences. In particular, we found the widespread initiation of antisense transcription from within the gene coding sequences is a common feature among all Entamoeba species. Interestingly, we observed the enrichment of antisense transcription in genes involved in several processes that are common to species infecting the human intestine, e.g., the metabolism of phospholipids. These results suggest a potentially conserved and compact gene regulatory system in Entamoeba

Quantification of stochastic noise of splicing and polyadenylation in Entamoeba histolytica

Alternative splicing and polyadenylation were observed pervasively in eukaryotic messenger RNAs. These alternative isoforms could either be consequences of physiological regulation or stochastic noise of RNA processing. To quantify the extent of stochastic noise in splicing and polyadenylation, we analyzed the alternative usage of splicing and polyadenylation sites in Entamoeba histolytica using RNA-Seq. First, we identified a large number of rarely spliced alternative junctions and then showed that the occurrence of these alternative splicing events is correlated with splicing site sequence, occurrence of constitutive splicing events and messenger RNA abundance. Our results implied the majority of these alternative splicing events are likely to be stochastic error of splicing machineries, and we estimated the corresponding error rates. Second, we observed extensive microheterogeneity of polyadenylation cleavage sites, and the extent of such microheterogeneity is correlated with the occurrence of constitutive cleavage events, suggesting most of such microheterogeneity is likely to be stochastic. Overall, we only observed a small fraction of alternative splicing and polyadenylation isoforms that are unlikely to be solely stochastic, implying the functional relevance of alternative splicing and polyadenylation in E. histolytica is limited. Lastly, we revised the gene models and annotated their 3′UTR in AmoebaDB, providing valuable resources to the community

Supervised Multiple Kernel Learning approaches for multi-omics data integration

Advances in high-throughput technologies have originated an ever-increasing availability of omics datasets. The integration of multiple heterogeneous data sources is currently an issue for biology and bioinformatics. Multiple kernel learning (MKL) has shown to be a flexible and valid approach to consider the diverse nature of multi-omics inputs, despite being an underused tool in genomic data mining.We provide novel MKL approaches based on different kernel fusion strategies.To learn from the meta-kernel of input kernels, we adaptedunsupervised integration algorithms for supervised tasks with support vector machines.We also tested deep learning architectures for kernel fusion and classification.The results show that MKL-based models can compete with more complex, state-of-the-art, supervised multi-omics integrative approaches. Multiple kernel learning offers a natural framework for predictive models in multi-omics genomic data. Our results offer a direction for bio-data mining research and further development of methods for heterogeneous data integration

Genomic determinants for initiation and length of natural antisense transcripts in Entamoeba histolytica

International audienceNatural antisense transcripts (NAT) have been reported in prokaryotes and eukaryotes. While the functions of most reported NATs remain unknown, their potentials in regulating the transcription of their counterparts have been speculated. Entamoeba histolytica, which is a unicellular eukaryotic parasite, has a compact protein-coding genome with very short intronic and intergenic regions. The regulatory mechanisms of gene expression in this compact genome are under-described. In this study, by genome-wide mapping of RNA-Seq data in the genome of E. histolytica, we show that a substantial fraction of its protein-coding genes (28%) has significant transcription on their opposite strand (i.e. NAT). Intriguingly, we found the location of transcription start sites or polyadenylation sites of NAT are determined by the specific motifs encoded on the opposite strand of the gene coding sequences, thereby providing a compact regulatory system for gene transcription. Moreover, we demonstrated that NATs are globally up-regulated under various environmental conditions including temperature stress and pathogenicity. While NATs do not appear to be consequences of spurious transcription, they may play a role in regulating gene expression in E. histolytica, a hypothesis which needs to be tested

Extensive transcriptome analysis correlates the plasticity of Entamoeba histolytica pathogenesis to rapid phenotype changes depending on the environment

International audienceAmoebiasis is a human infectious disease due to the amoeba parasite Entamoeba histolytica. The disease appears in only 20% of the infections. Diversity in phenotypes may occur within the same infectious strain in the gut; for instance, parasites can be commensal (in the intestinal lumen) or pathogenic (inside the tissue). The degree of pathogenesis of clinical isolates varies greatly. These findings raise the hypothesis that genetic derivation may account for amoebic diverse phenotypes. The main goal of this study was to analyse gene expression changes of a single virulent amoebic strain in different environmental contexts where it exhibit different degrees of virulence, namely isolated from humans and maintained through animal liver passages, in contact with the human colon and short or prolonged in vitro culture. The study reveals major transcriptome changes in virulent parasites upon contact with human colon explants, including genes related to sugar metabolism, cytoskeleton rearrangement, stress responses and DNA repair. Furthermore, in long-term cultured parasites, drastic changes in gene expression for proteins with functions for proteasome and tRNA activities were found. Globally we conclude that rapid changes in gene expression rather than genetic derivation can sustain the invasive phenotype of a single virulent isolate of E. histolytica

GH16 and GH81 family β-(1,3)-glucanases in Aspergillus fumigatus are essential for conidial cell wall morphogenesis.

International audienceThe fungal cell wall is a rigid structure because of fibrillar and branched β-(1,3)-glucan linked to chitin. Softening of the cell wall is an essential phenomenon during fungal morphogenesis, wherein rigid cell wall structures are cleaved by glycosylhydrolases. During the search for glycosylhydrolases acting on β-(1,3)-glucan, we identified seven genes in the Aspergillus fumigatus genome coding for potential endo-β-(1,3)-glucanase. ENG1 (previously characterized and named ENGL1, Mouyna et al., ), belongs to the Glycoside-Hydrolase 81 (GH81) family, while ENG2 to ENG7, to GH16 family. ENG1 and four GH16 genes (ENG2-5) were expressed in the resting conidia as well as during germination, suggesting an essential role during A. fumigatus morphogenesis. Here, we report the effect of sequential deletion of AfENG2-5 (GH16) followed by AfENG1 (GH81) deletion in the Δeng2,3,4,5 mutant. The Δeng1,2,3,4,5 mutant showed conidial defects, with linear chains of conidia unable to separate while the germination rate was not affected. These results show, for the first time in a filamentous fungus, that endo β-(1,3)-glucanases are essential for proper conidial cell wall assembly and thus segregation of conidia during conidiation

Comprehensive evaluation of normalization methods for High-throughput rnaSequencing Data analysis

Comprehensive evaluation of normalization methods for High-throughput rnaSequencing Data analysis. 4 Workshop of StatSeq COST Action TD080

Impact of the gut microbiota on the m 6 A epitranscriptome of mouse cecum and liver

International audienceThe intestinal microbiota modulates host physiology and gene expression via mechanisms that are not fully understood. Here we examine whether host epitranscriptomic marks are affected by the gut microbiota. We use methylated RNA-immunoprecipitation and sequen-cing (MeRIP-seq) to identify N6-methyladenosine (m 6 A) modifications in mRNA of mice carrying conventional, modified, or no microbiota. We find that variations in the gut micro-biota correlate with m 6 A modifications in the cecum, and to a lesser extent in the liver, affecting pathways related to metabolism, inflammation and antimicrobial responses. We analyze expression levels of several known writer and eraser enzymes, and find that the methyltransferase Mettl16 is downregulated in absence of a microbiota, and one of its target mRNAs, encoding S-adenosylmethionine synthase Mat2a, is less methylated. We furthermore show that Akkermansia muciniphila and Lactobacillus plantarum affect specific m 6 A modifications in mono-associated mice. Our results highlight epitranscriptomic modifications as an additional level of interaction between commensal bacteria and their host

Massive activation of archaeal defense genes during viral infection

Archaeal viruses display unusually high genetic and morphologic diversity. Studies of these viruses proved to be instrumental for the expansion of knowledge on viral diversity and evolution. The Sulfolobus islandicus Rod Shaped Virus 2 (SIRV2) is a model to study virus-host interactions in Archaea. It is a lytic virus that exploits a unique egress mechanism based on formation of remarkable pyramidal structures on the host cell envelope. Using whole transcriptome sequencing we present here a global map defining host and viral gene expression during the infection cycle of SIRV2 in its hyperthermophilic host S. islandicus LAL14/1. This information was used, in combination with a yeast two hybrid analysis of SIRV2 protein interactions, to advance current understanding of viral gene functions. As a consequence of SIRV2 infection, transcription of more than one third of S. islandicus genes was differentially regulated. While expression decreased of genes involved in cell division, those playing a role in anti-viral defense were activated on large scale. Expression of genes belonging to Toxin-Antitoxin and Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-Cas systems was specifically pronounced. The observed different degree of activation of various CRISPR-Cas systems highlights the specialized functions they perform. The information on individual gene expression and activation of anti-viral defense systems is expected to aid future studies aiming at detailed understanding of functions and interplay of these systems in vivo.status: publishe

Differential transcriptomic response of Anopheles arabiensis to Plasmodium vivax and Plasmodium falciparum infection

Publié dans BioRxiv le 28 mai 2021Plasmodium vivax malaria is now recognized as the second most dangerous parasitic threat to human health with the regular decrease of Plasmodium falciparum worldwide over recent decades. A very limited numbers of studies address the interaction of P. vivax with its Anopheles mosquito vectors. Those studies were conducted in P. vivax endemic countries with P.vivax local major vectors for which limited genomic and genetic tools are available. Despite the presence of P. vivax in several African countries and increasing reports on its occurrence in many others, there is virtually no data on the molecular responses of Anopheles arabiensis, a major African mosquito vector, to P. vivax, which limits the development of further mosquito-targeted interventions aimed at reducing P. vivax transmission. Taking advantage of the situation of Madagascar where P. falciparum, P. vivax and An. arabiensis are present, we explore the molecular responses of An. arabiensis towards these two human malaria parasites. RNA sequencing on RNAs isolated from mosquito midguts dissected at the early stage of infection (24 hours) was performed using mosquitoes fed on the blood of P. vivax and P. falciparum gametocyte carriers in a field station. From a de novo assembly of An. arabiensis midgut total RNA transcriptome, the comparative analysis revealed that a greater number of genes were differentially expressed in the mosquito midgut in response to P. vivax (209) than to P. falciparum (81). Among these, 15 common genes were identified to be significantly expressed in mosquito midgut 24 hours after ingesting P. vivax and P. falciparum gametocytes, including immune responsive genes and genes involved in amino-acid detoxification pathways. Importantly, working with both wild mosquitoes and field circulating parasites, our analysis revealed a strong mosquito genotype by parasite genotype interaction. Our study also identified 51 putative long non-coding RNAs differentially expressed in An. arabiensis mosquito infected midgut. Among these, several mapped to the published An. arabiensis genome at genes coding immune responsive genes such as gambicin 1, leucine-rich repeat containing genes, either on sense or antisense strands. This study constitutes the first comparison of An. arabiensis molecular interaction with P. vivax and P. falciparum, investigating both coding and long non-coding RNAs for the identification of potential transcripts, that could lead to the development of novel approaches to simultaneously block the transmission of vivax and falciparum malaria