Primary transcriptome analysis reveals importance of IS elements for the shaping of the transcriptional landscape of Bordetella pertussis

Combined transcriptomic and ChIPseq analyses of the Bordetella pertussis RisA regulon

The regulation of Bordetella pertussis virulence is mediated by the two-component system BvgA/S, which activates the transcription of virulence-activated genes (vags). In the avirulent phase, the vags are not expressed, but instead, virulence-repressed genes (vrgs) are expressed, under the control of another two-component system, RisA/K. Here, we combined transcriptomic and chromatin immunoprecipitation sequencing (ChIPseq) data to examine the RisA/K regulon. We performed RNAseq analyses of RisA-deficient and RisA-phosphoablative B. pertussis mutants cultivated in virulent and avirulent conditions. We confirmed that the expression of most vrgs is regulated by phosphorylated RisA. However, the expression of some, including those involved in flagellum biosynthesis and chemotaxis, requires RisA independently of phosphorylation. Many RisA-regulated genes encode proteins with regulatory functions, suggesting multiple RisA regulation cascades. By ChIPseq analyses, we identified 430 RisA-binding sites, 208 within promoter regions, 201 within open reading frames, and 21 in non-coding regions. RisA binding was demonstrated in the promoter regions of most vrgs and, surprisingly, of some vags, as well as for other genes not identified as vags or vrgs. Unexpectedly, many genes, including some vags, like prn, brpL, bipA, and cyaA, contain a BvgA-binding site and a RisA-binding site, which increases the complexity of the RisAK/BvgAS network in B. pertussis virulence regulation.IMPORTANCEThe expression of virulence-activated genes (vags) of Bordetella pertussis, the etiological agent of whooping cough, is under the transcriptional control of the two-component system BvgA/S, which allows the bacterium to switch between virulent and avirulent phases. In addition, the more recently identified two-component system RisA/K is required for the expression of B. pertussis genes, collectively named vrgs, that are repressed during the virulent phase but activated during the avirulent phase. We have characterized the RisA/K regulon by combined transcriptomic and chromatin immunoprecipitation sequencing analyses. We identified more than 400 RisA-binding sites. Many of them are localized in promoter regions, especially vrgs, but some were found within open reading frames and in non-coding regions. Surprisingly, RisA-binding sites were also found in promoter regions of some vags, illustrating the previously underappreciated complexity of virulence regulation in B. pertussis.

Investigating the impact of insertion sequences and transposons in the genomes of the most significant phytopathogenic bacteria

Genetic variability in phytopathogens is one of the main problems encountered for effective plant disease control. This fact may be related to the presence of transposable elements (TEs), but little is known about their role in host genomes. Here, we performed the most comprehensive analysis of insertion sequences (ISs) and transposons (Tns) in the genomes of the most important bacterial plant pathogens. A total of 35 692 ISs and 71 transposons were identified in 270 complete genomes. The level of pathogen-host specialization was found to be a significant determinant of the element distribution among the species. Some Tns were identified as carrying virulence factors, such as genes encoding effector proteins of the type III secretion system and resistance genes for the antimicrobial streptomycin. Evidence for IS-mediated ectopic recombination was identified in Xanthomonas genomes. Moreover, we found that IS elements tend to be inserted in regions near virulence and fitness genes, such ISs disrupting avirulence genes in X. oryzae genomes. In addition, transcriptome analysis under different stress conditions revealed differences in the expression of genes encoding transposases in the Ralstonia solanacearum, X. oryzae, and P. syringae species. Lastly, we also investigated the role of Tns in regulation via small noncoding regulatory RNAs and found these elements may target plant-cell transcriptional activators. Taken together, the results indicate that TEs may have a fundamental role in variability and virulence in plant pathogenic bacteria.

Interference of flagellar rotation up-regulates the expression of small RNA contributing to Bordetella pertussis infection

Bacterial small RNAs (sRNAs) posttranscriptionally regulate gene expressions involved in various biological processes, including pathogenicity. Our previous study identified sRNAs, the expression of which was up-regulated in Bordetella pertussis, the causative agent of whooping cough, upon tracheal colonization of the bacteria; however, their roles in bacterial infection remain unknown. Here, we found that one sRNA, Bpr4, contributes to B. pertussis infection by posttranscriptionally up-regulating filamentous hemagglutinin (FHA), a major adhesin of the bacteria. Bpr4 bound to the 5' untranslated region of fhaB mRNA encoding FHA and inhibited its degradation mediated by RNaseE. Our results demonstrated that Bpr4 up-regulation was triggered by the interference of flagellar rotation, which caused the disengagement of MotA, a flagellar stator. Subsequently, MotA activated a diguanylate cyclase to generate cyclic di-GMP, which plays a role in Bpr4 up-regulation through the RisK/RisA two-component system. Our findings indicate that a flagellum-triggered sensory system contributes to B. pertussis infection.

Insertion Sequence (IS)-Excision Enhancer (IEE)-Mediated IS Excision from the lacZ Gene Restores the Lactose Utilization Defect of Shiga Toxin-Producing Escherichia coli O121:H19 Strains and Is Responsible for Their Delayed Lactose Utilization Phenotype

Lactose utilization is one of the general biochemical characteristics of Escherichia coli, and the lac operon is responsible for this phenotype, which can be detected on lactose-containing media, such as MacConkey agar, after 24 h of incubation. However, some Shiga toxin-producing E. coli (STEC) O121:H19 strains exhibit an unusual phenotype called delayed lactose utilization (DLU), in which lactose utilization can be detected after 48 h of cultivation but not after only 24 h of cultivation. Insertion of an insertion sequence (IS), IS600, into the lacZ gene appears to be responsible for the DLU phenotype, and exposure to lactose has been reported to be necessary to observe this phenotype, but the mechanism underlying these phenomena remains to be elucidated. Here, we performed detailed analyses of the lactose utilization abilities of a set of O121:H19 strains and their mutants and found that IS-excision enhancer (IEE)-mediated excision of IS600 reactivates the lacZ gene and that the selective proliferation of IS-cured subclones in lactose-supplemented culture medium is responsible for the expression of the DLU phenotype. In addition, we analyzed the patterns of IS insertion into the lacZ and iee genes in the global O121:H19 population and revealed that while there are O121:H19 strains or lineage/sublineages that contain the IS insertion into iee or intact lacZ and thus do not show the DLU phenotype, most currently circulating O121:H19 strains contain IS600-inserted lacZ and intact iee and thus exhibit this phenotype.

IMPORTANCE Insertion sequences (ISs) can modulate gene expression by gene inactivation or activation. While phenotypic changes due to IS insertion/transposition are frequently observed, gene reactivation by precise or simple IS excision rarely occurs. In this study, we show that IS600 is excised from the lacZ gene by IS-excision enhancer (IEE) during the cultivation of Shiga toxin-producing Escherichia coli (STEC) O121:H19 strains that show an unusual phenotype called delayed lactose utilization (DLU). This excision rescued their lactose utilization defect, and the subsequent selective proliferation of IS-cured subclones in lactose-containing medium resulted in the expression of the DLU phenotype. As we also show that most currently circulating O121:H19 strains exhibit this phenotype, this study not only provides information helpful for the isolation and identification of O121:H19 STEC but also offers novel insights into the roles of IS and IEE in the generation of phenotypic variation in bacterial populations.

Towards comprehensive understanding of bacterial genetic diversity: large-scale amplifications in Bordetella pertussis and Mycobacterium tuberculosis

Bacterial genetic diversity is often described solely using base-pair changes despite a wide variety of other mutation types likely being major contributors. Tandem duplication/amplifications are thought to be widespread among bacteria but due to their often-intractable size and instability, comprehensive studies of these mutations are rare. We define a methodology to investigate amplifications in bacterial genomes based on read depth of genome sequence data as a proxy for copy number. We demonstrate the approach with Bordetella pertussis, whose insertion sequence element-rich genome provides extensive scope for amplifications to occur. Analysis of data for 2430 B. pertussis isolates identified 272 putative amplifications, of which 94 % were located at 11 hotspot loci. We demonstrate limited phylogenetic connection for the occurrence of amplifications, suggesting unstable and sporadic characteristics. Genome instability was further described in vitro using long-read sequencing via the Nanopore platform, which revealed that clonally derived laboratory cultures produced heterogenous populations rapidly. We extended this research to analyse a population of 1000 isolates of another important pathogen, Mycobacterium tuberculosis. We found 590 amplifications in M. tuberculosis, and like B. pertussis, these occurred primarily at hotspots. Genes amplified in B. pertussis include those involved in motility and respiration, whilst in M. tuberuclosis, functions included intracellular growth and regulation of virulence. Using publicly available short-read data we predicted previously unrecognized, large amplifications in B. pertussis and M. tuberculosis. This reveals the unrecognized and dynamic genetic diversity of B. pertussis and M. tuberculosis, highlighting the need for a more holistic understanding of bacterial genetics.

Delayed lactose utilization among Shiga toxin-producing Escherichia coli of serogroup O121

Two outbreaks of Shiga toxin-producing Escherichia coli O121:H19 associated with wheat flour, in the United States of America and Canada, involved strains with an unusual phenotype, delayed lactose utilization (DLU). These strains do not ferment lactose when initially cultured on MacConkey agar (MAC), but lactose fermentation occurs following subculture to a second plate of MAC. The prevalence of DLU was determined by examining the β-galactosidase activity of 49 strains of E. coli O121, and of 37 other strains of E. coli. Twenty four of forty three O121:H19 and one O121:NM displayed DLU. Two strains (O121:NM and O145:H34) did not have detectable β-galactosidase activity. β-glucuronidase activity of O121 strains was also determined. All but six DLU strains had normal β-glucuronidase activity. β-glucuronidase activity was suppressed on MAC for 17 of 23 O121 non-DLU strains. Genomic analysis found that DLU strains possessed an insertion sequence, IS600 (1267 bp), between lacZ (β-galactosidase) and lacY (β-galactoside permease), that was not present in strains exhibiting normal lactose utilization. The insert might reduce the expression of β-galactoside permease, delaying import of lactose, resulting in the DLU phenotype. The high probability of DLU should be considered when using lactose-containing media for the isolation of STEC O121.

Maturation of 23S rRNA includes removal of helix H1 in many bacteria

In most bacteria, the three ribosomal RNAs (rRNAs) are encoded together in each of several near-identical operons. As soon as the nascent precursor rRNA emerges from RNA polymerase, ribosome assembly begins. This process entails ribosomal protein binding, rRNA folding, rRNA modification, and rRNA processing. In the model organisms Escherichia coli and Bacillus subtilis, rRNA processing results in similar mature rRNAs, despite substantial differences in the cohort of RNAses involved. A recent study of Flavobacterium johnsoniae, a member of the phylum Bacteroidota (formerly Bacteroidetes), revealed that helix H1 of 23S rRNA is absent from ribosomes, apparently a consequence of rRNA maturation. In this work, we mined RNA-seq data from 19 individual organisms and ocean metatranscriptomic samples to compare rRNA processing across diverse bacterial lineages. We found that mature ribosomes from multiple clades lack H1, and typically these ribosomes also lack an encoded H98. For all groups analysed, H1 is predicted to form in precursor rRNA as part of a longer leader-trailer helix. Hence, we infer that evolutionary loss of H98 sets the stage for H1 removal during 50S subunit maturation.

Identification of BvgA-Dependent and BvgA-Independent Small RNAs (sRNAs) in Bordetella pertussis Using the Prokaryotic sRNA Prediction Toolkit ANNOgesic

Noncoding small RNAs (sRNAs) are crucial for the posttranscriptional regulation of gene expression in all organisms and are known to be involved in the regulation of bacterial virulence. In the human pathogen Bordetella pertussis, which causes whooping cough, virulence is controlled primarily by the master two-component system BvgA (response regulator)/BvgS (sensor kinase). In this system, BvgA is phosphorylated (Bvg⁺ mode) or nonphosphorylated (Bvg^- mode), with global transcriptional differences between the two. B. pertussis also carries the bacterial sRNA chaperone Hfq, which has previously been shown to be required for virulence. Here, we conducted transcriptomic analyses to identify possible B. pertussis sRNAs and to determine their BvgAS dependence using transcriptome sequencing (RNA-seq) and the prokaryotic sRNA prediction program ANNOgesic. We identified 143 possible candidates (25 Bvg⁺ mode specific and 53 Bvg^- mode specific), of which 90 were previously unreported. Northern blot analyses confirmed all of the 10 ANNOgesic candidates that we tested. Homology searches demonstrated that 9 of the confirmed sRNAs are highly conserved among B. pertussis, Bordetella parapertussis, and Bordetella bronchiseptica, with one that also has homologues in other species of the Alcaligenaceae family. Using coimmunoprecipitation with a B. pertussis FLAG-tagged Hfq, we demonstrated that 3 of the sRNAs interact directly with Hfq, which is the first identification of sRNA binding to B. pertussis Hfq. Our study demonstrates that ANNOgesic is a highly useful tool for the identification of sRNAs in this system and that its combination with molecular techniques is a successful way to identify various BvgAS-dependent and Hfq-binding sRNAs. IMPORTANCE Noncoding small RNAs (sRNAs) are crucial for posttranscriptional regulation of gene expression in all organisms and are known to be involved in the regulation of bacterial virulence. We have investigated the presence of sRNAs in the obligate human pathogen B. pertussis, using transcriptome sequencing (RNA-seq) and the recently developed prokaryotic sRNA search program ANNOgesic. This analysis has identified 143 sRNA candidates (90 previously unreported). We have classified their dependence on the B. pertussis two-component system required for virulence, namely, BvgAS, based on their expression in the presence/absence of the phosphorylated response regulator BvgA, confirmed several by Northern analyses, and demonstrated that 3 bind directly to B. pertussis Hfq, the RNA chaperone involved in mediating sRNA effects. Our study demonstrates the utility of combining RNA-seq, ANNOgesic, and molecular techniques to identify various BvgAS-dependent and Hfq-binding sRNAs, which may unveil the roles of sRNAs in pertussis pathogenesis.

A Unique Reverse Adaptation Mechanism Assists Bordetella pertussis in Resistance to Both Scarcity and Toxicity of Manganese

The ability of bacterial pathogens to acquire essential micronutrients is critical for their survival in the host environment. Manganese plays a complex role in the virulence of a variety of pathogens due to its function as an antioxidant and enzymatic cofactor. Therefore, host cells deprive pathogens of manganese to prevent or attenuate infection. Here, we show that evolution of the human-restricted pathogen Bordetella pertussis has selected for an inhibitory duplication within a manganese exporter of the calcium:cation antiporter superfamily. Intriguingly, upon exposure to toxic levels of manganese, the nonfunctional exporter becomes operative in resister cells due to a unique reverse adaptation mechanism. However, compared with wild-type (wt) cells, the resisters carrying a functional copy of the exporter displayed strongly reduced intracellular levels of manganese and impaired growth under oxidative stress. Apparently, inactivation of the manganese exporter and the resulting accumulation of manganese in the cytosol benefited the pathogen by improving its survival under stress conditions. The inhibitory duplication within the exporter gene is highly conserved among B. pertussis strains, absent from all other Bordetella species and from a vast majority of organisms across all kingdoms of life. Therefore, we conclude that inactivation of the exporter gene represents an exceptional example of a flexible genome decay strategy employed by a human pathogen to adapt to its exclusive host.

RNase III and RNase E Influence Posttranscriptional Regulatory Networks Involved in Virulence Factor Production, Metabolism, and Regulatory RNA Processing in Bordetella pertussis

Bordetella pertussis has been shown to encode regulatory RNAs, yet the posttranscriptional regulatory circuits on which they act remain to be fully elucidated. We generated mutants lacking the endonucleases RNase III and RNase E and assessed their individual impact on the B. pertussis transcriptome. Transcriptome sequencing (RNA-Seq) analysis showed differential expression of ∼25% of the B. pertussis transcriptome in each mutant, with only 28% overlap between data sets. Both endonucleases exhibited substantial impact on genes involved in amino acid uptake (e.g., ABC transporters) and in virulence (e.g., the type III secretion system and the autotransporters vag8, tcfA, and brkA). Interestingly, mutations in RNase III and RNase E drove the stability of many transcripts, including those involved in virulence, in opposite directions, a result that was validated by qPCR and immunoblotting for tcfA and brkA. Of note, whereas similar mutations to RNase E in Escherichia coli have subtle effects on transcript stability, a striking >20-fold reduction in four gene transcripts, including tcfA and vag8, was observed in B. pertussis. We further compared our data set to the regulon controlled by the RNA chaperone Hfq to identify B. pertussis loci influenced by regulatory RNAs. This analysis identified ∼120 genes and 19 operons potentially regulated at the posttranscriptional level. Thus, our findings revealed how changes in RNase III- and RNase E-mediated RNA turnover influence pathways associated with virulence and cellular homeostasis. Moreover, we highlighted loci potentially influenced by regulatory RNAs, providing insights into the posttranscriptional regulatory networks involved in fine-tuning B. pertussis gene expression. IMPORTANCE Noncoding, regulatory RNAs in bacterial pathogens are critical components required for rapid changes in gene expression profiles. However, little is known about the role of regulatory RNAs in the growth and pathogenesis of Bordetella pertussis. To address this, mutants separately lacking ribonucleases central to regulatory RNA processing, RNase III and RNase E, were analyzed by RNA-Seq. Here, we detail the first transcriptomic analysis of the impact of altered RNA degradation in B. pertussis. Each mutant showed approximately 1,000 differentially expressed genes, with significant changes in the expression of pathways associated with metabolism, bacterial secretion, and virulence factor production. Our analysis suggests an important role for these ribonucleases during host colonization and provides insights into the breadth of posttranscriptional regulation in B. pertussis, further informing our understanding of B. pertussis pathogenesis.

Comparative Omics Analysis of Historic and Recent Isolates of Bordetella pertussis and Effects of Genome Rearrangements on Evolution

Despite high vaccination coverage, pertussis is increasing in many industrialized countries, including the Czech Republic. To better understand Bordetella pertussis resurgence, we analyzed historic strains and recent clinical isolates by using a comparative omics approach. Whole-genome sequencing showed that historic and recent isolates of B. pertussis have substantial variation in genome organization and form separate phylogenetic clusters. Subsequent RNA sequence analysis and liquid chromatography with mass tandem spectrometry analyses showed that these variations translated into discretely separated transcriptomic and proteomic profiles. When compared with historic strains, recent isolates showed increased expression of flagellar genes and genes involved in lipopolysaccharide biosynthesis and decreased expression of polysaccharide capsule genes. Compared with reference strain Tohama I, all strains had increased expression and production of the type III secretion system apparatus. We detected the potential link between observed effects and insertion sequence element–induced changes in gene context only for a few genes.

Radical genome remodelling accompanied the emergence of a novel host-restricted bacterial pathogen

The emergence of new pathogens is a major threat to public and veterinary health. Changes in bacterial habitat such as a switch in host or disease tropism are typically accompanied by genetic diversification. Staphylococcus aureus is a multi-host bacterial species associated with human and livestock infections. A microaerophilic subspecies, Staphylococcus aureus subsp. anaerobius, is responsible for Morel’s disease, a lymphadenitis restricted to sheep and goats. However, the evolutionary history of S. aureus subsp. anaerobius and its relatedness to S. aureus are unknown. Population genomic analyses of clinical S. aureus subsp. anaerobius isolates revealed a highly conserved clone that descended from a S. aureus progenitor about 1000 years ago before differentiating into distinct lineages that contain African and European isolates. S. aureus subsp. anaerobius has undergone limited clonal expansion, with a restricted population size, and an evolutionary rate 10-fold slower than S. aureus. The transition to its current restricted ecological niche involved acquisition of a pathogenicity island encoding a ruminant host-specific effector of abscess formation, large chromosomal re-arrangements, and the accumulation of at least 205 pseudogenes, resulting in a highly fastidious metabolism. Importantly, expansion of ~87 insertion sequences (IS) located largely in intergenic regions provided distinct mechanisms for the control of expression of flanking genes, including a novel mechanism associated with IS-mediated anti-anti-sense decoupling of ancestral gene repression. Our findings reveal the remarkable evolutionary trajectory of a host-restricted bacterial pathogen that resulted from extensive remodelling of the S. aureus genome through an array of diverse mechanisms in parallel.

The emergence of new pathogens is a major threat to public and veterinary health. Some bacteria such as Staphylococcus aureus, have the capacity to infect many different host species including humans and livestock while others such as the closely-related S. aureus subsp. anaerobius, associated with a single type of pathology called Morel’s disease in small ruminants, are highly niche-restricted. However, our understanding of the genetic basis for such differences in bacterial host-tropism is very limited. Here, we discovered that S. aureus subsp. anaerobius evolved from an S. aureus ancestor and underwent an array of extensive changes to its genome that accompanied the transition to its current restricted lifestyle. We observed genome decay involving loss of function of hundreds of genes, large intra-chromosomal rearrangements affecting most of the genome, acquisition of a pathogenicity island, and expansion of large numbers of insertion sequences that are inserted at intergenic sites around the genome. Importantly, we found that IS elements affect the expression of neighbouring genes in different ways including a novel mechanism of IS-enabled disruption of ancestral gene repression. Taken together, we provide a remarkable example of radical genomic changes associated with evolutionary transition from a multi-host to highly restricted host ecology.

Omics Analysis of Blood-Responsive Regulon in Bordetella pertussis Identifies a Novel Essential T3SS Substrate

Bacterial pathogens sense specific cues associated with different host niches and integrate these signals to appropriately adjust the global gene expression. Bordetella pertussis is a Gram-negative, strictly human pathogen of the respiratory tract and the etiological agent of whooping cough (pertussis). Though B. pertussis does not cause invasive infections, previous results indicated that this reemerging pathogen responds to blood exposure. Here, omics RNA-seq and LC–MS/MS techniques were applied to determine the blood-responsive regulon of B. pertussis. These analyses revealed that direct contact with blood rewired global gene expression profiles in B. pertussis as the expression of almost 20% of all genes was significantly modulated. However, upon loss of contact with blood, the majority of blood-specific effects vanished, with the exception of several genes encoding the T3SS-secreted substrates. For the first time, the T3SS regulator BtrA was identified in culture supernatants of B. pertussis. Furthermore, proteomic analysis identified BP2259 protein as a novel secreted T3SS substrate, which is required for T3SS functionality. Collectively, presented data indicate that contact with blood represents an important cue for B. pertussis cells.

A Mutation Upstream of the rplN-rpsD Ribosomal Operon Downregulates Bordetella pertussis Virulence Factor Production without Compromising Bacterial Survival within Human Macrophages

The BvgS/BvgA two-component system controls expression of ∼550 genes of Bordetella pertussis, of which, ∼245 virulence-related genes are positively regulated by the BvgS-phosphorylated transcriptional regulator protein BvgA (BvgA∼P). We found that a single G-to-T nucleotide transversion in the 5'-untranslated region (5'-UTR) of the rplN gene enhanced transcription of the ribosomal protein operon and of the rpoA gene and provoked global dysregulation of B. pertussis genome expression. This comprised overproduction of the alpha subunit (RpoA) of the DNA-dependent RNA polymerase, downregulated BvgA and BvgS protein production, and impaired production and secretion of virulence factors by the mutant. Nonetheless, the mutant survived like the parental bacteria for >2 weeks inside infected primary human macrophages and persisted within infected mouse lungs for a longer period than wild-type B. pertussis These observations suggest that downregulation of virulence factor production by bacteria internalized into host cells may enable persistence of the whooping cough agent in the airways.IMPORTANCE We show that a spontaneous mutation that upregulates transcription of an operon encoding ribosomal proteins and causes overproduction of the downstream-encoded α subunit (RpoA) of RNA polymerase causes global effects on gene expression levels and proteome composition of Bordetella pertussis Nevertheless, the resulting important downregulation of the BvgAS-controlled expression of virulence factors of the whooping cough agent did not compromise its capacity to persist for prolonged periods inside primary human macrophage cells, and it even enhanced its capacity to persist in infected mouse lungs. These observations suggest that the modulation of BvgAS-controlled expression of virulence factors may occur also during natural infections of human airways by Bordetella pertussis and may possibly account for long-term persistence of the pathogen within infected cells of the airways.

Impact of insertion sequences on convergent evolution of Shigella species

Shigella species are specialised lineages of Escherichia coli that have converged to become human-adapted and cause dysentery by invading human gut epithelial cells. Most studies of Shigella evolution have been restricted to comparisons of single representatives of each species; and population genomic studies of individual Shigella species have focused on genomic variation caused by single nucleotide variants and ignored the contribution of insertion sequences (IS) which are highly prevalent in Shigella genomes. Here, we investigate the distribution and evolutionary dynamics of IS within populations of Shigella dysenteriae Sd1, Shigella sonnei and Shigella flexneri. We find that five IS (IS1, IS2, IS4, IS600 and IS911) have undergone expansion in all Shigella species, creating substantial strain-to-strain variation within each population and contributing to convergent patterns of functional gene loss within and between species. We find that IS expansion and genome degradation are most advanced in S. dysenteriae and least advanced in S. sonnei; and using genome-scale models of metabolism we show that Shigella species display convergent loss of core E. coli metabolic capabilities, with S. sonnei and S. flexneri following a similar trajectory of metabolic streamlining to that of S. dysenteriae. This study highlights the importance of IS to the evolution of Shigella and provides a framework for the investigation of IS dynamics and metabolic reduction in other bacterial species.

Combined RNAseq and ChIPseq Analyses of the BvgA Virulence Regulator of Bordetella pertussis

Bordetella pertussis, the etiological agent of whooping cough, remains a major global health problem. Despite the global usage of whole-cell vaccines since the 1950s and of acellular vaccines in the 1990s, it still is one of the most prevalent vaccine-preventable diseases in industrialized countries. Virulence of B. pertussis is controlled by BvgA/S, a two-component system responsible for upregulation of virulence-activated genes (vags) and downregulation of virulence-repressed genes (vrgs). By transcriptome sequencing (RNAseq) analyses, we identified more than 270 vags or vrgs, and chromatin immunoprecipitation sequencing (ChIPseq) analyses revealed 148 BvgA-binding sites, 91 within putative promoter regions, 52 within open reading frames, and 5 in noncoding regions. Some vags, such as dnt and fhaL, do not contain a BvgA-binding site, suggesting indirect regulation. In contrast, several vrgs and some genes not identified by RNAseq analyses under laboratory conditions contain strong BvgA-binding sites, indicating previously unappreciated complexities of BvgA/S biology.

Bordetella pertussis regulates the production of its virulence factors by the two-component system BvgAS. In the virulence phase, BvgS phosphorylates BvgA, which then activates the transcription of virulence-activated genes (vags). In the avirulence phase, such as during growth in the presence of MgSO₄, BvgA is not phosphorylated and the vags are not expressed. Instead, a set of virulence-repressed genes (vrgs) is expressed. Here, we performed transcriptome sequencing (RNAseq) analyses on B. pertussis cultivated with or without MgSO₄ and on a BvgA-deficient Tohama I derivative. We observed that 146 genes were less expressed under modulating conditions or in the BvgA-deficient strain than under the nonmodulating condition, while 130 genes were more expressed. Some of the genes code for proteins with regulatory functions, suggesting a BvgA/S regulation cascade. To determine which genes are directly regulated by BvgA, we performed chromatin immunoprecipitation sequencing (ChIPseq) analyses. We identified 148 BvgA-binding sites, 91 within putative promoter regions, 52 within open reading frames, and 5 in noncoding regions. Among the former, 32 are in BvgA-regulated putative promoter regions. Some vags, such as dnt and fhaL, contain no BvgA-binding site, suggesting indirect BvgA regulation. Unexpectedly, BvgA also bound to some vrg putative promoter regions. Together, these observations indicate an unrecognized complexity of BvgA/S biology.

IMPORTANCEBordetella pertussis, the etiological agent of whooping cough, remains a major global health problem. Despite the global usage of whole-cell vaccines since the 1950s and of acellular vaccines in the 1990s, it still is one of the most prevalent vaccine-preventable diseases in industrialized countries. Virulence of B. pertussis is controlled by BvgA/S, a two-component system responsible for upregulation of virulence-activated genes (vags) and downregulation of virulence-repressed genes (vrgs). By transcriptome sequencing (RNAseq) analyses, we identified more than 270 vags or vrgs, and chromatin immunoprecipitation sequencing (ChIPseq) analyses revealed 148 BvgA-binding sites, 91 within putative promoter regions, 52 within open reading frames, and 5 in noncoding regions. Some vags, such as dnt and fhaL, do not contain a BvgA-binding site, suggesting indirect regulation. In contrast, several vrgs and some genes not identified by RNAseq analyses under laboratory conditions contain strong BvgA-binding sites, indicating previously unappreciated complexities of BvgA/S biology.

Expression of small RNAs of Bordetella pertussis colonizing murine tracheas

We performed RNA sequencing on Bordetella pertussis, the causative agent of whooping cough, and identified nine novel small RNAs (sRNAs) that were transcribed during the bacterial colonization of murine tracheas. Among them, four sRNAs were more strongly expressed in vivo than in vitro. Moreover, the expression of eight sRNAs was not regulated by the BvgAS two‐component system, which is the master regulator for the expression of genes contributing to the bacterial infection. The present results suggest a BvgAS‐independent gene regulatory system involving the sRNAs that is active during B. pertussis infection.

Comparative genomics of whole-cell pertussis vaccine strains from India

BACKGROUND

Despite high vaccination coverage using acellular (ACV) and whole-cell pertussis (WCV) vaccines, the resurgence of pertussis is observed globally. Genetic divergence in circulating strains of Bordetella pertussis has been reported as one of the contributing factors for the resurgence of the disease. Our current knowledge of B. pertussis genetic evolution in circulating strains is mostly based on studies conducted in countries using ACVs targeting only a few antigens used in the production of ACVs. To better understand the adaptation to vaccine-induced selection pressure, it will be essential to study B. pertussis populations in developing countries which are using WCVs. India is a significant user and global supplier of WCVs. We report here comparative genome analyses of vaccine and clinical isolates reported from India. Whole-genome sequences obtained from vaccine strains: WCV (J445, J446, J447 and J448), ACV (BP165) were compared with Tohama-I reference strain and recently reported clinical isolates from India (BPD1, BPD2). Core genome-based phylogenetic analysis was also performed using 166 isolates reported from countries using ACV.

RESULTS

Whole-genome analysis of vaccine and clinical isolates reported from India revealed high genetic similarity and conserved genome among strains. Phylogenetic analysis showed that clinical and vaccine strains share genetic closeness with reference strain Tohama-I. The allelic profile of vaccine strains (J445:ptxP1/ptxA2/prn1/fim2-1/fim3-1; J446: ptxP2/ptxA4/prn7/fim2-2/fim3-1; J447 and J448: ptxP1/ptxA1/ prn1/fim2-1/fim3-1), which matched entirely with clinical isolates (BPD1:ptxP1/ptxA1/prn1/fim2-1 and BPD2: ptxP1/ptxA1/prn1/fim2-1) reported from India. Multi-locus sequence typing (MLST) demonstrated the presence of dominant sequence types ST2 and primitive ST1 in vaccine strains which will allow better coverage against circulating strains of B. pertussis.

CONCLUSIONS

The study provides a detailed characterization of vaccine and clinical strains reported from India, which will further facilitate epidemiological studies on genetic shifts in countries which are using WCVs in their immunization programs.

Identification of RNA 3´ ends and termination sites in Haloferax volcanii

Archaeal genomes are densely packed; thus, correct transcription termination is an important factor for orchestrated gene expression. A systematic analysis of RNA 3´ termini, to identify transcription termination sites (TTS) using RNAseq data has hitherto only been performed in two archaea, Methanosarcina mazei and Sulfolobus acidocaldarius. In this study, only regions directly downstream of annotated genes were analysed, and thus, only part of the genome had been investigated. Here, we developed a novel algorithm (Internal Enrichment-Peak Calling) that allows an unbiased, genome-wide identification of RNA 3´ termini independent of annotation. In an RNA fraction enriched for primary transcripts by terminator exonuclease (TEX) treatment we identified 1,543 RNA 3´ termini. Approximately half of these were located in intergenic regions, and the remainder were found in coding regions. A strong sequence signature consistent with known termination events at intergenic loci indicates a clear enrichment for native TTS among them. Using these data we determined distinct putative termination motifs for intergenic (a T stretch) and coding regions (AGATC). In vivo reporter gene tests of selected TTS confirmed termination at these sites, which exemplify the different motifs. For several genes, more than one termination site was detected, resulting in transcripts with different lengths of the 3´ untranslated region (3´ UTR).

The world of asRNAs in Gram-negative and Gram-positive bacteria

Bacteria exhibit an amazing diversity of mechanisms controlling gene expression to both maintain essential functions and modulate accessory functions in response to environmental cues. Over the years, it has become clear that bacterial regulation of gene expression is still far from fully understood. This review focuses on antisense RNAs (asRNAs), a class of RNA regulators defined by their location in cis and their perfect complementarity with their targets, as opposed to small RNAs (sRNAs) which act in trans with only short regions of complementarity. For a long time, only few functional asRNAs in bacteria were known and were almost exclusively found on mobile genetic elements (MGEs), thus, their importance among the other regulators was underestimated. However, the extensive application of transcriptomic approaches has revealed the ubiquity of asRNAs in bacteria. This review aims to present the landscape of studied asRNAs in bacteria by comparing 67 characterized asRNAs from both Gram-positive and Gram-negative bacteria. First we describe the inherent ambiguity in the existence of asRNAs in bacteria, second, we highlight their diversity and their involvement in all aspects of bacterial life. Finally we compare their location and potential mode of action toward their target between Gram-negative and Gram-positive bacteria and present tendencies and exceptions that could lead to a better understanding of asRNA functions.

iRAPs curb antisense transcription in E. coli

RNA polymerase-binding RNA aptamers (RAPs) are natural RNA elements that control transcription in cis by directly contacting RNA polymerase. Many RAPs inhibit transcription by inducing Rho-dependent termination in Escherichia coli. Here, we studied the role of inhibitory RAPs (iRAPs) in modulation of antisense transcription (AT) using in silico and in vivo approaches. We revisited the antisense transcriptome in cells with impaired AT regulators (Rho, H-NS and RNaseIII) and searched for the presence of RAPs within antisense RNAs. Many of these RAPs were found at key genomic positions where they terminate AT. By exploring the activity of several RAPs both in a reporter system and in their natural genomic context, we confirmed their significant role in AT regulation. RAPs coordinate Rho activity at the antisense strand and terminate antisense transcripts. In some cases, they stimulated sense expression by alleviating ongoing transcriptional interference. Essentially, our data postulate RAPs as key determinants of Rho-mediated AT regulation in E. coli.

A high resolution A-to-I editing map in the mouse identifies editing events controlled by pre-mRNA splicing

Pre-mRNA-splicing and adenosine to inosine (A-to-I) RNA-editing occur mostly cotranscriptionally. During A-to-I editing, a genomically encoded adenosine is deaminated to inosine by adenosine deaminases acting on RNA (ADARs). Editing-competent stems are frequently formed between exons and introns. Consistently, studies using reporter assays have shown that splicing efficiency can affect editing levels. Here, we use Nascent-seq and identify ∼90,000 novel A-to-I editing events in the mouse brain transcriptome. Most novel sites are located in intronic regions. Unlike previously assumed, we show that both ADAR (ADAR1) and ADARB1 (ADAR2) can edit repeat elements and regular transcripts to the same extent. We find that inhibition of splicing primarily increases editing levels at hundreds of sites, suggesting that reduced splicing efficiency extends the exposure of intronic and exonic sequences to ADAR enzymes. Lack of splicing factors NOVA1 or NOVA2 changes global editing levels, demonstrating that alternative splicing factors can modulate RNA editing. Finally, we show that intron retention rates correlate with editing levels across different brain tissues. We therefore demonstrate that splicing efficiency is a major factor controlling tissue-specific differences in editing levels.

Comparative Integrated Omics Analysis of the Hfq Regulon in Bordetella pertussis

Bordetella pertussis is a Gram-negative strictly human pathogen of the respiratory tract and the etiological agent of whooping cough (pertussis). Previously, we have shown that RNA chaperone Hfq is required for virulence of B. pertussis. Furthermore, microarray analysis revealed that a large number of genes are affected by the lack of Hfq. This study represents the first attempt to characterize the Hfq regulon in bacterial pathogen using an integrative omics approach. Gene expression profiles were analyzed by RNA-seq and protein amounts in cell-associated and cell-free fractions were determined by LC-MS/MS technique. Comparative analysis of transcriptomic and proteomic data revealed solid correlation (r² = 0.4) considering the role of Hfq in post-transcriptional control of gene expression. Importantly, our study confirms and further enlightens the role of Hfq in pathogenicity of B. pertussis as it shows that Δhfq strain displays strongly impaired secretion of substrates of Type III secretion system (T3SS) and substantially reduced resistance to serum killing. On the other hand, significantly increased production of proteins implicated in transport of important metabolites and essential nutrients observed in the mutant seems to compensate for the physiological defect introduced by the deletion of the hfq gene.

Physiological roles of antisense RNAs in prokaryotes

Prokaryotes encounter constant and often brutal modifications to their environment. In order to survive, they need to maintain fitness, which includes adapting their protein expression patterns. Many factors control gene expression but this review focuses on just one, namely antisense RNAs (asRNAs), a class of non-coding RNAs (ncRNAs) characterized by their location in cis and their perfect complementarity with their targets. asRNAs were considered for a long time to be trivial and only to be found on mobile genetic elements. However, recent advances in methodology have revealed that their abundance and potential activities have been underestimated. This review aims to illustrate the role of asRNA in various physiologically crucial functions in both archaea and bacteria, which can be regrouped in three categories: cell maintenance, horizontal gene transfer and virulence. A literature survey of asRNAs demonstrates the difficulties to characterize and assign a role to asRNAs. With the aim of facilitating this task, we describe recent technological advances that could be of interest to identify new asRNAs and to discover their function.

How Genomics Is Changing What We Know About the Evolution and Genome of Bordetella pertussis

The evolution of Bordetella pertussis from a common ancestor similar to Bordetella bronchiseptica has occurred through large-scale gene loss, inactivation and rearrangements, largely driven by the spread of insertion sequence element repeats throughout the genome. B. pertussis is widely considered to be monomorphic, and recent evolution of the B. pertussis genome appears to, at least in part, be driven by vaccine-based selection. Given the recent global resurgence of whooping cough despite the wide-spread use of vaccination, a more thorough understanding of B. pertussis genomics could be highly informative. In this chapter we discuss the evolution of B. pertussis, including how vaccination is changing the circulating B. pertussis population at the gene-level, and how new sequencing technologies are revealing previously unknown levels of inter- and intra-strain variation at the genome-level.

Signal transduction-dependent small regulatory RNA is involved in glutamate metabolism of the human pathogen Bordetella pertussis

Bordetella pertussis is the causative agent of human whooping cough, a highly contagious respiratory disease which despite vaccination programs remains the major cause of infant morbidity and mortality. The requirement of the RNA chaperone Hfq for virulence of B. pertussis suggested that Hfq-dependent small regulatory RNAs are involved in the modulation of gene expression. High-throughput RNA sequencing revealed hundreds of putative noncoding RNAs including the RgtA sRNA. Abundance of RgtA is strongly decreased in the absence of the Hfq protein and its expression is modulated by the activities of the two-component regulatory system BvgAS and another response regulator RisA. Whereas RgtA levels were elevated under modulatory conditions or in the absence of bvg genes, deletion of the risA gene completely abolished RgtA expression. Profiling of the ΔrgtA mutant in the ΔbvgA genetic background identified the BP3831 gene encoding a periplasmic amino acid-binding protein of an ABC transporter as a possible target gene. The results of site-directed mutagenesis and in silico analysis indicate that RgtA base-pairs with the region upstream of the start codon of the BP3831 mRNA and thereby weakens the BP3831 protein production. Furthermore, our data suggest that the function of the BP3831 protein is related to transport of glutamate, an important metabolite in the B. pertussis physiology. We propose that the BvgAS/RisA interplay regulates the expression of RgtA which upon infection, when glutamate might be scarce, attenuates translation of the glutamate transporter and thereby assists in adaptation of the pathogen to other sources of energy.

A Bordetella pertussis MgtC homolog plays a role in the intracellular survival

Bordetella pertussis, the causative agent of whooping cough, has the capability to survive inside the host cells. This process requires efficient adaptation of the pathogen to the intracellular environment and the associated stress. Among the proteins produced by the intracellular B. pertussis we identified a protein (BP0414) that shares homology with MgtC, a protein which was previously shown to be involved in the intracellular survival of other pathogens. To explore if BP0414 plays a role in B. pertussis intracellular survival a mutant strain defective in the production of this protein was constructed. Using standard in vitro growth conditions we found that BP0414 is required for B. pertussis growth under low magnesium availability or low pH, two environmental conditions that this pathogen might face within the host cell. Intracellular survival studies showed that MgtC is indeed involved in B. pertussis viability inside the macrophages. The use of bafilomycin A1, which inhibits phagosome acidification, abolished the survival defect of the mgtC deficient mutant strain suggesting that in intracellular B. pertussis the role of MgtC protein is mainly related to the bacterial adaptation to the acidic conditions found inside the of phagosomes. Overall, this work provides an insight into the importance of MgtC in B. pertussis pathogenesis and its contribution to bacterial survival within immune cells.