Identifying sequence variants in the human mitochondrial genome using high-resolution melt (HRM) profiling

DNA melting analysis

Melting is a fundamental property of DNA that can be monitored by absorbance or fluorescence. PCR conveniently produces enough DNA to be directly monitored on real-time instruments with fluorescently labeled probes or dyes. Dyes monitor the entire PCR product, while probes focus on a specific locus within the amplicon. Advances in amplicon melting include high resolution instruments, saturating DNA dyes that better reveal multiple products, prediction programs for domain melting, barcode taxonomic identification, high speed microfluidic melting, and highly parallel digital melting. Most single base variants and small insertions or deletions can be genotyped by high resolution amplicon melting. High resolution melting also enables heterozygote scanning for any variant within a PCR product. A web application (uMelt, http://www.dna-utah.org) predicts amplicon melting curves with multiple domains, a useful tool for verifying intended products. Additional applications include methylation assessment, copy number determination and verification of sequence identity. When amplicon melting does not provide sufficient detail, unlabeled probes or snapback primers can be used instead of covalently labeled probes. DNA melting is a simple, inexpensive, and powerful tool with many research applications that is beginning to make its mark in clinical diagnostics.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

Next-Generation Sequencing to Characterize Mitochondrial Genomic DNA Heteroplasmy

Mitochondria play a very important role in many crucial cellular functions. Each eukaryotic cell contains hundreds of mitochondria with hundreds of mitochondrial genomes. Mutant and wild-type mitochondrial DNA (mtDNA) may co-exist as heteroplasmy and cause human disease. The purpose of the protocols in this article is to simultaneously determine the mtDNA sequence and quantify the heteroplasmy level using parallel sequencing. The protocols include mitochondrial genomic DNA PCR amplification of two full-length products using two distinct sets of PCR primers. The PCR products are mixed at an equimolar ratio, and the samples are then barcoded and sequenced with high-throughput next-generation sequencing technology. This technology is highly sensitive, specific, and accurate in determining mtDNA mutations and the degree/level of heteroplasmy. © 2022 Wiley Periodicals LLC. Basic Protocol 1: PCR amplification of mitochondrial DNA Basic Protocol 2: Analysis of next-generation sequencing of mitochondrial DNA Basic Protocol 3: Mutect2 pipeline for automated sample processing and large-scale data analysis.

Investigation on the Prevalence of Canine Microfilaremia in Thailand Using a Novel Microfluidic Device in Combination with Real-Time PCR

We conducted a survey of canine microfilaraemia in 768 dogs in Chanthaburi, Samut Sakhon, and Narathiwat provinces of Thailand using a novel semi-automated, microfluidic device that is easy and rapid to perform. Microfilariae species were identified using High Resolution Melting real-time PCR (HRM real-time PCR). The prevalence of canine microfilaremia was 16.2% (45/278) in Chanthaburi and 5.5% (12/217) in Samut Sakhon. The prevalence of canine microfilaremia in Narathiwat was 22.7% (67/273). Brugia pahangi and Dirofilaria immitis were the predominant species of filariae found in the infected dogs from Chanthaburi and Narathiwat, respectively. The low prevalence of canine microfilaremia of Samut Sakhon may reflect the success of the Soi Dog foundation’s efforts and the establishment of veterinary control programs. An effective disease control and prevention strategies is needed in Chanthaburi and Narathiwat to reduce the risks of zoonotic transmission of the parasites. An appropriate drug treatment should be given to infected dogs and prophylactic drugs are suggested to be given to dogs age ≤1-year-old to prevent filarial infection. The novel microfluidic device could be implemented for surveillance of filariae infection in other animals.

Application of High Resolution Melt analysis (HRM) for screening haplotype variation in a non-model plant genus: Cyclopia (Honeybush)

Aim

This study has three broad aims: to (a) develop genus-specific primers for High Resolution Melt analysis (HRM) of members of Cyclopia Vent., (b) test the haplotype discrimination of HRM compared to Sanger sequencing, and (c) provide an example of using HRM to detect novel haplotype variation in wild C. subternata Vogel. populations.

Location

The Cape Floristic Region (CFR), located along the southern Cape of South Africa.

Methods

Polymorphic loci were detected through a screening process of sequencing 12 non-coding chloroplast DNA segments across 14 Cyclopia species. Twelve genus-specific primer combinations were designed around variable cpDNA loci, four of which failed to amplify under PCR; the eight remaining were applied to test the specificity, sensitivity and accuracy of HRM. The three top performing HRM Primer combinations were then applied to detect novel haplotypes in wild C. subternata populations, and phylogeographic patterns of C. subternata were explored.

Results

We present a framework for applying HRM to non-model systems. HRM accuracy varied across the PCR products screened using the genus-specific primers developed, ranging between 56 and 100%. The nucleotide variation failing to produce distinct melt curves is discussed. The top three performing regions, having 100% specificity (i.e. different haplotypes were never grouped into the same cluster, no false negatives), were able to detect novel haplotypes in wild C. subternata populations with high accuracy (96%). Sensitivity below 100% (i.e. a single haplotype being clustered into multiple unique groups during HRM curve analysis, false positives) was resolved through sequence confirmation of each cluster resulting in a final accuracy of 100%. Phylogeographic analyses revealed that wild C. subternata populations tend to exhibit phylogeographic structuring across mountain ranges (accounting for 73.8% of genetic variation base on an AMOVA), and genetic differentiation between populations increases with distance (p < 0.05 for IBD analyses).

Conclusions

After screening for regions with high HRM clustering specificity-akin to the screening process associated with most PCR based markers-the technology was found to be a high throughput tool for detecting genetic variation in non-model plants.

Genome sequencing and CRISPR/Cas9 gene editing of an early flowering Mini-Citrus (Fortunella hindsii)

Hongkong kumquat (Fortunella hindsii) is a wild citrus species characterized by dwarf plant height and early flowering. Here, we identified the monoembryonic F. hindsii (designated as 'Mini-Citrus') for the first time and constructed its selfing lines. This germplasm constitutes an ideal model for the genetic and functional genomics studies of citrus, which have been severely hindered by the long juvenility and inherent apomixes of citrus. F. hindsii showed a very short juvenile period (~8 months) and stable monoembryonic phenotype under cultivation. We report the first de novo assembled 373.6 Mb genome sequences (Contig-N50 2.2 Mb and Scaffold-N50 5.2 Mb) for F. hindsii. In total, 32 257 protein-coding genes were annotated, 96.9% of which had homologues in other eight Citrinae species. The phylogenomic analysis revealed a close relationship of F. hindsii with cultivated citrus varieties, especially with mandarin. Furthermore, the CRISPR/Cas9 system was demonstrated to be an efficient strategy to generate target mutagenesis on F. hindsii. The modifications of target genes in the CRISPR-modified F. hindsii were predominantly 1-bp insertions or small deletions. This genetic transformation system based on F. hindsii could shorten the whole process from explant to T1 mutant to about 15 months. Overall, due to its short juvenility, monoembryony, close genetic background to cultivated citrus and applicability of CRISPR, F. hindsii shows unprecedented potentials to be used as a model species for citrus research.

An Effective Strategy to Eliminate Inherent Cross-Contamination in mtDNA Next-Generation Sequencing of Multiple Samples

Heteroplasmic mutations in mitochondrial DNA (mtDNA) play critical roles in mitochondrial disease, aging, and cancer. Recently, next-generation sequencing (NGS) has been widely used to detect mtDNA mutations for diagnosis and monitoring of the above-mentioned diseases. However, little attention is paid on inherent cross-contamination generated during mtDNA capture and sequencing of mixed samples, which may seriously reduce the detection accuracy of mtDNA heteroplasmic mutations. In this study, a novel sequencing strategy based on a unique double-barcode design was established. The results showed that when single barcode-based analysis strategy was used, cross-contamination level of 20 DNA samples ranged from 0.27% to 11.90% on HiSeq 2500 and from 0.93% to 17.70% on HiSeq X ten, whereas double barcode-based strategy could effectively eliminate cross-contamination. Moreover, the data indicated that cross-contamination was mainly derived from capture process and was significantly affected by different NGS platforms. In addition, contamination level was negatively related to sequencing depth. Moreover, cross-contamination significantly increased the false-positive calling of mtDNA heteroplasmic mutations and remarkably affected the heteroplasmy level of mtDNA mutations. In contrast, cross-contamination had no notable effect on classification of mtDNA haplogroup. Taken together, our novel double barcode-based sequencing strategy is effective in eliminating cross-contamination, enhancing the detection accuracy of mtDNA NGS, and improving its application in diagnosis or monitoring of diseases associated with mtDNA mutations.

Diagnosis of feline filariasis assisted by a novel semi-automated microfluidic device in combination with high resolution melting real-time PCR

BACKGROUND

The diagnosis of filariasis traditionally relies on the detection of circulating microfilariae (mf) using Giemsa-stained thick blood smears. This approach has several limitations. We developed a semi-automated microfluidic device to improve and simplify the detection of filarial nematodes.

METHODS

The efficiency and repeatability of the microfluidic device was evaluated. Human EDTA blood samples were 'spiked' with B. malayi mf at high, moderate, and low levels, and subsequently tested 10 times. The device was also used for a field survey of feline filariasis in 383 domesticated cats in an area of Narathiwat Province, Thailand, the endemic area of Brugia malayi infection.

RESULTS

In the control blood arbitrarily spiked with mf, the high level, moderate level and low level mf-positive controls yielded coefficient variation (CV) values of 4.44, 4.16 and 4.66%, respectively, at the optimized flow rate of 6 µl/min. During the field survey of feline filariasis in Narathiwat Province, the device detected mf in the blood of 34 of 383 cats (8.9%) whereas mf were detected in 28 (7.3%) cats using the blood smear test. Genomic DNA was extracted from mf trapped in the device after which high-resolution melting (HRM) real-time PCR assay was carried out, which enabled the simultaneous diagnosis of filarial species. Among the 34 mf-positive samples, 12 were identified as B. malayi, 15 as Dirofilaria immitis and 7 as| D. repens.

CONCLUSIONS

We developed a semi-automated microfluidic device to detect mf of filarial parasites that could be used to diagnose lymphatic filariasis in human populations. This novel device facilitates rapid, higher-throughput detection and identification of infection with filariae in blood samples.

Mitochondrial and nuclear disease panel (Mito-aND-Panel): Combined sequencing of mitochondrial and nuclear DNA by a cost-effective and sensitive NGS-based method

Background

The diagnosis of mitochondrial disorders is challenging because of the clinical variability and genetic heterogeneity of these conditions. Next‐Generation Sequencing (NGS) technology offers a robust high‐throughput platform for nuclear and mitochondrial DNA (mtDNA) analyses.

Method

We developed a custom Agilent SureSelect Mitochondrial and Nuclear Disease Panel (Mito‐aND‐Panel) capture kit that allows parallel enrichment for subsequent NGS‐based sequence analysis of nuclear mitochondrial disease‐related genes and the complete mtDNA genome. Sequencing of enriched mtDNA simultaneously with nuclear genes was compared with the separated sequencing of the mitochondrial genome and whole exome sequencing (WES).

Results

The Mito‐aND‐Panel permits accurate detection of low‐level mtDNA heteroplasmy due to a very high sequencing depth compared to standard diagnostic procedures using Sanger sequencing/SNaPshot and WES which is crucial to identify maternally inherited mitochondrial disorders.

Conclusion

We established a NGS‐based method with combined sequencing of the complete mtDNA and nuclear genes which enables a more sensitive heteroplasmy detection of mtDNA mutations compared to traditional methods. Because the method promotes the analysis of mtDNA variants in large cohorts, it is cost‐effective and simple to setup, we anticipate this is a highly relevant method for sequence‐based genetic diagnosis in clinical diagnostic applications.

High-Resolution Melting Analysis for Rapid Detection of Sequence Type 131 Escherichia coli

Escherichia coli isolates belonging to the sequence type 131 (ST131) clonal complex have been associated with the global distribution of fluoroquinolone and β-lactam resistance. Whole-genome sequencing and multilocus sequence typing identify sequence type but are expensive when evaluating large numbers of samples. This study was designed to develop a cost-effective screening tool using high-resolution melting (HRM) analysis to differentiate ST131 from non-ST131 E. coli in large sample populations in the absence of sequence analysis. The method was optimized using DNA from 12 E. coli isolates. Singleplex PCR was performed using 10 ng of DNA, Type-it HRM buffer, and multilocus sequence typing primers and was followed by multiplex PCR. The amplicon sizes ranged from 630 to 737 bp. Melt temperature peaks were determined by performing HRM analysis at 0.1°C resolution from 50 to 95°C on a Rotor-Gene Q 5-plex HRM system. Derivative melt curves were compared between sequence types and analyzed by principal component analysis. A blinded study of 191 E. coli isolates of ST131 and unknown sequence types validated this methodology. This methodology returned 99.2% specificity (124 true negatives and 1 false positive) and 100% sensitivity (66 true positives and 0 false negatives). This HRM methodology distinguishes ST131 from non-ST131 E. coli without sequence analysis. The analysis can be accomplished in about 3 h in any laboratory with an HRM-capable instrument and principal component analysis software. Therefore, this assay is a fast and cost-effective alternative to sequencing-based ST131 identification.

The ability of human nuclear DNA to cause false positive low-abundance heteroplasmy calls varies across the mitochondrial genome

Background

Low-abundance mutations in mitochondrial populations (mutations with minor allele frequency ≤ 1%), are associated with cancer, aging, and neurodegenerative disorders. While recent progress in high-throughput sequencing technology has significantly improved the heteroplasmy identification process, the ability of this technology to detect low-abundance mutations can be affected by the presence of similar sequences originating from nuclear DNA (nDNA). To determine to what extent nDNA can cause false positive low-abundance heteroplasmy calls, we have identified mitochondrial locations of all subsequences that are common or similar (one mismatch allowed) between nDNA and mitochondrial DNA (mtDNA).

Results

Performed analysis revealed up to a 25-fold variation in the lengths of longest common and longest similar (one mismatch allowed) subsequences across the mitochondrial genome. The size of the longest subsequences shared between nDNA and mtDNA in several regions of the mitochondrial genome were found to be as low as 11 bases, which not only allows using these regions to design new, very specific PCR primers, but also supports the hypothesis of the non-random introduction of mtDNA into the human nuclear DNA.

Conclusion

Analysis of the mitochondrial locations of the subsequences shared between nDNA and mtDNA suggested that even very short (36 bases) single-end sequencing reads can be used to identify low-abundance variation in 20.4% of the mitochondrial genome. For longer (76 and 150 bases) reads, the proportion of the mitochondrial genome where nDNA presence will not interfere found to be 44.5 and 67.9%, when low-abundance mutations at 100% of locations can be identified using 417 bases long single reads. This observation suggests that the analysis of low-abundance variations in mitochondria population can be extended to a variety of large data collections such as NCBI Sequence Read Archive, European Nucleotide Archive, The Cancer Genome Atlas, and International Cancer Genome Consortium.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3375-x) contains supplementary material, which is available to authorized users.

Mitochondrial respiratory chain disorders in the Old Order Amish population

The Old Order Amish populations in the US are one of the Plain People groups and are descendants of the Swiss Anabaptist immigrants who came to North America in the early eighteenth century. They live in numerous small endogamous demes that have resulted in reduced genetic diversity along with a high prevalence of specific genetic disorders, many of them autosomal recessive. Mitochondrial respiratory chain deficiencies arising from mitochondrial or nuclear DNA mutations have not previously been reported in the Plain populations. Here we present four different Amish families with mitochondrial respiratory chain disorders. Mutations in two mitochondrial encoded genes leading to mitochondrial respiratory chain disorder were identified in two patients. In the first case, MELAS syndrome caused by a mitochondrial DNA (mtDNA) mutation (m.3243A>G) was identified in an extended Amish pedigree following a presentation of metabolic strokes in the proband. Characterization of the extended family of the proband by a high resolution melting assay identified the same mutation in many previously undiagnosed family members with a wide range of clinical symptoms. A MELAS/Leigh syndrome phenotype caused by a mtDNA mutation [m.13513G>A; p.Asp393Asn] in the ND5 gene encoding the ND5 subunit of respiratory chain complex I was identified in a patient in a second family. Mutations in two nuclear encoded genes leading to mitochondrial respiratory chain disorder were also identified in two patients. One patient presented with Leigh syndrome and had a homozygous deletion in the NDUFAF2 gene, while the second patient had a homozygous mutation in the POLG gene, [c.1399G>A; p.Ala467Thr]. Our findings identify mitochondrial respiratory chain deficiency as a cause of disease in the Old Order Amish that must be considered in the context of otherwise unexplained systemic disease, especially if neuromuscular symptoms are present.

High-resolution melt and morphological analyses of mealybugs (Hemiptera: Pseudococcidae) from cacao: tools for the control of Cacao swollen shoot virus spread

BACKGROUND

Mealybugs (Hemiptera: Coccoidea: Pseudococcidae) are key vectors of badnaviruses, including Cacao swollen shoot virus (CSSV), the most damaging virus affecting cacao (Theobroma cacao L.). The effectiveness of mealybugs as virus vectors is species dependent, and it is therefore vital that CSSV resistance breeding programmes in cacao incorporate accurate mealybug identification. In this work, the efficacy of a CO1-based DNA barcoding approach to species identification was evaluated by screening a range of mealybugs collected from cacao in seven countries.

RESULTS

Morphologically similar adult females were characterised by scanning electron microscopy, and then, following DNA extraction, were screened with CO1 barcoding markers. A high degree of CO1 sequence homology was observed for all 11 individual haplotypes, including those accessions from distinct geographical regions. This has allowed the design of a high-resolution melt (HRM) assay capable of rapid identification of the commonly encountered mealybug pests of cacao.

CONCLUSIONS

HRM analysis readily differentiated between mealybug pests of cacao that cannot necessarily be identified by conventional morphological analysis. This new approach, therefore, has potential to facilitate breeding for resistance to CSSV and other mealybug-transmitted diseases.

Rapid High-Resolution Melt Analysis of Cytauxzoon felis Cytochrome b To Aid in the Prognosis of Cytauxzoonosis

Cytauxzoon felis is a virulent, tick-transmitted, protozoan parasite that infects felines. Cytauxzoonosis was previously thought to be uniformly fatal in domestic cats. Treatment combining atovaquone and azithromycin (A&A) has been associated with survival rates of over 60%. Atovaquone, a ubiquinone analogue, targets C. felis cytochrome b (cytb), of which 30 unique genotypes have been identified. The C. felis cytb genotype cytb1 is associated with increased survival rates in cats treated with A&A. The purpose of this study was to design a PCR panel that could distinguish C. felis cytb1 from other cytochrome b genotypes. Primer pairs were designed to span five different nucleotide positions at which single-nucleotide polymorphisms in the C. felis cytb gene had been identified. Through the use of high-resolution melt analysis, this panel was predicted to distinguish cytb1 from other cytb genotypes. Assays were validated using samples from 69 cats with cytauxzoonosis for which the C. felis cytb genotypes had been characterized previously. The PCR panel identified C. felis cytb1 with 100% sensitivity and 98.2% specificity. High-resolution melt analysis can rapidly provide prognostic information for clients considering A&A treatment in cats with cytauxzoonosis.

Understanding structure, function, and mutations in the mitochondrial ATP synthase

The mitochondrial ATP synthase is a multimeric enzyme complex with an overall molecular weight of about 600,000 Da. The ATP synthase is a molecular motor composed of two separable parts: F₁ and F_o. The F₁ portion contains the catalytic sites for ATP synthesis and protrudes into the mitochondrial matrix. F_o forms a proton turbine that is embedded in the inner membrane and connected to the rotor of F₁. The flux of protons flowing down a potential gradient powers the rotation of the rotor driving the synthesis of ATP. Thus, the flow of protons though F_o is coupled to the synthesis of ATP. This review will discuss the structure/function relationship in the ATP synthase as determined by biochemical, crystallographic, and genetic studies. An emphasis will be placed on linking the structure/function relationship with understanding how disease causing mutations or putative single nucleotide polymorphisms (SNPs) in genes encoding the subunits of the ATP synthase, will affect the function of the enzyme and the health of the individual. The review will start by summarizing the current understanding of the subunit composition of the enzyme and the role of the subunits followed by a discussion on known mutations and their effect on the activity of the ATP synthase. The review will conclude with a summary of mutations in genes encoding subunits of the ATP synthase that are known to be responsible for human disease, and a brief discussion on SNPs.

Assessment of high resolution melting analysis as a potential SNP genotyping technique in forensic casework

High resolution melting (HRM) analysis is a simple, cost effective, closed tube SNP genotyping technique with high throughput potential. The effectiveness of HRM for forensic SNP genotyping was assessed with five commercially available HRM kits evaluated on the ViiA™ 7 Real Time PCR instrument. Four kits performed satisfactorily against forensically relevant criteria. One was further assessed to determine the sensitivity, reproducibility, and accuracy of HRM SNP genotyping. The manufacturer's protocol using 0.5 ng input DNA and 45 PCR cycles produced accurate and reproducible results for 17 of the 19 SNPs examined. Problematic SNPs had GC rich flanking regions which introduced additional melting domains into the melting curve (rs1800407) or included homozygotes that were difficult to distinguish reliably (rs16891982; a G to C SNP). A proof of concept multiplexing experiment revealed that multiplexing a small number of SNPs may be possible after further investigation. HRM enables genotyping of a number of SNPs in a large number of samples without extensive optimization. However, it requires more genomic DNA as template in comparison to SNaPshot®. Furthermore, suitably modifying pre-existing forensic intelligence SNP panels for HRM analysis may pose difficulties due to the properties of some SNPs.

Quantitative assessment of heteroplasmy of mitochondrial genome: perspectives in diagnostics and methodological pitfalls

The role of alterations of mitochondrial DNA (mtDNA) in the development of human pathologies is not understood well. Most of mitochondrial mutations are characterized by the phenomenon of heteroplasmy which is defined as the presence of a mixture of more than one type of an organellar genome within a cell or tissue. The level of heteroplasmy varies in wide range, and the expression of disease is dependent on the percent of alleles bearing mutations, thus allowing consumption that an upper threshold level may exist beyond which the mitochondrial function collapses. Recent findings have demonstrated that some mtDNA heteroplasmic mutations are associated with widely spread chronic diseases, including atherosclerosis and cancer. Actually, each etiological mtDNA mutation has its own heteroplasmy threshold that needs to be measured. Therefore, quantitative evaluation of a mutant allele of mitochondrial genome is an obvious methodological challenge, since it may be a keystone for diagnostics of individual genetic predisposition to the disease. This review provides a comprehensive comparison of methods applicable to the measurement of heteroplasmy level of mitochondrial mutations associated with the development of pathology, in particular, in atherosclerosis and its clinical manifestations.

Rapid and inexpensive body fluid identification by RNA profiling-based multiplex High Resolution Melt (HRM) analysis

Positive identification of the nature of biological material present on evidentiary items can be crucial for understanding the circumstances surrounding a crime. However, traditional protein-based methods do not permit the identification of all body fluids and tissues, and thus molecular based strategies for the conclusive identification of all forensically relevant biological fluids and tissues need to be developed. Messenger RNA (mRNA) profiling is an example of such a molecular-based approach. Current mRNA body fluid identification assays involve capillary electrophoresis (CE) or quantitative RT-PCR (qRT-PCR) platforms, each with its own limitations. Both platforms require the use of expensive fluorescently labeled primers or probes. CE-based assays require separate amplification and detection steps thus increasing the analysis time. For qRT-PCR assays, only 3-4 markers can be included in a single reaction since each requires a different fluorescent dye.

To simplify mRNA profiling assays, and reduce the time and cost of analysis, we have developed single- and multiplex body fluid High Resolution Melt (HRM) assays for the identification of common forensically relevant biological fluids and tissues. The incorporated biomarkers include IL19 (vaginal secretions), IL1F7 (skin), ALAS2 (blood), MMP10 (menstrual blood), HTN3 (saliva) and TGM4 (semen).  The HRM assays require only unlabeled PCR primers and a single saturating intercalating fluorescent dye (Eva Green). Each body-fluid-specific marker can easily be identified by the presence of a distinct melt peak. Usually, HRM assays are used to detect variants or isoforms for a single gene target. However, we have uniquely developed duplex and triplex HRM assays to permit the simultaneous detection of multiple targets per reaction. Here we describe the development and initial performance evaluation of the developed HRM assays. The results demonstrate the potential use of HRM assays for rapid, and relatively inexpensive, screening of biological evidence.

Identification and biochemical characterization of the novel mutation m.8839G>C in the mitochondrial ATP6 gene associated with NARP syndrome

Mutations in the ATP6 gene are reported to be associated with Leber hereditary optic neuropathy, bilateral striatal necrosis, coronary atherosclerosis risk and neuropathy, ataxia and retinitis pigmentosa (NARP)/maternally inherited Leigh syndromes. Here, we present a patient with NARP syndrome, in whom a previously undescribed mutation was detected in the ATP6 gene: m.8839G>C. Several observations support the concept that m.8839G>C is pathogenically involved in the clinical phenotype of this patient: (1) the mutation was heteroplasmic in muscle; (2) mutation load was higher in the symptomatic patient than in the asymptomatic carriers; (3) cybrids carrying this mutation presented lower cell proliferation, increased mitochondrial DNA (mtDNA) copy number, increased steady-state OxPhos protein levels and decreased mitochondrial membrane potential with respect to isogenic wild-type cybrids; (4) this change was not observed in 2959 human mtDNAs from different mitochondrial haplogroups; (5) the affected amino acid was conserved in all the ATP6 sequences analyzed; and (6) using in silico prediction, the mutation was classified as 'probably damaging'. However, measurement of ATP synthesis showed no differences between wild-type and mutated cybrids. Thus, we suggest that m.8839G>C may lower the efficiency between proton translocation within F0 and F1 rotation, required for ATP synthesis. Further experiments are needed to fully characterize the molecular mechanisms involved in m.8839G>C pathogenicity.

High-resolution melt analysis to detect sequence variations in highly homologous gene regions: application to CYP2B6

High-resolution melt (HRM) analysis using 'release-on-demand' dyes, such as EvaGreen(®) has the potential to resolve complex genotypes in situations where genotype interpretation is complicated by the presence of pseudogenes or allelic variants in close proximity to the locus of interest. We explored the utility of HRM to genotype a SNP (785A>G, K262R, rs2279343) that is located within exon 5 of the CYP2B6 gene, which contributes to the metabolism of a number of clinically used drugs. Testing of 785A>G is challenging, but crucial for accurate genotype determination. This SNP is part of multiple known CYP2B6 haplotypes and located in a region that is identical to CYP2B7, a nonfunctional pseudogene. Because small CYP2B6-specific PCR amplicons bracketing 785A>G cannot be generated, we simultaneously amplified both genes. A panel of 235 liver tissue DNAs and five Coriell samples were assessed. Eight CYP2B6/CYP2B7 diplotype combinations were found and a novel variant 769G>A (D257N) was discovered. The frequency of 785G corresponded to those reported for Caucasians and African-Americans. Assay performance was confirmed by CYP2B6 and/or CYP2B7 sequence analysis in a subset of samples, using a preamplified CYP2B6-specific long-range-PCR amplicon as HRM template. Inclusion rather than exclusion of a homologous pseudogene allowed us to devise a sensitive, reliable and affordable assay to test this CYP2B6 SNP. This assay design may be utilized to overcome the challenges and limitations of other methods. Owing to the flexibility of HRM, this assay design can easily be adapted to other gene loci of interest.

Primer effect in the detection of mitochondrial DNA point heteroplasmy by automated sequencing

The correct detection of mitochondrial DNA (mtDNA) heteroplasmy by automated sequencing presents methodological constraints. The main goals of this study are to investigate the effect of sense and distance of primers in heteroplasmy detection and to test if there are differences in the accurate determination of heteroplasmy involving transitions or transversions. A gradient of the heteroplasmy levels was generated for mtDNA positions 9477 (transition G/A) and 15,452 (transversion C/A). Amplification and subsequent sequencing with forward and reverse primers, situated at 550 and 150 bp from the heteroplasmic positions, were performed. Our data provide evidence that there is a significant difference between the use of forward and reverse primers. The forward primer is the primer that seems to give a better approximation to the real proportion of the variants. No significant differences were found concerning the distance at which the sequencing primers were placed neither between the analysis of transitions and transversions. The data collected in this study are a starting point that allows to glimpse the importance of the sequencing primers in the accurate detection of point heteroplasmy, providing additional insight into the overall automated sequencing strategy.

Comprehensive next-generation sequence analyses of the entire mitochondrial genome reveal new insights into the molecular diagnosis of mitochondrial DNA disorders

PURPOSE

The application of massively parallel sequencing technology to the analysis of the mitochondrial genome has demonstrated great improvement in the molecular diagnosis of mitochondrial DNA-related disorders. The objective of this study was to investigate the performance characteristics and to gain new insights into the analysis of the mitochondrial genome.

METHODS

The entire mitochondrial genome was analyzed as a single amplicon using a long-range PCR-based enrichment approach coupled with massively parallel sequencing. The interference of the nuclear mitochondrial DNA homologs was distinguished from the actual mitochondrial DNA sequences by comparison with the results obtained from conventional PCR-based Sanger sequencing using multiple pairs of primers.

RESULTS

Our results demonstrated the uniform coverage of the entire mitochondrial genome. Massively parallel sequencing of the single amplicon revealed the presence of single-nucleotide polymorphisms and nuclear homologs of mtDNA sequences that cause the erroneous and inaccurate variant calls when PCR/Sanger sequencing approach was used. This single amplicon massively parallel sequencing strategy provides an accurate quantification of mutation heteroplasmy as well as the detection and mapping of mitochondrial DNA deletions.

CONCLUSION

The ability to quantitatively and qualitatively evaluate every single base of the entire mitochondrial genome is indispensible to the accurate molecular diagnosis and genetic counseling of mitochondrial DNA-related disorders. This new approach may be considered as first-line testing for comprehensive analysis of the mitochondrial genome.Genet Med 2013:15(5):388-394.

Rapid detection and identification of Brugia malayi, B. pahangi, and Dirofilaria immitis by high-resolution melting assay

Human lymphatic filariasis is caused by filarial worms such as Brugia malayi for which the major reservoir is domestic cats. However, domestic cats or dogs also carry nonhuman filaria such as Brugia pahangi and Dirofilaria immitis. We have developed a single-tube, real-time PCR with a high-resolution melting (HRM) analysis assay for detection and identification of B. malayi, B. pahangi, and D. immitis in blood samples. The designated primer pair in the PCR can amplify a 114-bp region of mitochondrial 12S rRNA genes of these filarial worms. Subsequently, the HRM assay showed a specific melting temperature for each species. The assay showed the highest sensitivity and specificity in comparison with DNA sequences after assessment with 34 cat and 14 dog blood samples. This assay could be helpful for epidemiological studies of reservoirs and vectors.

Effect of glucocorticoid receptor gene polymorphisms on asthma phenotypes

The clinical presentation of asthma results from complex gene-gene and gene-environment interactions. The natural variability of the DNA sequence within the NR3C1 gene affects the activity of glucocorticoid receptors (GCRs). The NR3C1 gene is localized on chromosome 5q31–q32. The gene coding for the GCR comprises nine exons. The structural domains of the GCR determine the biological functions of the functional domains. The observed resistance to glucocorticosteroids and the normal metabolic profile of Tth111I single nucleotide polymorphism (SNP) carriers is due to the ER22/23EK polymorphism that is present in them. BclI polymorphism significantly affects the process of alternative NR3C1 gene splicing and within that mechanism increases the sensitivity to glucocorticoids (GCs). A total of 451 subjects were enrolled in the present study, including 235 qualified to the group of bronchial asthma patients. A group of 216 healthy participants with no history of asthma or atopic conditions was qualified for the study. Genotyping was accomplished using the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and PCR-high resolution melting (HRM) methods. No statistically significant differences were observed in the frequency of Tth111I, BclI and ER22/23EK polymorphisms of the NR3C1 gene when comparing mild, moderate and severe asthma vs. the control group. Investigative analyses demonstrated statistically significant correlations for alleles and genotypes of Tth111I polymorphism of the NR3C1 gene between healthy subjects and patients with severe asthma characterized by a control profile corresponding to an Asthma Control Test (ACT)™ score ≥20. It was established that only the Tth111I polymorphism of the NR3C1 gene plays an important role in the pathogenesis of chronic bronchitis leading to the development of asthma with both allergic and non-allergic etiology.

Applications of the method of high resolution melting analysis for diagnosis of Leber's disease and the three primary mutation spectrum of LHON in the Han Chinese population

Current screening methods, such as single strand conformational polymorphism (SSCP), denaturing high performance liquid chromatography (dHPLC) and direct DNA sequencing that are used for detecting mutation in Leber's hereditary optic neuropathy (LHON) subjects are time consuming and costly. Here we tested high-resolution melt (HRM) analysis for mtDNA primary mutations in LHON patients. In this study, we applied the high resolution melting (HRM) technology to screen mtDNA primary mutations in 50 LHON patients from their peripheral blood. In order to evaluate the reliability of this technique, we compared the results obtained by HRM and direct mtDNA sequencing. We also investigated the spectrum of three most common mtDNA mutations implicated in LHON in the Han Chinese population. The results showed HRM analysis differentiated all of the mtDNA primary mutations and identified 4 additional mtDNA mutations from 50 patients in the blind study. The prevalence of three primary mutations were 11778G>A (87.9%), 14484T>C (6.5%) and 3460G>A (1.7%) in the Han Chinese population. In conclusion, HRM analysis is a rapid, reliable, and low-cost tool for detecting mtDNA primary mutations and has practical applications in molecular genetics.

Synthetic spike-in standards improve run-specific systematic error analysis for DNA and RNA sequencing

While the importance of random sequencing errors decreases at higher DNA or RNA sequencing depths, systematic sequencing errors (SSEs) dominate at high sequencing depths and can be difficult to distinguish from biological variants. These SSEs can cause base quality scores to underestimate the probability of error at certain genomic positions, resulting in false positive variant calls, particularly in mixtures such as samples with RNA editing, tumors, circulating tumor cells, bacteria, mitochondrial heteroplasmy, or pooled DNA. Most algorithms proposed for correction of SSEs require a data set used to calculate association of SSEs with various features in the reads and sequence context. This data set is typically either from a part of the data set being “recalibrated” (Genome Analysis ToolKit, or GATK) or from a separate data set with special characteristics (SysCall). Here, we combine the advantages of these approaches by adding synthetic RNA spike-in standards to human RNA, and use GATK to recalibrate base quality scores with reads mapped to the spike-in standards. Compared to conventional GATK recalibration that uses reads mapped to the genome, spike-ins improve the accuracy of Illumina base quality scores by a mean of 5 Phred-scaled quality score units, and by as much as 13 units at CpG sites. In addition, since the spike-in data used for recalibration are independent of the genome being sequenced, our method allows run-specific recalibration even for the many species without a comprehensive and accurate SNP database. We also use GATK with the spike-in standards to demonstrate that the Illumina RNA sequencing runs overestimate quality scores for AC, CC, GC, GG, and TC dinucleotides, while SOLiD has less dinucleotide SSEs but more SSEs for certain cycles. We conclude that using these DNA and RNA spike-in standards with GATK improves base quality score recalibration.

Deep sequencing unearths nuclear mitochondrial sequences under Leber's hereditary optic neuropathy-associated false heteroplasmic mitochondrial DNA variants

Leber's hereditary optic neuropathy (LHON) is associated with mitochondrial DNA (mtDNA) ND mutations that are mostly homoplasmic. However, these mutations are not sufficient to explain the peculiar features of penetrance and the tissue-specific expression of the disease and are believed to be causative in association with unknown environmental or other genetic factors. Discerning between clear-cut pathogenetic variants, such as those that appear to be heteroplasmic, and less penetrant variants, such as the homoplasmic, remains a challenging issue that we have addressed here using next-generation sequencing approach. We set up a protocol to quantify MTND5 heteroplasmy levels in a family in which the proband manifests a LHON phenotype. Furthermore, to study this mtDNA haplotype, we applied the cybridization protocol. The results demonstrate that the mutations are mostly homoplasmic, whereas the suspected heteroplasmic feature of the observed mutations is due to the co-amplification of Nuclear mitochondrial Sequences.

Newborn screening for spinal muscular atrophy by calibrated short-amplicon melt profiling

BACKGROUND

The management options for the autosomal recessive neurodegenerative disorder spinal muscular atrophy (SMA) are evolving; however, their efficacy may require presymptom diagnosis and continuous treatment. To identify presymptomatic SMA patients, we created a DNA-based newborn screening assay to identify the homozygous deletions of the SMN1 (survival of motor neuron 1, telomeric) gene observed in 95%-98% of affected patients.

METHODS

We developed primers that amplify a 52-bp PCR product from homologous regions in the SMN1 and SMN2 (survival of motor neuron 2, centromeric) genes that flank a divergent site at site c.840. Post-PCR high-resolution melt profiling assessed the amplification product, and we used a unique means of melt calibration to normalize profiles. Samples that we had previously characterized for the numbers of SMN1 and SMN2 copies established genotypes associated with particular profiles. The system was evaluated with approximately 1000 purified DNA samples, 100 self-created dried blood spots, and >1200 dried blood spots from newborn screening tests.

RESULTS

Homozygous deletion of SMN1 exon 7 produced a distinctive melt profile that identified SMA patients. Samples with different numbers of SMN1 and SMN2 copies were resolved by their profiles. All samples with homozygous deletions were unambiguously recognized, and no normal sample was misidentified as a positive.

CONCLUSIONS

This assay has characteristics suitable for population-based screening. A reliable screening test will facilitate the identification of an SMA-affected cohort to receive early intervention to maximize the benefit from treatment. A prospective screening trial will allow the efficacy of treatment options to be assessed, which may justify the inclusion of SMA as a target for population screening.

Rapid genotyping of APOA5 -1131T>C polymorphism using high resolution melting analysis with unlabeled probes

The APOA5 -1131 T/C polymorphism (rs662799) exhibits a very strong association with elevated TG levels in different racial groups. High resolution melting (HRM) analysis with the use of unlabeled probes has shown to be a convenient and reliable tool to genotyping, but not yet been used for detecting rs662799 polymorphism. We applied the unlabeled probe HRM analysis and direct DNA sequencing to assay the -1131T>C SNP in 130 cases DNA samples blindly. This HRM analysis can be completed in <3 min for each sample. The two melting peaks were displayed at 66.1±0.4°C for CC homozygote and 68.7±0.2°C for TT homozygote; TC heterozygote showed the both melting peaks. The genotyping results by HRM method were completely concordant with direct DNA sequencing. The distribution of CC, TC, and TT genotypes for the -1131T>C SNP was 9.2, 49.2, and 41.5%, respectively. This assay was sensitive enough to detect C allele down to 20% and 10% for T allele. The limit of detection for C and T allele was 6.2 and 2.5 ng/μL DNA, respectively. The developed unlabeled probe HRM method provides an alternative mean to detect ApoA5 -1131T>C SNP rapidly and accurately.

Next generation sequencing to characterize mitochondrial genomic DNA heteroplasmy

This protocol describes the methodology to characterize mitochondrial DNA (mtDNA) heteroplasmy by parallel sequencing. Mitochondria play an important role in essential cellular functions. Each eukaryotic cell contains hundreds of mitochondria with hundreds of mitochondria genomes. Mutant and wild-type mtDNA may co-exist as heteroplasmy, and cause human disease. The purpose of this protocol is to simultaneously determine mtDNA sequence and quantify the heteroplasmic level. This protocol includes a two-fragment mitochondrial genome DNA PCR amplification. The PCR product is then mixed at an equimolar ratio. The samples are then barcoded and sequenced with high-throughput, next-generation sequencing technology. This technology is highly sensitive, specific, and accurate in determining mtDNA mutations and the level of heteroplasmy.

Large scale mtDNA sequencing reveals sequence and functional conservation as major determinants of homoplasmic mtDNA variant distribution

The mitochondrial DNA (mtDNA) is highly variable, containing large numbers of pathogenic mutations and neutral polymorphisms. The spectrum of homoplasmic mtDNA variation was characterized in 730 subjects and compared with known pathogenic sites. The frequency and distribution of variants in protein coding genes were inversely correlated with conservation at the amino acid level. Analysis of tRNA secondary structures indicated a preference of variants for the loops and some acceptor stem positions. This comprehensive overview of mtDNA variants distinguishes between regions and positions which are likely not critical, mainly conserved regions with pathogenic mutations and essential regions containing no mutations at all.

Rapid and cost-effective screening of newly identified microsatellite loci by high-resolution melting analysis

This study describes a new method for identifying microsatellite loci that will reliably amplify and show high degree of polymorphism in a given species. Microsatellites are the most powerful codominant markers available today, but the development of novel loci remains a labour-intensive and expensive process. In de novo isolation, approaches using next generation sequencing (NGS) are gradually replacing ones using Escherichia coli libraries, resulting in unparalleled numbers of candidate loci available. We present a systematic review of published microsatellite primer notes and show that, on average, about half of all candidate loci are lost due to insufficient PCR amplification, monomorphism or multicopy status in the genome, no matter what isolation strategy is used. Thus, the screening of candidate loci remains a major step in marker development. Re-assessing capillary-electrophoresis genotyped loci via high-resolution melting analysis (HRM), we evaluate the usefulness of HRM for this step. We demonstrate its applicability in a genotyping case study and introduce a fast, HRM-based workflow for the screening of microsatellite loci. This workflow may readily be applied to NGS-based marker development and has the potential to cut the costs of traditional testing by half to three quarters.

Searching for a needle in the haystack: comparing six methods to evaluate heteroplasmy in difficult sequence context

Mitochondrial DNA (mtDNA) mutations have been involved in disease, aging and cancer and furthermore exploited for evolutionary and forensic investigation. When investigating mtDNA mutations the peculiar aspects of mitochondrial genetics, such as heteroplasmy and threshold effect, require suitable approaches which must be sensitive enough to detect low-level heteroplasmy and, precise enough to quantify the exact mutational load. In order to establish the optimal approach for the evaluation of heteroplasmy, six methods were experimentally compared for their capacity to reveal and quantify mtDNA variants. Drawbacks and advantages of cloning, Fluorescent PCR (F-PCR), denaturing High Performance Liquid Chromatography (dHPLC), quantitative Real-Time PCR (qRTPCR), High Resolution Melting (HRM) and 454 pyrosequencing were determined. In particular, detection and quantification of a mutation in a difficult sequence context were investigated, through analysis of an insertion in a homopolymeric stretch (m.3571insC).

Quantitative PCR coupled with melt curve analysis for detection of selected pseudo-nitzschia spp. (Bacillariophyceae) from the Northwestern Mediterranean sea

The frequency and intensity of Pseudo-nitzschia spp. blooms along the coast of Catalonia have been increasing over the past 20 years. As species from this genus that are documented as toxigenic have been found in local waters, with both toxic and nontoxic species cooccurring in the same bloom, there is a need to develop management tools for discriminating the difference. Currently, differentiation of toxic and nontoxic species requires time-consuming electron microscopy to distinguish taxonomic features that would allow identification as to species, and cryptic species can still remain misidentified. In this study, cells of Pseudo-nitzschia from clonal cultures isolated from seawater were characterized to their species identity using scanning electron microscopy, and subsamples of each culture were used to create an internal transcribed spacer 1 (ITS-1), 5.8S, and ITS-2 ribosomal DNA database for development of species-specific quantitative PCR (qPCR) assays. Once developed, these qPCR assays were applied to field samples collected over a 2-year period in Alfaques Bay in the northwestern Mediterranean Sea to evaluate the possibility of a comprehensive surveillance for all Pseudo-nitzschia spp. using molecular methods to supplement optical microscopy, which can discern taxonomy only to the genus level within this taxon. Total Pseudo-nitzschia cell density was determined by optical microscopy from water samples collected weekly and compared to results obtained from the sum of eight Pseudo-nitzschia species-specific qPCR assays using duplicate samples. Species-specific qPCR followed by melt curve analysis allowed differentiation of amplicons and identification of false positives, and results correlated well with the total Pseudo-nitzschia cell counts from optical microscopy.

Detection of hemi/homozygotes through heteroduplex formation in high-resolution melting analysis

Heteroduplex formation, required for the complete detection of hemi/homozygotes using high-resolution melting analysis, can be induced either by pre-PCR mixing of genomic DNAs or by post-PCR mixing of PCR products from unknown and reference samples. This study investigates the effects of both methods using two single nucleotide polymorphisms in X-linked DMD gene. The results show that both methods resulted in the same effect when mixing samples with the same gene copy number. Mixing samples with different gene copy numbers has not been previously explored and we show that post-PCR mixing is insensitive to gene copy number differences as compared to pre-PCR mixing.

Paternal mosaicism of an STXBP1 mutation in OS

Ohtahara syndrome (OS) is one of the most severe and earliest forms of epilepsy. We have recently identified that the de novo mutations of STXBP1 are important causes for OS. Here we report a paternal somatic mosaicism of an STXBP1 mutation. The affected daughter had onset of spasms at 1 month of age, and interictal electroencephalogram showed suppression-burst pattern, leading to the diagnosis of OS. She had a heterozygous c.902+5G>A mutation of STXBP1, which affects donor splicing of exon 10, resulting in 138-bp insertion of intron 10 sequences in the transcript. The mutant transcript had a premature stop codon, and was degraded by nonsense-mediated mRNA decay in lymphoblastoid cells derived from the patient. High-resolution melting analysis of clinically unaffected parental DNAs suggested that the father was somatic mosaic for the mutation, which was also suggested by sequencing. Cloning of PCR products amplified with the paternal DNA samples extracted from blood, saliva, buccal cells, and nails suggested that 5.3%, 8.7%, 11.9%, and 16.9% of alleles harbored the mutation, respectively. This is a first report of somatic mosaicism of an STXBP1 mutation, which has implications in genetic counseling of OS.

COLD PCR HRM: a highly sensitive detection method for IDH1 mutations

The p.Arg132His mutation of isocitrate dehydrogenase 1 (IDH1(R132H) ) is a frequent alteration and a major prognostic marker in gliomas. However, direct sequencing of highly contaminated tumor samples may fail to detect this mutation. Our objective was to evaluated the sensitivity of a newly described amplification method, coamplification at lower temperature-PCR (COLD PCR), combined with high-resolution melting (HRM) for the detection of the IDH1(R132H) mutation. To this end, we used serial dilutions of mutant DNA with wild-type DNA. PCR-HRM assay detects IDH1(R132H) at an abundance of 25%, similar to the detection limit of direct Sanger sequencing. Introducing a run of COLD PCR allows the detection of 2% mutant DNA. Using two consecutive runs of COLD PCR, we detected 0.25% mutant DNA in a background of wild-type DNA, that mimics a tumor sample highly contaminated by normal DNA. We then analyzed 10 biopsies of tumor edges, considered free of tumor cells by histological analysis, and showed that immunohistochemistry of IDH1(R132H) was positive in three cases (30%), whereas double COLD PCR HRM was positive in the 10 cases studied (100%). In summary, COLD PCR HRM analysis is 100-fold more sensitive than Sanger sequencing, rendering this rapid and powerful strategy particularly useful for samples highly contaminated with normal tissue.

Mutational analysis of whole mitochondrial DNA in patients with MELAS and MERRF diseases

Mitochondrial diseases are clinically and genetically heterogeneous disorders, which make the exact diagnosis and classification difficult. The purpose of this study was to identify pathogenic mtDNA mutations in 61 Korean unrelated families (or isolated patients) with MELAS or MERRF. In particular, the mtDNA sequences were completely determined for 49 patients. From the mutational analysis of mtDNA obtained from blood, 5 confirmed pathogenic mutations were identified in 17 families, and 4 unreported pathogenically suspected mutations were identified in 4 families. The m.3243A>G in the tRNA(Leu(UUR))was predominantly observed in 10 MELAS families, and followed by m.8344A>G in the tRNA(Lys) of 4 MERRF families. Most pathogenic mutations showed heteroplasmy, and the rates were considerably different within the familial members. Patients with a higher rate of mutations showed a tendency of having more severe clinical phenotypes, but not in all cases. This study will be helpful for the molecular diagnosis of mitochondrial diseases, as well as establishment of mtDNA database in Koreans.

Detecting heteroplasmy from high-throughput sequencing of complete human mitochondrial DNA genomes

Heteroplasmy, the existence of multiple mtDNA types within an individual, has been previously detected by using mostly indirect methods and focusing largely on just the hypervariable segments of the control region. Next-generation sequencing technologies should enable studies of heteroplasmy across the entire mtDNA genome at much higher resolution, because many independent reads are generated for each position. However, the higher error rate associated with these technologies must be taken into consideration to avoid false detection of heteroplasmy. We used simulations and phiX174 sequence data to design criteria for accurate detection of heteroplasmy with the Illumina Genome Analyzer platform, and we used artificial mixtures and replicate data to test and refine the criteria. We then applied these criteria to mtDNA sequence reads for 131 individuals from five Eurasian populations that had been generated via a parallel tagged approach. We identified 37 heteroplasmies at 10% frequency or higher at 34 sites in 32 individuals. The mutational spectrum does not differ between heteroplasmic mutations and polymorphisms in the same individuals, but the relative mutation rate at heteroplasmic mutations is significantly higher than that estimated for all mutable sites in the human mtDNA genome. Moreover, there is also a significant excess of nonsynonymous mutations observed among heteroplasmies, compared to polymorphism data from the same individuals. Both mutation-drift and negative selection influence the fate of heteroplasmies to determine the polymorphism spectrum in humans. With appropriate criteria for avoiding false positives due to sequencing errors, next-generation technologies can provide novel insights into genome-wide aspects of mtDNA heteroplasmy.

High-resolution DNA melting analysis in clinical research and diagnostics

Among nucleic acid analytical methods, high-resolution melting analysis is gaining more and more attention. High-resolution melting provides simple, homogeneous solutions for variant scanning and genotyping, addressing the needs of today's overburdened laboratories with rapid turnaround times and minimal cost. The flexibility of the technique has allowed it to be adopted by a wide range of disciplines for a variety of applications. In this review we examine the broad use of high-resolution melting analysis, including gene scanning, genotyping (including small amplicon, unlabeled probe and snapback primers), sequence matching and methylation analysis. Four major application arenas are examined to demonstrate the methods and approaches commonly used in particular fields. The appropriate usage of high-resolution melting analysis is discussed in the context of known constraints, such as sample quality and quantity, with a particular focus placed on proper experimental design in order to produce successful results.

Genotyping mitochondrial DNA single nucleotide polymorphisms by PCR ligase detection reactions

BACKGROUND

The identification of human mitochondrial DNA (mtDNA) sequence variations, especially single nucleotide polymorphisms (SNPs), is important for many applications. The PCR-ligase detection reaction (LDR) method can reduce false-positives and eliminate the need for both post-PCR and post-ligation purifications in SNP analyses. In addition, it has been successfully employed to detect point mutations in various nuclear genes. In this study, we used the PCR-LDR platform to characterize mtDNA SNPs.

METHODS

Multiplex PCR-LDRs were used to genotype 19 mtDNA single nucleotide polymorphic sites from 812 samples. Performance of the method was assessed by direct sequencing of 44 samples.

RESULTS

We established an overall 97.4% success rate with 99.2% accuracy using the multiplex PCR-LDR methodology.

CONCLUSIONS

The PCR-LDR mtDNA genotyping technique is simple, highly accurate, has high-throughput, and is cost-effective. Therefore, this method is applicable to mtDNA haplotyping in various applications.

High-resolution melting analysis, a simple and effective method for reliable mutation scanning and frequency studies in the ACADVL gene

Expanded newborn screening uses tandem mass spectrometry (MS/MS) to identify patients affected with fatty acid oxidation defects by the presence of pathological acylcarnitine esters. A caveat to MS/MS assessment is that cut-off values for disease-specific acylcarnitines does not always clearly discriminate affected patients from carriers and healthy individuals. Diagnostic evaluation of screening-positive samples is required to confirm a metabolic deficiency. With MS/MS newborn screening becoming established in a growing number of countries, streamlined means for time- and -effective follow-on diagnostic evaluation is essential. Moreover, studies to evaluate the diagnostic accuracy of MS/MS newborn screening are needed for determination and adjustment of precise cut-off values for the disease-specific acylcarnitines. In the current study, we use the fatty acid oxidation disorder very-long-chain acyl-CoA dehydrogenase deficiency (VLCADD), the second most common fatty acid oxidation disorder detected by expanded newborn screening, to demonstrate accurate and fast diagnostic evaluation of the ACADVL gene utilizing DNA extracted from the newborn screening dried blood spot and high resolution melt (HRM) profiling. We also demonstrate that HRM is a very effective means to determine carrier frequency of prevalent ACADVL mutations in the general population. Based on estimates of the expected disease incidence, we discuss the diagnostic accuracy of MS/MS-based newborn screening to identify VLCADD in Denmark.