1. Introduction
Long QT syndrome is a rare disorder of cardiac ion channels, characterised by a prolonged QT interval and T-wave abnormalities on electrocardiogram (ECG) and the occurrence of the ventricular tachycardia t
In New Zealand, the majority of specimens referred for Long QT syndrome diagnostic testing are retrieved after death. Postmortem specimens are often difficult to handle as they are usually either tissue samples or severely haemolysed blood. The extracted DNA is frequently of low quality, due to the presence of unwanted material such as short fragments produced by degradation and chemical modifications from oxidation and hydrolysis processes [7]. As a result, only short sequences can be reliably amplified [7]. Moreover, capillary-based sequencing, the gold standard first-line diagnostic test for Long QT syndrome, is very sensitive to the presence of contaminants, such as proteins, RNA and residual salt. The presence of such contaminants leads to poor quality electropherograms for analysis, which tend to be compromised by the presence of dye blobs, C-shoulders and a variable degree of baseline noise (Figure 1). Although the QC requirements for array comparative genomic hybridisation (aCGH), which can be used to detect whole exon deletion and duplication mutations, are less stringent than those for sequencing, contamination or degradation of sample DNA can lead to suboptimal efficiency of labelling and hybridisation. Such difficulties arising in the practical procedure mean that the analysis of postmortem DNA can be time-consuming and challenging, and obtaining high quality data within a reasonable timeframe can be extremely difficult.
Long QT syndrome affects 1 in 2000 individuals [8] and contributes to 15-25% of sudden unexplained death in 1-40 year olds [9,10] and 10% of sudden unexplained deaths in infancy (SIDs) [11]. As a consequence, it is an important differential diagnosis to be considered in all cases of autopsy negative young sudden death. Molecular genetic testing is essential to make a postmortem diagnosis, given that screening for the electrocardiogram characteristics is no longer possible after death. Historically the turnaround time for diagnosis has been up to six months, due to the large number of genes to be analysed and the difficulties inherent in dealing with postmortem specimens (as detailed above). However, the demand by coroners for diagnostic laboratories to undertake more rapid analysis has been steadily increasing. A protocol tailored specifically to the treatment of postmortem specimens is necessary to meet this demand. Within our laboratory, we have successfully developed a robust process for sequence and dosage analysis of postmortem samples and have achieved an excellent turnaround time of 6-12 weeks. Here, we address the critical QC parameters that should be considered in order to obtain high quality data for rapid, accurate analysis. The discussion presented below concerns both sequence analysis and dosage analysis. The former uses conventional sequencing technology, while the latter involves the simultaneous high resolution screening of deletion and duplication mutations in multiple cardiac genes as opposed to the more conventional multiplex ligation dependent probe amplification (MLPA) technique, which many diagnostic laboratories still use [6].
2. Materials and methods
2.1. Sequencing
2.1.1. Primer design
We have used two approaches to design primers flanking each of the coding exons of the genes
Following the above approach, we developed an alternate design protocol [12,13]. Each mRNA sequence of interest was identified through the public UCSC genome browser page
2.1.2. DNA extraction
Genomic DNA (gDNA) was extracted from peripheral blood leucocytes (EDTA blood samples) using the Gentra Puregene DNA Extraction kit (Qiagen), according to the manufacturer's instructions.
A standard phenol/chloroform protease protocol was used to extract gDNA from postmortem tissue specimens. A small section (2mm x 2mm x 4mm) is usually cut from frozen tissue and diced as finely as possible using a scalpel blade. The tissue is placed into a 1.5ml microcentrifuge tube with 450µl of 1x TES Buffer (1M NaCl, 0.5M Tris-HCl, 10mM EDTA), 60µl of 20µg/µl Proteinase K (Roche) and 10µl of 10% SDS, and incubated overnight with vigorous shaking. Following digestion, an equal volume of phenol is added and the sample is vortexed vigorously. Once homogenous, the sample is centrifuged to separate the layers and the top aqueous layer is removed and transferred to a fresh 1.5ml centrifuge tube. An equal volume of chloroform is then added to the aqueous layer. This is vortexed, centrifuged, and the aqueous (top) layer again transferred to a fresh tube. A 2x volume of 100% ethanol is added to the aqueous layer to precipitate the DNA, followed by centrifugation at 13,000rpm and the supernatant is removed. The pellet is then washed with 70% ethanol and re-centrifuged for 2 minutes at 13,000rpm. The supernatant is again removed, and the pellet air dried prior to re-suspension in TE buffer (10mM Tris-HCl, 1mM EDTA, pH 7.0-8.0).
The quality and quantity of extracted gDNA is measured using a NanoDrop ND-1000 Spectrophotometer.
2.1.3. PCR
PCR amplification is performed in a final 25µl reaction volume with the following reagents: Faststart buffer (Roche Applied Science), 2mM MgCl2 (Roche Applied Science), M13-tailed forward and reverse primers at 0.8µM each (synthesised by Integrated DNA Technologies Inc), 0.4mM dNTPs (GE Healthcare Ltd), 1 unit Faststart Taq DNA polymerase (Roche Applied Science) and 5µl GC-rich solution (Roche Applied Science). 50ng of gDNA is included in each reaction. PCR amplification is carried out with the following conditions: denaturation at 95°C for 5 minutes, followed by 35 cycles of 94°C for 45 seconds, 60°C for 30 seconds and 72°C for 30 seconds, with a final extension of 72°C for 10 minutes.
2.1.4. Sequencing (Figure 2)
5µL of each PCR is cleaned with ExoSAP-IT (Affymetrix, USB) prior to bidirectional DNA sequencing using M13 forward and reverse primers and Big-Dye Terminator v3.0 (Applied Biosystems Ltd). 20µl of each sequenced product is manually purified using the CleanSEQ Sequencing Purification System (Agencourt Bioscience). Four different drying times prior to elution (20 minutes, 24 hours, four days and seven days) were assessed to establish an optimal drying time for generating high quality sequencing data. 15µL of purified product was then subjected to capillary electrophoresis using the Applied Biosystems model 3130xl Genetic Analyzer.
The analysis of sequence traces is performed using Variant Reporter v1.0 (Applied Biosystems). Variant Reporter uses advanced algorithms and quality metrics to automate the detection of variants and to streamline the analysis process.
2.2. Array comparative genomic hybridization (aCGH)
A Roche NimbleGen 12x135K Custom CGH Array was used for dosage analysis. This bespoke CGH array has been designed to screen for dosage changes within the genes responsible for LQT1-12 (
Two hundred and fifty nanograms of gDNA are processed according to the manufacturer’s instructions (NimbleGen Array User’s Guide: CGH and CNV Arrays v8.0;
3. Results
In order to overcome the historical difficulties faced when performing Long QT syndrome testing using postmortem specimens and meet the demand by coroners for rapid results, we addressed the following parameters:
3.1. DNA purity
Since poor quality gDNA leads to suboptimal PCR amplification affecting downstream applications, the purity of gDNA is an important criterion for success in generating high quality sequencing data [14,15]. A NanoDrop ND-1000 spectrophotometer was used to measure the quality and quantity of the extracted gDNA. The ratio of absorbance at 260nm and 280nm (A260/280 ratio) is used to assess the purity of gDNA, which should be in the range of 1.8 to 2.0 to be accepted as pure gDNA. The ratio of absorbance at 260nm and 230nm (A260/230 ratio) is used as a secondary measure of nucleic acid purity, and for pure gDNA should be in the range of 2.0-2.2 [16]. Postmortem gDNA isolated using a standard phenol/chloroform protease protocol may contain residual phenol, chloroform or ethanol. These contaminants inhibit the activity of DNA polymerase in downstream applications (Figures 3-5) [14], so the purity of the gDNA must be strictly monitored. If a suboptimal A260/280 or A260/230 ratio indicates that the gDNA is of low quality, a secondary cleanup process should be considered. In our laboratory, we perform this cleanup using another phenol extraction and ethanol precipitation to further purify the gDNA sample [17].
3.2. Amount of DNA template used in a sequencing reaction
The extent of dilution of the cleaned PCR product prior to sequencing determines the amount of DNA template used in the sequencing reaction, which can affect the data quality for analysis [15]. The presence of excessive DNA template in a sequencing reaction will lead to rapidly progressive signal loss on the electropherogram (Figure 6), while using insufficient DNA template in a sequencing reaction will result in weak signal strength and a loss of peak shape (Figure 7), causing difficulties in basecalling and accurate analysis [15]. In order to obtain high quality data with good signal strength, the amplified product should be diluted in ddH20 to as close as possible to 2.5ng per 100bp.
3.3. CleanSEQ treatment of sequencing reactions
The sequenced product is purified using an Agencourt CleanSEQ system, which uses SPRI (Solid Phase Reversible Immobilization) magnetic bead-based technology. According to the manufacturer’s recommendation, the sequencing product should be cleaned using undiluted CleanSEQ reagent, but in our experience, postmortem DNA samples should be cleaned using CleanSEQ reagent diluted 1:2 in ddH20 (figures 8,9).
3.4. Drying of DNA-bound beads prior to elution
The length of drying time following ethanol washing and prior to elution of the sequencing product from the CleanSEQ beads is one of the critical parameters for obtaining high quality sequence data. We assessed 4 different drying periods: 20 minutes, 24 hours, four days and seven days.
A drying time of 20 minutes, as recommended by the manufacturer, is frequently associated with the presence of large ‘dye blobs’, most probably as a result of residual ethanol in the eluted product (Figure 10); this problem resolves if drying time is extended to 24 hours (Figure 11). However, variable baseline noise on the electropherogram will be evident if the drying time exceeds four days (Figure 12); further increased background noise and loss of resolution occur when the drying time exceeds seven days (Figure 13).
4. Discussion
4.1. DNA quality
4.1.1. Sequence analysis
The gDNA used in a PCR should be as pure as possible in order to optimise the quality of the template that will be used in downstream applications. The purity of DNA is assessed in our laboratory using a NanoDrop ND-1000 spectrophotometer. It is important to consider both A260/280 and A260/230 as poor DNA quality will affect downstream applications [16,18]. An unacceptable A260/280 or A260/230 ratio indicates the presence of contaminants in the DNA; an abnormal A260/230 indicates the presence of residual phenol or other chemical from the extraction process, while an abnormal A260/280 most frequently indicates the presence of protein [16,18]. If any of the ratios appear to be abnormal, DNA purification should be considered before processing any further.
4.1.2. Dosage analysis (aCGH)
The quality requirements for gDNA used in an aCGH assay are not as stringent as those for a sequencing assay. An awareness of the presence of significant degradation is important; however, partial compensation for this can be made by increasing the volume of gDNA used in the amplification and labelling step. In our laboratory, we have found that doubling the volume of gDNA when degradation is visible on the 'check gel' (present in lanes 5-8 of Figure 14) is an effective counter-measure. The adequacy of amplification is measured using a NanoDrop ND-1000 spectrophotometer in order to allow the hybridisation of equimolar amounts of test and control DNA to the array slide. A post-amplification concentration of greater than 2500ng/µl is ideal for further processing; a sample with a post-amplification concentration of less than 1500ng/µl is unlikely to produce good quality data for analysis. Therefore, if there is sufficient volume of gDNA available, a repeat amplification with an increased volume of template gDNA should be undertaken.
4.2. Amount of DNA template used in a sequencing reaction
The DNA template should be diluted to 2.5ng per 100bp prior to sequencing. Incorrect quantification and dilution will alter the amount of input DNA template in a sequencing reaction, which can lead to problematic data for analysis. In our laboratory, we use Variant ReporterTM (Applied Biosystems Ltd) for automated analysis of sequence data. The signal strength of the sequence data is automatically assessed. Data with good signal strength is above the value of 200 in the Variant ReporterTM software. Sequence data with weak signal strength, a value below 100, indicates insufficient DNA template has been used in the sequencing reaction. Noisy background and a loss of sharp peak shape are also commonly evident. A repeat treatment with EXOSAP-IT and less extensive dilution of product is necessary to avoid inaccurate analysis using sequence data with unacceptable quality. In contrast, when excess DNA template is used in a sequencing reaction, a rapidly progressive signal loss is seen on the electropherogram. In this instance, trace quality can be improved by diluting the eluted product with water and reloading the sample on the capillary sequencing platform (we use an Applied Biosystems model 3130xl Genetic Analyzer).
4.3. CleanSEQ treatment of sequencing reactions
Agencourt CleanSEQ is routinely used to purify the sequenced products. According to the manufacturer’s recommendation, 20µl of sequencing product should be cleaned with 10µl of CleanSEQ reagent. However, we found that this leads to an excess of CleanSEQ beads in relation to the amount of sequence product, therefore unwanted short fragments (e.g unincorporated dNTPs and excess primers) were able to bind to beads during purification, and the quality of sequence data generated was consequently reduced. We found that diluting the CleanSEQ beads 1:2 with ddH20 leads to much higher quality data.
4.4. Drying of DNA-bound beads prior to elution
This drying step is a critical QC consideration that significantly affects sequence quality. When the purified product is dried for 24 hours prior to elution, the subsequent electropherograms are free of dye blobs, C-shoulders and baseline noise, and consistently pass all QC metrics in Variant Reporter.
A drying time of 20 minutes, as suggested by manufacturer, was found to be problematic, causing the presence of dye blobs. Dye blobs are a common artefact due to either the presence of excess unincorporated dye or residual ethanol following product purification [15]. Previously, operators have focused on purifying sequenced products using an accurate final concentration of ethanol in an attempt to achieve consistency of evaporation/drying between sequencing runs. Ethanol is highly volatile and it is therefore difficult to achieve consistency - even the short period of time in which the vessel containing the ethanol is open to allow access can result in a decrease in concentration. Here, we showed that the issue of dye blobs can be resolved by simply extending the drying time to 24 hours, allowing adequate time for all ethanol to evaporate completely, despite any minor differences in the concentration of ethanol used.
On the other hand, we found that the quality of sequence data is adversely affected if the drying time is extended to more than 4 days. This suggests that the stability of purified sequencing products is another factor that affects data quality. High quality sequence data was obtained with a drying time of 24 hours; however, the sequences were suboptimal when elution and electrophoresis were carried out at day 4 or day 7 of drying. This indicates that purified sequencing product is most stable for the first 24 hours and that consideration of the sample stability should be taken into account when aiming for high quality data.
4.5. Analysis parameters
4.5.1. Sequencing
All sequence traces are analysed in our laboratory using Variant Reporter Software v1.0 (Applied Biosystems). The current CMGS best practice guidelines for Sanger sequence analysis in diagnostic laboratories [19] recommend a PHRED score of at least 20 for bidirectional data (corresponding to 99% confidence that the base is called correctly), and a PHRED score of at least 30 for unidirectional data (99.9% confidence that each base is called correctly). The procedure we describe above produces bidirectional sequence data that meet our laboratory's even more stringent analysis criteria: a Variant Reporter trace score of 35, corresponding to a false base call rate of 0.031%.
4.5.2. Dosage analysis
In order to provide simultaneous dosage analysis of all the Long QT syndrome genes, we developed a custom designed Roche NimbleGen 12x135K CGH array. As part of the validation of this array we analysed twenty patients with known copy number abnormalities [20]. Analysis of the data from these individuals with known copy number changes allowed precise threshold criteria to be developed. It was determined that a log2ratio ≤-0.4 over 6 contiguous probes is indicative of a deletion, and a log2ratio of ≥0.4 over 15 contiguous probes indicative of a duplication.
4.6. Other parameters
It is important to monitor other QC parameters when it comes to sequence analysis of Long QT specimens: the freshness of the Polymer (POP-7) and the usage of the capillaries in the Applied Biosystems model 3130xl Genetic Analyzer.
Polymer (POP-7) is used to separate DNA fragments on genetic analyzers and the polymer remains stable for up to 7 days [15]. The capillary in the genetic analyzer should be replaced after 1000 injections. Within our laboratory, we perform a regular weekly maintenance of the analyzer, and the capillary is replaced when 600 injections are reached. From past experience, both the freshness of polymer and the usage of the capillary play a vital role in the data quality. Electropherograms with reduced resolution and peak shape are produced when either of these two parameters is suboptimal. In order to achieve high quality data, both freshness of polymer and usage of the capillary should be closely monitored.
5. Conclusions
We have successfully established a robust method for processing postmortem specimens for Long QT diagnostic testing in a timely manner. The electropherograms in Figure 15 are indicative of the high quality data routinely produced, despite the limitations inherent in the types of specimens that are referred. The blood sample used to extract gDNA for sequencing in this example was heavily haemolysed on arrival at the laboratory.
The QC parameters described above should be monitored closely in order to consistently achieve optimal results (Figure 16). Although the basic procedure is essentially the same as that routinely used for sequence and aCGH analysis, the poor quality of postmortem specimens as a source of template DNA mean that particular attention needs to be paid to each step, in particular the critical initial assessment of DNA quality, the addition of the appropriate volume of template DNA to the sequencing reaction, and the drying time of the beads used in the purification procedure prior to elution.
Acknowledgments
We acknowledge the assistance of Dr Andrew Dodd in the initial design of primers against the coding regions of the LQT genes, and Mr Daniel Lai for his assistance in optimizing aspects of our amplicon clean-up.
References
- 1.
Hunter J. D. Sharma P. Rathi S. 2008 Long QT syndrome 8 67 70 - 2.
Modell S. M. Lehmann M. H. 2006 The long QT syndrome family of cardiac ion channelopathies: a HuGE review. 8 3 143 55 - 3.
Alders M. Mannens M. M. A. M. 2003 Romano-Ward Syndrome. Feb 20 [Updated 2012 May 31]. 1993Available from: http://www.ncbi.nlm.nih.gov/books/NBK1129/ - 4.
Yang Y. Yang Y. Liang B. Liu J. Li J. Grunnet M. Olesen S. P. Rasmussen H. B. Ellinor P. T. Gao L. Lin X. Li L. Wang L. Junjie J. X. Liu Y. Y. Liu Y. Y. Zhang S. S. Lian D. D. Peng L. Y. Jespersen T. Chen Y. H. 2010 Identification of a Kir3.4 mutation in congenital long QT syndrome. 86 872 880 - 5.
Dan M. Roden M. D. 2008 Long-QT Syndrome. 358 169 176 - 6.
Eddy C. A. MacCormick J. M. Chung S. K. Crawford J. R. Love D. R. Rees M. I. et al. 2008 Identification of large gene deletions and duplications in KCNQ1 and KCNH2 in patients with long QT syndrome. 5 9 1275 1281 - 7.
Bonin S. Petrera F. Niccolini B. Stanta G. 2003 PCR analysis in archival post mortem tissues. 56 184 186 - 8.
Schwartz P. J. Stramba-Badiale M. Crotti L. Pedrazzini M. Besana A. Bosi G. et al. 2009 Prevalence of the congenital long-QT syndrome. 120 18 1761 1767 - 9.
Skinner J. R. Crawford J. Smith W. Aitken A. Heaven D. Evans CA 2011 Prospective, population-based long QT molecular autopsy study of postmortem negative sudden death in 1 to 40 year olds. 8 3 412 419 - 10.
Gladding P. A. Evans C. A. Crawford J. Chung S. K. Vaughan A. Webster D. et al. 2010 Posthumous diagnosis of long QT syndrome from neonatal screening cards. 7 4 481 486 - 11.
Skinner J. R. 2010 Sudden Unexplained Death in Infancy and Long QT syndrome. 6 48 55 - 12.
Doherty E. Marquis-Nicholson R. Love J. M. Brookes C. Prosser D. Love D. R. 2011 Primer Design to Sequence Analysis- a Pipeline for a Molecular Genetic Diagnostic Laboratory, Applications and Experiences of Quality Control, Ognyan Ivanov (Ed.) 978-9-53307-236-4 InTech, Available from http://www.intechopen.com/books/applications-and-experiences-of-quality-control/primer-design-to-sequence-analysis-a-pipeline-for-a-molecular-genetic-diagnostic-laboratory - 13.
Lai D. Love D. R. 2012 Automation of a primer design and evaluation pipeline for subsequent sequencing of the coding regions of all human Refseq genes. 8 8 363 366 - 14.
Verhagen O. J. Wijkhuijs A. J. van der Sluijs-Gelling A. J. Szczepanski T. van der Linden-Schrever B. E. Pongers-Willemse van M. J. Wering E. R. van Dongen J. J. van der Schoot C. E. 1999 Suitable DNA isolation method for the detection of minimal residual disease by PCR techniques. 13 8 1298 1299 - 15.
Applied Biosystems Chemistry Guide, second edition 2009 DNA Sequencing by Capillary Electrophoresis. Available from: http://www3.appliedbiosystems.com/cms/groups/mcb_support/documents/generaldocuments/cms_041003.pdf - 16.
T009 - Technical Support Bulletin and 260/230Ratios. 2008 Available from: http://www.nanodrop.com/Library/T009-NanoDrop%201000-&-NanoDrop%208000-Nucleic-Acid-Purity-Ratios.pdf - 17.
Powell R. Gannon F. 2002 Purification of DNA by phenol extraction and ethanol precipitation. Available from: http://fds.oup.com/www.oup.co.uk/pdf/pas/9v1-7 -3.pdf - 18.
T042 - Technical Bulletin 2011 Assessment of Nucleic Acid Purity. Available from: http://www.nanodrop.com/Library/T042-NanoDrop-Spectrophotometers-Nucleic-Acid-Purity-Ratios.pdf - 19.
Ellard S. Charlton R. Yau M. Gokhale D. Taylor G. Wallace A. Ramsden S. C. 2009 Practice guidelines for Sanger Sequencing Analysis and Interpretation. Available from: http://www.cmgs.org/BPGs/pdfs%20current%20bpgs/Sequencingv2.pdf - 20.
Marquis-Nicholson R. Doherty E. Thrush A. Love J. M. Lan-C C. George A. M. Love D. R. 2012 Array-based identification of copy number changes: simultaneous gene-focused and low resolution whole genome analysis. submitted