Length Heteroplasmy Research Articles

High heteroplasmy is associated with low mitochondrial copy number and selection against non-synonymous mutations in the snail Cepaea nemoralis

Molluscan mitochondrial genomes are unusual because they show wide variation in size, radical genome rearrangements and frequently show high variation (> 10%) within species. As progress in understanding this variation has been limited, we used whole genome sequencing of a six-generation matriline of the terrestrial snail Cepaea nemoralis, as well as whole genome sequences from wild-collected C. nemoralis, the sister species C. hortensis, and multiple other snail species to explore the origins of mitochondrial DNA (mtDNA) variation. The main finding is that a high rate of SNP heteroplasmy in somatic tissue was negatively correlated with mtDNA copy number in both Cepaea species. In individuals with under ten mtDNA copies per nuclear genome, more than 10% of all positions were heteroplasmic, with evidence for transmission of this heteroplasmy through the germline. Further analyses showed evidence for purifying selection acting on non-synonymous mutations, even at low frequency of the rare allele, especially in cytochrome oxidase subunit 1 and cytochrome b. The mtDNA of some individuals of Cepaea nemoralis contained a length heteroplasmy, including up to 12 direct repeat copies of tRNA-Val, with 24 copies in another snail, Candidula rugosiuscula, and repeats of tRNA-Thr in C. hortensis. These repeats likely arise due to error prone replication but are not correlated with mitochondrial copy number in C. nemoralis. Overall, the findings provide key insights into mechanisms of replication, mutation and evolution in molluscan mtDNA, and so will inform wider studies on the biology and evolution of mtDNA across animal phyla.

Population analysis of complete mitogenomes for 334 samples from El Salvador

The use of mitochondrial DNA (mtDNA) in the field of forensic genetics is widely spread mainly due to its advantages when identifying highly degraded samples. In this sense, massive parallel sequencing has made the analysis of the whole mitogenome more accessible, noticeably increasing the informativeness of mtDNA haplotypes. The civil war (1980–1992) in El Salvador caused many deaths and disappearances (including children) all across the country and the economic and social instability after the war forced many people to emigration. For this reason, different organizations have collected DNA samples from relatives with the aim of identifying missing people. Thus, we present a dataset containing 334 complete mitogenomes from the Salvadoran general population. To the best of our knowledge, this is the first publication of a nationwide forensic-quality complete mitogenome database of any Latin American country. We found 293 different haplotypes, with a random match probability of 0.0041 and 26.6 mean pairwise differences, which is similar to other Latin American populations, and which represent a marked improvement from the values obtained with just control region sequences. These haplotypes belong to 54 different haplogroups, being 91% of them of Native American origin. Over a third (35.9%) of the individuals carried at least a heteroplasmic site (excluding length heteroplasmies). Ultimately, the present database aims to represent mtDNA haplotype diversity in the general Salvadoran populations as a basis for the identification of people that disappeared during or after the civil war.

Whole Mitochondrial Genome Detection and Analysis of Two- to Four-Generation Maternal Pedigrees Using a New Massively Parallel Sequencing Panel.

Mitochondrial DNA (mtDNA) is an effective genetic marker in forensic practice, especially for aged bones and hair shafts. Detection of the whole mitochondrial genome (mtGenome) using traditional Sanger-type sequencing is laborious and time-consuming. Additionally, its ability to distinguish point heteroplasmy (PHP) and length heteroplasmy (LHP) is limited. The application of massively parallel sequencing in mtDNA detection helps researchers to study the mtGenome in-depth. The ForenSeq mtDNA Whole Genome Kit, which contains a total of 245 short amplicons, is one of the multiplex library preparation kits for the mtGenome. We used this system to detect the mtGenome in the blood samples and hair shafts of thirty-three individuals from eight two-generation pedigrees, one three-generation pedigree, and one four-generation pedigree. High-quality sequencing results were obtained. Ten unique mtGenome haplotypes were observed in the mothers from the ten pedigrees. A total of 26 PHPs were observed using the interpretation threshold of 6%. Eleven types of LHPs in six regions were evaluated in detail. When considering homoplasmic variants only, consistent mtGenome haplotypes were observed between the twice-sequenced libraries and between the blood and hair shafts from the same individual and among maternal relatives in the pedigrees. Four inherited PHPs were observed, and the remainder were de novo/disappearing PHPs in the pedigrees. Our results demonstrate the effective capability of the ForenSeq mtDNA Whole Genome Kit to generate the complete mtGenome in blood and hair shafts, as well as the complexity of mtDNA haplotype comparisons between different types of maternal relatives when heteroplasmy is considered.

Forensic evaluation of mitochondrial DNA heteroplasmy in Gujarat population, India

Background Owing to its high copy number and its small size, mtDNA analysis is the most reliable choice when biological materials from crime scenes are degraded or have mixed STR profiles. Aim To examine the occurrence of heteroplasmy along with its frequency and pattern in both HV1 and HV2 regions of the mtDNA among unrelated individuals from India. Subjects and methods Mitochondrial DNA control region [hypervariable region one (HV1) and hypervariable region two (HV2)] were analysed in blood and buccal tissues of 104 unrelated individuals from the Indian state of Gujarat. Results A high frequency of point heteroplasmy (PH) and length heteroplasmy (LH) was revealed. PH was detected in 7.69% of the population, with a higher frequency observed in blood than in buccal samples. However, there were no statistically significant differences in PH between the two tissues (Chi-square = 0.552, p ≥ 0.05). A total of six PH positions were detected: three at HV1, and another three at HV2. The studied population showed 46.15% LH in the HV1 and HV2 regions of both tissues. The LH positions observed in the Gujarat population were the same as those previously reported at HV1 np16184–16193 and HV2 np303–315. Conclusions Our findings suggest that differences in the pattern of heteroplasmy found in different tissues can complicate the forensic analysis, on the other hand, the probability of a match between the questioned and reference samples increases when the heteroplasmy is identical in both tissues. Variability of PH among persons and even within tissues recommends analysing multiple tissue samples before drawing a conclusion in forensic mtDNA analyses.

Development of a multiplex assay for detection of autosomal and Y-chromosomal STRs, assessment of the degradation state of mitochondrial DNA and presence of mitochondrial length heteroplasmies

The current focus in most routine forensic casework is detection of autosomal or gonosomal Short Tandem Repeats (STRs). With increasing degradation, STR analysis tends to be less successful up to complete failure. For challenging samples such as telogen hair roots and shafts, touch DNA samples or skeletal remains, mitochondrial DNA (mtDNA) analysis provides a powerful tool. Determination of DNA quantity is an important part in the casework workflow. Several ready-to-use kits are commercially available for nuclear DNA targets. However, quantification of mtDNA targets requires the establishment of an in-house method. Some assays even contain assessment of degradation, which alleviates the choice of target enrichment for sequencing through medium or small amplicons. As Sanger-type Sequencing (STS) still remains the golden standard in many laboratories, identification of heteroplasmies in C-tract regions prior to the sequencing reaction is advantageous. Firstly, primer selection can be expanded with primers binding near the C-tract and secondly, determination of the dominant variant is straightforward. All those quantity (nuclear and mtDNA) and quality (degradation and length heteroplasmies) evaluations usually require at least two separate reactions. Therefore, the aim of this project was the combination of all these targets in one multiplex assay using capillary electrophoresis to spare valuable sample extract. Amplification of representative autosomal and Y-chromosomal STRs allows estimate of success of (Y-)STR analysis. Simultaneously, five length heteroplasmies in the mitochondrial control region are targeted as well as three conservative regions of differing fragment lengths for assessment of the mitochondrial degradation state. Based on the outcome of this assay, forensic examiners can decide if STR analysis may be suitable. In case of absent STR peaks, appropriate proceeding of mtDNA sequencing can be determined.

Platinum-Quality Mitogenome Haplotypes from United States Populations.

A total of 1327 platinum-quality mitochondrial DNA haplotypes from United States (U.S.) populations were generated using a robust, semi-automated next-generation sequencing (NGS) workflow with rigorous quality control (QC). The laboratory workflow involved long-range PCR to minimize the co-amplification of nuclear mitochondrial DNA segments (NUMTs), PCR-free library preparation to reduce amplification bias, and high-coverage Illumina MiSeq sequencing to produce an average per-sample read depth of 1000 × for low-frequency (5%) variant detection. Point heteroplasmies below 10% frequency were confirmed through replicate amplification, and length heteroplasmy was quantitatively assessed using a custom read count analysis tool. Data analysis involved a redundant, dual-analyst review to minimize errors in haplotype reporting with additional QC checks performed by EMPOP. Applying these methods, eight sample sets were processed from five U.S. metapopulations (African American, Caucasian, Hispanic, Asian American, and Native American) corresponding to self-reported identity at the time of sample collection. Population analyses (e.g., haplotype frequencies, random match probabilities, and genetic distance estimates) were performed to evaluate the eight datasets, with over 95% of haplotypes unique per dataset. The platinum-quality mitogenome haplotypes presented in this study will enable forensic statistical calculations and thereby support the usage of mitogenome sequencing in forensic laboratories.

The mitochondrial genome sequence of the BS-C-1 cell line is at odds with the reported derivation from Cercopithecus aethiops

The BS-C-1 cell line widely used in virological research was reportedly derived from the African green monkey Cercopithecus (Chlorocebus) aethiops. We used Sanger sequencing to determine the full nucleotide sequence of mtDNA in BS-C-1 cells. The mitochondrial genome in this cell line is 16,456 bp long and has an organization typical of that in other mammalian genomes with 37 genes, including 2 rRNAs, 22 tRNAs, and 13 polypeptide-encoding genes. Surprisingly, this genome is most closely related to that of Chlorocebus pygerythrus, which does not agree with the reported derivation of this cell line from C. aethiops. Another unexpected observation is the presence of mtDNA length heteroplasmy within the MT-ND6 gene, with the variant containing a stretch of 8 Cs encoding a functional gene, whereas the variant containing 9 Cs encodes a frameshifted gene. It is unclear whether the latter variant is nonfunctional or whether it is corrected by programmed translational ribosomal −1 frameshifting. The availability of the full mtDNA sequence for the BS-C-1 cell line should increase its utility by enabling studies on mtDNA transcription and replication.

Unusual mtDNA Control Region Length Heteroplasmy in the COS-7 Cell Line.

The COS-7 cell line is a workhorse of virology research. To expand this cell line’s utility and to enable studies on mitochondrial DNA (mtDNA) transcription and replication, we determined the complete nucleotide sequence of its mitochondrial genome by Sanger sequencing. In contrast to other available mtDNA sequences from Chlorocebus aethiops, the mtDNA of the COS-7 cell line was found to contain a variable number of perfect copies of a 108 bp unit tandemly repeated in the control region. We established that COS-7 cells are heteroplasmic with at least two variants being present: with four and five repeat units. The analysis of the mitochondrial genome sequences from other primates revealed that tandem repeats are absent from examined mtDNA control regions of humans and great apes, but appear in lower primates, where they are present in a homoplasmic state. To our knowledge, this is the first report of mtDNA length heteroplasmy in primates.

Impact of the sequencing method on the detection and interpretation of mitochondrial DNA length heteroplasmy

Advancements in sequencing technologies allow for rapid and efficient analysis of mitochondrial DNA (mtDNA) in forensic laboratories, which is particularly beneficial for specimens with limited nuclear DNA. Next generation sequencing (NGS) offers higher throughput and sensitivity over traditional Sanger-type sequencing (STS) as well as the ability to quantitatively analyze the data. Changes in sample preparation, sequencing method and analysis required for NGS may alter the mtDNA haplotypes compared to previously generated STS data. Thus, the present study aimed to characterize the impact of different sequencing workflows on the detection and interpretation of length heteroplasmy (LHP), a particularly complicated aspect of mtDNA analysis. Whole mtDNA genome (mitogenome) data were generated for 16 high-quality samples using well-established Illumina and Ion methods, and the NGS data were compared to previously-generated STS mtDNA control region data. Although the mitogenome haplotypes were concordant with the exception of length and low-level variants (<30 % variant frequency), LHP in the hypervariable segment (HVS) polycytosine regions (C-tracts) differed across sequencing methods. Consistent with previous studies, LHP in HVS1 was observed in samples with nine or more consecutive cytosines (Cs) and eight Cs in the HVS2 region in the STS data. The Illumina data produced a similar pattern of LHP as the STS data, whereas the Ion data were noticeably different. More complex LHP (i.e. more length molecules) was observed in the Ion data, as length variation occurred in multiple homopolymer stretches within the targeted HVS regions. Further, the STS dominant or major molecule (MM) differed from the Ion MM in 11 (37 %) of the 30 regions evaluated and six instances (20 %) in Illumina data. This is of particular interest, as the MM is used by many forensic laboratories to report the HVS C-tract in the mtDNA haplotype. In general, the STS MMs were longer than the Illumina MMs, while the Ion MMs were the shortest. The higher rate of homopolymer indels in Ion data likely contributed to these differences. Supplemental analysis with alternative approaches demonstrated that the LHP pattern may also be altered by the bioinformatic tool and workflow used for data interpretation. The broader application of NGS in forensic laboratories will undoubtedly result in the use of varying sample preparation and sequencing methods. Based on these findings, minor LHP differences are expected across sequencing workflows, and it will be important that C-tract indels continue to be ignored for forensic queries and comparisons.

Length heteroplasmy in the predominate mitochondrial DNA haplogroups in the Croatian population

Mitochondrial control region represents the most variable segment of the mitochondrial genome. The frequency and pattern of heteroplasmy has been described in several studies; however, none of the reports documented the Croatian population. In the present study, we screened the control region (1122 bp) of 95 individuals belonging to two predominant mitochondrial phylogenetic branches in the Croatian population, haplogroups H and U. Length heteroplasmy occurred in polycytosine (poly-C) tracts within three hypervariable segments of the control region with the following frequencies: HVSI - 26.3%, HVSII - 52.6% and HVSIII - 7.4%. Furthermore, the association between certain polymorphisms in HVSI and length heteroplasmy was investigated. Our results indicate that only polymorphisms located in the poly-C tract are associated with HVSI length heteroplasmy. The T to C transition at np 16189 is significantly associated with the occurrence of length heteroplasmy (p&amp;lt;0.0001). This effect was even stronger if the C insertion was present in the position 16193. The data support the hypothesis that an uninterrupted poly-C tract of more than eight cytosines leads to length heteroplasmy. Length heteroplasmy associated with the T to C substitution in np 16189 was predominantly found in haplogroup U.

Spectrum of mitochondrial genome instability and implication of mitochondrial haplogroups in Korean patients with acute myeloid leukemia

BackgroundMitochondrial DNA (mtDNA) mutations may regulate the progression and chemosensitivity of leukemia. Few studies regarding mitochondrial aberrations and haplogroups in acute myeloid leukemia (AML) and their clinical impacts have been reported. Therefore, we focused on the mtDNA length heteroplasmies minisatellite instability (MSI), copy number alterations, and distribution of mitochondrial haplogroups in Korean patients with AML.MethodsThis study investigated 74 adult patients with AML and 70 controls to evaluate mtDNA sequence alterations, MSI, mtDNA copy number, haplogroups, and their clinical implications. The hypervariable (HV) control regions (HV1 and HV2), tRNAleu1gene, and cytochrome b gene of mtDNA were analyzed. Two mtDNA minisatellite markers, 16189 poly-C (16184CCCCCTCCCC16193, 5CT4C) and 303 poly-C (303CCCCCCCTCCCCC315, 7CT5C), were used to examine the mtDNA MSI.ResultsIn AML, most mtDNA sequence variants were single nucleotide substitutions, but there were no significant differences compared to those in controls. The number of mtMSI patterns increased in AML. The mean mtDNA copy number of AML patients increased approximately 9-fold compared to that of controls (P<0.0001). Haplogroup D4 was found in AML with a higher frequency compared to that in controls (31.0% vs. 15.7%, P=0.046). None of the aforementioned factors showed significant impacts on the outcomes.ConclusionAML cells disclosed more heterogeneous patterns with the mtMSI markers and had increased mtDNA copy numbers. These findings implicate mitochondrial genome instability in primary AML cells. Therefore, mtDNA haplogroup D4 might be associated with AML risk among Koreans.

The complete mitochondrial genome of the grooved carpet shell, Ruditapes decussatus (Bivalvia, Veneridae).

Despite the large number of animal complete mitochondrial genomes currently available in public databases, knowledge about mitochondrial genomics in invertebrates is uneven. This paper reports, for the first time, the complete mitochondrial genome of the grooved carpet shell, Ruditapes decussatus, also known as the European clam. Ruditapes decussatus is morphologically and ecologically similar to the Manila clam Ruditapes philippinarum, which has been recently introduced for aquaculture in the very same habitats of Ruditapes decussatus, and that is replacing the native species. Currently the production of the European clam is almost insignificant, nonetheless it is considered a high value product, and therefore it is an economically important species, especially in Portugal, Spain and Italy. In this work we: (i) assembled Ruditapes decussatus mitochondrial genome from RNA-Seq data, and validated it by Sanger sequencing; (ii) analyzed and characterized the Ruditapes decussatus mitochondrial genome, comparing its features with those of other venerid bivalves; (iii) assessed mitochondrial sequence polymorphism (SP) and copy number variation (CNV) of tandem repeats across 26 samples. Despite using high-throughput approaches we did not find evidence for the presence of two sex-linked mitochondrial genomes, typical of the doubly uniparental inheritance of mitochondria, a phenomenon known in ∼100 bivalve species. According to our analyses, Ruditapes decussatus is more genetically similar to species of the Genus Paphia than to the congeneric Ruditapes philippinarum, a finding that bolsters the already-proposed need of a taxonomic revision. We also found a quite low genetic variability across the examined samples, with few SPs and little variability of the sequences flanking the control region (Largest Unassigned Regions (LURs). Strikingly, although we found low nucleotide variability along the entire mitochondrial genome, we observed high levels of length polymorphism in the LUR due to CNV of tandem repeats, and even a LUR length heteroplasmy in two samples. It is not clear if the lack of genetic variability in the mitochondrial genome of Ruditapes decussatus is a cause or an effect of the ongoing replacement of Ruditapes decussatus with the invasive Ruditapes philippinarum, and more analyses, especially on nuclear sequences, are required to assess this point.

Investigation of length heteroplasmy in mitochondrial DNA control region by massively parallel sequencing

Accurate sequencing of the control region of the mitochondrial genome is notoriously difficult due to the presence of polycytosine bases, termed C-tracts. The precise number of bases that constitute a C-tract and the bases beyond the poly cytosines may not be accurately defined when analyzing Sanger sequencing data separated by capillary electrophoresis. Massively parallel sequencing has the potential to resolve such poor definition and provides the opportunity to discover variants due to length heteroplasmy. In this study, the control region of mitochondrial genomes from 20 samples was sequenced using both standard Sanger methods with separation by capillary electrophoresis and also using massively parallel DNA sequencing technology. After comparison of the two sets of generated sequence, with the exception of the C-tracts where length heteroplasmy was observed, all sequences were concordant. Sequences of three segments 16184–16193, 303–315 and 568–573 with C-tracts in HVI, II and III can be clearly defined from the massively parallel sequencing data using the program SEQ Mapper. Multiple sequence variants were observed in the length of C-tracts longer than 7 bases. Our report illustrates the accurate designation of all the length variants leading to heteroplasmy in the control region of the mitochondrial genome that can be determined by SEQ Mapper based on data generated by massively parallel DNA sequencing.

Problems in Mitochondrial DNA forensics: while interpreting length heteroplasmy conundrum of various Sindhi and Baluchi ethnic groups of Pakistan

The insight heterodox genetics of mtDNA infer new perspectives at the level of human mitochondrial control region heteroplasmy, which is substantial in evolutionary as well as forensic interpretation. The main goal of this study is to interrogate the recurrence and resolve the ambiguity of blurry spectrum of heteroplasmy in the human mtDNA control region of 50 Baluchi and 116 Sindhi unrelated individuals. Sanger sequencing was employed classically, that was further investigated by minisequencing. Only 20% Baluchi and 25.8% Sindhi were homoplasmic, whereas rest of 80% Baluchi and 74.1% Sindhi exhibited at least one heteroplasmy within the specimen. In total, 166 individuals have length heteroplasmy (LH) found at positions 16189, 303–315, 568–573, and 514–524, whilst point mutation heteroplasmy (PMH) was detected at positions 73, 16093, 16189, and 16234, respectively. Overall LH was observed albeit high frequency in Sindhi ethnic group (82%) rather than Baluchi’s (37%), whereas PMH accumulation was relatively extensive (24%) in Baluchi’s than Sindhi’s (11.2%). The obtained results ascertained that growing knowledge of heteroplasmy assisted to develop consciences in the forensic community that heteroplasmy plays a pivotal role in the legal interpretation on a regular basis and knowledge of its biological underpinnings has a vital niche in the forensic science. Limited studies have focused on heteroplasmy, yet scientific attention should be given, in order to determine its magnitude in different ethnic boundaries.

High occurrence of length heteroplasmy in domestic Bactrian camel (Camelus bactrianus)

Heteroplasmy is the presence of more than one mitochondrial DNA (mtDNA) variant within a cell, tissue, or individual. In this study, sequence variation was investigated in the control region (CR) of mitochondrial DNA (mtDNA) from 135 individuals belonging to five primary domestic Bactrian camel breeds in China and Mongolia. Due to variation of the repeat unit G(T/C)(AC)n, length heteroplasmy was detected within each camel by direct sequencing and fragment analysis. A high occurrence of mtDNA heteroplasmy, up to 100 percentages was observed in five camel populations. Our study provides the first evidence of extensive length heteroplasmy in Bactrian camels.

Defect of the mitochondrial DNA hypervariable region as a risk factor for canine mammary tumour

The aim of this study was to identify mutations in the hypervariable region of mitochondrial DNA in canine mammary tumours and to determine their association with the process of neoplastic transformation. A total of 93 biological samples, including blood as well as normal and neoplastic tissue samples from 31 dogs with diagnosed malignant canine mammary tumours were analysed. DNA extraction, amplification and sequencing of the D-loop as well as bioinformatic and statistical analyses were performed. In the mitochondrial D-loop sequence, 26 polymorphic loci and 5 mutations were identified. For the first time, D-loop length heteroplasmy was detected in dogs with mammary tumours. The malignancy grade exerted no effect on the presence of nucleotide changes. A statistically significant association between the presence of mutations and polymorphisms and the size of dogs was demonstrated. The 100% frequency of length heteroplasmy may imply that this is a hotspot mutation of canine mammary tumour.

Analysis of heteroplasmy in the major noncoding region of mitochondrial DNA in the blood and atherosclerotic plaques of carotid arteries

For identification of somatic mitochondrial DNA (mtDNA) mutations, the mtDNA major noncoding region (D-loop) sequence in blood samples and carotid atherosclerosis plaques from patients with atherosclerosis was analyzed. Five point heteroplasmic positions were observed in 4 of 23 individuals (17%). Only in two cases could heteroplasmy have resulted from somatic mutation, whereas three heteroplasmic positions were found in both vascular tissue and blood. In addition, length heteroplasmy in a polycytosine stretches was registered at nucleotide positions 303-315 in 16 individuals, and also in the 16184-16193 region--in four patients. The results suggest that somatic mtDNA mutations can occur during atherosclerosis, but some heteroplasmic mutations may appear in all tissues, possibly being inherited.

Investigation into length heteroplasmy in the mitochondrial DNA control region after treatment with bisulfite.

We report on a method to analyze length heteroplasmy within the human mitochondrial genome in which there are polycytosine [poly(C)] stretches. These poly(C) tracts induce heteroplasmy with the resultant inherent problems of accurate sequence designations. In this study, 20 samples that exhibited length heteroplasmy due to variation in the C-tracts within hypervariable region I (HVI) were treated with bisulfite, and one or more cytosine bases in these C-tracts were converted randomly to uracil. This resulted in an accurate sequence designation for nearly all samples. The only exceptions in which the DNA sequence could still not be determined occurred when there was total conversion, or a lack of conversion, of the cytosine bases. Replicate tests on the same samples showed that individual cytosine bases were randomly converted to uracil. This simple method was useful for investigating length heteroplasmy due to 16189C and 310C transitions in the mitochondrial-DNA control region. It is valuable for medical and forensic investigations.

Peripheral blood mitochondrial DNA copy number, length heteroplasmy and breast cancer risk: a replication study.

Oxidative stress has consistently been linked to breast carcinogenesis, and mitochondria play a significant role in regulating reactive oxygen species generation. In our previous study, we found that increased levels of mitochondrial DNA (mtDNA) copy number and the presence of mitochondrial length heteroplasmies in the hypervariable (HV) regions 1 and 2 (HV1 and HV2) in peripheral blood are associated with increased risk of breast cancer. In current study with 1000 breast cancer cases and 1000 healthy controls, we intended to replicate our previous findings. Overall, levels of mtDNA copy number were significantly higher in breast cancer cases than healthy controls (mean: 1.17 versus 0.94, P < 0.001). In the multivariate linear regression analysis, increased mtDNA copy number levels were associated with a 1.32-fold increased risk of breast cancer [adjusted odds ratio (OR) = 1.32, 95% confidence interval (CI) = 1.15-1.67]. Breast cancer cases were more likely to have HV1 and HV2 region length heteroplasmies than healthy controls (P < 0.001, respectively). The existence of HV1 and HV2 length heteroplasmies was associated with 2.01- and 1.63-folds increased risk of breast cancer (for HV1: OR = 2.01, 95% CI = 1.66-2.42; for HV2: OR = 1.63, 95% CI = 1.34-1.92). Additionally, joint effects among mtDNA copy number, HV1 and HV2 length heteroplasmies were observed. Our results are consistent with our previous findings and further support the roles of mtDNA copy number and mtDNA length heteroplasmies that may play in the development of breast cancer.

Length heteroplasmy of the polyC-polyT-polyC stretch in the dog mtDNA control region.

Previously, the mitochondrial control region of 214 Belgian dogs was sequenced. Analysis of this data indicated length heteroplasmy of the polyT stretch in the polyC-polyT-polyC stretch from positions 16661 to 16674. Nine polyC-polyT-polyC haplotype combinations were observed, consisting of seven major haplotypes (highest signal intensity) combined with minor haplotypes (lower signal intensity) one T shorter than the major haplotype in all but three dogs. The longer the polyT stretch, the smaller was the difference in signal intensity between the major and minor haplotype peaks. Additional sequencing, cloning, and PCR trap experiments were performed to further study the intra-individual variation of this mitochondrial DNA (mtDNA) region. Cloning experiments demonstrated that the proportion of clones displaying the minor haplotypes also increased with the length of the polyT stretch. Clone amplification showed that in vitro polymerase errors might contribute to the length heteroplasmy of polyT stretches with at least 10 Ts. Although major and minor polyC-polyT-polyC haplotypes did not differ intra-individually within and between tissues in this study, interpretation of polyT stretch variation should be handled with care in forensic casework.