Difference between revisions of "Team:Bielefeld-CeBiTec/Project/unnatural base pair/preservation system"

(references)
 
(26 intermediate revisions by 5 users not shown)
Line 4: Line 4:
 
<body>
 
<body>
 
<div class="container">
 
<div class="container">
 +
<div id="title" style="background-image: url(https://static.igem.org/mediawiki/2017/5/5b/T--Bielefeld-CeBiTec--nanoporesequenzer.jpeg);">
 +
<img src="https://static.igem.org/mediawiki/2017/7/74/T--Bielefeld-CeBiTec--title-img-centrifuge.jpg">
 +
<div id="title-bg">
 +
<div id="title-text">
 +
Retention and Preservation System
 +
</div>
 +
</div>
 +
</div>
 +
 +
<div class="contentbox">
 +
<div class="content">
 +
<h3>Short Summary</h3>
 +
<article>
 +
The retention and preservation of unnatural nucleoside triphosphates and unnatural bases in the DNA requires a cellular retention and preservation system acting on multiple cellular levels. On one hand, degradation of the unnatural nucleoside triphosphates by enzymes like cytosine deaminase must be prohibited for efficient integration into the DNA. On the other hand, once the unnatural bases are incorporated into the DNA, loss or change of the unnatural bases have to be prevented. We want to achieve this by adapting CRISPR/Cas9 to detect and prevent mutations of the unnatural bases.
 +
</article>
 +
</div>
 +
</div>
  
 
<div class="contentbox">
 
<div class="contentbox">
Line 10: Line 27:
 
 
 
<!-- Ueberschriften -->
 
<!-- Ueberschriften -->
<h2> Preservation system </h2>
 
 
<h3> Preservation system using Cas9 </h3>
 
<h3> Preservation system using Cas9 </h3>
 
 
Line 16: Line 32:
 
<!-- Normaler Text -->
 
<!-- Normaler Text -->
 
<article>
 
<article>
Due to tautomerisation of isoG and hydrolysis of isoC<sup>m</sup> there is a need for a system, to preserve the unnatural base pair (UBP) on the plasmid. In 2017, Zhang et al . successfully deployed a CRISPR (<i>clustered regularly interspaced short palindromic repeat</i>)-Cas9 system for retention of a UBP. We adapted this conservation system to our UBP and thus used the bacterial immune response to eliminate all plasmid DNA that had lost the UBP.  
+
Due to tautomerisation of isoG and hydrolysis of isoC<sup>m</sup> and the resulting loss of the unnatural base pair (UBP), there is a need for a system, to preserve the UBP on the plasmid. In 2017, Zhang&nbsp;<i>et&nbsp;al</i>. successfully deployed a CRISPR (<i>clustered regularly interspaced short palindromic repeat</i>)-Cas9 system for retention of a UBP. We adapted this conservation system to our UBP and thus used CRISPR/Cas9 to eliminate all plasmid DNA that had lost the UBP.  
 
<br>
 
<br>
The nuclease Cas9 is part of the adaptive immune system of <i>Streptococcus pyogenes</i>, where it induces double strand breaks in the genomic DNA. This enzyme is recruited  by a CRISPR RNA (crRNA). A crRNA consists of direct repeats interspaced by variable sequences called protospacer. Those protospacers are derived from foreign DNA and encode the Cas9 guiding sequence (guide RNA). An auxiliary transactivating crRNA (tracrRNA) helps processing the precursor crRNA array into an active crRNA that contains the 20 nucleotide guide RNA. The guide RNA binds to the complementary genomic DNA sequence via Watson-Crick base pairing. For this binding, the genomic DNA sequence needs to be located upstream of a CRISPR type II specific 5’-NGG protospacer adjacent motif (PAM). To combine crRNA and tracrRNA a chimeric single stranded guide RNA (sgRNA) was designed. Therefore, only the 20 nucleotide guiding sequence needs to be exchanged for targeting any genomic sequence followed by a PAM sequence (Ran et al, 2013 a, b). A double strand break introduced by Cas9 leads to DNA degradation by exonucleases in prokaryotic cells (Simmon and Lederberg, 1972).
+
The nuclease Cas9 is part of the adaptive immune system of <i>Streptococcus pyogenes</i>, where it induces double strand breaks in the genomic DNA. This enzyme is recruited  by a CRISPR&nbsp;RNA&nbsp;(crRNA). A crRNA consists of direct repeats interspaced by variable sequences called protospacer. Those protospacers are derived from foreign DNA and encode the Cas9 guiding sequence&nbsp;(guide RNA). An auxiliary transactivating crRNA (tracrRNA) helps processing the precursor crRNA array into an active crRNA that contains the 20 nucleotide guide RNA. The guide RNA binds to the complementary genomic DNA sequence via Watson&#x2011;Crick base pairing. For this binding, the genomic DNA sequence needs to be located upstream of a CRISPR type II specific 5’&#x2001;NGG&nbsp;protospacer&nbsp;adjacent&nbsp;motif&nbsp;(PAM). Synthetically chimeric single stranded guide RNA (sgRNA) was designed by combining crRNA and tracrRNA. In the sgRNA, only the 20 nucleotide guiding sequence needs to be exchanged for targeting any genomic sequence followed by a PAM sequence (Ran&nbsp;<i>et al</i>., 2013 a, b). The resulting double strand break introduced by Cas9 leads to exonucleolytic degradation of the DNA in prokaryotic cells (Simmon and Lederberg, 1972).
 
<br>
 
<br>
The UBP isoG and isoC<sup>m</sup> is an orthogonal system. UBP inside the target DNA causes a mismatch to the sgRNA generally reducing the cleavage activity of Cas9 (Zhang et al., 2017). Accordingly, Cas9 can be programmed by sgRNAs to cleave all plasmids, which had lost the UBP due to point mutations (Figure 1). Consequently, this leads to a retention of the UBP in the plasmids.
+
In our case we envision a retention system, where Cas9 cleaves Plasmids at sites where the UBP is absent. This works by using a sgRNA complementary to the DNA sequence without the UBP. In plasmids with the UBP present, the mismatch between isoG/isoC<sup>m</sup> and sgRNA greatly decreases Cas9 activity (Zhang <i>et al.</i>, 2017). In the event of UBP loss, the sgRNA now binds perfectly to the mutated site and restores Cas9 activity which leads to degradation of the mutated plasmid. Consequently, this leads to a retention of the UBP in the plasmids.
 
+
 
<!-- Mittleres zentriertes Bild -->
 
<!-- Mittleres zentriertes Bild -->
 
<div class="figure large">
 
<div class="figure large">
<img class="figure image" src="https://static.igem.org/mediawiki/2017/3/35/T--Bielefeld-CeBiTec--Cas9_conservation_system_from_poster.png">
+
<img class="figure image" src="https://static.igem.org/mediawiki/2017/6/67/T--Bielefeld-CeBiTec--Cas9_conservation_system_from_poster2.png">
<p class="figure subtitle"><b>Figure 1: UBP conservation system using Cas9.</b><br> sgRNAs are targeted against every possible DNA sequence that had lost the UBP, which was incorporated on a plasmid. A: The UBP gets lost, which leads to a point mutation. One of the sgRNAs can bind to the DNA target sequence. Cas9 is recruited and cleaves the plasmid, which is followed by its degradation. B: Plasmids that contain a UBP in the DNA target sequence lead to a mismatch with every sgRNA. Cas9 does not cleave the plasmid, leading to a retention of the UBP.</p>
+
<p class="figure subtitle"><b>Figure 1: UBP conservation system using Cas9.</b><br> sgRNAs are targeted against every DNA sequence emerging from UBP loss on a plasmid. A: Loss of the UBP leads to a point mutation. Now a sgRNA can bind to the DNA target sequence. Cas9 is recruited and cleaves the plasmid, which is followed by its degradation. B: Plasmids that contain a UBP in the DNA target sequence lead to a mismatch with every sgRNA. Cas9 does not cleave the plasmid, leading to retention of the UBP.</p>
 
</div>
 
</div>
For this conservation system, we integrated the <i>cas9</i> gene into the genome of <i>E. coli</i> BL21(DE3). This new strain is called X1. To be able to strictly regulate the expression of <i>cas9</i> we used the optimized IPTG-inducible promoter <i>P<sub>lac-tight</sub></i>. The genomic knock in was done according to the protocol by Cobb et al., 2014 using the pCRISPomyces plasmid system. A sgRNA was designed searching for a reverse sequence with the constraint N(16)R(4)NGG (Cobb et al., 2004) inside the coding sequence of the arsB gene (Zhang et al., 2017) of the <i>E. coli</i> genome. <i>arsB</i> is coding for an arsenic efflux pump membrane protein. As a result, 5‘&#x2011;TATTGTTCATAATAGAAGAGAGG &#x2011; 3‘ turned out to be the suggested guiding sequence with the highest on-target activity score with being unique within the complete genome. Furthermore, a repair template is necessary for the knock in. Therefore, a composite BioBrick was assembled using 1 kb long flanking sequences, two times the terminator BBa_B0015 flanking the <i>P<sub>lac-tight</sub></i> and the <i>cas9</i> coding sequence (Figure 2).
+
Because the iGEM competition is based on plasmids that can be submitted to the parts registry we refrained from creating a novel strain for the UBP retention. Consequently, our retention system consists of two plasmids. The first plasmid (BBa_K2201010) contains the nucleotide transporter gene <i>PtNTT2</i> and the <i>cas9</i> on a pSB1K3 backbone. This plasmid was transformed into <i>E. coli</i> BL21(DE3). After producing chemically competent cells, another transformation was done with the second plasmid. The second plasmid carries the UBP and the sgRNAs for recruiting Cas9 on a pSB3C5 plasmid. Since two large plasmids produce a great metabolic stress for an organism, a more sufficient UBP retention could be achieved with a newly created <i>E. coli</i> strain, that has the genes <i>PtNTT2</i> and <i>cas9</i> in its genomic DNA. We provide the parts that are required to create a repair template for the knock in strategy using a CRISPR/Cas9 system. These parts were designed to perform a genomic integration into <i>E. coli</i> BL21(DE3). The genomic knock in can be done according to the protocol by Cobb <i>et al.</i>, 2014 using the pCRISPomyces plasmid system. A sgRNA was designed searching for a reverse sequence with the constraint N(16)R(4)NGG (Cobb <i>et al.</i>, 2004) inside the coding sequence of the <i>arsB</i> gene (Zhang <i>et al.</i>, 2017) of the <i>E. coli</i> genome. <i>arsB</i> is coding for an arsenic efflux pump membrane protein. As a result, 5‘ TATTGTTCATAATAGAAGAG 3‘ turned out to be the suggested guiding sequence with the highest on target activity score with being unique within the complete genome. The required repair template is designed with 1 kb long flanking sequences and two terminators BBa_B0015 flanking the gene of interest for example <i>cas9</i> (Figure 2). To be able to strictly regulate the expression of <i>cas9</i> we rationally designed an optimized IPTG inducible promoter <i>P<sub>lacO&#x2011;tight1</sub></i>.
 
 
 
<!-- Mittleres zentriertes Bild -->
 
<!-- Mittleres zentriertes Bild -->
 
<div class="figure large">
 
<div class="figure large">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/7/7c/T--Bielefeld-CeBiTec--repair_template_picture.png">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/7/7c/T--Bielefeld-CeBiTec--repair_template_picture.png">
<p class="figure subtitle"><b>Figure 2: The repair template for genomic integration into <i>E. coli</i> BL21(DE3).</b><br> The 1 kb left flanking sequence (LFS) and right flanking sequence (RFS) are taken from the genome of <i>E. coli</i> BL21(DE3). Inside the genome <i>LFS</i> and <i>RFS</i> are directly flanking the coding sequence of <i>arsB</i>. The strong terminators from BBa_B0015 were used to avO<sub>id</sub> basal expression and to stop transcription right after <i>cas9</i>. The coding sequence of <i>cas9</i> is taken from the pCRISPomyces plasmid system by Cobb et al., 2014 and is originated from <i>S. pyogenes</i>. It is negatively regulated by the IPTG-inducible promoter <i>P<sub>lac-tight</sub></i>.</p>
+
<p class="figure subtitle"><b>Figure 2: The repair template (BBa_K2201028) for genomic integration of <i>cas9</i> into <i>E.&nbsp;coli</i>&nbsp;BL21(DE3).</b><br> The 1&nbsp;kb left flanking sequence&nbsp;(LFS: BBa_K2201021) and right flanking sequence&nbsp;(RFS: K2201022) are taken from the genome of <i>E.&nbsp;coli</i>&nbsp;BL21(DE3). Inside the genome <i>LFS</i> and <i>RFS</i> are directly flanking the coding sequence of <i>arsB</i>. The strong terminators from BBa_B0015 were used to avoid basal expression and to stop transcription right after <i>cas9</i>. There are composite parts the consist of <i>LFS</i> + terminator + <i>P<sub>lacO&#x2011;tight1</sub></i> (BBa_K2201024) and terminator + <i>RFS</i> (BBa_K2201025) that can be assembled with any coding sequence of interest to create a repair template. The coding sequence of <i>cas9</i> is taken from the pCRISPomyces plasmid system by Cobb&nbsp;<i>et&nbsp;al</i>.,&nbsp;2014 and is originated from <i>S.&nbsp;pyogenes</i>. It is negatively regulated by the IPTG&#x2011;inducible promoter <i>P<sub>lacO&#x2011;tight1</sub></i>.</p>
 
</div>
 
</div>
This composite BioBrick from Figure 2 needs to be restricted with <i>Not</i>I to separate the linear repair template from the pSB1C3 backbone. Together with the target plasmid containing the sgRNA and the pCRISPomyces plasmid containing the cas9 all three elements were co-transfected into <i>E. coli</i> BL21(DE3). The genomic integration was verified by sequencing.
+
This composite BioBrick from Figure&nbsp;2 needs to be restricted with <i>Not</i>I to separate the linear repair template from the pSB1C3 backbone. Together with the target plasmid containing the sgRNA and the pCRISPomyces plasmid containing the <i>cas9</i> all three elements were co&#x2011;transfected into <i>E.&nbsp;coli</i>&nbsp;BL21(DE3). The genomic integration was verified by sequencing.
 
 
 
</article>
 
</article>
Line 42: Line 57:
 
<article>
 
<article>
 
 
We designed a <i>lac</i> operon for a tight repression called <i>P<sub>lac-tight</sub></i> to achieve a low transcription rate. The induction of the wild type lac operon increases the level of &#946; &#x2011; galactosidase 1000-fold (Müller et al., 1996). As part of the <i>lac</i> operon the <i>lac</i> repressor was the first repressor isolated and sequenced in 1966 by Gilbert and Müller-Hill. The wild type <i>lac</i> operon consists of the genes <i>lacZ</i>, <i>lacY</i> and <i>lacA</i> that are transcribed from the lac promoter <i>P<sub>lac</sub></i> into a polycistronic mRNA (Figure 4). These genes code for the proteins &#946; &#x2011; galactosidase, Lac permease and Lac transacetylase (Oehler et al., 1994).
+
We designed a tightly repressed <i>lac</i> operon called <i>P<sub>lacO&#x2011;tight1</sub></i> to achieve a low basal transcription rate. The <i>lac</i> operon was extensively researched starting from the midst of last century. This showed the role of the <i>lac</i> repressor and the inducible nature of the <i>lac</i> promoter, which results in 1000-fold increased expression in its activated state (Müller <i>et al.</i>, 1996). The wild type <i>lac</i> operon consists of the genes <i>lacZ</i>, <i>lacY</i> and <i>lacA</i> that are transcribed from the lac promoter <i>P<sub>lac</sub></i> into a polycistronic mRNA (Figure&nbsp;4). These genes code for the proteins &#946;&#x2011;galactosidase, Lac permease and Lac transacetylase (Oehler&nbsp;<i>et&nbsp;al</i>.,&nbsp;1994).
 
 
 
<!-- Mittleres zentriertes Bild -->
 
<!-- Mittleres zentriertes Bild -->
 
<div class="figure large">
 
<div class="figure large">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/c/c9/T--Bielefeld-CeBiTec--polycistronische_mRNA_lacO.png">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/c/c9/T--Bielefeld-CeBiTec--polycistronische_mRNA_lacO.png">
<p class="figure subtitle"><b>Figure 3: DNA map of the wild type <i>lac</i> operon and its transcriptional and translational products (Oehler et al., 1994).</b></p>
+
<p class="figure subtitle"><b>Figure&nbsp;3: DNA map of the wild type <i>lac</i> operon and its transcriptional and translational products (Oehler&nbsp;<i>et&nbsp;al</i>.,&nbsp;1994).</b></p>
 
</div>
 
</div>
 
 
The transcription is constitutively activated by the CAP protein. <i>lacI</i> codes for the tetrameric Lac repressor and is expressed by the promoter <i>P<sub>i</sub></i>. The Lac repressor can bind to the lac operators O1, O2 or O3. O3 is located 92 bp upstream of O1 and O2 401 bp downstream of O1. The inter-operator distances are counted from the center of O1 to the center of the distal operator. By binding simultaniously to O1 and O2 or to O1 and O3, the Lac repressor forms a DNA loop that negatively controlls the expression of <i>P<sub>lac</sub></i> (Oehler et al., 1994).  
+
The transcription is constitutively activated by the CAP&nbsp;protein. <i>lacI</i> codes for the tetrameric Lac repressor and is expressed by the promoter <i>P<sub>i</sub></i>. The Lac repressor can bind to the <i>lac</i> operators O1, O2 or O3. O3 is located 92&nbsp;bp upstream of O1 and O2 401&nbsp;bp downstream of O1. The inter&#x2011;operator distances are counted from the center of O1 to the center of the distal operator. By binding simultaniously to O1 and O2 or to O1 and O3, the Lac repressor forms a DNA loop that negatively controlls the expression of <i>P<sub>lac</sub></i> (Oehler&nbsp;<i>et&nbsp;al</i>.,&nbsp;1994).  
 
 
 
<!-- Mittleres zentriertes Bild -->
 
<!-- Mittleres zentriertes Bild -->
 
<div class="figure large">
 
<div class="figure large">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/f/fd/T--Bielefeld-CeBiTec--DNA_fold_Lac-repressor.png">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/f/fd/T--Bielefeld-CeBiTec--DNA_fold_Lac-repressor.png">
<p class="figure subtitle"><b>Figure 4: The wild type lac operators.</b><br>The tetrameric Lac repressor can bind either the lac operators O1 and O3 or O1 and O2 to form a DNA loop. The DNA loops efficiently inhibit the transcription by the CAP protein (Oehler et al., 1994).</p>
+
<p class="figure subtitle"><b>Figure&nbsp;4: The wild type <i>lac</i> operators and its regulatory structures.</b><br>The tetrameric Lac repressor can bind either the <i>lac</i> operators O1 and O3 or O1 and O2 to form a DNA loop. The DNA loops efficiently inhibit the transcription by the CAP&nbsp;protein (Oehler&nbsp;<i>et&nbsp;al</i>.,&nbsp;1994).</p>
 
</div>
 
</div>
 
 
In 1994, Oehler et al. showed that an inactivated O2 in its natural position does not decrease the repression by low amounts of tetrameric Lac repressor. Based on this result, we designed our <i>P<sub>lac-tight</sub></i> without the lac operator O2. It was also shown that two weak operators result in a tighter repression than a single strong operator. This can be explained by the thermodynamic concept that a second operator increases the local concentration of the Lac repressor for the neighboring operator. As a consequence there is a higher probability of occupation for two operators by the Lac repressor leading to a tighter repression (Oehler et al., 1996). In 1983, Sadler et al. proposed a ideal lac operator O<sub>id</sub> that binds the Lac repressor 10-fold tighter than the natural strong lac operator O1. O<sub>id</sub> is a inverted repeat of the left half of O1 (Figure 5).  
+
In 1994, Oehler&nbsp;<i>et&nbsp;al</i>. showed that an inactivated O2 in its natural position does not decrease the repression by low amounts of tetrameric Lac repressor. Based on this result, we designed our <i>P<sub>lacO&#x2011;tight1</sub></i> without the <i>lac</i> operator O2. It was also shown that two weak operators result in a tighter repression than a single strong operator. This can be explained by the thermodynamic concept that a second operator increases the local concentration of the Lac repressor for the neighboring operator. As a consequence there is a higher probability of occupation for two operators by the Lac repressor leading to a tighter repression (Oehler&nbsp;<i>et&nbsp;al</i>.,&nbsp;1996). In 1983, Sadler&nbsp;<i>et&nbsp;al</i>. proposed a ideal <i>lac</i> operator O<sub>id</sub> that binds the Lac repressor 10&#x2011;fold tighter than the natural strong <i>lac</i> operator O1. O<sub>id</sub> is an inverted repeat of the left half of O1 (Figure&nbsp;5).  
 
 
 
<!-- Mittleres zentriertes Bild -->
 
<!-- Mittleres zentriertes Bild -->
 
<div class="figure medium">
 
<div class="figure medium">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/c/ce/T--Bielefeld-CeBiTec--Oid_and_O1_compared2.png">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/c/ce/T--Bielefeld-CeBiTec--Oid_and_O1_compared2.png">
<p class="figure subtitle"><b>Figure 5: The <i>lac</i> operator O<sub>id</sub>.</b><br>The inversion is indicated by the arrow for the perfectly symmetric lac operator O<sub>id</sub>. It is the inverted repeat of the left half of O1 (Sadler et al. 1983).</p>
+
<p class="figure subtitle"><b>Figure 5: The <i>lac</i> operator O<sub>id</sub>.</b><br>The inversion is indicated by the arrow for the perfectly symmetric <i>lac</i> operator&nbsp;O<sub>id</sub>. It is the inverted repeat of the left half of O1 (Sadler&nbsp;<i>et&nbsp;al</i>.&nbsp;1983).</p>
 
</div>
 
</div>
 
 
The operator-DNA-operator complex requires energy for the bending process in order to form a DNA loop. Additional energy for a torsion is required when the two <i>lac</i> operators lay on opposite sites of the helical DNA surface. Therefore, the DNA loop formation is energetically favoured for lac operators in phase (Müller et al., 1996). In 1996, Mueller et al. investigated the strength of repression for an inter-operator distance of O<sub>id</sub> and O1 from 57.5 bp up to 1493.5 bp. The repression values were compared to the repression by a single O1 at its natural position. A shorter spacing than 57 .5 bp could not be examined due to the  35 box of the promoter. Phase dependency for the repression was observed for a spacing around 200 bp. That leads to the observation of periodically maxima for repression values (Figure 6).
+
The operator&#x2011;DNA&#x2011;operator complex requires energy for the bending process to form a DNA loop. Additional energy for a torsion is required when the two <i>lac</i> operators lay on opposite sites of the helical DNA surface. Therefore, the DNA loop formation is energetically favoured for <i>lac</i> operators in phase (Müller&nbsp;<i>et&nbsp;al</i>.,&nbsp;1996). In 1996, Mueller&nbsp;<i>et&nbsp;al</i>. investigated the strength of repression for an inter&#x2011;operator distance of O<sub>id</sub> and O1 from 57.5&nbsp;bp up to 1493.5&nbsp;bp. The repression values were compared to the repression by a single O1 at its natural position. A shorter spacing than 57.5&nbsp;bp could not be examined due to the  35&nbsp;box of the promoter. Phase dependency for the repression was observed for a spacing around 200&nbsp;bp. That leads to the observation of periodically maxima for repression values (Figure&nbsp;6).
 
 
 
<!-- Mittleres zentriertes Bild -->
 
<!-- Mittleres zentriertes Bild -->
 
<div class="figure medium">
 
<div class="figure medium">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/c/c5/T--Bielefeld-CeBiTec--Oid_o1_bp_distance.png">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/c/c5/T--Bielefeld-CeBiTec--Oid_o1_bp_distance.png">
<p class="figure subtitle"><b>Figure 6: Repression values dependent on inter-operator distances between O1 and O<sub>id</sub>.</b><br>The repression values refer to the repression of the chromosomal <i>lacZ</i> gene under the control of O1 at its natural positon and O<sub>id</sub> at the indicated position. With 50 tetrameric Lac repressors per cell the repression value is calculated by the specific activity of &#946 &#x2011galactosidase in absence of active Lac repressor devided by the specific activity of &#946 &#x2011galactosidase in the presence of active Lac repressor. The dashed line shows the repression value for a single natural O1 operator (Mueller et al., 1996).</p>
+
<p class="figure subtitle"><b>Figure&nbsp;6: Repression values dependent on inter&#x2011;operator distances between O1 and O<sub>id</sub>.</b><br>The repression values refer to the repression of the chromosomal <i>lacZ</i> gene under the control of O1 at its natural position and O<sub>id</sub> at the indicated position. With 50 tetrameric Lac repressors per cell the repression value is calculated by the specific activity of &#946;&#x2011;galactosidase in absence of active Lac repressor divided by the specific activity of &#946;&#x2011;galactosidase in the presence of active Lac repressor. The dashed line shows the repression value for a single natural O1 operator (Mueller&nbsp;<i>et&nbsp;al</i>.,&nbsp;1996).</p>
 
</div>
 
</div>
 
 
The distance of 70 .5 bp showed the strongest repression value which is 50-fold higher than the natural repression. Repression drops sharply to 15-fold at a 150.5 bp spacing and to threefold at around 600 bp. All inter-operator distances beyond 600 bp kept a twofold increased repression value (Müller et al., 1996).
+
The distance of 70.5&nbsp;bp showed the strongest repression value which is 50&#x2011;fold higher than the natural repression. Repression drops sharply to 15&#x2011;fold at a 150.5&nbsp;bp spacing and to threefold at around 600&nbsp;bp. All inter&#x2011;operator distances beyond 600&nbsp;bp kept a twofold increased repression value (Müller&nbsp;<i>et&nbsp;al</i>.,&nbsp;1996).
 
<br>
 
<br>
According to these results, our <i>P<sub>lac-tight</sub></i> consists of the auxiliary operator O<sub>id</sub> with a 70.5 bp spacing to O1 at its natural position. The residuray sequence like the P lac was kept as the natural lac operon taken from <i>E. coli</i> BL21(DE3) (Figure 7).
+
According to these results, our <i>P<sub>lacO&#x2011;tight1</sub></i> consists of the auxiliary operator O<sub>id</sub> with a 70.5&nbsp;bp spacing to O1 at its natural position. The residuray sequence like the P<sub><i>lac</i></sub> was kept as in the natural <i>lac</i> operon taken from <i>E.&nbsp;coli</i>&nbsp;BL21(DE3) (Figure&nbsp;7).
  
 
<!-- Mittleres zentriertes Bild -->
 
<!-- Mittleres zentriertes Bild -->
 
<div class="figure large">
 
<div class="figure large">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/7/78/T--Bielefeld-CeBiTec--tight_lac_operon_geneious.png">
 
<img class="figure image" src="https://static.igem.org/mediawiki/2017/7/78/T--Bielefeld-CeBiTec--tight_lac_operon_geneious.png">
<p class="figure subtitle"><b>Figure 7: Tight <i>lac</i> operon <i>P<sub>lac-tight</sub></i>.</b><br>The figure with its annotations was created with the software Geneious 10.0.8.</p>
+
<p class="figure subtitle"><b>Figure&nbsp;7: Tight <i>lac</i> operon <i>P<sub>lacO&#x2011;tight1</sub></i>.</b><br>The figure with its annotations was created with the software Geneious 10.0.8.</p>
 
</div>
 
</div>
 
 
Line 93: Line 108:
 
</div>
 
</div>
  
</div>
+
</div>
 +
 
 +
<div class="contentbox">
 +
<div class="bevel tr"></div>
 +
<div class="content">
 +
<h3> Deletion of <i>codA</i> </h3>
 +
 
 +
<article>
 +
To retain the unnatural base pair and keep a sufficient level of unnatural nucleoside triphosphate in the cell, the degradation of the unnatural nucleoside triphosphates must be minimized. In <i>E. coli</i>, the gene <i>codA</i> codes for the cytosine deaminase, an enzyme of the pyrimidine metabolism. The cytosine deaminase catalyzes the reaction of cytosine to uracil. Furthermore, it can catalyze the deamination of isoguanosione and isocytosine. Isoguanosine is formed during oxidative stress by reaction of the radical oxygen species (ROS) •OH with adenine. Other products from reactions of adenine with ROS include 8-oxoadenine and 6-N-hydroxyaminopurine.
 +
 
 +
<div class="figure fourty">
 +
<img class="figure image" src="https://static.igem.org/mediawiki/2017/6/61/T--Bielefeld-CeBiTec--codA_reactions_jpg.jpeg">
 +
<p class="figure subtitle"><b>Figure 8: Reactions catalyzed by the cytosine deaminase.</b><br> A) The conversion of cytosine to uracil is a normal reaction step within the pyrimidine metabolism. B) Isocytosine is converted into uracil by the cytosine deaminase. C) Isoguanine is formed during oxidative stress. Cytosine deaminase catalyzes the reaction to the non-mutagenic xanthine.</p>
 +
</div>
 +
 
 +
In <i>E. coli</i>, the <i>codA</i> gene is part of the <i>codBA</i> operon. <i>CodB</i> codes for the cytosine permease,
 +
while <i>codA</i> encodes cytosine deaminase.The cytosine deaminase hydrolyses cytosine to ammonia and uracil by hydrolytic deamination.
 +
This reaction poses the only way how cytosine can the metabolized in <i>E.coli</i>(Danielsen <i>et al.</i>, 1992). When isoguanosine and isocytosine are provided to the cell,
 +
they are converted to uracil in the wildtype. Therefore, it is necessary to knock out or delete the codA gene to enable the stable
 +
availability of isoguanosine and isocytosine.
 +
 
 +
<div class="figure large">
 +
<img class="figure image" src="https://static.igem.org/mediawiki/2017/9/90/T--Bielefeld-CeBiTec--codBA-operon.jpeg">
 +
<p class="figure subtitle"><b>Figure 9: Arrangement of the <i>codBA</i> operon. </b><br> <i>codA</i> is located downstream of <i>codB</i> and 1.3 kb in size. Both genes overlap by 11 bases. <i>codB</i> codes for cytosine permease, while <i>codA</i> codes for cytosine deaminase.</p>
 +
</div>
 +
 
 +
<i>codA</i> is 1281 bp in size, resulting in a protein composed of 427 amino acids after translation.
 +
The protein CodA is located in the cytosol and has an atomic mass of 47.5 kD.
 +
 
 +
<div class="figure large">
 +
<img class="figure image" src="https://static.igem.org/mediawiki/2017/4/44/T--Bielefeld-CeBiTec--codA-crystal-structure.jpeg">
 +
<p class="figure subtitle"><b>Figure 10: Crystal structure of cytosine deaminase from <i>Escherichia coli</i> complexed with zinc and phosphono-cytosine. </b><br> The structure was determined by X-Ray crystallography with a resolution of 1.71 Å (Hall et al., 2011).</p>
 +
</div>
 +
 
 +
</article>
 +
 
 +
</div>
 +
<div class="bevel bl"></div>
 +
</div>
 +
 
 +
<div class="container">
 +
<div class="contentbox">
 +
<div class="bevel tr"></div>
 +
<div class="content">
 +
            <h2> References </h2>
 +
<article>
 +
<div class="article">
 +
<b>Cobb, R.E., Wang, Y., and Zhao, H.</b> (2015). High-efficiency multiplex genome editing of Streptomyces species using an engineered CRISPR/Cas system. ACS Synth. Biol. <b>4</b>: 723–8.
 +
</div>
 +
<div class="article">
 +
<b>Danielsen, S., Kilstrup, M., Barilla, K., Jochimsen, B., and Neuhard, J.</b> (1992). Characterization of the Escherichia coli codBA operon encoding cytosine permease and cytosine deaminase. Mol. Microbiol. <b>6</b>: 1335–1344.
 +
</div>
 +
<div class="article">
 +
<b>Hall, R.S., Fedorov, A.A., Xu, C., Fedorov, E. V., Almo, S.C., and Raushel, F.M.</b> (2011). Three-Dimensional Structure and Catalytic Mechanism of Cytosine Deaminase. Biochemistry <b>50</b>: 5077–5085
 +
</div>
 +
<div class="article">
 +
<b>Müller, J., Oehler, S., and Müller-Hill, B.</b> (1996). Repression of lac Promoter as a Function of Distance, Phase and Quality of an Auxiliary lac Operator. J. Mol. Biol. <b>257</b>: 21–29.
 +
</div>
 +
<div class="article">
 +
<b>Oehler, S., Amouyal, M., Kolkhof, P., von Wilcken-Bergmann, B., and Müller-Hill, B.</b> (1994). Quality and position of the three lac operators of E. coli define efficiency of repression. EMBO J. <b>13</b>: 3348–3355.
 +
</div>
 +
<div class="article">
 +
<b>Ran, F.A., Hsu, P.D., Lin, C., Gootenberg, J.S., Konermann, S., Trevino, A.E., Scott, D. a, Inoue, A., Matoba, S., Zhang, Y., and Zhang, F.</b> (2013). Double nicking by RNA-guided CRISPR Cas9 for enhanced genome editing specificity. Cell <b>154</b>: 1380–9.
 +
</div>
 +
<div class="article">
 +
<b>Ran, F.A., Hsu, P.D., Wright, J., Agarwala, V., Scott, D.A., and Zhang, F.</b> (2013). Genome engineering using the CRISPR-Cas9 system. Nat. Protoc. <b>8</b>: 2281–2308.
 +
</div>
 +
<div class="article">
 +
<b>Sadler, J.R., Sasmor, H., and Betz, J.L.</b> (1983). A perfectly symmetric lac operator binds the lac repressor very tightly. Proc. Natl. Acad. Sci. U. S. A. <b>80</b>: 6785–9.
 +
</div>
 +
<div class="article">
 +
<b>Simmon, V.F. and Lederberg, S.</b> (1972). Degradation of bacteriophage lambda deoxyribonucleic acid after restriction by Escherichia coli K-12. J. Bacteriol. <b>112</b>: 161–9.
 +
</div>
 +
<div class="article">
 +
<b>Zhang, Y., Lamb, B.M., Feldman, A.W., Zhou, A.X., Lavergne, T., Li, L., and Romesberg, F.E.</b> (2017). A semisynthetic organism engineered for the stable expansion of the genetic alphabet. Proc. Natl. Acad. Sci. <b>114</b>: 1317–1322.
 +
</div>
 +
</article>
 +
 
 +
 
 +
        </div>
 +
<div class="bevel bl"></div>
 +
    </div>
 +
 
</body>
 
</body>
 +
<script>
 +
$("#project").addClass("active");
 +
$("#project-ubp").addClass("active");
 +
$("#project-ubp-preservation").addClass("active");
 +
</script>
 
</html>
 
</html>
 
{{Team:Bielefeld-CeBiTec/Footer}}
 
{{Team:Bielefeld-CeBiTec/Footer}}

Latest revision as of 02:33, 2 November 2017

Retention and Preservation System

Short Summary

The retention and preservation of unnatural nucleoside triphosphates and unnatural bases in the DNA requires a cellular retention and preservation system acting on multiple cellular levels. On one hand, degradation of the unnatural nucleoside triphosphates by enzymes like cytosine deaminase must be prohibited for efficient integration into the DNA. On the other hand, once the unnatural bases are incorporated into the DNA, loss or change of the unnatural bases have to be prevented. We want to achieve this by adapting CRISPR/Cas9 to detect and prevent mutations of the unnatural bases.

Preservation system using Cas9

Due to tautomerisation of isoG and hydrolysis of isoCm and the resulting loss of the unnatural base pair (UBP), there is a need for a system, to preserve the UBP on the plasmid. In 2017, Zhang et al. successfully deployed a CRISPR (clustered regularly interspaced short palindromic repeat)-Cas9 system for retention of a UBP. We adapted this conservation system to our UBP and thus used CRISPR/Cas9 to eliminate all plasmid DNA that had lost the UBP.
The nuclease Cas9 is part of the adaptive immune system of Streptococcus pyogenes, where it induces double strand breaks in the genomic DNA. This enzyme is recruited by a CRISPR RNA (crRNA). A crRNA consists of direct repeats interspaced by variable sequences called protospacer. Those protospacers are derived from foreign DNA and encode the Cas9 guiding sequence (guide RNA). An auxiliary transactivating crRNA (tracrRNA) helps processing the precursor crRNA array into an active crRNA that contains the 20 nucleotide guide RNA. The guide RNA binds to the complementary genomic DNA sequence via Watson‑Crick base pairing. For this binding, the genomic DNA sequence needs to be located upstream of a CRISPR type II specific 5’ NGG protospacer adjacent motif (PAM). Synthetically chimeric single stranded guide RNA (sgRNA) was designed by combining crRNA and tracrRNA. In the sgRNA, only the 20 nucleotide guiding sequence needs to be exchanged for targeting any genomic sequence followed by a PAM sequence (Ran et al., 2013 a, b). The resulting double strand break introduced by Cas9 leads to exonucleolytic degradation of the DNA in prokaryotic cells (Simmon and Lederberg, 1972).
In our case we envision a retention system, where Cas9 cleaves Plasmids at sites where the UBP is absent. This works by using a sgRNA complementary to the DNA sequence without the UBP. In plasmids with the UBP present, the mismatch between isoG/isoCm and sgRNA greatly decreases Cas9 activity (Zhang et al., 2017). In the event of UBP loss, the sgRNA now binds perfectly to the mutated site and restores Cas9 activity which leads to degradation of the mutated plasmid. Consequently, this leads to a retention of the UBP in the plasmids.

Figure 1: UBP conservation system using Cas9.
sgRNAs are targeted against every DNA sequence emerging from UBP loss on a plasmid. A: Loss of the UBP leads to a point mutation. Now a sgRNA can bind to the DNA target sequence. Cas9 is recruited and cleaves the plasmid, which is followed by its degradation. B: Plasmids that contain a UBP in the DNA target sequence lead to a mismatch with every sgRNA. Cas9 does not cleave the plasmid, leading to retention of the UBP.

Because the iGEM competition is based on plasmids that can be submitted to the parts registry we refrained from creating a novel strain for the UBP retention. Consequently, our retention system consists of two plasmids. The first plasmid (BBa_K2201010) contains the nucleotide transporter gene PtNTT2 and the cas9 on a pSB1K3 backbone. This plasmid was transformed into E. coli BL21(DE3). After producing chemically competent cells, another transformation was done with the second plasmid. The second plasmid carries the UBP and the sgRNAs for recruiting Cas9 on a pSB3C5 plasmid. Since two large plasmids produce a great metabolic stress for an organism, a more sufficient UBP retention could be achieved with a newly created E. coli strain, that has the genes PtNTT2 and cas9 in its genomic DNA. We provide the parts that are required to create a repair template for the knock in strategy using a CRISPR/Cas9 system. These parts were designed to perform a genomic integration into E. coli BL21(DE3). The genomic knock in can be done according to the protocol by Cobb et al., 2014 using the pCRISPomyces plasmid system. A sgRNA was designed searching for a reverse sequence with the constraint N(16)R(4)NGG (Cobb et al., 2004) inside the coding sequence of the arsB gene (Zhang et al., 2017) of the E. coli genome. arsB is coding for an arsenic efflux pump membrane protein. As a result, 5‘ TATTGTTCATAATAGAAGAG 3‘ turned out to be the suggested guiding sequence with the highest on target activity score with being unique within the complete genome. The required repair template is designed with 1 kb long flanking sequences and two terminators BBa_B0015 flanking the gene of interest for example cas9 (Figure 2). To be able to strictly regulate the expression of cas9 we rationally designed an optimized IPTG inducible promoter PlacO‑tight1.

Figure 2: The repair template (BBa_K2201028) for genomic integration of cas9 into E. coli BL21(DE3).
The 1 kb left flanking sequence (LFS: BBa_K2201021) and right flanking sequence (RFS: K2201022) are taken from the genome of E. coli BL21(DE3). Inside the genome LFS and RFS are directly flanking the coding sequence of arsB. The strong terminators from BBa_B0015 were used to avoid basal expression and to stop transcription right after cas9. There are composite parts the consist of LFS + terminator + PlacO‑tight1 (BBa_K2201024) and terminator + RFS (BBa_K2201025) that can be assembled with any coding sequence of interest to create a repair template. The coding sequence of cas9 is taken from the pCRISPomyces plasmid system by Cobb et al., 2014 and is originated from S. pyogenes. It is negatively regulated by the IPTG‑inducible promoter PlacO‑tight1.

This composite BioBrick from Figure 2 needs to be restricted with NotI to separate the linear repair template from the pSB1C3 backbone. Together with the target plasmid containing the sgRNA and the pCRISPomyces plasmid containing the cas9 all three elements were co‑transfected into E. coli BL21(DE3). The genomic integration was verified by sequencing.

Optimization of the negatively regulating promoter Plac

We designed a tightly repressed lac operon called PlacO‑tight1 to achieve a low basal transcription rate. The lac operon was extensively researched starting from the midst of last century. This showed the role of the lac repressor and the inducible nature of the lac promoter, which results in 1000-fold increased expression in its activated state (Müller et al., 1996). The wild type lac operon consists of the genes lacZ, lacY and lacA that are transcribed from the lac promoter Plac into a polycistronic mRNA (Figure 4). These genes code for the proteins β‑galactosidase, Lac permease and Lac transacetylase (Oehler et al., 1994).

Figure 3: DNA map of the wild type lac operon and its transcriptional and translational products (Oehler et al., 1994).

The transcription is constitutively activated by the CAP protein. lacI codes for the tetrameric Lac repressor and is expressed by the promoter Pi. The Lac repressor can bind to the lac operators O1, O2 or O3. O3 is located 92 bp upstream of O1 and O2 401 bp downstream of O1. The inter‑operator distances are counted from the center of O1 to the center of the distal operator. By binding simultaniously to O1 and O2 or to O1 and O3, the Lac repressor forms a DNA loop that negatively controlls the expression of Plac (Oehler et al., 1994).

Figure 4: The wild type lac operators and its regulatory structures.
The tetrameric Lac repressor can bind either the lac operators O1 and O3 or O1 and O2 to form a DNA loop. The DNA loops efficiently inhibit the transcription by the CAP protein (Oehler et al., 1994).

In 1994, Oehler et al. showed that an inactivated O2 in its natural position does not decrease the repression by low amounts of tetrameric Lac repressor. Based on this result, we designed our PlacO‑tight1 without the lac operator O2. It was also shown that two weak operators result in a tighter repression than a single strong operator. This can be explained by the thermodynamic concept that a second operator increases the local concentration of the Lac repressor for the neighboring operator. As a consequence there is a higher probability of occupation for two operators by the Lac repressor leading to a tighter repression (Oehler et al., 1996). In 1983, Sadler et al. proposed a ideal lac operator Oid that binds the Lac repressor 10‑fold tighter than the natural strong lac operator O1. Oid is an inverted repeat of the left half of O1 (Figure 5).

Figure 5: The lac operator Oid.
The inversion is indicated by the arrow for the perfectly symmetric lac operator Oid. It is the inverted repeat of the left half of O1 (Sadler et al. 1983).

The operator‑DNA‑operator complex requires energy for the bending process to form a DNA loop. Additional energy for a torsion is required when the two lac operators lay on opposite sites of the helical DNA surface. Therefore, the DNA loop formation is energetically favoured for lac operators in phase (Müller et al., 1996). In 1996, Mueller et al. investigated the strength of repression for an inter‑operator distance of Oid and O1 from 57.5 bp up to 1493.5 bp. The repression values were compared to the repression by a single O1 at its natural position. A shorter spacing than 57.5 bp could not be examined due to the 35 box of the promoter. Phase dependency for the repression was observed for a spacing around 200 bp. That leads to the observation of periodically maxima for repression values (Figure 6).

Figure 6: Repression values dependent on inter‑operator distances between O1 and Oid.
The repression values refer to the repression of the chromosomal lacZ gene under the control of O1 at its natural position and Oid at the indicated position. With 50 tetrameric Lac repressors per cell the repression value is calculated by the specific activity of β‑galactosidase in absence of active Lac repressor divided by the specific activity of β‑galactosidase in the presence of active Lac repressor. The dashed line shows the repression value for a single natural O1 operator (Mueller et al., 1996).

The distance of 70.5 bp showed the strongest repression value which is 50‑fold higher than the natural repression. Repression drops sharply to 15‑fold at a 150.5 bp spacing and to threefold at around 600 bp. All inter‑operator distances beyond 600 bp kept a twofold increased repression value (Müller et al., 1996).
According to these results, our PlacO‑tight1 consists of the auxiliary operator Oid with a 70.5 bp spacing to O1 at its natural position. The residuray sequence like the Plac was kept as in the natural lac operon taken from E. coli BL21(DE3) (Figure 7).

Figure 7: Tight lac operon PlacO‑tight1.
The figure with its annotations was created with the software Geneious 10.0.8.

Deletion of codA

To retain the unnatural base pair and keep a sufficient level of unnatural nucleoside triphosphate in the cell, the degradation of the unnatural nucleoside triphosphates must be minimized. In E. coli, the gene codA codes for the cytosine deaminase, an enzyme of the pyrimidine metabolism. The cytosine deaminase catalyzes the reaction of cytosine to uracil. Furthermore, it can catalyze the deamination of isoguanosione and isocytosine. Isoguanosine is formed during oxidative stress by reaction of the radical oxygen species (ROS) •OH with adenine. Other products from reactions of adenine with ROS include 8-oxoadenine and 6-N-hydroxyaminopurine.

Figure 8: Reactions catalyzed by the cytosine deaminase.
A) The conversion of cytosine to uracil is a normal reaction step within the pyrimidine metabolism. B) Isocytosine is converted into uracil by the cytosine deaminase. C) Isoguanine is formed during oxidative stress. Cytosine deaminase catalyzes the reaction to the non-mutagenic xanthine.

In E. coli, the codA gene is part of the codBA operon. CodB codes for the cytosine permease, while codA encodes cytosine deaminase.The cytosine deaminase hydrolyses cytosine to ammonia and uracil by hydrolytic deamination. This reaction poses the only way how cytosine can the metabolized in E.coli(Danielsen et al., 1992). When isoguanosine and isocytosine are provided to the cell, they are converted to uracil in the wildtype. Therefore, it is necessary to knock out or delete the codA gene to enable the stable availability of isoguanosine and isocytosine.

Figure 9: Arrangement of the codBA operon.
codA is located downstream of codB and 1.3 kb in size. Both genes overlap by 11 bases. codB codes for cytosine permease, while codA codes for cytosine deaminase.

codA is 1281 bp in size, resulting in a protein composed of 427 amino acids after translation. The protein CodA is located in the cytosol and has an atomic mass of 47.5 kD.

Figure 10: Crystal structure of cytosine deaminase from Escherichia coli complexed with zinc and phosphono-cytosine.
The structure was determined by X-Ray crystallography with a resolution of 1.71 Å (Hall et al., 2011).

References

Cobb, R.E., Wang, Y., and Zhao, H. (2015). High-efficiency multiplex genome editing of Streptomyces species using an engineered CRISPR/Cas system. ACS Synth. Biol. 4: 723–8.
Danielsen, S., Kilstrup, M., Barilla, K., Jochimsen, B., and Neuhard, J. (1992). Characterization of the Escherichia coli codBA operon encoding cytosine permease and cytosine deaminase. Mol. Microbiol. 6: 1335–1344.
Hall, R.S., Fedorov, A.A., Xu, C., Fedorov, E. V., Almo, S.C., and Raushel, F.M. (2011). Three-Dimensional Structure and Catalytic Mechanism of Cytosine Deaminase. Biochemistry 50: 5077–5085
Müller, J., Oehler, S., and Müller-Hill, B. (1996). Repression of lac Promoter as a Function of Distance, Phase and Quality of an Auxiliary lac Operator. J. Mol. Biol. 257: 21–29.
Oehler, S., Amouyal, M., Kolkhof, P., von Wilcken-Bergmann, B., and Müller-Hill, B. (1994). Quality and position of the three lac operators of E. coli define efficiency of repression. EMBO J. 13: 3348–3355.
Ran, F.A., Hsu, P.D., Lin, C., Gootenberg, J.S., Konermann, S., Trevino, A.E., Scott, D. a, Inoue, A., Matoba, S., Zhang, Y., and Zhang, F. (2013). Double nicking by RNA-guided CRISPR Cas9 for enhanced genome editing specificity. Cell 154: 1380–9.
Ran, F.A., Hsu, P.D., Wright, J., Agarwala, V., Scott, D.A., and Zhang, F. (2013). Genome engineering using the CRISPR-Cas9 system. Nat. Protoc. 8: 2281–2308.
Sadler, J.R., Sasmor, H., and Betz, J.L. (1983). A perfectly symmetric lac operator binds the lac repressor very tightly. Proc. Natl. Acad. Sci. U. S. A. 80: 6785–9.
Simmon, V.F. and Lederberg, S. (1972). Degradation of bacteriophage lambda deoxyribonucleic acid after restriction by Escherichia coli K-12. J. Bacteriol. 112: 161–9.
Zhang, Y., Lamb, B.M., Feldman, A.W., Zhou, A.X., Lavergne, T., Li, L., and Romesberg, F.E. (2017). A semisynthetic organism engineered for the stable expansion of the genetic alphabet. Proc. Natl. Acad. Sci. 114: 1317–1322.