Flow cell alignment methods and systems

Information

  • Patent Grant
  • 10249038
  • Patent Number
    10,249,038
  • Date Filed
    Monday, January 23, 2017
    7 years ago
  • Date Issued
    Tuesday, April 2, 2019
    5 years ago
Abstract
A DNA flow cell processing method including positioning the flow cell on a stage at a predetermined location relative to a camera, illuminating the flow cell from a side with a first light source to reflect light off the DNA fragment bead locations, obtaining a first image of the flow cell and identifying a first reference pattern of bead locations in the first image, moving at least one of the flow cell and the stage relative to the camera, attempting to reposition the stage at the predetermined location, obtaining a second image of the flow cell, identifying the first reference pattern in the second image, and evaluating a first offset, relative to the camera, between the first reference pattern in the first image and the first reference pattern in the second image.
Description

This application claims the benefit of and priority to U.S. application Ser. No. 13/832,509, filed on Mar. 15, 2013, which issued as U.S. Pat. No. 9,554,095 on Jan. 24, 2017, the contents of which are incorporated herein by reference.


BACKGROUND OF THE INVENTION

Over the past 25 years, the amount of DNA sequence information generated and deposited into Genbank has grown exponentially. Many of the next-generation sequencing technologies use a form of sequencing by synthesis (SBS), wherein specially designed nucleotides and DNA polymerases are used to read the sequence of chip-bound, single-stranded DNA templates in a controlled manner. Other next-generation sequencing technologies may use native nucleotides and/or polymerases or labeled oligonucleotides and ligation enzymes to determine nucleic acid sequences. To attain high throughput, many millions of such template spots, each being either single or multiple molecules, are arrayed across a sequencing chip and their sequence is independently read out and recorded. The desire to perform high throughput sequencing stems from the need for faster processing and reduced costs. However, commercial high throughput systems, while reducing the cost of large scale sequencing (e.g., 10-100 gigabases), make smaller scale sequencing (e.g., 100 megabases to 1 gigabase) costly and inconvenient.


Recently, instruments have been developed to perform sequencing on a much smaller scale than conventional devices. Exemplary apparatus and methods that may be used for performing smaller scale sequencing operations are described in U.S. Patent Publication Nos. 2010/0323350 (application Ser. No. 12/719,469, field Mar. 8, 2010), 2010/0152050 (application Ser. No. 12/704,842, field Feb. 12, 2010), and 2009/0298131 (application Ser. No. 12/370,125, field Feb. 12, 2009). The foregoing are incorporated herein by reference. Such instruments use an “assembly line” type system, which may be arranged in the form of a carousel, to simultaneously process a number of relatively small flow cells containing single-stranded DNA templates. During operation, each flow cell is physically moved through a series of processing stations. Some of these processing stations purge the flow cell and fill it with a new reagent, while others are used for imaging the flow cell, or as idle stations where the flow cell is held without substantive processing. Other processing stations may also be provided. These instruments provide high throughput SBS operations, while offering significant savings in reagents and other processing costs. This new generation of instruments is expected to expand the public's access to SBS operations to use for various purposes, at a reduced cost, and with more rapid turnaround than earlier devices could offer.


There continues to be a need to advance the state of the art of sequencing instruments, and particularly those that use movable flow cells for small-scale sequencing operations.


SUMMARY

The present invention is embodied in imagers and alignment methods for use by imagers imaging deoxyribonucleic acid (DNA) fragments on a flow cell are disclosed. The imagers capture intensity values at DNA fragment bead locations in tiles with each tile having a reference location in the flow cell.


Flow cells may be aligned by obtaining a dark field image of each tile during a first imaging session, identifying a first dark field constellation of bead locations within a first tile and a second dark field constellation of bead locations within a second tile during the first imaging session, identifying constellations corresponding to the first and second dark field constellations during a second imaging session, altering the reference location of at least one tile during the second imaging session to correct for a linear offset in the corresponding first dark field constellations, and applying at least one correction factor for reading out intensity values from the imager for the bead locations in the flow cell to correct for an angular offset determined from offsets in the corresponding first and second dark field constellations.


An imager may include an xy stage supporting the flow cell, a camera configured to capture images of flow cell tiles, and a processor coupled to the xy stage and the camera. The processor may be configured to position the XY stage and control the camera to obtaining a dark field image of each tile during a first imaging session, identify a first dark field constellation of bead locations within a first tile and a second dark field constellation of bead locations within a second tile during the first imaging session, identify constellations corresponding to the first and second dark field constellations during a second imaging session, alter the reference location of at least one tile during the second imaging session when positioning the XY stage to correct for a linear offset in the corresponding first dark field constellations, and apply at least one correction factor for reading out intensity values from the camera for the bead locations in the flow cell to correct for an angular offset determined from offsets in the corresponding first and second dark field constellations.





BRIEF DESCRIPTION OF THE DRAWINGS

The invention is best understood from the following detailed description when read in connection with the accompanying drawings, with like elements having the same reference numerals. When a plurality of similar elements are present, a single reference numeral may be assigned to the plurality of similar elements with a small letter designation referring to specific elements. When referring to the elements collectively or to a non-specific one or more of the elements, the small letter designation may be dropped. Lines without arrows connecting components may represent a bi-directional exchange between these components. This emphasizes that according to common practice, the various features of the drawings are not drawn to scale. On the contrary, the dimensions of the various features are arbitrarily expanded or reduced for clarity. Included in the drawings are the following figures:



FIG. 1 is a block diagram of an imager in accordance with aspects of the present invention;



FIG. 2 a perspective view of an imager in accordance with aspects of the present invention;



FIG. 3 is a side view of a flow cell within the imager of FIG. 2 illustrating DNA fragment beads reflecting dark field LED light;



FIG. 4 is a top view of a flow cell within the imager of FIG. 2 illustrating virtual tiles for imaging the flow cell;



FIG. 5 is a timing diagram for imaging a flow cell in accordance with aspects of the present invention;



FIG. 6 is a flow chart of steps for imaging a flow cell in accordance with aspects of the present invention; and



FIG. 7 is a flow chart of steps for altering tile reference locations and registering DNA fragment bead locations in accordance with aspects of the present invention.





DETAILED DESCRIPTION OF THE INVENTION


FIG. 1 depicts an imager 100 for imaging bead locations of deoxyribonucleic acid (DNA) fragments on a flow cell of a chip 102 in accordance with one embodiment of the present invention. As will be discussed in further detail below, the imager 100 images the flow cell in tiles with each tile having a reference location in the flow cell. A transport 104 and an xy stage 106 position the chip 102 for imaging by the imager 100. The transport 104 moves the XY stage 106 to remove the chip from a carousel or other type of conveyor system and position the chip for imaging within the imager 100. The xy stage 106 positions the chip 102 in a horizontal plane under an objective lens 108 for imaging of individual tiles.


The imager 100 includes one or more dark field LEDs 110 and at least one bank of fluorescent LEDs 112 to illuminate the DNA fragments on the flow cell. In the illustrated embodiment, the dark field LED 110 illuminates the chip 102 from the side and the fluorescent field LEDs 112 illuminate the chip 102 from above through a lens 113, an excitation filter 114, and a dichroic mirror 116. The dichroic mirror 116 reflects light in one or more frequency bands and passes light not in those frequency band(s). The dark field LED 110 may be a red LED. The fluorescent field LEDs 112 and corresponding filters 114 may be selected to illuminate the chip 102 with desired colors of light.


When the dark field LEDs 110 are used to illuminate the chip, the mirror 116 may be removed from the pathway such that light reflected by the bead locations in the flow cell along an imaging axis 111 passes though the objective lens 108, an emission filter 109, and a lens system 109. The light then impinges upon a camera 118 where it is captured. When the fluorescent field LEDs 112 are used, light from the LEDs 112 pass through the filter 114 and are reflected by the mirror 116 through the objective lens 108 toward the chip 102. The bead locations emit light in response to being illuminated by light from the fluorescent field LEDs 112 (with the emitted light having a different wavelength). The emitted light may be produced by fluorescing dyes or other compositions associated with particular nucleotides, as described in U.S. application Ser. No. 12/719,469. The emitted light passes through the objective lens 108, the mirror 116, the emission filter 107, and the lens system 109 to impinge on the camera 118. The camera may be 3296 by 2472 charge coupled device (CCD). Although one bank of fluorescent field LEDs 112 is illustrated, additional banks of LEDs and corresponding filters may be employed to illuminate the chip with different colors of light. For example, four different colors of light may be produced, e.g., blue, green, yellow, and red.


A processor 120 is coupled to the various components to control the positioning of the chip and the imaging of the flow cell on the chip. A memory 122 is coupled to the processor 120. The processor 120 may control the imaging system 100 to implement one or more of the steps described herein. The memory is a non-transitory computer readable medium that may store instructions that, when carried out by the processor 120, implement one or more of the steps described herein. Additionally, the memory may store image information from camera 118 such as intensity values at bead locations obtained during imaging.



FIG. 2 is a perspective view of the imaging system 100 depicting the xy stage 106 with the chip 102 thereon in greater detail. The chip 102 includes a flow cell 200 containing DNA fragments. The xy stage 106 includes an x-stage 202 and corresponding controller 204 and a y-stage 206 and corresponding controller 208. The processor 120 (FIG. 1) is coupled to the x-controller 204 and the y-controller 208. The processor 120 controls the x-controller 204 to move the x-stage 200 (and thereby the flow cell) along a first axis (referred to herein as the x-axis) perpendicular to an imaging axis 111. The processor 120 additionally controls the y-controller 208 to move the y-stage 206 (and thereby the flow cell 200) along a second axis (referred to herein as the y-axis) perpendicular to the imaging axis 111 and to the first axis.



FIG. 3 depicts a cutaway of a portion a chip 102 being illuminated by the dark field LED 110. The flow cell 200 includes a transparent layer 300 including DNA fragments thereon. The transparent layer 300 may be glass. The DNA fragments reflect light from the dark field LED 110 through the objective lens 108 (FIG. 1) toward the camera 118, which identifies their location. These locations of DNA fragments are referred to herein as bead locations. FIG. 3 depicts four bead locations (302a, b, c, d) although in actual use there will be many orders of magnitude more bead locations per tile. As illustrated in FIG. 3, the dark field LED 110 may be tilted to optimize illumination of the beads. In an embodiment, the objective lens 108 and the dark field LED 112 are stationary in the x,y plane so that the bead illumination incidence angle does not change.



FIG. 4 depicts a top view of a chip 102. The chip includes a flow cell 200. During imaging, the flow cell 200 is imaged in tiles 400. In the illustrated embodiment, the flow cell 200 is imaged using 72 tiles (three rows of 24 tiles). Each tile has a reference location corresponding to the position of the tile on the chip. As will be described in further detail below, at least two tiles, e.g., tile 400a and 400b, may be selected for use in aligning the flow cell 200 and registering image intensities captured for individual tiles 400 within the flow cell 200 during imaging. The at least two tiles may be located essentially anywhere within the flow cell. In an embodiment, the tiles are in a center row and are at least 5 tiles away from one another, e.g., 9 tiles away from one another.


A reference constellation of bead locations within each tile, e.g., 402a and 402b, may be used for alignment and registration purposes. In one embodiment, one reference constellation 402a is used to determine an offset for the tile reference locations when the flow cell has been removed from the imager 100 and added back in for a subsequent imaging session. Two or more reference constellations may be used in other embodiments. Additionally, the difference in the offsets for two constellations may be used to calculate an angular adjustment for bead locations when reading out intensity values.



FIG. 5 depicts a timing diagram 500 and FIG. 6 depicts a flow chart 600 for imaging a flow cell in accordance with aspects of the present invention. Although the timing diagram and flow chart are described with reference to FIGS. 1-4, other suitable apparatus for carrying out the steps described herein will be understood by one of skill in the art from the description herein and are considered within the scope of the present invention. Additionally, one or more of the steps described herein may be omitted or performed in a different order without departing from the spirit and scope of the present invention. The total time for imaging all tiles in a flow cell may be on the order of 6 minutes.


At step 602, a flow cell (FC) is lifted from a flow cell carrier such as a carousel. The flow cell may be lifted by a transport 104 into position under an objective lens 108 of an imager 100 for imaging during an imaging sequence. It may take approximately 7 seconds for the transport to move the flow cell from the carrier to the imager 100 for imaging.


At block 604, the flow cell is aligned within the imager. If it is the first time the flow cell is positioned within the imager, a dark field image of a first tile, e.g., tile 400a, and a second tile, e.g., tile 400b, is obtained and a subset of the beads in each image is identified (referred to herein as a constellation; e.g., constellations 402a and 402b). Thereafter, in subsequent imaging sessions, a dark field image of the first tile and the second tile is obtained and the dark field images from the current session are compared to the dark field images from the first imaging session to identify matching constellations and to determine a corresponding difference in position. The difference in position between the constellations from at least one tile may then be used to alter the reference positions of each tile for the current imaging session. Additionally, dark field images for all the tiles may be obtained at this stage. It may take approximately 29 seconds or more to align the stage and obtain the dark field images.


At block 606, a filter and mirror are positioned for fluorescent field imaging. The processor 102 may control the positioning of the filter 114 and mirror 116, e.g., though an additional transport device (not shown). It may take approximately 4 seconds to move the filter and the mirror.


At block 608, the flow cell is positioned for imaging of a particular tile, e.g., tile 0 located in the upper left hand corner of the tiles 400 (FIG. 4). The processor 102 may position the flow cell for imaging through the use of the x-controller 204 and the y-controller 208 (FIG. 2) to move the flow cell to the reference position corresponding to the desired tile as altered in block 604. It may take approximately 7 seconds to position the flow cell at the first tile. Thereafter, it may take a little over 0.5 second to move to the next tile to be imaged.


At block 610, the imaging system 100 autofocuses (ATF) to create a sharp image of the flow cell tile on the camera. The processor 120 may perform autofocussing by adjusting the position of the camera 118, objective lens 108, or the chip 102 (e.g., via the transport 104 (FIG. 1)) along the imaging axis 111. It may take approximately 0.25 seconds to perform autofocussing.


At block 612, the fluorescent LEDs are turned on. The processor 120 may turn on the fluorescent LEDs 112. It may take approximately 120 milliseconds to turn on the LEDs 112.


At block 614, the camera is triggered (snap) and, at block 616, the imaging system is synchronized. The processor 120 may trigger the camera 118 to take an image and synchronize itself with the camera 118. It may take approximately 110 milliseconds to trigger the camera 118 and approximately 200 milliseconds to synchronize.


At block 618, the camera images the fluorescent light from the flow cell tile (exposure time). The camera 118 may have an exposure period of 350 milliseconds during which light from the flow cell is captured.


At block 620, an optional conversion step is performed to covert the imaged from a 1-dimensional (1D) array to a 2-dimensional (2D) array. Processor 120 may perform the conversion. This step may be performed as a programming convenience and may be omitted. It may take approximately 43 milliseconds to perform the conversion.


At block 622, a decision is made regarding tiles to be processed. If there are additional tiles to process, processing proceeds at block 608 with the next tile and, in parallel, processing proceeds at block 624 with further processing of the light captured from the current tile. If there are no additional tiles, processing proceeds only at block 624.


At block 624, the image is downloaded from the camera. The processor 120 may download the image from the camera 118, e.g., into random access memory (RAM) associated with the processor. It may take approximately 700 milliseconds to download the image.


At block 626, post processing of the downloaded image is performed. Processor 120 may perform the processing and store the result in memory 122. The post processing may involve registration of the captured images with the bead locations determined from the dark field image as adjusted based on offsets determined at block 604. It may take approximately 365 milliseconds to register the images with the bead locations and approximately 25 milliseconds to extract the pixels, e.g., for storage in memory 122.



FIG. 7 depicts a flow chart 700 of steps describing the alignment of flow cells and the registration of captured images in further detail in accordance with aspects of the present invention. Although flow chart 700 is described with reference to FIGS. 1-4, other suitable apparatus for carrying out the steps described herein will be understood by one of skill in the art from the description herein and are considered within the scope of the present invention. Additionally, one or more of the steps described herein may be omitted or performed in a different order without departing from the spirit and scope of the present invention.


At block 702, a flow cell is positioned for a first imaging session. A flow cell 200 on a chip 102 may be lifted by a transport 104 into position under an objective lens 108 of an imager 100 under control of processor 120 for imaging during the first imaging session.


At block 704, a dark field image is obtained and intensity values are read out for each tile in the flow cell to be imaged. In one embodiment, a dark field image is first obtained for two reference tiles 400 within the flow cell 200 and a dark field image for all tiles 400 is subsequently obtained. In another embodiment, a dark field image for all tiles 400 is initially obtained without separately obtaining another dark field image for the two reference tiles.


The dark field image for each tile 400 may be obtained by moving the flow cell 200 to a first reference location associated with a first tile using an xy-stage 106 under control of the processor 120, illuminating the side of the flow cell 200 with a dark field LED 110, capturing light reflected by the DNA fragments with a camera 118, and storing intensity values in a memory 122 (with relatively high intensity values representing the locations of the DNA fragment bead locations). Subsequent tiles within the flow cell may then be imaged by iteratively moving the flow cell to a respective reference location associated with each of those tiles and repeating the capturing process. During the first imaging session for each flow cell, the reference locations for the tiles may be nominal locations established based on a standard flow geometry and xy-stage 106 positions. For example, the xy-stage 106 may start at a “home” location at one limit of the x and y travel range, and each nominal location may be predetermined as an x and y distance from the home location that is expected to place the center of each tile directly under the objective lens 108. In this example, tile 0 may have a nominal location stored in a lookup table as (x0, y0), where x0 is a distance from the home x location, and y0 is a distance from the home y location, and each other tile may have its own nominal location stored in a similar manner (distances may be measured as distances, motor movement increments, and so on).


At block 706, bead locations are determined. The locations within the dark field image having a relatively high intensity value, represent the bead locations within the flow cells. The bead locations (representing the location of the DNA fragments) may be determined by analyzing the intensity values from the dark field image with the processor 120 and selecting those with an intensity value above a value indicating the presence of a bead. After each bead is located in the dark field, the x and y locations of the bead are stored and used in later imaging steps as expected bead locations. For example, an array of x,y locations corresponding to each bead may be saved for each tile. Having identified the bead locations in terms of x and y locations, the full-resolution image itself may optionally be discarded.


At block 708, dark field reference constellations are identified within at least two tiles. In an embodiment, the processor 120 identifies two dark field reference constellations—a first dark field reference constellation within a first tile and a second dark field reference constellation within a second tile. Each dark field reference constellation includes the identified bead locations located within at least a portion of its respective tile. In an embodiment, the dark field reference constellation of a tile is the bead locations in an approximately 500 by 500 pixel area in the center of the tile. The dark field reference constellations may be determined at the time the dark field image is obtained for the two reference tiles (i.e., in step 704), or they may be identified in a separate imaging step.


At block 710, a fluorescent field image is obtained for each tile in the flow cell to be imaged. In an embodiment, a fluorescent field image is obtained by moving the flow cell 200 to the first nominal reference location associated with a first tile using an xy-stage 106 under control of the processor 120, illuminating the flow cell 200 with a filtered fluorescent field LED 112, and capturing light emitted by the DNA fragments with a camera 118. Subsequent tiles within the flow cell may then be imaged by iteratively moving the flow cell to the respective nominal reference location associated with each of those tiles, and repeating the capturing process.


At block 712, correction factors are determined and applied. The processor 120 may determine the correction factors and apply them to the expected bead locations from the dark field image to obtain corrected bead locations.


As the xy-stage moves the flow cell for imaging the individual tiles, error may be introduced by the physical movement from tile to tile. This error may be accommodated using software to compare a test region near the center of the fluorescent field image of each tile with the dark field image of the same tile to determine an offset, and incorporating the offset into a correction factor when reading out the intensity values at the bead locations. The test region may be a 500 by 500 pixel area at the center of the fluorescent field tile image. The process may be conducted by identifying a number of x, y bead locations determined in the dark field imaging step, and adding up the intensities of the same locations in the 500 by 500 pixel fluorescent field test region to obtain a total intensity value. Next, the 500 by 500 pixel test region is iteratively shifted by a predetermined amount, e.g., ±10 pixels in the x and y directions, and the intensities at the x and y bead locations are once again measured and summed up for each iteration. Ultimately, the test region will be shifted to a point where the intensities of the fluorescent light emitted from the beads will overlie the x,y bead locations identified in the dark field imaging step, and the sum of the intensities will reach a peak value. At this point, the offset is determined as the total x and y shift distances necessary to obtain the peak intensity. These x and y values are used as a correction factor to shift the fluorescent field tile image data and thereby align the intensities across the entire fluorescent field image with the bead locations from the dark field image.


Additional correction factors may be determined. These other factors may include, for example, “pincushion” correction (correcting to account for a different diffraction angle of the particular wavelength as compared to the dark field LED), predetermined distortion corrections for known lens aberrations, and so on.


At block 714, the intensity values of the corrected fluorescent field image for the tile at the corrected bead locations are read out. The processor 120 may read out the intensity values at the corrected bead locations.


If fluorescent field images for other colors are to be determined, processing may return to block 710 with steps 710-714 repeated for each additional color. The flow cell is thereafter removed from the imaging station.


At block 716, the flow cell is repositioned for a second imaging session. A flow cell 200 on the chip 102 may be lifted by the transport 104 into position under an objective lens 108 of an imager 100 for imaging during the second imaging session.


At block 718, dark field constellations corresponding to the dark field reference constellations identified for the flow cell during the first imaging session (block 708) are determined. In an embodiment, the processor 120 identifies two dark field constellations corresponding to the dark field referenced constellations identified during the first imaging session—a first dark field constellation within the first tile and a second dark field constellation with a second tile. To do so, the xy-stage 106 moves the flow cell 200 to the nominal tile locations for each of the two tiles in which the constellations should be found and takes a dark field image of each. The dark field constellation within each tile may be the bead locations located within a 500 by 500 pixel area at the geometric center of the tile.


At block 720, offsets between the dark field constellations in the second imaging session and the dark field reference constellations in the first imaging session are determined. The offset between the constellations for a tile may be determined by iteratively shifting the dark field constellation in the second imaging session by a predetermined amount, e.g., ±10 pixels in the x and y direction, and evaluating the intensities of this dark field image at the x,y locations corresponding to the bead locations found in the original reference constellation. The offset can be determined by summing the intensities of the dark field image at the expected x,y bead positions determined from the dark field reference constellation at each shift iteration, with the shift iteration producing the maximum intensity at the expected bead positions representing the offset. This may be essentially the same process as described above in step 712.


In one embodiment, a single dark field constellation is compared to a single reference constellation, and the offset between the two may be used to establish a linear offset of the flow cell in the x and/or y directions introduced when the flow cell was repositioned for a subsequent imaging. In embodiments using a second constellation as a second reference, the system can also calculate an angular offset that was introduced when the flow cell was repositioned for the subsequent imaging. For example, on each subsequent placement in the imaging station, the new location of the first constellation 402a may be determined and used as a first offset reference point, and the new location of the second constellation 402b may be determined and used as a second offset reference point. By comparing the first and second offset reference points to the original constellation locations, the total x,y offset and angular rotation of the flow cell are determined.


At block 722, new reference locations for the tiles are calculated to accommodate for any shift in the flow cell's position as compared to its position in the first round of imaging. In one embodiment, the processor 120 uses the x,y offset from one reference constellation to calculate a geometric correction factor for each of the nominal tile center locations that were used in step 704. For example, if the constellation found in step 718 was offset by +50 x pixels and +100 y pixels from the referenced constellation location, each nominal tile center location would be offset a corresponding −50 x and −100 y pixels, and the xy-stage would be driven to these recalculated locations to place the center of each tile directly (or nearly directly) under the objective lens for further imaging steps.


In another embodiment, the processor 120 uses the x,y offset from two reference constellations to calculate a geometric correction factor for each of the nominal tile center locations that were used in step 704. In this case, the processor 120 may calculate x, y and angular correction factors for each nominal tile center location to accurately position the center of each tile under the objective lens 108. The shape of the flow cell remains constant between successive visits to the imaging system, and therefore the calculation of the geometric correction factors is a matter of basic geometry that need not be described herein. This process is expected to be more accurate than simply using an x,y offset as in the prior embodiment, but at the cost of greater processing time to evaluate a second constellation. As with the prior embodiment, the new reference locations are used to drive the xy-stage 106 to place the center of each tile directly (or nearly directly) under the objective lens 108 for the following steps.


At block 724, fluorescent field images are obtained for each tile during the subsequent imaging session. The fluorescent field images may be determined as described above with respect to block 710 using the reference locations as altered at block 722.


At block 726, correction factor(s) are applied to the fluorescent field image. The corrections factors applied in step 726 may include corrections as described above with reference to block 712 (e.g., registering the image to match the dark field image, pincushion and distortion correction, etc.). In addition, a correction factor also may be applied to rotate the fluorescent field image to correct for any angular offset identified in step 720. Following the corrections, the fluorescent field tile image should be properly shifted in the x,y direction and rotated to match the illuminated points and regions to the x,y locations of the beads determined from the original dark field imaging step for the respective tile.


At block 728, the intensity values at the corrected bead locations are read out. The processor 120 may read out the intensity values at the corrected bead locations.


If fluorescent field images for other colors are to be determined, processing may return to block 724 with steps 724-728 repeated for each additional color. The flow cell is thereafter removed from the imaging station. If the flow cell is to be imaged an additional time after further processing, processing may resume at block 716.


It will be appreciated from the foregoing example that embodiments of the invention can be used to help resolve a significant issue relating to sequencing by synthesis (“SBS”) instruments that use multiple flow cells that are repeatedly mounted on and dismounted from an imaging system. By its nature, the SBS process requires repeated examination of exactly the same nucleotide strand to properly identify the nucleotides as they join the sequence. In conventional systems, fixed flow cells are used. Such systems have very little free play between the flow cell and the camera, so it is relatively easy to place the camera at almost exactly the same location during each successive imaging step. In contrast, in systems that use flow cells that are repeatedly dismounted and mounted on the imaging platform, there is much more free play. The free play in such devices may derive from clearances necessary to automate the mounting and dismounting process, manufacturing variances between one flow cell and the next, and so on. In these systems, it has been found that it is very difficult to properly align the flow cell in exactly the same location during successive visits to the imaging station. Since the beads being examined may be only a few hundred microns in diameter (or smaller), even a small amount of misalignment will make it impossible to correlate the data obtained from successive imaging steps, potentially rendering the device useless. The above process and versions thereof can be used to overcome this problem, and provide rapid image acquisition and processing in SBS instruments that use removable flow cells.


Although the invention is illustrated and described herein with reference to specific embodiments, the invention is not intended to be limited to the details shown. Rather, various modifications may be made in the details within the scope and range of equivalents of the claims and without departing from the invention.

Claims
  • 1. A method for processing flow cells containing deoxyribonucleic acid (DNA) fragments, the method comprising: positioning a flow cell containing DNA fragment beads on an xy stage with the flow cell facing an imaging axis of a camera, wherein the xy stage is movable relative to the camera along an x-axis and a y-axis, the y-axis is perpendicular to the x-axis, and the imaging axis is perpendicular to the x-axis and the y-axis;positioning the xy stage at a first predetermined nominal location along the x-axis and the y-axis with respect to the camera;illuminating the flow cell from a side of and above the flow cell with a dark field light source to cause a first light to strike the DNA fragment beads at an angle relative to the imaging axis, and reflect off the DNA fragment bead locations;obtaining a first image of the flow cell;evaluating the first light reflected from the DNA fragment beads to identify, in the first image, a first reference pattern of bead locations;moving at least one of the flow cell and the xy stage relative to the camera;moving the xy stage to reposition the xy stage at a second position substantially at the first predetermined nominal location along the x-axis and the y-axis with respect to the camera;obtaining a second image of the flow cell;identifying, in the second image, the first reference pattern of bead locations;evaluating a first offset, relative to the camera, between the first reference pattern in the first image and the first reference pattern in the second image; andusing the first offset to perform a first correction.
  • 2. The method of claim 1, further comprising using the first offset to perform a first correction comprising at least one of: moving the xy stage to reposition the flow cell relative to the camera, andapplying a correction factor for reading out intensity values obtained during a fluorescent imaging process.
  • 3. The method of claim 2, wherein moving at least one of the flow cell and the xy stage relative to the camera comprises removing the flow cell from the xy stage and replacing the flow cell on the xy stage.
  • 4. The method of claim 2, wherein moving at least one of the flow cell and the xy stage relative to the camera comprises moving the xy stage in one or both of the x-axis and the y-axis.
  • 5. The method of claim 2, wherein: the first offset comprises a linear offset between a first position of the flow cell relative to the camera during the step of positioning the xy stage at the first predetermined nominal location, and the second position; andperforming the first correction further comprises correcting for the linear offset by moving the xy stage to place the flow cell at the first position.
  • 6. The method of claim 2, wherein the first offset comprises a linear offset, and applying a correction factor for reading out intensity values obtained during a fluorescent imaging process comprises correcting for the linear offset.
  • 7. The method of claim 2, wherein the first offset comprises an angular offset, and applying a correction factor for reading out intensity values obtained during a fluorescent imaging process comprises correcting for the angular offset.
  • 8. The method of claim 2, further comprising: positioning the xy stage at a second predetermined nominal location along the x-axis and the y-axis with respect to the camera;obtaining a third image of the flow cell;evaluating the first light reflected from the DNA fragment beads to identify, in the third image, a second reference pattern of bead locations;after the step of moving at least one of the flow cell and the xy stage relative to the camera, moving the xy stage to reposition the xy stage at a fourth position along the x-axis and the y-axis with respect to the camera;obtaining a fourth image of the flow cell;identifying, in the fourth image, the second reference pattern of bead locations; andevaluating a second offset, relative to the camera, between the second reference pattern in the third image and the second reference pattern in the fourth image.
  • 9. The method of claim 8, further comprising using the second offset to perform a second correction comprising at least one of: moving the xy stage to reposition the flow cell relative to the camera, andapplying a second correction factor for reading out intensity values obtained during the fluorescent imaging process.
  • 10. The method of claim 9, wherein the first correction comprises a linear correction and the second correction comprises an angular correction.
  • 11. The method of claim 8, wherein the first offset comprises a linear offset, and the second offset comprises an angular offset.
  • 12. The method of claim 8, wherein the first offset comprises a first linear offset, the second offset comprises a second linear offset, and the method further comprises: evaluating the first linear offset and the second linear offset to determine an angular offset.
  • 13. The method of claim 12, further comprising using the angular offset to perform a second correction comprising at least one of: moving the xy stage to reposition the flow cell relative to the camera, andapplying an angular correction factor for reading out intensity values obtained during the fluorescent imaging process.
  • 14. The method of claim 8, wherein: the flow cell comprises a plurality of tiles, each tile having a respective reference location;the first image and the second image are made at a first tile; andthe third image and the fourth image are made at a second tile.
  • 15. The method of claim 14, wherein the second tile is spaced from the first tile.
  • 16. The method of claim 15, wherein the first offset comprises a first linear offset, the second offset comprises a second linear offset, and the method further comprises: evaluating the first linear offset and the second linear offset to determine an angular offset.
  • 17. The method of claim 16, further comprising using the angular offset to perform a second correction comprising at least one of: moving the xy stage to reposition the flow cell relative to the camera, andapplying an angular correction factor for reading out intensity values obtained during the fluorescent imaging process.
  • 18. The method of claim 1, further comprising: obtaining a third image of the flow cell;evaluating the first light reflected from the DNA fragment beads to identify, in the third image, a second reference pattern of bead locations;obtaining a fourth image of the flow cell;identifying, in the fourth image, the second reference pattern of bead locations; andevaluating a second offset, relative to the camera, between the second reference pattern in the third image and the second reference pattern in the fourth image.
  • 19. The method of claim 18, further comprising evaluating the first offset and the second offset to determine an angular offset.
  • 20. The method of claim 19, further comprising: using the first offset to perform a first correction comprising at least one of: moving the xy stage to reposition the flow cell relative to the camera, andapplying a correction factor for reading out intensity values obtained during a fluorescent imaging process; andusing the angular offset to perform a second correction comprising at least one of: moving the xy stage to reposition the flow cell relative to the camera, andapplying an angular correction factor for reading out intensity values obtained during a fluorescent imaging process.
  • 21. The method of claim 18, further comprising: using the first offset to perform a first linear correction comprising at least one of: moving the xy stage to reposition the flow cell relative to the camera, andapplying a linear correction factor for reading out intensity values obtained during a fluorescent imaging process; andusing the second offset to perform an angular correction comprising applying an angular correction factor for reading out intensity values obtained during a fluorescent imaging process.
US Referenced Citations (300)
Number Name Date Kind
5374395 Robinson Dec 1994 A
5502773 Tibbetts Mar 1996 A
5695934 Brenner Dec 1997 A
5720923 Haff Feb 1998 A
5756334 Perler May 1998 A
5863722 Brenner Jan 1999 A
5900481 Lough May 1999 A
5981956 Stern Nov 1999 A
6136543 Anazawa Oct 2000 A
6185030 Overbeck Feb 2001 B1
6207960 Stern Mar 2001 B1
6222624 Yonezawa Apr 2001 B1
6223128 Allex Apr 2001 B1
6242193 Anazawa Jun 2001 B1
6255475 Kwiatkowski Jul 2001 B1
6309836 Kwiatkowski Oct 2001 B1
6337188 Head Jan 2002 B1
6404907 Gilchrist Jun 2002 B1
6448066 Wheatcroft Sep 2002 B1
6576424 Fodor Jun 2003 B2
6596483 Choong Jul 2003 B1
6597000 Stern Jul 2003 B2
6611767 Fiekowsky Aug 2003 B1
6613513 Parce Sep 2003 B1
6639088 Kwiatkowski Oct 2003 B2
6646243 Pirrung Nov 2003 B2
6664079 Ju Dec 2003 B2
6867851 Blumenfeld Mar 2005 B2
6890741 Fan May 2005 B2
6929944 Matson Aug 2005 B2
6995841 Scott et al. Feb 2006 B2
7033754 Chee Apr 2006 B2
7037687 Williams May 2006 B2
7052847 Korlach May 2006 B2
7057026 Barnes Jun 2006 B2
7085651 Yasuda Aug 2006 B2
7110618 Bobrov Sep 2006 B2
7145645 Blumenfeld Dec 2006 B2
7157228 Hashmi Jan 2007 B2
7169560 Lapidus Jan 2007 B2
7200254 Kira Apr 2007 B2
7233393 Tomaney Jun 2007 B2
7270951 Stemple Sep 2007 B1
7279563 Kwiatkowski Oct 2007 B2
7282337 Harris Oct 2007 B1
7306918 Hashmi Dec 2007 B2
7315637 Xi Jan 2008 B2
7329496 Dower Feb 2008 B2
7344865 Parce Mar 2008 B2
7386399 Izmailov Jun 2008 B1
7387891 Boege Jun 2008 B2
7405823 Tomaney Jul 2008 B2
7414116 Milton Aug 2008 B2
7425416 Hashmi Sep 2008 B2
7427673 Balasubramanian Sep 2008 B2
7476504 Turner Jan 2009 B2
7482120 Buzby Jan 2009 B2
7482167 Sammak Jan 2009 B2
7491498 Lapidus Feb 2009 B2
7508516 Sugiyama Mar 2009 B2
7541444 Milton Jun 2009 B2
7544794 Benner Jun 2009 B1
7566537 Balasubramanian Jul 2009 B2
7566538 Parce Jul 2009 B2
7567695 Frost Jul 2009 B2
7592435 Milton Sep 2009 B2
7635562 Harris Dec 2009 B2
7635565 Hashmi Dec 2009 B2
7635578 Ju Dec 2009 B2
7635588 King Dec 2009 B2
7666593 Lapidus Feb 2010 B2
7680790 Indeck Mar 2010 B2
7687260 Gutekunst Mar 2010 B2
7713698 Ju May 2010 B2
7714303 Lundquist May 2010 B2
7745129 Schatz Jun 2010 B1
7771973 Milton Aug 2010 B2
7785796 Balasubramanian Aug 2010 B2
7795424 Liu Sep 2010 B2
7816503 Milton Oct 2010 B2
7835871 Kain Nov 2010 B2
7843567 Moon Nov 2010 B2
7883869 Ju Feb 2011 B2
7888073 Densham Feb 2011 B2
7897345 Lapidus Mar 2011 B2
7993895 Eid Aug 2011 B2
7994304 Siddiqi Aug 2011 B2
7998717 Eid Aug 2011 B2
8012690 Berka Sep 2011 B2
8018593 Tan Sep 2011 B2
8022194 Piepenburg Sep 2011 B2
8030466 Shin Oct 2011 B2
8031918 Roth Oct 2011 B2
8032310 Stenger Oct 2011 B2
8053192 Bignell Nov 2011 B2
8058414 Menchen Nov 2011 B2
8071346 Eid Dec 2011 B2
8071739 Milton Dec 2011 B2
8084590 Milton Dec 2011 B2
8088575 Ju Jan 2012 B2
8094312 Ulmer Jan 2012 B2
8116988 Glick Feb 2012 B2
8143599 Feng Mar 2012 B2
8148503 Litosh Apr 2012 B2
8158346 Balasubramanian Apr 2012 B2
8178360 Barnes May 2012 B2
8182993 Tomaney May 2012 B2
8209130 Kennedy Jun 2012 B1
8212020 Benner Jul 2012 B2
8236532 Ronaghi Aug 2012 B2
8241573 Banerjee Aug 2012 B2
8247177 Smith Aug 2012 B2
8263364 Williams Sep 2012 B2
8263365 Williams Sep 2012 B2
8271206 Liu Sep 2012 B2
8280640 Levin Oct 2012 B2
8298792 Ju Oct 2012 B2
8299226 Piepenburg Oct 2012 B2
8300971 Osher Oct 2012 B2
8301394 Chen Oct 2012 B2
8370079 Sorenson Feb 2013 B2
8383039 Zhou Feb 2013 B2
8394586 Balasubramanian Mar 2013 B2
8399188 Zhao Mar 2013 B2
8399196 Hoser Mar 2013 B2
8407012 Erlich Mar 2013 B2
8407554 Kermani Mar 2013 B2
8427637 Gao Apr 2013 B2
8429108 Eusterbrock Apr 2013 B2
8460865 Chee Jun 2013 B2
8463718 Ben-Hur Jun 2013 B2
8481259 Gordon Jul 2013 B2
8481903 Triener Jul 2013 B2
8594951 Homer Nov 2013 B2
8612161 Gordon Dec 2013 B2
8666678 Davey Mar 2014 B2
8703422 Tomaney Apr 2014 B2
8725425 Heiner May 2014 B2
8728399 Jaffe May 2014 B2
8738300 Porreca May 2014 B2
8748789 Triener Jun 2014 B2
8774494 Staker Jul 2014 B2
8795971 Kersey Aug 2014 B2
8812268 Scoullar Aug 2014 B2
8831316 Tyurina Sep 2014 B2
8852910 Smith Oct 2014 B2
8886535 Chong Nov 2014 B2
8900810 Gordon Dec 2014 B2
8929630 Fu Jan 2015 B2
8934098 Cox Jan 2015 B2
8940481 Gordon Jan 2015 B2
8965076 Garcia Feb 2015 B2
20010030290 Stern Oct 2001 A1
20020004204 OKeefe Jan 2002 A1
20020110899 Wheatcroft Aug 2002 A1
20030087289 Zuzan May 2003 A1
20030128894 Bobrov et al. Jul 2003 A1
20030215816 Sundararajan Nov 2003 A1
20030224419 Corcoran Dec 2003 A1
20040002073 Li Jan 2004 A1
20040029126 Bloecker Feb 2004 A1
20040036731 Ready Feb 2004 A1
20040142347 Stockwell Jul 2004 A1
20040152108 Keith Aug 2004 A1
20040175718 Chee Sep 2004 A1
20040185466 Ju Sep 2004 A1
20040229269 Hashmi Nov 2004 A1
20050019902 Mathies Jan 2005 A1
20050026163 Sundararajan Feb 2005 A1
20050052509 Gilligan Mar 2005 A1
20050064469 Schulz Mar 2005 A1
20050170367 Quake Aug 2005 A1
20050211874 Takeyama Sep 2005 A1
20050227231 Tcherkassov Oct 2005 A1
20050239085 Buzby Oct 2005 A1
20050244827 Olsson Nov 2005 A1
20060034543 Bacus Feb 2006 A1
20060073493 Fasulo Apr 2006 A1
20060105349 Ekenberg May 2006 A1
20060147935 Linnarsson Jul 2006 A1
20060166203 Tooke Jul 2006 A1
20060240439 Smith Oct 2006 A1
20060275782 Gunderson Dec 2006 A1
20060281109 BarrOst Dec 2006 A1
20070048748 Williams Mar 2007 A1
20070070349 Harris Mar 2007 A1
20070076983 Doerrer Apr 2007 A1
20070082358 Fuerst Apr 2007 A1
20070117103 Buzby May 2007 A1
20070117104 Buzby May 2007 A1
20070177799 Tyurina Aug 2007 A1
20070196832 Efcavitch Aug 2007 A1
20070196846 Hanzel Aug 2007 A1
20080038839 Linder Feb 2008 A1
20080087826 Harris Apr 2008 A1
20080088823 Harris Apr 2008 A1
20080103053 Siddiqi May 2008 A1
20080138804 Buzby Jun 2008 A1
20080161194 Turner Jul 2008 A1
20080161195 Turner Jul 2008 A1
20080176241 Eid Jul 2008 A1
20080182757 Heiner Jul 2008 A1
20080246949 Harris Oct 2008 A1
20080280773 Fedurco Nov 2008 A1
20080293071 Gelfand Nov 2008 A1
20090056390 Falkensson Mar 2009 A1
20090099027 Greiner Apr 2009 A1
20090155793 Oliphant Jun 2009 A1
20090186771 Siddiqi Jul 2009 A1
20090263791 Ju Oct 2009 A1
20090298131 Gordon Dec 2009 A1
20090325260 Otto Dec 2009 A1
20100009871 Reed Jan 2010 A1
20100015611 Webster Jan 2010 A1
20100035253 Gordon Feb 2010 A1
20100035273 Axen Feb 2010 A1
20100036110 Xie Feb 2010 A1
20100041029 Ju Feb 2010 A1
20100063742 Hart Mar 2010 A1
20100063743 Gordon Mar 2010 A1
20100087327 Weng Apr 2010 A1
20100121582 Pan May 2010 A1
20100143932 Lapidus Jun 2010 A1
20100152050 Gordon Jun 2010 A1
20100153420 Yang Jun 2010 A1
20100159531 Gordon Jun 2010 A1
20100167413 Lundquist Jul 2010 A1
20100227327 Xie Sep 2010 A1
20100235105 Volkov Sep 2010 A1
20100246977 Fu Sep 2010 A1
20100311061 Korlach Dec 2010 A1
20100323350 Gordon Dec 2010 A1
20100330569 Olejnik Dec 2010 A1
20110014611 Ju Jan 2011 A1
20110045489 Gardner Feb 2011 A1
20110049389 Kishima Mar 2011 A1
20110052446 Hirano Mar 2011 A1
20110080487 Venkataraman et al. Apr 2011 A1
20110081647 Siddiqi Apr 2011 A1
20110087016 Suo Apr 2011 A1
20110124054 Olejnik May 2011 A1
20110212437 Emig Sep 2011 A1
20110220775 Triener Sep 2011 A1
20110251828 Scoullar Oct 2011 A1
20110254969 Tyurina Oct 2011 A1
20110256631 Tomaney Oct 2011 A1
20110257889 Klammer Oct 2011 A1
20110268347 Staker Nov 2011 A1
20110300534 Chiou Dec 2011 A1
20120020537 Garcia Jan 2012 A1
20120041727 Mishra Feb 2012 A1
20120052490 Eid Mar 2012 A1
20120092477 Kawano Apr 2012 A1
20120095201 Milton Apr 2012 A1
20120109598 Davey May 2012 A1
20120116688 Mishra May 2012 A1
20120156680 Ju Jun 2012 A1
20120173159 Davey Jul 2012 A1
20120197623 Homer Aug 2012 A1
20120252010 Balasubramanian Oct 2012 A1
20120258869 Kersey Oct 2012 A1
20120270305 Reed Oct 2012 A1
20120282708 Corbett Nov 2012 A1
20120330566 Chaisson Dec 2012 A1
20130053280 Hamasaki Feb 2013 A1
20130054171 Chen Feb 2013 A1
20130060482 Sikora Mar 2013 A1
20130096015 Ju Apr 2013 A1
20130116128 Shen May 2013 A1
20130116130 Fu May 2013 A1
20130122485 Hong May 2013 A1
20130137091 Gordon May 2013 A1
20130137110 Kraihanzel May 2013 A1
20130138358 Tang May 2013 A1
20130217006 Sorenson Aug 2013 A1
20130235388 Segale Sep 2013 A1
20130260372 Buermann Oct 2013 A1
20130281321 Mann Oct 2013 A1
20130289921 Gopalan Oct 2013 A1
20130301888 Gordon Nov 2013 A1
20130303384 Rearick Nov 2013 A1
20140051584 Davey Feb 2014 A1
20140073514 Shen Mar 2014 A1
20140073517 Zhou Mar 2014 A1
20140220558 Homer Aug 2014 A1
20140222399 Davey Aug 2014 A1
20140235461 Yin Aug 2014 A1
20140248618 Shaikh Sep 2014 A1
20140267669 Stoops Sep 2014 A1
20140303016 Tomaney Oct 2014 A1
20140309143 Kersey Oct 2014 A1
20140315328 Lee Oct 2014 A1
20140329694 Buermann Nov 2014 A1
20140356884 Mittal Dec 2014 A1
20150038339 Gordon Feb 2015 A1
20150045234 Stone Feb 2015 A1
20150080231 Staker Mar 2015 A1
20150080232 Ju Mar 2015 A1
20150099642 Barany Apr 2015 A1
20150297949 Aman et al. Oct 2015 A1
Foreign Referenced Citations (20)
Number Date Country
1507217 Feb 2005 EP
1647327 Apr 2006 EP
2005164815 Jun 2005 JP
2006110523 Apr 2006 JP
2012152219 Aug 2012 JP
0027521 May 2000 WO
0058507 Oct 2000 WO
02093144 Nov 2002 WO
2004018493 Mar 2004 WO
2004018497 Mar 2004 WO
2004042403 May 2004 WO
2005042763 May 2005 WO
2006074351 Jul 2006 WO
2008037568 Apr 2008 WO
2008134867 Nov 2008 WO
2009085328 Jul 2009 WO
2011011738 Jan 2011 WO
2012047889 Apr 2012 WO
2014055963 Apr 2014 WO
2014165554 Oct 2014 WO
Non-Patent Literature Citations (17)
Entry
International Search Report and Written Opinion of the International Searching Authority for International Application No. PCT/IB2016/001318, dated Mar. 14, 2017, 18 pages.
Chen et al., “High-throughout platform for real-time monitoring of biological processes by multicolor single-molecule fluorescence”, Proceedings of the National Academy of Sciences, vol. III, No. 2, pp. 664-669.
Chen et al., Suppporting Information for “High-throughout platform for real-time monitoring of biological processes by multicolor single-molecule fluroescence”, Proceedings of the National Academy of Sciences, vol. III, No. 2, Dec. 30, 2013, 13 pages.
Eid et al., “Real-Time DNA Sequencing from Single Polymerase Molecules”, Science, vol. 323, No. 5910, Nov. 20, 2008, pp. 133-138.
Eid et al., Supporting Online Material for “REal-Time DNA Sequencing from Single Polymerase Molecules”, Science, vol. 323, No. 5910, Nov. 20, 2008, 21 pages.
Guo et al., “CMOS Time-Resolved, Contact and Multispectral Fluorescence Imaging for DNA Molecular Diagnostics”, Sensors, vol. 14, No. 11, Jan. 2, 2009, pp. 20602-20619.
Lerose et al., “CMOS-integrated geometrically tunable optical filters”, Applied Optics, vol. 52, No. 8, Mar. 10, 2013, pp. 1655-1662.
Sobhy et al., “Versatile single-molecule multi-color excitation and detection fluroescence setup for studying biomolecular dynamics”, Review of Scientific Instruments, vol. 82, Nov. 1, 2011, 11 pages.
Entire patent prosecution history of U.S. Appl. No. 13/802,957, filed Mar. 14, 2013, entitled, “Apparatus and Methods for Purging Flow Cells in Nucleic Acid Sequencing Instruments.”
Entire patent prosecution history of U.S. Appl. No. 13/832,509, filed Mar. 15, 2013, entitled, “Flow Cell Alignment Methods and Systems.”
Extended European Search Report for EP Application No. 14767413.9, dated Nov. 18, 2016, 8 pagaes.
International Search Report and Written Opinion for PCT International Application No. PCT/US2014/23474 dated Jul. 14, 2014.
International Search Report and Written Opinion of PCT International Application No. PCT/US2014/23490 dated Jun. 17, 2014.
Notice of Reasons for Rejection for Japanese Application No. 2016-501266, dated Feb. 20, 2018, including English translation, 18 pages.
Chinese Office Action for Chinese Application No. 201480028007.2, dated Apr. 2, 2018, including English translation, 10 pages.
Non Final Office Action for U.S. Appl. No. 15/253,015, dated Sep. 18, 2018, 14 pages.
Australian Examination Report for Australian Application No. 2014237067, dated Nov. 7, 2018, 4 pages.
Related Publications (1)
Number Date Country
20170132787 A1 May 2017 US
Continuations (1)
Number Date Country
Parent 13832509 Mar 2013 US
Child 15412743 US