Improvement of variables interpretability in kernel PCA

Kernel methods have been proven to be a powerful tool for the integration and analysis of highthroughput 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 KPCA representation of the data. The proposed method, kernel PCA Interpretable Gradient (KPCA-IG), provides a datadriven 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. 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

Temporal transcriptomic response during arsenic stress in Herminiimonas arsenicoxydans

Background: Arsenic is present in numerous ecosystems and microorganisms have developed various mechanisms to live in such hostile environments. Herminiimonas arsenicoxydans, a bacterium isolated from arsenic contaminated sludge, has acquired remarkable capabilities to cope with arsenic. In particular our previous studies have suggested the existence of a temporal induction of arsenite oxidase, a key enzyme in arsenic metabolism, in the presence of As(III). Results: Microarrays were designed to compare gene transcription profiles under a temporal As(III) exposure. Transcriptome kinetic analysis demonstrated the existence of two phases in arsenic response. The expression of approximatively 14% of the whole genome was significantly affected by an As(III) early stress and 4% by an As(III) late exposure. The early response was characterized by arsenic resistance, oxidative stress, chaperone synthesis and sulfur metabolism. The late response was characterized by arsenic metabolism and associated mechanisms such as phosphate transport and motility. The major metabolic changes were confirmed by chemical, transcriptional, physiological and biochemical experiments. These early and late responses were defined as general stress response and specific response to As(III), respectively. Conclusion: Gene expression patterns suggest that the exposure to As(III) induces an acute response to rapidly minimize the immediate effects of As(III). Upon a longer arsenic exposure, a broad metabolic response was induced. These data allowed to propose for the first time a kinetic model of the As(III) response in bacteria

Transcriptome analysis of antigenic variation in Plasmodium falciparum - var silencing is not dependent on antisense RNA

BACKGROUND: Plasmodium falciparum, the causative agent of the most severe form of malaria, undergoes antigenic variation through successive presentation of a family of antigens on the surface of parasitized erythrocytes. These antigens, known as Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) proteins, are subject to a mutually exclusive expression system, and are encoded by the multigene var family. The mechanism whereby inactive var genes are silenced is poorly understood. To investigate transcriptional features of this mechanism, we conducted a microarray analysis of parasites that were selected to express different var genes by adhesion to chondroitin sulfate A (CSA) or CD36. RESULTS: In addition to oligonucleotides for all predicted protein-coding genes, oligonucleotide probes specific to each known var gene of the FCR3 background were designed and added to the microarray, as well as tiled sense and antisense probes for a subset of var genes. In parasites selected for adhesion to CSA, one full-length var gene (var2csa) was strongly upregulated, as were sense RNA molecules emanating from the 3' end of a limited subset of other var genes. No global relationship between sense and antisense production of var genes was observed, but notably, some var genes had coincident high levels of both antisense and sense transcript. CONCLUSION: Mutually exclusive expression of PfEMP1 proteins results from transcriptional silencing of non-expressed var genes. The distribution of steady-state sense and antisense RNA at var loci are not consistent with a silencing mechanism based on antisense silencing of inactive var genes. Silencing of var loci is also associated with altered regulation of genes distal to var loci

Growth of Yersinia pseudotuberculosis in human plasma: impacts on virulence and metabolic gene expression

<p>Abstract</p> <p>Background</p> <p>In man, infection by the Gram-negative enteropathogen <it>Yersinia pseudotuberculosis </it>is usually limited to the terminal ileum. However, in immunocompromised patients, the microorganism may disseminate from the digestive tract and thus cause a systemic infection with septicemia.</p> <p>Results</p> <p>To gain insight into the metabolic pathways and virulence factors expressed by the bacterium at the blood stage of pseudotuberculosis, we compared the overall gene transcription patterns (the transcriptome) of bacterial cells cultured in either human plasma or Luria-Bertani medium. The most marked plasma-triggered metabolic consequence in <it>Y. pseudotuberculosis </it>was the switch to high glucose consumption, which is reminiscent of the acetogenic pathway (known as "glucose overflow") in <it>Escherichia coli</it>. However, upregulation of the glyoxylate shunt enzymes suggests that (in contrast to <it>E. coli</it>) acetate may be further metabolized in <it>Y. pseudotuberculosis</it>. Our data also indicate that the bloodstream environment can regulate major virulence genes (positively or negatively); the <it>yadA </it>adhesin gene and most of the transcriptional units of the pYV-encoded type III secretion apparatus were found to be upregulated, whereas transcription of the pH6 antigen locus was strongly repressed.</p> <p>Conclusion</p> <p>Our results suggest that plasma growth of <it>Y. pseudotuberculosis </it>is responsible for major transcriptional regulatory events and prompts key metabolic reorientations within the bacterium, which may in turn have an impact on virulence.</p

Endoplasmic Reticulum Stress-Sensing Mechanism Is Activated in Entamoeba histolytica upon Treatment with Nitric Oxide

The Endoplasmic Reticulum stores calcium and is a site of protein synthesis and modification. Changes in ER homeostasis lead to stress responses with an activation of the unfolded protein response (UPR). The Entamoeba histolytica endomembrane system is simple compared to those of higher eukaryotes, as a canonical ER is not observed. During amoebiasis, an infection of the human intestine and liver by E. histolytica, nitric oxide (NO) triggers an apoptotic-like event preceded by an impairment of energy production and a loss of important parasite pathogenic features. We address the question of how this ancient eukaryote responds to stress induced by immune components (i.e. NO) and whether stress leads to ER changes and subsequently to an UPR. Gene expression analysis suggested that NO triggers stress responses marked by (i) dramatic up-regulation of hsp genes although a bona fide UPR is absent; (ii) induction of DNA repair and redox gene expression and iii) up-regulation of glycolysis-related gene expression. Enzymology approaches demonstrate that NO directly inhibits glycolysis and enhance cysteine synthase activity. Using live imaging and confocal microscopy we found that NO dramatically provokes extensive ER fragmentation. ER fission in E. histolytica appears as a protective response against stress, as it has been recently proposed for neuron self-defense during neurologic disorders. Chronic ER stress is also involved in metabolic diseases including diabetes, where NO production reduces ER calcium levels and activates cell death. Our data highlighted unique cellular responses of interest to understand the mechanisms of parasite death during amoebiasis

Dynamic RNA profiling in Plasmodium falciparum synchronized blood stages exposed to lethal doses of artesunate

<p>Abstract</p> <p>Background</p> <p>Translation of the genome sequence of <it>Plasmodium sp</it>. into biologically relevant information relies on high through-put genomics technology which includes transcriptome analysis. However, few studies to date have used this powerful approach to explore transcriptome alterations of <it>P. falciparum </it>parasites exposed to antimalarial drugs.</p> <p>Results</p> <p>The rapid action of artesunate allowed us to study dynamic changes of the parasite transcriptome in synchronous parasite cultures exposed to the drug for 90 minutes and 3 hours. Developmentally regulated genes were filtered out, leaving 398 genes which presented altered transcript levels reflecting drug-exposure. Few genes related to metabolic pathways, most encoded chaperones, transporters, kinases, Zn-finger proteins, transcription activating proteins, proteins involved in proteasome degradation, in oxidative stress and in cell cycle regulation. A positive bias was observed for over-expressed genes presenting a subtelomeric location, allelic polymorphism and encoding proteins with potential export sequences, which often belonged to subtelomeric multi-gene families. This pointed to the mobilization of processes shaping the interface between the parasite and its environment. In parallel, pathways were engaged which could lead to parasite death, such as interference with purine/pyrimidine metabolism, the mitochondrial electron transport chain, proteasome-dependent protein degradation or the integrity of the food vacuole.</p> <p>Conclusion</p> <p>The high proportion of over-expressed genes encoding proteins exported from the parasite highlight the importance of extra-parasitic compartments as fields for exploration in drug research which, to date, has mostly focused on the parasite itself rather than on its intra and extra erythrocytic environment. Further work is needed to clarify which transcriptome alterations observed reflect a specific response to overcome artesunate toxicity or more general perturbations on the path to cellular death.</p

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

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

Inflammatory control in AIDS-resistant non human primates

International audienceAfrican non human primates are natural hosts of SIV. The infection is non-pathogenic despite plasma viral load levels similar to those in HIV-1 infected humans and SIVmac-infected macaques (MAC) progressing towards AIDS. The most striking difference between non-pathogenic SIV and pathogenic HIV-1/SIVmac infections is the lack of chronic T cell activation in natural hosts. In HIV and SIVmac infections, chronic T cell activation is known to drive CD4+T cell depletion. Intense research efforts are worldwide put on the search of the mechanisms that can control chronic T cell activation in HIV/SIV infections. Innate immune responses play a determinant role in the regulation of T cell activation profiles. Type I interferons (IFN-I) are part of the first-wave response of the innate immune system in viral infections. We compared the IFN-I responses between pathogenic (MAC) and non-pathogenic SIV infections (African Green monkey, AGM) at the level of blood and lymph nodes (LN) during the early and chronic stage of infection. During the acute SIVagm infection, we detected high amounts of IFN-α in the plasma of AGMs, although the mean levels at the peak were three times lower than in MAC. The microarray data revealed a rapid and strong up-regulation of type I Interferon-Stimulated Genes (ISG) in AGMs during acute SIVagm infection. ISGs denote the in vivo activity of IFN-I. Using a functional assay, we demonstrated that low IFN-α concentrations (50 times lower than the IFN-α levels in plasma at the peak) were sufficient to induce strong ISG responses in AGM and MAC cells. Surprisingly, our direct comparison of blood and LNs showed that ISG induction was broader in blood of AGMs than in MAC, while in LN, it was the contrary. Thus, in AGMs, less ISG were induced in LNs as compared to MAC already during the acute phase of infection. Moreover, our tight kinetic analysis showed that this ISG expression was efficiently controlled after day 28 post-infection in AGMs, while in MAC the ISGs expression remained uncontrolled. Finally, we identified genes that were differentially expressed between the two species and which might be involved in the discriminating responses. Altogether, this shows that AGMs are capable to mount a well coordinated and efficient regulative response to innate immune activation

From array-based hybridization of Helicobacter pylori isolates to the complete genome sequence of an isolate associated with MALT lymphoma

<p>Abstract</p> <p>Background</p> <p><it>elicobacter pylori </it>infection is associated with several gastro-duodenal inflammatory diseases of various levels of severity. To determine whether certain combinations of genetic markers can be used to predict the clinical source of the infection, we analyzed well documented and geographically homogenous clinical isolates using a comparative genomics approach.</p> <p>Results</p> <p>A set of 254 <it>H. pylori </it>genes was used to perform array-based comparative genomic hybridization among 120 French <it>H. pylori </it>strains associated with chronic gastritis (n = 33), duodenal ulcers (n = 27), intestinal metaplasia (n = 17) or gastric extra-nodal marginal zone B-cell MALT lymphoma (n = 43). Hierarchical cluster analyses of the DNA hybridization values allowed us to identify a homogeneous subpopulation of strains that clustered exclusively with <it>cag</it>PAI minus MALT lymphoma isolates. The genome sequence of B38, a representative of this MALT lymphoma strain-cluster, was completed, fully annotated, and compared with the six previously released <it>H. pylori </it>genomes (i.e. J99, 26695, HPAG1, P12, G27 and Shi470). B38 has the smallest <it>H. pylori </it>genome described thus far (1,576,758 base pairs containing 1,528 CDSs); it contains the <it>vacA</it>s2m2 allele and lacks the genes encoding the major virulence factors (absence of <it>cag</it>PAI, <it>bab</it>B, <it>bab</it>C, <it>sab</it>B, and <it>hom</it>B). Comparative genomics led to the identification of very few sequences that are unique to the B38 strain (9 intact CDSs and 7 pseudogenes). Pair-wise genomic synteny comparisons between B38 and the 6 <it>H. pylori </it>sequenced genomes revealed an almost complete co-linearity, never seen before between the genomes of strain Shi470 (a Peruvian isolate) and B38.</p> <p>Conclusion</p> <p>These isolates are deprived of the main <it>H. pylori </it>virulence factors characterized previously, but are nonetheless associated with gastric neoplasia.</p