Differentiating between monozygotic twins through DNA methylation-specific high-resolution melt curve analysis

Although short tandem repeat profiling is extremely powerful in identifying individuals from crime scene stains, it is unable to differentiate between monozygotic (MZ) twins. Efforts to address this include mutation analysis through whole genome sequencing and through DNA methylation studies. Methylation of DNA is affected by environmental factors; thus, as MZ twins age, their DNA methylation patterns change. This can be characterized by bisulfite treatment followed by pyrosequencing. However, this can be time-consuming and expensive; thus, it is unlikely to be widely used by investigators. If the sequences are different, then in theory the melting temperature should be different. Thus, the aim of this study was to assess whether high-resolution melt curve analysis can be used to differentiate between MZ twins. Five sets of MZ twins provided buccal swabs that underwent extraction, quantification, bisulfite treatment, polymerase chain reaction amplification and high-resolution melting curve analysis targeting two markers, Alu-E2F3 and Alu-SP. Significant differences were observed between all MZ twins targeting Alu-E2F3 and in four of five MZ twins targeting Alu-SP (P < 0.05). Thus, it has been demonstrated that bisulfite treatment followed by high-resolution melting curve analysis could be used to differentiate between MZ twins.     

A reappraisal of complete mtDNA variation in East Asian families with hearing impairment

In a number of recent studies, we summarized the obvious errors and shortcomings that can be spotted in many (if not most) mitochondrial DNA (mtDNA) data sets published in medical genetics. We have reanalyzed here the complete mtDNA genome data published in various recent reports of East Asian families with hearing impairment, using a phylogenetic approach, in order to demonstrate the persistence of lab-specific mistakes in mtDNA genome sequencing in cases where those caveats were (deliberately) neglected. A phylogenetic reappraisal of complete mtDNAs with mutation A1555G (or G11778A) indeed supports the suggested lack of association between haplogroup background and phenotypic presentation of these mutations in East Asians. In contrast, the claimed pathogenicity of mutation T1095C in Chinese families with hearing impairment seems unsupported, basically because this mutation is rather basal in the mtDNA phylogeny, being specific to haplogroup M11 in East Asia. The roles of other haplogroup specific or associated variants, such as A827G, T961C, T1005C, in East Asian subjects with aminoglycoside-induced and non-syndromic hearing loss are also unclear in view of the known mtDNA phylogeny.     

Capturing single-copy nuclear genes,organellar genomes,and nuclear ribosomal DNA from deep genome skimming data for plant phylogenetics: A case study in Vitaceae

With the decreasing cost and availability of many newly developed bioinformatics pipelines, next-generation sequencing (NGS) has revolutionized plant systematics in recent years. Genome skimming has been widely used to obtain high-copy fractions of the genomes, including plastomes, mitochondrial DNA (mtDNA), and nuclear ribosomal DNA (nrDNA). In this study, through simulations, we evaluated the optimal (minimum) sequencing depth and performance for recovering single-copy nuclear genes (SCNs) from genome skimming data, by subsampling genome resequencing data and generating 10 data sets with different sequencing coverage in silico. We tested the performance of four data sets (plastome, nrDNA, mtDNA, and SCNs) obtained from genome skimming based on phylogenetic analyses of the Vitis clade at the genus level and Vitaceae at the family level, respectively. Our results showed that optimal minimum sequencing depth for high-quality SCNs assembly via genome skimming was about 10× coverage. Without the steps of synthesizing baits and enrichment experiments, coupled with incredibly low sequencing costs, we showcase that deep genome skimming (DGS) is as effective for capturing large data sets of SCNs as the widely used Hyb-Seq approach, in addition to capturing plastomes, mtDNA, and entire nrDNA repeats. DGS may serve as an efficient and economical alternative and may be superior to the popular target enrichment/Hyb-Seq approach.     

Unique mitochondrial DNA in highly inbred feral cattle

The Chillingham herd of wild Northumbrian cattle remains viable despite over 300 years of in-breeding and a near-homozygous nuclear genome. Here we report the complete mitochondrial DNA sequence using ultra-deep next generation sequencing. Random population sampling of ~10% of the extant herd identified a single mtDNA haplotype harbouring a unique bovine variant present in all other higher mammals (m.11789C/Y421H) which may contribute to their survival.     

Very Low-Level Heteroplasmy mtDNA Variations Are Inherited in Humans

Little is known about the inheritance of very low heteroplasmy mitochondria DNA （mtDNA） variations. Even with the development of new next-generation sequencing methods, the practical lower limit of measured heteroplasmy is still about 1% due to the inherent noise level of the sequencing. In this study, we sequenced the mitochondrial genome of 44 individuals using Illumina high-throughput sequencing technology and obtained high-coverage mitochondria sequencing data. Our study population contains many mother-offspring pairs. This unique study design allows us to bypass the usual heteroplasmy limitation by analyzing the correlation of mutation levels at each position in the mtDNA sequence between maternally related pairs and non-related pairs. The study showed that very low heteroplasmy variants, down to almost 0.1%, are inherited maternally and that this inheritance begins to decrease at about 0.5%, cor- resnondin to abottleneck of about 200 mtDNA.     

Analysis of mitochondrial DNA in microfluidic systems

Abnormalities in mitochondrial function play a major role in many human diseases. It is often of critical importance to ascertain what proportion of the mitochondria within a cell, or cells, bear a given mutation (the mitochondrial "demographics"). In this work, a rapid, novel, on-chip procedure was used, in which a restriction enzyme was employed to excise a mitochondrial DNA (mtDNA) sequence from plasmid DNA that acted as a prototypical mitochondrial genome. The DNA was then denatured, reassembled to form duplexes, fluorescently labelled and analysed. This method was able to differentiate between a homogeneous population and a heterogeneous population. Using a microfluidic chip, the method could be performed in about 45 min, even without robotics or multiplexed operation, whereas conventional methods of analysis require days to perform. This method may ultimately form the basis for a means of characterizing the mitochondrial demographics of a single cell.     

Regionally Specific and Genome-Wide Analyses Conclusively Demonstrate the Absence of CpG Methylation in Human Mitochondrial DNA

Although CpG methylation clearly distributes genome-wide in vertebrate nuclear DNA, the state of methylation in the vertebrate mitochondrial genome has been unclear. Several recent reports using immunoprecipitation, mass spectrometry, and enzyme-linked immunosorbent assay methods concluded that human mitochondrial DNA (mtDNA) has much more than the 2 to 5% CpG methylation previously estimated. However, these methods do not provide information as to the sites or frequency of methylation at each CpG site. Here, we have used the more definitive bisulfite genomic sequencing method to examine CpG methylation in HCT116 human cells and primary human cells to independently answer these two questions. We found no evidence of CpG methylation at a biologically significant level in these regions of the human mitochondrial genome. Furthermore, unbiased next-generation sequencing of sodium bisulfite treated total DNA from HCT116 cells and analysis of genome-wide sodium bisulfite sequencing data sets from several other DNA sources confirmed this absence of CpG methylation in mtDNA. Based on our findings using regionally specific and genome-wide approaches with multiple human cell sources, we can definitively conclude that CpG methylation is absent in mtDNA. It is highly unlikely that CpG methylation plays any role in direct control of mitochondrial function.     

The complete nucleotide sequence of Chinese hamster (Cricetulus griseus) mitochondrial DNA.

Mammalian mitochondria contain their own approximately 16.5 kb circular genome. Mitochondrial DNA (mtDNA) encodes for a subset of the proteins involved in the electron transport chain and depletion or mutation of the sequence is implicated in a number of human disease processes. The recent finding is that mitochondrial damage mediates genotoxicity after exposure to chemical carcinogens has focused attention on the role of mtDNA mutations in the development of cancer. Although the entire genome has been sequenced for a number of mammals, only a small fraction of the mtDNA sequence is available for hamsters. We have obtained here the entire 16,284 bp sequence of the Chinese hamster mitochondrial genome, which will enable detailed analysis of mtDNA mutations caused by exposure to mutagens in hamster-derived cell lines.     

Analysis of the complete human mtDNA genome: methodology and inferences for human evolution

The analysis of mitochondrial DNA (mtDNA) sequences has been a potent tool in our understanding of human evolution. However, almost all studies of human evolution based on mtDNA sequencing have focused on the control region, which constitutes less than 7% of the mitochondrial genome. The rapid development of technology for automated DNA sequencing has made it possible to study the complete mtDNA genomes in large numbers of individuals, opening the field of mitochondrial population genomics. Here we describe a suitable methodology for determining the complete human mitochondrial sequence and the global mtDNA diversity in humans. Also, we discuss the implications of the results with respect to the different hypotheses for the evolution of modern humans.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

A topoisomerase II cleavage site is associated with a novel mitochondrial DNA deletion

Mitochondrial myopathies and encephalopathies can be caused by nucleotide substitutions, deletions or duplications of the mitochondrial DNA (mtDNA). In one such disorder, Kearns-Sayre Syndrome (KSS), large-scale hetero-plasmic mtDNA deletions are often found. We describe a 14-year-old boy with clinical features of KSS, plus some additional features. Analysis of the entire mitochondrial genome by the polymerase chain reaction and Southern blotting revealed a 7864-bp mtDNA deletion, heteroplasmic in its tissue distribution. DNA sequencing established that the deletion was between nucleotides 6238 and 14103, and flanked by a 4-bp (TCCT) direct repeat sequence. Deletions between direct repeats have been hypothesised to occur by a slipped-mismatching or illegitimate recombination event, or following the DNA cleavage action of topoisomerase II. Analysis of the gene sequence in the region surrounding the mtDNA deletion breakpoint in this patient revealed the presence of putative vertebrate topoisomerase II sites. We suggest that direct repeat sequences, together with putative topoisomerase II sites, may predispose certain regions of the mitochondrial genome to deletions.     

Rapid identification of unknown heteroplasmic mutations across the entire human mitochondrial genome with mismatch-specific Surveyor Nuclease

Mitochondrial DNA (mtDNA) mutations are responsible for mitochondrial diseases in numerous patients. But, until now, no rapid method was available for the identification of unknown deleterious point mutations. Here, we describe a new strategy for the rapid identification of heteroplasmic mtDNA mutations using mismatch-specific Surveyor Nuclease. This protocol involves the following three steps: (i) PCR amplification of the entire human mitochondrial genome in 17 overlapping fragments; (ii) localization of mtDNA mismatch(es) after digestion of the 17 amplicons by Surveyor Nuclease; and (iii) identification of the mutation by sequencing the region containing the mismatch. This Surveyor Nuclease-based strategy allows a systematic screening of the entire mtDNA to identify a mutation within 2 days. It represents an important diagnostic approach for mitochondrial diseases that can be routinely used in molecular diagnostic laboratories.     

High mutation rates in the mitochondrial genomes of Daphnia pulex

Despite the great utility of mitochondrial DNA (mtDNA) sequence data in population genetics and phylogenetics, key parameters describing the process of mitochondrial mutation (e.g., the rate and spectrum of mutational change) are based on few direct estimates. Furthermore, the variation in the mtDNA mutation process within species or between lineages with contrasting reproductive strategies remains poorly understood. In this study, we directly estimate the mtDNA mutation rate and spectrum using Daphnia pulex mutation-accumulation (MA) lines derived from sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) lineages. The nearly complete mitochondrial genome sequences of 82 sexual and 47 asexual MA lines reveal high mtDNA mutation rate of 1.37 × 10(-7) and 1.73 × 10(-7) per nucleotide per generation, respectively. The Daphnia mtDNA mutation rate is among the highest in eukaryotes, and its spectrum is dominated by insertions and deletions (70%), largely due to the presence of mutational hotspots at homopolymeric nucleotide stretches. Maximum likelihood estimates of the Daphnia mitochondrial effective population size reveal that between five and ten copies of mitochondrial genomes are transmitted per female per generation. Comparison between sexual and asexual lineages reveals no statistically different mutation rates and highly similar mutation spectra.     

Modulating heteroplasmy

Patients with mitochondrial DNA (mtDNA) disease usually harbor a mixture of mutant and wild-type mtDNA (a state termed heteroplasmy), and the clinical features of the disease depend on the percentage of mutant mtDNA (the "mutation load") in vulnerable tissues. Factors that modulate the mutation load are poorly understood, but recent work has started to unravel the mechanisms. In certain circumstances heteroplasmy might be regulated at the level of the individual mitochondrial genome.     

The whole mitochondrial genome sequence of the zebra finch (Taeniopygia guttata)

Here we describe the complete nucleotide sequence of the mitochondrial genome (16 583/4 bp) of the zebra finch (Taeniopygia guttata). Primers were designed based on highly conserved regions of an alignment of three passerine complete mitochondrial DNA (mtDNA) sequences. A combination of overlapping long polymerase chain reaction (PCR) purification, followed by fully nested PCR and sequencing was used to determine the complete mtDNA genome. Six birds, from distinct maternal lineages of a pedigreed population were sequenced. Five novel haplotypes were identified. These sequences provide the first data for sequence variation across the whole mitochondrial genome of a passerine bird species.     

The mitochondrial genome sequence and molecular phylogeny of the turkey, Meleagris gallopavo

The mitochondrial genome (mtGenome) has been little studied in the turkey ( Meleagris gallopavo ), a species for which there is no publicly available mtGenome sequence. Here, we used PCR-based methods with 19 pairs of primers designed from the chicken and other species to develop a complete turkey mtGenome sequence. The entire sequence (16 717 bp) of the turkey mtGenome was obtained, and it exhibited 85% similarity to the chicken mtGenome sequence. Thirteen genes and 24 RNAs (22 tRNAs and 2 rRNAs) were annotated. An mtGenome-based phylogenetic analysis indicated that the turkey is most closely related to the chicken, Gallus gallus , and quail, Corturnix japonica . Given the importance of the mtGenome, the present work adds to the growing genomic resources needed to define the genetic mechanisms that underlie some economically significant traits in the turkey.     

The complete sequence of the mitochondrial genome of Saccharomyces cerevisiae

The currently available yeast mitochondrial DNA (mtDNA) sequence is incomplete, contains many errors and is derived from several polymorphic strains. Here, we report that the mtDNA sequence of the strain used for nuclear genome sequencing assembles into a circular map of 85 779 bp which includes 10 kb of new sequence. We give a list of seven small hypothetical open reading frames (ORFs). Hot spots of point mutations are found in exons near the insertion sites of optional mobile group I intron-related sequences. Our data suggest that shuffling of mobile elements plays an important role in the remodelling of the yeast mitochondrial genome.     

Application of high-resolution melting for genotyping bovine mitochondrial DNA

Recent studies have demonstrated that mitochondrial DNA (mtDNA) haplotype has a significant impact on the efficiency of bovine somatic cell nuclear transfer. Conventional methods for detecting mtDNA variations and haplotypes, such as restriction fragment length polymorphism (RFLP), temporal temperature gradient gel electrophoresis, dHPLC and sequencing, are labor intensive or expensive and have low sensitivity. High-resolution melting (HRM) analysis is a new technique for mutation detection and has the advantages of speed, cost, and accuracy. Here, we describe the genotyping of bovine mtDNA using HRM analysis. DNA samples containing mtDNA were extracted from 75 Holstein cows and subjected to rapid-cycle (<20 min) PCR of small amplicons (<120 bp) using specific primer sets. Capillaries containing the PCR products were then subjected to HRM analysis; data were acquired in 2 min and analyzed using the instrument's software. Five common bovine mtDNA single nucleotide polymorphisms were identified: 9602 G>A, 169 A>G, 166A>G with 173A>G, and 363C>G. These results agree with both sequencing and RFLP analysis. In addition, a very small amount of heteroplasmic variants (<5%) was sufficiently to be distinguished by HRM analysis that would be very useful to differentiate heteroplasmy vs. homoplasmy. HRM analysis thus provides a new approach to genotyping bovine mtDNA sequence variations and has many advantages over other methods, including speed of analysis, cost, and accuracy. We believe this will be a valuable technique for determining the efficiency of nuclear transfer in cloned embryos and for studying maternal effects on nuclear-cytoplasm interactions.     

Concerted action of two novel tRNA mtDNA point mutations in chronic progressive external ophthalmoplegia

CPEO (chronic progressive external ophthalmoplegia) is a common mitochondrial disease phenotype in adults which is due to mtDNA (mitochondrial DNA) point mutations in a subset of patients. Attributing pathogenicity to novel tRNA mtDNA mutations still poses a challenge, particularly when several mtDNA sequence variants are present. In the present study we report a CPEO patient for whom sequencing of the mitochondrial genome revealed three novel tRNA mtDNA mutations: G5835A, del4315A, T1658C in tRNATyr, tRNAIle and tRNAVal genes. In skeletal muscle, the tRNAVal and tRNAIle mutations were homoplasmic, whereas the tRNATyr mutation was heteroplasmic. To address the pathogenic relevance, we performed two types of functional tests: (i) single skeletal muscle fibre analysis comparing G5835A mutation loads and biochemical phenotypes of corresponding fibres, and (ii) Northern-blot analyses of mitochondrial tRNATyr, tRNAIle and tRNAVal. We demonstrated that both the G5835A tRNATyr and del4315A tRNAIle mutation have serious functional consequences. Single-fibre analyses displayed a high threshold of the tRNATyr mutation load for biochemical phenotypic expression at the single-cell level, indicating a rather mild pathogenic effect. In contrast, skeletal muscle tissue showed a severe decrease in respiratory-chain activities, a reduced overall COX (cytochrome c oxidase) staining intensity and abundant COX-negative fibres. Northern-blot analyses showed a dramatic reduction of tRNATyr and tRNAIle levels in muscle, with impaired charging of tRNAIle, whereas tRNAVal levels were only slightly decreased, with amino-acylation unaffected. Our findings suggest that the heteroplasmic tRNATyr and homoplasmic tRNAIle mutation act together, resulting in a concerted effect on the biochemical and histological phenotype. Thus homoplasmic mutations may influence the functional consequences of pathogenic heteroplasmic mtDNA mutations.