Nanobiopsy investigation of the subcellular mtDNA heteroplasmy in human tissues.

Mitochondrial function is critical to continued cellular vitality and is an important contributor to a growing number of human diseases. Mitochondrial dysfunction is typically heterogeneous, mediated through the clonal expansion of mitochondrial DNA (mtDNA) variants in a subset of cells in a given tissue. To date, our understanding of the dynamics of clonal expansion of mtDNA variants has been technically limited to the single cell-level. Here, we report the use of nanobiopsy for subcellular sampling from human tissues, combined with next-generation sequencing to assess subcellular mtDNA mutation load in human tissue from mitochondrial disease patients. The ability to map mitochondrial mutation loads within individual cells of diseased tissue samples will further our understanding of mitochondrial genetic diseases.

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.

Mitochondrial DNA: Inherent Complexities Relevant to Genetic Analyses.

Mitochondrial DNA (mtDNA) exhibits distinct characteristics distinguishing it from the nuclear genome, necessitating specific analytical methods in genetic studies. This comprehensive review explores the complex role of mtDNA in a variety of genetic studies, including genome-wide, epigenome-wide, and phenome-wide association studies, with a focus on its implications for human traits and diseases. Here, we discuss the structure and gene-encoding properties of mtDNA, along with the influence of environmental factors and epigenetic modifications on its function and variability. Particularly significant are the challenges posed by mtDNA's high mutation rate, heteroplasmy, and copy number variations, and their impact on disease susceptibility and population genetic analyses. The review also highlights recent advances in methodological approaches that enhance our understanding of mtDNA associations, advocating for refined genetic research techniques that accommodate its complexities. By providing a comprehensive overview of the intricacies of mtDNA, this paper underscores the need for an integrated approach to genetic studies that considers the unique properties of mitochondrial genetics. Our findings aim to inform future research and encourage the development of innovative methodologies to better interpret the broad implications of mtDNA in human health and disease.

Single-cell mtDNA dynamics in tumors is driven by coregulation of nuclear and mitochondrial genomes.

The extent of cell-to-cell variation in tumor mitochondrial DNA (mtDNA) copy number and genotype, and the phenotypic and evolutionary consequences of such variation, are poorly characterized. Here we use amplification-free single-cell whole-genome sequencing (Direct Library Prep (DLP+)) to simultaneously assay mtDNA copy number and nuclear DNA (nuDNA) in 72,275 single cells derived from immortalized cell lines, patient-derived xenografts and primary human tumors. Cells typically contained thousands of mtDNA copies, but variation in mtDNA copy number was extensive and strongly associated with cell size. Pervasive whole-genome doubling events in nuDNA associated with stoichiometrically balanced adaptations in mtDNA copy number, implying that mtDNA-to-nuDNA ratio, rather than mtDNA copy number itself, mediated downstream phenotypes. Finally, multimodal analysis of DLP+ and single-cell RNA sequencing identified both somatic loss-of-function and germline noncoding variants in mtDNA linked to heteroplasmy-dependent changes in mtDNA copy number and mitochondrial transcription, revealing phenotypic adaptations to disrupted nuclear/mitochondrial balance.

Single-cell analysis reveals context-dependent, cell-level selection of mtDNA.

Heteroplasmy occurs when wild-type and mutant mitochondrial DNA (mtDNA) molecules co-exist in single cells1. Heteroplasmy levels change dynamically in development, disease and ageing2,3, but it is unclear whether these shifts are caused by selection or drift, and whether they occur at the level of cells or intracellularly. Here we investigate heteroplasmy dynamics in dividing cells by combining precise mtDNA base editing (DdCBE)4 with a new method, SCI-LITE (single-cell combinatorial indexing leveraged to interrogate targeted expression), which tracks single-cell heteroplasmy with ultra-high throughput. We engineered cells to have synonymous or nonsynonymous complex I mtDNA mutations and found that cell populations in standard culture conditions purge nonsynonymous mtDNA variants, whereas synonymous variants are maintained. This suggests that selection dominates over simple drift in shaping population heteroplasmy. We simultaneously tracked single-cell mtDNA heteroplasmy and ancestry, and found that, although the population heteroplasmy shifts, the heteroplasmy of individual cell lineages remains stable, arguing that selection acts at the level of cell fitness in dividing cells. Using these insights, we show that we can force cells to accumulate high levels of truncating complex I mtDNA heteroplasmy by placing them in environments where loss of biochemical complex I activity has been reported to benefit cell fitness. We conclude that in dividing cells, a given nonsynonymous mtDNA heteroplasmy can be harmful, neutral or even beneficial to cell fitness, but that the 'sign' of the effect is wholly dependent on the environment.

AMPK Suppression Due to Obesity Drives Oocyte mtDNA Heteroplasmy via ATF5-POLG Axis.

Due to the exclusive maternal transmission, oocyte mitochondrial dysfunction reduces fertility rates, affects embryonic development, and programs offspring to metabolic diseases. However, mitochondrial DNA (mtDNA) are vulnerable to mutations during oocyte maturation, leading to mitochondrial nucleotide variations (mtSNVs) within a single oocyte, referring to mtDNA heteroplasmy. Obesity (OB) accounts for more than 40% of women at the reproductive age in the USA, but little is known about impacts of OB on mtSNVs in mature oocytes. It is found that OB reduces mtDNA content and increases mtSNVs in mature oocytes, which impairs mitochondrial energetic functions and oocyte quality. In mature oocytes, OB suppresses AMPK activity, aligned with an increased binding affinity of the ATF5-POLG protein complex to mutated mtDNA D-loop and protein-coding regions. Similarly, AMPK knockout increases the binding affinity of ATF5-POLG proteins to mutated mtDNA, leading to the replication of heteroplasmic mtDNA and impairing oocyte quality. Consistently, AMPK activation blocks the detrimental impacts of OB by preventing ATF5-POLG protein recruitment, improving oocyte maturation and mitochondrial energetics. Overall, the data uncover key features of AMPK activation in suppressing mtSNVs, and improving mitochondrial biogenesis and oocyte maturation in obese females.

Evidence for hybridization-driven heteroplasmy maintained across generations in a ricefish endemic to a Wallacean ancient lake.

Heteroplasmy, the presence of multiple mitochondrial DNA (mtDNA) haplotypes within cells of an individual, is caused by mutation or paternal leakage. However, heteroplasmy is usually resolved to homoplasmy within a few generations because of germ-line bottlenecks; therefore, instances of heteroplasmy are limited in nature. Here, we report heteroplasmy in the ricefish species Oryzias matanensis, endemic to Lake Matano, an ancient lake in Sulawesi Island, in which one individual was known to have many heterozygous sites in the mitochondrial NADH dehydrogenase subunit 2 (ND2) gene. In this study, we cloned the ND2 gene for some additional individuals with heterozygous sites and demonstrated that they are truly heteroplasmic. Phylogenetic analysis revealed that the extra haplotype within the heteroplasmic O. matanensis individuals clustered with haplotypes of O. marmoratus, a congeneric species inhabiting adjacent lakes. This indicated that the heteroplasmy originated from paternal leakage due to interspecific hybridization. The extra haplotype was unique and contained two non-synonymous substitutions. These findings demonstrate that this hybridization-driven heteroplasmy was maintained across generations for a long time to the extent that the extra mitochondria evolved within the new host.

Clinical Characteristics of Diabetes in People with Mitochondrial DNA 3243A>G Mutation in Korea.

Maternally inherited diabetes and deafness (MIDD) is a rare mitochondrial disorder primarily resulting from m.3243A>G mutation. The clinical characteristics of MIDD exhibit significant heterogeneity. Our study aims to delineate these characteristics and determine the potential correlation with m.3243A>G heteroplasmy levels. This retrospective, descriptive study encompassed patients with confirmed m.3243A>G mutation and diabetes mellitus at Seoul National University Hospital. Our cohort comprises 40 patients with MIDD, with a mean age at study enrollment of 33.3±12.9 years and an average % of heteroplasmy of 30.0%± 14.6% in the peripheral blood. The most prevalent comorbidity was hearing loss (90%), followed by albuminuria (61%), seizure (38%), and stroke (33%). We observed a significant negative correlation between % of heteroplasmy and age at diabetes diagnosis. These clinical features can aid in the suspicion of MIDD and further consideration of genetic testing for m.3243A>G mutation.

Mitochondrial point heteroplasmy: insights from deep-sequencing of human replicate samples.

BACKGROUND: Human mitochondrial heteroplasmy is an extensively investigated phenomenon in the context of medical diagnostics, forensic identification and molecular evolution. However, technical limitations of high-throughput sequencing hinder reliable determination of point heteroplasmies (PHPs) with minor allele frequencies (MAFs) within the noise threshold. RESULTS: To investigate the PHP landscape at an MAF threshold down to 0.1%, we sequenced whole mitochondrial genomes at approximately 7.700x coverage, in multiple technical and biological replicates of longitudinal blood and buccal swab samples from 11 human donors (159 libraries in total). The results obtained by two independent sequencing platforms and bioinformatics pipelines indicate distinctive PHP patterns below and above the 1% MAF cut-off. We found a high inter-individual prevalence of low-level PHPs (MAF < 1%) at polymorphic positions of the mitochondrial DNA control region (CR), their tissue preference, and a tissue-specific minor allele linkage. We also established the position-dependent potential of minor allele expansion in PHPs, and short-term PHP instability in a mitotically active tissue. We demonstrate that the increase in sensitivity of PHP detection to minor allele frequencies below 1% within a robust experimental and analytical pipeline, provides new information with potential applicative value. CONCLUSIONS: Our findings reliably show different mutational loads between tissues at sub-1% allele frequencies, which may serve as an informative medical biomarker of time-dependent, tissue-specific mutational burden, or help discriminate forensically relevant tissues in a single person, close maternal relatives or unrelated individuals of similar phylogenetic background.

Quantifying mitochondrial heteroplasmy diversity: A computational approach.

Biodiversity plays a pivotal role in sustaining ecosystem processes, encompassing diverse biological species, genetic types and the intricacies of ecosystem composition. However, the precise definition of biodiversity at the individual level remains a challenging endeavour. Hill numbers, derived from Rényi's entropy, have emerged as a popular measure of diversity, with a recent unified framework extending their application across various levels, from genetics to ecosystems. In this study, we employ a computational approach to exploring the diversity of mitochondrial heteroplasmy using real-world data. By adopting Hill numbers with q = 2, we demonstrate the feasibility of quantifying mitochondrial heteroplasmy diversity within and between individuals and populations. Furthermore, we investigate the alpha diversity of mitochondrial heteroplasmy among different species, revealing heterogeneity at multiple levels, including mitogenome components and protein-coding genes (PCGs). Our analysis explores large-scale mitochondrial heteroplasmy data in humans, examining the relationship between alpha diversity at the mitogenome components and PCGs level. Notably, we do not find a significant correlation between these two levels. Additionally, we observe significant correlations in alpha diversity between mothers and children in blood samples, exceeding the reported R2 value for allele frequency correlations. Moreover, our investigation of beta diversity and local overlay similarity demonstrates that heteroplasmy variant distributions in different tissues of children more closely resemble those of their mothers. Through systematic quantification and analysis of mitochondrial heteroplasmy diversity, this study enhances our understanding of heterogeneity at multiple levels, from individuals to populations, providing new insights into this fundamental dimension of biodiversity.

Whole mitogenome sequencing uncovers a relation between mitochondrial heteroplasmy and leprosy severity.

BACKGROUND: In recent years, the mitochondria/immune system interaction has been proposed, so that variants of mitochondrial genome and levels of heteroplasmy might deregulate important metabolic processes in fighting infections, such as leprosy. METHODS: We sequenced the whole mitochondrial genome to investigate variants and heteroplasmy levels, considering patients with different clinical forms of leprosy and household contacts. After sequencing, a specific pipeline was used for preparation and bioinformatics analysis to select heteroplasmic variants. RESULTS: We found 116 variants in at least two of the subtypes of the case group (Borderline Tuberculoid, Borderline Lepromatous, Lepromatous), suggesting a possible clinical significance to these variants. Notably, 15 variants were exclusively found in these three clinical forms, of which five variants stand out for being missense (m.3791T > C in MT-ND1, m.5317C > A in MT-ND2, m.8545G > A in MT-ATP8, m.9044T > C in MT-ATP6 and m.15837T > C in MT-CYB). In addition, we found 26 variants shared only by leprosy poles, of which two are characterized as missense (m.4248T > C in MT-ND1 and m.8027G > A in MT-CO2). CONCLUSION: We found a significant number of variants and heteroplasmy levels in the leprosy patients from our cohort, as well as six genes that may influence leprosy susceptibility, suggesting for the first time that the mitogenome might be involved with the leprosy process, distinction of clinical forms and severity. Thus, future studies are needed to help understand the genetic consequences of these variants.

Mitochondrial DNA Variants at Low-Level Heteroplasmy and Decreased Copy Numbers in Chronic Kidney Disease (CKD) Tissues with Kidney Cancer.

The human mitochondrial genome (mtDNA) is a circular DNA molecule with a length of 16.6 kb, which contains a total of 37 genes. Somatic mtDNA mutations accumulate with age and environmental exposure, and some types of mtDNA variants may play a role in carcinogenesis. Recent studies observed mtDNA variants not only in kidney tumors but also in adjacent kidney tissues, and mtDNA dysfunction results in kidney injury, including chronic kidney disease (CKD). To investigate whether a relationship exists between heteroplasmic mtDNA variants and kidney function, we performed ultra-deep sequencing (30,000×) based on long-range PCR of DNA from 77 non-tumor kidney tissues of kidney cancer patients with CKD (stages G1 to G5). In total, this analysis detected 697 single-nucleotide variants (SNVs) and 504 indels as heteroplasmic (0.5% ≤ variant allele frequency (VAF) < 95%), and the total number of detected SNVs/indels did not differ between CKD stages. However, the number of deleterious low-level heteroplasmic variants (pathogenic missense, nonsense, frameshift and tRNA) significantly increased with CKD progression (p < 0.01). In addition, mtDNA copy numbers (mtDNA-CNs) decreased with CKD progression (p < 0.001). This study demonstrates that mtDNA damage, which affects mitochondrial genes, may be involved in reductions in mitochondrial mass and associated with CKD progression and kidney dysfunction.

Vanadium Pentoxide Exposure Causes Strain-Dependent Changes in Mitochondrial DNA Heteroplasmy, Copy Number, and Lesions, but Not Nuclear DNA Lesions.

Interstitial lung diseases (ILDs) are lethal lung diseases characterized by pulmonary inflammation and progressive lung interstitial scarring. We previously developed a mouse model of ILD using vanadium pentoxide (V2O5) and identified several gene candidates on chromosome 4 associated with pulmonary fibrosis. While these data indicated a significant genetic contribution to ILD susceptibility, they did not include any potential associations and interactions with the mitochondrial genome that might influence disease risk. To conduct this pilot work, we selected the two divergent strains we previously categorized as V2O5-resistant C57BL6J (B6) and -responsive DBA/2J (D2) and compared their mitochondrial genome characteristics, including DNA variants, heteroplasmy, lesions, and copy numbers at 14- and 112-days post-exposure. While we did not find changes in the mitochondrial genome at 14 days post-exposure, at 112 days, we found that the responsive D2 strain exhibited significantly fewer mtDNA copies and more lesions than control animals. Alongside these findings, mtDNA heteroplasmy frequency decreased. These data suggest that mice previously shown to exhibit increased susceptibility to pulmonary fibrosis and inflammation sustain damage to the mitochondrial genome that is evident at 112 days post-V2O5 exposure.

Negative correlation between organ heteroplasmy, particularly hepatic heteroplasmy, and age at death revealed by post-mortem studies of m.3243A > G cases.

Mitochondrial DNA m.3243A > G mutation causes mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and its associated multi-organ disorders, including diabetes. To clarify associations between m.3243A > G organ heteroplasmy and clinical phenotypes, including the age at death, we combined genetic and pathological examinations from seven unreported and 36 literature cases of autopsied subjects. Clinical characteristics of subjects were as follows: male, 13; female, 28; unknown, 2; the age at death, 36.9 ± 20.2 [4-82] years; BMI, 16.0 ± 2.9 [13.0-22.3]; diabetes, N = 21 (49%), diabetes onset age 38.6 ± 14.2 years; deafness, N = 27 (63%); stroke-like episodes (StLEp), N = 25 (58%); congestive heart failure (CHF), N = 15 (35%); CHF onset age, 51.3 ± 14.5 years. Causes of death (N = 32) were as follows: cardiac, N = 13 (41%); infection, N = 8 (25%); StLEp, N = 4 (13%); gastrointestinal, N = 4 (13%); renal, N = 2 (6%); hepatic, N = 1 (2%). High and low heteroplasmies were confirmed in non-regenerative and regenerative organs, respectively. Heteroplasmy of the liver, spleen, leukocytes, and kidney for all subjects was significantly associated with the age at death. Furthermore, the age at death was related to juvenile-onset (any m.3243A > G-related symptoms appeared before 20) and stroke-like episodes. Multiple linear regression analysis with the age at death as an objective variable showed the significant contribution of liver heteroplasty and juvenile-onset to the age at death. m.3243A > G organ heteroplasmy levels, particularly hepatic heteroplasmy, are significantly associated with the age at death in deceased cases.

Deleterious heteroplasmic mitochondrial mutations are associated with an increased risk of overall and cancer-specific mortality.

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia.

PHFinder: assisted detection of point heteroplasmy in Sanger sequencing chromatograms.

Heteroplasmy is the presence of two or more organellar genomes (mitochondrial or plastid DNA) in an organism, tissue, cell or organelle. Heteroplasmy can be detected by visual inspection of Sanger sequencing chromatograms, where it appears as multiple peaks of fluorescence at a single nucleotide position. Visual inspection of chromatograms is both consuming and highly subjective, as heteroplasmy is difficult to differentiate from background noise. Few software solutions are available to automate the detection of point heteroplasmies, and those that are available are typically proprietary, lack customization or are unsuitable for automated heteroplasmy assessment in large datasets. Here, we present PHFinder, a Python-based, open-source tool to assist in the detection of point heteroplasmies in large numbers of Sanger chromatograms. PHFinder automatically identifies point heteroplasmies directly from the chromatogram trace data. The program was tested with Sanger sequencing data from 100 humpback whales (Megaptera novaeangliae) tissue samples with known heteroplasmies. PHFinder detected most (90%) of the known heteroplasmies thereby greatly reducing the amount of visual inspection required. PHFinder is flexible and enables explicit specification of key parameters to infer double peaks (i.e., heteroplasmies).

CmVCall: An automated and adjustable nanopore analysis pipeline for heteroplasmy detection of the control region in human mitochondrial genome.

Genetic associations between human mitochondrial DNA (mtDNA) heteroplasmy and mitochondrial diseases, aging, and cancer have been elaborated, contributing a lot to the further understanding of mtDNA polymorphic spectrum in anthropology, population, and forensic genetics. In the past decade, heteroplasmy detection using Sanger sequencing and next generation sequencing (NGS) was hampered by the former's inefficiency and the latter's inherent bias due to amplification and mapping of short reads, respectively. Nanopore sequencing stands out for its ability to yield long contiguous segments of DNA, providing a new insight into heterogeneity authentication. In addition to MinION from Oxford Nanopore Technologies, an alternative nanopore sequencer QNome (Qitan Technology) has also been applied to various biological research and the forensic applicability of this platform has been proved recently. In this study, we evaluated the performance of four commonly used variant callers in the heterogeneity authentication of the control region of human mtDNA based on simulations of different ratios generated by mixing QNome nanopore sequencing reads of two synthetic sequences. Then, an open-source and python-based nanopore analytics pipeline, CmVCall was developed and incorporated multiple programs including reads filtering, removal of nuclear mitochondrial sequences (NUMTs), alignment, optional 'Correction' mode, and heterogeneity identification. CmVCall can achieve high precision, accuracy, and recall of 100%, 99.9%, and 92.3% with a 5% heteroplasmy level in 'Correction' mode. Moreover, blood, saliva, and hair shaft samples from monozygotic (MZ) twins were used for heterogeneity evaluation and comparison with the NGS data. Results of MZ twin samples showed that CmVCall could identify more point heteroplasmy sites, revealing significant levels of inter- and intra-individual mtDNA polymorphism. In conclusion, we believe that this analysis pipeline will lay a solid foundation for the development of a comprehensive nanopore analysis pipeline targeting the whole mitochondrial genome.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans.

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

Genetic diversity of thrips populations on Allium species around the world.

Thrips are a serious pest in many crops. In onion cultivation, Thrips tabaci is the most important, but not the only thrips species causing damage. We investigated which thrips species affects onion and related species worldwide, how much genetic variation there is within T. tabaci populations, and how this evolves. Furthermore, we determined the reproductive mode and the correlation between the genetic and geographic distances. Thrips samples from infested onions or related species were obtained from 14 different locations worldwide. Species and haplotypes were determined through DNA barcoding with the mitochondrial Cytochrome Oxidase subunit I (COI) gene. Thrips tabaci was the most commonly observed species, but Scirtothrips dorsalis, Thrips palmi, Frankliniella intonsa, Frankliniella occidentalis and Frankliniella tenuicornis were also found, especially at the beginning of the growing seasons and depending on the location. The Nei's genetic distance within T. tabaci was less than 5% and the haplotypes were clustered into two phylogenetic groups, each linked to a specific mode of reproduction, thelytokous or arrhenotokous. Thelytokous thrips were more common and more widely distributed than arrhenotokous thrips. A high percentage of heteroplasmy was detected in the arrhenotokous group. Heteroplasmic thrips were only found in populations where thelytokous and arrhenotokous were present in sympatry. Some T. tabaci haplotypes were present in high frequency at several sampled locations. No correlation was found between the genetic and geographic distances, which points to anthropic activities spreading thrips haplotypes throughout the world.