Improvement of variables interpretability in kernel PCA.

BACKGROUND: Kernel methods have been proven to be a powerful tool for the integration and analysis of high-throughput technologies generated data. Kernels offer a nonlinear version of any linear algorithm solely based on dot products. The kernelized version of principal component analysis is a valid nonlinear alternative to tackle the nonlinearity of biological sample spaces. This paper proposes a novel methodology to obtain a data-driven feature importance based on the kernel PCA representation of the data. RESULTS: The proposed method, kernel PCA Interpretable Gradient (KPCA-IG), provides a data-driven feature importance that is computationally fast and based solely on linear algebra calculations. It has been compared with existing methods on three benchmark datasets. The accuracy obtained using KPCA-IG selected features is equal to or greater than the other methods' average. Also, the computational complexity required demonstrates the high efficiency of the method. An exhaustive literature search has been conducted on the selected genes from a publicly available Hepatocellular carcinoma dataset to validate the retained features from a biological point of view. The results once again remark on the appropriateness of the computed ranking. CONCLUSIONS: The black-box nature of kernel PCA needs new methods to interpret the original features. Our proposed methodology KPCA-IG proved to be a valid alternative to select influential variables in high-dimensional high-throughput datasets, potentially unravelling new biological and medical biomarkers.

Lassa virus activates myeloid dendritic cells but suppresses their ability to stimulate T cells.

Lassa virus (LASV) is responsible for a viral hemorrhagic fever in humans and the death of 3,000 to 5,000 people every year. The immune response to LASV is poorly understood, but type I interferon (IFN-I) and T-cell responses appear to be critical for the host. We studied the response of myeloid dendritic cells (mDC) to LASV, as mDCs are involved in both IFN-I production and T-cell activation. We compared the response of primary human mDCs to LASV and Mopeia virus (MOPV), which is similar to LASV, but non-pathogenic. We showed that mDCs produced substantial amounts of IFN-I in response to both LASV and MOPV. However, only MOPV-infected mDCs were able to activate T cells. More surprisingly, coculture with T cells completely inhibited the activation of LASV-infected mDCs. These differences between LASV and MOPV were mostly due to the LASV nucleoprotein, which has major immunosuppressive properties, but the glycoprotein was also involved. Overall, these results suggest that mDCs may be important for the global response to LASV and play a role in the outcome of Lassa fever.

Tri-mannose grafting of chitosan nanocarriers remodels the macrophage response to bacterial infection.

BACKGROUND: Infectious diseases are still a leading cause of death and, with the emergence of drug resistance, pose a great threat to human health. New drugs and strategies are thus urgently needed to improve treatment efficacy and limit drug-associated side effects. Nanotechnology-based drug delivery systems are promising approaches, offering hope in the fight against drug resistant bacteria. However, how nanocarriers influence the response of innate immune cells to bacterial infection is mostly unknown. RESULTS: Here, we used Mycobacterium tuberculosis as a model of bacterial infection to examine the impact of mannose functionalization of chitosan nanocarriers (CS-NCs) on the human macrophage response. Both ungrafted and grafted CS-NCs were similarly internalized by macrophages, via an actin cytoskeleton-dependent process. Although tri-mannose ligands did not modify the capacity of CS-NCs to escape lysosomal degradation, they profoundly remodeled the response of M. tuberculosis-infected macrophages. mRNA sequencing showed nearly 900 genes to be differentially expressed due to tri-mannose grafting. Unexpectedly, the set of modulated genes was enriched for pathways involved in cell metabolism, particularly oxidative phosphorylation and sugar metabolism. CONCLUSIONS: The ability to modulate cell metabolism by grafting ligands at the surface of nanoparticles may thus be a promising strategy to reprogram immune cells and improve the efficacy of encapsulated drugs.

Bacteriocin from epidemic Listeria strains alters the host intestinal microbiota to favor infection.

Listeria monocytogenes is responsible for gastroenteritis in healthy individuals and for a severe invasive disease in immunocompromised patients. Among the three identified L. monocytogenes evolutionary lineages, lineage I strains are overrepresented in epidemic listeriosis outbreaks, but the mechanisms underlying the higher virulence potential of strains of this lineage remain elusive. Here, we demonstrate that Listeriolysin S (LLS), a virulence factor only present in a subset of lineage I strains, is a bacteriocin highly expressed in the intestine of orally infected mice that alters the host intestinal microbiota and promotes intestinal colonization by L. monocytogenes, as well as deeper organ infection. To our knowledge, these results therefore identify LLS as the first bacteriocin described in L. monocytogenes and associate modulation of host microbiota by L. monocytogenes epidemic strains to increased virulence.

Expansion of the SOS regulon of Vibrio cholerae through extensive transcriptome analysis and experimental validation.

BACKGROUND: The SOS response is an almost ubiquitous response of cells to genotoxic stresses. The full complement of genes in the SOS regulon for Vibrio species has only been addressed through bioinformatic analyses predicting LexA binding box consensus and in vitro validation. Here, we perform whole transcriptome sequencing from Vibrio cholerae treated with mitomycin C as an SOS inducer to characterize the SOS regulon and other pathways affected by this treatment. RESULTS: Comprehensive transcriptional profiling allowed us to define the full landscape of promoters and transcripts active in V. cholerae. We performed extensive transcription start site (TSS) mapping as well as detection/quantification of the coding and non-coding RNA (ncRNA) repertoire in strain N16961. To improve TSS detection, we developed a new technique to treat RNA extracted from cells grown in various conditions. This allowed for identification of 3078 TSSs with an average 5'UTR of 116 nucleotides, and peak distribution between 16 and 64 nucleotides; as well as 629 ncRNAs. Mitomycin C treatment induced transcription of 737 genes and 28 ncRNAs at least 2 fold, while it repressed 231 genes and 17 ncRNAs. Data analysis revealed that in addition to the core genes known to integrate the SOS regulon, several metabolic pathways were induced. This study allowed for expansion of the Vibrio SOS regulon, as twelve genes (ubiEJB, tatABC, smpA, cep, VC0091, VC1190, VC1369-1370) were found to be co-induced with their adjacent canonical SOS regulon gene(s), through transcriptional read-through. Characterization of UV and mitomycin C susceptibility for mutants of these newly identified SOS regulon genes and other highly induced genes and ncRNAs confirmed their role in DNA damage rescue and protection. CONCLUSIONS: We show that genotoxic stress induces a pervasive transcriptional response, affecting almost 20% of the V. cholerae genes. We also demonstrate that the SOS regulon is larger than previously known, and its syntenic organization is conserved among Vibrio species. Furthermore, this specific co-localization is found in other γ-proteobacteria for genes recN-smpA and rmuC-tatABC, suggesting SOS regulon conservation in this phylum. Finally, we comment on the limitations of widespread NGS approaches for identification of all RNA species in bacteria.

MEMHDX: an interactive tool to expedite the statistical validation and visualization of large HDX-MS datasets.

MOTIVATION: With the continued improvement of requisite mass spectrometers and UHPLC systems, Hydrogen/Deuterium eXchange Mass Spectrometry (HDX-MS) workflows are rapidly evolving towards the investigation of more challenging biological systems, including large protein complexes and membrane proteins. The analysis of such extensive systems results in very large HDX-MS datasets for which specific analysis tools are required to speed up data validation and interpretation. RESULTS: We introduce a web application and a new R-package named 'MEMHDX' to help users analyze, validate and visualize large HDX-MS datasets. MEMHDX is composed of two elements. A statistical tool aids in the validation of the results by applying a mixed-effects model for each peptide, in each experimental condition, and at each time point, taking into account the time dependency of the HDX reaction and number of independent replicates. Two adjusted P-values are generated per peptide, one for the 'Change in dynamics' and one for the 'Magnitude of ΔD', and are used to classify the data by means of a 'Logit' representation. A user-friendly interface developed with Shiny by RStudio facilitates the use of the package. This interactive tool allows the user to easily and rapidly validate, visualize and compare the relative deuterium incorporation on the amino acid sequence and 3D structure, providing both spatial and temporal information. AVAILABILITY AND IMPLEMENTATION: MEMHDX is freely available as a web tool at the project home page http://memhdx.c3bi.pasteur.fr CONTACT: marie-agnes.dillies@pasteur.fr or sebastien.brier@pasteur.frSupplementary information: Supplementary data is available at Bioinformatics online.

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

The 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.

Repressor activity of the RpoS/σS-dependent RNA polymerase requires DNA binding.

The RpoS/σ(S) sigma subunit of RNA polymerase (RNAP) activates transcription of stationary phase genes in many Gram-negative bacteria and controls adaptive functions, including stress resistance, biofilm formation and virulence. In this study, we address an important but poorly understood aspect of σ(S)-dependent control, that of a repressor. Negative regulation by σ(S) has been proposed to result largely from competition between σ(S) and other σ factors for binding to a limited amount of core RNAP (E). To assess whether σ(S) binding to E alone results in significant downregulation of gene expression by other σ factors, we characterized an rpoS mutant of Salmonella enterica serovar Typhimurium producing a σ(S) protein proficient for Eσ(S) complex formation but deficient in promoter DNA binding. Genome expression profiling and physiological assays revealed that this mutant was defective for negative regulation, indicating that gene repression by σ(S) requires its binding to DNA. Although the mechanisms of repression by σ(S) are likely specific to individual genes and environmental conditions, the study of transcription downregulation of the succinate dehydrogenase operon suggests that σ competition at the promoter DNA level plays an important role in gene repression by Eσ(S).

Bacteriophages to reduce gut carriage of antibiotic resistant uropathogens with low impact on microbiota composition.

Uropathogenic Escherichia coli (UPEC) is the leading cause of urinary tract infections (UTIs) worldwide, causing over 150 million clinical cases annually. There is currently no specific treatment addressing the asymptomatic carriage in the gut of UPEC before they initiate UTIs. This study investigates the efficacy of virulent bacteriophages to decrease carriage of gut pathogens. Three virulent bacteriophages infecting an antibiotic-resistant UPEC strain were isolated and characterized both in vitro and in vivo. A new experimental murine model of gut carriage of E. coli was elaborated and the impact of virulent bacteriophages on colonization levels and microbiota diversity was assessed. A single dose of a cocktail of the three bacteriophages led to a sharp decrease in E. coli levels throughout the gut. We also observed that microbiota diversity was much less affected by bacteriophages than by antibiotics. Therefore, virulent bacteriophages can efficiently target UPEC strains residing in the gut, with potentially profound public health and economic impacts. These results open a new area with the possibility to manipulate specifically the microbiota using virulent bacteriophages, which could have broad applications in many gut-related disorders/diseases and beyond.

Single nucleotide resolution RNA-seq uncovers new regulatory mechanisms in the opportunistic pathogen Streptococcus agalactiae.

BACKGROUND: Streptococcus agalactiae, or Group B Streptococcus, is a leading cause of neonatal infections and an increasing cause of infections in adults with underlying diseases. In an effort to reconstruct the transcriptional networks involved in S. agalactiae physiology and pathogenesis, we performed an extensive and robust characterization of its transcriptome through a combination of differential RNA-sequencing in eight different growth conditions or genetic backgrounds and strand-specific RNA-sequencing. RESULTS: Our study identified 1,210 transcription start sites (TSSs) and 655 transcript ends as well as 39 riboswitches and cis-regulatory regions, 39 cis-antisense non-coding RNAs and 47 small RNAs potentially acting in trans. Among these putative regulatory RNAs, ten were differentially expressed in response to an acid stress and two riboswitches sensed directly or indirectly the pH modification. Strikingly, 15% of the TSSs identified were associated with the incorporation of pseudo-templated nucleotides, showing that reiterative transcription is a pervasive process in S. agalactiae. In particular, 40% of the TSSs upstream genes involved in nucleotide metabolism show reiterative transcription potentially regulating gene expression, as exemplified for pyrG and thyA encoding the CTP synthase and the thymidylate synthase respectively. CONCLUSIONS: This comprehensive map of the transcriptome at the single nucleotide resolution led to the discovery of new regulatory mechanisms in S. agalactiae. It also provides the basis for in depth analyses of transcriptional networks in S. agalactiae and of the regulatory role of reiterative transcription following variations of intra-cellular nucleotide pools.

A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis.

During the last 3 years, a number of approaches for the normalization of RNA sequencing data have emerged in the literature, differing both in the type of bias adjustment and in the statistical strategy adopted. However, as data continue to accumulate, there has been no clear consensus on the appropriate normalization method to be used or the impact of a chosen method on the downstream analysis. In this work, we focus on a comprehensive comparison of seven recently proposed normalization methods for the differential analysis of RNA-seq data, with an emphasis on the use of varied real and simulated datasets involving different species and experimental designs to represent data characteristics commonly observed in practice. Based on this comparison study, we propose practical recommendations on the appropriate normalization method to be used and its impact on the differential analysis of RNA-seq data.

Identification of the virulence landscape essential for Entamoeba histolytica invasion of the human colon.

Entamoeba histolytica is the pathogenic amoeba responsible for amoebiasis, an infectious disease targeting human tissues. Amoebiasis arises when virulent trophozoites start to destroy the muco-epithelial barrier by first crossing the mucus, then killing host cells, triggering inflammation and subsequently causing dysentery. The main goal of this study was to analyse pathophysiology and gene expression changes related to virulent (i.e. HM1:IMSS) and non-virulent (i.e. Rahman) strains when they are in contact with the human colon. Transcriptome comparisons between the two strains, both in culture conditions and upon contact with human colon explants, provide a global view of gene expression changes that might contribute to the observed phenotypic differences. The most remarkable feature of the virulent phenotype resides in the up-regulation of genes implicated in carbohydrate metabolism and processing of glycosylated residues. Consequently, inhibition of gene expression by RNA interference of a glycoside hydrolase (ß-amylase absent from humans) abolishes mucus depletion and tissue invasion by HM1:IMSS. In summary, our data suggest a potential role of carbohydrate metabolism in colon invasion by virulent E. histolytica.

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.

A putative regulatory genetic locus modulates virulence in the pathogen Leptospira interrogans.

Limited research has been conducted on the role of transcriptional regulators in relation to virulence in Leptospira interrogans, the etiological agent of leptospirosis. Here, we identify an L. interrogans locus that encodes a sensor protein, an anti-sigma factor antagonist, and two genes encoding proteins of unknown function. Transposon insertion into the gene encoding the sensor protein led to dampened transcription of the other 3 genes in this locus. This lb139 insertion mutant (the lb139(-) mutant) displayed attenuated virulence in the hamster model of infection and reduced motility in vitro. Whole-transcriptome analyses using RNA sequencing revealed the downregulation of 115 genes and the upregulation of 28 genes, with an overrepresentation of gene products functioning in motility and signal transduction and numerous gene products with unknown functions, predicted to be localized to the extracellular space. Another significant finding encompassed suppressed expression of the majority of the genes previously demonstrated to be upregulated at physiological osmolarity, including the sphingomyelinase C precursor Sph2 and LigB. We provide insight into a possible requirement for transcriptional regulation as it relates to leptospiral virulence and suggest various biological processes that are affected due to the loss of native expression of this genetic locus.

Identification and characterization of the genetic changes responsible for the characteristic smooth-to-rough morphotype alterations of clinically persistent Mycobacterium abscessus.

Mycobacterium abscessus is an emerging pathogen that is increasingly recognized as a relevant cause of human lung infection in cystic fibrosis patients. This highly antibiotic-resistant mycobacterium is an exception within the rapidly growing mycobacteria, which are mainly saprophytic and non-pathogenic organisms. M. abscessus manifests as either a smooth (S) or a rough (R) colony morphotype, which is of clinical importance as R morphotypes are associated with more severe and persistent infections. To better understand the molecular mechanisms behind the S/R alterations, we analysed S and R variants of three isogenic M. abscessus S/R pairs using an unbiased approach involving genome and transcriptome analyses, transcriptional fusions and integrating constructs. This revealed different small insertions, deletions (indels) or single nucleotide polymorphisms within the non-ribosomal peptide synthase gene cluster mps1-mps2-gap or mmpl4b in the three R variants, consistent with the transcriptional differences identified within this genomic locus that is implicated in the synthesis and transport of Glyco-Peptido-Lipids (GPL). In contrast to previous reports, the identification of clearly defined genetic lesions responsible for the loss of GPL-production or transport makes a frequent switching back-and-forth between smooth and rough morphologies in M. abscessus highly unlikely, which is important for our understanding of persistent M. abscessus infections.

Massive activation of archaeal defense genes during viral infection.

Archaeal viruses display unusually high genetic and morphological 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 the 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 of genes involved in cell division decreased, those genes playing a role in antiviral defense were activated on a 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 antiviral defense systems is expected to aid future studies aimed at detailed understanding of the functions and interplay of these systems in vivo.

Eoulsan: a cloud computing-based framework facilitating high throughput sequencing analyses.

UNLABELLED: We developed a modular and scalable framework called Eoulsan, based on the Hadoop implementation of the MapReduce algorithm dedicated to high-throughput sequencing data analysis. Eoulsan allows users to easily set up a cloud computing cluster and automate the analysis of several samples at once using various software solutions available. Our tests with Amazon Web Services demonstrated that the computation cost is linear with the number of instances booked as is the running time with the increasing amounts of data. AVAILABILITY AND IMPLEMENTATION: Eoulsan is implemented in Java, supported on Linux systems and distributed under the LGPL License at: http://transcriptome.ens.fr/eoulsan/

Identification of RNA partners of viral proteins in infected cells.

RNA viruses exhibit small-sized genomes encoding few proteins, but still establish complex networks of protein-protein and RNA-protein interactions within a cell to achieve efficient replication and spreading. Deciphering these interactions is essential to reach a comprehensive understanding of the viral infection process. To study RNA-protein complexes directly in infected cells, we developed a new approach based on recombinant viruses expressing tagged viral proteins that were purified together with their specific RNA partners. High-throughput sequencing was then used to identify these RNA molecules. As a proof of principle, this method was applied to measles virus nucleoprotein (MV-N). It revealed that in addition to full-length genomes, MV-N specifically interacted with a unique population of 5' copy-back defective interfering RNA genomes that we characterized. Such RNA molecules were able to induce strong activation of interferon-stimulated response element promoter preferentially via the cytoplasmic pattern recognition receptor RIG-I protein, demonstrating their biological functionality. Thus, this method provides a new platform to explore biologically active RNA-protein networks that viruses establish within infected cells.

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.