ADAPTIVE BASE CALLING SYSTEMS AND METHODS
Methods for updating a system comprising a sequencer are described herein. In some exemplary methods, the system is updated through generating sequencing data for a plurality of nucleic acid molecule colonies, selecting sequencing data for a subset of the nucleic acid molecule colonies, calling preliminary sequences for the subset of the nucleic acid colonies, mapping the called preliminary sequences to a known reference sequence, and updating the pre-trained sequencer-specific machine-learning model. Also described herein are systems for carrying out such methods and computer readable memory for storing such methods.
1 . A method of updating a system comprising a sequencer, the method comprising:
receiving, at one or more processors, sequencing data for a plurality of nucleic acid molecule colonies comprising nucleic acid molecules derived from a selected species, wherein the sequencing data was generated according to a method comprising extending sequencing primers hybridized to the nucleic acid molecules using a plurality of sequencing flow steps, each sequencing flow step comprising combining the plurality of nucleic acid molecule colonies with nucleotides, wherein at least a portion of the nucleotides are labeled, and measuring, for each nucleic acid molecule colony, a signal intensity value indicating nucleotide incorporation into the sequencing primers hybridized to said nucleic acid molecules, wherein the sequencing data comprises, for each nucleic acid molecule colony, a signal intensity value at each sequencing flow step;
selecting, using the one or more processors, sequencing data for a subset of the nucleic acid molecule colonies;
calling, using the one or more processors, preliminary sequences for the subset of the nucleic acid molecule colonies, comprising inputting the selected sequencing data into a pre-trained sequencer-specific machine-learning model configured to call a homopolymer length or a homopolymer length likelihood for each sequencing flow step based on the signal intensity values, wherein the pre-trained sequencer-specific machine-learning model was previously trained based on sequencing data previously generated using the same sequencer and nucleic acid molecules from the same selected species;
mapping, using the one or more processors, the called preliminary sequences to a known reference sequence to identify corresponding reference sequence fragments for the called preliminary sequences; and
updating, using the one or more processors, the pre-trained sequencer-specific machine-learning model based on a training data set comprising the selected sequencing data and the corresponding reference sequence fragments.
2 . The method of claim 1 , comprising generating, using the sequencer, the sequencing data.
3 . The method of claim 1 , wherein the pre-trained sequencer-specific machine-learning model was previously updated based on penultimate sequencing data generated using the same sequencer and penultimate nucleic acid molecules from the same selected species.
4 . The method of claim 3 , wherein the pre-trained sequencer-specific machine-learning model was previously updated by a method comprising:
generating the penultimate sequencing data for a plurality of penultimate nucleic acid molecule colonies comprising the penultimate nucleic acid molecules, comprising extending sequencing primers hybridized to the nucleic acid molecules using a plurality of sequencing flow steps, each sequencing flow step comprising combining the plurality of penultimate nucleic acid molecule colonies with nucleotides, wherein at least a portion of the nucleotides are labeled, and measuring, for each penultimate nucleic acid molecule colony, a signal intensity value indicating nucleotide incorporation into the sequencing primers hybridized to said nucleic acid molecules, wherein the penultimate sequencing data comprises, for each penultimate nucleic acid molecule colony, a signal intensity value at each sequencing flow step;
selecting penultimate sequencing data for a subset of the penultimate nucleic acid molecule colonies;
calling penultimate preliminary sequences for the subset of the penultimate nucleic acid molecule colonies, comprising inputting the selected penultimate sequencing data into a penultimate pre-trained sequencer-specific machine-learning model configured to call a homopolymer length or a homopolymer length likelihood for each sequencing flow step based on the signal intensity values, wherein the penultimate pre-trained sequencer-specific machine-learning model was previously trained based on antepenultimate sequencing data previously generated using the same sequencer and antepenultimate nucleic acid molecules from the same selected species;
mapping the called penultimate preliminary sequences to the known reference sequence to identify corresponding penultimate reference sequence fragments for the called penultimate preliminary sequences; and
updating the penultimate pre-trained sequencer-specific machine-learning model based on a penultimate training data set comprising the selected penultimate sequencing data and the corresponding penultimate reference sequence fragments.
5 . The method of claim 1 , wherein the pre-trained sequencer-specific machine-learning model is selected from a plurality of pre-trained sequencer-specific machine-learning models based on a quality score, wherein the plurality of pre-training sequencer-specific machine-learning models were each trained based on sequencing data generated using the same sequencer and nucleic acid molecules from the same selected species.
6 . The method of claim 1 , wherein the pre-trained sequencer-specific machine-learning model was previously initialized using sequencing data previously generated using the same sequencer and nucleic acid molecules from a different selected species.
7 . The method of claim 6 , wherein the different selected species has a smaller genome than the selected species.
8 . The method of claim 6 , wherein the different selected species is a bacterial species or a viral species.
9 . The method of claim 6 , wherein the different selected species is Escherichia coli.
10 . The method of claim 1 , wherein the selected species is a primate.
11 . The method of claim 1 , wherein the selected species is a human.
12 . The method of claim 1 , wherein the sequencer-specific machine-learning model is a neural network.
13 . The method of claim 1 , wherein the sequencer-specific machine-learning model is a convoluted neural network.
14 . The method of claim 1 , wherein the sequencing data comprises, for each nucleic acid molecule colony, a vector comprising a signal intensity value at each sequencing flow step.
15 . The method of claim 1 , wherein updating the pre-trained sequencer-specific machine-learning model based on the training data set comprises iteratively updating the pre-trained sequencer-specific machine-learning model using the same training data set until a predetermined quality control threshold is met or surpassed.
16 . The method of claim 15 , wherein the predetermined quality control threshold is a convergence threshold.
17 . The method of claim 1 , wherein updating the pre-trained sequencer-specific machine-learning model based on the training data set comprises iteratively updating the pre-trained sequencer-specific machine-learning model using the same training data set until a change in a quality control value between iterations is below a predetermined threshold.
18 . The method of claim 15 , wherein the predetermined threshold is a convergence threshold.
19 . The method of claim 1 , wherein the selected sequencing data for the subset of the nucleic acid molecule colonies is randomly selected.
20 . A method of determining a sequence of a target nucleic acid molecule, comprising:
updating a system according to the method of claim 1 , wherein the plurality of nucleic acid molecule colonies comprises a colony comprising the target nucleic acid molecule;
inputting, using the one or more processors, the sequencing data for the colony comprising the target nucleic acid molecule into the updated sequencer-specific machine-learning model; and
calling, using the one or more processors, for the target nucleic acid molecule, a homopolymer length for each sequencing flow step using the updated sequencer-specific machine-learning model.
21 - 56 . (canceled)