Use of stereotypical mutational motifs to define resolution limits for the ultra-deep resequencing of mitochondrial DNA

Massively parallel resequencing of mitochondrial DNA (mtDNA) has led to significant advances in the study of heteroplasmic mtDNA variants in health and disease, but confident resolution of very low-level variants (<2% heteroplasmy) remains challenging due to the difficulty in distinguishing signal from noise at this depth. However, it is likely that such variants are precisely those of greatest interest in the study of somatic (acquired) mtDNA mutations. Previous approaches to this issue have included the use of controls such as phage DNA and mtDNA clones, both of which may not accurately recapitulate natural mtDNA. We have therefore explored a novel approach, taking advantage of mtDNA with a known stereotyped mutational motif (nAT>C, from patient with MNGIE, mitochondrial neurogastrointestinal encephalomyopathy) and comparing mutational pattern distribution with healthy mtDNA by ligation-mediated deep resequencing (Applied Biosystems SOLiD). We empirically derived mtDNA-mutant heteroplasmy detection limits, demonstrating that the presence of stereotypical mutational motif could be statistically validated for heteroplasmy thresholds ≥0.22% (P=0.034). We therefore provide empirical evidence from biological samples that very low-level mtDNA mutants can be meaningfully resolved by massively parallel resequencing, confirming the utility of the approach for studying somatic mtDNA mutation in health and disease. Our approach could also usefully be employed in other settings to derive platform-specific deep resequencing resolution limits.     

Analyses of the mitochondrial mutations in the Chinese patients with sporadic Creutzfeldt–Jakob disease

Pathogenic mitochondrial DNA (mtDNA) mutations leading to mitochondrial dysfunction can cause a variety of chronic diseases in central nervous system (CNS). However, the role of mtDNA mutations in sporadic Creutzfeldt–Jakob disease (sCJD) has still been unknown. In this study, we comparatively analyzed complete mtDNA sequences of 31 Chinese sCJD patients and 32 controls. Using MITOMASTER and PhyloTree, we characterized 520 variants in sCJD patients and 507 variants in control by haplogroup and allele frequencies. We classified the mtDNAs into 40 sub-haplogroups of 5 haplogroups, most of them being Asian-specific haplogroups. Haplogroup U, an European-specific haplogroups mtDNA, was found only in sCJD. The analysis to control region (CR) revealed a 31% increase in the frequency of mtDNA CR mutations in sCJD versus controls. In functional elements of the mtDNA CR, six CR mutations were in conserved sequence blocks I (CSBI) in sCJD, while only one in control (P<0.05). More mutants in transfer ribonucleic acid-Leu (tRNA-Leu) were detected in sCJD. The frequencies of two synonymous amino-acid changes, m.11467A>G, p.(=) in NADH dehydrogenase subunit 4 (ND4) and m.12372G>A, p.(=) in NADH dehydrogenase subunit 5 (ND5), in sCJD patients were higher than that of controls. Our study, for the first time, screened the variations of mtDNA of Chinese sCJD patients and identified some potential disease-related mutations for further investigations.     

Effect of EME1 exon variant Ile350Thr on risk and early onset of breast cancer in southern Chinese women

Essential meiotic endonuclease 1 homolog 1 (EME1) is a key DNA repair protein that participates in the recognition and repair of DNA double-strand breaks. Deficiency of the EME1 gene can lead to spontaneous genomic instability and thus contribute to tumorgenesis. We hypothesized that the exon variants of EME1 confer genetic susceptibility to breast cancer. In a case-control study of 748 breast cancer patients and 778 normal controls, we analyzed the association between two exon variants of EME1 (i.e.,Ile350Thr: rs12450550T > C and Glu69Asp: rs3760413T > G) and breast cancer risk. We found that compared to the common Ile/Ile genotype, the Thr variant genotypes (Thr/Ile + Thr/Thr) conferred a 1.47-fold increased risk of breast cancer (OR=1.47, 95% CI=1.13-1.92). The variant Ile350Thr was also associated with early onset of breast cancer (r = -0.116, P = 0.002). The mean age of onset was 44.4 years for Thr/Thr genotype carriers and 46.5 years for Thr/Ile genotype carriers, which was significantly lower than that (49.4 years) for Ile/Ile genotype carriers (P = 0.006). Moreover, no significant association was observed between the Glu69Asp variant and breast cancer risk. Our findings suggest that the EME1 variant Ile350Thr contributes to an increased risk and early onset of breast cancer.     

Mutations in the rpoB gene of Mycobacterium tuberculosis that interfere with PCR-single-strand conformation polymorphism analysis for rifampin susceptibility testing.

Rifampin susceptibility of 32 rifampin-resistant and 26 rifampin-susceptible Mycobacterium tuberculosis strains was analyzed by PCR-single-strand conformation polymorphism (SSCP) and DNA sequencing within the 157-bp region of the rpoB gene (Ala500 to Val550). Two false-positive PCR-SSCP results were observed among the susceptible strains due to the silent mutation Gln513 (CAA-->CAG) and the deletion mutation Thr508 and Ser509. Another silent mutation [Leu511 (CTG-->CTA)], combined with the mutation Ser531-->Leu, was observed in a resistant strain. These results suggest that to rule out false-positive PCR-SSCP results, sequencing of the target DNA is required.     

线粒体耳聋与核修饰基因

线粒体突变是导致耳聋的重要原因,可以引起综合征和非综合征耳聋,这些突变主要位于线粒体12S rRNA和tRNA基因上。引起非综合征耳聋的突变主要以同质性形式存在,并且表型呈现高度异质性,说明其它因素参与了疾病的形成,包括环境因子,线粒体单倍型和核修饰基因。本文介绍了影响线粒体耳聋表型的4个候选核修饰基因:MTO1、GTPBP3、TFB1M和TRMU及其研究方法。     

Linkage disequilibrium between TAP2 variants and HLA class II alleles; no primary association between TAP2 variants and insulin-dependent diabetes mellitus

The TAP1 and TAP2 genes, located in the HLA class II region, encode subunits of a peptide transporter. Both genes display limited genetic variability; four different nucleotide substitutions have been found in theTAP2 gene. Here studies on linkage disequilibrium between TAP2 variants and HLA class II alleles are reported, in an attempt to evaluate whether TAP2 variants are associated with insulin-dependent diabetes mellitus (IDDM). As reported previously, a significant decrease of homozygosity for TAP2 alleles encoding alanine at residue 665 (665 Ala) and glut amine at 687 (687 Gin) paralleled by an increase in homozygosity for TAP2 alleles encoding threonine at residue 665 (665 Thr) and a stop codon at 687 (687 Stop), was found in both Finnish and Norwegian IDDM patients compared to random controls. However, a strong linkage disequilibrium between these TAP2 polymorphisms and given HLA-DR and -DQ genes was observed among healthy controls. The frequent 665 Thr and 687 Stop variants were in linkage disequilibrium both with the DR4-DQ8 and the DR3-DQ2 haplotypes, haplotypes which are strongly associated with IDDM. In contrast, the DR1-DQ5 and DR13-DQ6 (e.g. DQB1*0603) haplotypes, which are decreased among IDDM patients, were associated with the 665 Ala and 687 Gin variants. Thus, when DR- and DQ-matched patients and controls were compared, associations of the investigated TAP2 variants and IDDM were no longer detectable. These data, therefore, indicate that the associations previously found between certain TAP2 variants and IDDM are secondary to a primary association between this disease and particular DQαβ heterodimers.     

Mitochondrial DNA variants in inclusion body myositis characterized by deep sequencing

Muscle pathology in inclusion body myositis (IBM) typically includes inflammatory cell infiltration, muscle fibers with rimmed vacuoles and cytochrome c oxidase (COX)‐deficient fibers. Previous studies have revealed clonal expansion of large mitochondrial DNA (mtDNA) deletions in the COX‐deficient muscle fibers. Technical limitations have prevented complete investigations of the mtDNA deletions and other mtDNA variants. Detailed characterization by deep sequencing of mtDNA in muscle samples from 21 IBM patients and 10 age‐matched controls was performed after whole genome sequencing with a mean depth of mtDNA coverage of 46,000x. Multiple large mtDNA deletions and duplications were identified in all IBM and control muscle samples. In general, the IBM muscles demonstrated a larger number of deletions and duplications with a mean heteroplasmy level of 10% (range 1%‐35%) compared to controls (1%, range 0.2%‐3%). There was also a small increase in the number of somatic single nucleotide variants in IBM muscle. More than 200 rearrangements were recurrent in at least two or more IBM muscles while 26 were found in both IBM and control muscles. The deletions and duplications, with a high recurrence rate, were mainly observed in three mtDNA regions, m.534‐4429, m.6330‐13993, and m.8636‐16072, where some were flanked by repetitive sequences. The mtDNA copy number in IBM muscle was reduced to 42% of controls. Immunohistochemical and western blot analyses of IBM muscle revealed combined complex I and complex IV deficiency affecting the COX‐deficient fibers. In conclusion, deep sequencing and quantitation of mtDNA variants revealed that IBM muscles had markedly increased levels of large deletions and duplications, and there were also indications of increased somatic single nucleotide variants and reduced mtDNA copy numbers compared to age‐matched controls. The distribution and type of variants were similar in IBM muscle and controls indicating an accelerated aging process in IBM muscle, possibly associated with chronic inflammation.     

Fine mapping analysis confirms and strengthens linkage of four chromosomal regions in familial hypospadias

Hypospadias is a common male genital malformation and is regarded as a complex disease affected by multiple genetic as well as environmental factors. In a previous genome-wide scan for familial hypospadias, we reported suggestive linkage in nine chromosomal regions. We have extended this analysis by including new families and additional markers using non-parametric linkage. The fine mapping analysis displayed an increased LOD score on chromosome 8q24.1 and 10p15 in altogether 82 families. On chromosome 10p15, with the highest LOD score, we further studied AKR1C2, AKR1C3 and AKR1C4 involved in steroid metabolism, as well as KLF6 expressed in preputial tissue from hypospadias patients. Mutation analysis of the AKR1C3 gene showed a new mutation, c.643G>A (p.(Ala215Thr)), in a boy with penile hypospadias. This mutation is predicted to have an impact on protein function and structure and was not found in controls. Altogether, we homed in on four chromosomal regions likely to harbor genes for hypospadias. Future studies will aim for studying regulatory sequence variants in these regions.     

Variation and association to diabetes in 2000 full mtDNA sequences mined from an exome study in a Danish population

In this paper, we mine full mtDNA sequences from an exome capture data set of 2000 Danes, showing that it is possible to get high-quality full-genome sequences of the mitochondrion from this resource. The sample includes 1000 individuals with type 2 diabetes and 1000 controls. We characterise the variation found in the mtDNA sequence in Danes and relate the variation to diabetes risk as well as to several blood phenotypes of the controls but find no significant associations. We report 2025 polymorphisms, of which 393 have not been reported previously. These 393 mutations are both very rare and estimated to be caused by very recent mutations but individuals with type 2 diabetes do not possess more of these variants. Population genetics analysis using Bayesian skyline plot shows a recent history of rapid population growth in the Danish population in accordance with the fact that >40% of variable sites are observed as singletons.     

Candidacidal activity of recombinant human salivary histatin-5 and variants.

Human salivary histatins possess fungicidal and bactericidal activities. The current investigation evaluates the structure-function relationship of histatins with regard to their candidacidal activity by using recombinant histatin-5 and its variants produced in Escherichia coli. The purified recombinant histatins were examined for their candidacidal activity and secondary structure. The m21 (with Lys-13 replaced by Thr [Lys-13-->Thr]) and m71 (Lys-13-->Glu) variants are significantly less effective than recombinant histatin-5 in killing Candida albicans, suggesting that Lys-13 is critical for candidacidal activity. The m68 (Lys-13-->Glu and Arg-22-->Gly) variant is significantly less potent than the recombinant histatin-5 as well as m71, indicating that Arg-22 is crucial for the cidal activity. The candidacidal activities of m1 (Arg-12-->Ile), m2 (Arg-12-->Ile and Lys-17-->Asp), m12 (Arg-12-->Lys and His-21-->Leu), and m70 (His-19-->Pro and His-21-->Arg) variants, however, are comparable to that of recombinant histatin-5, indicating that Arg-12, Lys-17, His-19, and His-21 are not functionally important. The conformational preferences of histatin-5 and variants were determined by circular dichroism. The results indicate that all proteins have a strong tendency to adopt alpha-helical conformation in trifluoroethanol. Previously, we have shown that the alpha-helical conformation is one of the important structural requirements for eliciting appreciable candidacidal activity. Collectively, the data suggest that in addition to the helical conformation, specific residues such as Lys-13 and Arg-22 in the sequence of histatin-5 are, indeed, important for candidacidal activity.     

Polymorphisms of the TAP1 and TAP2 genes in human alveolar echinococcosis

We postulated that TAP genes may influence the susceptibility of some individuals to Echinococcus multilocularis infection. Six coding region variants (codons 333 and 637 in TAP1, and 379, 565, 651 and 665 in TAP2) were typed in 94 patients and 100 controls. Thr/Thr homozygosity at TAP2/665 was more prevalent in patients than in controls [64% vs. 45%, respectively; odds ratio (OR) = 2.1 (95% confidence interval (CI) 1.1; 2.7)] and Thr/Ala heterozygozity was less prevalent (32% vs. 50%, respectively) (P = 0.014). Of the 38 patients with progressive lesions, 76% were Thr/Thr, as compared with 55% of patients without progressive lesions and 45% of controls (P = 0.058 and 0.02, respectively), independent of HLA status. To determine whether this association is functionally relevant, functional analyses and/or confirmation in distinct populations of patients with alveolar echinococcosis would be required.     

Phylogenetic analysis of the mitochondrial genome indicates significant differences between patients with Alzheimer disease and controls in a French-Canadian founder population.

The activity of cytochrome oxidase (CO), the terminal enzyme of the mitochondrial electron transport chain, has been reported to be lower in the brains of Alzheimer disease (AD) patients. This suggests that a modification of mitochondrial DNA (mtDNA) may be responsible for this decrease of CO activity. Many mtDNA variants were found by different studies at a higher frequency in AD patients, suggesting that mtDNA variants could confer a genetic susceptibility to AD. In this study, we sequenced the entire mitochondrial genome region that encompasses the three CO genes and the 22 mitochondrial tRNA in 69 AD patients and 83 age-matched controls. We detected a total of 95 mtDNA variants. The allele frequencies of the majority of these variants were similar in patients and controls. However, a haplotype composed of three different modifications (positions: 5633, 7476, and 15812) was present in three of the 69 late-onset AD patients (4.3%) and also in 1 of 16 early-onset AD patients (6.2%) but not in control individuals. Given that one of these variants (15812) has already been shown to be associated with another neurodegenerative disease and that all three modifications are relatively conserved and their frequencies in the general population is only 0.1%, our data suggest that the presence of this haplotype may represent a risk factor for AD. We also found a significant association (P < 0.05) of two other variants at positions 709 (rRNA 12S) and 15928 (tRNA(Thr)). These two mtDNA variants are three times more frequent in control individuals compared with AD patients, suggesting that they may be protective against AD.