Visualization and quality assessment of de novo genome assemblies.

Recent technological progress has greatly facilitated de novo genome sequencing. However, de novo assemblies consist in many pieces of contiguous sequence (contigs) arranged in thousands of scaffolds instead of small numbers of chromosomes. Confirming and improving the quality of such assemblies is critical for subsequent analysis. We present a method to evaluate genome scaffolding by aligning independently obtained transcriptome sequences to the genome and visually summarizing the alignments using the Cytoscape software. Applying this method to the genome of the red fire ant Solenopsis invicta allowed us to identify inconsistencies in 7%, confirm contig order in 20% and extend 16% of scaffolds.Scripts that generate tables for visualization in Cytoscape from FASTA sequence and scaffolding information files are publicly available at https://github.com/ksanao/TGNet.

Test of four colon cancer risk-scores in formalin fixed paraffin embedded microarray gene expression data.

BACKGROUND: Prognosis prediction for resected primary colon cancer is based on the T-stage Node Metastasis (TNM) staging system. We investigated if four well-documented gene expression risk scores can improve patient stratification. METHODS: Microarray-based versions of risk-scores were applied to a large independent cohort of 688 stage II/III tumors from the PETACC-3 trial. Prognostic value for relapse-free survival (RFS), survival after relapse (SAR), and overall survival (OS) was assessed by regression analysis. To assess improvement over a reference, prognostic model was assessed with the area under curve (AUC) of receiver operating characteristic (ROC) curves. All statistical tests were two-sided, except the AUC increase. RESULTS: All four risk scores (RSs) showed a statistically significant association (single-test, P < .0167) with OS or RFS in univariate models, but with HRs below 1.38 per interquartile range. Three scores were predictors of shorter RFS, one of shorter SAR. Each RS could only marginally improve an RFS or OS model with the known factors T-stage, N-stage, and microsatellite instability (MSI) status (AUC gains < 0.025 units). The pairwise interscore discordance was never high (maximal Spearman correlation = 0.563) A combined score showed a trend to higher prognostic value and higher AUC increase for OS (HR = 1.74, 95% confidence interval [CI] = 1.44 to 2.10, P < .001, AUC from 0.6918 to 0.7321) and RFS (HR = 1.56, 95% CI = 1.33 to 1.84, P < .001, AUC from 0.6723 to 0.6945) than any single score. CONCLUSIONS: The four tested gene expression-based risk scores provide prognostic information but contribute only marginally to improving models based on established risk factors. A combination of the risk scores might provide more robust information. Predictors of RFS and SAR might need to be different.

LRF-mediated Dll4 repression in erythroblasts is necessary for hematopoietic stem cell maintenance.

Hematopoietic stem cells (HSCs) are the most primitive cells in the hematopoietic system and are under tight regulation for self-renewal and differentiation. Notch signals are essential for the emergence of definitive hematopoiesis in mouse embryos and are critical regulators of lymphoid lineage fate determination. However, it remains unclear how Notch regulates the balance between HSC self-renewal and differentiation in the adult bone marrow (BM). Here we report a novel mechanism that prevents HSCs from undergoing premature lymphoid differentiation in BM. Using a series of in vivo mouse models and functional HSC assays, we show that leukemia/lymphoma related factor (LRF) is necessary for HSC maintenance by functioning as an erythroid-specific repressor of Delta-like 4 (Dll4) expression. Lrf deletion in erythroblasts promoted up-regulation of Dll4 in erythroblasts, sensitizing HSCs to T-cell instructive signals in the BM. Our study reveals novel cross-talk between HSCs and erythroblasts, and sheds a new light on the regulatory mechanisms regulating the balance between HSC self-renewal and differentiation.

Mice carrying ubiquitin-specific protease 2 (Usp2) gene inactivation maintain normal sodium balance and blood pressure.

Ubiquitylation plays an important role in the control of Na⁺ homeostasis by the kidney. It is well established that the epithelial Na⁺ channel ENaC is regulated by the ubiquitin-protein ligase NEDD4-2, limiting ENaC cell surface expression and activity. Ubiquitylation can be reversed by the action of deubiquitylating enzymes (DUBs). One such DUB, USP2-45, was identified previously as an aldosterone-induced protein in the kidney and is also a circadian output gene. In heterologous expression systems, USP2-45 binds to ENaC, deubiquitylates it, and enhances channel density and activity at the cell surface. Because the role of USP2-45 in renal Na⁺ transport had not been studied in vivo, we investigated here the effect of Usp2 gene inactivation in this process. We demonstrate first that USP2-45 protein has a rhythmic expression with a peak at ZT12. Usp2-KO mice did not show any differences from wild-type littermates with respect to the diurnal control of Na⁺ or K⁺ urinary excretion and plasma levels either on a standard diet or after acute and chronic changes to low- and high-Na⁺ diets, respectively. Moreover, they had similar aldosterone levels on either a low- or high-Na⁺ diet. Blood pressure measurements using telemetry did not reveal variations compared with control mice. Usp2-KO mice did not display alterations in expression of genes involved in sodium homeostasis or the ubiquitin system, as evidenced by transcriptome analysis in the kidney. Our data suggest that USP2 does not play a primary role in the control of Na⁺ balance or blood pressure.

The evolution of gene expression levels in mammalian organs.

Changes in gene expression are thought to underlie many of the phenotypic differences between species. However, large-scale analyses of gene expression evolution were until recently prevented by technological limitations. Here we report the sequencing of polyadenylated RNA from six organs across ten species that represent all major mammalian lineages (placentals, marsupials and monotremes) and birds (the evolutionary outgroup), with the goal of understanding the dynamics of mammalian transcriptome evolution. We show that the rate of gene expression evolution varies among organs, lineages and chromosomes, owing to differences in selective pressures: transcriptome change was slow in nervous tissues and rapid in testes, slower in rodents than in apes and monotremes, and rapid for the X chromosome right after its formation. Although gene expression evolution in mammals was strongly shaped by purifying selection, we identify numerous potentially selectively driven expression switches, which occurred at different rates across lineages and tissues and which probably contributed to the specific organ biology of various mammals.

Phenotypic switching in Pseudomonas brassicacearum involves GacS- and GacA-dependent Rsm small RNAs.

The plant-beneficial bacterium Pseudomonas brassicacearum forms phenotypic variants in vitro as well as in planta during root colonization under natural conditions. Transcriptome analysis of typical phenotypic variants using microarrays containing coding as well as noncoding DNA fragments showed differential expression of several genes relevant to secondary metabolism and of the small RNA (sRNA) genes rsmX, rsmY, and rsmZ. Naturally occurring mutations in the gacS-gacA system accounted for phenotypic switching, which was characterized by downregulation of antifungal secondary metabolites (2,4-diacetylphloroglucinol and cyanide), indoleacetate, exoenzymes (lipase and protease), and three different N-acyl-homoserine lactone molecules. Moreover, in addition to abrogating these biocontrol traits, gacS and gacA mutations resulted in reduced expression of the type VI secretion machinery, alginate biosynthesis, and biofilm formation. In a gacA mutant, the expression of rsmX was completely abolished, unlike that of rsmY and rsmZ. Overexpression of any of the three sRNAs in the gacA mutant overruled the pleiotropic changes and restored the wild-type phenotypes, suggesting functional redundancy of these sRNAs. In conclusion, our data show that phenotypic switching in P. brassicacearum results from mutations in the gacS-gacA system.

KAP1 regulates gene networks controlling T-cell development and responsiveness.

Chromatin remodeling at specific genomic loci controls lymphoid differentiation. Here, we investigated the role played in this process by Kruppel-associated box (KRAB)-associated protein 1 (KAP1), the universal cofactor of KRAB-zinc finger proteins (ZFPs), a tetrapod-restricted family of transcriptional repressors. T-cell-specific Kap1-deleted mice displayed a significant expansion of immature thymocytes, imbalances in CD4(+)/CD8(+) cell ratios, and altered responses to TCR and TGFβ stimulation when compared to littermate KAP1 control mice. Transcriptome and chromatin studies revealed that KAP1 binds T-cell-specific cis-acting regulatory elements marked by the H3K9me3 repressive mark and enriched in Ikaros/NuRD complexes. Also, KAP1 directly controls the expression of several genes involved in TCR and cytokine signaling. Among these, regulation of FoxO1 seems to play a major role in this system. Likely responsible for tethering KAP1 to at least part of its genomic targets, a small number of KRAB-ZFPs are selectively expressed in T-lymphoid cells. These results reveal the so far unsuspected yet important role of KAP1-mediated epigenetic regulation in T-lymphocyte differentiation and activation.

MALT1 auto-proteolysis is essential for NF-κB-dependent gene transcription in activated lymphocytes.

Mucosa-associated lymphoid tissue 1 (MALT1) controls antigen receptor-mediated signalling to nuclear factor κB (NF-κB) through both its adaptor and protease function. Upon antigen stimulation, MALT1 forms a complex with BCL10 and CARMA1, which is essential for initial IκBα phosphorylation and NF-κB nuclear translocation. Parallel induction of MALT1 protease activity serves to inactivate negative regulators of NF-κB signalling, such as A20 and RELB. Here we demonstrate a key role for auto-proteolytic MALT1 cleavage in B- and T-cell receptor signalling. MALT1 cleavage occurred after Arginine 149, between the N-terminal death domain and the first immunoglobulin-like region, and did not affect its proteolytic activity. Jurkat T cells expressing an un-cleavable MALT1-R149A mutant showed unaltered initial IκBα phosphorylation and normal nuclear accumulation of NF-κB subunits. Nevertheless, MALT1 cleavage was required for optimal activation of NF-κB reporter genes and expression of the NF-κB targets IL-2 and CSF2. Transcriptome analysis confirmed that MALT1 cleavage after R149 was required to induce NF-κB transcriptional activity in Jurkat T cells. Collectively, these data demonstrate that auto-proteolytic MALT1 cleavage controls antigen receptor-induced expression of NF-κB target genes downstream of nuclear NF-κB accumulation.

Evolutionary dynamics of coding and non-coding transcriptomes.

Gene expression changes may underlie much of phenotypic evolution. The development of high-throughput RNA sequencing protocols has opened the door to unprecedented large-scale and cross-species transcriptome comparisons by allowing accurate and sensitive assessments of transcript sequences and expression levels. Here, we review the initial wave of the new generation of comparative transcriptomic studies in mammals and vertebrate outgroup species in the context of earlier work. Together with various large-scale genomic and epigenomic data, these studies have unveiled commonalities and differences in the dynamics of gene expression evolution for various types of coding and non-coding genes across mammalian lineages, organs, developmental stages, chromosomes and sexes. They have also provided intriguing new clues to the regulatory basis and phenotypic implications of evolutionary gene expression changes.

Bacterial subfamily of LuxR regulators that respond to plant compounds.

Pseudomonas fluorescens are rhizobacteria known for their biocontrol properties. Several antimicrobial functions are crucial for this process, and the experiments described here investigate the modulation of their expression during the plant-bacterium interaction. The role of a LuxR family regulator in interkingdom signaling has been investigated using genome-scale transcriptome analysis, gene promoter studies in vivo and in vitro, biocontrol assays, and response to plant compounds. PsoR, a LuxR solo or orphan regulator of P. fluorescens, was identified. PsoR is solubilized and activates a lux-box-containing promoter only in the presence of macerated plants, suggesting the presence of a plant molecule(s) that most likely binds to PsoR. Gene expression profiles revealed that genes involved in the inhibition of plant pathogens were affected by PsoR, including a chitinase gene, iron metabolism genes, and biosynthetic genes of antifungal compounds. 2,4-Diacetylphloroglucinol production is PsoR dependent both in vitro and in vivo. psoR mutants were significantly reduced for their ability to protect wheat plants from root rot, and damping-off caused by Pythium ultimum infection. PsoR most likely senses a molecule(s) in the plant and modulates expression of genes that have a role in biocontrol. PsoR and related proteins form a subfamily of LuxR family regulators in plant-associated bacteria.

State of genomics and epigenomics research in the perspective of HIV cure.

PURPOSE OF REVIEW: One of the seven key scientific priorities identified in the road map on HIV cure research is to 'determine the host mechanisms that control HIV replication in the absence of therapy'. This review summarizes the recent work in genomics and in epigenetic control of viral replication that is relevant for this mission. RECENT FINDINGS: New technologies allow the joint analysis of host and viral transcripts. They identify the patterns of antisense transcription of the viral genome and its role in gene regulation. High-throughput studies facilitate the assessment of integration at the genome scale. Integration site, orientation and host genomic context modulate the transcription and should also be assessed at the level of single cells. The various models of latency in primary cells can be followed using dynamic study designs to acquire transcriptome and proteome data of the process of entry, maintenance and reactivation of latency. Dynamic studies can be applied to the study of transcription factors and chromatin modifications in latency and upon reactivation. SUMMARY: The convergence of primary cell models of latency, new high-throughput quantitative technologies applied to the study of time series and the identification of compounds that reactivate viral transcription bring unprecedented precision to the study of viral latency.

Molecular and Functional Characterization of Neuroblastoma-Initiating Cells

Neuroblastoma (NB) is the most common extracranial malignant tumor in young children and arises at any site of the sympathetic nervous system. The disease exhibits a remarkable phenotypic diversity ranging from spontaneous regression to fatal disease. Poor outcome results from a rapidly progressive, metastatic and drug-resistant disease. Recent studies have suggested that solid tumors may arise from a minor population of cancer stem cells (CSCs) with stem cell markers and typical properties such as self-renewal ability, asymmetric division and drug resistance. In this model, CSCs possess the exclusive ability to initiate and maintain the tumor, and to produce distant metastases. Tumor cell subpopulations with stem-like phenotypes have indeed been identified in several cancer including leukemia, breast, brain and colon cancers. CSC hypothesis still needs to be validated in the other cancers including NB.NB originates from neural crest-derived malignant sympatho-adrenal cells. We have identified rare cells that express markers in conformity with neural crest stem cells and their derived lineages within primary NB tissue and cell lines, leading us to postulate the existence of CSCs in NB tumors.In the absence of specific markers to isolate CSCs, we adapted to NB tumor cells the sphere functional assay, based on the ability of stem cells to grow as spheres in non-adherent conditions. By serial passages of spheres from bone marrow NB metastases, a subset of cells was gradually selected and its specific gene expression profile identified by micro-array time-course analysis. The differentially expressed genes in spheres are enriched in genes implicated in development including CD133, ABC-transporters, WNT and NOTCH genes, identified in others solid cancers as CSCs markers, and other new markers, all referred by us as the Neurosphere Expression Profile (NEP). We confirmed the presence of a cell subpopulation expressing a combination of the NEP markers within a few primary NB samples.The tumorigenic potential of NB spheres was assayed by in vivo tumor growth analyses using orthotopic (adrenal glands) implantations of tumor cells into immune-compromised mice. Tumors derived from the sphere cells were significantly more frequent and were detected earlier compared to whole tumor cells. However, NB cells expressing the neurosphere-associated genes and isolated from the bulk tumors did not recapitulate the CSC-like phenotype in the orthotopic model. In addition, the NB sphere cells lost their higher tumorigenic potential when implanted in a subcutaneous heterotopic in vivo model.These results highlighted the complex behavior of CSC functions and led us to consider the stem-like NB cells as a dynamic and heterogeneous cell population influenced by microenvironment signals.Our approach identified for the first time candidate genes that may be associated with NB self-renewal and tumorigenicity and therefore would establish specific functional targets for more effective therapies in aggressive NB.

Douleur chronique et mécanismes moléculaires

La douleur est définie par l'International association for the study of pain (IASP) comme une expérience sensorielle et émotionnelle désagréable, associée à une lésion tissulaire réelle ou à une lésion potentielle, ou décrite en des termes évoquant une telle lésion. Sa fonction est de signaler une menace potentielle pour l'intégrité de l'organisme. Mais ce ressenti peut persister ou être présent en l'absence d'une telle menace. Il s'agit dans ce cas d'une douleur pathologique dont les mécanismes étiologiques et physiopathologiques ne sont pas encore bien compris.¦Ce travail de maîtrise a pour objectif l'étude d'une mutation génétique responsable d'un syndrome douloureux chronique. Cette mutation génétique a été décelée chez une famille de patients lausannois atteints de PEPD (paroxysmal extreme pain disorder) et touche le canal sodique Nav1.7. Ce canal est exprimé principalement dans les neurones des ganglions sensitifs et des ganglions sympathiques. Il est considéré comme responsable de la transmission de la douleur vers le SNC car des mutations « perte de fonction » de ce canal sont à l'origine d'une insensibilité congénitale à la douleur. Les mutations « gain de fonction » de ce canal sont à l'origine de syndromes douloureux chroniques tel le syndrome PEPD ou l'érythromélalgie. La mutation présente chez les patients lausannois est située entre les segments transmembranaires S4 et S5 sur le 4ème domaine de la sous-unité α du canal sodique Nav1.7. Cette mutation ne touche qu'un seul acide aminé, en position 1612 où une leucine est remplacée par une proline.¦Les méthodes employées dans ce travail sont la mutagenèse pour générer des plasmides contenant le gène SCN9A muté (T4835C) codant pour le canal Nav1.7 muté (L1612P), l'amplification de ces plasmides dans des bactéries et la transfection de cellules HEK293 avec les plasmides contenant le gène SCN9A muté (T4835C). Nous avons ainsi pu induire l'expression du canal muté dans des cellules HEK293. Ces cellules pourront être utilisées par la suite pour enregistrer les courants ioniques transitant à travers les canaux exprimés à la membrane plasmique. Cela permettra de comparer les propriétés électrophysiologiques du canal Nav1.7 muté L1612P avec celles du canal non muté. Nous avons également recherché l'expression d'ARNm codant pour les composants des canaux sodiques (sous-unités α et β) dans les cellules HEK293 non transfectées par la technique de qRT-PCR. Ceci afin de répertorier les composants des canaux sodiques exprimés constitutivement par les cellules HEK293 qui pourraient avoir une influence sur les mesures électrophysiologiques qui seront effectuées sur ces cellules.¦Ce travail a permis de générer des plasmides contenant le gène SCN9A muté T4835C qui sont des outils nécessaires à la réalisation d'études plus détaillées sur le fonctionnement et les propriétés du canal Nav1.7 muté L1612P. Ce travail a également permis, par la méthode de la qRT-PCR, une analyse de l'expression d'ARNm codant pour la sous-unité α du canal Nav1.7 et des sous-unités β 1 à 4 par les cellules HEK293. Les résultats ainsi obtenus permettent de mieux caractériser le transcriptome des cellules HEK293 et seront utiles pour interpréter avec plus de précisions les expérimentations électrophysiologiques utilisant ces cellules. L'étude électrophysiologique des cellules HEK293 exprimant le canal Nav1.7 muté par la technique du patch clamp est en cours. Elle sera poursuivie dans le cadre d'un travail de recherche dépassant le cadre de ce travail de maîtrise. Les cellules HEK293 exprimant le canal Nav1.7 muté (L1612P) pourront être également utilisées pour tester l'effet de divers médicaments sur ce canal. Ce qui pourrait d'une part permettre d'optimiser le traitement des patients souffrant de PEPD et d'autre part être utile pour tout traitement à but antalgique.

Genome-wide transposon insertion scanning of environmental survival functions in the polycyclic aromatic hydrocarbon degrading bacterium Sphingomonas wittichii RW1.

Sphingomonas wittichii RW1 is a dibenzofuran and dibenzodioxin-degrading bacterium with potentially interesting properties for bioaugmentation of contaminated sites. In order to understand the capacity of the microorganism to survive in the environment we used a genome-wide transposon scanning approach. RW1 transposon libraries were generated with around 22 000 independent insertions. Libraries were grown for an average of 50 generations (five successive passages in batch liquid medium) with salicylate as sole carbon and energy source in presence or absence of salt stress at -1.5 MPa. Alternatively, libraries were grown in sand with salicylate, at 50% water holding capacity, for 4 and 10 days (equivalent to 7 generations). Library DNA was recovered from the different growth conditions and scanned by ultrahigh throughput sequencing for the positions and numbers of inserted transposed kanamycin resistance gene. No transposon reads were recovered in 579 genes (10% of all annotated genes in the RW1 genome) in any of the libraries, suggesting those to be essential for survival under the used conditions. Libraries recovered from sand differed strongly from those incubated in liquid batch medium. In particular, important functions for survival of cells in sand at the short term concerned nutrient scavenging, energy metabolism and motility. In contrast to this, fatty acid metabolism and oxidative stress response were essential for longer term survival of cells in sand. Comparison to transcriptome data suggested important functions in sand for flagellar movement, pili synthesis, trehalose and polysaccharide synthesis and putative cell surface antigen proteins. Interestingly, a variety of genes were also identified, interruption of which cause significant increase in fitness during growth on salicylate. One of these was an Lrp family transcription regulator and mutants in this gene covered more than 90% of the total library after 50 generations of growth on salicylate. Our results demonstrate the power of genome-wide transposon scanning approaches for analysis of complex traits.

First-generation linkage map for the European tree frog (Hyla arborea) with utility in congeneric species.

BACKGROUND: Western Palearctic tree frogs (Hyla arborea group) represent a strong potential for evolutionary and conservation genetic research, so far underexploited due to limited molecular resources. New microsatellite markers have recently been developed for Hyla arborea, with high cross-species utility across the entire circum-Mediterranean radiation. Here we conduct sibship analyses to map available markers for use in future population genetic applications. FINDINGS: We characterized eight linkage groups, including one sex-linked, all showing drastically reduced recombination in males compared to females, as previously documented in this species. Mapping of the new 15 markers to the ~200 My diverged Xenopus tropicalis genome suggests a generally conserved synteny with only one confirmed major chromosome rearrangement. CONCLUSIONS: The new microsatellites are representative of several chromosomes of H. arborea that are likely to be conserved across closely-related species. Our linkage map provides an important resource for genetic research in European Hylids, notably for studies of speciation, genome evolution and conservation.