<figure><img src="image_1.png" /></figure>

HAL
open science

# Genome-Wide Transcription Start Site Mapping and Promoter Assignments to a Sigma Factor in the Human Enteropathogen Clostridioides difficile

Olga Soutourina, Thomas Dubois, Marc Monot, Pavel Shelyakin, Laure

Saujet, Pierre Boudry, Mikhail Gelfand, Bruno Dupuy, Isabelle

Martin-Verstraete

► To cite this version:

Olga Soutourina, Thomas Dubois, Marc Monot, Pavel Shelyakin, Laure Saujet, et al. Genome-Wide Transcription Start Site Mapping and Promoter Assignments to a Sigma Factor in the Human Enteropathogen Clostridioides difficile. Frontiers in Microbiology, 2020, 11, pp.1939. 10.3389/fmicb.2020.01939. [hal-02921326](https://hal.inrae.fr/hal-02921326)

HAL Id: hal-02921326

https://hal.inrae.fr/hal-02921326

Submitted on 15 Oct 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

OPEN ACCESS

**Edited by:**
Miroslav Patek,
Academy of Sciences of the Czech
Republic, Czechia

**Reviewed by:**
Björn Voß,
University of Stuttgart, Germany
Masaya Fujita,
University of Houston, United States

*Correspondence:
Olga Soutourina
olga.soutourina@i2bc.paris-saclay.fr
Isabelle Martin-Verstraete
isabelle.martin-verstraete@pasteur.fr

†These authors have contributed
equally to this work

‡Present address:
Thomas Dubois,
UMR UMET, INRA, CNRS, Univ. Lille
1, Villeneuve d'Ascq, France
Laure Saujet,
Pherecydes Pharma, Romainville,
France
Pierre Boudry,
Université Paris-Saclay, CEA, CNRS,
Institute for Integrative Biology of the
Cell (I2BC), Gif-sur-Yvette, France

**Specialty section:**
This article was submitted to
Microbial Physiology and Metabolism,
a section of the journal
Frontiers in Microbiology

Received: 28 April 2020

Accepted: 23 July 2020

Published: 13 August 2020

**Citation:**
Soutourina O, Dubois T, Monot M, Shelyakin PV, Saujet L, Boudry P, Gelfand MS, Dupuy B and Martin-Verstraete I (2020) Genome-Wide Transcription Start Site Mapping and Promoter Assignments to a Sigma Factor in the Human Enteropathogen Clostridioides difficile. *Front. Microbiol.* 11:1939. doi: 10.3389/fmicb.2020.01939

# Genome-Wide Transcription Start Site Mapping and Promoter Assignments to a Sigma Factor in the Human Enteropathogen Clostridioides difficile

Olga Soutourina<fnref n="1" />,<fnref n="2" />,³*†, Thomas Dubois<fnref n="1" />†*, Marc Monot<fnref n="1" />†, Pavel V. Shelyakin⁴,
Laure Saujet<fnref n="1" />‡, Pierre Boudry<fnref n="1" />‡, Mikhail S. Gelfand⁴,⁵, Bruno Dupuy<fnref n="1" /> and
Isabelle Martin-Verstraete<fnref n="1" />,<fnref n="2" />*

The emerging human enteropathogen *Clostridioides difficile* is the main cause of diarrhea associated with antibiotherapy. Regulatory pathways underlying the adaptive responses remain understudied and the global view of *C. difficile* promoter structure is still missing. In the genome of *C. difficile* 630, 22 genes encoding sigma factors are present suggesting a complex pattern of transcription in this bacterium. We present here the first transcriptional map of the *C. difficile* genome resulting from the identification of transcriptional start sites (TSS), promoter motifs and operon structures. By 5'-end RNA-seq approach, we mapped more than 1000 TSS upstream of genes. In addition to these primary TSS, this analysis revealed complex structure of transcriptional units such as alternative and internal promoters, potential RNA processing events and 5' untranslated regions. By following an *in silico* iterative strategy that used as an input previously published consensus sequences and transcriptomic analysis, we identified candidate promoters upstream of most of protein-coding and non-coding RNAs genes. This strategy also led to refine consensus sequences of promoters recognized by major sigma factors of *C. difficile*. Detailed analysis focuses on the transcription in the pathogenicity locus and regulatory genes, as well as regulons of transition phase and sporulation sigma factors as important components of *C. difficile* regulatory network governing toxin gene expression and spore formation. Among the still uncharacterized regulons of the major sigma factors of *C. difficile*, we defined the SigL regulon by combining transcriptome and *in silico* analyses. We showed that the SigL regulon is largely involved in amino-acid degradation, a metabolism crucial for *C. difficile* gut colonization. Finally, we combined our TSS mapping, *in silico* identification of promoters and RNA-seq data to improve gene annotation and to suggest operon organization

<footnote n="1">¹ Laboratoire Pathogenèses des Bactéries Anaérobies, Institut Pasteur, UMR CNRS 2001, Université de Paris, Paris, France,</footnote>

<footnote n="2">² Institut Universitaire de France, Paris, France, ³ Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the</footnote>

<footnote>Cell (I2BC), Gif-sur-Yvette, France, ⁴ Institute for Information Transmission Problems, Moscow, Russia, ⁵ Skolkovo Institute</footnote>

<footnote>of Science and Technology, Moscow, Russia</footnote>
in *C. difficile*. These data will considerably improve our knowledge of global regulatory circuits controlling gene expression in *C. difficile* and will serve as a useful rich resource for scientific community both for the detailed analysis of specific genes and systems biology studies.

**Keywords:** transcription initiation, transcription unit architecture, sigma factors, sigma 54, amino acid catabolism

INTRODUCTION

*Clostridioides difficile* (formerly *Clostridium difficile*) is an emerging human enteropathogen causing nosocomial antibiotic-associated diarrhea in adults (Carroll and Bartlett, 2011). This pathogen became a key public health issue worldwide. Alarming incidence of *C. difficile* infections was further accentuated by the recent emergence of antibiotic resistance, of hypervirulent epidemic strains broadening the population at risk and severity of disease, high rate of recurrent infection as well as an overall aging of population in industrial countries (Rupnik et al., 2009; Banawas, 2018). This anaerobic, spore-forming, Gram-positive bacterium can be found in soil and aquatic environments as well as in intestinal tracts of humans and animals (Keessen et al., 2011). In humans, the *C. difficile* carriage may be asymptomatic, but when symptoms appear, they can vary from mild diarrhea to life-threatening pseudomembranous colitis. The main *C. difficile* virulence factors are two toxins, TcdA and TcdB, produced by toxigenic strains (Vedantam et al., 2012). A binary toxin CDT is also present in some isolates, as well as additional factors that help colonization, like adhesins, pili, and flagella (Janoir, 2016). Despite recent efforts of research community, a lot of questions remains unanswered on the regulatory processes responsible for the adaptation of *C. difficile* inside the host. This information is largely missing and urgently needed for our understanding of the success of this pathogen and development of new diagnostic and therapeutic strategies.

Recent advances in high throughput approaches resulted in the accumulation of new invaluable data on the regulatory strategies of pathogenic bacteria at genomic, transcriptomic and metabolomics levels (Lo et al., 2017; Saliba et al., 2017; Hor et al., 2018). RNA-seq studies notably revealed high complexity of the transcriptome landscape in bacteria (Sorek and Cossart, 2010; Saliba et al., 2017). The precise exploration and interpretation of this large amount of data is essential to improve our understanding of the *C. difficile* infection cycle. One of the first steps towards the establishment of *C. difficile* regulatory pathways is the definition of the transcriptional map at a genome-wide level. In particular, the analysis of transcriptome data combined with *in silico* predictions could provide key information on molecular details of regulatory mechanisms including promoter sequences, type of sigma factor associated to the RNA polymerase (RNAP) involved in the initiation of transcription, as well as other regulatory elements. The genome-wide transcriptional start site (TSS) mapping allows the determination of transcriptional start positions at single-nucleotide resolution (Wurtzel et al., 2010; Sharma and Vogel, 2014). Two methods developed in 2010

have proven their validity for TSS mapping. Differential RNA-seq (dRNA-seq) method largely used in prokaryotes includes two enzymatic steps, i.e., terminator 5'-phosphate-dependent exonuclease (TEX) followed by tobacco acid pyrophosphatase (TAP) treatment for adapter ligation, while 5'-end RNA-seq compares two samples treated or not with TAP enzyme to distinguish between primary 5'-triphosphate and processed 5'-monophosphate transcripts (Sharma et al., 2010; Jager et al., 2014; Sharma and Vogel, 2014; Papenfort et al., 2015; Babski et al., 2016). The 5'-end RNA-seq approach (Wurtzel et al., 2010) that we are using in this study has been successfully implemented for the TSS mapping in Achaea and bacteria including pathogenic and non-pathogenic Listeria (Wurtzel et al., 2012a), *Pseudomonas aeruginosa* (Wurtzel et al., 2012b), *Streptococcus agalactiae* (Rosinski-Chupin et al., 2015), *Streptococcus pyogenes* (Rosinski-Chupin et al., 2019). The detection of TSS implies the presence of a promoter in its upstream region. This promoter element will define the site directing the RNAP for the initiation of transcription and represents a key element to understand the regulation of gene expression in bacteria. Promoters differ at their consensus sequences depending on the interchangeable RNAP sigma factor, which provides DNA recognition specificity (Burgess and Anthony, 2001).

During its infection cycle, *C. difficile* have to face changing conditions including variations in pH, oxygen content, osmolarity and exposure to various antimicrobial compounds (Abt et al., 2016). The *C. difficile* 630 genome carries 14 genes encoding sigma factors including two copies of housekeeping SigA and several alternative sigma factors allowing the initiation of transcription under various physiological conditions (Sebaihia et al., 2006) (Table 1). These sigma factors belong to the sigma 70 family except for one of them, SigL, belonging to the sigma 54 family (Gruber and Gross, 2003). Four specific sigma factors of sporulation, SigF and SigG in the forespore and SigE and SigK in the mother cell and the sigma factor of transition phase, SigH, which is involved in the initiation of sporulation are present in *C. difficile* and their inactivation leads to an asporogenous phenotype (Saujet et al., 2011, 2013; Fimlaid et al., 2013). The sigma factor of the general stress response, SigB, is involved in adaptive strategy during gut colonization and plays a role in resistance to low pH, to antimicrobial peptides, to reactive oxygen species and nitric oxide as well as in oxygen tolerance (Kint et al., 2017). TcdR is an alternative sigma factor required for toxin genes expression (*tcdA* and *tcdB*) while *tcdR* itself is transcribed by RNAP with SigD, which controls flagellar synthesis and motility (Mani and Dupuy, 2001; El Meouche et al., 2013). 3 extracytoplasmic function (ECF) sigma factors, SigV, CsfU and CsfT are also
<table><caption>**TABLE 1** | Sigma factors encoded in the genome of C. difficile strain 630.</caption><thead><tr><th>Gene ID</th><th>Gene name</th><th>Function</th><th>Group*</th><th>% identity / similarity with sigma factors of B. subtilis</th><th>References</th></tr></thead><tbody><tr><td>CD1455</td><td>sigA1</td><td>SigA, exponential growth</td><td>1</td><td>65/77</td><td>Saujet et al., 2011</td></tr><tr><td>CD1498</td><td>sigA2</td><td>SigA, stationary phase</td><td>1</td><td>56/74</td><td>Saujet et al., 2011</td></tr><tr><td>CD0011</td><td>sigB</td><td>SigB, general stress response</td><td>2</td><td>34/54</td><td>Kint et al., 2017</td></tr><tr><td>CD0266</td><td>sigD</td><td>SigD, flagella formation</td><td>3</td><td>30/53</td><td>El Meouche et al., 2013</td></tr><tr><td>CD3176</td><td>sigL</td><td>SigL, sigma 54 factor</td><td>NA</td><td>29/50</td><td>Dubois et al., 2016</td></tr><tr><td>CD0057</td><td>sigH</td><td>SigH factor of transition phase</td><td>3</td><td>63/78</td><td>Saujet et al., 2011</td></tr><tr><td>CD2643</td><td>sigE</td><td>SigE factor of sporulation</td><td>3</td><td>64/77</td><td>Pereira et al., 2013; Saujet et al., 2013</td></tr><tr><td>CD0772</td><td>sigF</td><td>SigF factor of sporulation</td><td>3</td><td>49/72</td><td>Pereira et al., 2013; Saujet et al., 2013</td></tr><tr><td>CD2642</td><td>sigG</td><td>SigG factor of sporulation</td><td>3</td><td>68/86</td><td>Pereira et al., 2013; Saujet et al., 2013</td></tr><tr><td>CD1230</td><td>sigK</td><td>SigK factor of sporulation</td><td>3</td><td>56/77</td><td>Pereira et al., 2013; Saujet et al., 2013</td></tr><tr><td>CD1558</td><td>csfV/sigV</td><td>SigV, ECF sigma factor, resistance to lyzozyme</td><td>4</td><td>57/76</td><td>Ho and Ellermeier, 2011</td></tr><tr><td>CD1887</td><td>csfU</td><td>ECF sigma factor SigW-like</td><td>4</td><td>30/48</td><td>Ho and Ellermeier, 2011</td></tr><tr><td>CD0677</td><td>csfT</td><td>ECF sigma factor</td><td>4</td><td>-</td><td>Ho and Ellermeier, 2011</td></tr><tr><td>CD0659</td><td>tcdR</td><td>Sigma factor of toxin genes</td><td>5</td><td>-</td><td>Mani and Dupuy, 2001</td></tr></tbody></table>

Eight additional Sigma factors are annotated in the genome of C. difficile strain 630 (Sebaiha et al., 2006). The genes CD0359, CD0435, CD0510, CD1094, CD1878, CD2223, CD3329, and CD3371 encoding these sigma factors are associated with conjugative transposon. Group* of σ<sup>70</sup> family sigma factors according to the classification by Gruber and Gross (2003). "NA" not applicable.

present (Ho and Ellermeier, 2011; Sineva et al., 2017). The regulons of several RNAP sigma factors have been recently defined in C. difficile including sporulation specific sigma factors, the general stress-response sigma factor, SigB and the motility-associated sigma factor, SigD (Saujet et al., 2011; El Meouche et al., 2013; Pereira et al., 2013; Saujet et al., 2013; Kint et al., 2017). Only two studies have previously shown a role of SigL, belonging to the Sigma 54 family, in the degradation of cysteine associated with a control of toxin production (Dubois et al., 2016; Gu et al., 2017; Nie et al., 2019). Interestingly, the genes for 24 transcriptional activators (named EBP for enhancer binding proteins activating SigL-dependent promoters) containing a AAA<sup>+</sup> domain, which is responsible for ATP hydrolysis and their interaction with SigL (Francke et al., 2011), are present in the genome of C. difficile. Only two EBPs, CsdR and PrdR, controlling cysteine and proline catabolism, respectively, have been experimentally characterized in C. difficile (Bouillaut et al., 2013; Gu et al., 2017; Gu et al., 2018).

By combining *in silico* analysis, RNA-seq and genome-wide promoter mapping, we have recently identified more than 200 non-coding RNAs (ncRNAs) in C. difficile from different functional classes including riboswitches, trans- and cis-acting antisense RNAs (Soutourina et al., 2013). However, the global view of C. difficile promoter structure is still missing. To fill the gap in our current understanding of C. difficile genes and transcriptional unit structure and regulation, we present here a transcriptional map of the C. difficile genome including the identification of TSSs, operon structures and promoter motifs. We also defined the SigL regulon. These data would considerably improve our knowledge on the regulons of the major sigma factors and on the global regulatory circuits that control gene expression in C. difficile. This work will be essential for genome-wide and

gene-specific studies of the regulatory mechanisms in this emerging pathogen.

## MATERIALS AND METHODS

### Bacterial Strains and Growth Conditions

*Clostridioides difficile* strains were grown anaerobically (5% H<sub>2</sub>, 5% CO<sub>2</sub>, and 90% N<sub>2</sub>) in TY (Bacto tryptone 30 g.l<sup>-1</sup>, yeast extract 20 g.l<sup>-1</sup>, pH 7.4) or Brain Heart Infusion (BHI, Difco). For solid media, agar was added to a final concentration of 17 g.L<sup>-1</sup>. When necessary, cefoxitin (Cfx; 25 µg/ml), erythromycin (Erm; 2.5 µg/ml) and thiamphenicol (Tm; 15 µg/ml) were added to C. difficile cultures. Strains and plasmids used in this study are listed in Supplementary Table S1.

### Volatile Fatty Acid Analysis

The end products of fermentation were detected by gas-liquid chromatography. Strain 630Δ*erm* and the *sigL::erm* mutant were grown in TY for 48 h at 37°C. After centrifugation, supernatants were recovered and mixed with sulfuric acid and ether to extract the volatile fatty acids. For each sample, 5 µl of the supernatants was injected in a gas chromatograph (CP-3380; Varian) equipped with a flame ionization detector and connected to an integrator (model C-55A; Shimadzu). A glass column (2 m by 4 mm) packed with 10% SP 1000 plus 1% H<sub>3</sub>PO<sub>4</sub> on Chromosorb W AW (100/120 mesh) was used. The instrument was operated at 170°C for 18 min. The operating conditions were as follows: injector temperature, 200°C; detector temperature, 200°C; carrier gas (nitrogen); flow rate, 30 ml min<sup>-1</sup>. Each peak of the GLC patterns was identified by the retention time compared with that obtained with standard (2-Methylpentanoic acid at 10 mM). The amounts of fatty acids were calculated by comparison with an internal standard (Carlier and Sellier, 1989).
RNA Extraction, Quantitative Real-Time
PCR and 5' RACE

For the 5'-end RNA-seq experiment, total RNA was isolated
from C. difficile 630Δerm strain grown in TY medium either
after 4 h or 10 h of growth or under starvation conditions that
correspond to a 1 h resuspension of exponentially grown cells (6 h
of growth) into PBS buffer for 1 h at 37°C as previously described
(Andre et al., 2008; Soutourina et al., 2013). 12 ml of culture for
each condition have been used yielding at least 100 µg of total
RNA. The analysis of genes controlled by SigL was performed
with RNA extracted from cells of strain 630Δerm or 630Δerm
sigL::erm mutant (Dubois et al., 2016) after 4 h of growth in
TY. After centrifugation, the culture pellets were resuspended
in RNAproTM solution (MP Biomedicals) and RNA extracted
using the FastRNA Pro Blue Kit, according to the manufacturer's
instructions. The RNA quality was determined using RNA 6000
Nano Reagents (Agilent). Quantitative real-time PCR analysis
was performed as previously described (Saujet et al., 2011). In
each sample, the quantity of cDNAs of a gene was normalized
to the quantity of cDNAs of the dnaF gene (CD1305) encoding
DNA polymerase III. The relative change in gene expression
was recorded as the ratio of normalized target concentrations
(ΔΔCt) (Livak and Schmittgen, 2001). 5' RACE experiments
were performed as previously described (Soutourina et al., 2013).

RNA-Seq and 5'-End RNA-Seq

Non-orientated library for whole transcript analysis by RNA-
seq was realized on a RNA sample extracted from *C. difficile*
630Δ*erm* strain grown at the late-exponential phase (6 h of
growth) as previously described (Wurtzel et al., 2010). For 5'-
end RNA-seq, two strand-specific cDNA libraries were generated
from mixed RNA samples depleted for rRNAs and treated
or not with Tobacco Acid Pyrophosphatase (TAP +/−), as
previously described (Soutourina et al., 2013). TAP converts
5'-triphosphates into 5'-monophosphates allowing sequencing
adaptor ligation and thus the enrichment with primary transcript
reads in the library treated with TAP. The TAP+/− library
construction for high-throughput sequencing (5'-end RNA-seq)
was realized on mixed sample combining RNAs extracted from
three different growth conditions including exponential growth
(4 h), entry to a stationary phase (10 h) and nutriment starvation
(1h incubation in PBS buffer). We prioritized depth of coverage
by combining several conditions in a single 5'-end RNA-seq to
maximize the number of genes expressed and thus the number
of detectable TSS. 15 µg of total RNA treated with TurboDNAse
(Ambion) was used for depletion of ribosomal RNA with
the MicrobExpress kit (Ambion) following the manufacturer
instructions. RNA depleted for rRNA was divided into two
similar fractions and 1500 ng of this partially purified RNA was
used for each library preparation. To convert the 5'PPP ends in
5'P ends, RNA was denatured 10 min at 65°C, placed on ice and
incubated 1 h at 37°C with 10 units of TAP (Epicenter) (TAP+
library). For TAP- library construction, RNA depleted for rRNA
was incubated with buffer alone under the same conditions.
Products were purified by phenol/chloroform extraction and
ethanol precipitation. cDNA library construction for Illumina

sequencing was performed as previously described (Sahr et al.,
2012). The Illumina reads were first scanned using Tagdust
for adaptor removal. To eliminate reads matching the rDNA,
sequence reads were mapped to the rDNA operon sequences
of C. difficile 630 strain using the Bowtie software. Remaining
sequencing reads were mapped to the C. difficile genome as
previously described (Soutourina et al., 2013). We have set up
a semi-automatic procedure for data analysis. The data were
visualized using COV2HTML (Monot et al., 2014) at a strand-
specific manner (for 5'-end RNA-seq libraries) or as a whole
transcript coverage map (for RNA-seq): http://mmonot.eu/
COV2HTML/visualisation.php?str_id=-44. This interface allows
to visualize the accumulation of the first bases of the reads,
which is the signal generated by the 5'-end RNA-seq approach.
The sequences were localized on the genome and formed peaks
(shown in green, see Figures 1-3, 6) that were easily detectable
because the background noise was low. All transcriptional start
signals detected by 5'-end RNA-seq were inspected manually to
identify potential TSS and cleavage sites helped with a score that
is the first base coverage ratio at the given position between
the two conditions TAP+/TAP- (normalized by the number
of total reads). The analysis consisted of scanning the genome
systematically and each time a green peak was detected, we
zoomed in the region to access the surrounding sequences for
detection of promoter motifs. Following criteria have been taken
into consideration for the TSS validation: TAP+/TAP- ratio
(cut-off value of 1.5 in 90% of cases), the identification of
promoter motifs at defined distance upstream of TSS (in silico
promoter prediction score of at least 2 with -10 and -35
boxes positions from TSS separated by 16-18 bp or -12 and
-24 consensus motif positions for the SigL promoters), the TSS
location with respect to the annotated gene (upstream of the RBS
for protein-coding genes) and whole transcript RNA-seq data
coverage. We expected a distribution of RNA-seq coverage signal
from TSS position to the 3'-end of the transcript provided that the
gene was expressed at sufficient level under tested conditions. In
case of previously studied sigma factors, an additional criterion
has been included on the differential expression in comparative
transcriptome analysis between strain inactivated for a given
sigma factor and wild type strain [SigH (Saujet et al., 2011), SigE-
F-G-K (Saujet et al., 2013), SigD (El Meouche et al., 2013), SigB
(Kint et al., 2017), Sigma 54 (present study)].

Then, we also used the PhageTerm software (Garneau et al., 2017) to automatically scan for significant peaks through a statistic module included (Supplementary Table S2). PhageTerm has been developed to statistically detect sequencing bias due to the way phages are packaged. One of these packaging mode (Pac) creates an accumulation of the first bases of sequences identical to the TSS signal. A transcript can be thought of as a small linear genome like a phage and then processed by PhageTerm. First, each gene locus is segmented according to coverage using a regression tree. Then, a gamma distribution is fitted to starting position coverage for each segment and an adjusted *p*-value is computed for each position. This two-parameter method is quite similar the one used in TSSAR, which uses Poisson distribution with a second parameter that depends on the local strength of expression (Amman et al., 2014). The analysis was limited to a
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Visualization of TSS mapping of dual (A) and internal (B) promoters. Representative examples of 5'-end RNA-seq (TAP-/TAP+ profile comparison) and RNA-seq data for dual tandem TSS and internal TSS inside the coding sequences are shown in panel (A) and (B), respectively. Cov2HTML (Monot et al., 2014) was used for the visualization. On a RNA-seq and 5'-end RNA-seq sequence read mapping visualization, coding sequences are indicated by blue arrows. The 5'-end RNA-seq data for either positive "strand +" or negative "strand -" strands are presented in the panels. The TSS identified by 5'-end RNA-seq are indicated by red broken arrows and potential processing sites are indicated by scissors mark. Sigma factor consensus associated with a given TSS is indicated. The TSS corresponds to a position with significantly greater number of reads in TAP+ sample, potential cleavage site corresponds to position with large number of reads in both TAP- and TAP+ samples. 5'-end RNA-seq data show 51-bp reads matching to the 5'-transcript ends, while RNA-seq data show reads covering whole transcript.</figcaption></figure>

region of 500 bases before each gene for process time constraints (160 computing days to scan 3000 genes). PhageTerm detected 856 TSS (Supplementary Figure S1), 74% (633) were already determined by semi-automatic procedure and 26% (233) were new. Of these 233, only 31 contained grounds for inclusion in the list according to TSS criteria described above. PhageTerm analysis results thus overlapped largely with the predictions of our semi-automatic analyses. Most of the remainder were potential cleavage sites. The complete experimental TSS raw data set was deposited in the SRA database with an accession number PRJNA626554.

## Microarray Design, DNA-Array Hybridization and Transcriptome Analysis

The microarray of *C. difficile* 630 genome was designed as previously described (Saujet et al., 2011) (GEO database accession number GPL10556). Transcriptome was performed using four

different RNA preparations for each condition (630Δ*erm* and the *sigL::erm* mutant). Labeled DNA hybridization to microarrays and array scanning were done as previously described (Saujet et al., 2011). The complete experimental data set was deposited in the GEO database with a superseries record accession number GSE149245. All the slides were analyzed using R and limma software (Linear Model for Microarray Data) from Bioconductor project<fnref n="1" />. For each slide, we corrected background with the ‘normexp’ method, resulting in strictly positive values and reducing variability in the log ratios for genes with low levels of hybridization signal. Then, we normalized each slide with the ‘loess’ method (Smyth and Speed, 2003). To test for differential expression, we used the bayesian adjusted t-statistics and we performed a multiple testing correction of Benjamini and Hochberg based on the false discovery rate (FDR) (Benjamini and Hochberg, 1995). For the transcriptome

<footnote n="1">¹www.bioconductor.org</footnote>
<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Examples of cleavage sites detected by TSS-mapping. Representative examples of 5'-end RNA-seq (TAP-/TAP+ profile comparison) and RNA-seq data for potential processing profiles are shown. The RNA-seq and 5'-end RNA-seq data visualization is presented as in **Figure 1**. Potential cleavage site shown by scissors mark corresponds to a position with large number of reads in both TAP- and TAP+ samples. RBS are shown by green boxes to highlight the positions of potential cleavage sites in the proximity or inside RBS.</figcaption></figure>

data, a gene is considered as differentially expressed when the *p*-value is < 0.05.

## Strategy of Identification of Promoters Upstream of TSS

Identification of the consensus sequences for each of the sigma factors corresponding to transcription start sites (TSSs) was made using a modified procedure from Saujet et al. (2013). Positional weight matrices (PWMs) for the different sigma factors SigK, SigH, SigF, SigG, SigD, SigE and SigB were made based on a set of experimentally determined binding sites in *C. difficile* genome (Eckweiler et al., 2018) and for SigA and SigL based on *Bacillus subtilis* sites from DBTBS (Sierro et al., 2008). These PWMs, in addition to nucleotides in the proximal and distant promoter boxes, assign weights to distances between these boxes. Genome regions upstream of TSSs were scanned with these PWMs, the maximal allowed distance from a TSS to the proximal box was set to 15 nucleotides. For each sigma factor the highest scoring candidate site was retained, yielding a matrix of scores TSS × PWM. To account for the fact that PWMs have different distributions of scores, precluding direct comparison, all scores were transformed to Z-scores by subtracting the mean and

dividing by the standard deviation of scores produced by the respective PWM on a set of random sequences. To account for uneven GC-content, for each upstream region, the set of random sequences was constructed by shuffling of its sequence 10000 times. In the case of sigma 54 (SigL), raw scores of sites that did not include known important dinucleotides GG and GC in the distal and proximal boxes, respectively (Nie et al., 2019), were artificially decreased by the maximal achievable score for this matrix for each mismatch. We considered independently two variants of PWMs for SigA, the standard two-box one and the extended one without the distal box, but with invariant TG preceding the proximal box. Finally, for each TSS we ranged sigma factors by their Z-scores and retained two sigma factors with maximal Z-scores as most probable regulators of the TSS.

To validate the assignment of sigma factors to promoters, we used the data on gene expression levels in mutants with knocked out sigma factors: SigH (Saujet et al., 2011), SigE-F-G-K (Saujet et al., 2013), SigD (El Meouche et al., 2013), SigB (Kint et al., 2017), Sigma 54 (present study). Quality of prediction was assessed using a confusion matrix. Rows in the matrix correspond to the actual sigma factors and columns to the *in silico* predicted sigma factors, so the correct predictions are located at the diagonal. A gene was assumed to be actually regulated by
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Examples of promoters controlling regulatory genes detected by TSS-mapping. Representative examples of 5'-end RNA-seq (TAP-/TAP+ profile comparison) data for the identification of TSS for genes encoding important transcriptional regulators are shown. The 5'-end RNA-seq data visualization is presented as in **Figure 1**. The sequence of promoter region is shown upstream of TSS with the -35 and -10 promoter elements indicated in blue and TSS indicated in red.</figcaption></figure>

a sigma factor if its expression level changed by at least 10% in knockout mutants. For each TSS, candidate regulating sigma factors were predicted and ranged as described above. A value in a diagonal cell counts cases when the actual sigma factor was among two best ones predicted (or among three best ones if the two best corresponded to SigA and extended SigA motif). A value in an off-diagonal cell is the count of incorrect predictions with the column corresponding to the predicted sigma factor with the highest score, and the row corresponding to the actual sigma factor (or sigma factors if a gene changed its expression in several mutants).

# RESULTS AND DISCUSSION

## Genome Wide Mapping of Transcriptional Start Sites by 5'-End RNA-Seq

A total of 3684 genes have been annotated in the genome of *C. difficile* strain 630 (Sebaihia et al., 2006; Monot et al., 2011). The main objective of this study was to provide accurate identification of transcriptional start sites (TSS) for a large

number of *C. difficile* 630 transcriptional units. For genome-wide detection of TSS, we performed a differential 5'-end sequencing approach (5'-end RNA-seq). To increase the number of TSS identified, the 5'-end RNA-seq analysis was realized with mixed RNA samples extracted from *C. difficile* 630 Δ*erm* cells harvested during exponential growth phase (4 h), at the onset of stationary phase (10 h) and under nutrient starvation conditions. We pooled together three samples from different growth conditions before sequencing to cover the *C. difficile* transcription map, a strategy largely used for similar studies (Rosinski-Chupin et al., 2015, 2019). The goal was to identify the TSS position and not to analyze the differential gene expression. The lack of replicates could be a problem for genes expressed exclusively in one of tested conditions, however, based on our previous transcriptome analysis and gene-specific studies, we assume that the great part of genes is expressed at least in two or three samples. This study is not exhaustive but provide valuable information about TSS positions for a part of *C. difficile* genes that is useful for further investigations of this pathogen.

This approach allowed us to identify using both the PhageTerm software and manual inspection of sequence data a total of 1562 potential TSS upstream of *C. difficile* strain 630 genes including 274 TSS upstream of non-coding RNA genes
(Soutourina et al., 2013) (**Supplementary Tables S2, S3**). Among TSS associated with ncRNA genes, 66 are located upstream riboswitches, 7 upstream rRNA and 24 upstream tRNA genes (**Supplementary Table S3**). For the cluster of highly expressed tRNA genes, a number of suggested TSS could correspond to cleavage sites representing processing events during transcript maturation process.

Together with automatic *in silico* approaches (Amman et al., 2014; Jorjani and Zavolan, 2014), manual data curation still remains valuable for gene annotation and data interpretation and validation. For the TSS mapping data analysis we have used a semi-automatic procedure. From the 5-end RNA-seq data the accumulation of the first bases of the reads has been visualized as green peaks in our visualization interface. The genome has been then systematically scanned and each green peak has been inspected for surrounding sequences to detect transcription markers. We have used as a basis an automatic TAP+/TAP - ratio score at the target position representing the first base coverage of TAP+ reads divided by the TAP- coverage of the first base (normalized by the number of total reads) (**Supplementary Tables S2, S3**). We then combined several additional criteria for validation of a potential TSS. These include the presence of well-located promoter motifs upstream of TSS with *in silico* promoter prediction score (see paragraph on *in silico* prediction of promoter motifs associated with different sigma factors), whole transcript RNA-seq data for transcript coverage, differential expression in comparative transcriptome analysis between strain inactivated for a given sigma factor and wild type strain when available.

A total number of 1288 protein-coding genes could be associated with a potential TSS through our analysis. The genomic position of all identified TSSs with corresponding genes is listed in **Supplementary Table S2**. A large proportion of these TSS corresponds to the first gene of *C. difficile* operons. We have then included the data on possible operon structure for transcriptional unit organization analysis from available gene annotation (Vallenet et al., 2020) and whole RNA-seq expression profiles leading to a total number of more than 2000 genes covered by our TSS mapping (**Supplementary Table S2**). Additional genes missing in our analysis could be expressed under conditions different from those used in this study, and their TSS identification would require specific expression conditions.

As an example of a large gene cluster, **Supplementary Figure S2** depicts the operon map and the promoters of the genes coding for flagella biosynthesis and function. Both one SigA-dependent and three alternative sigma factor SigD-specific promoters could be associated with TSSs for flagella genes (**Supplementary Table S2** and **Supplementary Figure S2**). In the *C. difficile* 630 genome, three loci encode flagellum-associated proteins : late stage flagellar genes CD0226-CD0240, flagellar glycosylation genes CD0241-CD0244 and early stage flagellar genes CD0245-CD0272 (Stabler et al., 2006; Aubry et al., 2012; Stevenson et al., 2015). The hierarchy of flagellar transcription starts with early stage CD0245/flgB flagellar operon that contains genes for assembly of the basal body and for the flagellar alternative sigma factor SigD that activates the transcription of the late stage operons. These late stage flagellar genes are involved

in assembly of the flagellar hook, filament, and cap, and for post-
translational modification of the flagellar filament. We compared
our results of TSSs and associated promoter identification with
previously suggested transcriptional organization of flagellar
cluster and completed its TSS map (**Supplementary Figure S2**
and **Supplementary Table S2**, El Meouche et al., 2013; Stevenson
et al., 2015). In accordance with previous data, we identified
a TSS associated with a SigA-dependent promoter for early
stage *flgB* operon and TSSs associated with SigD-dependent
promoters for late stage genes CD0226 encoding putative lytic
transglycosylase, *flgM* (CD0229) encoding flagellar anti-sigma
factor and *fliC* (CD0239) encoding flagellin (**Supplementary
Table S2**). Additional TSS could be defined for CD0241,
*fliQ* (CD0261) and *flgG* (CD0269) genes. Similarly, to the
previously reported antisense TSS in the flagellar biosynthesis
cluster of *Legionella pneumophila* (Sahr et al., 2012), several
TSS could be found in antisense orientation to flagellar genes
corresponding to CD630_n00050 – CD630_n00120 antisense
RNAs (Soutourina et al., 2013, **Supplementary Figure S2** and
**Supplementary Table S3**).

To further validate our approach for TSS mapping, we
compared our data with about 35 TSSs previously mapped by
5' rapid amplification of cDNA ends (5'RACE) or 3'/5' RACE
approach (Table 2). We also performed 5'RACE experiments
that unambiguously identified TSSs for two additional genes
(sinR encoding a transcriptional regulator and tcdC encoding
protein negatively controlling toxin gene expression). Almost all
these TSSs were in agreement with the TSSs identified by 5'-
end RNA-seq with only minor deviations in case of multiple
possible TSS and/or potential cleavage sites detected. These
results and the TSSs already characterized upstream potential
ncRNA genes confirmed that we have obtained a robust dataset
that describes the transcriptional map of the C. difficile strain 630
(Supplementary Table S3, Soutourina et al., 2013).

The analysis of the nucleotide composition of TSS identified
in this study revealed a strong selection for purine with the
majority of A (75%) generally required for efficient initiation
of transcription by RNAP and to a lower extent of G (14%)
and T (9%) (Supplementary Figure S3). The TSS mapping
also allowed us to clarify the position of translational start
sites for 10 genes. Indeed, the translational start sites of several
genes are located upstream of the TSSs suggesting a mis-
annotation for translation initiation. We carefully checked the
CDS and identified start codons (ATG, TTG or GTG) with
a ribosome-binding site located upstream. For example, we
modified the translational start sites for the CD0341, CD1088,
CD2055, CD2064, CD2112, CD3271 genes encoding conserved
proteins of unknown function, CD1015 and CD1099 genes
encoding two-component system response regulator or histidine
kinase, CD1234 encoding a small protein contributing to skin
element excision during sporulation and msrAB (CD2166) gene
encoding peptide methionine sulfoxide reductase. The modified
annotations are available on MaGe MicroScope platform
(Vallenet et al., 2020).

The 5' untranslated region (5'UTR) also called mRNA leader
region located between TSS and translation initiation codon
constitutes often a target for important regulatory processes.
<table><caption>**TABLE 2** | Comparison of transcriptional start positions identified by RACE and TSS mapping by 5'-end RNA-seq.</caption><thead><tr><th>Gene</th><th>5' RACE position</th><th>5'-end RNA-seq<br>5' position</th><th>Orientation</th><th>Promoter</th><th>References</th></tr></thead><tbody><tr><td>tcdC CD630_06640</td><td>T 805066<br>A 805037<br>T 805031<br>C 805030<br>A 805028<br>A 805027<br>A 805022<br>T 805505<br>T 805503<br>A 805481</td><td>T 805066</td><td>-</td><td></td><td>This work</td></tr><tr><td>tcdR CD630_06590</td><td>A 784644<br>A 786446<br>T 786504<br>A 786505<br>A 786379<br>A 786239</td><td>A 786444</td><td>+</td><td>P-SigD</td><td>El Meouche et al., 2013</td></tr><tr><td>sinR CD630_22140</td><td>G 2566722</td><td>G 2566722</td><td>+</td><td>P-SigA</td><td>This work</td></tr><tr><td>sigh CD630_00570</td><td>A 82971</td><td>A 82971</td><td>+</td><td>P-SigA</td><td>Saujet et al., 2011</td></tr><tr><td>spo0A CD630_12140</td><td>T 1412457<br>G 1412532</td><td>G 1412458<br>A 1412531</td><td>+</td><td>P-SigA</td><td>Saujet et al., 2011</td></tr><tr><td>spollAA CD630_07700</td><td>T 942481</td><td>G 942480</td><td>+</td><td>P-SigH</td><td>Saujet et al., 2011</td></tr><tr><td>CD630_24920</td><td>A 2877271</td><td>not mapped</td><td>-</td><td>P-SigH</td><td>Saujet et al., 2011</td></tr><tr><td>dnaG CD630_14540</td><td>A 1682862</td><td>not mapped</td><td>+</td><td>P-SigA</td><td>Saujet et al., 2011</td></tr><tr><td>sigA2 CD630_14980</td><td>G 1734857</td><td>A 1734856</td><td>+</td><td>P-SigH</td><td>Saujet et al., 2011</td></tr><tr><td>sigG CD630_26420</td><td>T 3050599</td><td>A 3050598</td><td>-</td><td>P-SigF</td><td>Saujet et al., 2011</td></tr><tr><td>CD630_02260</td><td>G 292850</td><td>G 292850</td><td>+</td><td>P-SigD</td><td>El Meouche et al., 2013</td></tr><tr><td>flgm CD630_02290</td><td>A 294568</td><td>A 294568</td><td>+</td><td>P-SigD</td><td>El Meouche et al., 2013</td></tr><tr><td>fliC CD630_02390</td><td>G 300874</td><td>G 300874</td><td>+</td><td>P-SigD</td><td>El Meouche et al., 2013</td></tr><tr><td>CD630_35270</td><td>G 4122554</td><td>G 4122554</td><td>-</td><td>P-SigD</td><td>El Meouche et al., 2013</td></tr><tr><td>CD630_23440</td><td>A 2713978</td><td>A 2713981</td><td>-</td><td>P-SigA</td><td>Dineen et al., 2010</td></tr><tr><td>glgC CD630_08820</td><td>A 1059909</td><td>T 1060072</td><td>+</td><td>P-SigA</td><td>Dineen et al., 2010</td></tr><tr><td>hisZ CD630_15470</td><td>around 1795492</td><td>A 1795496</td><td>+</td><td>P-SigA</td><td>Dineen et al., 2010</td></tr><tr><td>cdsB CD630_32320</td><td>A 3784442</td><td>A 3784442</td><td>-</td><td>P-SigL</td><td>Gu et al., 2018</td></tr><tr><td>cwpV CD630_05140</td><td>A 607231</td><td>A 607231</td><td>+</td><td>P-SigA</td><td>Emerson et al., 2009</td></tr><tr><td>dltD CD630_28540</td><td>A 3337889<br>A 3337884</td><td>A 3337889</td><td>-</td><td>P-SigA/P-SigV<br>P-SigA/P-SigV</td><td>Woods et al., 2016</td></tr><tr><td>CD630_25171</td><td>T 2908422<br>A 2908421</td><td>T 2908422</td><td>-</td><td>P-SigA/P-SigB</td><td>Maikova et al., 2018</td></tr><tr><td>CD630_29071</td><td>A 3398599</td><td>A 3398599</td><td>-</td><td>P-SigA/P-SigB</td><td>Maikova et al., 2018</td></tr><tr><td>CD630_09562</td><td>A 1124042</td><td>A 1124042</td><td>+</td><td>P-SigA/P-SigB</td><td>Maikova et al., 2018</td></tr><tr><td>SQ1781 (RCd8)</td><td>T 2907991<br>T 2908066<br>A 2907896</td><td>A 2908006<br>T 2908013</td><td>+</td><td></td><td>Maikova et al., 2018</td></tr><tr><td>CD630_n00370 (RCd10)</td><td>A 1124339</td><td>A 1124339</td><td>-</td><td>P-SigA/P-SigB</td><td>Maikova et al., 2018</td></tr><tr><td>CD630_n01000 (RCd9)</td><td>A 3398302</td><td>A 3398302</td><td>+</td><td>P-SigA/P-SigB</td><td>Maikova et al., 2018</td></tr><tr><td>CD630_n00860</td><td>A 2908058</td><td>A 2908058</td><td>-</td><td></td><td>Maikova et al., 2018</td></tr><tr><td>SQ2025</td><td>A 3306816<br>T 3306807<br>A 3306797<br>G 3306788</td><td>cleavage T 3306807</td><td>-</td><td></td><td>Soutourina et al., 2013</td></tr><tr><td>CD630_n00210 (RCd4)</td><td>A 655066<br>A 655075</td><td>A 655072</td><td>+</td><td></td><td>Soutourina et al., 2013</td></tr></tbody></table>

(Continued)
<table><caption>**TABLE 2** | Continued</caption><thead><tr><th>Gene</th><th>5' RACE position</th><th>5'-end RNA-seq<br>5' position</th><th>Orientation</th><th>Promoter</th><th>References</th></tr></thead><tbody><tr><td>CD630_n00680 (RCd5)</td><td>A 2285913</td><td>G 655119<br>A 2285913</td><td>+</td><td></td><td>Soutourina et al., 2013</td></tr><tr><td>SQ1002</td><td>A 1761105<br>T 1761106</td><td>A 1761105<br>T 1760987<br>T 1761212</td><td>-</td><td></td><td>Soutourina et al., 2013</td></tr><tr><td>SQ173</td><td>T 308770</td><td>T 308776</td><td>+</td><td></td><td>Soutourina et al., 2013</td></tr><tr><td>SQ1498</td><td>C 2441928<br>A 2441933</td><td>T 2441927<br>A 241078</td><td>-</td><td></td><td>Soutourina et al., 2013</td></tr><tr><td>CD630_n00030 (RCd2)</td><td>T 241079<br>A 241065<br>T 241067</td><td>A 241078</td><td>-</td><td></td><td>Soutourina et al., 2013</td></tr><tr><td>CD630_n00170 (RCd6)</td><td>A 560340<br>C 560320</td><td>A 560340</td><td>-</td><td></td><td>Soutourina et al., 2013</td></tr></tbody></table>

This region could carry conserved motifs for RNA-based regulations or binding sites for regulatory proteins that could be involved in various transcriptional and post-transcriptional regulatory mechanisms including repression or attenuation of transcription, mRNA stabilization and degradation, riboswitch-based premature termination of transcription and control of initiation of translation. The length of these regions is usually related to their functional importance in gene expression regulatory processes. We analyzed the length distribution of the 5'UTR of the genes with identified TSS (Supplementary Figure S4). For the majority of *C. difficile* mRNAs (75% of the TSSs mapped), the first predicted translational start lies within 100-bp region downstream from the TSS. A total of 392 genes had a TSS located within 20-40 bp from the translational start site, and only few potentially leaderless genes (15) could be found. The presence of leaderless mRNAs has been reported in prokaryotes associated with atypical mechanisms of translation initiation (Jager et al., 2009; Nakagawa et al., 2010). By contrast, 137 mRNAs had leaders longer than 300 bp that could be associated with the presence of particular regulatory motifs. Most of these genes are regulated by specific riboswitches. For example, a number of long-leader containing genes are associated with specific T-boxes responsive to the level of tRNA aminoacylation, among them are genes encoding specific tRNA synthetases, e.g., *serS1*, *metG*, *proS*, *valS* and *asnC*, as well as amino acid biosynthesis and transport genes, e.g., *trpP*, *leuA*, *thrC* and *argH*. The large *flgB* flagella operon is associated with a 496-bp 5'UTR region carrying a c-di-GMP-responsive riboswitch (Sudarsan et al., 2008; Soutourina et al., 2013) and a “flagellar genetic switch” that controls the phase variable production of flagella (Anjuwon-Foster and Tamayo, 2017).

## Identification of Two Promoters Upstream Genes and Internal Promoters

For a total of 83 genes, we detected the possible existence of two TSSs (Supplementary Table S2 column A highlighted in green, Figure 1A, Supplementary Figure S5A). In the

majority of cases, these additional TSSs were associated with lower signal from deep-sequencing data and either the same or a different sigma factor recognition motif could be found upstream of the TSSs. We could hypothesize that for these genes two alternative promoters could be used depending on the environmental factors and growth conditions. For example, *rpmH* (*CD3680*) gene encoding 50S ribosomal protein L34 is associated with one *SigF*- and one *SigA*-dependent promoter having different deep-sequencing signal intensity for both 5-end RNA-seq and whole transcript RNA-seq (Figure 1A). Likewise, for *cbpA* (*CD3145*) gene encoding surface-exposed adhesin, a first TSS potentially associated with *SigA*-dependent promoter is detected with lower signal intensity while a second *SigA*-dependent TSS is revealed with higher sequence reads number (Figure 1A). Interestingly, at least two potential cleavage sites could be found in 5'-end of *cbpA* CDS with high intensity signal from both TAP-treated and non-treated samples. In addition, potential internal TSS could be detected inside the coding part of this gene (Figure 1A). Two different promoters associated with different sigma factors were detected upstream of the *CD0761* gene encoding a putative ATP-dependent RNA helicase and *mreB* (*CD1145*) gene encoding rod-shape determining protein *MreB* (Supplementary Figure S5A). In these cases, either a *SigF*- or a *SigH*-dependent promoter is present in addition to a *SigA*-dependent promoter (Supplementary Figure S5A and Supplementary Table S2).

A number of internal TSSs in the coding regions of genes was also detected, which further demonstrated a complex transcriptional architecture of *C. difficile* genome (Supplementary Table S2, column I highlighted in blue, Figure 1B). In many cases, these internal TSS map to the genes with identified primary TSS in accordance with the mechanisms for internal transcription initiation by elongating RNAP complexes suggested in bacteria (Shao et al., 2014; Cuklina et al., 2016; Harden et al., 2016). Moreover, internal transcription initiation could be at the origin of 3'-end derived ncRNA as recently demonstrated in other bacterial species
(Chao et al., 2012; Guo et al., 2014). In addition to a primary TSS identified upstream of the CD0156 gene, we could find a TSS associated with -10 box for SigA-dependent promoter in the 3'-end of this gene encoding a putative membrane protein while no known downstream gene could be identified in sense orientation (Figure 1B). Similarly, the TSS for coaE (CD1129) gene encoding dephospho-CoA kinase lies in the middle of the upstream polA (CD1128) gene encoding DNA polymerase I (Supplementary Figure S5B). In addition to primary TSS upstream of CD3045 gene encoding putative ATPase, two internal promoters could be identified inside this gene that could drive the transcription of downstream CD3044 gene encoding a transcriptional regulator of the RpiR family (Supplementary Figure S5B). This is an example of internal TSS that could be found inside the upstream gene co-transcribed with the downstream gene within the same operon. The presence of these internal promoters might allow to change the ratio of protein production of co-transcribed genes, a differential expression of genes inside operons or more complex regulatory processes related to alternative operon structure and co-transcriptional relationship between genes (Sorek and Cossart, 2010; Shao et al., 2014).

## Identification of Cleavage Sites Within Transcripts

In addition to TSS mapping, the comparison of RNA-seq signal intensity between TAP+ and TAP- samples allowed us to suggest several potential cleavage sites corresponding to positions with large number of reads in both samples treated or not with TAP enzyme (TAP+/TAP- ratio close to 1), located downstream from identified TSS and not associated with promoter consensus sequences from bioinformatics search. Examples of these sites potentially cleaved by cellular ribonucleases are shown in Figure 2 and Supplementary Figure S6. Different ribonucleases such as RNase III, RNase Y or RNase J could be involved in these potential processing events and could recognize particular secondary structures rather than specific sequence motifs (Trinquier et al., 2020). In the majority of cases, the potential cleavage sites are situated just downstream from the TSS at 23-62 bp and up to 87 bp distance and from 0 to 54 bp upstream of translation initiation RBS site (Figure 2). The processed mRNA would thus retain in the majority of cases the RBS and codon for translation initiation. For the CD0564 gene encoding ATP-dependent protease, the modA (CD0869) gene encoding molybdenum ABC-type transporter, the gcvTPA (CD1657) glycine decarboxylase operon and the CD2396 gene encoding conserved protein of unknown function the potential cleavage site is situated at 0 to 3 bp from RBS (Figure 2). By contrast, for the CD1392 gene encoding a putative ribonuclease, the potential cleavage site lies at "A" position inside the "GGAGG" RBS (Supplementary Figure S6). Sometimes, it is difficult to distinguish between potential cleavage sites and the presence of additional TSS based on TSS mapping data visualization (Figure 2, Supplementary Figure S6 and data not shown). Conversely, a part of potential TSS (Supplementary Tables S2, S3) could correspond to

cleavage sites. A complex profile (Supplementary Figure S6) could be observed with several TSS and potential cleavage sites associated with the same gene, e.g., for the CD3145 adhesin-encoding gene (Figure 1A). CD1329 gene encoding RNaseY, CD1392 gene encoding a putative ribonuclease and ptsH (CD2756) gene encoding HPr protein of PTS system also presented a complex profile with several potential cleavage sites inside the gene (Supplementary Figure S6). Interestingly, a potential cleavage site could be also detected in the 3'-part of the agrB gene encoding accessory gene regulator of a quorum-sensing system (Supplementary Figure S6). This gene is co-transcribed with autoinducer prepeptide agrD gene. The cleavage of the bicistronic agrB-agrD mRNA will result in an agrD transcript and could lead to the stabilization of the processed transcript.

## Transcription in the Pathogenicity Locus and Key Regulatory Genes

In *C. difficile*, the main virulence factors, toxins TcdA and TcdB, as well as toxin expression regulatory components TcdR, TcdC and TcdE are encoded in the 19.6 kb chromosomal region, called pathogenicity locus (PaLoc) (Braun et al., 1996; Martin-Verstraete et al., 2016). TcdR is an alternative sigma factor (Table 1) driving the transcription of tcdA, tcdB and tcdR genes, while TcdC negatively controls toxin gene expression by interfering with transcription initiation by TcdR-containing RNAP (Mani and Dupuy, 2001; Matamouros et al., 2007). Despite a relatively low level of PaLoc region gene expression, the analysis of sequencing data allowed us to define a TSS for both tcdR and tcdC regulatory genes, that were in accordance with RACE analysis previously published for tcdR (El Meouche et al., 2013) or performed in this study for tcdC (Table 2). The absence of detection of TSSs upstream of tcdA and tcdB that have been previously mapped (Dupuy and Sonenshein, 1998) might be due to the low level of transcription of these genes under the conditions used.

**Figure 3 and Supplementary Figure S7** show representative examples of TSS mapping for several genes encoding regulators such as CcpA, CodY, SigH, RstA, LuxS, LexA, and SinR known to directly or indirectly control toxin gene expression (Martin-Verstraete et al., 2016). Most of these regulatory genes are driven by SigA-dependent promoters. In addition to primary TSS upstream of the codY gene, a second TSS, which is potentially associated with a SigF-dependent promoter could be identified inside the upstream gene (Supplementary Figure S7). The comparison of TAP-treated and non-treated samples also suggested the presence of potential cleavages sites just downstream from the TSS of ccpA and lexA genes, and in the 3'-end of codY and fur genes. Interestingly, a cleavage site could be detected between cggR (CD3175) and gapA (CD3174) genes as observed for the polycistronic cggR-gapA operon maturated by RNase Y in B. subtilis and S. aureus (DeLoughery et al., 2018) (Figure 3). This processing event could also lead to a differential stability and thus discordant
abundance of co-transcribed genes encoding glycolytic enzyme and its transcriptional regulator in C. difficile.

## In silico Prediction of Promoter Motifs Associated With Different Sigma Factors

14 genes encoding sigma factors (two copies of SigA, SigH, SigF, SigE, SigG, SigK, SigB, SigD, TcdR, SigV, Csfu, Csft and SigL) are present in the C. difficile 630 genome (Table 1). In addition, eight potential sigma factors associated with mobile genetic elements mainly in conjugative transposons of Tn916, Tn1549 or Tn5397 family are annotated in the genome (Sebaihia et al., 2006; Brouwer et al., 2011). These proteins might be involved in the transcription of genes inside these mobile genetic elements but we cannot exclude that they also contribute to expression of genes outside the transposons. The relatively high number of sigma factors suggests a very complex pattern of transcription in C. difficile. To characterize potential promoter consensus sequences, we analyzed the nucleotide sequences upstream of the TSSs. We extracted the 100 nucleotides (from -99 to +1) upstream of the TSS for protein-coding genes and ncRNA genes and analyzed them for the presence of motifs previously identified for a smaller set of C. difficile promoters under the control of SigB, SigH, SigE, SigF, SigG, SigK, TcdR and SigD or B. subtilis promoters (SigA, SigL) (Mani et al., 2002; Saujet et al., 2011, 2013; El Meouche et al., 2013; Kint et al., 2017). By contrast, consensus sequences for the three ECF sigma factors and for the putative sigma factors associated with mobile genetic elements in C. difficile remain to be identified and we cannot search for the corresponding promoters upstream of TSS.

The automatic search for promoters upstream of TSS is known to be difficult due to variations in the distance between -10 and -35 boxes or between the TSS and the -10 element and sometimes degenerated consensus sequences. The AT-rich nature of C. difficile genome especially in intergenic regions (Supplementary Figure S3) also complicates the search for promoters. To identify sigma factors corresponding to TSSs as reliably as possible, positional weight matrices (PWMs) for sigma factors were made on the basis of a set of binding sites in C. difficile genome that are identified by comparing expression profiles of WT strain versus mutants inactivated for each sigma factor. Data were previously obtained for SigK, SigH, SigF, SigG, SigD, SigE, SigB and TcdR (Mani et al., 2002; Saujet et al., 2011, 2013; El Meouche et al., 2013; Fimlaid et al., 2013; Kint et al., 2017). Concerning SigL, we performed a transcriptomic analysis (see below). For SigA, the consensus was based on B. subtilis sites from database of transcriptional regulation in B. subtilis (DBTBS) (Sierro et al., 2008). After training (Supplementary Figure S8), we applied the obtained positional weight matrices to genes that have demonstrated differential expression in mutants for each sigma factor. The results are shown in Figure 4 and Supplementary Table S4. The predictions for sporulation specific sigma factors (SigF, SigE, SigG, SigK) largely matched the experimental data but allowed to slightly extend each regulon (see below). The automatic predictions for sigma factors with larger regulons (SigL, SigB, SigH) were relatively less reliable. In particular, SigA

was predicted as the main factor for many of them. It could be caused by two reasons. Firstly, these promoters could be recognized by two sigma factors, SigA and another sigma factor. Secondly, some genes could indeed be transcribed only by SigA, and their expression could change in the other sig mutants due to downstream global effects on other genes involving regulatory cascades.

The motifs for SigA-dependent promoter could be found closely upstream of TSS for 743 genes in C. difficile genome as best prediction (Supplementary Tables S2, S3). An extended -10 box with consensus sequence TGNTATAAT could be identified upstream of 245 TSSs (Supplementary Tables S2, S3) (Browning and Busby, 2016). For these genes, no -35 consensus or only a weak consensus could be usually found in promoter region. Extended -10 promoters do not require a strong -35 sequence for efficient transcription (Keilty and Rosenberg, 1987; Gruber and Gross, 2003). Indeed, the extra TG dinucleotide stabilizes open complex formation by providing critical contacts with region 3.0 of RNAP (Gruber and Gross, 2003; Murakami and Darst, 2003). In addition, for several promoters, it is difficult to discriminate between extended -10 promoters and SigB-dependent promoters using the consensus recently defined in C. difficile (Kint et al., 2017) (Supplementary Tables S2, S3, S4). Finally, some genes controlled by transcriptional activators have only weak consensus motifs at -35 position or even no -35 box and we found a series of promoters with only a -10 box (Supplementary Tables S2, S3, S4).

Based on the PWMs, previously reported transcriptomic data and manual inspection of promoters, we assigned a corresponding sigma factor involved in the transcription for most of the genes (Supplementary Tables S2, S3) and we refined the respective consensus sequences of promoters recognized by individual RNAP associated with specific sigma factors (Figure 5). We will discuss below in more detail the results obtained for SigH, the sigma factors specific of sporulation and SigL.

In accordance to our previous analysis (Soutourina et al., 2013), the majority of ncRNA genes could be associated with SigA-dependent promoters. Upstream of a few ncRNA TSS, we predicted the presence of promoters probably recognized by SigD, SigB or sporulation-specific sigma factors (Supplementary Tables S3, S4).

## Genes Transcribed by Sigma Factors of the Transition Phase and of Sporulation

Transition from exponential growth toward stationary phase represents a key checkpoint in C. difficile development associated with the induction of toxin production and decisions to initiate sporulation or biofilm formation. By comparing the transcriptome of the C. difficile strain 630Δ*erm* and its isogenic sigH mutant, we have previously identified 286 genes, which are positively controlled by SigH (Saujet et al., 2011). To identify genes likely transcribed by RNAP associated with SigH, we combined the transcriptome analysis previously published (Saujet et al., 2011) and the *in silico* analysis of promoters
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Confusion matrix for the prediction of the sites recognized by the main sigma factors of C. difficile. Rows: actual sigma factors, based on changes in gene expression in knockout mutants. Columns: predicted sigma factors. A prediction was listed as correct if the actual sigma factor was among two best ones, and incorrect otherwise, in the latter case the top-scoring sigma factor was listed as the prediction.</figcaption></figure>

(Table 3, Supplementary Tables S2, S4). We can then propose a list of 41 genes down-regulated in transcriptome in the sigH mutant with a consensus for SigH-dependent promoters identified as best hit (Table 3). Two promoters having a lower score of promoter identification were not used to refine the consensus (Figure 5). 12 additional genes with a consensus for SigH as best hit are not controlled by SigH in transcriptome in our conditions, two of these genes having a second mapped promoter. In addition to genes involved in sporulation (spo0A, spoIIAA, spoVD, spoVS, spoVG and CD0572) (Saujet et al., 2011) and cell division (minC, ftsZ, CD2650 and CD3673) (Saujet et al., 2011), we found genes involved in translation (fusA1, frr) (Sebaihia et al., 2006), in envelope-related processes (cwp27, CD0573, lplA, CD3458) (Sebaihia et al., 2006), in metabolism (CD0865, fdhF, CD0999, CD1484 and CD2989) (Neumann-Schaal et al., 2019) and genes encoding 11 conserved hypothetical proteins of unknown function. Interestingly, two genes of the CD0999 and CD2989 operons encoding sulfonate ABC transporters as well as a gene CD1484 encoding a sulfonate binding protein were expressed under the control of a P-SigH promoter. This suggests that these sulfur-containing compounds play a key role at the onset of stationary phase as alternative sulfur-sources. It is worth noting that ssuA is the last gene of an operon encoding another sulfonate ABC transporter. The first two genes of this operon were not controlled

by SigH. We detected an internal promoter located within CD1483 with a consensus recognized by SigH in agreement with the specific transcriptional control of ssuA by SigH and the increased transcript level observed for this gene in RNA-seq (Figure 6A). An internal promoter recognized by SigH is also present upstream of CD0789 in the CD0788-CD0789 operon (Figure 6B).

Spore formation is a tightly controlled process that constitutes an essential step in C. difficile life cycle for its dissemination and survival. We have previously compared our TSS mapping data (Soutourina et al., 2013) with transcriptomic data done with microarrays comparing wild-type strain with mutants inactivated for sigF, sigE, sigG or sigK genes encoding major sporulation sigma factors at different times during sporulation in liquid media (Saujet et al., 2013). We now extended this comparison to RNA-seq data obtained using RNA extracted after 18 h of growth on plates by Aimee Shen and co-workers (Fimlaid et al., 2013; Pishdadian et al., 2015). Even if the microarray and RNA-seq data largely overlapped, this new analysis allowed us to propose a more complete list of genes transcribed by the sigma factors of the forespore, SigF and SigG (Supplementary Table S5) or the mother cell, SigE and SigK (Supplementary Table S6) and to refine the consensus of the promoters recognized by each of these sigma factors (Figure 5). We identified 10 new SigE-dependent promoters
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Consensus sequences for promoters recognized by the different sigma factors. These consensus were generated by WebLogo using the data presented in **Supplementary Tables S2, S4–S6** and **Table 3**. The consensus for TcdR (**Table 1**) was obtained from characterized promoters (Mani et al., 2002; Martin-Verstraete et al., 2016).</figcaption></figure>

(CD0557, CD0629, CD1022, CD1043, CD1403, CD1575, CD1629, CD2316, CD2443, CD3466) while only two additional SigK-dependent promoters (CD2167 and CD2720) can be proposed (**Supplementary Table S6**). For *glmS* and *CD1340* that are controlled by SigK, both a SigK- (best hit) and a SigA-dependent promoter are detected upstream of these genes. In addition, for four genes (*CD1511*, *CD1930*, *CD2443* and *CD3350*), the best hit detected in *silico* did not fit with the transcriptome data. In the forespore (**Supplementary Table S5**), we identified 56 promoters recognized by SigF and/or SigG, including 12 new ones compared to our previous analysis (Saujet et al., 2013). 19 and 16 genes can be assigned to the SigF or the SigG regulon, respectively. 15 additional promoters have the key motifs recognized by SigF (**Supplementary Table S2**) but we failed to detect forespore-dependent control (SigF- or SigG-dependent) of their expression. For six of them, this absence of control might be due to the presence of a second promoter upstream of the gene. For 21 additional forespore-controlled genes, our data are ambiguous (**Supplementary Table S5**). Although a consensus closer to those recognized by SigG than by SigF was identified upstream of the TSS, the expression of CD2266, CD2599 and CD2856 was only controlled by SigF. Unexpectedly, CD1789, CD2245.1 and *sigG* were controlled by SigG but not by SigF despite the fact

that a consensus recognized by SigF is proposed upstream of these genes (**Supplementary Tables S5, S4**). In *C. difficile*, *sigG* is transcribed from at least two promoters, one in front of *spoIIGA* recognized by SigA, and the other just upstream of *sigG* recognized by forespore specific sigma factors (Saujet et al., 2013). Interestingly, *CD2245.1* is located downstream from the *cspBAC* operon, a member of the SigE regulon and is also less expressed in the *sigE* mutant than in the WT strain in transcriptome contrary to most forespore-controlled genes (Saujet et al., 2013). It is tempting to speculate that a more complex profile of expression with a possible readthrough from *cspBAC* also exists for *CD2245.1*. Finally, 15 genes controlled both by SigF and SigG as expected for members of the SigG regulon have a consensus closer to those of SigF-dependent promoters than SigG-dependent promoters. This included genes encoding 3 SASPs (*CD1290*, *CD3220.1*, *CD3249*), a catalase (*CD1567*), a dipicolinate transporter (*spoVAC* operon) and the SpoVT regulator. The results obtained for *spoVT* strongly suggest that this gene is transcribed by both SigF and SigG (Fimlaid et al., 2013; Saujet et al., 2013). The existence of a residual transcription of *spoIIR* in the *sigF* mutant in the forespore also suggests that SigG might allow *spoIIR* transcription. Due to the promiscuity of the consensus of SigF and SigG-controlled promoters (**Figure 5**),
<table><caption>TABLE 3 | Promoters transcribed by SigH associated to the RNAP.</caption><thead><tr><th>Gene</th><th>Function</th><th>Expression ratio<sup>a</sup> sigH/630Δerm</th><th>Promoter recognized by SigH</th><th>Score promoter</th></tr></thead><tbody><tr><td>CD0022 fusA1</td><td>Elongation factor G</td><td>0.04</td><td>AAAAGAGGACAA TTA C C T C C A GAT GTAGAAAT A A A G C T A T G A</td><td>4.29</td></tr><tr><td>CD0142#</td><td>RNA-binding protein</td><td>0.14</td><td>AAAAGAGGATTA T G A G A T T C G T GTAGAATATATTAT</td><td>4.30</td></tr><tr><td>CD0148</td><td>Conserved hypothetical protein</td><td>0.42</td><td>GAAAAAAGGTATCTACAATATACAGAGCGAATTTATAATA</td><td>4.16</td></tr><tr><td>CD0440 cwp27</td><td>Cell wall binding protein</td><td>0.25</td><td>CATAAAGGAAAA TAT T C T T T T T AT GTAGAATTA T A T A T TTAAGTAA</td><td>5.82</td></tr><tr><td>CD0572#</td><td>Sporulation protein</td><td>0.49</td><td>GATAAAGGTATTTGGAAAAGTATGAAGAATACAGAAATAGA</td><td>3.64</td></tr><tr><td>CD0573*</td><td>Membrane protein</td><td>0.45</td><td>ATAAAGCTATCCATATAAAATTAAAAAAATTCTTTTTA</td><td>1.57</td></tr><tr><td>CD0770 spoIIA</td><td>Anti-SigF factor</td><td>0.03</td><td>ATTGAAGGAATAAAAAATAAATTATAGAATTGATTTAAAG</td><td>4.93</td></tr><tr><td>CD0789</td><td>Conserved hypothetical protein</td><td>0.3</td><td>AATGCAGGAAAA T C A C C C T G T T G A C G A A T T A A T A A G A</td><td>4.02</td></tr><tr><td>CD0838</td><td>DNA-binding protein</td><td>0.26</td><td>GTACAAGGATTTTAAAGATAAATATAGAAATTAATGTTT</td><td>4</td></tr><tr><td>CD0865</td><td>ADP-ribose binding protein</td><td>0.08</td><td>TTA ACAGGATTTATGTAGGTGTTTATAGAAATAAATAAAT</td><td>4.7</td></tr><tr><td>CD0999</td><td>ABC-type transport system, sulfonate/taurine</td><td>0.35</td><td>AAAAGAGGATTA TTC GATAAATTACAGAAATTAGTATAA</td><td>4.3</td></tr><tr><td>CD1149 minC</td><td>Cell division regulator</td><td>0.09</td><td>ATAAAGGTTTTAAAGCGTATTTGAAGAATATAAATGA</td><td>3.48</td></tr><tr><td>CD1214 spo0A#</td><td>Stage 0 sporulation protein A</td><td>0.2</td><td>TTAGGAGGAATATAATTTTGGAGTGTCGAATATGCTTTA</td><td>5.79</td></tr><tr><td>CD1221</td><td>Membrane protein</td><td>0.5</td><td>ATAAATGAATAACCCTTTTGTCCTAAAGCATACTATTTTATA</td><td>3.41</td></tr><tr><td>CD1264</td><td>Conserved hypothetical protein</td><td>0.08</td><td>AAAAAGAGGAAAA AT C A T T T T T AAT GTAGAATAATTGTACA</td><td>5.2</td></tr><tr><td>CD1317</td><td>Conserved hypothetical protein</td><td>0.42</td><td>TTAAAGGAAAGACAGGATAAATATAGAAATTTTATAACA</td><td>4.55</td></tr><tr><td>CD1484 ssuA</td><td>ABC-type transport system, alkanesulfonates family</td><td>0.08</td><td>ATAAAGGATTA A A C A T T T A G T G C G A A G T A G A T A T A T A</td><td>4.42</td></tr><tr><td>CD1498 sigA2</td><td>RNA polymerase sigma factor</td><td>0.05</td><td>ATAAAGGAGGAT T A T G C T G T T A G A A G T A A T A T A T T T A</td><td>4.74</td></tr><tr><td>CD1543.1</td><td>Conserved hypothetical protein</td><td>0.07</td><td>TATAAAGGAAAA ACCTCTTTTAATGTAGAAACTTATATTG</td><td>5.4</td></tr><tr><td>CD1622#</td><td>Conserved hypothetical protein</td><td>0.46</td><td>AAAAGAGGCTTGTTTCCATTTTGAGACAGTATAATTTTTATGT</td><td>2.09</td></tr><tr><td>CD1654 lplA</td><td>Lipoate-protein ligase</td><td>0.34</td><td>TTTAAAGGAGTTTTTATATTATTGTCGAATATTAAAAATT</td><td>5.43</td></tr><tr><td>CD1878</td><td>Sigma factor Tn1549-like</td><td>0.45</td><td>GATGGAGGAAACTTTATGGCAAAAGAGTATTATCTTTATATCA</td><td>2.73</td></tr><tr><td>CD1900</td><td>Conserved hypothetical protein</td><td>0.12</td><td>ATAAAGGAATTTACAAGTTATGTGTTGTATA GATATTGT</td><td>2.54</td></tr><tr><td>CD1935 spoVS</td><td>Stage V sporulation protein S</td><td>0.03</td><td>AATAAAGGTTTTCTTAAACGATTATAGAAGTATATTTTTA</td><td>4.24</td></tr><tr><td>CD1941</td><td>Conserved hypothetical protein</td><td>0.05</td><td>TAACAAGGAAACAGAACCCTCTC ATAGAAATTAATTATTA</td><td>3.68</td></tr><tr><td>CD1967</td><td>Conserved hypothetical protein</td><td>0.03</td><td>TTAGAAGGAAAA TAGCTTTTATC ATCAAATTATAAATT</td><td>5.08</td></tr><tr><td>CD2063</td><td>Conserved hypothetical protein</td><td>0.45</td><td>GATAGAGGATTTATAAGTGTTTAAAGGTGAATTAAATAAAT</td><td>3.48</td></tr><tr><td>CD2137 frr</td><td>Ribosome-recycling factor</td><td>0.41</td><td>AATAAAGGTATTTGAGCTTACAACAGAGAATATAATAAGA</td><td>3.62</td></tr><tr><td>CD2447</td><td>Putative histidine triad protein</td><td>0.06</td><td>GTAAGAAGGAATTTTGCTATAACATGTAGAAATTATAATATTTA</td><td>4.83</td></tr><tr><td>CD2646 ftsZ</td><td>Cell division protein FtsZ</td><td>0.34</td><td>TAAAAAGGAAAA TTTACGTTTTTGTTGAATATGTTACTT</td><td>4.57</td></tr><tr><td>CD2650</td><td>Cell division protein Fts-Q type</td><td>0.34</td><td>TTTGTAGGAAAAACAAGCTCTAAAGGTGTATTATTAACCA</td><td>2.66</td></tr><tr><td>CD2656 spoVD</td><td>Stage V sporulation protein D</td><td>0.15</td><td>TCTAAATGAATATAAAAAATAAAAAAAAGAATATAATTATAAAA</td><td>3.58</td></tr><tr><td>CD2657*</td><td>Conserved hypothetical protein</td><td>0.36</td><td>TTTAGAGGAAGA AATATGATTAATACAAATATAGTATA</td><td>1.82</td></tr><tr><td>CD2989 ssuA2</td><td>Sulfonate-family ABC-type transporter</td><td>0.39</td><td>AAAAGAGGATTAAGTATGGATGATGTGGAATTTGTTAATA</td><td>4.6</td></tr><tr><td>CD3221</td><td>Peptidase, M20D family</td><td>0.16</td><td>GAAACAGGGGAATTATTATTAATGGTGAATTTATTTAATA</td><td>2.85</td></tr><tr><td>CD3290</td><td>Conserved hypothetical protein</td><td>0.19</td><td>ACAAGAGGGATTGTTGATGATTTTATCGAAATCCTAATT</td><td>5.02</td></tr><tr><td>CD3317 fdhF</td><td>Formate dehydrogenase-H</td><td>0.31</td><td>TTAAGAGGAATTGTGAGAAAAATTGTTGAATTTATAAGATA</td><td>3.76</td></tr><tr><td>CD3458</td><td>Membrane protein</td><td>0.07</td><td>GTTAAAGGAAATTGTAGGTAGTTTATCGAATTTGTTAATCA</td><td>5.13</td></tr><tr><td>CD3516 spoVG</td><td>Regulator of spore cortex synthesis</td><td>0.01</td><td>AAAAGAGGATATCCCTAGTTGTTCCATAGAATTATTTA</td><td>4.89</td></tr><tr><td>CD3673</td><td>DNA-binding protein SpoOJ-like</td><td>0.21</td><td>TTCTCAGGAAATAATTAAGCTTACTGTAAAAAAACAA</td><td>2.16</td></tr></tbody></table>

<footnote># means that two promoters are mapped upstream of the gene. *promoters with a lower score not used to obtain the consensus, which was determined with the promoters with a score > 2. & as determined by transcriptome analysis (Saujet et al., 2011). The genes underlined have been previously mapped by RACE (Saujet et al., 2011). The TTSs mapped are underlined and the –10 and –35 boxes are indicated in bold.</footnote>
we can propose a probable existence of an overlapping SigF and SigG regulons.

## SigL-Dependent Promoters in *C. difficile*

SigL is a sigma 54 type sigma factor. In the genome of *C. difficile*, an important number of EBPs is present (Nie et al., 2019). Since EBPs activate SigL-dependent promoters, we can expect a large set of SigL-controlled genes in this enteropathogen. However, SigL and its regulon remain poorly characterized in *C. difficile* despite its possible important role in the control of metabolism. Among the TSS mapped by 5-end RNA-seq and identified as potential SigL targets in silico (Supplementary Table S2 and Figure 4), we first searched for the presence of the highly conserved sequence GGC at position -24 and GC at position -12 (Francke et al., 2011). We identified 13 promoters with these conserved motifs upstream of TSSs (Supplementary Table S4 and Table 4). The alignment of these "-24,-12" promoters led to the consensus TGGCA-N<sup>6</sup>-[A,T]TGCT[A,T] (Figure 5). We used this consensus to check for additional "-24,-12" promoters in the 200 pb upstream of start codons of all *C. difficile* genes. By this approach, we identified 30 putative SigL-dependent promoters corresponding to about 95 controlled genes (Table 4). As described in other bacteria (Danson et al., 2019; Nie et al., 2019), it is worth noting that 23 genes or operons with a "-24,-12" promoter out of 30 are adjacent on the chromosome to genes encoding a sigma 54-dependent activator (EBP), which is probably involved in their control.

## SigL-Mediated Control of Gene Expression in *C. difficile*

To identify genes expressed under the control of SigL, we performed a transcriptome analysis comparing the expression in the strain 630Δ*erm* and the *sigL::erm* mutant after 4 h of growth in TY. Approximately 7.5 % of the *C. difficile* genes were found to be differentially expressed between these two strains. 165 genes were up-regulated and 124 genes were down-regulated in the *sigL::erm* mutant (Supplementary Table S7). We confirmed these results by qRT-PCR analyses for 10 genes (Table 4). Only 27 genes identified as containing a "-24,-12" promoter and very likely controlled by SigL were down-regulated according to the transcriptomic data obtained. We were able to detect 4 additional operons (13 genes) as controlled by SigL using qRT-PCR (Table 4). These results are not surprising since most of the genes transcribed by RNAP associated to SigL respond to specific inducers through their associated EBP activators sensing these signals (Francke et al., 2011; Nie et al., 2019). For example, cysteine or proline strongly induces the expression of *cdsB* and *prdC* genes and the *prdA* operon, respectively, in a CdsR- or a PrdR-dependent manner (Bouillaut et al., 2013; Gu et al., 2017). We did not detect a SigL-dependent control of *cdsB* expression as previously observed (Gu et al., 2017) and only a two-fold decrease of expression was detected for three genes of the *prdA* operon in the *sigL::erm* mutant (Supplementary Table S7). For most of the genes very likely transcribed by SigL, the

inducer is probably absent in TY medium and we were unable to detect a control by this sigma factor on their basal level of expression.

## Role of SigL in the Physiology of *C. difficile*

To complete our view of the role of SigL in the physiology of *C. difficile*, we combined the data on the control of expression by SigL with those obtained on the identification of "-24,-12" promoters. As described in other firmicutes (Deutscher et al., 2006; Francke et al., 2011; Nie et al., 2019), we observed SigL-dependent promoters upstream of 7 operons encoding phosphotransferase systems (PTS) belonging to the mannose (*CD0284*, *CD0490*, *CD3279*), cellobiose (*CD0860*) or mannitol/galactitol (*CD0040*, *CD2283*, *CD2327*) sub-families (Table 4 and Supplementary Figure S9). These PTS operons are associated with 5 LevR-type activators (*CD0040*, *CD0283*, *CD2283*, *CD2328*, *CD3280*) and a sixth one, a pseudogen, *CD0858* (Stulke et al., 1998). Three operons encoding PTS (*CD0284*, *CD0490* and *CD2327*) were down-regulated in transcriptome (TY 4h of growth) or in qRT-PCR experiments (Table 4). These PTS are likely second-line systems used when glucose, fructose or mannitol are absent or depleted (Deutscher et al., 2006).

Interestingly, we also observed a large set of genes involved in amino-acid degradation or encoding peptidases and amino-acid permeases that are controlled either directly or indirectly by SigL and/or transcribed by the RNAP associated to SigL (Table 4 and Supplementary Table S7). *C. difficile* can use some amino acids as energy source through Stickland reactions (Neumann-Schaal et al., 2019). These reactions consist in the coupled fermentation of two amino acids in which one is oxidatively deaminated or decarboxylated and the other is reductively deaminated or reduced. The Stickland donors are valine, leucine, isoleucine and alanine while the acceptors are proline, glycine and leucine (Jackson et al., 2006). The *prdA* operon and the *prdC* gene, as well as the *hadA* operon required for the proline and the leucine reductive branch were transcribed by SigL and controlled by SigL-dependent EBPs, PrdR (Bouillaut et al., 2013) and probably LeuR (*CD0402*), respectively (Table 4 and Figure 7) (Nie et al., 2019). The expression of the *hadA* operon was strongly down-regulated in the *sigL::erm* mutant while the *prdA* operon was only weakly controlled by SigL in TY medium (Table 4). By contrast, the expression of the *grd* operon involved in glycine reduction was strongly up-regulated in the *sigL::erm* mutant compared to the wild-type strain (8 to 60-fold) (Supplementary Table S7). In *C. difficile*, SigL likely plays a key role in the hierarchy of amino acid utilization through Stickland reactions favoring the reductive degradation of proline and leucine as an energy source in detriment of glycine. A "-24,-12" promoter recognized by SigL was also mapped upstream of the *ord* operon involved in the oxidative degradation of ornithine (Fonknechten et al., 2009) associated with the EBP-OrdR. In addition, SigL seemed to control directly or indirectly the expression of genes involved in the uptake and degradation of aromatic amino acids (tyrosine
<table><caption>TABLE 4 | SigL-dependent promoters in C. difficile 630 based on TSS mapping, transcriptome data and in silico analysis.</caption><thead><tr><th>Gene</th><th>Operons</th><th>Functions</th><th>Fold change<br>sigL::erm/ 630Δerm</th><th>Promoter -24,-12</th><th>SigL-dependent<br>associated regulator</th></tr></thead><tbody><tr><td>CD0040</td><td>CD0040-CD0043</td><td>Activator, PTS Galactitol family</td><td></td><td>ATGGCATATAAGTTGCTAT</td><td>CD0040, LevR-type</td></tr><tr><td>CD0166</td><td>CD0166-CD0165</td><td>Peptidase, amino acid transporter</td><td>0.17*</td><td>ATGGCAATAATTTGCTTA</td><td>CD0167, GamR-type</td></tr><tr><td>CD0284</td><td>CD0284-CD0289</td><td>PTS Mannose/fructose/sorbose family</td><td>0.2 to 0.3</td><td>TTGGCACGGCAATTGCTTA</td><td>CD0283, LevR-type</td></tr><tr><td>CD0395/ hadA#</td><td>hadAIBC-acdB-etfBA1</td><td>Leucine utilization</td><td>&lt; 0.001 (0.00001*)</td><td>TTGGCACGATTTATGCTTT</td><td>CD0402/LeuR</td></tr><tr><td>CD0442/ ord#</td><td>ord-ortAB-oraSEF-orr-nhaC</td><td>Ornithine degradation</td><td>NR</td><td>TTGGCACGATTTATGCTTT</td><td>CD0441/OrdR</td></tr><tr><td>CD0490#</td><td>CD0490-CD0494</td><td>Sugar-P-dehydrogenase, PTS mannose/fructose/sorbose family</td><td>0.24*</td><td>TTGGCATGAAAGTTGCTTT</td><td>CD0516? LevR-type</td></tr><tr><td>CD0800/ crt1#</td><td>crt1-CD0801-catB-bcd-etfBA2</td><td>Crotonase, permease, CoA transferase, acyl-CoA dehydrogenase, EtfBA</td><td>NR</td><td>TTGGCATAGTACTTTGCTAT</td><td>CD0806/YctR</td></tr><tr><td>CD0860</td><td>CD0860-CD0863 -malH1</td><td>PTS lactose/cellobiose family- Maltose-6"P glucosidase</td><td></td><td>CTGGCAATAACTTTGCTTA</td><td>CD0858°</td></tr><tr><td>CD1187#</td><td>CD1187-CD1189</td><td>CHP, γ-glutamyl-γ-aminobutyrate hydrolase, amino acid permease</td><td>NR</td><td>TTGGCATACATATTGCTAA</td><td>CD1186 GamR-type</td></tr><tr><td>CD1413/ rhaT#</td><td>CD1413</td><td>Membrane protein (RhaT)</td><td>0.31 (0.19*)</td><td>ATGGCATAGTTTTGCTTA</td><td>CD1412/XhaQ</td></tr><tr><td>CD1555#</td><td>CD1555</td><td>Putative serine/threonine exchanger</td><td>0.36 (0.5*)</td><td>TTGGCATATAATATGCTTA</td><td>?</td></tr><tr><td>CD1740#</td><td>CD1740-CD1741</td><td>Sarcosine reductase complex</td><td>NR</td><td>TTGGCATAGAAAATGCTTT</td><td>CD1739/SarR, TCS SarRS</td></tr><tr><td>CD2091 #</td><td>CD2091-CD2089</td><td>Putative Xanthine/uracile permease (PbuX), adenosine derivate deaminase</td><td>NR</td><td>TTGGCATATTATAATTGCTTC</td><td>CD2092-DioR1</td></tr><tr><td>CD2279</td><td>CD2279</td><td>Sugar-P dehydrogenase</td><td></td><td>ATGGCATAGATATTGCTAT</td><td>CD2283?</td></tr><tr><td>CD2283</td><td>CD2283-CD2280</td><td>Activator, PTS fructose/mannitol family</td><td>NR</td><td>ATGGCATGATAGTTGCTTA</td><td>CD2283 LevR-type</td></tr><tr><td>CD2327#</td><td>CD2327-CD2323</td><td>Arabitol/xylitol PTS, sugar-P dehydrogenases</td><td>0.03*</td><td>TTGGCACACAACATGCTTT</td><td>CD2328 LevR-type</td></tr><tr><td>CD2382#</td><td>CD2382-iorAB-butK</td><td>Aromatic aminotransferase, Indole-pyruvate oxidoreductase, butyrate kinase</td><td>0.2 to 0.02 (0.01*)</td><td>TTGGCATAGTAAATTGCTTA</td><td>CD2383/ZypR</td></tr><tr><td>CD2699#</td><td>CD2699-CD2697</td><td>Membrane proteins, peptidase</td><td>0.03 to 0.09 (0.01*)</td><td>TTGGCATAGTTTTGCTTA</td><td>CD2700</td></tr><tr><td>CD2733</td><td>CD2733-CD2734.1</td><td>PLP-dependent transferase, Na+/H+ antiporter, membrane protein</td><td></td><td>TTGGCACGTTGTTTGCTTA</td><td>CD2732</td></tr><tr><td>CD2862#</td><td>CD2862-CD2860</td><td>dipeptidase, membrane proteins</td><td></td><td>TTGGCATCAATATTGCTAC</td><td>CD2863</td></tr><tr><td>CD2870/ kdgT1</td><td>kdgT1-uxaA°</td><td>2-keto 3-deoxygluconate permease, altronate dehydratase</td><td>0.12 to 0.23</td><td>TTGGCATAGTAAATTGCTTT</td><td>CD2869/XduR</td></tr><tr><td>CD3093#</td><td>CD3093</td><td>γ-glutamyl-γ-aminobutyrate hydrolase</td><td></td><td>TTGGTATGCTACTTTGCTCT</td><td>CD3094 GamR-type</td></tr><tr><td>CD3094</td><td>CD3094</td><td>Sigma-54 dependent regulator</td><td></td><td>TTGGCAATTTTTGCTTT</td><td>CD3094 GamR-type</td></tr><tr><td>CD3184/ dpaL2</td><td>dpaL2</td><td>Diaminopropionate ammonia lyase</td><td>0.3*</td><td>TTGGCACGGTAATTGCTTT</td><td>CD3186/DioR2</td></tr><tr><td>CD3187</td><td>tdcF</td><td>Putative regulatory endoribonuclease</td><td></td><td>TTGGCACGTTAATTGCTTT</td><td>CD3186/DioR2</td></tr><tr><td>CD2085</td><td>dpaL1</td><td>Diaminopropionate ammonia lyase</td><td></td><td>TTGGCATGTTAATTGCTTA</td><td>CD3186/DioR2</td></tr><tr><td>CD3232/ cdsB #</td><td>cdsB</td><td>Cysteine desulfidase</td><td>NR°</td><td>ATGGCATGTATTTTGCTAT</td><td>CD3233/CdsR°</td></tr><tr><td>CD3244/ prdA</td><td>prdA-CD3243-prdBDE-prdE2F<br>CD3236</td><td>Proline utilization</td><td>0.4 to 0.5</td><td>TTGGCATAGGAATTGCTTA</td><td>CD3245/PrdR</td></tr><tr><td>CD3247/ prdC</td><td>prdC</td><td>Proline utilization</td><td></td><td>TTGGCATAGAAAATTGCTTT</td><td>CD3245/PrdR</td></tr><tr><td>CD3279</td><td>CD3279-CD3275</td><td>PTS Mannose/fructose/sorbose family, sugar-P isomerase</td><td></td><td>TTGGCATACTTTTTGCTTT</td><td>CD3280 LevR-type</td></tr></tbody></table>

<footnote># the transcriptional start sites were mapped by 5'-end RNA-seq (see Supplementary Table S2). * as determined by qRT-PCR. °pseudogen in strain 630. % cdsB is controlled by SigL in TY in the presence of cysteine (Gu et al., 2018). The CdsR regulator of C. difficile is of type 2 (Nie et al., 2019). NR = non regulated by SigL in qRT-PCR in the conditions tested. For most of them, the transcriptional regulators and genes were renamed as proposed by Nie et al. (2019) or by Dannheim et al. (2017).</footnote>
<figure><img src="image_9.png" /><figcaption>**FIGURE 6** | Promoters controlled by SigH in *C. difficile*. Examples of 5'-end RNA-seq (TAP-/TAP+ profile comparison) and RNA-seq data for dual tandem TSSs and/or internal TSSs inside coding sequences corresponding to SigH-dependent promoters are shown. Panel (A): *CD1482-CD1484* operon and panel (B): *CD0788-CD0789* operon.</figcaption></figure>

and phenylalanine) (Bradshaw et al., 2019). SigL was also found to control the expression of genes involved in the degradation of serine (sdaB) and cysteine (cdsB and maybe also cysK) and encoding a probable serine/threonine exchanger (CD1555) or a transporter of cysteine/cystine (CD2174 operon) (Dubois et al., 2016; Gu et al., 2017). Finally, a "-24,-12" promoter was mapped upstream of the CD1740-CD1741 operon encoding a sarcosine reductase associated with a two-component system, SarR-SarS. In conclusion, SigL plays a crucial role in the control of peptide and amino acid catabolism.

## Impact of SigL Inactivation on Growth and End-Fermentation Products

To confirm the role of SigL in *C. difficile*, we tested the growth of the WT and *sigL::erm* mutant strains (Dubois et al., 2016) in TY medium (Figure 8A). The mutant inactivated for SigL showed a reduced growth rate as compared to the strain 630Δ*erm* (Figure 8A). This result is in agreement with the proposed key role of SigL in the degradation of amino acids that are probably used as energy sources by *C. difficile* in TY. We then tested the effect on growth

of the addition of glucose or glycine since the *grd* operon encoding the glycine reductase is not transcribed from a "-24,-12" promoter. Interestingly, the addition of glucose to TY partially restored the growth of the *sigL::erm* mutant (Figure 8A) but this was not the case for the addition of glycine (data not shown). These results suggested that the addition of glucose rerouting the metabolism (Antunes et al., 2012) allows to compensate the reduced utilization of peptides while glycine is not sufficient to restore growth even if the *grd* operon is strongly up-regulated (8- to 60-fold) in the *sigL::erm* mutant (Supplementary Table S7).

To confirm the impact of SigL inactivation on the fermentation processes, we also analyzed the end products of fermentation in the 630Δ*erm* and of the *sigL::erm* mutant by gas-liquid chromatography after 48 h of growth in TY medium. While the total concentration of volatile acids slightly increased, we observed a drastic decrease of non-volatile acids in the *sigL::erm* mutant compared to the 630Δ*erm* strain (Figures 8B,C). These results confirmed that SigL inactivation led to important changes in metabolism and metabolite production. We observed a drastic depletion of pyruvic acid (Figure 8C) that can be associated at least partly to a drop-in
<figure><img src="image_10.png" /></figure>

expression of the genes involved in cysteine and serine uptake and degradation (Figure 7) and to an increased expression of *ldh* and maybe also *adhE* (Supplementary Figure S9). In the *sigL* mutant, the strong decreased production of the end-product of leucine degradation, isocaproic acid, (Figure 8B) correlates with the drastic reduction of expression of the *hadA* operon (Figure 7). Decreased leucine degradation through the reductive Stickland reactions could redirect leucine catabolism to the oxidative Stickland reactions explaining the significant increase of isovaleric acid in the supernatant of the *sigL::erm* mutant (Figure 8B) (Neumann-Schaal et al., 2015). The apparent depletion of succinate might result from the strong up-regulation of the pathway of conversion of succinate to crotonyl-CoA (Supplementary Figure S9). Globally, the effects observed on metabolites could also result from a change in the metabolism balance because of the lack of SigL, to an indirect influence of SigL on the fermentation pathways through the modulation of intracellular concentration of metabolites or a cascade of regulation by uncharacterized mechanisms. It is worth noting that a large set of transcriptional regulators are up- or down-regulated in the *sigL::erm* mutant (Supplementary Table S7) allowing indirect metabolic controls.

CONCLUSION

In the present study, we provide the first transcriptional map of the *C. difficile* genome demonstrating a complex structure of transcriptional units and operon organization in

this pathogen. We have applied the combination of *in silico* and experimental strategies to establish the list of proposed TSS that will serve as a start point for further studies in global and gene-specific scale. In addition to primary TSSs, this genome-wide TSS mapping revealed the presence of tandem and internal promoters suggesting alternative ways to accommodate gene expression changes during *C. difficile* development. This pathogen uses its large arsenal of sigma factors to determine the promoter selectivity. *C. difficile* has also two housekeeping SigA encoding genes, one transcribed from a SigH-dependent promoter (Saujet et al., 2011) and the second one likely in operon with *dnaG* and an unexpected number of extended –10 boxes. This work extends our knowledge on the expression program during stationary phase (mediated by stationary phase sigma factors, SigH and SigB) and sporulation. The promiscuity for the sequence recognition for the SigF and SigG sigma factors suggests the existence of a less controlled forespore program compared to *B. subtilis* as previously proposed (Saujet et al., 2013, 2014; Fimlaid et al., 2013). Finally, we identify several regulators likely expressed under the control of sporulation specific sigma factors (CD0629 and CD2316 by SigE, CD2599 by SigF or SigG) in addition to conserved SpoIIID and SpoVT transcriptional regulators (Saujet et al., 2014) suggesting that additional uncharacterized mechanisms of sporulation regulation might exist in *C. difficile*. Focusing on one of important alternative sigma factors, we characterized here the SigL regulon largely involved in amino-acid degradation pathways as a key trait for the colonization of the gut. In addition to determination of TSS, valuable information
<figure><img src="image_11.png" /><figcaption>**FIGURE 8** | Growth and fermentation products of the WT and *sigL::erm* mutant strains. (A) Growth curves of the 630Δ*erm* and the *sigL::erm* strains in TY and TY + 0.5% glucose. Growth curves are representative of at least three independent experiments. Concentrations of volatile (B) and non-volatile (C) fermentations end products of the 630Δ*erm* and the *sigL::erm* strains grown 48 h in TY. Gas-liquid chromatography analysis was performed from the supernatant of the culture of both strains as described in materials and methods. The concentration of the fermentation end products was standardized on the OD₆₀₀nm of the cultures after 48 h of growth. We could not detect butanol nor ethanol by these assays. The error bars represent the standard deviation of the mean. Asterisks indicate statistical significance (t-test or t-test Welch, ***p* ≤ 0.001, **p* ≤ 0.01, *p* ≤ 0.05).</figcaption></figure>

could be deduced from 5'-end RNA-seq on the potential
RNA processing events impacting the RNA stability, as
well as on the 5' untranslated regions harboring important
regulatory motifs for transcriptional and post-transcriptional
regulatory mechanisms. These data complete out previous
genome-wide identification of ncRNAs in *C. difficile* linking
together different facets of global regulatory networks in this
pathogen. This study thus constitutes the first essential step
towards better understanding of the complex transcriptional

and post-transcriptional regulations governing the infection
cycle and adaptation of this successful enteropathogen to
host environments.

DATA AVAILABILITY STATEMENT

The datasets presented in this study can be found in
online repositories. The names of the repository/repositories
and accession number(s) can be found in the article/
**Supplementary Material.**

## AUTHOR CONTRIBUTIONS

OS and IM-V co-designed the study, analyzed the data and wrote the manuscript. OS, TD, and LS performed the experiments. MM, PS, and MG collected and analyzed the data in *silico*. OS, IM-V, LS, PB, and BD manually analyzed the TSS. TD collected and analyzed the transcriptomic data. BD supervised the experimental work. TD, MM, MG, and BD revised the manuscript. All authors contributed to the article and approved the submitted version.

## FUNDING

This work was supported by the Institut Pasteur, the Université de Paris, the Université Paris-Saclay, the Institute for Integrative Biology of the Cell, the Institut Universitaire de France (to OS and to IM-V), the Agence Nationale de la Recherche (“CloSTARn”, ANR-13-JSV3-0005-01 to OS), the DIM-1HEALTH regional Ile-de-France program (LSP Grant No. 164466), the CNRS-RFBR PRC 2019 (Grant No. 288426 N° 19-54-15003) to OS. The computational analysis was supported by a grant from the Russian Science Foundation (18-14-00358).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01939/full#supplementary-material

**FIGURE S1** | Representative example of TSS identification by PhageTerm software. Cov2HTML (Monot et al., 2014) visualization of 5'-end RNA-seq data for TAP-/TAP+ profile comparison is presented in the left panel and PhageTerm (Garneau et al., 2017). TSS identification is shown in the right panel.

**FIGURE S2** | Operon organization of the genes coding the proteins involved in flagella biosynthesis and function associated with TSS mapping. The flagella genes are organized in several transcriptional units driven by RNAP in complex with alternative flagellum SigD or housekeeping SigA sigma factors. The TSS are shown by broken arrows, SigD-associated promoters are indicated in blue, SigA-associated promoters are shown in green and undefined promoters are depicted in black. Three boxed images show a zoom to the TSS mapping visualization for selected promoters. The genomic region visualization was made with MAGE platform (Vallenet et al., 2020), coding sequences are indicated by ORF in red and previously identified noncoding RNAs are shown by blue-green boxes.

**FIGURE S3** | Nucleotide composition of TSS and frequency of each initiation nucleotide in *C. difficile*. 100 nucleotide regions (from -99 to +1) upstream of the TSS for *C. difficile* strain 630 genes have been aligned and the deduced consensus sequence generated with Weblogo program is presented.

+1" indicates TSS position, "-10" indicates the position of -10 promoter element.

**FIGURE S4** | Length distribution of 5'UTR regions in the *C. difficile* strain 630 genome. The graph shows the length of 5'UTR (distance from the TSS to the translational start) and the number of genes with corresponding 5'UTR length range.

**FIGURE S5** | TSS mapping of dual (A) and internal (B) promoters. Representative examples of 5'-end RNA-seq (TAP-/TAP+ profile comparison) and RNA-seq data for dual tandem TSS and internal TSS inside the coding sequences are shown in panels (A,B), respectively. Cov2HTML (Monot et al., 2014) was used for the visualization. On a RNA-seq and 5'-end RNA-seq sequence read mapping visualization, coding sequences are indicated by blue arrows. The 5'-end RNA-seq data for either positive "strand +" or negative "strand -" strands are presented in the panels. The TSS identified by 5'-end RNA-seq are indicated by red broken arrows and potential processing sites are indicated by scissors mark. Sigma factor consensus associated with a given TSS is indicated. The TSS corresponds to a position with significantly greater number of reads in TAP+ sample, potential cleavage site corresponds to position with large number of reads in both TAP- and TAP+ samples. 5'-end RNA-seq data show 51-bp reads matching to the 5'-transcript ends, while RNA-seq data show reads covering whole transcript.

**FIGURE S6** | Examples of cleavage sites detected by TSS-mapping. Representative examples of 5'-end RNA-seq (TAP-/TAP+ profile comparison) and RNA-seq data for complex potential processing profiles are shown. The RNA-seq and 5'-end RNA-seq data visualization is presented as in **Figure 1**. Potential cleavage site shown by scissors mark corresponds to a position with large number of reads in both TAP- and TAP+ samples. RBS are shown by green boxes to highlight the positions of potential cleavage sites in the proximity or inside RBS.

**FIGURE S7** | Additional examples of promoters controlling regulatory genes detected by TSS-mapping. Representative examples of 5'-end RNA-seq (TAP-/TAP+ profile comparison) data for the identification of TSS for genes encoding important transcriptional regulators are shown. The 5'-end RNA-seq data visualization is presented as in **Figure 1**. The sequence of promoter region is shown upstream of TSS with the -35 and -10 promoter elements indicated in blue and TSS indicated in red.

**FIGURE S8** | Confusion matrix on training with the standard and strict criterion. Confusion matrices for the training set. (A) Listing scheme as in **Figure 4**, (B) Strict listing scheme with one top-scoring sigma factor taken as the prediction.

**FIGURE S9** | Carbon metabolism genes controlled by SigL in *C. difficile*. Upstream of genes indicated in green, a "-24, -12" promoter was mapped or identified in *silico* (Table 4). * indicated genes with a "-24, -12" promoter (Table 4). Genes in green are positively controlled by SigL in transcriptome or in qRT-PCR experiments as indicated by #. Genes indicated in red are negatively controlled by SigL in transcriptome (Supplementary Table S7). Green arrow indicates a compound less detected by Gas-liquid chromatography analysis (Figure 8).

**TABLE S1** | Strains and plasmids used in this study.

**TABLE S2** | Complete list of TSS and predicted promoters for *C. difficile* 630 CDS genes.

**TABLE S3** | Complete list of TSS and predicted promoters for *C. difficile* 630 ncRNA genes.

**TABLE S4** | *In silico* association of TSS with known sigma factors in *C. difficile*.

**TABLE S5** | Identification of promoters recognized by SigF and SigG.

**TABLE S6** | Identification of promoters recognized by SigE and/or SigK.

**TABLE S7** | Transcriptomics results for *sigL::erm* mutant.
REFERENCES

Abt, M. C., Mckenney, P. T., and Pamer, E. G. (2016). *Clostridium difficile* colitis: pathogenesis and host defence. *Nat. Rev. Microbiol.* 14, 609–620. doi: 10.1038/nrmicro.2016.108

Amman, F., Wolfinger, M. T., Lorenz, R., Hofacker, I. L., Stadler, P. F., and Findeiss, S. (2014). TSSAR: TSS annotation regime for dRNA-seq data. *BMC Bioinform.* 15:89. doi: 10.1186/1471-2105-15-89

Andre, G., Even, S., Putzer, H., Burguiere, P., Croux, C., Danchin, A., et al. (2008). S-box and T-box riboswitches and antisense RNA control a sulfur metabolic operon of *Clostridium acetobutylicum*. *Nucleic Acids Res.* 36, 5955–5969. doi: 10.1093/nar/gkn601

Anjuwon-Foster, B. R., and Tamayo, R. (2017). A genetic switch controls the production of flagella and toxins in *Clostridium difficile*. *PLoS Genet.* 13:e1006701. doi: 10.1371/journal.pgen.1006701

Antunes, A., Camiade, E., Monot, M., Courtois, E., Barbut, F., Sernova, N. V., et al. (2012). Global transcriptional control by glucose and carbon regulator CcpA in *Clostridium difficile*. *Nucleic Acids Res.* 40, 10701–10718. doi: 10.1093/nar/gks864

Aubry, A., Hussack, G., Chen, W., Kuolee, R., Twine, S. M., Fulton, K. M., et al. (2012). Modulation of toxin production by the flagellar regulon in *Clostridium difficile*. *Infect. Immun.* 80, 3521–3532. doi: 10.1128/iai.00224-12

Babski, J., Haas, K. A., Nather-Schindler, D., Pfeiffer, F., Forstner, K. U., Hammelmann, M., et al. (2016). Genome-wide identification of transcriptional start sites in the haloarchaeon *Haloferax volcanii* based on differential RNA-Seq (dRNA-Seq). *BMC Genom.* 17:629. doi: 10.1186/1471-2105-15-629

Banawas, S. S. (2018). *Clostridium difficile* infections: a global overview of drug sensitivity and resistance mechanisms. *Biomed. Res. Int.* 2018:8414257.

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. *J. R. Stat. Soc. Ser. B* 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

Bouillaut, L., Self, W. T., and Sonenshein, A. L. (2013). Proline-dependent regulation of *Clostridium difficile* Stickland metabolism. *J. Bacteriol.* 195, 844–854. doi: 10.1128/jb.01492-12

Bradshaw, W. J., Bruxelle, J. F., Kovacs-Simon, A., Harmer, N. J., Janoir, C., Pechine, S., et al. (2019). Molecular features of lipoprotein CD0873: a potential vaccine against the human pathogen *Clostridioides difficile*. *J. Biol. Chem.* 294, 15850–15861. doi: 10.1074/jbc.ra119.010120

Braun, V., Hundsberger, T., Leukel, P., Sauerborn, M., and Von Eichel-Streiber, C. (1996). Definition of the single integration site of the pathogenicity locus in *Clostridium difficile*. *Gene* 181, 29–38. doi: 10.1016/s0378-1119(96)00398-8

Brouwer, M. S., Warburton, P. J., Roberts, A. P., Mullany, P., and Allan, E. (2011). Genetic organisation, mobility and predicted functions of genes on integrated, mobile genetic elements in sequenced strains of *Clostridium difficile*. *PLoS One* 6:e23014. doi: 10.1371/journal.pgen.0023014

Browning, D. F., and Busby, S. J. (2016). Local and global regulation of transcription initiation in bacteria. *Nat. Rev. Microbiol.* 14, 638–650. doi: 10.1038/nrmicro.2016.103

Burgess, R. R., and Anthony, L. (2001). How sigma docks to RNA polymerase and what sigma does. *Curr. Opin. Microbiol.* 4, 126–131. doi: 10.1016/s1369-5274(00)00177-6

Carlier, J. P., and Sellier, N. (1989). Gas chromatographic-mass spectral studies after methylation of metabolites produced by some anaerobic bacteria in spent media. *J. Chromatogr.* 493, 257–273. doi: 10.1016/S0378-4347(00)82733-4

Carroll, K. C., and Bartlett, J. G. (2011). Biology of *Clostridium difficile*: implications for epidemiology and diagnosis. *Annu. Rev. Microbiol.* 65, 501–521.

Chao, Y., Papenfort, K., Reinhardt, R., Sharma, C. M., and Vogel, J. (2012). An atlas of Hfq-bound transcripts reveals 3' UTRs as a genomic reservoir of regulatory small RNAs. *EMBO J.* 31, 4005–4019. doi: 10.1038/emboj.2012.229

Cuklina, J., Hahn, J., Imakaev, M., Omasits, U., Forstner, K. U., Ljubimov, N., et al. (2016). Genome-wide transcription start site mapping of *Bradyrhizobium japonicum* grown free-living or in symbiosis - a rich resource to identify new transcripts, proteins and to study gene regulation. *BMC Genomics* 17:302. doi: 10.1186/1471-2105-15-302

Dannheim, H., Riedel, T., Neumann-Schaal, M., Bunk, B., Schober, I., Sproer, C., et al. (2017). Manual curation and reannotation of the genomes of *Clostridium difficile* 630Deltaerm and *C. difficile* 630. *J. Med. Microbiol.* 66, 286–293. doi: 10.1099/jmm.0.000427

Danson, A. E., Jovanovic, M., Buck, M., and Zhang, X. (2019). Mechanisms of sigma(54)-dependent transcription initiation and regulation. *J. Mol. Biol.* 431, 3960–3974. doi: doi: 10.1016/j.jmb.2019.04.022

DeLoughery, A., Lalanne, J. B., Losick, R., and Li, G. W. (2018). Maturation of polycistronic mRNAs by the endoribonuclease RNase Y and its associated Y-complex in *Bacillus subtilis*. *Proc. Natl. Acad. Sci. U.S.A.* 115, E5585–E5594.

Deutscher, J., Francke, C., and Postma, P. W. (2006). How phosphotransferase system-related protein phosphorylation regulates carbohydrate metabolism in bacteria. *Microbiol. Mol. Biol. Rev.* 70, 939–1031. doi: 10.1128/mmbr.00204-06

Dineen, S. S., Mcbride, S. M., and Sonenshein, A. L. (2010). Integration of metabolism and virulence by *Clostridium difficile* CodY. *J. Bacteriol.* 192, 5350–5362. doi: 10.1128/jb.00341-10

Dubois, T., Dancer-Thibonnier, M., Monot, M., Hamiot, A., Bouillaut, L., Soutourina, O., et al. (2016). Control of *Clostridium difficile* physiopathology in response to cysteine availability. *Infect. Immun.* 84, 2389–2405. doi: 10.1128/iai.00121-16

Dupuy, B., and Sonenshein, A. L. (1998). Regulated transcription of *Clostridium difficile* toxin genes. *Mol. Microbiol.* 27, 107–120.

Eckweiler, D., Dudek, C. A., Hartlich, J., Brotje, D., and Jahn, D. (2018). PRODORIC2: the bacterial gene regulation database in 2018. *Nucleic Acids Res.* 46, D320–D326.

El Meouche, I., Peltier, J., Monot, M., Soutourina, O., Pestel-Caron, M., Dupuy, B., et al. (2013). Characterization of the SigD regulon of *C. difficile* and its positive control of toxin production through the regulation of tcdR. *PLoS One* 8:e83748. doi: 10.1371/journal.pgen.0083748

Emerson, J. E., Reynolds, C. B., Fagan, R. P., Shaw, H. A., Goulding, D., and Fairweather, N. F. (2009). A novel genetic switch controls phase variable expression of CwpV, a *Clostridium difficile* cell wall protein. *Mol. Microbiol.* 74, 541–556. doi: 10.1111/j.1365-2958.2009.06812.x

Fimlaid, K. A., Bond, J. P., Schutz, K. C., Putnam, E. E., Leung, J. M., Lawley, T. D., et al. (2013). Global analysis of the sporulation pathway of *Clostridium difficile*. *PLoS Genet.* 9:e1003660. doi: 10.1371/journal.pgen.1003660

Fonknechten, N., Perret, A., Perchat, N., Tricot, S., Lechaplais, C., Vallenet, D., et al. (2009). A conserved gene cluster rules anaerobic oxidative degradation of L-ornithine. *J. Bacteriol.* 191, 3162–3167. doi: 10.1128/jb.01777-08

Francke, C., Groot Kormelink, T., Hagemeijer, Y., Overmars, L., Sluijter, V., Moezelaar, R., et al. (2011). Comparative analyses imply that the enigmatic Sigma factor 54 is a central controller of the bacterial exterior. *BMC Genomics* 12:385. doi: 10.1186/1471-2105-15-385

Garneau, J. R., Depardieu, F., Fortier, L. C., Bikard, D., and Monot, M. (2017). PhageTerm: a tool for fast and accurate determination of phage termini and packaging mechanism using next-generation sequencing data. *Sci. Rep.* 7:8292.

Gruber, T. M., and Gross, C. A. (2003). Multiple sigma subunits and the partitioning of bacterial transcription space. *Annu. Rev. Microbiol.* 57, 441–466. doi: 10.1146/annurev.micro.57.030502.090913

Gu, H., Shi, K., Liao, Z., Qi, H., Chen, S., Wang, H., et al. (2018). Time-resolved transcriptome analysis of *Clostridium difficile* R20291 response to cysteine. *Microbiol. Res.* 215, 114–125. doi: 10.1016/j.micres.2018.07.003

Gu, H., Yang, Y., Wang, M., Chen, S., Wang, H., Li, S., et al. (2017). Novel cysteine desulfidase CdsB involved in releasing cysteine repression of toxin synthesis in *Clostridium difficile*. *Front. Cell Infect. Microbiol.* 7:531. doi: 10.3389/fcimb.2017.00531

Guo, M. S., Updegrove, T. B., Gogol, E. B., Shabalina, S. A., Gross, C. A., and Storz, G. (2014). MicL, a new sigmaE-dependent sRNA, combats envelope stress by repressing synthesis of Lpp, the major outer membrane lipoprotein. *Genes Dev.* 28, 1620–1634. doi: 10.1101/gad.243485.114

Harden, T. T., Wells, C. D., Friedman, L. J., Landick, R., Hochschild, A., Kondev, J., et al. (2016). Bacterial RNA polymerase can retain sigma70 throughout transcription. *Proc. Natl. Acad. Sci. U.S.A.* 113, 602–607. doi: 10.1073/pnas.1513899113

Ho, T. D., and Ellermeier, C. D. (2011). PrsW is required for colonization, resistance to antimicrobial peptides, and expression of extracytoplasmic
function sigma factors in *Clostridium difficile*. *Infect. Immun.* 79, 3229–3238. doi: 10.1128/iai.000019-11

Hor, J., Gorski, S. A., and Vogel, J. (2018). Bacterial RNA biology on a genome scale. *Mol. Cell.* 70, 785–799. doi: 10.1016/j.molcel.2017.12.023

Jackson, S., Calos, M., Myers, A., and Self, W. T. (2006). Analysis of proline reduction in the nosocomial pathogen *Clostridium difficile*. *J. Bacteriol.* 188, 8487–8495. doi: 10.1128/jb.01370-06

Jager, D., Forstner, K. U., Sharma, C. M., Santangelo, T. J., and Reeve, J. N. (2014). Primary transcriptome map of the hyperthermophilic archaeon *Thermococcus kodakarensis*. *BMC Genom.* 15:684. doi: 10.1186/1471-2105-15-684

Jager, D., Sharma, C. M., Thomsen, J., Ehlers, C., Vogel, J., and Schmitz, R. A. (2009). Deep sequencing analysis of the methanosarcina mazei Go1 transcriptome in response to nitrogen availability. *Proc. Natl. Acad. Sci. U.S.A.* 106, 21878–21882. doi: 10.1073/pnas.0909051106

Janoir, C. (2016). Virulence factors of *Clostridium difficile* and their role during infection. *Anaerobe* 37, 13–24. doi: 10.1016/j.anaerobe.2015.10.009

Jorjani, H., and Zavolan, M. (2014). TSSer: an automated method to identify transcription start sites in prokaryotic genomes from differential RNA sequencing data. *Bioinformatics* 30, 971–974. doi: 10.1093/bioinformatics/btt752

Keessen, E. C., Gaastra, W., and Lipman, L. J. (2011). *Clostridium difficile* infection in humans and animals, differences and similarities. *Vet. Microbiol.* 153, 205–217. doi: 10.1016/j.vetmic.2011.03.020

Keilty, S., and Rosenberg, M. (1987). Constitutive function of a positively regulated promoter reveals new sequences essential for activity. *J. Biol. Chem.* 262, 6389–6395.

Kint, N., Janoir, C., Monot, M., Hoys, S., Soutourina, O., Dupuy, B., et al. (2017). The alternative sigma factor sigmaB plays a crucial role in adaptive strategies of *Clostridium difficile* during gut infection. *Environ. Microbiol.* 19, 1933–1958. doi: 10.1111/1462-2920.13696

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. *Methods* 25, 402–408. doi: 10.1006/meth.2001.1262

Lo, A. W., Moriel, D. G., Phan, M. D., Schulz, B. L., Kidd, T. J., Beatson, S. A., et al. (2017). 'Omic' approaches to study uropathogenic *Escherichia coli* virulence. *Trends Microbiol.* 25, 729–740. doi: 10.1016/j.tim.2017.04.006

Maikova, A., Peltier, J., Boudry, P., Hajnsdorf, E., Kint, N., Monot, M., et al. (2018). Discovery of new type I toxin-antitoxin systems adjacent to CRISPR arrays in *Clostridium difficile*. *Nucleic Acids Res.* 46, 4733–4751. doi: 10.1093/nar/gky124

Mani, N., and Dupuy, B. (2001). Regulation of toxin synthesis in *Clostridium difficile* by an alternative RNA polymerase sigma factor. *Proc. Natl. Acad. Sci. U.S.A.* 98, 5844–5849. doi: 10.1073/pnas.101126598

Mani, N., Lyras, D., Barroso, L., Howarth, P., Wilkins, T., Rood, J. I., et al. (2002). Environmental response and autoregulation of *Clostridium difficile* TxeR, a sigma factor for toxin gene expression. *J. Bacteriol.* 184, 5971–5978. doi: 10.1128/jb.184.21.5971-5978.2002

Martin-Verstraete, I., Peltier, J., and Dupuy, B. (2016). The regulatory networks that control *Clostridium difficile* toxin synthesis. *Toxins* 8:153. doi: 10.3390/toxins8050153

Matamouros, S., England, P., and Dupuy, B. (2007). *Clostridium difficile* toxin expression is inhibited by the novel regulator TcdC. *Mol. Microbiol.* 64, 1274–1288. doi: 10.1111/j.1365-2958.2007.05739.x

Monot, M., Boursaux-Eude, C., Thibonnier, M., Vallenet, D., Moszer, I., Medigue, C., et al. (2011). Reannotation of the genome sequence of *Clostridium difficile* strain 630. *J. Med. Microbiol.* 60, 1193–1199. doi: 10.1099/jmm.0.030452-0

Monot, M., Orgeur, M., Camiade, E., Brehier, C., and Dupuy, B. (2014). COV2HTML: a visualization and analysis tool of bacterial next generation sequencing (NGS) data for postgenomics life scientists. *OMICS* 18, 184–195. doi: 10.1089/omi.2013.0119

Murakami, K. S., and Darst, S. A. (2003). Bacterial RNA polymerases: the whole story. *Curr. Opin. Struct. Biol.* 13, 31–39. doi: 10.1016/s0959-440x(02)00005-2

Nakagawa, S., Niimura, Y., Miura, K., and Gojobori, T. (2010). Dynamic evolution of translation initiation mechanisms in prokaryotes. *Proc. Natl. Acad. Sci. U.S.A.* 107, 6382–6387. doi: 10.1073/pnas.1002036107

Neumann-Schaal, M., Hofmann, J. D., Will, S. E., and Schomburg, D. (2015). Time-resolved amino acid uptake of *Clostridium difficile* 630Deltaerm and concomitant fermentation product and toxin formation. *BMC Microbiol.* 15:281. doi: 10.1186/1471-2105-15-281

Neumann-Schaal, M., Jahn, D., and Schmidt-Hohagen, K. (2019). Metabolism the difference way: the key to the success of the pathogen *Clostridioides difficile*. *Front. Microbiol.* 10:219. doi: 10.3389/fcimb.2017.00219

Nie, X., Dong, W., and Yang, C. (2019). Genomic reconstruction of sigma(54) regulons in *Clostridiales*. *BMC Genomics* 20:565. doi: 10.1186/1471-2105-15-565

Papenfort, K., Forstner, K. U., Cong, J. P., Sharma, C. M., and Bassler, B. L. (2015). Differential RNA-seq of *Vibrio cholerae* identifies the VqmR small RNA as a regulator of biofilm formation. *Proc. Natl. Acad. Sci. U.S.A.* 112, E766–E775.

Pereira, F. C., Saujet, L., Tome, A. R., Serrano, M., Monot, M., Couture-Tosi, E., et al. (2013). The spore differentiation pathway in the enteric pathogen *Clostridium difficile*. *PLoS Genet.* 9:e1003782. doi: 10.1371/journal.pgen.1003782

Pishdadian, K., Fimlaid, K. A., and Shen, A. (2015). SpoIID-mediated regulation of sigmaK function during *Clostridium difficile* sporulation. *Mol. Microbiol.* 95, 189–208. doi: 10.1111/mmi.12856

Rosinski-Chupin, I., Sauvage, E., Fouet, A., Poyart, C., and Glaser, P. (2019). Conserved and specific features of *Streptococcus pyogenes* and *Streptococcus agalactiae* transcriptional landscapes. *BMC Genomics* 20:236. doi: 10.1186/1471-2105-15-236

Rosinski-Chupin, I., Sauvage, E., Sismeiro, O., Villain, A., Da Cunha, V., Caliot, M. E., et al. (2015). Single nucleotide resolution RNA-seq uncovers new regulatory mechanisms in the opportunistic pathogen *Streptococcus agalactiae*. *BMC Genomics* 16:419. doi: 10.1186/1471-2105-15-419

Rupnik, M., Wilcox, M. H., and Gerding, D. N. (2009). *Clostridium difficile* infection: new developments in epidemiology and pathogenesis. *Nat. Rev. Microbiol.* 7, 526–536. doi: 10.1038/nrmicro2164

Sahr, T., Rusniok, C., Dervins-Ravault, D., Sismeiro, O., Coppee, J. Y., and Buchrieser, C. (2012). Deep sequencing defines the transcriptional map of *L. pneumophila* and identifies growth phase-dependent regulated ncRNAs implicated in virulence. *RNA Biol.* 9, 503–519. doi: 10.4161/rna.20270

Saliba, A. E., Santos, S. C., and Vogel, J. (2017). New RNA-seq approaches for the study of bacterial pathogens. *Curr. Opin. Microbiol.* 35, 78–87. doi: 10.1016/j.mib.2017.01.001

Saujet, L., Monot, M., Dupuy, B., Soutourina, O., and Martin-Verstraete, I. (2011). The key sigma factor of transition phase, SigH, controls sporulation, metabolism, and virulence factor expression in *Clostridium difficile*. *J. Bacteriol.* 193, 3186–3196. doi: 10.1128/jb.00272-11

Saujet, L., Pereira, F. C., Henriques, A. O., and Martin-Verstraete, I. (2014). The regulatory network controlling spore formation in *Clostridium difficile*. *FEMS Microbiol. Lett.* 358, 1–10. doi: 10.1111/1574-6968.12540

Saujet, L., Pereira, F. C., Serrano, M., Soutourina, O., Monot, M., Shelyakin, P. V., et al. (2013). Genome-wide analysis of cell type-specific gene transcription during spore formation in *Clostridium difficile*. *PLoS Genet.* 9:e1003756. doi: 10.1371/journal.pgen.1003756

Sebaihia, M., Wren, B. W., Mullany, P., Fairweather, N. F., Minton, N., Stabler, R., et al. (2006). The multidrug-resistant human pathogen *Clostridium difficile* has a highly mobile, mosaic genome. *Nat. Genet.* 38, 779–786. doi: 10.1038/ng1830

Shao, W., Price, M. N., Deutschbauer, A. M., Romine, M. F., and Arkin, A. P. (2014). Conservation of transcription start sites within genes across a bacterial genus. *mBio* 5:e1398-14.

Sharma, C. M., Hoffmann, S., Darfeuille, F., Reignier, J., Findeiss, S., Sittka, A., et al. (2010). The primary transcriptome of the major human pathogen *Helicobacter pylori*. *Nature* 464, 250–255. doi: 10.1038/nature08756

Sharma, C. M., and Vogel, J. (2014). Differential RNA-seq: the approach behind and the biological insight gained. *Curr. Opin. Microbiol.* 19, 97–105. doi: 10.1016/j.mib.2014.06.010

Sierro, N., Makita, Y., De Hoon, M., and Nakai, K. (2008). DBTBS: a database of transcriptional regulation in *Bacillus subtilis* containing upstream intergenic conservation information. *Nucleic Acids Res.* 36, D93–D96.

Sineva, E., Savkina, M., and Ades, S. E. (2017). Themes and variations in gene regulation by extracytoplasmic function (ECF) sigma factors. *Curr. Opin. Microbiol.* 36, 128–137. doi: 10.1016/j.mib.2017.05.004

Smyth, G. K., and Speed, T. (2003). Normalization of cDNA microarray data. *Methods* 31, 265–273. doi: 10.1016/s1046-2023(03)00155-5

Sorek, R., and Cossart, P. (2010). Prokaryotic transcriptomics: a new view on regulation, physiology and pathogenicity. *Nat. Rev. Genet.* 11, 9–16. doi: 10.1038/nrg2695
Soutourina, O. A., Monot, M., Boudry, P., Saujet, L., Pichon, C., Sismeiro, O., et al. (2013). Genome-wide identification of regulatory RNAs in the human pathogen *Clostridium difficile*. *PLoS Genet.* 9:e1003493. doi: 10.1371/journal.pgen.1003493

Stabler, R. A., Gerding, D. N., Songer, J. G., Drudy, D., Brazier, J. S., Trinh, H. T., et al. (2006). Comparative phylogenomics of *Clostridium difficile* reveals clade specificity and microevolution of hypervirulent strains. *J. Bacteriol.* 188, 7297-7305. doi: 10.1128/jb.00664-06

Stevenson, E., Minton, N. P., and Kuehne, S. A. (2015). The role of flagella in *Clostridium difficile* pathogenicity. *Trends Microbiol* 23, 275-282. doi: 10.1016/j.tim.2015.01.004

Stulke, J., Arnaud, M., Rapoport, G., and Martin-Verstraete, I. (1998). PRD-a protein domain involved in PTS-dependent induction and carbon catabolic repression of catabolic operons in bacteria. *Mol. Microbiol.* 28, 865-874. doi: 10.1046/j.1365-2958.1998.00839.x

Sudarsan, N., Lee, E. R., Weinberg, Z., Moy, R. H., Kim, J. N., Link, K. H., et al. (2008). Riboswitches in eubacteria sense the second messenger cyclic di-GMP. *Science* 321, 411-413.

Trinquier, A., Durand, S., Braun, F., and Condon, C. (2020). Regulation of RNA processing and degradation in bacteria. *Biochim. Biophys. Acta Gene Regul. Mech.* 1863:194505. doi: 10.1016/j.bbagrm.2020.194505

Vallenet, D., Calteau, A., Dubois, M., Amours, P., Bazin, A., Beuvin, M., et al. (2020). MicroScope: an integrated platform for the annotation and exploration of microbial gene functions through genomic, pangenomic and metabolic comparative analysis. *Nucleic Acids Res.* 48, D579-D589.

Vedantam, G., Clark, A., Chu, M., Mcquade, R., Mallozzi, M., and Viswanathan, V. (2012). *Clostridium difficile* infection: toxins and non-toxin virulence factors,

and their contributions to disease establishment and host response. *Gut Microb.* 3, 121-134. doi: 10.4161/gmic.19399

Woods, E. C., Nawrocki, K. L., Suarez, J. M., and Mcbride, S. M. (2016). The *Clostridium difficile* Dlt pathway is controlled by the extracytoplasmic function sigma Factor sigmaV in response to Lysozyme. *Infect. Immun.* 84, 1902-1916. doi: 10.1128/iai.00207-16

Wurtzel, O., Sapra, R., Chen, F., Zhu, Y., Simmons, B. A., and Sorek, R. (2010). A single-base resolution map of an archaeal transcriptome. *Genome Res.* 20, 133-141. doi: 10.1101/gr.100396.109

Wurtzel, O., Sesto, N., Mellin, J. R., Karunker, I., Edelheit, S., Becavin, C., et al. (2012a). Comparative transcriptomics of pathogenic and non-pathogenic *Listeria* species. *Mol. Syst. Biol.* 8:583. doi: 10.1038/msb.2012.11

Wurtzel, O., Yoder-Himes, D. R., Han, K., Dandekar, A. A., Edelheit, S., Greenberg E. P., et al. (2012b). The single-nucleotide resolution transcriptome of *Pseudomonas aeruginosa* grown in body temperature. *PLoS Pathog.* 8:e1002945. doi: 10.1371/journal.pgen.1002945

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Soutourina, Dubois, Monot, Shelyakin, Saujet, Boudry, Gelfand, Dupuy and Martin-Verstraete. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.