Sec. Infectious Agents and Disease
An Optimized Method for Extracting Bacterial RNA from Mouse Skin Tissue Colonized by Mycobacterium ulcerans
- 1Center for Research in Cancerology and Immunology Nantes-Angers, Institut National de la Santé et de la Recherche Médicale, Université de Nantes, Université d’Angers, Angers, France
- 2Equipe Atip-Avenir, Center for Research in Cancerology and Immunology Nantes-Angers, Institut National de la Santé et de la Recherche Médicale, Centre Hospitalier Universitaire et Université d’Angers, Angers, France
- 3Transcriptome and Epigenome Platform, Biomics, Center for Innovation and Technological Research, Institut Pasteur, Paris, France
Bacterial transcriptome analyses during host colonization are essential to decipher the complexity of the relationship between the bacterium and its host. RNA sequencing (RNA-seq) is a promising approach providing valuable information about bacterial adaptation, the host response and, in some cases, mutual tolerance underlying crosstalk, as recently observed in the context of Mycobacterium ulcerans infection. Buruli ulcer is caused by M. ulcerans. This neglected disease is the third most common mycobacterial disease worldwide. Without treatment, M. ulcerans provokes massive skin ulcers. A healing process may be observed in 5% of Buruli ulcer patients several months after the initiation of disease. This spontaneous healing process suggests that some hosts can counteract the development of the lesions caused by M. ulcerans. Deciphering the mechanisms involved in this process should open up new treatment possibilities. To this end, we recently developed the first mouse model for studies of the spontaneous healing process. We have shown that the healing process is based on mutual tolerance between the bacterium and its host. In this context, RNA-seq seems to be the most appropriate method for deciphering bacterial adaptation. However, due to the low bacterial load in host tissues, the isolation of mycobacterial RNA from skin tissue for RNA-seq analysis remains challenging. We developed a method for extracting and purifying mycobacterial RNA whilst minimizing the amount of host RNA in the sample. This approach was based on the extraction of bacterial RNA by a differential lysis method. The challenge in the development of this method was the choice of a lysis system favoring the removal of host RNA without damage to the bacterial cells. We made use of the thick, resistant cell wall of M. ulcerans to achieve this end.
RNA sequencing (RNA-seq), a next-generation sequencing technique, opens up unique opportunities for deciphering interactions between microorganisms and their hosts. It provides information about the relative levels of expression of the various genes and support for proteomic results. RNA-seq can also be used to identify the regulatory networks controlled by non-coding RNA, as reported for Mycobacterium tuberculosis (Arnvig et al., 2011).
RNA sequencing can be used to investigate and characterize the different facets of bacterial life, including, in particular, the host/microorganism tolerance underlying crosstalk between the two species (La et al., 2008; Skvortsov and Azhikina, 2010; Westermann et al., 2012; Nalpas et al., 2013; Szafranska et al., 2014; Amorim-Vaz et al., 2015). In this context, M. ulcerans is a fascinating microorganism, with a complex biology due to the different facets of its life cycle. M. ulcerans is the causal agent of Buruli ulcer, a severe cutaneous infection (Vincent et al., 2014) and the third most frequent mycobacterial disease worldwide, after tuberculosis and leprosy (Asiedu et al., 2000). M. ulcerans has developed sophisticated strategies for colonizing various hosts, from aquatic organisms (aquatic plants, insects, etc.) to humans, suggesting a “parasite lifestyle” (Portaels et al., 1999, 2001; Marsollier et al., 2002, 2004, 2005, 2007a,b; Johnson et al., 2007; Merritt et al., 2010; Garchitorena et al., 2014, 2015; Marion et al., 2014a, 2016a; Zogo et al., 2015; Sanhueza et al., 2016).
Mycobacterium ulcerans colonizes human tissues in several phases. Following the inoculation of the dermis with M. ulcerans, there is an intracellular phase of infection, in which the bacterium remains within macrophages and neutrophils, allowing it to evade immune system recognition (Torrado et al., 2007). M. ulcerans then kills the host macrophage by producing mycolactone, a lipid toxin, initiating an extracellular stage, in which local mycolactone concentrations increase considerably, leading to massive host tissue destruction. During these two stages, mycolactone is not only cytotoxic, it also modulates the immune system, modifying cytokine production and acting on the peripheral nervous system to induce the formation of a painless lesion (George et al., 1999; Coutanceau et al., 2005; Oliveira et al., 2005; Torrado et al., 2007, 2010; Silva et al., 2009; Fraga et al., 2010, 2012; Marion et al., 2014b). These pleiotropic effects of mycolactone facilitate host colonization by this bacillus. This toxin is a distinctive feature of M. ulcerans and seems to play a key role in its eco-epidemiology and pathogenesis.
We recently showed that 5% of Buruli ulcer patients display spontaneous healing without treatment (Marion et al., 2016a). This clinically relevant observation demonstrates that patients can develop responses that counteract the effects of M. ulcerans and its toxin. Deciphering the mechanisms involved in this process will open up new therapeutic strategies.
We have developed the first dedicated mouse model for studies of the spontaneous healing process (Marion et al., 2016b). During the characterization of this model, we made an interesting discovery concerning the dynamics of viable bacterial load in healed tissues: the load of cultivable bacilli was found to be both high and stable in the long term (Marion et al., 2016b). We then demonstrated that mycolactone synthesis was inhibited in healed tissues. Surprisingly, transcriptomic studies based on RT-qPCR showed that the bacteria in these tissues were not dormant. Paradoxically, transcription levels for the principal genes involved directly in toxin synthesis were unaffected, suggesting that mycolactone synthesis was regulated upstream, as already shown in vitro (Deshayes et al., 2013). Our previous findings suggest that RNA-seq is the most appropriate approach for deciphering the regulation of mycolactone synthesis in vivo.
This approach requires the isolation of large amounts of high-quality bacterial RNA, which is challenging in studies performed in vivo, because host RNA is much more abundant than bacterial RNA in samples. It was therefore necessary to optimize the method for extracting RNA from tissues, so as to minimize the amount of host RNA in the sample whilst ensuring the isolation of sufficient quantities of high-quality mycobacterial RNA. We present here this optimized method, based on differential lysis for the analysis of the whole transcriptome of M. ulcerans.
Materials and Methods
Ethics Statement for Animal Experiments
All animal experiments were performed in accordance with national guidelines (articles R214–87 to R214–90 from the French “rural code”) and European guidelines (directive 2010/63/EU of the European Parliament and of the council of September 22, 2010 on the protection of animals used for scientific purposes). All protocols were approved by the ethics committee of the Pays de la Loire region, under protocol nos. CEEA 2009.14 and CEEA 2012.145, and performed at the required biosafety level. Animals were maintained under specific pathogen-free conditions in the animal house facility of Angers University Hospital, France (agreement A 49 007 002).
Bacterial Strains and Inoculation
Mycobacterium ulcerans strain 01G897 was originally isolated from patients from French Guiana (De Gentile et al., 1992). Bacterial suspensions were prepared as previously described (Marsollier et al., 2007b; Marion et al., 2016b), with adjustment to 2 × 105 acid-fast bacilli/ml for inoculation (50 μl) into the tail of 6-week-old females of the inbred FVB/N mouse strain (Charles River Laboratories, Saint-Germain-Nuelles, France).
RNA Extraction and Purification
RNA was extracted from infected tail skins 30 days post-infection, by the Trizol/chloroform method (Method 1), which co-extracts host and bacterial RNA, or by the differential lysis method (Method 2), optimized and adapted from that described by (Rustad et al., 2009), for the isolation of bacterial RNA alone.
Method 1: Total RNA Extraction
(i) Sample preparation: Tail skin was excised from infected mice and immediately placed in a Petri dish containing a mixture of 1 ml Trizol (Ambion) and 1 ml RLT buffer (Qiagen) supplemented with 1% β-mercaptoethanol. Skin tissues were cut into smaller pieces, transferred to round-bottomed tubes and broken up with a TissueRuptor (Qiagen).
(ii) RNA extraction and purification: Samples were transferred to two bead beating tubes (0.1 mm glass beads, MoBio) and were shaken with TissueLyser (Qiagen) at 4°C for 5 min at 30 Hz. The samples were immediately placed on ice and centrifuged at 10,000 × g for 5 min at 4°C to remove cell debris. The supernatant was transferred to a 15 ml tube containing 1 ml Trizol. 200 μl chloroform/isoamyl alcohol (24:1) was added to the tube, which was then repeatedly inverted to mix and centrifuged at 10,000 × g for 5 min. The aqueous phase was transferred to a clean tube containing 400 μl chloroform/isoamyl alcohol and centrifuged again, as in the previous step. The aqueous phase was transferred to a clean tube containing 1 volume of 70% ethanol, and the tube was repeatedly gently inverted to mix. Total RNA was purified with the RNeasy Midi kit (Qiagen), with DNase treatment, according to the manufacturer’s protocol, and eluted in 100 μl RNase- and DNase-free water.
(iii) DNase treatment: Contaminating DNA was removed by retreating the RNA with DNase, for 45 min at 37°C, with the TURBO DNA-free kit (Ambion), according to the manufacturer’s protocol.
(iv) Eukaryotic RNA removal: To reduce the levels of contaminating host RNA from the samples, MICROBEnrich kit (Ambion) was used according to the manufacturer’s protocol.
Method 2: Bacterial RNA Isolation by Differential Lysis Method
(i) Sample preparation: Tail skin was excised from infected mice and immediately placed in a Petri dish containing 2 ml Tris-EDTA (TE) buffer (10 mM Tris-HCl, 1 mM EDTA). Skin tissues were cut into smaller pieces, transferred to round-bottomed tubes and broken up with a TissueRuptor (Qiagen). 4 ml TE buffer was added and the tissue homogenates were digested with 1 ml of a 20 mg/ml proteinase K solution (Qiagen) for 10 min at 55°C (without shaking). Another 6 ml of TE buffer was added and the samples were centrifuged at 3,200 × g for 15 min at 4°C. The pellet, which contained the bacterial cells, was resuspended in 300 μl Tri-reagent (Zymo Research) and 300 μl RLT buffer (Qiagen) supplemented with 1% β-mercaptoethanol.
(ii) RNA extraction and purification: The samples were transferred into a bead beating tube (0.1 mm glass beads, MoBio) and shaken with TissueLyser (Qiagen) at 4°C for 5 min at 30 Hz. The samples were immediately placed on ice and centrifuged at 10,000 × g for 5 min at 4°C to remove cell debris. The supernatant was transferred into a clean tube containing 1 volume of 100% ethanol, with which it was mixed by repeated gentle inversion. The RNA was purified and treated with DNase with the Direct-zol RNA MiniPrep kit (Zymo Research), according to the manufacturer’s protocol, and eluted in 50 μl of RNase- and DNase-free water.
(iii) DNase treatment: Contaminating DNA was removed by retreating the RNA with DNase for 45 min at 37°C with the TURBO DNA-free kit (Ambion), according to the manufacturer’s protocol.
(i) Quantification and purity analysis: Total RNA concentration and purity (A260nm/A280nm) were assessed with 1 μl of the RNA preparation, on a NanoDrop 1000TM (Thermo Scientific) spectrophotometer. A ratio greater than 1.8 is usually considered to indicate satisfactory RNA purity (Imbeaud et al., 2005).
(ii) Quality and integrity analysis: RNA quality and integrity were assessed with the Experion automated electrophoresis system (Bio-Rad). The total RNA sample (Method 1) was diluted 10-fold and the Experion RNA StdSens analysis chip was used (quantification of 5–500 ng/μl RNA). With enriched bacterial RNA preparations (Method 2), the RNA was diluted twofold and the Experion RNA HighSens analysis chip was used (quantification of 100–5,000 pg/μl RNA). The RNA quality indicator (RQI) method returns a number between 1 (highly degraded RNA) and 10 (intact RNA) for each RNA sample (Imbeaud et al., 2005; Schroeder et al., 2006).
Transcriptional Analysis by RT-qPCR
RT-qPCR targeting the M. ulcerans ppk gene was performed to detect mycobacterial RNA transcripts in RNA samples. RT-qPCR targeting the mouse gapdh gene was performed to evaluate the contamination of RNA samples with host RNA. The ppk and gapdh genes were selected as housekeeping genes for the RT-qPCR analyses for M. ulcerans and mouse, respectively. Amplification efficiency (Table 1) was determined from the slope of a standard curve of cDNA serial dilutions.
(i) Reverse transcription: The first-strand cDNA was synthesized in a reaction volume of 20 μl containing 2 μl of total RNA, 500 ng of random primers (Invitrogen) and the M-MLV reverse transcriptase (Invitrogen). Contaminating DNA in RNA sample was checked by performing a negative control with no reverse transcriptase (RT-) for each sample.
(ii) Quantitative real-time PCR: qPCR was performed in a reaction volume of 10 μl containing Absolute Blue qPCR mix (Thermo Scientific), 300 nM primers, 100 nM Taqman probe (Table 1) and 2.5 μl of diluted twofold dilution of cDNA/RT-. The sequences of the primers and probes used are provided in Table 1. Reactions were run on a AriaMx Thermocycler (Agilent), with the following program: 10 min at 95°C and 40 cycles of 10 s at 95°C and 1 min at 60°C. Each sample was analyzed in duplicate.
(i) Ribosomal RNA (rRNA) removal: rRNA was depleted from the mycobacterial total RNA preparations with the RiboZero Epidemiology Illumina kit, which removes eukaryote and prokaryote rRNA in a single step.
(ii) Preparation of RNA-seq libraries: The RNA-seq libraries were prepared with the TruSeq Stranded Total RNA LT Sample Prep kit (Illumina). The quality of all libraries was checked with the DNA-1000 kit (Agilent) on a 2100 Bioanalyzer and quantification was performed with Quant-It assays on a Qubit 1.0 fluorometer (Invitrogen).
(iii) RNA-seq: Clusters were generated for the resulting libraries, with Illumina HiSeq SR Cluster Kit v4 reagents. Sequencing was performed with the Illumina HiSeq 2500 system and HiSeq SBS kit v4 reagents. Runs were carried out over 65 cycles, including seven indexing cycles, to obtain 65-bp single-end reads. Sequencing data were processed with Illumina Pipeline software (Casava version 1.9). All 65-bp reads were aligned against the complete genome and plasmid sequences of M. ulcerans (Agy99 strain) obtained from the Burulist database with Bowtie software, and against the mouse genome (GRCm38) obtained from the Ensembl database with STAR software. Data were normalized and analyzed in R, with the Bioconductor packages.
We recently developed a model for studies of the spontaneous healing process. Our studies in this model revealed that the host response modulated toxin synthesis, providing evidence for crosstalk between M. ulcerans and the host. We therefore decided to develop a method for extracting bacterial RNA from host tissue in conditions suitable for high-throughput RNA-seq, to make it possible to decipher the regulation of toxin production. We evaluated two methods for isolating M. ulcerans RNA from host tissue (tail) for RNA-seq analysis. In this study, we inoculated the tails of mice that were then killed for analysis 30 days post-infection. All infected mice presented clinical edemas and bacterial load was estimated at about 106 CFU (8.9 × 105–2.1 × 106 CFU). Each of the RNA extraction methods was evaluated on three samples (taken from three mice).
Method 1: Extraction of Total RNA from Tail Skin Infected with M. ulcerans
For studies of the interactions between the host and the bacteria during the infection of mice with M. ulcerans, we initially tried to extract bacterial and host RNA together, from the tail skin of mice infected with M. ulcerans. The infected tissues were broken up with a TissueRuptor and subjected to chemical and mechanical lysis, to release the bacterial and eukaryotic RNA. The RNA was then subjected to Trizol/chloroform extraction and purified on a Qiagen midi column which allows the purification of up to 1 mg RNA (Figure 1).
FIGURE 1. Scheme for total RNA extraction (Method 1) and comparison with the differential lysis method (Method 2).
With this approach, we obtained 200 to 300 μg of total RNA per sample (Figure 2C). Electrophoretic analysis revealed strong bands corresponding to the 18S and 28S rRNA of eukaryotic cells, with no bacterial RNA bands corresponding to 16S and 23S rRNA in samples (Figure 2A and Supplementary Figure S1). The electropherogram trace confirmed this result, as it contained peaks only for mouse rRNA (Figure 2B and Supplementary Figure S1). The mycobacterial RNA transcripts in samples were assessed by a reverse transcription-polymerase chain reaction (RT-qPCR) targeting the M. ulcerans housekeeping gene, ppk. The samples contained only small amounts of mycobacterial RNA, with Ct values of 32–33 (Table 2). The A260nm/A280nm ratio exceeded 1.8, confirming the purity of the RNA in each sample (Figure 2C), and electrophoresis showed a complete absence of RNA degradation, with an RQI value greater than 9 (Figure 2C). However, we were unable to assess the integrity of bacterial RNA, because no bacterial RNA was visualized.
FIGURE 2. Evaluation of the coextraction method (Method 1). (A) Image of the electrophoresis gel for three independent RNA samples extracted from mice infected with M. ulcerans. L, RNA ladder. (B) Electropherogram: example of sample RNA derived from lane 1. (C) Table with RNA concentrations (ng/μl), indicator of purity (ratio A260/A280), integrity of RNA samples (RQI, RNA Quality Indicator), and total quantity (μg) of RNA extracted with Method 1. (D) Electropherogram: example of sample RNA before and after MicrobEnrich (ME) treatment. The 18S and 28S (mouse), and 16S and 23S (bacterial) rRNA bands are indicated in red and green, respectively.
TABLE 2. Comparison of Ct values between the coextraction method (Method 1) and the differential lysis method (Method 2).
To eliminate the host RNA from the samples, MICROBEnrich kit (Ambion) was used which selectively removes 18S, 28S rRNA, and polyadenylated mRNA. Following MICROBEnrich treatment, electrophoresis was performed and showed this treatment to have been poorly effective in our hands, because strong bands corresponding to the 18S and 28S rRNAs were still detected (Figure 2D). These data confirmed the predominance of eukaryotic ribosomal RNA in samples, and suggested that mouse mRNA would probably also predominate over bacterial transcripts. However, based on our qPCR demonstrating the presence of mycobacterial RNA in samples, we decided to attempt a deep sequencing approach (RNA-seq) on a sample before and after bacterial RNA enrichment with the MICROBEnrich kit. The obtained RNA-seq data showed an absence of reads aligned with the M. ulcerans genome, even though mycobacterial RNA transcripts were detected by RT-qPCR, a more sensitive technique. In this case, sequencing depth was insufficient to obtain reads mapping to the M. ulcerans genome, confirming that the removal of mouse RNA was not effective enough to achieve an enrichment of the sample in mycobacterial RNA. In conclusion, these data demonstrate that studies of the M. ulcerans transcriptome in vivo will require a method capable of significantly decreasing eukaryotic RNA levels, to facilitate the detection of mycobacterial transcripts by RNA-seq.
Method 2: Bacterial RNA Isolation by the Differential Lysis Method
We developed an alternative technique for the purification of bacterial RNA from infected host tissues by differential lysis, to achieve our goal. The principal challenge in this method was the choice of a lysis system resulting in host-cell lysis without damage to the bacterial cells. We recently demonstrated that the mechanical disruption of tissues had no effect on bacterial viability (Marion et al., 2016b). We therefore disrupted the host cells by mechanical lysis and treated them with proteinase K, to digest and degrade the mouse tissue whilst leaving the bacterial cells intact. The intact bacterial cells were then separated from the lysate containing host RNA by centrifugation, and the pellet containing the bacteria was resuspended in lysis buffer and subjected to bead beating. The bacterial RNA was purified on a Zymo column, making it possible to purify up to 100 μg of RNA directly from samples in Tri-reagent, a much less time-consuming approach (Figure 1).
We obtained 150 to 500 ng of total RNA per sample (Figure 3C), about one-thousandth the amount obtained with Method 1 (Figure 2C). However, gel electrophoresis on the total RNA obtained with this method revealed the presence of bands corresponding to the 16S and 23S bacterial rRNA whereas the mouse 28S and 18S rRNA were not detected (Figure 3A and Supplementary Figure S1). This finding was confirmed by the electropherogram trace (Figure 3B and Supplementary Figure S1).
FIGURE 3. Evaluation of the differential lysis method (Method 2). (A) Image of the electrophoresis gel for three independent RNA samples extracted from mice infected with M. ulcerans. L, RNA ladder. The 16S and 23S (bacterial) rRNA bands are indicated in green. (B) Electropherogram: example of sample RNA derived from lane 3. (C) Table with RNA concentrations (ng/μl), indicator of purity (ratio A260/A280), integrity of RNA samples (RQI, RNA quality indicator), and total quantity (ng) of RNA extracted with Method 2.
The enrichment of samples in mycobacterial RNA by the differential lysis method (Method 2) was evaluated by RT-qPCR. The efficacy of eukaryotic RNA removal during sample preparation with Method 2 was evaluated by comparing equivalent fractions of enriched (Method 2) and non-enriched (Method 1) RNA samples. We quantified the expression of both bacterial and mouse genes in the samples obtained with the two methods.
No PCR product was obtained if the amplification was performed without prior reverse transcription (RT-). There was therefore no detectable DNA contamination. A comparison of ppk transcription levels showed that the transcript of this gene was present in both cases, with similar Ct values (30–33). An analysis of murine gapdh gene transcription revealed a difference in Ct of 15 between Method 1 and Method 2 (Table 2), corresponding to the presence of five orders of magnitude less host RNA with Method 2 than with Method 1. These results clearly demonstrate the efficacy of Method 2 for removing eukaryotic RNA. This method was, thus, considered suitable for the enrichment of samples in M. ulcerans RNA. The purity and integrity of bacterial RNA were assessed by determining the A260nm/A280nm ratio and the RQI value, respectively. The A260nm/A280nm ratio was greater than 1.8 (Figure 3C) and RQI values ranged from 7.3 to 9 (Figure 3C), indicating that the RNA was of sufficiently high quality for use in subsequent experiments, including transcriptional analyses.
Given the promising nature of these findings, we performed deep sequencing (RNA-seq) on the enriched samples. The removal of rRNA is a crucial step in RNA-seq, because rRNA signals can prevent adequate coverage of the bacterial transcriptome. We therefore used the Ribo-Zero gold rRNA Removal Kit (Epidemiology, Epicentre) to remove both mouse and bacterial rRNA before producing the cDNA library. The RNA-seq data confirmed that less than 5% of the reads mapped to rRNA, demonstrating the efficacy of the rRNA removal kit. The datasets for samples M2-1, M2-2 and M2-3 obtained by the optimized method (Method 2) contained 9.5, 7.6 and 7.9 million reads aligning with the M. ulcerans genome, respectively, whereas no reads aligned with this genome sequence were detected with the samples obtained by Method 1 (Figure 4). This enrichment procedure therefore provides sufficient coverage of sequences mapping to the M. ulcerans, because a minimum of 2 to 5 million reads from a ribosomal RNA-depleted library is required to provide adequate coverage of the gene expression profiles of bacteria in RNA-seq experiments (Rienksma et al., 2015).
FIGURE 4. Percentage of reads aligned with the M. ulcerans and mouse genomes. M1, Total RNA sample obtained with Method 1; M1_ME, Total RNA sample obtained with Method 1 after enrichment with MICROBEnrich (ME); M2, Bacterial RNA-enriched sample obtained by differential lysis (Method 2). M2-1, M2-2, M2-3 represent three replicates for Method 2.
Finally, two biological replicates, M2-1 and M2-2, were used to evaluate the reproducibility of the gene expression profiles obtained with RNA-seq technology. Spearman’s coefficient of correlation between samples (r = 0.9813) indicated a similar overall pattern of relative gene expression in the biological replicates, indicating that Method 2 was reproducible (Figure 5). This optimized method can therefore be used to study the mycobacterial transcriptome in a mouse model.
FIGURE 5. Reproducibility of Method 2 for gene expression profiles obtained by RNA-seq. Two biological replicates (M2-1 and M2-2) prepared with the same procedure (Method 2) were used to evaluate reproducibility in Spearman’s correlation tests. Each dot corresponds to a single gene. The x and y axes indicate the level of expression of a single gene in two different samples (M2-1 and M2-2). Spearman’s correlation coefficient = 0.9813.
Over the last few years, RNA-seq has become a powerful tool for studies of complex interactions between microorganisms and their hosts (La et al., 2008; Skvortsov and Azhikina, 2010; Westermann et al., 2012; Szafranska et al., 2014; Amorim-Vaz et al., 2015). This method is highly suitable for studies aiming to decipher the complexity of regulation during interactions between host and bacteria.
RNA sequencing provides ready access to information about the host response during the different phases of colonization by the microorganism, but it is much more difficult to investigate transcription patterns in the microorganism. This is the principal reason for which most reports have focused essentially on host gene expression (Jenner and Young, 2005; Hossain et al., 2006; La et al., 2008). Bacterial transcriptomic signatures are much less widely documented (Jansen and Yu, 2006; La et al., 2008) because bacterial mRNA yields are much lower than those for host RNA (Lucchini et al., 2001; Hinton et al., 2004; Ferreira-Machado et al., 2015). It is therefore necessary to overcome this technical difficulty to decipher, through transcriptomic approaches, the basis of microbial attack and the ways in which the microbe escapes host defenses.
The spontaneous healing process observed during M. ulcerans infection highlights the ability of the host to counteract the lesion development caused by the bacterium. We recently developed the first mouse model for investigating this phenomenon (Marion et al., 2016b). Our investigations have revealed a strong inhibition of mycolactone synthesis during spontaneous healing. In conclusion, this model provides a unique opportunity to understand the regulation of mycolactone synthesis and to identify possible drug targets. To this end, RNA-seq seems to be the most appropriate method. However, due to the low bacterial load in host tissues, the isolation of mycobacterial RNA from skin tissue for RNA-seq analysis remains challenging.
In this context, we initially attempted to co-extract host and bacterial RNA together, for studies of the total transcriptome of M. ulcerans in vivo. This approach is hindered by a technical issue for bacterial transcriptome analysis: lower yields for bacterial RNA than for host RNA. The high abundance of host RNA reduced the coverage of sequences mapping to the M. ulcerans genome. Commercial kits, such as the MicrobEnrich kit, have been developed to overcome this problem, and these kits have been successfully used in some studies (Bergman et al., 2007; Audia et al., 2008; Bielecki et al., 2011; Date et al., 2014). However, in our context, use of the MicrobEnrich kit did not provide a sufficiently high level of enrichment in mycobacterial RNA, because no mycobacterial transcripts were detected by RNA-seq. This may reflect the large proportion of host RNA in sample. We therefore developed an alternative method for studying the whole transcriptome of M. ulcerans in vivo by RNA-seq, in which mycobacterial enrichment was increased during the sample preparation process.
This method is based on the extraction of bacterial RNA by a differential lysis method. The principal challenge in this method is the choice of a system capable of lysing the host cells without damaging the bacterial cells. Chemical or mechanical differential lysis methods are generally used for RNA isolation in vivo (Schnappinger et al., 2003; Talaat et al., 2004, 2007; Tuanyok et al., 2006). Our method took advantage of the thick, resistant wall of M. ulcerans, which is generally considered problematic for the extraction of RNA, DNA and proteins. This resistance of the bacterial wall makes it possible (i) to lyse eukaryotic cells without damaging M. ulcerans cells, (ii) to remove large amounts of eukaryotic RNA by simple centrifugation and, (iii) to obtain high-quality bacterial RNA for RNA-seq analysis.
We have developed the first simple protocol for the selective extraction of M. ulcerans RNA from host tissues. This method represents a significant improvement as it provides RNA of sufficiently high quality for RNA-seq analysis. This strategy will make it possible to perform in vivo studies of the interactions between mycobacteria and their host, improving our understanding of the molecular mechanisms underlying infection. Finally, it may be possible to adapt this method for the isolation of other mycobacteria in the host colonization context.
MR-S, LM, and EM conceived and designed the experiments. MR-S, JB, OS, and EM performed the experiments. MR-S, OS, LM and EM performed data analysis. MR-S and LM wrote the paper.
This work was supported by Institut Nationale de la Santé et de la Recherche Médicale (INSERM), Atip Avenir program, Fondation Raoul Follereau-France, Agence Nationale de la Recherche (ANR) BU_SPONT_HEAL project (Infect-Era project and MYCOPARADOX ANR project), Angers University REPAIR project and Société Française de Dermatologie, Région Pays de la Loire (ARMINA project). The Transcriptome and EpiGenome Platform is a member of the France Génomique consortium (ANR10-NBS-09-08).
Conflict of Interest Statement
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.
We thank all members of the SCAHU platform for their technical assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2017.00512/full#supplementary-material
FIGURE S1 | RNA profile of a mixed sample containing bacterial and mouse RNA (M+B), a mouse sample (M) and a M. ulcerans sample (B). (A) Image of the electrophoresis gel for the RNA samples. (B) Electropherogram of RNA samples. The 18S and 28S (mouse), and 16S and 23S (bacterial) rRNA bands are indicated in red and green, respectively.
Amorim-Vaz, S., Tran Vdu, T., Pradervand, S., Pagni, M., Coste, A. T., and Sanglard, D. (2015). RNA enrichment method for quantitative transcriptional analysis of pathogens In Vivo applied to the fungus Candida albicans. MBio 6, e00942–15. doi: 10.1128/mBio.00942-15
Arnvig, K. B., Comas, I., Thomson, N. R., Houghton, J., Boshoff, H. I., Croucher, N. J., et al. (2011). Sequence-based analysis uncovers an abundance of non-coding RNA in the total transcriptome of Mycobacterium tuberculosis. PLoS Pathog. 7:e1002342. doi: 10.1371/journal.ppat.1002342
Audia, J. P., Patton, M. C., and Winkler, H. H. (2008). DNA microarray analysis of the heat shock transcriptome of the obligate intracytoplasmic pathogen Rickettsia prowazekii. Appl. Environ. Microbiol. 74, 7809–7812. doi: 10.1128/AEM.00896-08
Bergman, N. H., Anderson, E. C., Swenson, E. E., Janes, B. K., Fisher, N., Niemeyer, M. M., et al. (2007). Transcriptional profiling of Bacillus anthracis during infection of host macrophages. Infect. Immun. 75, 3434–3444. doi: 10.1128/IAI.01345-06
Bielecki, P., Puchalka, J., Wos-Oxley, M. L., Loessner, H., Glik, J., Kawecki, M., et al. (2011). In-vivo expression profiling of Pseudomonas aeruginosa infections reveals niche-specific and strain-independent transcriptional programs. PLoS ONE 6:e24235. doi: 10.1371/journal.pone.0024235
Coutanceau, E., Marsollier, L., Brosch, R., Perret, E., Goossens, P., Tanguy, M., et al. (2005). Modulation of the host immune response by a transient intracellular stage of Mycobacterium ulcerans: the contribution of endogenous mycolactone toxin. Cell Microbiol. 7, 1187–1196. doi: 10.1111/j.1462-5822.2005.00546.x
Date, S. V., Modrusan, Z., Lawrence, M., Morisaki, J. H., Toy, K., Shah, I. M., et al. (2014). Global gene expression of methicillin-resistant Staphylococcus aureus USA300 during human and mouse infection. J. Infect. Dis. 209, 1542–1550. doi: 10.1093/infdis/jit668
De Gentile, P. L., Mahaza, C., Rolland, F., Carbonnelle, B., Verret, J. L., and Chabasse, D. (1992). Cutaneous ulcer from Mycobacterium ulcerans. Apropos of 1 case in French Guiana. Bull. Soc. Pathol. Exot. 85, 212–214.
Deshayes, C., Angala, S. K., Marion, E., Brandli, I., Babonneau, J., Preisser, L., et al. (2013). Regulation of mycolactone, the Mycobacterium ulcerans toxin, depends on nutrient source. PLoS Negl. Trop. Dis. 7:e2502. doi: 10.1371/journal.pntd.0002502
Ferreira-Machado, A. B., Freitas, M. C., Saji, G. R., Rezende, A. B., Almeida, P. E., Cesar, D. E., et al. (2015). Integrity of prokaryotic mRNA isolated from complex samples for in vivo bacterial transcriptome analysis. Genet. Mol. Res. 14, 14752–14759. doi: 10.4238/2015.November.18.40
Fraga, A. G., Cruz, A., Martins, T. G., Torrado, E., Saraiva, M., Pereira, D. R., et al. (2010). Mycobacterium ulcerans triggers T-cell immunity followed by local and regional but not systemic immunosuppression. Infect. Immun. 79, 421–430. doi: 10.1128/IAI.00820-10
Fraga, A. G., Martins, T. G., Torrado, E., Huygen, K., Portaels, F., Silva, M. T., et al. (2012). Cellular immunity confers transient protection in experimental Buruli ulcer following BCG or mycolactone-negative Mycobacterium ulcerans vaccination. PLoS ONE 7:e33406. doi: 10.1371/journal.pone.0033406
Garchitorena, A., Guegan, J. F., Leger, L., Eyangoh, S., Marsollier, L., and Roche, B. (2015). Mycobacterium ulcerans dynamics in aquatic ecosystems are driven by a complex interplay of abiotic and biotic factors. Elife 4:e07616. doi: 10.7554/eLife.07616
Garchitorena, A., Roche, B., Kamgang, R., Ossomba, J., Babonneau, J., Landier, J., et al. (2014). Mycobacterium ulcerans ecological dynamics and its association with freshwater ecosystems and aquatic communities: results from a 12-month environmental survey in Cameroon. PLoS Negl. Trop. Dis. 8:e2879. doi: 10.1371/journal.pntd.0002879
George, K. M., Chatterjee, D., Gunawardana, G., Welty, D., Hayman, J., Lee, R., et al. (1999). Mycolactone: a polyketide toxin from Mycobacterium ulcerans required for virulence. Science 283, 854–857. doi: 10.1126/science.283.5403.854
Hinton, J. C., Hautefort, I., Eriksson, S., Thompson, A., and Rhen, M. (2004). Benefits and pitfalls of using microarrays to monitor bacterial gene expression during infection. Curr. Opin. Microbiol. 7, 277–282. doi: 10.1016/j.mib.2004.04.009
Imbeaud, S., Graudens, E., Boulanger, V., Barlet, X., Zaborski, P., Eveno, E., et al. (2005). Towards standardization of RNA quality assessment using user-independent classifiers of microcapillary electrophoresis traces. Nucleic Acids Res. 33:e56. doi: 10.1093/nar/gni054
Johnson, P. D., Azuolas, J., Lavender, C. J., Wishart, E., Stinear, T. P., Hayman, J. A., et al. (2007). Mycobacterium ulcerans in mosquitoes captured during outbreak of Buruli ulcer, southeastern Australia. Emerg. Infect. Dis. 13, 1653–1660. doi: 10.3201/eid1311.061369
Marion, E., Chauty, A., Kempf, M., Le Corre, Y., Delneste, Y., Croue, A., et al. (2016a). Clinical features of spontaneous partial healing during Mycobacterium ulcerans infection. Open Forum Infect. Dis. 3:ofw013. doi: 10.1093/ofid/ofw013
Marion, E., Chauty, A., Yeramian, E., Babonneau, J., Kempf, M., and Marsollier, L. (2014a). A case of guilt by association: water bug bite incriminated in M. ulcerans infection. Int. J. Mycobacteriol. 3, 158–161. doi: 10.1016/j.ijmyco.2014.01.004
Marion, E., Jarry, U., Cano, C., Savary, C., Beauvillain, C., Robbe-Saule, M., et al. (2016b). FVB/N mice spontaneously heal ulcerative lesions induced by Mycobacterium ulcerans and switch M. ulcerans into a Low Mycolactone Producer. J. Immunol. 196, 2690–2698. doi: 10.4049/jimmunol.1502194
Marion, E., Song, O. R., Christophe, T., Babonneau, J., Fenistein, D., Eyer, J., et al. (2014b). Mycobacterial toxin induces analgesia in buruli ulcer by targeting the angiotensin pathways. Cell 157, 1565–1576. doi: 10.1016/j.cell.2014.04.040
Marsollier, L., Aubry, J., Coutanceau, E., Andre, J. P., Small, P. L., Milon, G., et al. (2005). Colonization of the salivary glands of Naucoris cimicoides by Mycobacterium ulcerans requires host plasmatocytes and a macrolide toxin, mycolactone. Cell Microbiol. 7, 935–943. doi: 10.1111/j.1462-5822.2005.00521.x
Marsollier, L., Brodin, P., Jackson, M., Kordulakova, J., Tafelmeyer, P., Carbonnelle, E., et al. (2007a). Impact of Mycobacterium ulcerans biofilm on transmissibility to ecological niches and Buruli ulcer pathogenesis. PLoS Pathog. 3:e62. doi: 10.1371/journal.ppat.0030062
Marsollier, L., Deniaux, E., Brodin, P., Marot, A., Wondje, C. M., Saint-Andre, J. P., et al. (2007b). Protection against Mycobacterium ulcerans lesion development by exposure to aquatic insect saliva. PLoS Med. 4:e64. doi: 10.1371/journal.pmed.0040064
Marsollier, L., Robert, R., Aubry, J., Saint Andre, J. P., Kouakou, H., Legras, P., et al. (2002). Aquatic insects as a vector for Mycobacterium ulcerans. Appl. Environ. Microbiol. 68, 4623–4628. doi: 10.1128/AEM.68.9.4623-4628.2002
Marsollier, L., Stinear, T., Aubry, J., Saint Andre, J. P., Robert, R., Legras, P., et al. (2004). Aquatic plants stimulate the growth of and biofilm formation by Mycobacterium ulcerans in axenic culture and harbor these bacteria in the environment. Appl. Environ. Microbiol. 70, 1097–1103. doi: 10.1128/AEM.70.2.1097-1103.2004
Merritt, R. W., Walker, E. D., Small, P. L., Wallace, J. R., Johnson, P. D., Benbow, M. E., et al. (2010). Ecology and transmission of Buruli ulcer disease: a systematic review. PLoS Negl. Trop. Dis. 4:e911. doi: 10.1371/journal.pntd.0000911
Nalpas, N. C., Park, S. D., Magee, D. A., Taraktsoglou, M., Browne, J. A., Conlon, K. M., et al. (2013). Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro. BMC Genomics 14:230. doi: 10.1186/1471-2164-14-230
Oliveira, M. S., Fraga, A. G., Torrado, E., Castro, A. G., Pereira, J. P., Filho, A. L., et al. (2005). Infection with Mycobacterium ulcerans induces persistent inflammatory responses in mice. Infect. Immun. 73, 6299–6310. doi: 10.1128/IAI.73.10.6299-6310.2005
Portaels, F., Elsen, P., Guimaraes-Peres, A., Fonteyne, P. A., and Meyers, W. M. (1999). Insects in the transmission of Mycobacterium ulcerans infection. Lancet 353, 986. doi: 10.1016/S0140-6736(98)05177-0
Rienksma, R. A., Suarez-Diez, M., Mollenkopf, H. J., Dolganov, G. M., Dorhoi, A., Schoolnik, G. K., et al. (2015). Comprehensive insights into transcriptional adaptation of intracellular mycobacteria by microbe-enriched dual RNA sequencing. BMC Genomics 16:34. doi: 10.1186/s12864-014-1197-2
Sanhueza, D., Chevillon, C., Colwell, R., Babonneau, J., Marion, E., Marsollier, L., et al. (2016). Chitin promotes Mycobacterium ulcerans growth. FEMS Microbiol. Ecol. 92:fiw067. doi: 10.1093/femsec/fiw067
Schnappinger, D., Ehrt, S., Voskuil, M. I., Liu, Y., Mangan, J. A., Monahan, I. M., et al. (2003). Transcriptional adaptation of Mycobacterium tuberculosis within macrophages: insights into the phagosomal environment. J. Exp. Med. 198, 693–704. doi: 10.1084/jem.20030846
Schroeder, A., Mueller, O., Stocker, S., Salowsky, R., Leiber, M., Gassmann, M., et al. (2006). The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol. Biol. 7:3. doi: 10.1186/1471-2199-7-3
Silva, M. T., Portaels, F., and Pedrosa, J. (2009). Pathogenetic mechanisms of the intracellular parasite Mycobacterium ulcerans leading to Buruli ulcer. Lancet Infect. Dis. 9, 699–710. doi: 10.1016/S1473-3099(09)70234-8
Szafranska, A. K., Oxley, A. P., Chaves-Moreno, D., Horst, S. A., Rosslenbroich, S., Peters, G., et al. (2014). High-resolution transcriptomic analysis of the adaptive response of Staphylococcus aureus during acute and chronic phases of osteomyelitis. MBio 5:e01775–14. doi: 10.1128/mBio.01775-14
Talaat, A. M., Lyons, R., Howard, S. T., and Johnston, S. A. (2004). The temporal expression profile of Mycobacterium tuberculosis infection in mice. Proc. Natl. Acad. Sci. U.S.A. 101, 4602–4607. doi: 10.1073/pnas.0306023101
Talaat, A. M., Ward, S. K., Wu, C. W., Rondon, E., Tavano, C., Bannantine, J. P., et al. (2007). Mycobacterial bacilli are metabolically active during chronic tuberculosis in murine lungs: insights from genome-wide transcriptional profiling. J. Bacteriol. 189, 4265–4274. doi: 10.1128/JB.00011-07
Torrado, E., Fraga, A. G., Castro, A. G., Stragier, P., Meyers, W. M., Portaels, F., et al. (2007). Evidence for an intramacrophage growth phase of Mycobacterium ulcerans. Infect. Immun. 75, 977–987. doi: 10.1128/IAI.00889-06
Torrado, E., Fraga, A. G., Logarinho, E., Martins, T. G., Carmona, J. A., Gama, J. B., et al. (2010). IFN-gamma-dependent activation of macrophages during experimental infections by Mycobacterium ulcerans is impaired by the toxin mycolactone. J. Immunol. 184, 947–955. doi: 10.4049/jimmunol.0902717
Tuanyok, A., Tom, M., Dunbar, J., and Woods, D. E. (2006). Genome-wide expression analysis of Burkholderia pseudomallei infection in a hamster model of acute melioidosis. Infect. Immun. 74, 5465–5476. doi: 10.1128/IAI.00737-06
Vincent, Q. B., Ardant, M. F., Adeye, A., Goundote, A., Saint-Andre, J. P., Cottin, J., et al. (2014). Clinical epidemiology of laboratory-confirmed Buruli ulcer in Benin: a cohort study. Lancet Glob. Health 2, e422–e430. doi: 10.1016/S2214-109X(14)70223-2
Zogo, B., Djenontin, A., Carolan, K., Babonneau, J., Guegan, J. F., Eyangoh, S., et al. (2015). A field study in benin to investigate the role of mosquitoes and other flying insects in the ecology of Mycobacterium ulcerans. PLoS Negl. Trop. Dis. 9:e0003941. doi: 10.1371/journal.pntd.0003941
Keywords: RNA purification, RNA-seq, Mycobacterium ulcerans, cross-talk, host-bacteria interaction
Citation: Robbe-Saule M, Babonneau J, Sismeiro O, Marsollier L and Marion E (2017) An Optimized Method for Extracting Bacterial RNA from Mouse Skin Tissue Colonized by Mycobacterium ulcerans. Front. Microbiol. 8:512. doi: 10.3389/fmicb.2017.00512
Received: 12 December 2016; Accepted: 13 March 2017;
Published: 24 March 2017.
Edited by:Dongsheng Zhou, Beijing Institute of Microbiology and Epidemiology, China
Reviewed by:Jorge Pedrosa, University of Minho, Portugal
Tianyu Zhang, Guangzhou Institutes of Biomedicine and Health (CAS), China
Copyright © 2017 Robbe-Saule, Babonneau, Sismeiro, Marsollier and Marion. 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) or licensor 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.
*Correspondence: Estelle Marion, firstname.lastname@example.org