ORIGINAL RESEARCH article

Front. Microbiol., 28 February 2020

Sec. Physiology and Metabolism of Microorganisms

Volume 11 - 2020 | https://doi.org/10.3389/fmicb.2020.00282

Regulation of Microbial Metabolic Rates Using CRISPR Interference With Expanded PAM Sequences

  • Department of Systems Biotechnology, Chung-Ang University, Anseong, South Korea

Abstract

Genome-editing CRISPR/Cas9 technology has led to the development of artificial transcriptional repressors, also known as CRISPR interference (CRISPRi). The deactivated Cas9 (dCas9) protein guided by crRNA can specifically bind to target DNA sequences, including promoters and operators, without DNA cleavage. Protospacer adjacent motif (PAM) sequence dependence may be disadvantageous in the design of target-specific CRISPRi, as the PAM sequence is essential for DNA cleavage by the CRISPR/Cas9 system. We constructed a chromosomally integrated dCas9 system (ΔaraBAD:dcas9) in Escherichia coli under the control of the L-arabinose-inducible PBAD promoter. Plasmids carrying various crRNAs with target sequences specific for the gal promoter (−10 region), and the galETK structural genes in the gal operon, were transformed into dCas9-expressing E. coli. Cellular growth and/or galactose metabolic rates were monitored in the presence or absence of gratuitous L-arabinose. D-galactose consumption and cell growth rates were partially retarded by targeting transcriptional elongation but were fully inhibited by targeting transcriptional initiation. Moreover, RT-qPCR analysis showed that CRISPRi with several modified PAM sequences can repress the transcription of target DNAs. These results indicate that cellular metabolic rates and cell growth can be controlled by targeting structural genes or regulatory regions using CRISPRi; also, a loose PAM sequence dependence can expand the DNA targets of CRISPRi.

Introduction

The biological role of CRISPR has been identified as a prokaryotic adaptive immune system (Mojica et al., 2005; Pourcel et al., 2005). Since then, the application of CRISPR has been expanding explosively in the field of biotechnology (Doudna and Charpentier, 2014), especially because the components of CRISPR nucleases (i.e. crRNAs and Cas proteins) are separated functionally in the recognition of target nucleotide sequences and breakage of phosphodiester bonds in nucleic acids (Dagdas et al., 2017). Moreover, a specific protospacer adjacent motif (PAM) (∼3–5 bp) sequence must be located near the target sequence for the proper operation of the CRISPR/Cas system (Rath et al., 2015). In addition, the 5′-NGG sequence should be located right after target DNA sequences in the case of CRISPR/Cas9 (O’connell et al., 2014).

The CRISPR/Cas9 nuclease has been widely used in gene knockout studies in eukaryotic cells by the highly specific digestion of target DNA sequences followed by non-homologous end joining, resulting in the insertion or deletion of target genes (Zhang et al., 2014b). In addition, site-specific genome editing can be achieved by the introduction of single-stranded oligonucleotides and/or double-stranded DNA fragments flanking homologous regions into prokaryotic or eukaryotic cells, respectively, which can be negatively selected by the CRISPR/Cas9 nuclease (Knott and Doudna, 2018).

In addition to genome-editing tools, deactivated or dead Cas9 (dCas9) was created to be applied to artificial gene regulators (Larson et al., 2013). dCas9 has inactivated HNH and RuvC domains responsible for nuclease activity. It has been known that the target DNA sequence cannot be digested but is occupied by dCas9-crRNA complexes, which can affect the expression of target genes, and this is called CRISPR interference (CRISPRi) (Marraffini and Sontheimer, 2010a). When the dCas9-crRNA complex binds to the promoter of target genes, the initiation of transcription can be inhibited. If the complex binds to the coding region of the target gene, the elongation of transcription can be blocked (Qi et al., 2013). In prokaryotic cells, gene regulation is hardly achieved by RNA interference. Therefore, CRISPRi is very helpful in the engineering of microbial cells (Qiu et al., 2018; Xu and Qi, 2019).

Since the bacterial adaptive immune system, CRISPR, has evolved to distinguish self and non-self target sequences (Marraffini and Sontheimer, 2010b), the presence of the PAM sequence near the target DNA is essential for the DNA double-strand breakage generated by the Cas9-crRNA complex (Anders et al., 2014). Such PAM dependence can narrow down the possible target sequences in the Cas9-crRNA complex. However, only binding to the target is necessary for the purpose of dCas9-crRNA, which suggests that PAM dependence is lower than Cas9-crRNA dependence.

In this study, we evaluated how strongly the dCas9-crRNA complex inhibits initiation and elongation in the transcription of the gal operon. Also, we tested whether dCas9-crRNA can recognize DNA targets with modified PAM sequences (NNG and NGN). Our results demonstrate that cellular metabolic rates can be controlled with CRISPRi and that a loose PAM sequence dependence can expand the DNA targets of dCas9-crRNA.

Materials and Methods

Strains and Culture Conditions

Escherichia coli K-12 DH5α and MG1655 strains were used as a cloning host and a dCas9 expression host, respectively (Table 1). E. coli cells were grown routinely in LB medium. When needed, antibiotics such as ampicillin (50 μg/ml), kanamycin (25 μg/ml), and spectinomycin (75 μg/ml) were used.

TABLE 1

NameCharacteristicsSource/reference
Strain
DH5αfhuA2 lacΔU169 phoA glnV44 Φ80lacZΔM15 gyrA96 recA1 relA1 endA1 thi-1 hsdR17Laboratory stock
MG1655FilvG rfb-50 rph-1S. Adhya
HK1060MG1655, ΔaraBAD:PBAD-dcas9-KmRThis study
Plasmid
pKD46pSC101orits, araC, λ red genes, AmpRDatsenko and Wanner (2000)
pdCas9P15A ori, dcas9 gene, CmRAddgene #46569
pTargetFpBR322 ori, crRNA scaffold, SpRAddgene #62226
pHK459PAM(TGG), pBR322 ori, crRNA target 498AGGCTGTAACTGCGGGATCA517 in galK), SpRThis study
pBJ003PAM(CGG), pBR322 ori, crRNA target 573CGAACAGAAATCACCAATGC592 in galT), SpRThis study
pBJ004PAM(TGG), pBR322 ori, crRNA target 406TACGTTGAAAGCTTCCCGAC425 in galE), SpRThis study
pBJ005PAM(TGG), pBR322 ori, crRNA target –30TTCGCATCTTTGTTATGCTA–11 in Pgal), SpRThis study
pBJ021PAM(CAG), pBR322 ori, crRNA target 503GTAACTGCGGGATCATGGAT522 in galK), SpRThis study
pBJ022PAM(GCG), pBR322 ori, crRNA target 489CCAGTTTGTAGGCTGTAACT508 in galK), SpRThis study
pBJ023PAM(ATG), pBR322 ori, crRNA target 497TAGGCTGTAACTGCGGGATC516 in galK), SpRThis study
pBJ024PAM(GGA), pBR322 ori, crRNA target 499GGCTGTAACTGCGGGATCAT518 in galK), SpRThis study
pBJ025PAM(AGC), pBR322 ori, crRNA target 504TAACTGCGGGATCATGGATC523 in galK), SpRThis study
pBJ026PAM(TGT), pBR322 ori, crRNA target 482CAGAAAACCAGTTTGTAGGC501 in galK), SpRThis study
pBJ027PAM(AAG), pBR322 ori, crRNA target –1GAAATAACCATAGCATAACA–20 in Pgal), SpRThis study
pBJ028PAM(GCG), pBR322 ori, crRNA target –6AACCATAGCATAACAAAGAT–25 in Pgal), SpRThis study
pBJ029PAM(ATG), pBR322 ori, crRNA target –31TTTCGCATCTTTGTTATGCT–12 in Pgal), SpRThis study
pBJ030PAM(AGA), pBR322 ori, crRNA target –2AAATAACCATAGCATAACAA–21 in Pgal), SpRThis study
pBJ031PAM(TGC), pBR322 ori, crRNA target –5TAACCATAGCATAACAAAGA–24 in Pgal), SpRThis study
pBJ032PAM(GGT), pBR322 ori, crRNA target –29TCGCATCTTTGTTATGCTAT–10 in Pgal), SpRThis study

Strains and plasmids used in this study.

Genomic Integration

The primers used in this study are listed in Table 2. In order to make a dcas9-KmR cassette, the dcas9 gene was amplified using plasmid pdCas9 [a gift from Luciano Marraffini, Addgene plasmid #46569 (Bikard et al., 2013)] as a template and fused with a kanamycin resistance gene by overlap PCR. The PCR products of the dcas9-KmR cassette with a homologous region for recombineering at the araBAD promoter in the chromosome were electroporated into E. coli K-12 MG1655 cells harboring plasmid pKD46 (Datsenko and Wanner, 2000) to generate E. coli HK1060 (ΔaraBAD:PBAD-dcas9-KmR). Then, plasmid pKD46 was cured at 42°C. The disruption of the araBAD genes and the integration of the dcas9 gene in HK1060 cells were confirmed by PCR.

TABLE 2

NameSequence (5′→3′)Description
Cas-AraC-F1GCTTTTTATCGCAACTCTCTACTGTTTCTCCATACCCGTTTTTTTGGConstruction of cas9-KmR cassette
ATGGAGTGAAACGATGGATAAGAAATACTCAATAGGCT
KmR-AraD-R4GTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTConstruction of cas9-KmR cassette
GGCGTGGCGCGAGTGTAGGCTGGAGCTGCTTCGAAGTT
Cas-KmR-OR2GACGGATCCCCGGAATTCAGTCACCTCCTAGCTGACTCAAAConstruction of cas9-KmR cassette
Cas-KmR-OF3GCTAGGAGGTGACTGAATTCCGGGGATCCGTCGACCTGCAGConstruction of cas9-KmR cassette
Ara_F1AGCAGCTCCGAATAGCGCCCTTCCCCTTGCConfirmation of ΔaraBAD
AraD_500dnCCAGCCAGAAGGAGACTTCTGTCCCTTGConfirmation of ΔaraBAD
Sm_ATG_RGATACTGGGCCGGCAGGCGCTCCATTGCCCConstruction of crRNA plasmids
Sm_TAA_FGCAATGGAGCGCCTGCCGGCCCAGTATCAGConstruction of crRNA plasmids
pTarget_seqAACGCCTGGTATCTTTATAGTCCTGTCGSequencing of crRNA plasmids
Cas-galK-WFAGGCTGTAACTGCGGGATCAGTTTTAGAGCTAGAAATAGCAAGConstruction of pHK459 plasmid
Cas-galK-WRTGATCCCGCAGTTACAGCCTACTAGTATTATACCTAGGACTGConstruction of pHK459 plasmid
galT_Target_FCGAACAGAAATCACCAATGCGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ003 plasmid
galT_Target_RGCATTGGTGATTTCTGTTCGACTAGTATTATACCTAGGACTGConstruction of pBJ003 plasmid
galE Target_FTACGTTGAAAGCTTCCCGACGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ004 plasmid
galE_Target_RGTCGGGAAGCTTTCAACGTAACTAGTATTATACCTAGGACTGConstruction of pBJ004 plasmid
Pgal Target_FTTCGCATCTTTGTTATGCTAGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ005 plasmid
Pgal Target_RTAGCATAACAAAGATGCGAAACTAGTATTATACCTAGGACTGConstruction of pBJ005 plasmid
galK_CAG_FGTAACTGCGGGATCATGGATGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ021 plasmid
galK_CAG_RATCCATGATCCCGCAGTTACACTAGTATTATACCTAGGACTGConstruction of pBJ021 plasmid
galK_GCG_FCCAGTTTGTAGGCTGTAACTGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ022 plasmid
galK_GCG_RAGTTACAGCCTACAAACTGGACTAGTATTATACCTAGGACTGConstruction of pBJ022 plasmid
galK_ATG_FTAGGCTGTAACTGCGGGATCGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ023 plasmid
galK_ATG_RGATCCCGCAGTTACAGCCTAACTAGTATTATACCTAGGACTGConstruction of pBJ023 plasmid
galK_GGA_FGGCTGTAACTGCGGGATCATGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ024 plasmid
galK_GGA_RATGATCCCGCAGTTACAGCCACTAGTATTATACCTAGGACTGConstruction of pBJ024 plasmid
galK_AGC_FTAACTGCGGGATCATGGATCGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ025 plasmid
galK_AGC_RGATCCATGATCCCGCAGTTAACTAGTATTATACCTAGGACTGConstruction of pBJ025 plasmid
galK_TGT_FCAGAAAACCAGTTTGTAGGCGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ026 plasmid
galK_TGT_RGCCTACAAACTGGTTTTCTGACTAGTATTATACCTAGGACTGConstruction of pBJ026 plasmid
Pgal_AAG_FGAAATAACCATAGCATAACAGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ027 plasmid
Pgal_AAG_RTGTTATGCTATGGTTATTTCACTAGTATTATACCTAGGACTGConstruction of pBJ027 plasmid
Pgal_GCG_FAACCATAGCATAACAAAGATGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ028 plasmid
Pgal_GCG_RATCTTTGTTATGCTATGGTTACTAGTATTATACCTAGGACTGConstruction of pBJ028 plasmid
Pgal_ATG_FTTTCGCATCTTTGTTATGCTGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ029 plasmid
Pgal_ATG_RAGCATAACAAAGATGCGAAAACTAGTATTATACCTAGGACTGConstruction of pBJ029 plasmid
Pgal_AGA_FAAATAACCATAGCATAACAAGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ030 plasmid
Pgal_AGA_RTTGTTATGCTATGGTTATTTACTAGTATTATACCTAGGACTGConstruction of pBJ030 plasmid
Pgal_TGC_FTAACCATAGCATAACAAAGAGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ031 plasmid
Pgal_TGC_RTCTTTGTTATGCTATGGTTAACTAGTATTATACCTAGGACTGConstruction of pBJ031 plasmid
Pgal_GGT_FTCGCATCTTTGTTATGCTATGTTTTAGAGCTAGAAATAGCAAGConstruction of pBJ032 plasmid
Pgal_GGT_RATAGCATAACAAAGATGCGAACTAGTATTATACCTAGGACTGConstruction of pBJ032 plasmid
galE-RTLCCACCGTTTATGGCGATCAGRT-qPCR for transcripts of galE gene
galE-RTRGTTCCACCATCAGCTTGCTTRT-qPCR for transcripts of galE gene
galK-RTLGCCTCTATGCGCGATGATTTRT-qPCR for transcripts of galK gene
galK-RTRCACCTTTGTCGCCAATCACART-qPCR for transcripts of galK gene
16S-RTLCAGCAGCCGCGGTAATACRT-qPCR for 16S rRNA
16S-RTRACCAGGGTATCTAATCCTGTRT-qPCR for 16S rRNA

Primers used in this study.

Plasmid Construction

To generate the crRNA expression plasmids, pTargetF (a gift from Sheng Yang, Addgene plasmid #62226) was used as a template to amplify 0.9 and 1.2-kb fragments with primer pairs of Cas-galK-WF and Sm_ATG_R and Sm_TAA_F and Cas-galK-WR, respectively. These two fragments were purified for isothermal assembly using Gibson Assembly® Master Mix (NEB, MA, United States) to make pHK459 plasmid (Table 1). Other crRNA plasmids carrying the spectinomycin resistance gene as a selection marker were constructed by Gibson Assembly in the same manner. Target sequences recognized by crRNAs were designed in the DNA sequences (i.e. the gal promoter and the galE, galT, and galK genes) in the gal operon. The extended −10 region was selected as the DNA target in the gal promoter. For the experiments with modified PAM sequences, we constructed additional crRNA plasmids targeting the gal promoter and the galK gene. All the crRNA plasmids were confirmed by Sanger sequencing and transformed into HK1060 cells for the growth assay.

Growth Assay

HK1060 cells carrying each crRNA plasmid were grown in LB broth containing spectinomycin at 37°C overnight as starter cultures. Next, 1.0% of cells were inoculated in 25 ml of M9 minimal medium containing sodium succinate (final 20 mM) as a carbon source and spectinomycin in 250-ml flasks. crRNAs were expressed constitutively, and gratuitous L-arabinose (final 1 mM) was added for the expression of dCas9 proteins to form the dCas9-crRNA complex in HK1060 cells carrying crRNA plasmids. To determine whether the transcription of the gal operon could be negatively regulated by the dCas9-crRNA complex, D-galactose (final 20 mM) was added as an additional carbon source 6 h after the beginning of flask cultures. Optical density at 600 nm was measured every 3 h to monitor cellular growth using an Ultrospec 8000 spectrophotometer (GE Healthcare, Uppsala, Sweden).

Transcript Analysis

Expression levels of galE and galK genes were monitored using RT-qPCR. HK1060 cells carrying each crRNA plasmid were grown M9 minimal medium containing sodium succinate (final 20 mM) as a carbon source with or without L-arabinose (final 1 mM). D-galactose (final 20 mM) was added 6 h after the beginning of flask cultures. RNAs were isolated using the RNeasy® Mini kit (Cat. No. 74104; Qiagen, Hilden, Germany) 3 h after the addition of D-galactose. PCR primer sequences for the target genes were designed at the Universal Probe Library Assay Design Center1 and listed in Table 2. Quantitative real-time PCR (RT-qPCR) reactions were carried out on a LightCycler 96 (Roche Diagnostics, Mannheim, Germany) using the RealHelixTM qPCR kit (Nanohelix, South Korea). Five nanograms of each total RNA were used in RT-qPCR reactions under the following conditions: cDNA synthesis (50°C, 40 min); denaturation (95°C, 12 min); amplification for 50 cycles (95°C, 20 s; 60°C, 1 min). The raw fluorescence data were normalized against the expression level of 16S ribosomal RNA. In order to calculate the relative abundance of galE (374–453 base region from +1 start codon) and galK (895–963 base region) transcripts, mRNA levels of galE and galK genes in the presence of L-arabinose were divided by mRNA levels of the corresponding genes in the absence of L-arabinose.

Metabolite Analysis

The concentrations of D-galactose in the culture medium were determined by high-performance liquid chromatography (HPLC, RID-10A RI monitor, Shimadzu, Japan) with an Aminex HPX-87H column (300 × 7.8 mm, Hercules, BioRad), as described previously (Kim et al., 2013). After centrifugation of the cell culture broth, the supernatant was filtered through a 0.2-μm syringe filter. The column was isocratically eluted at 47°C with a flow rate of 0.5 ml min–1 using 0.01 N H2SO4.

Results

Repression of gal Operon by CRISPRi

Cellular growth and D-galactose utilization were monitored to assess how well CRISPRi interferes with transcription of the gal operon. Since CRISPRi can inhibit the initiation and elongation of transcription, the gal promoter and the gal structural genes (galE, galT, and galK) were used as the DNA targets of the dCas9-crRNA complex (Figure 1A). Cells carrying the L-arabinose-inducible dcas9 gene were transformed with various crRNA plasmids. The transformed cells were grown in succinate minimal media in the presence or absence of L-arabinose, which can serve as a gratuitous inducer. Succinate was used as a carbon source to prevent carbon catabolite repression. At the early exponential phase (∼OD600 of 0.3), D-galactose was added to the cell culture at 6 h after inoculation. If the gal transcription was not inhibited by dCas9-crRNAs, the increase of cellular growth and the consumption of D-galactose could be observed.

FIGURE 1

HK1060 cells were transformed with crRNA plasmids expressing crRNAs recognizing the DNA targets with the original PAM sequence (5′-NGG). As a result, in the absence of L-arabinose, all the transformant cells reached an OD600 of 4, and consumed D-galactose (total 20 mM) within 18 h (Figures 1B,C). However, cell growth rates were decreased in the presence of non-metabolizable L-arabinose. This means L-arabinose-induced dCas9 can combine with crRNA to make dCas9-crRNA, which can interfere with the transcription of D-galactose metabolic genes.

In the case of the gal promoter, D-galactose was not consumed at all until the end of the experiment (Figure 1C). The gal promoter is composed of two overlapped extended −10 promoters (Lewis et al., 2015), which may be covered thoroughly by dCas9-crRNA complex. These results indicate that transcriptional initiation can be completely inhibited by CRISPRi.

In the case of structural genes as the inhibition targets, the cell growth and D-galactose consumption were differentially retarded. This means that the expression of gal genes was repressed partially by the blocking of transcriptional elongation. In the profiles of HPLC, the level of gratuitous L-arabinose was not changed in any of the samples, because the ara operon was disrupted by the integration of the dcas9 gene in HK1060 cells. It was confirmed that our model CRISPRi system operated fully and partially in the promoter and structural genes of D-galactose metabolism, respectively.

Repression of gal Promoter by CRISPRi With Expanded PAM Sequences

In order to determine whether CRISPRi follows the PAM sequence rule (5′-NGG) of CRISPR/Cas9, various crRNA plasmids were designed, as shown in Table 1. Modified PAM sequences (5′-NGN and 5′-NNG) were tested in the gal promoter as the DNA target of the dCas9-crRNA complex (Figures 2A,B).

FIGURE 2

In addition to in NGG (TGG), the growth retardation was observed in the case of NAG (AAG), NTG (ATG), and NGA (AGA) (Figure 2C). Moreover, D-galactose consumption was not observed with the expanded PAM sequences until the end of the experiments (Figure 2D). The D-galactose consumption rates (between 9 and 12 h in graphs) in the cases of NCG (GCG), NGC (TGC), and NGT (GGT) were not changed in the presence (dCas9 expressed) or absence (dCas9 not expressed) of L-arabinose (Figure 2E). Transcript analysis showed that the differential repression of the gal promoter by dCas9-crRNA complexes with the modified PAM sequences. The order of repression was as follows: NGG (TGG) > NAG (AAG) > NGA (AGA) > NTG (ATG) (Figure 3A). In the cases of NCG (GCG), NGC (TGC), and NGT (GGT), the repression of gal promoter could not be observed in the RT-qPCR experiment. These data showed that even weakly attached dCas9-crRNA complex with modified PAM sequences did not result in sufficient gal transcription for galactose metabolism.

FIGURE 3

Exploration of galK Structural Gene as the Target of CRISPRi With Expanded PAM Sequences

The same set of modified PAM sequences (5′-NGN and 5′-NNG) was tested in the galK structural gene as the DNA target of the dCas9-crRNA complex (Figures 4A,B). As a result, cell growth was affected slightly in the cases of NGG (TGG), NTG (ATG), and NGA (GGA) as PAM sequences (Figure 4C). According to the D-galactose consumption rate, the order of interference was as follows: NGG (TGG) > NGA (GGA) > NTG (ATG) (Figures 4D,E). The D-galactose consumption rate (between 9 and 12 h in graphs) in the case of ATG was decreased from ∼3 mM/h (dCas9 not expressed due to the absence of L-arabinose) to ∼1 mM/h (dCas9 expressed by L-arabinose).

FIGURE 4

We analyzed the abundance of mRNA of the galK gene in the presence or absence of dCas9-crRNA complexes (Figure 3B). The galK transcripts in downstream region of dCas9-crRNA complex-binding sites were monitored by RT-qPCR. In case of NGG (TGG) as a PAM sequence, the transcriptional elongation of the galK gene was affected by CRISPRi. However, any significant interference in the transcriptional elongation was not clearly observed in other PAM sequences using RT-qPCR analysis.

Discussion

CRISPR interference can regulate the expression of DNA targets by inhibiting the transcription process. Designed target recognition sequences of modular crRNAs can be used to bind any DNA target to dCas9-crRNA complexes (Hawkins et al., 2015). While nonsense point mutation can affect cells permanently by deleting the function of the target gene, dCas9-mediated gene expression control allows the selection of desired time points (Quebatte and Dehio, 2017). Endogenous competing pathway genes can be repressed by CRISPRi in E. coli for metabolic engineering (Kim et al., 2017). Even if the genes are essential for growth, metabolic engineering with CRISPRi provides the necessary metabolites without disturbing cell growth (Cress et al., 2017).

Gene regulation using CRISPRi is carried out in two ways. First, binding of the dCas9-crRNA complex in the promoter region can fully repress transcriptional initiation. Also, elongation of transcription can be blocked by positioning the dCas9-crRNA complex in the middle of the open reading frame of the target gene (Cao et al., 2017). In our study, the gal operon and galactose metabolism were used as a model to demonstrate how dCas9-crRNA can precisely regulate the expression of target genes (gal promoter and structural galE, galT, and galK genes) and the corresponding metabolism (Figure 1A).

Non-metabolizable L-arabinose-induced dCas9 proteins can combine with target-specific crRNAs to make the dCas9-crRNA complex, which can bind to the DNA targets in the galactose promoter, preventing cellular consumption of D-galactose (Figure 1D). When the gal structural genes (galE, galT, and galK) were designed as the DNA targets of the dCas9-crRNA complex, we observed differential D-galactose consumption rates. RT-qPCR analysis showed that the initiation of galE transcription was significantly repressed, and the elongation of galK transcription was slightly blocked the dCas9-crRNAs with TGG as PAM sequence (Figure 3). These results might be explained by that RNA polymerase cannot bind to the gal promoter occupied by the dCas9-crRNA complex, but RNA polymerase can pass through some of the roadblocks of dCas9-crRNA during transcriptional elongation. Although all enzymes encoded by galE, galT, and galK genes are essential in galactose metabolism, the most retarded cell growth was observed in galE-targeted cells (Figure 1B), presumably because galE mutant cells cannot grow well due to the depletion of pyrimidine nucleotides (Lee et al., 2009).

The PAM sequence should be considered in CRISPR technologies when the target DNA sequence is designed. The biological role of the PAM sequence is to discriminate between non-self DNA from outside cells and self DNA stored in the chromosome in bacterial adaptive immune system. It has been known that the PAM sequence is indispensable for DNA cleavage by the Cas9-crRNA complex (Mojica et al., 2009). It has been reported that the hydrogen bonds between nucleotides in the PAM sequence (5′-NGG) and the two arginine residues (Arg1333 and Arg1335) in the Cas9 protein are important in DNA cleavage (Anders et al., 2014). The PAM dependence in the CRISPR/Cas9 system has been studied for expanding the target sequence and reducing off-target activity. Reportedly, the 5′-NGA sequence can function as PAM, although the modified PAM sequence exhibits less gfp gene cleavage activity (Zhang et al., 2014a). It is still unclear whether the PAM sequence is essential in the recognition of target sequences in CRISPR/dCas9 in bacterial cells.

The PAM dependence in CRISPRi was investigated using regulation of the gal promoter (Figure 2) and the galK gene (Figure 4) by dCas9-crRNA complex. The differential repressions of gal promoter by modified PAM sequences were observed by RT-qPCR analysis (Figure 3A). The order of gal promoter repression showed that target binding affinities of dCas9-crRNAs with modified PAM sequences might be weaker than those with the original PAM sequence (NGG). However, the CRISPRi-induced growth retardation and totally impaired D-galactose consumption suggest that NAG, NTG, and NGA can function as the original PAM sequence (NGG) in the gal transcription initiation (Figure 2). This means that the gal promoter even weakly occupied by dCas9-crRNA complex cannot be properly recognized by the RNA polymerase in the cells.

In the case of the galK structural gene as the DNA target of CRISPRi, transcriptional elongation can be clearly modulated by CRISPRi with the original PAM sequence (NGG) based on the results of cellular growth, metabolic rates, and transcript abundance. Since the original PAM sequence occurs frequently in the open reading frame of common genes, the expression of essential genes can be attenuated by CRISPRi with the original PAM sequence.

These results showed that the DNA target sequence for CRISPRi can be expanded, allowing the design of more specific target sequences such as narrow promoter regions. Our results also show that cellular growth and metabolic rates can be controlled by dCas9-crRNAs. The information about CRISPRi with expanded PAM sequences may be useful in the design of DNA targets and crRNAs for fine gene regulation in microbial biotechnology.

Statements

Data availability statement

All datasets generated for this study are included in the article/supplementary material.

Author contributions

BK, HK, and SL designed the research, analyzed the data and wrote the manuscript. BK and HK performed the experiments.

Funding

This study was supported by the National Research Foundation of Korea (NRF-2019R1A4A1024764, 2018R1A6A3A11051083, 2017R1E1A1A01075124, and 2015R1A2A2A0105402).

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.

References

  • 1

    AndersC.NiewoehnerO.DuerstA.JinekM. (2014). Structural basis of PAM-dependent target DNA recognition by the Cas9 endonuclease.Nature513569573. 10.1038/nature13579

  • 2

    BikardD.JiangW.SamaiP.HochschildA.ZhangF.MarraffiniL. A. (2013). Programmable repression and activation of bacterial gene expression using an engineered CRISPR-Cas system.Nucleic Acids Res.4174297437. 10.1093/nar/gkt520

  • 3

    CaoY.LiX.LiF.SongH. (2017). CRISPRi-sRNA: transcriptional-translational regulation of extracellular electron transfer in Shewanella oneidensis.ACS Synth. Biol.616791690. 10.1021/acssynbio.6b00374

  • 4

    CressB. F.LeitzQ. D.KimD. C.AmoreT. D.SuzukiJ. Y.LinhardtR. J.et al (2017). CRISPRi-mediated metabolic engineering of E. coli for O-methylated anthocyanin production.Microb. Cell Fact.16:10. 10.1186/s12934-016-0623-3

  • 5

    DagdasY. S.ChenJ. S.SternbergS. H.DoudnaJ. A.YildizA. (2017). A conformational checkpoint between DNA binding and cleavage by CRISPR-Cas9.Sci. Adv.3:eaao0027. 10.1126/sciadv.aao0027

  • 6

    DatsenkoK. A.WannerB. L. (2000). One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products.Proc. Natl. Acad. Sci. U.S.A.9766406645. 10.1073/pnas.120163297

  • 7

    DoudnaJ. A.CharpentierE. (2014). Genome editing. The new frontier of genome engineering with CRISPR-Cas9.Science346;1258096. 10.1126/science.1258096

  • 8

    HawkinsJ. S.WongS.PetersJ. M.AlmeidaR.QiL. S. (2015). Targeted transcriptional repression in bacteria using CRISPR interference (CRISPRi).Methods Mol. Biol.1311349362. 10.1007/978-1-4939-2687-9_23

  • 9

    KimH. J.HouB. K.LeeS. G.KimJ. S.LeeD. W.LeeS. J. (2013). Genome-wide analysis of redox reactions reveals metabolic engineering targets for D-lactate overproduction in Escherichia coli.Metab. Eng.184452. 10.1016/j.ymben.2013.03.004

  • 10

    KimS. K.SeongW.HanG. H.LeeD. H.LeeS. G. (2017). CRISPR interference-guided multiplex repression of endogenous competing pathway genes for redirecting metabolic flux in Escherichia coli.Microb. Cell Fact.16:188. 10.1186/s12934-017-0802-x

  • 11

    KnottG. J.DoudnaJ. A. (2018). CRISPR-Cas guides the future of genetic engineering.Science361866869. 10.1126/science.aat5011

  • 12

    LarsonM. H.GilbertL. A.WangX.LimW. A.WeissmanJ. S.QiL. S. (2013). CRISPR interference (CRISPRi) for sequence-specific control of gene expression.Nat. Protoc.821802196. 10.1038/nprot.2013.132

  • 13

    LeeS. J.TrostelA.LeP.HarinarayananR.FitzgeraldP. C.AdhyaS. (2009). Cellular stress created by intermediary metabolite imbalances.Proc. Natl. Acad. Sci. U.S.A.1061951519520. 10.1073/pnas.0910586106

  • 14

    LewisD. E. A.LeP.AdhyaS. (2015). Differential role of base pairs on gal promoters strength.J. Mol. Biol.427792806. 10.1016/j.jmb.2014.12.010

  • 15

    MarraffiniL. A.SontheimerE. J. (2010a). CRISPR interference: RNA-directed adaptive immunity in bacteria and archaea.Nat. Rev. Genet.11181190. 10.1038/nrg2749

  • 16

    MarraffiniL. A.SontheimerE. J. (2010b). Self versus non-self discrimination during CRISPR RNA-directed immunity.Nature463568571. 10.1038/nature08703

  • 17

    MojicaF. J.Diez-VillasenorC.Garcia-MartinezJ.AlmendrosC. (2009). Short motif sequences determine the targets of the prokaryotic CRISPR defence system.Microbiology155733740. 10.1099/mic.0.023960-0

  • 18

    MojicaF. J.Diez-VillasenorC.Garcia-MartinezJ.SoriaE. (2005). Intervening sequences of regularly spaced prokaryotic repeats derive from foreign genetic elements.J. Mol. Evol.60174182. 10.1007/s00239-004-0046-3

  • 19

    O’connellM. R.OakesB. L.SternbergS. H.East-SeletskyA.KaplanM.DoudnaJ. A. (2014). Programmable RNA recognition and cleavage by CRISPR/Cas9.Nature516263266. 10.1038/nature13769

  • 20

    PourcelC.SalvignolG.VergnaudG. (2005). CRISPR elements in Yersinia pestis acquire new repeats by preferential uptake of bacteriophage DNA, and provide additional tools for evolutionary studies.Microbiology151653663. 10.1099/mic.0.27437-0

  • 21

    QiL. S.LarsonM. H.GilbertL. A.DoudnaJ. A.WeissmanJ. S.ArkinA. P. (2013). Repurposing CRISPR as an RNA-guided platform for sequence-specific control of gene expression.Cell15211731183. 10.1016/j.cell.2013.02.022

  • 22

    QiuH.GongJ.ButayeP.LuG.HuangK.ZhuG.et al (2018). CRISPR/Cas9/sgRNA-mediated targeted gene modification confirms the cause-effect relationship between gyrA mutation and quinolone resistance in Escherichia coli.FEMS Microbiol. Lett.365:fny127. 10.1093/femsle/fny127

  • 23

    QuebatteM.DehioC. (2017). Systems-level interference strategies to decipher host factors involved in bacterial pathogen interaction: from RNAi to CRISPRi.Curr. Opin. Microbiol.393441. 10.1016/j.mib.2017.08.002

  • 24

    RathD.AmlingerL.RathA.LundgrenM. (2015). The CRISPR-Cas immune system: biology, mechanisms and applications.Biochimie117119128. 10.1016/j.biochi.2015.03.025

  • 25

    XuX.QiL. S. (2019). A CRISPR-dCas toolbox for genetic engineering and synthetic biology.J. Mol. Biol.4313447. 10.1016/j.jmb.2018.06.037

  • 26

    ZhangY.GeX.YangF.ZhangL.ZhengJ.TanX.et al (2014a). Comparison of non-canonical PAMs for CRISPR/Cas9-mediated DNA cleavage in human cells.Sci. Rep.4:5405. 10.1038/srep05405

  • 27

    ZhangF.WenY.GuoX. (2014b). CRISPR/Cas9 for genome editing: progress, implications and challenges.Hum. Mol. Genet.23R40R46.

Summary

Keywords

metabolic regulation, PAM dependence, interference, dCas9, CRISPR

Citation

Kim B, Kim HJ and Lee SJ (2020) Regulation of Microbial Metabolic Rates Using CRISPR Interference With Expanded PAM Sequences. Front. Microbiol. 11:282. doi: 10.3389/fmicb.2020.00282

Received

07 October 2019

Accepted

07 February 2020

Published

28 February 2020

Volume

11 - 2020

Edited by

Weiwen Zhang, Tianjin University, China

Reviewed by

Soo Rin Kim, Kyungpook National University, South Korea; Yinhua Lu, Shanghai Normal University, China

Updates

Copyright

*Correspondence: Sang Jun Lee,

These authors have contributed equally to this work

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

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics