41J Blog

Reading through some of Stratos’ more recent patents I came across a non-expandomer sequencing approach which I found quite interesting. The patent shows that Stratos had been thinking about other approaches to nanopore based sequencing. This suggests that the Roche acquisition may not just have been for their expandomer IP…

The Approach

My understanding of the approach from [1] is summarized in the schematic above. Essentially you construct a polymerase [2] that has a “tether” attached to it. The tether is composed of a PEG repeat, which threads through a nanopore [3]. The PEG region has a short oligo on the end. Once it’s threaded through the pore, another oligo can be hybridized to it. This secures the tether in the pore.

With polymerase-tether complex in the pore, the system should look something like the diagram above. The polymerase is secured on the top of the pore. Due to its charge, I guess the polymerase would be pulled toward the pore, but it’s too big to pass through.

As in a standard Ionic/Protein nanopore setup, there’s a bias voltage, and ionic current passing through the pore. The polymerase is now blocking the pore. This causes a significant reduction in current flow. One of the figures in the patent illustrates this (and looks like experimental data):

Now, a template to be sequenced is introduced. Single stranded DNA enters the polymerase, and synthesis of the complementary strand occurs. However, the template DNA doesn’t interact with the pore directly.

The idea is that as the polymerase incorporates bases it will undergo confirmational changes. The patent suggests that these confirmational changes can be up to about 1nm.

The idea here is that different confirmational changes will cause the polymerase to block the pore with varying efficiency. So, for example when incorporating a base it might go through a series of confirmational changes, which result in less current flow, then more current flow, then back to baseline.

Ideally, each base would induce a distinct conformational change, and current blockage. The current trace can then be used to infer the sequence of the DNA template.

However, if you’re only able to detect incorporation/non-incorporation then nucleotides could be sequentially introduced.

The approach is somewhat reminiscent of that proposed by Roswell. Both of which seem to be trying to detect confirmational changes in the polymerase. In comparison to the Roswell approach, this seems like a simpler setup. The polymerase is also closer to the “sensor” which might make detecting confirmation changes easier.

Contrasting this with other ionic nanopore approaches, it’s kind of nice that the strand doesn’t go through the pore. Ideally in the Stratos approach there would be at most 4 different signal types, one for each base. Rather than the signal being some combination of all the nucleotides currently in the pore. Theoretically this could result in lower error rates.

Overall the idea seems at least plausible, and I’ll be curious to see how it plays out.

References and Notes

[1] http://www.freepatentsonline.com/20170159115.pdf

[2] A number of enzymes that process DNA could work (for example an exonuclease). And the patent states this, but a polymerase seems the most likely option.

[3] They talk about and show various protein nanopores. But I imagine a solid state nanopore may also be viable, particularly as the dimensions of the pore maybe slightly less critical in this approach.

[4] “The tethers were constructed of three domains (i.e., “segments”): 1) a polyethylene glycol (PEG) repeat region, located proximal to the polymerase and designed to span the nanopore channel; 2) a short oligonucleotide, designed to hybridize to a single-stranded oligonucleotide on the opposite side of the nanopore relative to the polymerase to anchor the assembly; and 3) a negatively charged phosphoramidite tail, located most distal to the polymerase and designed to facilitate threading of the tether through the nanopore. FIG. 9 is a SDS/PAGE gel that shows the size of the unmodified KF polymerase (lane 1), the KF-tether 1 conjugate (lane 2) and the KF-tether 2 conjugate (lane 3). As expected, the conjugates show an increase in mass compared to the unmodified polymerase.”

[5] “When this polymerase complexes with a nucleotide that is the complement to the template base in the next extension position the polymerase reconfigures into what is referred to in the art as a “closed” conformation. At a more detailed structural level, the transition from the open to closed conformation is characterized by relative movement within the polymerase resulting in the “thumb” domain and “fingers” domain being closer to each other. In the open conformation the thumb domain is further from the fingers domain, akin to the opening and closing of the palm of a hand. In various polymerases, the distance between the tip of the finger and the thumb can change up to 10 angstroms between the “open” and “closed” conformations. The distance between the tip of the finger and the rest of the protein domains can also change up to 10 Angstroms. It will be understood that this change will be exploited in a method set forth herein.”

[6] This figure in the patent is supposed to show the threading/polymerase. But I find it a bit confusing.

Other quotes….

“The present disclosure relates to methods and constructs for single molecule electronic sequencing of template nucleic acids. The constructs are molecular sensor complexes which comprise a processive nucleic acid processing enzyme localized to a nanopore. Conformational changes in the enzyme induced by single nucleic acid processing events are transduced into electric signals by the nanopore, which are used to identify individual nucleotides. The methods can include the steps of providing a membrane with the nanopore and the enzyme complexed with a template nucleic acid localized proximal to an opening in the pore, contacting the enzyme with an ion conductive reaction mixture including the reagents required for nucleic acid processing, providing a voltage drop across the pore that induces ion current through the pore that is modulated by conformational changes in the enzyme, measuring current through the pore over time to detect nucleotide-dependent conformational changes in the enzyme, and identifying the type of nucleotide processed by the enzyme using current modulation characteristics, thus determining sequencing information about the nucleic acid molecule.”

“The polymerase is secured to the pore by hybridizing a short oligonucleotide anchor to the tether construct on the distal side of the nanopore.”

“These results indicate that a tether and a tether-polymerase conjugate can be anchored to a nanopore and, moreover, that the resulting complex can generate reproducible electrical signals. Polymerase-nanopore complexes are thus capable of modulating current flow through the pore and show promise as useful sensors to transduct mechanical events into electrical signals.”

“One or more of the transitions that a polymerase undergoes when adding a nucleotide to a nucleic acid can be detected using a molecular sensor complex as described herein.”

“FIG. 4C depicts the polymerase in a second, ™, “closed” configuration, which is induced, e.g., by binding of incoming nucleotide 605 to form a correct base pair with the template nucleic acid. In this second configuration, the degree to which the enzyme physically occludes the pore is reduced, and consequently the flow of current through the pore will increase. Such modulation of current flow generates an electronic signal specific for nucleotide species”

“In one embodiment, each of the four nucleotides induces a different polymerase conformation, as illustrated in FIG. 4C. . The movement of the polymerase during the incorporation of a nucleotide will modulate the ion current through the pore in a characteristic and reproducible manner, generating a signature electric signal.”

“In another embodiment, the average amplitude of the current modulation doesn’t change, but rather the noise in the current modulation changes as a single nucleotide is bound and incorporated. In yet another embodiment, the current modulation system only indicates an incorporation event but does not discriminate the base type. In this embodiment, the sequence information about a nucleic acid is obtained by sequentially flooding the senor complex with one of four reaction mixtures containing one of the four nucleotides and detecting the presence or absence of an electric signal.”

“For proper function of the molecular sensor complexes of the present invention, it is necessary that the enzyme be stably localized to the pore in sufficiently close proximity to reliably influence, or modulate, current flow through the pore. Several alternative localization and/or attachment structures or compositions are contemplated by the present invention, some which are illustrated schematically in FIGS. 5-8. FIG. 5A depicts one embodiment in which enzyme 500 is localized to pore 220 by covalent attachment to tethering structure 325, herein referred to simply as a “tether”. Tethers may be designed to thread through the lumen of the pore, from one side of membrane 100 to the other. Tethers may comprise one or more structural domains, or “segments”, designed to perform one or more functions.”

I’ve previously written about Apton Biosystems. When I wrote that post there wasn’t much to go on. However, a patent [2] has recently been published which reveals a bit more information.

The motivation stated in the patent is that “to reach a $10 30× genome”…”the amount of data per unit area needs to increase by 100 fold”. Elsewhere in the patent they mention that the prior art is a pitch of 1 micron. HiSeq wells were ~500nm. So they want decrease well size to ~100nm.

This premise, is slightly shaky as Illumina flowcells and reagents are sold at significant profit. I imagine a large part of Illumina’s costs are related to logistical issues, rather than consumables themselves.

In any case, the patent proposes a massive cost reduction by more densely packing DNA on the flowcell. The patent mostly refers to ordered arrays, and many examples refer to a single molecule approach. The basic chemistry however seems to be pretty standard Illumina style sequencing-by-synthesis.

The figure below shows a simulation of DNA attached to a surface, at varying pitch (spacing). The right-hand images are de-convoluted versions of the left. It’s clear that as the pitch gets smaller, the image gets more crowded, and it’s harder to identify individual spots.

When imaging using a standard optical microscope, you would expect your density to be diffraction limited. Essentially, you can’t clearly identify features smaller than the wavelength of light (~200nm)… normally.

However, a number of recent techniques have broken the diffraction limit. These have allowed optical microscopes to resolve features down to 10s of nanometers. In this patent, Apton apply some “super-resolution”-like approaches… but in a limited scope (we’ll revisit what Illumina might be doing here later).

A basic super-resolution approach is shown below (not from Apton):

Each “peak” in part A of the figure above is the signal from a single fluorophore. Because the peaks are well separated we can extract each one and look at its distribution. In part B we see a single distribution. This is a 2D Gaussian. If we just took the pixel of highest intensity as the location of the fluorophore our resolution would be diffraction limited to ~200nm. However, by performing a Gaussian fit over the distribution we can determine the location to sub-pixel resolution. In this case, they could identify fluorophore locations at a final resolution of 1.5nm.

The above approach only works because the flurophores are well separated. If the Gaussians overlapped, the fit wouldn’t work. In the image above you can see the FWHM of the Gaussian is about 3 pixels, this represents ~250nm on the surface. I’d imagine if flurophores were any closer than this you’d have issues.

In their patent, Apton use the above approach to identify positions of single DNA strands on the surface to a sub-diffraction limited resolution of “10 nm RMS or less”. Apton appear to use essentially the above approach. However they have a problem, they want to pack the molecules as closely as they can to improve density. This means they are not well separated like those in the figure above.

To get round this Apton seem to use a couple of approaches. The first is that they use signals from multiple cycles to identify molecule positions. If the oligos attached to the surface have a fairly random distribution of bases (like the human genome) this should help a lot. For each molecule, you can select a cycle where it is illuminated, but none of its neighbors are. This means there is no crosstalk at this position and you should be able to get a good estimate of its position.

The figure below shows what appear to be images from a real single molecule experiment:

The right image shows positions identified from different cycles piled up (super-resolved positions I assume). In the left image they’ve used these to create another Gaussian. I would guess they then take the peak of this second Gaussian to give the final location of the molecule. This way, they can incorporate information from multiple cycles to give themselves the best estimate of the molecules location.

Using the above process I’d imagine they can get pretty good spot locations. They also mention the use of a crosstalk correction algorithm during the location identification step. But I’d imagine just filtering out bad looking Gaussian would work reasonably well.

While they may now have good spot/molecule locations this doesn’t mean they can pack DNA at ~10nm on the surface. This is because in any given cycle, there will be adjacent molecules which are fluorescing. The resulting Gaussian PSFs as imaged will overlap meaning that spots can’t be resolved. This is essentially crosstalk between adjacent spots.

Apton appear to be trying to use there super-accurate spot locations as the input to their crosstalk correction algorithm. The crosstalk correction process isn’t described in detail. But I can see that with very accurate spot locations, you can parameterize a model to which you can fit your observed signal.

This sounds great, but crosstalk seems to increase exponentially:

The examples say “molecule locations are resolvable with molecule-to-molecule distances under 200 nm” and elsewhere they say “acceptable crosstalk levels” … “occurs for pitches at or above 210 nm”.

So it seems based on this, a pitch of ~200nm is viable, but it’s not clear that you can go lower than this. This seems unfortunate, as it’s only about a quarter of the size of Illumina’s wells.

What About Illumina?

As mentioned above, super-resolution has been around for a while. In fact, the Genome analyzer 2 used super-resolution-like techniques. Cluster locations could be identified to sub-pixel resolution. Rather than just picking the “brightest” pixel in a cluster, adjacent pixel intensities could be fitted to a PSF to give a more accurate cluster location.

Illumina appear to have now filed a bunch of patents on various approaches to increasing density. One patent uses a DNA-PAINT [3] approach, which they suggest can increase the packing density to “may be less than about 20 nm”. Another describes a STED approach [4] (200nm).

There was however one patent that I found quite fascinating. This approach appears to be for an iSeq-like [6] platform where clusters sit over a CMOS image sensor. The device incorporates additional electrodes which allow an electric field to be created under the cluster.

This field can then be use to electrically deactivate a fluorophore (in one example, by pulling a quencher down on to it):

This seems very neat. Essentially it would be that you could image one cluster while quenching all its neighbors. This would remove any crosstalk. Giving you good separation of signals while retaining density. Theoretically you could push clusters very close together.

While the Illumina patents describes some interesting approaches, I didn’t see anything that looked like a real experimental setup or real datasets. So, maybe much of this is theory at the moment. I guess we’ll have to wait an see!

References and Notes

[1] Myosin V Walks Hand-Over-Hand:Single Fluorophore Imaging with1.5-nm Localization.

[2] http://www.freepatentsonline.com/10378053.pdf

[3] http://www.freepatentsonline.com/y2019/0276886.html

[4] http://www.freepatentsonline.com/y2019/0219835.html

[5] http://www.freepatentsonline.com/9193998.html

[6] I don’t see any reason why a similar setup might not be used with a normal (patterned or otherwise) flowcell with embedded electrodes. But the patent seems to focus on a iSeq-like apporach.

I wanted to download a set of mutations in SARS-CoV-2. CoV-GLUE seems to be a reasonable database of mutations in SARS-CoV-2. However the web interface doesn’t seem to have an option to download a dataset. And there isn’t a published API. So I threw together some ugly bash/awk to get what I wanted. I don’t imagine this will work for long, as the website appears to be under active development. But here are my notes anyway.

The website works off a (undocumented?) JSON API. I used the follow JSON template to get replacements (non-synonymous substitutions) which occur in 2 or more sequences:

{"multi-render":{"tableName":"cov_replacement","allObjects":false,"whereClause":"(true) and  (((num_seqs >= 2)))","rendererModuleName":"covListReplacementsRenderer","pageSize":500,"fetchLimit":500,"fetchOffset":FETCHOFFSET,"sortProperties":"-num_seqs,+variation.featureLoc.feature.name,+codon_label_int,+replacement_aa"}}

The above goes in a file called templ. I then just modify “FETCHOFFSET” using sed and download the first 4500 mutations (at the time of writing there are 4000 odd mutations. You’d want to stick this all in a loop… but I didn’t bother:

rm all
rm index.html;cp templ c;sed -i 's/FETCHOFFSET/0/g' c;wget http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/ --post-file=./c --header "Cookie: _cle=accepted" --header "Content-Type: application/json";cat index.html >> all
rm index.html;cp templ c;sed -i 's/FETCHOFFSET/500/g' c;wget http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/ --post-file=./c --header "Cookie: _cle=accepted" --header "Content-Type: application/json";cat index.html >> all
rm index.html;cp templ c;sed -i 's/FETCHOFFSET/1000/g' c;wget http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/ --post-file=./c --header "Cookie: _cle=accepted" --header "Content-Type: application/json";cat index.html >> all
rm index.html;cp templ c;sed -i 's/FETCHOFFSET/1500/g' c;wget http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/ --post-file=./c --header "Cookie: _cle=accepted" --header "Content-Type: application/json";cat index.html >> all
rm index.html;cp templ c;sed -i 's/FETCHOFFSET/2000/g' c;wget http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/ --post-file=./c --header "Cookie: _cle=accepted" --header "Content-Type: application/json";cat index.html >> all
rm index.html;cp templ c;sed -i 's/FETCHOFFSET/2500/g' c;wget http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/ --post-file=./c --header "Cookie: _cle=accepted" --header "Content-Type: application/json";cat index.html >> all
rm index.html;cp templ c;sed -i 's/FETCHOFFSET/3000/g' c;wget http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/ --post-file=./c --header "Cookie: _cle=accepted" --header "Content-Type: application/json";cat index.html >> all
rm index.html;cp templ c;sed -i 's/FETCHOFFSET/3500/g' c;wget http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/ --post-file=./c --header "Cookie: _cle=accepted" --header "Content-Type: application/json";cat index.html >> all
rm index.html;cp templ c;sed -i 's/FETCHOFFSET/4000/g' c;wget http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/ --post-file=./c --header "Cookie: _cle=accepted" --header "Content-Type: application/json";cat index.html >> all

The we extract all the mutation IDs:

awk  'BEGIN{RS="id\":\"";FS="\""}{print $1}' all > allmuts

And then fetch them from the server, this will create 4000 odd files, which we can then parse further:

#awk '{print "wget --header \"Cookie: _cle=accepted\" --header \"Content-Type: application/json\" --post-file=./info.json http://cov-glue.cvr.gla.ac.uk/gluetools-ws/project/cov/custom-table-row/cov_replacement/" $1}' allmuts  > allmuts.get

The CoV-GLUE database seems like a great resource. I hope they add a feature to download sequences/results soon. I’ve seen database results presented in a few preprints. It would be nice it those papers could also include the raw data, otherwise they’re unfortunately going to end up being difficult to replicate…

UPDATE: If you’re interested in Singular, you should probably take a look at my review of their S-1.

Company

Singular Genomics Systems was founded in 2016. Pitchbook lists them as having raise 45.5MUSD and as being at Series B. They list Coatue Management, Domain Associates, F-Prime Capital, Revelation Partners and Arch Ventures [1] as investors [3]. They are a member of JLabs [2]. There appear to be ~60 employees listed on LinkedIn.

Technology

It appears that Singular Genomics have license technology from Jingyue Ju’s lab [4]. Jingyue Ju’s lab has generated a huge amount of IP around various approaches to DNA sequencing, so this doesn’t really narrow things down very much.

Singular’s patents also describe two different optical sequencing approaches. One is a single molecule “real time” sequencing approach (the closest similar commercial platform would be PacBio). The other is a ensemble approach (with an example showing amplified DNA on beads). I’ll briefly review these two patents below. But the main takeaway is that they appear to be working on a optical approach. I suspect it’s slightly more likely that they are working on an ensemble approach (as these are more common, and easier to get working).

The ensemble approach also shows the closest to what could be real data. So let’s look at this first:

Ensemble Approach

The above figure from [6] shows one innovation they describe on the basic sequencing-by-synthesis approach. Essentially what they’re suggesting is that after flowing in your standard labelled reversible terminators you flow in a “chasing mix”. This chasing mix is a set of all 4 nucleotides, with reversible terminators, but no labels. What this means is that you give all the strands a second chance to incorporate a nucleotide. This “second chance” doesn’t give you any more signal. But it does hopefully mean that unextended strands are kept “in step”.

So called “phasing” (strands failing to incorporate a base, and getting out of step) is a major source of error in sequencing-by-synthesis. I guess the idea here, is that an unlabelled nucleotide might incorporate with better efficiency than an labelled one.

Beyond this, the patent discusses methods of speeding up imaging, potentially by taking images during the “chasing step”. This is interesting in the sense that were otherwise the imaging time would be wasted, you can use it here to help extend unextended strands, without otherwise altering the signal.

The graph above shows the first 10 cycles of a fragment of PhiX. To properly understand this data it would need normalization, but there doesn’t seem to be much in the way of phasing. This work appears to have been preformed on beads. This seems to suggest that they have something up and running. Unfortunately it doesn’t tell us much about their proposed amplification/cluster/polony generation approach. I’d guess they are using a bead based platform to evaluate the chemistry and have other ideas around amplification. But it’s also possible that they are designing a bead based platform (like Ion Torrent/454).

Single Molecule Approach

Another patent [5] discusses a single molecule approach. In this approach they’re watching a polymerase incorporate nucleotides in realtime. Here they’re suggesting detection through FRET one option appears to be have a couple of FRET acceptor/donor sites on the polymerase. As the polymerase incorporates a nucleotide a conformational change occurs and your get a FRET. You also use a label on the nucleotide to then observe incorporation using FRET.

The paper suggests observation via grating style TIRF is desirable. Where the grating could be incorporated into the flowcell. There are a few other bits and pieces of interest in the patent, such as attachment methods, but nothing that looked like an experimental setup or data to me.

Overall these patents don’t give a clear picture of where Singular Genomics is heading. It seems that the approach is likely optical, and I suspect not single molecule based on the lack of experimental data. Will be interesting to see how things develop!

Notes

[1] https://www.archventure.com/portfolio/

[2] https://jlabs.jnjinnovation.com/sites/jlabs/files/JLABSPortfolioSocial.pdf

[3] https://pitchbook.com/profiles/company/226128-25#investors

[4] https://www.cbinsights.com/company/singular-genomics

[5] Single molecule approach: https://patents.google.com/patent/US20180258472A1

[6] http://www.freepatentsonline.com/y2020/0102609.html

In order to decrease SBS cycle times, in embodiments of the present disclosure the identity of distinguishable, blocked dNTP analogue incorporated into the labeled, blocked extension product(s) generated in the sequencing reaction of such cycle is assessed while the sequencing reaction is running, i.e., before the chasing reaction is initiated. In other embodiments, such assessment is conducted during the chasing reaction. In embodiments, such assessment is conducted about less than 60 seconds before termination of the sequencing reaction, about less than 60 seconds before initiation of the chasing reaction, about less than 300 seconds after initiation of the chasing reaction, or about less than 60 to about less than 10 seconds before termination of the chasing reaction. In embodiments, such assessment is conducted substantially simultaneously with initiation of the chasing reaction. In embodiments, such assessment is conducted at the conclusion of or after the chasing reaction.

chasing conditions (i.e., conditions under which an unlabeled, blocked dNTP analogue species can be incorporated into a primed template DNA molecule that was not extended to include a distinguishable, blocked dNTP analogue species), thereby forming the unlabeled, blocked extension product(s).

DNA Sequencing with Simultaneous Imaging and Chase Steps

Provided here is an example of the embodiment of a sequencing-by-synthesis method where the DNA bases were identified during the chase step. In this example, identical DNA fragments derived from the PhiX 174 genome were immobilized on 1 micron beads. The beads were tethered to a glass coverslip which was part of a flow cell. All necessary reagents for SBS were sequentially delivered into the flow cell. At first, four distinguishable, blocked dNTP analogues were presented into the flow cell. Each dNTP was labeled with a different fluorophore as follows: dCTP-Bodipy, dTTP-R6G, dATP-AF568, dGTP-AF647. A sequencing polymerase was used to incorporate these dNTPs into the complementary strand. A small volume of buffer was then used to remove any excess dye-labeled, blocked dNTPs. As a second step, dark, blocked dNTPs were introduced into the flow cell. During this second step, a set of four images was taken, one for each of the colors corresponding to each dye-labeled dNTP, while the dark dNTPs continued to be incorporated into any unextended DNA templates on the bead. The images were obtained using a Nikon microscope, with a 20×0.75 NA objective, and standard filter sets corresponding to each of the dyes. Note that the images were taken simultaneously with the chasing step, at a temperature of 60° C., demonstrating the compatibility of the two processes in terms of reaction conditions. The excess dark, blocked dNTPs were then washed out and a deprotection reagent was brought in. This reagent cleaved the blocking group and the dye from the incorporated dNTPs. The cycle was then repeated as many times as desirable. The results from the first 10 cycles obtained using this method are shown in FIG. 1. Each bar represents the fluorescent signal from the corresponding base. Spectral cross-talk correction, which compensates for the spectral bleed through from one emission channel into another, was applied to these data. No additional corrections have been applied. We have shown sequencing read lengths of >75 bases in this manner.

http://www.freepatentsonline.com/y2019/0077726.html

http://www.freepatentsonline.com/y2020/0102609.html

http://www.freepatentsonline.com/y2019/0352508.html

https://www.indeed.com/cmp/Singular-Genomics-1/reviews

https://patents.google.com/patent/US20180258472A1/en?inventor=eli+glezer&oq=eli+glezer&sort=new&page=2