Outbreak of Invasive Wound Mucormycosis in a Burn Unit Due to Multiple Strains of Mucor circinelloides f. circinelloides Resolved by Whole-Genome Sequencing – Garcia-Hermoso D, Criscuolo A, Lee SC, Legrand M, Chaouat M, Denis B, Lafaurie M, Rouveau M, Soler C, Schaal J-V, Mimoun M, Mebazaa A, Heitman J, Dromer F, Brisse S, Bretagne S, Alanio A

Mucorales are ubiquitous environmental molds responsible for mucormycosis in diabetic, immunocompromised, and severely burned patients. Small outbreaks of invasive wound mucormycosis (IWM) have already been reported in burn units without extensive microbiological investigations. We faced an outbreak of IWM in our center and investigated the clinical isolates with whole-genome sequencing (WGS) analysis. We analyzed M. circinelloides isolates from patients in our burn unit (BU1, Hôpital Saint-Louis, Paris, France) together with nonoutbreak isolates from Burn Unit 2 (BU2, Paris area) and from France over a 2-year period (2013 to 2015). A total of 21 isolates, including 14 isolates from six BU1 patients, were analyzed by whole-genome sequencing (WGS). Phylogenetic classification based on de novo assembly and assembly free approaches showed that the clinical isolates clustered in four highly divergent clades. Clade 1 contained at least one of the strains from the six epidemiologically linked BU1 patients. The clinical isolates were specific to each patient. Two patients were infected with more than two strains from different clades, suggesting that an environmental reservoir of clonally unrelated isolates was the source of contamination. Only two patients from BU1 shared one strain, which could correspond to direct transmission or contamination with the same environmental source. In conclusion, WGS of several isolates per patients coupled with precise epidemiological data revealed a complex situation combining potential cross-transmission between patients and multiple contaminations with a heterogeneous pool of strains from a cryptic environmental reservoir.

Parasites & vectors
Diverse laboratory colonies of Aedes aegypti harbor the same adult midgut bacterial microbiome – Dickson LB, Ghozlane A, Volant S, Bouchier C, Ma L, Vega-Rúa A, Dusfour I, Jiolle D, Paupy C, Mayanja MN, Kohl A, Lutwama JJ, Duong V, Lambrechts L

BACKGROUND: Host-associated microbes, collectively known as the microbiota, play an important role in the biology of multicellular organisms. In mosquito vectors of human pathogens, the gut bacterial microbiota influences vectorial capacity and has become the subject of intense study. In laboratory studies of vector biology, genetic effects are often inferred from differences between geographically and genetically diverse colonies of mosquitoes that are reared in the same insectary. It is unclear, however, to what extent genetic effects can be confounded by uncontrolled differences in the microbiota composition among mosquito colonies. To address this question, we used 16S metagenomics to compare the midgut bacterial microbiome of six laboratory colonies of Aedes aegypti recently derived from wild populations representing the geographical range and genetic diversity of the species. RESULTS: We found that the diversity, abundance, and community structure of the midgut bacterial microbiome was remarkably similar among the six different colonies of Ae. aegypti, regardless of their geographical origin. We also confirmed the relatively low complexity of bacterial communities inhabiting the mosquito midgut. CONCLUSIONS: Our finding that geographically diverse colonies of Ae. aegypti reared in the same insectary harbor a similar gut bacterial microbiome supports the conclusion that the gut microbiota of adult mosquitoes is environmentally determined regardless of the host genotype. Thus, uncontrolled differences in microbiota composition are unlikely to represent a significant confounding factor in genetic studies of vector biology.

Genome Announcement
Draft Genome Sequence of the Fish Pathogen Flavobacterium columnare Genomovar III Strain PH-97028 (=CIP 109753) – Criscuolo A, Chesneau O, Clermont D, Bizet C

Flavobacterium columnare strain PH-97028 (=CIP 109753) is a genomovar III reference strain that was isolated from a diseased Ayu fish in Japan. We report here the analysis of the first available genomovar III sequence of this species to aid in identification, epidemiological tracking, and virulence studies.

Sequanix: A Dynamic Graphical Interface for Snakemake Workflows. – Desvillechabrol Dimitri, Legendre Rachel, Rioualen Claire, Bouchier Christiane, van Helden Jacques, Kennedy Sean, Cokelaer Thomas

Summary We designed a PyQt graphical user interface – Sequanix – aimed at democratizing the use of Snakemake pipelines in the NGS space and beyond. By default, Sequanix includes Sequana NGS pipelines (Snakemake format) (, and is also capable of loading any external Snakemake pipeline. New users can easily, visually, edit configuration files of expert-validated pipelines and can interactively execute these production-ready workflows. Sequanix will be useful to both Snakemake developers in exposing their pipelines and to a wide audience of users. Availability: Source on, bio-containers on and Singularity hub (

PLOS Neglected Tropical Diseases
Individual co-variation between viral RNA load and gene expression reveals novel host factors during early dengue virus infection of the Aedes aegypti midgut – Vincent Raquin, Sarah Hélène Merkling, Valérie Gausson, Isabelle Moltini-Conclois, Lionel Frangeul, Hugo Varet, Marie-Agnès Dillies, Maria-Carla Saleh, Louis Lambrechts

Dengue virus (DENV) causes more human infections than any other mosquito-borne virus. The current lack of antiviral strategies has prompted genome-wide screens for host genes that are required for DENV infectivity. Earlier transcriptomic studies that identified DENV host factors in the primary vector Aedes aegypti used inbred laboratory colonies and/or pools of mosquitoes that erase individual variation. Here, we performed transcriptome sequencing on individual midguts in a field-derived Ae. aegypti population to identify new candidate host factors modulating DENV replication. We analyzed the transcriptomic data using an approach that accounts for individual co-variation between viral RNA load and gene expression. This approach generates a prediction about the agonist or antagonist effect of candidate genes on DENV replication based on the sign of the correlation between gene expression and viral RNA load. Using this method, we identified 39 candidate genes that went undetected by conventional pairwise comparison of gene expression levels between DENV-infected midguts and uninfected controls. Only four candidate genes were detected by both methods, emphasizing their complementarity. We demonstrated the value of our approach by functional validation of a candidate agonist gene encoding a sterol regulatory element-binding protein (SREBP), which was identified by correlation analysis but not by pairwise comparison. We confirmed that SREBP promotes DENV infection in the midgut by RNAi-mediated gene knockdown in vivo. We suggest that our approach for transcriptomic analysis can empower genome-wide screens for potential agonist or antagonist factors by leveraging inter-individual variation in gene expression. More generally, this method is applicable to a wide range of phenotypic traits displaying inter-individual variation.

Journal of Crohn’s and Colitis
Characteristics of Fecal Microbiota in Pediatric Crohn’s Disease and Their Dynamic Changes During Infliximab Therapy. – Yizhong Wang, Xuefeng Gao, Amine Ghozlane, Hui Hu, Xiaolu Li, Yongmei Xiao, Dan Li, Guangjun Yu, Ting Zhang

BACKGROUND: Crohn’s disease [CD] is known to be associated with gut microbial dysbiosis. Infliximab [IFX] is increasingly used to treat pediatric CD, however it is not clear how gut microbiota is modified during IFX treatment. The aim of this study was to characterize the fecal microbiota community composition in pediatric CD patients and to assess its dynamic changes during IFX therapy. METHODS: 16S rRNA sequencing approach was applied to determine the compositions of microbial communities in fecal samples. The composition and function of the fecal microbiota were compared between CD patients and health controls. RESULTS: Characteristics of fecal microbiome composition in pediatric CD patients prior to IFX treatment were represented by a lower biodiversity, a gain in Enterococcus, and a significant loss in multiple short-chain fatty acid [SCFA]-producing bacteria, including Anaerostipes, Blautia, Coprococcus, Faecalibacterium, Lachnospira, Odoribacter, Roseburia, Ruminococcus, and Sutterella. Additionally, alterations were observed in metabolic functions of the gut microbial community in CD. IFX treatment increased the biodiversity of gut microbiota and shifted its composition as well as its functional capabilities in the pediatric CD patients toward a healthy status. However, multiple SCFA-producing taxa were not significantly expanded. The sustained response of pediatric CD patients to IFX was associated with abundance of SCFA-producing bacteria. CONCLUSIONS: A lower biodiversity with alterations in the composition and function of fecal microbial community, characterizing gut microbial dysbiosis, were observed in the Chinese pediatric CD patients. IFX diminished the CD-associated gut microbial dysbiosis but was deficient in increasing certain SCFA-producing taxa.

Science Advances
Carryover effects of larval exposure to different environmental bacteria drive adult trait variation in a mosquito vector – Laura B. Dickson, Davy Jiolle, Guillaume Minard, Isabelle Moltini-Conclois, Stevenn Volant, Amine Ghozlane, Christiane Bouchier, Diego Ayala, Christophe Paupy, Claire Valiente Moro, Louis Lambrechts

Conditions experienced during larval development of holometabolous insects can affect adult traits, but whether differences in the bacterial communities of larval development sites contribute to variation in the ability of insect vectors to transmit human pathogens is unknown. We addressed this question in the mosquito Aedes aegypti, a major arbovirus vector breeding in both sylvatic and domestic habitats in Sub-Saharan Africa. Targeted metagenomics revealed differing bacterial communities in the water of natural breeding sites in Gabon. Experimental exposure to different native bacterial isolates during larval development resulted in significant differences in pupation rate and adult body size but not life span. Larval exposure to an Enterobacteriaceae isolate resulted in decreased antibacterial activity in adult hemolymph and reduced dengue virus dissemination titer. Together, these data provide the proof of concept that larval exposure to different bacteria can drive variation in adult traits underlying vectorial capacity. Our study establishes a functional link between larval ecology, environmental microbes, and adult phenotypic variation in a holometabolous insect vector.

PLOS Pathogens
LipL21 lipoprotein binding to peptidoglycan enables Leptospira interrogans to escape NOD1 and NOD2 recognition – Ratet G, Santecchia I, Fanton d’Andon M, Vernel-Pauillac F, Wheeler R, Lenormand P, Fischer F, Lechat P, Haake DA, Picardeau M, Boneca IG, Werts C

Leptospirosis is a widespread zoonosis, potentially severe in humans, caused by spirochetal bacteria, Leptospira interrogans (L. interrogans). Host defense mechanisms involved in leptospirosis are poorly understood. Recognition of lipopolysaccharide (LPS) and lipoproteins by Toll-Like Receptors (TLR)4 and TLR2 is crucial for clearance of leptospires in mice, yet the role of Nucleotide Oligomerization Domain (NOD)-like receptors (NOD)1 and NOD2, recognizing peptidoglycan (PG) fragments has not previously been examined. Here, we show that pathogenic leptospires escape from NOD1 and NOD2 recognition both in vitro and in vivo, in mice. We found that leptospiral PG is resistant to digestion by certain hydrolases and that a conserved outer membrane lipoprotein of unknown function, LipL21, specific for pathogenic leptospires, is tightly bound to the PG. Leptospiral PG prepared from a mutant not expressing LipL21 (lipl21-) was more readily digested than the parental or complemented strains. Muropeptides released from the PG of the lipl21- mutant, or prepared using a procedure to eliminate the LipL21 protein from the PG of the parental strain, were recognized in vitro by the human NOD1 (hNOD1) and NOD2 (hNOD2) receptors, suggesting that LipL21 protects PG from degradation into muropeptides. LipL21 expressed in E. coli also resulted in impaired PG digestion and NOD signaling. We found that murine NOD1 (mNOD1) did not recognize PG of L. interrogans. This result was confirmed by mass spectrometry showing that leptospiral PG was primarily composed of MurTriDAP, the natural agonist of hNOD1, and contained only trace amounts of the tetra muropeptide, the mNOD1 agonist. Finally, in transgenic mice expressing human NOD1 and deficient for the murine NOD1, we showed enhanced clearance of a lipl21- mutant compared to the complemented strain, or to what was observed in NOD1KO mice, suggesting that LipL21 facilitates escape from immune surveillance in humans. These novel mechanisms allowing L. interrogans to escape recognition by the NOD receptors may be important in circumventing innate host responses.

BMC Genomics
Unraveling the evolution and coevolution of small regulatory RNAs and coding genes in Listeria – Franck Cerutti, Ludovic Mallet, Anaïs Painset, Claire Hoede, Annick Moisan, Christophe Bécavin, Mélodie Duval, Olivier Dussurget, Pascale Cossart, Christine Gaspin and Hélène Chiapello

Background: Small regulatory RNAs (sRNAs) are widely found in bacteria and play key roles in many important physiological and adaptation processes. Studying their evolution and screening for events of coevolution with other genomic features is a powerful way to better understand their origin and assess a common functional or adaptive relationship between them. However, evolution and coevolution of sRNAs with coding genes have been sparsely investigated in bacterial pathogens. Results: We designed a robust and generic phylogenomics approach that detects correlated evolution between sRNAs and protein-coding genes using their observed and inferred patterns of presence-absence in a set of annotated genomes. We applied this approach on 79 complete genomes of the Listeria genus and identified fifty-two accessory sRNAs, of which most were present in the Listeria common ancestor and lost during Listeria evolution. We detected significant coevolution between 23 sRNA and 52 coding genes and inferred the Listeria sRNA-coding genes coevolution network. We characterized a main hub of 12 sRNAs that coevolved with genes encoding cell wall proteins and virulence factors. Among them, an sRNA specific to L. monocytogenes species, rli133, coevolved with genes involved either in pathogenicity or in interaction with host cells, possibly acting as a direct negative post-transcriptional regulation. Conclusions: Our approach allowed the identification of candidate sRNAs potentially involved in pathogenicity and host interaction, consistent with recent findings on known pathogenicity actors. We highlight four sRNAs coevolving with seven internalin genes, some of which being important virulence factors in Listeria.

Cell Reports
In Situ Fixation Redefines Quiescence and Early Activation of Skeletal Muscle Stem Cells – Léo Machado , Joana Esteves de Lima , Odile Fabre , Caroline Proux , Rachel Legendre , Anikó Szegedi , Hugo Varet , Lars Roed Ingerslev , Romain Barrès , Frédéric Relaix, Philippos Mourikis

State of the art techniques have been developed to isolate and analyze cells from various tissues, aiming to capture their in vivo state. However, the majority of cell isolation protocols involve lengthy mechanical and enzymatic dissociation steps followed by flow cytometry, exposing cells to stress and disrupting their physiological niche. Focusing on adult skeletal muscle stem cells, we have developed a protocol that circumvents the impact of isolation procedures and captures cells in their native quiescent state. We show that current isolation protocols induce major transcriptional changes accompanied by specific histone modifications while having negligible effects on DNA methylation. In addition to proposing a protocol to avoid isolation-induced artifacts, our study reveals previously undetected quiescence and early activation genes of potential biological interest.

The Journal of Infectious Diseases
Sustained malaria control over an eight-year period in Papua New Guinea: the challenge of low-density asymptomatic infections – Koepfli C, Ome-Kaius M, Jally S, Malau E, Maripal S, Ginny J, Timinao L, Kattenberg JH, Obadia T, White M, Rarau P, Senn N, Barry AE, Kazura JW, Mueller I, Robinson LJ

BACKGROUND: The scale-up of effective malaria control in the last decade has resulted in a substantial decline in the incidence of clinical malaria in many countries. The effects on the proportions of asymptomatic and submicroscopic infections, and on transmission potential are yet poorly understood. METHODS: In Papua New Guinea, vector control has been intensified since 2008, and improved diagnosis and treatment introduced in 2012. Cross-sectional surveys were conducted in Madang Province in 2006 (n=1280), 2010 (n=2117) and 2014 (n=2516). Infections were quantified by highly sensitive qPCR and gametocytes by RT-qPCR. RESULTS: P. falciparum prevalence by qPCR decreased from 42% in 2006 to 9% in 2014. P. vivax prevalence decreased from 42% in 2006 to 13% in 2010, but then increased to 20% in 2014. Parasite densities decreased 5-fold from 2006 to 2010; 72% of P. falciparum and 87% of P. vivax infections were submicroscopic in 2014. Gametocyte density and positivity correlated closely with parasitemia, and population gametocyte prevalence decreased 3-fold for P. falciparum and 29% for P. vivax from 2010 to 2014. CONCLUSIONS: Sustained control has resulted in reduced transmission potential but an increasing proportion of gametocyte carriers are asymptomatic and submicroscopic and represent a challenge to malaria control.

PLoS One
Infection of Ixodes ricinus by Borrelia burgdorferi sensu lato in peri-urban forests of France – Axelle Marchant , Alain Le Coupanec , Claire Joly , Emeline Perthame, Natacha Sertour, Martine Garnier, Vincent Godard, Elisabeth Ferquel , Valerie Choumet

Lyme borreliosis is the most common tick-borne disease in the northern hemisphere. In Europe, it is transmitted by Ixodes ticks that carry bacteria belonging to the Borrelia burgdorferi sensu lato complex. The objective of this work was to explore eco-epidemiological factors of Lyme borreliosis in peri-urban forests of France (Sénart, Notre-Dame and Rambouillet). We investigated whether the introduction of Tamias sibiricus in Sénart could alter the density of infected ticks. Moreover, the density and tick infection were investigated according to the tree species found in various patches of Sénart forest. For this purpose, ticks were sampled during 3 years. In the Sénart forest, the density of nymph and adult ticks showed no significant difference between 2008, 2009 and 2011. The nymph density varied significantly as a function of the month of collection. Regarding the nymphs, a higher rate of infection and infected density were found in 2009. Plots with chipmunks (C) presented a lower density of both nymphs and adult ticks than plots without chipmunks (NC) did. A higher rate of infection of nymphs with Borrelia was seen in C plots. The prevalence of the various species of Borrelia was also found to vary between C and NC plots with the year of the collect. The presence of chestnut trees positively influenced the density of both nymphs and adults. The infected nymph density showed a significant difference depending on the peri-urban forest studied, Sénart being higher than Rambouillet. The prevalence of Borrelia species also differed between the various forests studied. Concerning the putative role that Tamias sibiricus may play in the transmission of Borrelia, our results suggest that its presence is correlated with a higher rate of infection of questing ticks by Borrelia genospecies and if its population increases, it could play a significant role in the risk of transmission of Lyme borreliosis.

Systematic and Applied Microbiology
Manual and expert annotation of the nearly complete genome sequence of Staphylococcus sciuri strain ATCC 29059: A reference for the oxidase-positive staphylococci that supports the atypical phenotypic features of the species group – Christo-Foroux E, Vallaeys T, Loux V, Dassa E, Deutscher J, Wandersman C, Livernois A, Hot C, Criscuolo A, Dauga C, Clermont D, Chesneau O

Staphylococcus sciuri is considered to be one of the most ancestral species in the natural history of the Staphylococcus genus that consists of 48 validly described species. It belongs to the basal group of oxidase-positive and novobiocin-resistant staphylococci that diverged from macrococci approximately 250 million years ago. Contrary to other groups, the S. sciuri species group has not developed host-specific colonization strategies. Genome analysis of S. sciuri ATCC 29059 provides here the first genetic basis for atypical traits that would support the switch between the free-living style and the infective state in animals and humans. From among the most remarkable features, it was noticed in this extensive study that there were a number of phosphoenolpyruvate:carbohydrate phosphotransferase systems (PTS), almost twice as many as any other staphylococci, and the co-occurrence of mevalonate and non-mevalonate pathways for isoprenoid synthesis. The sequenced strain was devoid of the main virulence factors present in Staphylococcus aureus, although it exhibited numerous heme and iron acquisition systems, as well as crt and aldH genes necessary for gold pigment synthesis. The sensing and signaling networks, exemplified by a large and typical repertoire of two-component regulatory systems and a complete panel of master regulators, such as agr, rex, mgrA, rot, sarA and sarR genes, depict the background in which S. aureus virulence genes were later acquired. An additional sigma factor, a distinct set of electron transducer elements and many gene operons similar to those found in Bacillus spp. would constitute the most visible remnant links with Bacillaceae organisms.

Antonie Van Leeuwenhoek
Revisiting the taxonomy of the genus Elizabethkingia using whole-genome sequencing, optical mapping, and MALDI-TOF, along with proposal of three novel Elizabethkingia species: Elizabethkingia bruuniana sp. nov., Elizabethkingia ursingii sp. nov., and Elizabethkingia occulta sp. nov – Nicholson AC, Gulvik CA, Whitney AM, Humrighouse BW, Graziano J, Emery B, Bell M, Loparev V, Juieng P, Gartin J, Bizet C, Clermont D, Criscuolo A, Brisse S, McQuiston JR

The genus Elizabethkingia is genetically heterogeneous, and the phenotypic similarities between recognized species pose challenges in correct identification of clinically derived isolates. In addition to the type species Elizabethkingia meningoseptica, and more recently proposed Elizabethkingia miricola, Elizabethkingia anophelis and Elizabethkingia endophytica, four genomospecies have long been recognized. By comparing historic DNA-DNA hybridization results with whole genome sequences, optical maps, and MALDI-TOF mass spectra on a large and diverse set of strains, we propose a comprehensive taxonomic revision of this genus. Genomospecies 1 and 2 contain the type strains E. anophelis and E. miricola, respectively. Genomospecies 3 and 4 are herein proposed as novel species named as Elizabethkingia bruuniana sp. nov. (type strain, G0146(T) = DSM 2975(T) = CCUG 69503(T) = CIP 111191(T)) and Elizabethkingia ursingii sp. nov. (type strain, G4122(T) = DSM 2974(T) = CCUG 69496(T) = CIP 111192(T)), respectively. Finally, the new species Elizabethkingia occulta sp. nov. (type strain G4070(T) = DSM 2976(T) = CCUG 69505(T) = CIP 111193(T)), is proposed.

International Journal of Systematic and Evolutionary Microbiology
Psychrobacter pasteurii and Psychrobacter piechaudii sp. nov., two novel species within the genus Psychrobacter – Hurtado-Ortiz R, Nazimoudine A, Criscuolo A, Hugon P, Mornico D, Brisse S, Bizet C, Clermont D

Six Gram-negative, non-motile, non-spore-forming, non-pigmented, oxidase- and catalase-positive bacterial strains were deposited in 1972, in the Collection of the Institut Pasteur (CIP), Paris, France. The strains, previously identified as members of the genus Moraxella on the basis of their phenotypic and biochemical characteristics, were placed within the genus Psychrobacter based on the results from comparative 16S rRNA gene sequence studies. Their closest phylogenetic relatives were Psychrobacter sanguinis CIP 110993T, Psychrobacter phenylpyruvicus CIP 82.27T and Psychrobacter lutiphocae CIP 110018T. The DNA G+C contents were between 42.1 and 42.7 mol%. The predominant fatty acids were C18 : 1ω9c, C16 : 0, C12 : 0 3-OH, and C18 : 0. Average nucleotide identity between the six strains and their closest phylogenetic relatives, as well as their phenotypic characteristics, supported the assignment of these strains to two novel species within the genus Psychrobacter. The proposed names for these strains are Psychrobacter pasteurii sp. nov., for which the type strain is A1019T (=CIP 110853T=CECT 9184T), and Psychrobacter piechaudii sp. nov., for which the type strain is 1232T (=CIP110854T=CECT 9185T).

Emerging Infectious Diseases
Real-time whole-genome sequencing for surveillance of Listeria monocytogenes, France. – Moura A, Tourdjman M, Leclercq A, Hamelin E, Laurent E, Fredriksen N, Van Cauteren D, Bracq-Dieye H, Thouvenot P, Vales G, Tessaud-Rita N, Maury MM, Alexandru A, Criscuolo A, Quevillon E, Donguy MP, Enouf V, de Valk H, Brisse S, Lecuit M.

During 2015-2016, we evaluated the performance of whole-genome sequencing (WGS) as a routine typing tool. Its added value for microbiological and epidemiologic surveillance of listeriosis was compared with that for pulsed-field gel electrophoresis (PFGE), the current standard method. A total of 2,743 Listeria monocytogenes isolates collected as part of routine surveillance were characterized in parallel by PFGE and core genome multilocus sequence typing (cgMLST) extracted from WGS. We investigated PFGE and cgMLST clusters containing human isolates. Discrimination of isolates was significantly higher by cgMLST than by PFGE (p<0.001). cgMLST discriminated unrelated isolates that shared identical PFGE profiles and phylogenetically closely related isolates with distinct PFGE profiles. This procedure also refined epidemiologic investigations to include only phylogenetically closely related isolates, improved source identification, and facilitated epidemiologic investigations, enabling identification of more outbreaks at earlier stages. WGS-based typing should replace PFGE as the primary typing method for L. monocytogenes.

Nature Communication
Evolutionary dynamics and genomic features of the Elizabethkingia anophelis 2015 to 2016 Wisconsin outbreak strain. – Perrin A, Larsonneur E, Nicholson AC, Edwards DJ, Gundlach KM, Whitney AM, Gulvik CA, Bell ME, Rendueles O, Cury J, Hugon P, Clermont D, Enouf V, Loparev V, Juieng P, Monson T, Warshauer D, Elbadawi L, Walters MS, Crist MB, Noble-Wang J, Borlaug G, Rocha EPC, Criscuolo A, Touchon M, Davis JP, Holt KE, McQuiston JR, Brisse S

An atypically large outbreak of Elizabethkingia anophelis infections occurred in Wisconsin. Here we show that it was caused by a single strain with thirteen characteristic genomic regions. Strikingly, the outbreak isolates show an accelerated evolutionary rate and an atypical mutational spectrum. Six phylogenetic sub-clusters with distinctive temporal and geographic dynamics are revealed, and their last common ancestor existed approximately one year before the first recognized human infection. Unlike other E. anophelis, the outbreak strain had a disrupted DNA repair mutY gene caused by insertion of an integrative and conjugative element. This genomic change probably contributed to the high evolutionary rate of the outbreak strain and may have increased its adaptability, as many mutations in protein-coding genes occurred during the outbreak. This unique discovery of an outbreak caused by a naturally occurring mutator bacterial pathogen provides a dramatic example of the potential impact of pathogen evolutionary dynamics on infectious disease epidemiology.

Scientific Reports
Proteome remodelling by the stress sigma factor RpoS/σS in Salmonella: identification of small proteins and evidence for post-transcriptional regulation – Lago M, Monteil V, Douche T, Guglielmini J, Criscuolo A, Maufrais C, Matondo M, Norel F

The RpoS/σS sigma subunit of RNA polymerase is the master regulator of the general stress response in many Gram-negative bacteria. Extensive studies have been conducted on σS-regulated gene expression at the transcriptional level. In contrast, very limited information regarding the impact of σS on global protein production is available. In this study, we used a mass spectrometry-based proteomics approach to explore the wide σS-dependent proteome of the human pathogen Salmonella enterica serovar Typhimurium. Our present goals were twofold: (1) to survey the protein changes associated with the ΔrpoS mutation and (2) to assess the coding capacity of σS-dependent small RNAs. Our proteomics data, and complementary assays, unravelled the large impact of σS on the Salmonella proteome, and validated expression and σS regulation of twenty uncharacterized small proteins of 27 to 96 amino acids. Furthermore, a large number of genes regulated at the protein level only were identified, suggesting that post-transcriptional regulation is an important component of the σS response. Novel aspects of σS in the control of important catabolic pathways such as myo-inositol, L-fucose, propanediol, and ethanolamine were illuminated by this work, providing new insights into the physiological remodelling involved in bacterial adaptation to a non-actively growing state.

Pathogens and Disease
Molecular characterization of Chlamydia pneumoniae associated to atherosclerosis. – Yazouli LE, Criscuolo A, Hejaji H, Bouaaza M, Elmdaghri N, Alami AA, Amraoui A, Dakka N, Radouani F

Chlamydia pneumoniae is a respiratory pathogen associated with chronic inflammatory diseases such as asthma and atherosclerosis, and its detection in human carotid and coronary atheroma suggests some support for its involvement in atherogenesis. The main objective of our study was to evaluate the association between Chlamydia pneumoniae and atherosclerosis in Moroccan patients through a case/control approach and detected strain genotyping. A total of 137 cases and 124 controls were enrolled, nested PCR was performed for Chlamydia pneumoniae screening of the peripheral blood mononuclear cells (PBMCs) of both cases and controls as well as atheroma plaques from 37 cases, and positive samples were subjected to sequencing for genotyping and phylogenetic analysis. The results showed 54% and 18%, respectively, for positivity in cases and control PBMCs and 86.5% in atheroma plaques, the difference being significant between the two groups (p<0.001, ORa = 8.580, CI, 95% [3.273-22.491]). Strain sequence analyses showed more than 98% similarity with human reference strains, and revealed various genotypes. This study supports the involvement of Chlamydia pneumoniae in atherosclerosis in the studied population and genotyping revealed that detected strains were identical to human strains circulating worldwide.

Nature Microbiology
N-terminomics identifies Prli42 as a membrane miniprotein conserved in Firmicutes and critical for stressosome activation in Listeria monocytogenes – Impens F, Rolhion N, Radoshevich L, Bécavin C, Duval M, Mellin J, García Del Portillo F, Pucciarelli MG, Williams AH, Cossart P

To adapt to changing environments, bacteria have evolved numerous pathways that activate stress response genes. In Gram-positive bacteria, the stressosome, a cytoplasmic complex, relays external cues and activates the sigma B regulon. The stressosome is structurally well-characterized in Bacillus, but how it senses stress remains elusive. Here, we report a genome-wide N-terminomic approach in Listeria that strikingly led to the discovery of 19 internal translation initiation sites and 6 miniproteins, among which one, Prli42, is conserved in Firmicutes. Prli42 is membrane-anchored and interacts with orthologues of Bacillus stressosome components. We reconstituted the Listeria stressosome in vitro and visualized its supramolecular structure by electron microscopy. Analysis of a series of Prli42 mutants demonstrated that Prli42 is important for sigma B activation, bacterial growth following oxidative stress and for survival in macrophages. Taken together, our N-terminonic approach unveiled Prli42 as a long-sought link between stress and the stressosome.

Frontiers in Cellular and Infectectious Microbiology
Genome-Wide Transcriptional Start Site Mapping and sRNA Identification in the Pathogen Leptospira interrogans – Zhukova A, Fernandes LG, Hugon P, Pappas CJ, Sismeiro O, Coppée JY, Becavin C, Malabat C, Eshghi A, Zhang JJ, Yang FX, Picardeau M

Leptospira are emerging zoonotic pathogens transmitted from animals to humans typically through contaminated environmental sources of water and soil. Regulatory pathways of pathogenic Leptospira spp. underlying the adaptive response to different hosts and environmental conditions remains elusive. In this study, we provide the first global Transcriptional Start Site (TSS) map of a Leptospira species. RNA was obtained from the pathogen Leptospira interrogans grown at 30°C (optimal in vitro temperature) and 37°C (host temperature) and selectively enriched for 5′ ends of native transcripts. A total of 2865 and 2866 primary TSS (pTSS) were predicted in the genome of L. interrogans at 30 and 37°C, respectively. The majority of the pTSSs were located between 0 and 10 nucleotides from the translational start site, suggesting that leaderless transcripts are a common feature of the leptospiral translational landscape. Comparative differential RNA-sequencing (dRNA-seq) analysis revealed conservation of most pTSS at 30 and 37°C. Promoter prediction algorithms allow the identification of the binding sites of the alternative sigma factor sigma 54. However, other motifs were not identified indicating that Leptospira consensus promoter sequences are inherently different from the Escherichia coli model. RNA sequencing also identified 277 and 226 putative small regulatory RNAs (sRNAs) at 30 and 37°C, respectively, including eight validated sRNAs by Northern blots. These results provide the first global view of TSS and the repertoire of sRNAs in L. interrogans. These data will establish a foundation for future experimental work on gene regulation under various environmental conditions including those in the host.

Listeriomics: an Interactive Web Platform for Systems Biology of Listeria. – Bécavin C, Koutero M, Tchitchek N, Cerutti F, Lechat P, Maillet N, Hoede C, Chiapello H, Gaspin C, Cossart P

As for many model organisms, the amount of Listeria omics data produced has recently increased exponentially. There are now >80 published complete Listeria genomes, around 350 different transcriptomic data sets, and 25 proteomic data sets available. The analysis of these data sets through a systems biology approach and the generation of tools for biologists to browse these various data are a challenge for bioinformaticians. We have developed a web-based platform, named Listeriomics, that integrates different tools for omics data analyses, i.e., (i) an interactive genome viewer to display gene expression arrays, tiling arrays, and sequencing data sets along with proteomics and genomics data sets; (ii) an expression and protein atlas that connects every gene, small RNA, antisense RNA, or protein with the most relevant omics data; (iii) a specific tool for exploring protein conservation through the Listeria phylogenomic tree; and (iv) a coexpression network tool for the discovery of potential new regulations. Our platform integrates all the complete Listeria species genomes, transcriptomes, and proteomes published to date. This website allows navigation among all these data sets with enriched metadata in a user-friendly format and can be used as a central database for systems biology analysis. IMPORTANCE In the last decades, Listeria has become a key model organism for the study of host-pathogen interactions, noncoding RNA regulation, and bacterial adaptation to stress. To study these mechanisms, several genomics, transcriptomics, and proteomics data sets have been produced. We have developed Listeriomics, an interactive web platform to browse and correlate these heterogeneous sources of information. Our website will allow listeriologists and microbiologists to decipher key regulation mechanism by using a systems biology approach.

Scientific Reports
Extensive transcriptome analysis correlates the plasticity of Entamoeba histolytica pathogenesis to rapid phenotype changes depending on the environment. – Christian Weber, Mikael Koutero, Marie-Agnes Dillies, Hugo Varet, Cesar Lopez-Camarillo, Jean Yves Coppée, Chung-Chau Hon, Nancy Guillén

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

Promyelocytic Leukemia Protein (PML) Controls Listeria monocytogenes Infection – Ribet D, Lallemand-Breitenbach V, Ferhi O, Nahori M-A, Varet H, de Thé H, Cossart P.

The promyelocytic leukemia protein (PML) is the main organizer of stress-responsive subnuclear structures called PML nuclear bodies. These structures recruit multiple interactors and modulate their abundance or their posttranslational modifications, notably by the SUMO ubiquitin-like modifiers. The involvement of PML in antiviral responses is well established. In contrast, the role of PML in bacterial infection remains poorly characterized. Here, we show that PML restricts infection by the pathogenic bacterium Listeria monocytogenes but not by Salmonella enterica serovar Typhimurium. During infection, PML undergoes oxidation-mediated multimerization, associates with the nuclear matrix, and becomes de-SUMOylated due to the pore-forming activity of the Listeria toxin listeriolysin O (LLO). These events trigger an antibacterial response that is not observed during in vitro infection by an LLO-defective Listeria mutant, but which can be phenocopied by specific induction of PML de-SUMOylation. Using transcriptomic and proteomic microarrays, we also characterized a network of immunity genes and cytokines, which are regulated by PML in response to Listeria infection but independently from the listeriolysin O toxin. Our study thus highlights two mechanistically distinct complementary roles of PML in host responses against bacterial infection. The promyelocytic leukemia protein (PML) is a eukaryotic protein that can polymerize in discrete nuclear assemblies known as PML nuclear bodies (NBs) and plays essential roles in many different cellular processes. Key to its function, PML can be posttranslationally modified by SUMO, a ubiquitin-like modifier. Identification of the role of PML in antiviral defenses has been deeply documented. In contrast, the role of PML in antibacterial defenses remains elusive. Here, we identify two mechanistically distinct complementary roles of PML in antibacterial responses against pathogens such as Listeria: (i) we show that PML regulates the expression of immunity genes in response to bacterial infection, and (ii) we unveil the fact that modification of PML SUMOylation by bacterial pore-forming toxins is sensed as a danger signal, leading to a restriction of bacterial intracellular multiplication. Taken together, our data reinforce the concept that intranuclear bodies can dynamically regulate important processes, such as defense against invaders.

Nature Microbioloby
Whole genome-based population biology and epidemiological surveillance of Listeria monocytogenes – Alexixandra Moura, Alexis Criscuolo, Hannes Pouseele, Mylène M. Maury, Alexandre Leclercq, Cheryl Tarr, Jonas T. Björkman, Timothy Dallman, Aleisha Reimer, Vincent Enouf, Elise Larsonneur, Heather Carleton, Hélène Bracq-Dieye, Lee S. Katz, Louis Jones, Marie Touchon, Mathieu Tourdjman, Matthew Walker, Steven Stroika, Thomas Cantinelli, Viviane Chenal-Francisque, Zuzana Kucerova, Eduardo P. C. Rocha, Celine Nadon, Kathie Grant, Eva M. Nielsen, Bruno Pot, Peter Gerner-Smidt, Marc Lecuit, Sylvain Brisse

Listeria monocytogenes (Lm) is a major human foodborne pathogen. Numerous Lm outbreaks have been reported worldwide and associated with a high case fatality rate, reinforcing the need for strongly coordinated surveillance and outbreak control. We developed a universally applicable genome-wide strain genotyping approach and investigated the population diversity of Lm using 1,696 isolates from diverse sources and geographical locations. We define, with unprecedented precision, the population structure of Lm, demonstrate the occurrence of international circulation of strains and reveal the extent of heterogeneity in virulence and stress resistance genomic features among clinical and food isolates. Using historical isolates, we show that the evolutionary rate of Lm from lineage I and lineage II is low (∼2.5 × 10−7 substitutions per site per year, as inferred from the core genome) and that major sublineages (corresponding to so-called ‘epidemic clones’) are estimated to be at least 50–150 years old. This work demonstrates the urgent need to monitor Lm strains at the global level and provides the unified approach needed for global harmonization of Lm genome-based typing and population biology.

Concomitant emergence of the antisense protein gene of HIV-1 and of the pandemic – Elodie Cassan, Anne-Muriel Arigon-Chifolleau, Jean-Michel Mesnard, Antoine Gross, and Olivier Gascuel

HIV-1 is commonly assumed to have nine genes. However, in 1988 a 10th gene was suggested, overlapped by the env gene, but read on the antisense strand. The corresponding protein was named AntiSense Protein (ASP). Several pieces of evidence argue in favor of ASP expression in vivo, but its function is still unknown. We performed the first evolutionary study of ASP, using a very large number of HIV-1 and SIV (simian) sequences. Our results show that ASP is specific to group M of HIV-1, which is responsible for the pandemic. Moreover, we demonstrated that evolutionary forces act to maintain the asp gene within the M sequences and showed a striking correlation of asp with the spread of the pandemic.

Phylogenomic analysis supports the ancestral presence of LPS-outer membranes in the Firmicutes. – Antunes LC, Poppleton D, Klingl A, Criscuolo A, Dupuy B, Brochier-Armanet C, Beloin C, Gribaldo S

One of the major unanswered questions in evolutionary biology is when and how the transition between diderm (two membranes) and monoderm (one membrane) cell envelopes occurred in Bacteria. The Negativicutes and the Halanaerobiales belong to the classically monoderm Firmicutes, but possess outer membranes with lipopolysaccharide (LPS-OM). Here, we show that they form two phylogenetically distinct lineages, each close to different monoderm relatives. In contrast, their core LPS biosynthesis enzymes were inherited vertically, as in the majority of bacterial phyla. Finally, annotation of key OM systems in the Halanaerobiales and the Negativicutes shows a puzzling combination of monoderm and diderm features. Together, these results support the hypothesis that the LPS-OMs of Negativicutes and Halanaerobiales are remnants of an ancient diderm cell envelope that was present in the ancestor of the Firmicutes, and that the monoderm phenotype in this phylum is a derived character that arose multiple times independently through OM loss.

MEMHDX: An interactive tool to expedite the statistical validation and visualization of large HDX-MS datasets – Hourdel V, Volant S, O’Brien DP, Chenal A, Chamot-Rooke J, Dillies MA, Brier S

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: MEMHDX is freely available as a web tool at the project home page CONTACT:; SUPPLEMENTARY INFORMATION: Supplementary data is available at Bioinformatics online.

Key experimental evidence of chromosomal DNA transfer among selected tuberculosis-causing mycobacteria – Boritsch EC, Khanna V, Pawlik A, Honoré N, Navas VH, Ma L, Bouchier C, Seemann T, Supply P, Stinear TP, Brosch R

Horizontal gene transfer (HGT) is a major driving force of bacterial diversification and evolution. For tuberculosis-causing mycobacteria, the impact of HGT in the emergence and distribution of dominant lineages remains a matter of debate. Here, by using fluorescence-assisted mating assays and whole genome sequencing, we present unique experimental evidence of chromosomal DNA transfer between tubercle bacilli of the early-branching Mycobacterium canettii clade. We found that the obtained recombinants had received multiple donor-derived DNA fragments in the size range of 100 bp to 118 kbp, fragments large enough to contain whole operons. Although the transfer frequency between M. canettii strains was low and no transfer could be observed among classical Mycobacterium tuberculosis complex (MTBC) strains, our study provides the proof of concept for genetic exchange in tubercle bacilli. This outstanding, now experimentally validated phenomenon presumably played a key role in the early evolution of the MTBC toward pathogenicity. Moreover, our findings also provide important information for the risk evaluation of potential transfer of drug resistance and fitness mutations among clinically relevant mycobacterial strains.

Genome biology and evolution
Species- and strain-specific adaptation of the HSP70 super family in pathogenic Trypanosomatids – Drini S, Criscuolo A, Lechat P, Imamura H, Skalický T, Rachidi N, Lukeš J, Dujardin JC, Späth GF

All eukaryotic genomes encode multiple members of the heat shock protein 70 (HSP70) family, which evolved distinctive structural and functional features in response to specific environmental constraints. Phylogenetic analysis of this protein family thus can inform on genetic and molecular mechanisms that drive species-specific environmental adaptation. Here we use the eukaryotic pathogen Leishmania spp. as a model system to investigate the evolution of the HSP70 protein family in an early-branching eukaryote that is prone to gene amplification and adapts to cytotoxic host environments by stress-induced and chaperone-dependent stage differentiation. Combining phylogenetic and comparative analyses of trypanosomatid genomes, draft genome of Paratrypanosoma and recently published genome sequences of 204 L. donovani field isolates, we gained unique insight into the evolutionary dynamics of the Leishmania HSP70 protein family. We provide evidence for (i) significant evolutionary expansion of this protein family in Leishmania through gene amplification and functional specialization of highly conserved canonical HSP70 members, (ii) evolution of trypanosomatid-specific, non-canonical family members that likely gained ATPase-independent functions, and (iii) loss of one atypical HSP70 member in the Trypanosoma genus. Finally, we reveal considerable copy number variation of canonical cytoplasmic HSP70 in highly related L. donovani field isolates, thus identifying this locus as a potential hot spot of environment-genotype interaction. Our data draw a complex picture of the genetic history of HSP70 in trypanosomatids that is driven by the remarkable plasticity of the Leishmania genome to undergo massive intra-chromosomal gene amplification to compensate for the absence of regulated transcriptional control in these parasites.

Scientific reports
Genomic epidemiology and global diversity of the emerging bacterial pathogen Elizabethkingia anophelis – Breurec S, Criscuolo A, Diancourt L, Rendueles O, Vandenbogaert M, Passet V, Caro V, Rocha EP, Touchon M, Brisse S

Elizabethkingia anophelis is an emerging pathogen involved in human infections and outbreaks in distinct world regions. We investigated the phylogenetic relationships and pathogenesis-associated genomic features of two neonatal meningitis isolates isolated 5 years apart from one hospital in Central African Republic and compared them with Elizabethkingia from other regions and sources. Average nucleotide identity firmly confirmed that E. anophelis, E. meningoseptica and E. miricola represent demarcated genomic species. A core genome multilocus sequence typing scheme, broadly applicable to Elizabethkingia species, was developed and made publicly available ( Phylogenetic analysis revealed distinct E. anophelis sublineages and demonstrated high genetic relatedness between the African isolates, compatible with persistence of the strain in the hospital environment. CRISPR spacer variation between the African isolates was mirrored by the presence of a large mobile genetic element. The pan-genome of E. anophelis comprised 6,880 gene families, underlining genomic heterogeneity of this species. African isolates carried unique resistance genes acquired by horizontal transfer. We demonstrated the presence of extensive variation of the capsular polysaccharide synthesis gene cluster in E. anophelis. Our results demonstrate the dynamic evolution of this emerging pathogen and the power of genomic approaches for Elizabethkingia identification, population biology and epidemiology.

J Antimicrob Chemother
Bacteriophage LM33_P1, a fast-acting weapon against the pandemic ST131-O25b:H4 Escherichia coli clonal complex. – Dufour N, Clermont O, La Combe B, Messika J, Dion S, Khanna V, Denamur E, Ricard JD, Debarbieux L

Amongst the highly diverse Escherichia coli population, the ST131-O25b:H4 clonal complex is particularly worrisome as it is associated with a high level of antibiotic resistance. The lack of new antibiotics, the worldwide continuous increase of infections caused by MDR bacteria and the need for narrow-spectrum antimicrobial agents have revived interest in phage therapy. In this article, we describe a virulent bacteriophage, LM33_P1, which specifically infects O25b strains, and provide data related to its therapeutic potential. A large panel of E. coli strains (n = 283) was used to assess both the specificity and the activity of bacteriophage LM33_P1. Immunology, biochemistry and genetics-based methods confirmed this specificity. Virology methods and sequencing were used to characterize this bacteriophage in vitro, while three relevant mouse models were employed to show its in vivo efficacy. Bacteriophage LM33_P1 exclusively infects O25b E. coli strains with a 70% coverage on sequence types associated with high antibiotic resistance (ST131 and ST69). This specificity is due to an interaction with the LPS mediated by an original tail fibre. LM33_P1 also has exceptional intrinsic properties with a high adsorption constant and produces over 300 virions per cell in <10 min. Using animal pneumonia, septicaemia and urinary tract infection models, we showed the in vivo efficacy of LM33_P1 to reduce the bacterial load in several organs. Bacteriophage LM33_P1 represents the first weapon that specifically and quickly kills O25b E. coli strains. Therapeutic approaches derived from this bacteriophage could be developed to stop or slow down the spread of the ST131-O25b:H4 drug-resistant clonal complex in humans.

Nucleic Acids Research
DNA minicircles clarify the specific role of DNA structure on retroviral integration – Pasi Marco, Mornico Damien, Volant Stevenn, Juchet Anna, Batisse Julien, Bouchier Christiane, Parissi Vincent, Ruff Marc, Lavery Richard, Lavigne Marc

Chromatin regulates the selectivity of retroviral integration into the genome of infected cells. At the nucleosome level, both histones and DNA structure are involved in this regulation. We propose a strategy that allows to specifically study a single factor: the DNA distortion induced by the nucleosome. This strategy relies on mimicking this distortion using DNA minicircles (MCs) having a fixed rotational orientation of DNA curvature, coupled with atomic-resolution modelling. Contrasting MCs with linear DNA fragments having identical sequences enabled us to analyse the impact of DNA distortion on the efficiency and selectivity of integration. We observed a global enhancement of HIV-1 integration in minicircles and an enrichment of integration sites in the outward-facing DNA major grooves. Both of these changes are favoured by LEDGF/p75, revealing a new, histone-independent role of this integration cofactor. PFV integration is also enhanced in MCs, but is not associated with a periodic redistribution of integration sites, thus highlighting its distinct catalytic properties. MCs help to separate the roles of target DNA structure, histone modifications and integrase cofactors during retroviral integration and to reveal integrase-specific regulation mechanisms.

Cell Microbiol
GH16 and GH81 family β-(1,3)-glucanases in Aspergillus fumigatus are essential for conidial cell wall morphogenesis – Mouyna I, Aimanianda V, Hartl L, Prevost MC, Sismeiro O, Dillies MA, Jagla B, Legendre R, Coppée JY, Latgé JP.

The fungal cell wall is a rigid structure due to 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 A. fumigatus genome coding for potential endo-β-(1,3)-glucanase. ENG1, (previously characterized and named ENGL1, Mouyna et al., 2002), 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.

SARTools: A DESeq2- and EdgeR-Based R Pipeline for Comprehensive Differential Analysis of RNA-Seq Data – Hugo Varet, Loraine Brillet-Guégen, Jean-Yves Coppée, Marie-Agnès Dillies

Several R packages exist for the detection of differentially expressed genes from RNA-Seq data. The analysis process includes three main steps, namely normalization, dispersion estimation and test for differential expression. Quality control steps along this process are recommended but not mandatory, and failing to check the characteristics of the dataset may lead to spurious results. In addition, normalization methods and statistical models are not exchangeable across the packages without adequate transformations the users are often not aware of. Thus, dedicated analysis pipelines are needed to include systematic quality control steps and prevent errors from misusing the proposed methods. SARTools is an R pipeline for differential analysis of RNA-Seq count data. It can handle designs involving two or more conditions of a single biological factor with or without a blocking factor (such as a batch effect or a sample pairing). It is based on DESeq2 and edgeR and is composed of an R package and two R script templates (for DESeq2 and edgeR respectively). Tuning a small number of parameters and executing one of the R scripts, users have access to the full results of the analysis, including lists of differentially expressed genes and a HTML report that (i) displays diagnostic plots for quality control and model hypotheses checking and (ii) keeps track of the whole analysis process, parameter values and versions of the R packages used. SARTools provides systematic quality controls of the dataset as well as diagnostic plots that help to tune the model parameters. It gives access to the main parameters of DESeq2 and edgeR and prevents untrained users from misusing some functionalities of both packages. By keeping track of all the parameters of the analysis process it fits the requirements of reproducible research.

Bacteriocin from epidemic Listeria strains alters the host intestinal microbiota to favor infection – Juan J. Quereda, Olivier Dussurget, Marie-Anne Nahori, Amine Ghozlane, Stevenn Volant, Marie-Agnès Dillies, Béatrice Regnault, Sean Kennedy, Stanislas Mondot, Barbara Villoing, Pascale Cossart, and Javier Pizarro-Cerda

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.

Genome Research
Improved definition of the mouse transcriptome via targeted RNA sequencing – Giovanni Bussotti, Tommaso Leonardi, Michael B. Clark, Tim R. Mercer, Joanna Crawford, Lorenzo Malquori, Cedric Notredame, Marcel E. Dinger, John S. Mattick and Anton J. Enright

Targeted RNA sequencing (CaptureSeq) uses oligonucleotide probes to capture RNAs for sequencing, providing enriched read coverage, accurate measurement of gene expression, and quantitative expression data. We applied CaptureSeq to refine transcript annotations in the current murine GRCm38 assembly. More than 23,000 regions corresponding to putative or annotated long noncoding RNAs (lncRNAs) and 154,281 known splicing junction sites were selected for targeted sequencing across five mouse tissues and three brain subregions. The results illustrate that the mouse transcriptome is considerably more complex than previously thought. We assemble more complete transcript isoforms than GENCODE, expand transcript boundaries, and connect interspersed islands of mapped reads. We describe a novel filtering pipeline that identifies previously unannotated but high-quality transcript isoforms. In this set, 911 GENCODE neighboring genes are condensed into 400 expanded gene models. Additionally, 594 GENCODE lncRNAs acquire an open reading frame (ORF) when their structure is extended with CaptureSeq. Finally, we validate our observations using current FANTOM and Mouse ENCODE resources.

F1000 Research
DREAMTools: a Python package for scoring collaborative challenges – Thomas Cokelaer, Mukesh Bansal, Christopher Bare, Erhan Bilal, Brian M. Bot, Elias Chaibub Neto, Federica Eduati, Alberto de la Fuente, Mehmet Gönen, Steven M. Hill, Bruce Hoff, Jonathan R. Karr, Robert Küffner, Michael P. Menden, Pablo Meyer, Raquel Norel, Abhishek Pratap, Robert J. Prill, Matthew T. Weirauch, James C. Costello, Gustavo Stolovitzky, Julio Saez-Rodriguez

DREAM challenges are community competitions designed to advance computational methods and address fundamental questions in system biology and translational medicine. Each challenge asks participants to develop and apply computational methods to either predict unobserved outcomes or to identify unknown model parameters given a set of training data. Computational methods are evaluated using an automated scoring metric, scores are posted to a public leaderboard, and methods are published to facilitate community discussions on how to build improved methods. By engaging participants from a wide range of science and engineering backgrounds, DREAM challenges can comparatively evaluate a wide range of statistical, machine learning, and biophysical methods. Here, we describe DREAMTools, a Python package for evaluating DREAM challenge scoring metrics. DREAMTools provides a command line interface that enables researchers to test new methods on past challenges, as well as a framework for scoring new challenges. As of March 2016, DREAMTools includes more than 80% of completed DREAM challenges. DREAMTools complements the data, metadata, and software tools available at the DREAM website and on the Synapse platform at

Inferring causal molecular networks: empirical assessment through a community-based effort – Steven M Hill, Laura M Heiser, Thomas Cokelaer, Michael Unger, Nicole K Nesser, Daniel E Carlin, Yang Zhang , Artem Sokolov , Evan O Paull, Chris K Wong, Kiley Graim, Adrian Bivol , Haizhou Wang, Fan Zhu, Bahman Afsari, Ludmila V Danilova, Alexander V Favorov, Wai Shing Lee, Dane Taylor, Chenyue W Hu, Byron L Long, David P Noren, Alexander J Bisberg , HPN-DREAM Consortium, Gordon B Mills, Joe W Gray, Michael Kellen, Thea Norman, Stephen Friend, Amina A Qutub, Elana J Fertig, Yuanfang Guan, Mingzhou Song, Joshua M Stuart, Paul T Spellman, Heinz Koeppl, Gustavo Stolovitzky, Julio Saez-Rodriguez & Sach Mukherjee

It remains unclear whether causal, rather than merely correlational, relationships in molecular networks can be inferred in complex biological settings. Here we describe the HPN-DREAM network inference challenge, which focused on learning causal influences in signaling networks. We used phosphoprotein data from cancer cell lines as well as in silico data from a nonlinear dynamical model. Using the phosphoprotein data, we scored more than 2,000 networks submitted by challenge participants. The networks spanned 32 biological contexts and were scored in terms of causal validity with respect to unseen interventional data. A number of approaches were effective, and incorporating known biology was generally advantageous. Additional sub-challenges considered time-course prediction and visualization. Our results suggest that learning causal relationships may be feasible in complex settings such as disease states. Furthermore, our scoring approach provides a practical way to empirically assess inferred molecular networks in a causal sense

Environmental Microbiology
Bacteriophages to redu ce gut carriage of antibiotic resistant uropathogens with low impact on microbiota composition – Matthieu Galtier, Luisa De Sordi, Damien Maura, Harindra Arachchi, Stevenn Volant, Marie-Agnès Dillies and Laurent Debarbieux

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.

BMC Biology
Genome-wide replication landscape of Candida glabrata – Descorps-Declère S., Saguez C., Cournac A., Marbouty M., Rolland T., Ma L., Bouchier C., Moszer I., Dujon B., Koszul R., Richard GF.

BACKGROUND: The opportunistic pathogen Candida glabrata is a member of the Saccharomycetaceae yeasts. Like its close relative Saccharomyces cerevisiae, it underwent a whole-genome duplication followed by an extensive loss of genes. Its genome contains a large number of very long tandem repeats, called megasatellites. In order to determine the whole replication program of the C. glabrata genome and its general chromosomal organization, we used deep-sequencing and chromosome conformation capture experiments. RESULTS: We identified 253 replication fork origins, genome wide. Centromeres, HML and HMR loci, and most histone genes are replicated early, whereas natural chromosomal breakpoints are located in late-replicating regions. In addition, 275 autonomously replicating sequences (ARS) were identified during ARS-capture experiments, and their relative fitness was determined during growth competition. Analysis of ARSs allowed us to identify a 17-bp consensus, similar to the S. cerevisiae ARS consensus sequence but slightly more constrained. Megasatellites are not in close proximity to replication origins or termini. Using chromosome conformation capture, we also show that early origins tend to cluster whereas non-subtelomeric megasatellites do not cluster in the yeast nucleus. CONCLUSIONS: Despite a shorter cell cycle, the C. glabrata replication program shares unexpected striking similarities to S. cerevisiae, in spite of their large evolutionary distance and the presence of highly repetitive large tandem repeats in C. glabrata. No correlation could be found between the replication program and megasatellites, suggesting that their formation and propagation might not be directly caused by replication fork initiation or termination.

Nature Genetics
Uncovering Listeria monocytogenes hypervirulence by harnessing its biodiversity – Mylène M Maury, Yu-Huan Tsai, Caroline Charlier, Marie Touchon, Viviane Chenal-Francisque, Alexandre Leclercq, Alexis Criscuolo, Charlotte Gaultier, Sophie Roussel, Anne Brisabois, Olivier Disson, Eduardo PC Rocha, Sylvain Brisse, Marc Lecuit

Microbial pathogenesis studies are typically performed with reference strains, thereby overlooking within-species heterogeneity in microbial virulence. Here we integrated human epidemiological and clinical data with bacterial population genomics to harness the biodiversity of the model foodborne pathogen Listeria monocytogenes and decipher the basis of its neural and placental tropisms. Taking advantage of the clonal structure of this bacterial species, we identify clones epidemiologically associated either with food or with human central nervous system (CNS) or maternal-neonatal (MN) listeriosis. The latter clones are also most prevalent in patients without immunosuppressive comorbidities. Strikingly, CNS- and MN-associated clones are hypervirulent in a humanized mouse model of listeriosis. By integrating epidemiological data and comparative genomics, we have uncovered multiple new putative virulence factors and demonstrate experimentally the contribution of the first gene cluster mediating L. monocytogenes neural and placental tropisms. This study illustrates the exceptional power in harnessing microbial biodiversity to identify clinically relevant microbial virulence attributes.

Nucleic Acids Research
topIb, a phylogenetic hallmark gene of Thaumarchaeota encodes a functional eukaryote-like topoisomerase IB. – Narimane Dahmane, Danièle Gadelle, Stéphane Delmas, Alexis Criscuolo, Stephan Eberhard, Nicole Desnoues, Sylvie Collin, Hongliang Zhang, Yves Pommier, Patrick Forterre, Guennadi Sezonov

Type IB DNA topoisomerases can eliminate torsional stresses produced during replication and transcription. These enzymes are found in all eukaryotes and a short version is present in some bacteria and viruses. Among prokaryotes, the long eukaryotic version is only observed in archaea of the phylum Thaumarchaeota. However, the activities and the roles of these topoisomerases have remained an open question. Here, we demonstrate that all available thaumarchaeal genomes contain a topoisomerase IB gene that defines a monophyletic group closely related to the eukaryotic enzymes. We show that the topIB gene is expressed in the model thaumarchaeon Nitrososphaera viennensis and we purified the recombinant enzyme from the uncultivated thaumarchaeon Candidatus Caldiarchaeum subterraneum. This enzyme is active in vitro at high temperature, making it the first thermophilic topoisomerase IB characterized so far. We have compared this archaeal type IB enzyme to its human mitochondrial and nuclear counterparts. The archaeal enzyme relaxes both negatively and positively supercoiled DNA like the eukaryotic enzymes. However, its pattern of DNA cleavage specificity is different and it is resistant to camptothecins (CPTs) and non-CPT Top1 inhibitors, LMP744 and lamellarin D. This newly described thermostable topoisomerases IB should be a promising new model for evolutionary, mechanistic and structural studies.

G3: Genes, Genomes, Genetics
Massive Amplification at an Unselected Locus Accompanies Complex Chromosomal Rearrangements in Yeast – Agnès Thierry, Varun Khanna, Bernard Dujon

Gene amplification has been observed in different organisms in response to environmental constraints such as limited nutrients or exposure to a variety of toxic compounds, conferring them specific phenotypic adaptations by increased expression levels. But the presence of multiple gene copies has generally not been found in natural genomes in absence of specific functional selection. Here we show that the massive amplification of a chromosomal locus (up to 880 copies per cell) occurs in absence of any direct selection, associated with low-order amplifications of flanking segments in complex chromosomal alterations. These results were obtained in the mutants with restored phenotypes that spontaneously appeared from genetically engineered strains of the yeast Saccharomyces cerevisiae suffering from severe fitness reduction. Grossly extended chromosomes (macrotene) were formed, with complex structural alterations but sufficient stability to propagate unchanged over successive generations. Their detailed molecular analysis, including complete genome sequencing, identification of sequence breakpoints and comparisons between mutants revealed novel mechanisms to their formation whose combined action underlies the astonishing dynamics of eukaryotic chromosomes and its consequences.

Scientific Reports
Identification of protein secretion systems in bacterial genomes – Sophie S Abby, Jean Cury, Julien Guglielmini, Bertrand Néron, Marie Touchon, Eduardo PC Rocha

Bacteria with two cell membranes (diderms) have evolved complex systems for protein secretion. These systems were extensively studied in some model bacteria, but the characterisation of their diversity has lagged behind due to lack of standard annotation tools. We built models for accurate identification of protein secretion systems and related appendages in bacteria with LPS-containing outer membranes. They can be used with MacSyFinder (standalone program) or online ( They include protein profiles and information on the system’s composition and genetic organisation. They can be used to search for T1SS-T6SS, T9SS, and accessorily for flagella, Type IV and Tad pili. We identified ~10,000 systems in bacterial genomes, where T1SS and T5SS were by far the most abundant and widespread. The recently described T6SSiii and T9SS were restricted to Bacteroidetes, and T6SSii to Francisella. T2SS, T3SS, and T4SS were frequently encoded in single-copy in one locus, whereas most T1SS were encoded in two loci. The secretion systems of diderm Firmicutes were similar to those found in other diderms. Novel systems may remain to be discovered, since some clades of environmental bacteria lacked all known protein secretion systems. Our models can be fully customized, which should facilitate the identification of novel systems.