1Department of Mechanical Engineering, Texas A&M University, 2Department of Mechanical Engineering and Department of Nuclear Engineering, Texas A&M University, 3Department of Chemical Engineering, Texas A&M University
Muddu, R., Hassan, Y. A., Ugaz, V. M. Rapid PCR Thermocycling using Microscale Thermal Convection. J. Vis. Exp. (49), e2366, doi:10.3791/2366 (2011).
Many molecular biology assays depend in some way on the polymerase chain reaction (PCR) to amplify an initially dilute target DNA sample to a detectable concentration level. But the design of conventional PCR thermocycling hardware, predominantly based on massive metal heating blocks whose temperature is regulated by thermoelectric heaters, severely limits the achievable reaction speed1. Considerable electrical power is also required to repeatedly heat and cool the reagent mixture, limiting the ability to deploy these instruments in a portable format.
Thermal convection has emerged as a promising alternative thermocycling approach that has the potential to overcome these limitations2-9. Convective flows are an everyday occurrence in a diverse array of settings ranging from the Earth's atmosphere, oceans, and interior, to decorative and colorful lava lamps. Fluid motion is initiated in the same way in each case: a buoyancy driven instability arises when a confined volume of fluid is subjected to a spatial temperature gradient. These same phenomena offer an attractive way to perform PCR thermocycling. By applying a static temperature gradient across an appropriately designed reactor geometry, a continuous circulatory flow can be established that will repeatedly transport PCR reagents through temperature zones associated with the denaturing, annealing, and extension stages of the reaction (Figure 1). Thermocycling can therefore be actuated in a pseudo-isothermal manner by simply holding two opposing surfaces at fixed temperatures, completely eliminating the need to repeatedly heat and cool the instrument.
One of the main challenges facing design of convective thermocyclers is the need to precisely control the spatial velocity and temperature distributions within the reactor to ensure that the reagents sequentially occupy the correct temperature zones for a sufficient period of time10,11. Here we describe results of our efforts to probe the full 3-D velocity and temperature distributions in microscale convective thermocyclers12. Unexpectedly, we have discovered a subset of complex flow trajectories that are highly favorable for PCR due to a synergistic combination of (1) continuous exchange among flow paths that provides an enhanced opportunity for reagents to sample the full range of optimal temperature profiles, and (2) increased time spent within the extension temperature zone the rate limiting step of PCR. Extremely rapid DNA amplification times (under 10 min) are achievable in reactors designed to generate these flows.
1. Basic Design of the Convective Thermocycler
2. Preparation and Loading of the PCR Reaction Mixture
3. Running the Reaction in the Convective Thermocycler
4. Performing Gel Electrophoresis Analysis
5. Representative Results:
Optimal design of convective thermocyclers involves choosing the correct reactor geometry that will generate a circulatory flow capable of transporting reagents through the key temperatures involved in the PCR process. The geometric parameters that can be varied in the cylindrical reactors considered here are the height (h) and diameter (d), or equivalently the aspect ratio (h/d). We have explored the 3-D flow fields inside convective PCR reactors over a range of different aspect ratios using computational fluid dynamics (CFD), and found that unexpectedly complex patterns can arise. More importantly, our analysis has uncovered a subset of these complex flow fields that significantly accelerate the reaction12.
These geometric effects can be clearly seen by comparing the flow fields in reactors at high and low aspect ratios. In the high aspect ratio case (h/d = 9; a tall, skinny cylinder), fluid elements are advected along trajectories tracing out paths that are essentially closed loops, and they are locked-in to follow the same paths for long periods of time (Figure 3a). Consequently, there is little opportunity for exchange between flow trajectories that expose reagents to the optimal sequence of thermal conditions for PCR (oscillating between the annealing and denaturing extremes at the reactor's top and bottom surfaces) and the much larger ensemble of remaining trajectories that do not contribute to amplification (localized closer to the center of the reactor).
The flow field is much different at a smaller aspect ratio (h/d = 3; a shorter, wider cylinder), becoming more chaotic in the sense that the fluid element trajectories no longer follow closed paths (Figure 3b). Thus, even though reagents are subjected to a more complex temperature profile, fluid elements are able to explore a much wider range of trajectories so that more of the reagent volume has an opportunity to experience optimal thermal profiles. These results counterintuitively suggest that while individual flow trajectories at h/d = 9 may appear to be more favorable for PCR by producing temperature profiles that "look" similar to those employed in a conventional thermocycler, the chaotic nature of the flow field at h/d = 3 ultimately dominates on a global scale by promoting enhanced exchange so that reagents do not become trapped in unfavorable trajectories for very long.
To test this hypothesis and determine which reactor geometry was more favorable for PCR, we used both of them to perform PCR replication of a 295 bp target associated with the β-actin gene from a human genomic DNA template. Remarkably, we routinely achieved amplification of the correct target product in only 10 min at h/d = 3 (Figure 4a), while the same reaction required at least 20 min before detectable products were observed at h/d = 9 (Figure 4b). The reaction specificity is also much greater at h/d = 3 where a single PCR product is obtained, while multiple nonspecific products are generated at h/d = 9 where the flow field traps fluid elements in unfavorable thermal trajectories for extended periods of time.
Figure 1. Thermal convection in a cylindrical chamber whose top and bottom surfaces are maintained at different fixed temperatures. If the temperature at the bottom surface is higher than at the top, a vertical density gradient is established within the enclosed fluid that is capable of generating a circulatory flow pattern. With the right choice of geometric parameters (height h and diameter d), the convective flow field can be harnessed to actuate PCR thermocycling when the top and bottom surfaces are manintained near annealing and denaturing temperatures, respectively (gravity acts vertically downward).
Figure 2. (a) Key components of a convective flow PCR thermocyler. A reactor cartridge consists of a polycarbonate block containing an array of cylindrical chambers. The block is sandwiched between a top plate designed to be connected to a circulating water bath and a bottom plate incorporating embedded heaters. (b) The PCR reagents are loaded and sealed inside the reactor cartridge, after which the device is assembled to perform the reaction. Upon completion, the procucts can be easily pipetted out for subsequent analysis.
Figure 3. Computational simulations of 3-D convective flow fields generated inside cylindrical reactors with temperatures of 53 and 96 °C imposed at the top and bottom surfaces, respectively. Results are shown at aspect ratios of (a) h/d = 9 (38.2 μL reactor volume) and (b) h/d = 3 (18.5 μL reactor volume). A representative trajectory followed by a fluid element traveling through the 3D flow field is shown at the left along with top- and side view projections of the path; the corresponding plot of temperature versus time following a fluid element is shown at the right. The thermal profile at h/d = 9 appears close to that of a conventional thermocycler, but the chaotic flow field at h/d = 3 is counterintuitively more favorable for PCR.
Figure 4. DNA replication results obtained in the reactor geometries shown in Figure 3 (upper and lower surfaces were maintained at 53 and 96 °C, respectively). (a) PCR is greatly accelerated at h/d = 3, as evident by strong products visible after only 10 min of reaction time (M: 100 bp ladder, lanes 1-4: products from parallel reactions in 4 different cylindrical chambers). (b) The same reaction performed at h/d = 9 requires at least 20 min before visible products are observed, and multiple bands are evident indicating replication of non-specific target DNA in addition to the desired product (M: 100 bp ladder, lanes 1-2: products from parallel reactions in 2 different cylindrical chambers).
The interplay between flow and chemical reaction can change dramatically under the influence of chaotic advection. But these complex flow states are challenging to study, particularly in realistic geometries where 3D effects become important. Rayleigh-Bénard convection at h/d > 1 not only provides a convenient platform to explore these phenomena, their application to perform PCR also enables the coupling between flow and chemical reactions to be probed. This unique capability makes it possible to select optimal flow states where chaotic effects act (counterintuitively) to accelerate the rate of DNA replication.
No conflicts of interest declared.
We gratefully thank T. Ragucci, B. Watkins, S. Wong, and B. Wallek at Lynntech, Inc. for technical assistance and many helpful discussions. This work was supported by the US National Science Foundation (NSF) under grant CBET-0933688.
|KOD DNA Polymerase kit||Novagen, EMD Millipore||71085-3|
|TaqMan β-actin Control Reagents kit||Applied Biosystems||401846|
|Aluminum tape||Axygen Scientific||PCR-AS-200|
|Bovine serum albumin||Sigma-Aldrich||A2153|
|Rain-X Anti-Fog||SOPUS Products|
|Long Gel-loading Pipette Tips||Fisher Scientific||05-408-151|
|FEP Teflon tape||McMaster-Carr||7562A17|
|TAE running buffer||Bio-Rad||141-0743|
|SYBR Green I||Invitrogen||S7563|
|6x Orange Loading Dye||Fermentas||R0631|
|100 bp DNA ladder sizing marker||Bio-Rad||170-8202|