SNP calling from RNA-seq data without a reference genome: identification, quantification, differential analysis and impact on the protein sequence.

SNPs (Single Nucleotide Polymorphisms) are genetic markers whose precise identification is a prerequisite for association studies. Methods to identify them are currently well developed for model species, but rely on the availability of a (good) reference genome, and therefore cannot be applied to non-model species. They are also mostly tailored for whole genome (re-)sequencing experiments, whereas in many cases, transcriptome sequencing can be used as a cheaper alternative which already enables to identify SNPs located in transcribed regions. In this paper, we propose a method that identifies, quantifies and annotates SNPs without any reference genome, using RNA-seq data only. Individuals can be pooled prior to sequencing, if not enough material is available from one individual. Using pooled human RNA-seq data, we clarify the precision and recall of our method and discuss them with respect to other methods which use a reference genome or an assembled transcriptome. We then validate experimentally the predictions of our method using RNA-seq data from two non-model species. The method can be used for any species to annotate SNPs and predict their impact on the protein sequence. We further enable to test for the association of the identified SNPs with a phenotype of interest.

S1PR5 is pivotal for the homeostasis of patrolling monocytes.

Patrolling Ly6C(-) monocytes are blood-circulating cells that play a role in inflammation and in the defense against pathogens. Here, we show that similar to natural killer (NK) cells, patrolling monocytes express high levels of S1PR5, a G-coupled receptor for sphingosine-1 phosphate. We found that S1pr5(-/-) mice lack peripheral Ly6C(-) monocytes but have a normal number of these cells in the bone marrow (BM). Various lines of evidence exclude a direct contribution of S1PR5 in the survival of Ly6C(-) monocytes at the periphery. Rather, our data support a role for S1PR5 in the egress of Ly6C(-) monocytes from the BM. In particular, we observed a reduced frequency of patrolling monocytes in BM sinusoids of S1PR5 KO mice. Unexpectedly, S1P was not a chemoattractant for patrolling monocytes and had no significant effect on their viability in vitro. Moreover, the disruption of S1P gradients in vivo did not alter Ly6C(-) monocyte trafficking and viability. These data suggest that S1PR5 regulates the trafficking of monocytes via a mechanism independent of S1P gradients.

Negative regulation of NKG2D expression by IL-4 in memory CD8 T cells.

IL-4 is one of the main cytokines produced during Th2-inducing pathologies. This cytokine has been shown to affect a number of immune processes such as Th differentiation and innate immune responses. However, the impact of IL-4 on CD8 T cell responses remains unclear. In this study, we analyzed the effects of IL-4 on global gene expression profiles of Ag-induced memory CD8 T cells in the mouse. Gene ontology analysis of this signature revealed that IL-4 regulated most importantly genes associated with immune responses. Moreover, this IL-4 signature overlapped with the set of genes preferentially expressed by memory CD8 T cells over naive CD8 T cells. In particular, IL-4 downregulated in vitro and in vivo in a STAT6-dependent manner the memory-specific expression of NKG2D, thereby increasing the activation threshold of memory CD8 T cells. Furthermore, IL-4 impaired activation of memory cells as well as their differentiation into effector cells. This phenomenon could have an important clinical relevance as patients affected by Th2 pathologies such as parasitic infections or atopic dermatitis often suffer from viral-induced complications possibly linked to inefficient CD8 T cell responses.

Genomic analysis of the regulatory elements and links with intrinsic DNA structural properties in the shrunken genome of Buchnera.

BACKGROUND: Buchnera aphidicola is an obligate symbiotic bacterium, associated with most of the aphididae, whose genome has drastically shrunk during intracellular evolution. Gene regulation in Buchnera has been a matter of controversy in recent years as the combination of genomic information with the experimental results has been contradictory, refuting or arguing in favour of a functional and responsive transcription regulation in Buchnera.The goal of this study was to describe the gene transcription regulation capabilities of Buchnera based on the inventory of cis- and trans-regulators encoded in the genomes of five strains from different aphids (Acyrthosiphon pisum, Schizaphis graminum, Baizongia pistacea, Cinara cedri and Cinara tujafilina), as well as on the characterisation of some intrinsic structural properties of the DNA molecule in these bacteria. RESULTS: Interaction graph analysis shows that gene neighbourhoods are conserved between E. coli and Buchnera in structures called transcriptons, interactons and metabolons, indicating that selective pressures have acted on the evolution of transcriptional, protein-protein interaction and metabolic networks in Buchnera. The transcriptional regulatory network in Buchnera is composed of a few general DNA-topological regulators (Nucleoid Associated Proteins and topoisomerases), with the quasi-absence of any specific ones (except for multifunctional enzymes with a known gene expression regulatory role in Escherichia coli, such as AlaS, PepA and BolA, and the uncharacterized hypothetical regulators YchA and YrbA). The relative positioning of regulatory genes along the chromosome of Buchnera seems to have conserved its ancestral state, despite the genome erosion. Sigma-70 promoters with canonical thermodynamic sequence profiles were detected upstream of about 94% of the CDS of Buchnera in the different aphids. Based on Stress-Induced Duplex Destabilization (SIDD) measurements, unstable σ70 promoters were found specifically associated with the regulator and transporter genes. CONCLUSIONS: This genomic analysis provides supporting evidence of a selection of functional regulatory structures and it has enabled us to propose hypotheses concerning possible links between these regulatory elements and the DNA-topology (i.e., supercoiling, curvature, flexibility and base-pair stability) in the regulation of gene expression in the shrunken genome of Buchnera.

Structure and dynamics of the operon map of Buchnera aphidicola sp. strain APS.

BACKGROUND: Gene expression regulation is still poorly documented in bacteria with highly reduced genomes. Understanding the evolution and mechanisms underlying the regulation of gene transcription in Buchnera aphidicola, the primary endosymbiont of aphids, is expected both to enhance our understanding of this nutritionally based association and to provide an intriguing case-study of the evolution of gene expression regulation in a reduced bacterial genome. RESULTS: A Bayesian predictor was defined to infer the B. aphidicola transcription units, which were further validated using transcriptomic data and RT-PCR experiments. The characteristics of B. aphidicola predicted transcription units (TUs) were analyzed in order to evaluate the impact of operon map organization on the regulation of gene transcription.On average, B. aphidicola TUs contain more genes than those of E. coli. The global layout of B. aphidicola operon map was mainly shaped by the big reduction and the rearrangements events, which occurred at the early stage of the symbiosis. Our analysis suggests that this operon map may evolve further only by small reorganizations around the frontiers of B. aphidicola TUs, through promoter and/or terminator sequence modifications and/or by pseudogenization events. We also found that the need for specific transcription regulation exerts some pressure on gene conservation, but not on gene assembling in the operon map in Buchnera. Our analysis of the TUs spacing pointed out that a selection pressure is maintained on the length of the intergenic regions between divergent adjacent gene pairs. CONCLUSIONS: B. aphidicola can seemingly only evolve towards a more polycistronic operon map. This implies that gene transcription regulation is probably subject to weak selection pressure in Buchnera conserving operons composed of genes with unrelated functions.

Immune signatures of protective spleen memory CD8 T cells.

Memory CD8 T lymphocyte populations are remarkably heterogeneous and differ in their ability to protect the host. In order to identify the whole range of qualities uniquely associated with protective memory cells we compared the gene expression signatures of two qualities of memory CD8 T cells sharing the same antigenic-specificity: protective (Influenza-induced, Flu-TM) and non-protective (peptide-induced, TIM) spleen memory CD8 T cells. Although Flu-TM and TIM express classical phenotypic memory markers and are polyfunctional, only Flu-TM protects against a lethal viral challenge. Protective memory CD8 T cells express a unique set of genes involved in migration and survival that correlate with their unique capacity to rapidly migrate within the infected lung parenchyma in response to influenza infection. We also enlighten a new set of poised genes expressed by protective cells that is strongly enriched in cytokines and chemokines such as Ccl1, Ccl9 and Gm-csf. CCL1 and GM-CSF genes are also poised in human memory CD8 T cells. These immune signatures are also induced by two other pathogens (vaccinia virus and Listeria monocytogenes). The immune signatures associated with immune protection were identified on circulating cells, i.e. those that are easily accessible for immuno-monitoring and could help predict vaccines efficacy.

Overexpression of the transcription factor Sp1 activates the OAS-RNAse L-RIG-I pathway.

Deregulated expression of oncogenes or transcription factors such as specificity protein 1 (Sp1) is observed in many human cancers and plays a role in tumor maintenance. Paradoxically in untransformed cells, Sp1 overexpression induces late apoptosis but the early intrinsic response is poorly characterized. In the present work, we studied increased Sp1 level consequences in untransformed cells and showed that it turns on an early innate immune transcriptome. Sp1 overexpression does not activate known cellular stress pathways such as DNA damage response or endoplasmic reticulum stress, but induces the activation of the OAS-RNase L pathway and the generation of small self-RNAs, leading to the upregulation of genes of the antiviral RIG-I pathway at the transcriptional and translational levels. Finally, Sp1-induced intrinsic innate immune response leads to the production of the chemokine CXCL4 and to the recruitment of inflammatory cells in vitro and in vivo. Altogether our results showed that increased Sp1 level in untransformed cells constitutes a novel danger signal sensed by the OAS-RNase L axis leading to the activation of the RIG-I pathway. These results suggested that the OAS-RNase L-RIG-I pathway may be activated in sterile condition in absence of pathogen.

Overexpression of transcription factor Sp1 leads to gene expression perturbations and cell cycle inhibition.

BACKGROUND: The ubiquitous transcription factor Sp1 regulates the expression of a vast number of genes involved in many cellular functions ranging from differentiation to proliferation and apoptosis. Sp1 expression levels show a dramatic increase during transformation and this could play a critical role for tumour development or maintenance. Although Sp1 deregulation might be beneficial for tumour cells, its overexpression induces apoptosis of untransformed cells. Here we further characterised the functional and transcriptional responses of untransformed cells following Sp1 overexpression. METHODOLOGY AND PRINCIPAL FINDINGS: We made use of wild-type and DNA-binding-deficient Sp1 to demonstrate that the induction of apoptosis by Sp1 is dependent on its capacity to bind DNA. Genome-wide expression profiling identified genes involved in cancer, cell death and cell cycle as being enriched among differentially expressed genes following Sp1 overexpression. In silico search to determine the presence of Sp1 binding sites in the promoter region of modulated genes was conducted. Genes that contained Sp1 binding sites in their promoters were enriched among down-regulated genes. The endogenous sp1 gene is one of the most down-regulated suggesting a negative feedback loop induced by overexpressed Sp1. In contrast, genes containing Sp1 binding sites in their promoters were not enriched among up-regulated genes. These results suggest that the transcriptional response involves both direct Sp1-driven transcription and indirect mechanisms. Finally, we show that Sp1 overexpression led to a modified expression of G1/S transition regulatory genes such as the down-regulation of cyclin D2 and the up-regulation of cyclin G2 and cdkn2c/p18 expression. The biological significance of these modifications was confirmed by showing that the cells accumulated in the G1 phase of the cell cycle before the onset of apoptosis. CONCLUSION: This study shows that the binding to DNA of overexpressed Sp1 induces an inhibition of cell cycle progression that precedes apoptosis and a transcriptional response targeting genes containing Sp1 binding sites in their promoter or not suggesting both direct Sp1-driven transcription and indirect mechanisms.

Systemic analysis of the symbiotic function of Buchnera aphidicola, the primary endosymbiont of the pea aphid Acyrthosiphon pisum.

Buchnera aphidicola is the primary obligate intracellular symbiont of most aphid species. B. aphidicola and aphids have been evolving in parallel since their association started, about 150 Myr ago. Both partners have lost their autonomy, and aphid diversification has been confined to smaller ecological niches by this co-evolution. B. aphidicola has undergone major genomic and biochemical changes as a result of adapting to intracellular life. Several genomes of B. aphidicola from different aphid species have been sequenced in the last decade, making it possible to carry out analyses and comparative studies using system-level in silico methods. This review attempts to provide a systemic description of the symbiotic function of aphid endosymbionts, particularly of B. aphidicola from the pea aphid Acyrthosiphon pisum, by analyzing their structural genomic properties, as well as their genetic and metabolic networks.