You might know about the CRISPR-based technology called SHERLOCK (Specific High-sensitivity Enzymatic Reporter un-LOCKing) which aims to detect nucleic acid with a high specificity. We, the iGEM Montpellier team, developed Shell’lock which is based on this well known technique and that we adapted to the detection of oyster pathogens. In this page, you will find all the design details of Shell’lock.

Modules


Module 1

Guide and Target Design

Module 2

Paper Based Test

Module 3

Fuorescent Output

Module 1: Guide And Target Design


Overview

The SHERLOCK detection system is based on the specific recognition of a target RNA that is present in the sample by a specifically designed guide RNA (crRNA), which acts as a guide for the Cas13c protein. The specificity is easily programmable and relies on introducing in the crRNA a 28 nucleotide long sequence that is perfectly complementary to the target RNA . In our strategy, detecting a specific organism then consists in identifying a sequence that is unique to this organism. To obtain a proof of concept in the lab, we first identified suitable target sequences from the genome of pathogens and then worked with the corresponding in vitro translated RNAs.

Positive Control

As positive control, we used the synthetic sequence and matching guide RNA designed by Kellner et al. in “SHERLOCK: Nucleic acid detection with CRISPR nuclease” (1).

Genome mining and design of the target sequences

We designed multiple sequences from different genes of our organisms of interest. We focused on different type of genes, the detection of which can provide relevant information for each of them:

  • Toxicity genes that are directly linked to pathogenicity (vam (2) and toxR (3))
  • Gene that are specific to v.aestuarianus (dnaJ (4))
  • Gene for DNA gyrase (gyrB (5))
  • Gene involved in cell division, drug resistance and sensitivity (mreB (6))
  • Gene coding for uridylate kinase, involved in pyrimidine ribonucleotide biosynthesis (pyrH)
  • Gene involved in DNA replication, recombination and repair (recA and topA)
  • Ribosomal RNA gene (16s RNA(7))
  • Developmental genes that have been involved in pathogenicity in other vibrio species, but are not specific to V.aestuarianus (flaA (8) and vspR (9))

All the genes are listed in the table below.

Vibrio aestuarianus gene Strain Accession number
AY605667.1 01/32 AY605667.1
toxR 02/041 AM183570.1
flaA 02/041 Sequence extracted from mage genscope (10)
vspR 02/041 Sequence extracted from mage genscope (11)
dnaJ GTC2706 AB263018.1
16s rRNA gene 02/041 AJ845017.1
gyrB 02/041 AJ852513.1
mreB 02/041 Sequence extracted from mage genscope (12)
pyrH 02/041 Sequence extracted from mage genscope (13)
recA 02/041 Sequence extracted from mage genscope (14)
topA 02/041 Sequence extracted from mage genscope(15)
Table 1: Vibrio aestuarianus genes used as target.

Although Cas13 recognizes RNA, to lower the cost of our experiments, we decided to design DNA template sequences and then perform transcription to obtain target and guide RNAs. The template for all the target sequences is described in figure 1.

Figure 1: DNA target sequence template. The 28 nucleotides that are complementary to the guide sequence are flanked by 40 nucleotides on each side that constitute the natural context of the target in the reference genome sequence. This 108 nucleotides-long sequence is preceded by the T7 promoter and followed by a T7 terminator. This entire sequence is flanked by two spacers that provide standard priming sites for PCR amplification.

Procedure to obtain RNA target sequences:


To design the synthetic target sequences, we had to choose 28 nucleotides in the gene sequence. We randomly chose these 28 nucleotides but as the Cas13a protein has dinucleotide preferences (cuts preferentially polyU and A/U (1)) we tried to find regions enriched in these nucleotides. (1))) we tried to find regions enriched in these nucleotides.

Design of the guide RNA

The guide RNA is complementary to the 28 nucleotides of the coding strand of the target sequence. We also decided to design the guides as DNA to have enough material for a large number of experiments. The template used for all the guides is described in figure 2.

Figure 2: Guide DNA template. We designed primers that hybridize themselves (20 nucleotides of base pairing) during the PCR. The forward primer (in orange) is the same for all the PCR but the reverse primer (in pink) changes depending on the specific 28 nucleotides of the guide sequence. The PCR product consists of a T7 promoter that drives the expression of a scaffold as described by Kellner et al.(1), to which a 28 nucleotides sequence is perfectly complementary to the target sequence.

To filter the guide RNA and to choose the most specific one, we developed a program. The goal is to have an automatic pipeline that Blast a sequence against NCBI databases. We used it to blast our potential guide RNA candidates to check if there are alignments with other sequences. The analysis is run via python through the NCBI server. The principle is shown in figure 3:

Figure 3: Blast program principle. First, open the sequence and transform the file type as an object, then call the Blast function with adjusted parameters, finally, extract the results and filter them (we applied filters to only have results of organisms that can be found in water).

We identified some limits of this program:

  • The Expect Value Threshold default setting is reduced to 0,05
  • The maximum number of target sequences limit is no more than 5000
  • The maximum allowed query length for nucleotide queries (blastn, blastx and tblastx) is 1,000,000 and 100,000 for protein queries (blastp and tblastn) (16)
  • Submitting searches on off-hours (8 pm to 8 am EST) may provide better throughput.

Even if we identified these limits, we think that our program is a necessary tool to identify best guide candidates. In fact, it avoids the design of guide RNA that will recognize the wrong target sequence and thus produce false positives.

The script is available in the iGEM software GitLab.

We analyzed our results and decided to remove sequences that have a high probability to match with sequences in other organisms than Vibrio Aestuarianus. We end up with 9 guide RNAs that can target the sequences of interest: two guides targeting 16s, two guides targeting dnaJ, oneguide targeting flaA, one targeting vam, one targeting toxR and the last one targeting vspR.
All the DNA sequences used are reported in this table.

References

1. J.Kellner M, Koob J, S.Gootenberg J, O.Abudayyeh O, Zhang F. SHERLOCK: Nucleic acid detection with CRISPR nucleases. 2020 Mar; Available from: 10.1038/s41596-019-0210-2
2. Lebreuche Y, Le Roux F, Henry J, Zatylny C. Vibrio aestuarianus zinc metalloprotease causes lethality in the Pacific oyster Crassostrea gigas and impairs the host cellular immune defenses. 2010 Nov;753–8. 10.1016/j.fsi.2010.07.007
3. Garnier M, Lebreuche Y, Nicolas JL. Molecular and phenotypic characterization of Vibrio aestuarianus subsp. francensis subsp. nov., a pathogen of the oyster Crassostrea gigas. 2008 Oct; Available from: http://dx.doi.org/10.1016/j.syapm.2008.06.003
4. Saulnier D, De Decker S, Haffner P. Real-time PCR assay for rapid detection and quantification of Vibrio aestuarianus in oyster and seawater: A useful tool for epidemiologic studies. 2009 May; Available from: https://doi.org/10.1016/j.mimet.2009.01.021
5. DNA gyrase, subunit B [Internet]. Available from: https://www.ebi.ac.uk/interpro/entry/InterPro/IPR011557/
6. mreB gene [Internet]. Available from: https://mage.genoscope.cns.fr/microscope/mage/seq.php?id=12834278
7.Hoffman M, W Brown E, CH Feng P, Fischer M, R Monday S. PCR-based method for targeting 16S-23S rRNA intergenic spacer regions among Vibrio species. 2010 Mar; Available from: https://doi.org/10.1186/1471-2180-10-90
8. L.Milton D, O’Toole R, Hörstedt P. Flagellin A Is Essential for the Virulence of Vibrio anguillarum. 1996 Mar; Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC177804/pdf/1781310.pdf
9. H.Rashid M, Rajanna C, Zhang D, Pascuale V. Role of exopolysaccharide, the rugose phenotype and VpsR in the pathogenesis of epidemic Vibrio cholerae. 2004 Jan; Available from: https://doi.org/10.1016/S0378-1097(03)00879-6
10. flaA gene [Internet]. Available from: https://mage.genoscope.cns.fr/microscope/mage/getInfoLabel.php?id=12832451
11. vspR gene [Internet]. Available from: https://mage.genoscope.cns.fr/microscope/mage/getInfoLabel.php?id=20021630
12. mreB gene [Internet]. Available from: https://mage.genoscope.cns.fr/microscope/mage/seq.php?id=12834278
13. pyrH gene [Internet]. Available from: https://mage.genoscope.cns.fr/microscope/mage/seq.php?id=12832370
14.recA gene [Internet]. Available from: https://mage.genoscope.cns.fr/microscope/mage/seq.php?id=12834138
15. topA gene [Internet]. Available from: https://mage.genoscope.cns.fr/microscope/mage/seq.php?id=12833395
16.New BLAST default parameters and search limits coming in September [Internet]. Available from: https://ncbiinsights.ncbi.nlm.nih.gov/2020/06/18/new-blast-settings/

Module 2: Paper Based Test


In this section we present our lateral flow assay using Shell’ Lock to detect oyster pathogens. Current methods of pathogen detection require the use of conventional molecular biology techniques such as PCR (1), (2). To stay close to the first concern of our project, the oyster farmers, we wanted to develop an easy detection technique. In this section we describe the design principle behind using paper-strips and lateral flow assay to perform the Shell’Lock reaction.

Experimental Design:

The test can be decomposed in two main steps:

  • A. Main reaction
  • B. Paper based revealing

A. Main reaction

The first step of our assay is the Cas13 activation and probe cleavage. This can be preceded by an optional target amplification step. This step will be performed in a first test tube. A second test tube contains the enzyme (Cas13a), the guide RNA (crRNA) and the probe. The probe is a15 uridine RNA molecule conjugated to biotin in 5’ and FAM (6-fluorescein amidite) in 3’. In the reaction tube one can insert the product of the RPA or directly the tested sample (see results for detailed protocol).

B. Paper based revealing

To the reaction tube are added gold nanoparticles coupled with an anti-FITC antibody (that can bind to FAM) and reaction buffer (Milenia hybrid-detect kit, CITY AND COUNTRY). The strip contains two different types of capture molecules: a biotin ligand (control band) and an anti-rabbit antibody (test band). Both are fixed on the strip and serve to capture antibodies. The sample is deposited on a sample pad, bycapillarity the aqueous solution will migrate. The anti-FITC coupled gold nanoparticles bind to the RNA- coupled FITC molecules and co-migrate with the sample on the paper strip. In the case of a SHERLOCK reaction two different scenarios can occur:

  • In the case of a negative test, the probes will be under their native form (FAM-RNA-Biotin) and be captured by the control band (biotin ligand). The gold nanoparticles will recognize the FITC which is located at the control band and a visible readout can be done
  • In the case of a positive test, the probes will be cutted by the Cas13a. In that case, the biotin part will still be recognized by the biotin ligand molecule but the FITC will continue to migrate. The gold nanoparticle will still recognize the FITC part of the probe but will also be captured by the anti-rabbit antibody present at the test band. This will result in a visible band at the test position.

Design of all the components:

Figure 1: Functioning of the lateral flow detection assay

Semi-quantitative readout:

In order to have a semi-quantitative readout, we used the absorbance properties of the gold nanoparticles. As shown in figure 2, gold nanoparticles have a peak of absorbance around 500 nm. Using a gel imager we exposed our readouts to white light and obtained an RGB image of our test. Using the ImageJ software (3), we split the channels and inverted the green channel in order to have access to the absorbance measurement. We then measured the intensity of each band and the intensity of the background. Using a custom made python pipeline we subtracted the background intensity from the intensity of each band and plotted the ratio of intensity between the test and control band. An example of such measurement can be found in figure 3 . More details can be found in the results section.

Figure 2: Absorbance spectra of a Gold nanoparticle. Credits: Taku Ueki
Figure 3: Quantification of the SHERLOCK lateral flow assay. The SHERLOCK reaction was carried out as described before (link). NC: negative control, SHERLOCK reaction without LwCas13a.

Limits and troubleshooting of the test:

  • If the gold nanoparticles are in excess compared to the FITC epitope, they can bind to the test band, even if not associated with the cut probe. This causes a faint but visible band even in the negative conditions. We optimized the concentrations of our test accordingly to decrease as much as possible this “false positive”.
  • Quantifications depend on the time between the end of the experiment and imaging. To avoid these causes of variations we set up a standardized protocol to image.

Results

Results of lateral flow experiments can be found on this page.

References

F. S. Le Guyader et al., « Detection of Multiple Noroviruses Associated with an International Gastroenteritis Outbreak Linked to Oyster Consumption », J Clin Microbiol, vol. 44, no 11, p. 3878‑3882, nov. 2006, doi: 10.1128/JCM.01327-06
2. P. K. C. Cheng, D. K. K. Wong, T. W. H. Chung, et W. W. L. Lim, « Norovirus contamination found in oysters worldwide », J. Med. Virol., vol. 76, no 4, p. 593‑597, août 2005, doi:10.1002/jmv.20402
3. V. Séchet et al., « Characterization of toxin-producing strains of Dinophysis spp. (Dinophyceae) isolated from French coastal waters, with a particular focus on the D. acuminata-complex », Harmful Algae, vol. 107, p. 101974, juill. 2021, doi: 10.1016/j.hal.2021.101974. 10.1002/jmv.20402
4.C. Myhrvold et al., « Field-deployable viral diagnostics using CRISPR-Cas13 », Science, vol. 360, no 6387, p. 444‑448, avr. 2018, doi:10.1126/science.aas8836
5.C. A. Schneider, W. S. Rasband, et K. W. Eliceiri, « NIH Image to ImageJ: 25 years of image analysis », Nat Methods, vol. 9, no 7, Art. no 7, juill. 2012, doi:10.1038/nmeth.2089

Module 3: Fluorescent Output


Our diagnostic tool is based on the SHERLOCK (1) system, which itself relies on the Cas13a enzymatic reaction. Cas13a is an enzyme that possesses both a specific and unspecific nuclease activity (2) . In other words, Cas13a cleaves target RNA sequences that are complementary to a guide RNA and other non-specific and off-target RNA sequences. We are using that property of Cas13a coupled with a variety of RNA probes in order to detect the presence of certain pathogens in Thau’s lagoon water.

Our project was split in two different tests as explained in other sections:

  • An in vitro test used for the characterisation of the reaction.
  • A field test directly usable by oyster farmers. This test relies on the use of antibodies and RNA probes linked to molecules recognized by antibodies in a lateral flow assay.

Our diagnostic is intended to be easy to use and rapid, making it qualitative. In order to access more quantitative outputs and perform more in-depth analysis we developed an in vitro test.
The test relies on the use of a specific RNA probe linked to a fluorophore and a quencher. Thus using fluorescence output we developed a plate-reader-based Shell’lock detection. In this section we are going to detail every step of the experimental design of this approach.

Designing of the probe

The probe is the core of the fluorescent detection. We designed the cleavage reporter according to Gootemberg et al. (3) . This reporter is composed of the fluorophore 6-FAM (IDT), a poly U sequence and the quencher BHQ_1 (IDT).

Mechanism of quenching

Quenchers are a class of molecules that can, through different mechanisms suppress the fluorescence of a fluorophore. The most common form of quenching is through Förster resonance energy transfer (FRET). This mechanism is based on the overlapping of the excitation spectra of the quencher with the emission spectra of the fluorophore. The energy of the fluorescent emission by the fluorophore is transmitted to the quencher and can then be emitted at another wavelength or dissipated. In the case of Dark quencher the energy is dissipated which creates a decrease in the fluorescence detection.

Many mechanisms of quenching have been described in the litterature, involving energy transfer, dimerization of the fluorophore, chemical modification of the fluorophore… We decided to use Black Hole Quencher®-1 quencher from IDT. The mechanism of quenching of this particular fluorophore is the property of IDT.

Fluorophore

The choice of the fluorophore was adapted from the quencher. We chose a fluorophore that corresponded to the characteristics imposed by the quencher i.e with an overlap between the emission spectra of the fluorophore and the absorbance of the quencher. We decided to choose FAM (6-Carboxyfluorescein ).

Use of the Cas

For our experiments we followed the protocol published in the SHERLOCK publication by Kellner et al. (1) . For this we used the LwCas13a (Leptotrichia wadei). Multiple Cas13a have been tested and used for the SHERLOCK reaction (4) . The plasmid of expression of the Cas used can be found on Addgene.

Principle of the reaction

The reaction starts by the recognition of the guide RNA by the Cas protein. This step depends on the presence of a scaffold sequence in the guide in 5’ that is recognized by LwCas13a. This sequence was obtained by aligning varius guide RNA found in the literature using the same Cas (LwCas13a). After assembly of the Cas-guide complex, there is specific recognition of the target. This specificity comes from the “spacer” sequence in the guide (20nt) that is complementary to the target RNA. Upon recognition the Cas protein will become activated and will perform it’s nuclease activities:

  • The unspecific “collateral” nuclease activity, which will cut unspecifically RNA present in physical proximity with the complex. In our test this will lead to the cutting of the RNA probes and the beginning of fluorescence measurement
  • The specific nuclease activity, which will cleave the target RNA. Upon cleave of the target the Cas-guide will unbind from the target and be recycled for further use. This step will also stop the collateral activity.

A scheme of the SHERLOCK reaction can be seen in figure 1.

Figure 1: Schematic representation of the Shell’lock reaction

Experimental protocol

The components of the SHERLOCK mix can be found in our protocol. The concentration of the components were adjusted for our experimental setup after characterization experiments. We performed kinetic experiments to follow the evolution of the concentration of the cleaved probes. The experiments ran over 8h by measuring at 5 minutes intervals. The complete protocols for the different experiments we carried out can be found in our protocol.

Wet-lab

The Protocol is found here: Protocol