3-GD of first derivative of high-resolution melt curves for estimation of mutant percentage in genome-edited samples.

<p>gDNA was isolated from HEK293T cells transfected with F8-S2 targeting RGENs or CCR5 targeting TALENs and PCR amplified using corresponding primer pairs to obtain high resolution melt curve data (Materials and Methods). 3-GD curve fitting was done on first derivative melt curves using CurveExpert Professional and <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0190192#pone.0190192.e013" target="_blank">Eq 12</a> as described in Materials and Methods. The individual Gaussians-g1(x) (purple dashed line), g2(x) (brown dashed line) and g3(x) (green dashed line) and their sum- g1(x)+ g2(x) + g3(x) (red solid line) were overlaid over the first derivative melt curve (blue dots). GD of F8-S2 is shown in A and of CCR5 in B. Table C shows the parameters (weights, centers and SDs) of 3-GD. The parameters that were fixed from GD of mocks and those that were set free during 3-GD of edited samples are shown in the Comments column. The g1 weight (w₁) represents the mutation frequencies in the amplicons of genome-edited F8-S2 and CCR5 target sites, respectively.</p

Validation of GD method using predefined amplicon mixes.

<p>High-resolution melt curves of samples containing different proportions of F8-S3 amplicon (S3Wt) or F8-S2 mutant (S2Mt) in F8-S2 amplicon (S2Wt) were analyzed by GD as detailed in Materials and Methods. (A) Derivative melt curve data (blue dots) of indicated S3Wt-S2Wt mixes were fitted using 3-GD (red traces). The nominal percentage of S3Wt in the mix is shown below (indicated by S3Wt%) and the GD-estimated amount in the top left corner of each plot. (B) Derivative melt curve data (blue dots) of indicated S2Mt-S2Wt mixes were fitted using 3-GD (red traces). The nominal percentage of S2Mt in the mix is shown below (indicated by S2Mt%) and the GD-estimated amount in the top left corner of each plot. (C) Scatter plot of nominal F8-S3Wt% in mix (X-axis) vs 3-GD estimated F8-S3Wt% (Y-axis). (D) Scatter plot of nominal F8-S2Mt% in mix (X-axis) vs 3-GD estimated F8-S2Mt% (Y-axis). The equations of linear regression analysis of both types of dose-response curves and the correlation coefficients (R²) are shown. The samples were tested in duplicate (replicates ‘a’ and ‘b’).</p

Parameters determined by 3-GD of two different size amplicons from the CCR5-edited target site.

<p>Parameters determined by 3-GD of two different size amplicons from the CCR5-edited target site.</p

Temperature-dependent quenching of fluorescence of free and dsDNA-bound fluorophore and its correction.

<p>(A) Plot of first-order polynomial curve fit of raw RFU vs. temperature in no template controls (NTC). The equation shown in the plot is the mean ± SD of six different sample slopes and constants. (B) The unprocessed high-resolution melting profile (blue trace) and the extrapolation from first-order polynomial curve fitting of the post-melt curve region (red dashed line) from an amplicon of an unedited target site. (C) High-resolution melting profile of background subtracted RFU (BcRFU, blue trace) and that of ‘unquenched’ or fluorescence-compensated BcRFU (FcRFU, green trace) from an unedited target site. The red dashed line shows extrapolation of pre-melt region from first-order polynomial curve fitting of BcRFU and depicts the predicted BcRFU in the absence of DNA melting. D) Comparison of first-order polynomial curve fitting of post-melt and pre-melt portions of melting curves. Normalized data were used to enable plotting of the two sets of data.</p

Size of PCR product does not affect determination of mutant percentage by GD.

<p>The CCR5 target site in gDNA of unmodified or genome-edited cells were amplified using two pairs of primers designed to produce two distinct sizes of product (107 bp and 140 bp, respectively). The amplicons were subjected to high-resolution melting and then processed to correct for temperature-dependent quenching of fluorescence of free and dsDNA-bound fluorophore. The resulting melt curves of genome-edited (for clone pair L1R1) and unmodified controls (Mock) are shown (A & C). Corresponding first-derivatives of processed melt curves are shown in B and D. Replicates G1 and G2, A1 and A2 refer to gDNA samples amplified using primers that produce 107 bp amplicon, whereas G5 and G6, and A5 and A6 refer to gDNA samples amplified using primers that produce 140 bp amplicon. The derivative melt curves were decomposed using the 3-GD model to estimate the mutant frequency. The estimated mutant frequencies for both sizes of amplicons are shown in (E). Error bar = 1 SD.</p

Analysis of F8-S2 and CCR5 target sequence features and melting properties <i>in silico</i>.

<p>Sliding window analysis of percentage of AT (%AT) in F8-S2 (A) or CCR5 (B) sequences of target sites amplified by PCR. The percentage of As and Ts were determined in a sliding overlapping window of 10-mers. The shift was by 1 bp. These are shown as green dashes. The data was smoothed using running averages with a period of 5 (solid green line). The sum of free energies (∆Gs) in a sliding window of 10-mers and a shift of 1 bp is shown along the left y-axis in kJ/mol (blue dots). The running averages were calculated as for %AT traces and are shown as blue traces. Putative AT-rich domains are marked I-IV. (C- H) The F8-S2 and CCR5 target sequences were used as input in the UMelt web analysis tool (29). UMelt predicted derivative melt curve (C and D), "Dynamic Profile” of melting (E and F) using a sliding temperature control that was situated close to the predicted Tm for each sequence to identify portions of the target sequences (nucleotide position indicated on the x-axis) that may have melted earlier than the rest. The web tool also provided a "Melting Profile" analysis that shows potential regions that might show greater tendency to melt earlier (G and H).</p

Comparison of mutant percentage estimation by 2- and 3-GD.

<p>First derivatives of high-resolution melt curves from genome-edited samples were curve fitted using 2- or 3-GD models as described in Materials and Methods (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0190192#pone.0190192.e012" target="_blank">Eq 11</a> and <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0190192#pone.0190192.e013" target="_blank">Eq 12</a>, respectively). The mutant percentages estimated from curve fitting are shown along the y-axis for F8-S2 (A) and CCR5 (B). Two molecular clones (10 and 11) of dgRNAs targeting F8-S2 site and two pairs of TALENs (L1R1 and L2R2) targeting CCR5 site were tested. The mutant percentages were compared using Student’s t-test (two-tailed). The p-values of the pair-wise comparisons of 2-GD and 3-GD are shown above the bars.</p

Comparison of different methods of processing melt curve data for background and fluorescence quenching correction.

<p>Melt curve data from amplicons of unmodified or control samples from F8-S2 (A) or CCR5 target loci (B) were either unprocessed (-dF/dT, blue trace) or corrected using exponential background subtraction method of Palais and Wittwer (24) (-dF/dT-dB/dT, red dashes) or the method described in this study (-d(nFcRFU)/dT, green trace).</p

2-GD model shows better fit than 1-GD for derivative melt curve data of mocks.

<p>2-GD model shows better fit than 1-GD for derivative melt curve data of mocks.</p

Flow diagram of steps involved in processing high-resolution melt curve data for GD.

<p>(A-D) The steps needed to process the raw melt curve data to correct for background from free fluorophore (A-B) and to correct for temperature dependent quenching of bound-fluorophore (C-D) by dividing background-corrected RFU (BcRFU) by the efficiency (C). The normalized fluorescence corrected RFU (nFCRFU) of mocks is then differentiated (E) before curve fitting using 2-GD (F) to determine the parameters values for use as constants in 3-GD of genome-edited samples. (G) Processed melt curves of genome-edited samples using the same steps outlined in A-E are then curve-fitted by 3-GD as described in the text. Equation numbers refer to those in Materials and Methods.</p