Differential diagnosis between Tunga penetrans (L., 1758) and T. trimamillata Pampiglione et al., 2002 (Insecta, Siphonaptera), the two species of the genus Tunga parasitic in man

Of the ten currently known species of sand fleas, only two, Tunga penetrans and Tunga trimamillata, are known to be parasites of man, besides other warm blooded animals, most of which are peridomestic. The hosts of the other eight are limited to a few genera of wild mammals. T. trimamillata was only recently identified and differentiated from T. penetrans by features of the gravid female phase. In the present paper the different morphological characters of both for non-gravid females and males of the two species are described. In the non-gravid-females the distinguishing characters are: a) differences in the length of the first two segments of the maxillary palps (the first is the longest in T. trimamillata, whereas the second is the longest in T. penetrans), this feature is associated with another character i.e. the presence of short, thick spines in addition to the bristles on the surface of the segments only in T. trimamillata; b) the presence, only in T. trimamillata, of a row of spines on the antero-medial surface of the libia of the 3rd pair of legs; c) the last abdominal spiracle protrudes in T. trimamillata but not in T. penetrans; d) the hood of the hilla in spermatheca is surrounded by a papilla only in T. penetrans. The following morphological characters differentiate males of T. trimamillata and T. penetrans: a) T. trimamillata has a row of spines on the antero-medial surface of the tibia of the 3rd pair of legs; b) the diameter of the abdominal spiracles of T. trimamillata is smaller than that of T. penetrans and the edges of the spiracles are more regular; c) the claspers and aedeagus of T. trimamillata are shorter than those of T. penetrans; d) the profile of the oedeagal apodeme differs between the two species.     

Wolbachia infection in the newly described Ecuadorian sand flea, Tunga trimamillata

Wolbachia pipientis is an intracellular endosymbiont producing reproductive alterations in its hosts. This bacterium have been reported in many arthropods and nematodes. By PCR amplification and sequencing of the 16S rDNA and ftsZ genes we have identified a Wolbachia strain in the newly described sand-flea, Tunga trimamillata. Prevalence of this endosymbiont in the 26 individuals screened is equal to 35%. Sympatric and allopatric specimens of the related species Tunga penetrans were also analysed, but in contrast to literature data, Wolbachia appears absent in the presently analysed 24 specimens. Field studies evidence a female-biased sex-ratio in T. trimamillata, suggesting that Wolbachia may cause sex-ratio distortion in this species. By means of BLAST search and phylogenetic analysis we found that the Wolbachia strain from T. trimamillata pertains to the arthropod-infecting Wolbachia; this strain is highly differentiated from the Wolbachia strain of T. penetrans described in literature.     

Phylogenetic relationships of phrynosomatid lizards based on nuclear and mitochondrial data,and a revised phylogeny for Sceloporus

Phrynosomatid lizards are among the most common and diverse groups of reptiles in western North America, Mexico, and Central America. Phrynosomatidae includes 136 species in 10 genera. Phrynosomatids are used as model systems in many research programs in evolution and ecology, and much of this research has been undertaken in a comparative phylogenetic framework. However, relationships among many phrynosomatid genera are poorly supported and in conflict between recent studies. Further, previous studies based on mitochondrial DNA sequences suggested that the most species-rich genus (Sceloporus) is possibly paraphyletic with respect to as many as four other genera (Petrosaurus, Sator, Urosaurus, and Uta). Here, we collect new sequence data from five nuclear genes and combine them with published data from one additional nuclear gene and five mitochondrial gene regions. We compare trees from nuclear and mitochondrial data from 37 phrynosomatid taxa, including a “species tree” (from BEST) for the nuclear data. We also present a phylogeny for 122 phrynosomatid species based on maximum likelihood analysis of the combined data, which provides a strongly-supported hypothesis for relationships among most phrynosomatid genera and includes most phrynosomatid species. Our results strongly support the monophyly of Sceloporus (including Sator) and many of the relationships within it. We present a new classification for phrynosomatid lizards and the genus Sceloporus, and offer a new tree with branch lengths for use in comparative studies.     

Congruent mammalian trees from mitochondrial and nuclear genes using Bayesian methods

Analyses of mitochondrial and nuclear gene sequences have often produced different mammalian tree topologies, undermining confidence in the merit of molecular approaches with respect to "traditional" morphological classification. The recent sequencing of the complete mitochondrial genomes of two additional rodents (Spalax judaei and Jaculus jaculus) and one lagomorph (Ochotona princeps) has prompted us to reinvestigate the issue. Using Bayesian phylogenetics, we found phylogenetic relationships between mammalian species highly congruent with previous results based on nuclear genes. Our results show the existence of four primary lineages of placental mammals: Xenarthra, Afrotheria, Laurasiatheria, and Euarchontoglires. Relationships between and within these lineages strongly suggest that the gene trees may also be congruent with the underlying species phylogeny.     

Phylogeography of the pine processionary moth Thaumetopoea wilkinsoni in the Near East

Phylogeographic structure of the eastern pine processionary moth Thaumetopoea wilkinsoni was explored in this study by means of nested clade phylogeographic analyses of COI and COII sequences of mitochondrial DNA and Bayesian estimates of divergence times. Intraspecific relationships were inferred and hypotheses tested to understand historical spread patterns and spatial distribution of genetic variation. Analyses revealed that all T. wilkinsoni sequences were structured in three clades, which were associated with two major biogeographic events, the colonization of the island of Cyprus and the separation of southwestern and southeastern Anatolia during the Pleistocene. Genetic variation in populations of T. wilkinsoni was also investigated using amplified fragment length polymorphisms and four microsatellite loci. Contrasting nuclear with mitochondrial data revealed recurrent gene flow between Cyprus and the mainland, related to the long-distance male dispersal. In addition, a reduction in genetic variability was observed at both mitochondrial and nuclear markers at the expanding boundary of the range, consistent with a recent origin of these populations, founded by few individuals expanding from nearby localities. In contrast, several populations fixed for one single mitochondrial haplotype showed no reduction in nuclear variability, a pattern that can be explained by recurrent male gene flow or selective sweeps at the mitochondrial level. The use of both mitochondrial and nuclear markers was essential in understanding the spread patterns and the population genetic structure of T. wilkinsoni, and is recommended to study colonizing species characterized by sex-biased dispersal.     

TESTS OF PHYLOGEOGRAPHIC MODELS WITH NUCLEAR AND MITOCHONDRIAL DNA SEQUENCE VARIATION IN THE STONE CRABS,MENIPPE ADINA AND MENIPPE MERCENARIA

Evolutionary relationships among stone crabs (Menippe) from the Gulf of Mexico and western Atlantic were investigated by comparisons of restriction sites within anonymous nuclear DNA sequences and nucleotide sequences of both mitochondrial and a duplicated nuclear form of the mitochondrial large subunit ribosomal RNA (LSrDNA) gene. A survey of over 100 restriction sites by Southern blot analysis with 10 anonymous nuclear DNA sequence probes failed to reveal any differences between Menippe adina and M. mercenaria. Sequence comparisons of both mitochondrial and nuclear forms of the LSrDNA gene also did not distinguish these species. Although both LSrDNA gene sequences were variable, some haplotypes were shared by the two species, implying either incomplete gene lineage sorting or introgressive hybridization. Based on molecular clock calibrations, we estimate that all of the observed mitochondrial LSrDNA sequences share a common ancestor between 1.5 and 2.7 million years before present (M.Y.B.P.). However, because identical sequences are shared by the two species, these data are also compatible with a more recent common ancestry. These findings conflict with a previously proposed biogeographic scenario for North American Menippe, which featured a relict hybrid zone on the Atlantic Coast. We suggest an alternative scenario based on relatively recent events and ongoing, rather than historical, gene flow.     

Population genetic patterns revealed by microsatellite data challenge the mitochondrial DNA based taxonomy of Astyanax in Mexico (Characidae, Teleostei)

Astyanax has become an important model system for evolutionary studies of cave animals. We investigated correlations of population genetic patterns revealed by microsatellite data and phylogeographic patterns shown by mitochondrial DNA sequences in Mexican cave and surface fish of the genus Astyanax (Characidae, Teleostei) to improve the understanding of the colonization history of this neotropical fish in Central and North America and to assess a recent taxonomic classification. The distribution of nuclear genotypes is not congruent with that of the mitochondrial clades. Admixture analyses suggest there has been nuclear gene flow between populations defined by different mitochondrial clades. The microsatellite data indicate that there was mitochondrial capture of a cave population from adjacent populations. Furthermore, gene flow also occurred between populations belonging to different nuclear genotypic clusters. This indicates that neither the nuclear genotypic clusters nor the mitochondrial clades represent independent evolutionary units, although the mitochondrial divergences are high and in a range usually characteristic for different fish species. This conclusion is supported by the presence of morphologically intermediate forms. Our analyses show that the Trans-Mexican Volcanic Belt limited gene flow, but has been crossed by Astyanax several times. In Yucatán, where obvious geographic barriers are missing, the incongruence between the distribution of nuclear and mitochondrial markers reflects random colonization events caused by inundations or marine transgressions resulting in random phylogeographic breaks. Thus, conclusions about the phylogeographic history and even more about the delimitation of species should not be based on single genetic markers.     

A comprehensive analysis of mammalian mitochondrial genome base composition and improved phylogenetic methods

Phylogenetic analysis of mammalian species using mitochondrial protein genes has proved to be problematic in many previous studies. The high mutation rate of mitochondrial DNA and unusual base composition of several species has prompted us to conduct a detailed study of the composition of 69 mammalian mitochondrial genomes. Most major changes in base composition between lineages can be attributed to shifts between the proportions of C and T on the L-strand. These changes are significant at all codon positions and are shown to affect amino acid composition. Correlated changes in the base composition of the RNA loops and stems are also observed. Following up from previous studies, we investigate changes in the base composition of all 12 H-strand proteins and find that variability in proportions of C and T is correlated with location on the genome. Variation in base composition across genes and species is known to adversely affect the performance of phylogenetic inference methods. We have, therefore, developed a customized three-state general time-reversible DNA substitution model, implemented in the PHASE phylogenetic inference package, which lumps C and T into a composite pyrimidine state. We compare the phylogenetic tree obtained using the new three-state model with that obtained using a standard four-state model. Results using the three-state model are more congruent with recent studies using large sets of nuclear genes and help resolve some of the apparent conflicts between studies using nuclear and mitochondrial proteins.     

Speciation, introgressive hybridization and nonlinear rate of molecular evolution in flycatchers

Evolutionary history of Muscicapidae flycatchers is inferred from nuclear and mitochondrial DNA (mtDNA) sequence comparisons and population genetic analysis of nuclear and mtDNA markers. Phylogenetic reconstruction based on sequences from the two genomes yielded similar trees with respect to the order at which the species split off. However, the genetic distances fitted a nonlinear, polynomial model reflecting diminishing divergence rate of the mtDNA sequences compared to the nuclear DNA sequences. This could be explained by Haldane's rule because genetic isolation might evolve more rapidly on the mitochondrial rather than the nuclear genome in birds. This is because hybrid sterility of the heterogametic sex (females) would predate that of the homogametic sex (males), leading to sex biased introgression of nuclear genes. Analyses of present hybrid zones of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) may indicate a slight sexual bias in rate of introgression, but the introgression rates were too low to allow proper statistical analyses. It is suggested, however, that the observed deviation from linearity can be explained by a more rapid mutational saturation of the mtDNA sequences than of the nuclear DNA sequences, as supported by analyses of third codon position transversions at two protein coding mtDNA genes. A phylogeographic scenario for the black and white flycatcher species is suggested based on interpretation of the genetic data obtained. Four species appear to have diverged from a common ancestor relatively simultaneously during the Pleistocene. After the last glaciation period, pied and collared flycatchers expanded their breeding ranges and eventually came into secondary contact in Central and Eastern Europe and on the Baltic Isles.     

Population Structure of the Giant Garter Snake, Thamnophis gigas

The giant garter snake, Thamnophis gigas, is a threatened species endemic to California’s Central Valley. We tested the hypothesis that current watershed boundaries have caused genetic differentiation among populations of T. gigas. We sampled 14 populations throughout the current geographic range of T. gigas and amplified 859 bp from the mitochondrial gene ND4 and one nuclear microsatellite locus. DNA sequence variation from the mitochondrial gene indicates there is some genetic structuring of the populations, with high F_ST values and unique haplotypes occurring at high frequency in several populations. We found that clustering populations by watershed boundary results in significant between-region genetic variance for mtDNA. However, analysis of allele frequencies at the microsatellite locus NSU3 reveals very low F_ST values and little between-region variation in allele frequencies. The discordance found between mitochondrial and microsatellite data may be explained by aspects of molecular evolution and/or T. gigas life history characteristics. Differences in effective population size between mitochondrial and nuclear DNA, or male-biased gene flow, result in a lower migration rate of mitochondrial haplotypes relative to nuclear alleles. However, we cannot exclude homoplasy as one explanation for homogeneity found for the single microsatellite locus. The mitochondrial nucleotide sequence data supports conservation practices that identify separate management units for T. gigas.     

Morphology,molecules and taxonomy: extreme incongruence in pleurocerids (Gastropoda,Cerithioidea, Pleuroceridae)

Pronounced mitochondrial heterogeneity within putative species of Pleuroceridae has prevented meaningful systematic revisions of this critically imperilled freshwater family. Previous studies have demonstrated that this mitochondrial diversity often produces polyphyletic species on mitochondrial gene trees, but its significance is unclear. Hypotheses advanced to explain this pattern have included cryptic species, retained ancestral polymorphisms and introgression; other possible explanations such as doubly uniparental inheritance or the presence of pseudogenes have not been given due consideration. Previous analyses have not included adequate sampling, neither in terms of number of individuals nor in geographic coverage, to adequately test any of these hypotheses. To fully characterize mitochondrial heterogeneity in pleurocerids and robustly assess possible causal explanations, we collected 239 individuals representing four putative species from seven sites and sequenced the COI and 16S rRNA mitochondrial genes and the H3 nuclear gene for all individuals. We also used whole‐genome shotgun sequencing to construct and annotate a mitochondrial genome for one individual. Characters with demonstrated utility in morphospecies delineation of gastropods (head–foot coloration, shell and radular morphology, pallial oviduct anatomy) were examined for a subset of individuals to determine whether morphology co‐varied by haplotype clade. We found pronounced mitochondrial heterogeneity at both the population level and species level in three of the species examined, but our data reject paralogous nuclear copies of mitochondrial genes (NUMTs) and doubly uniparental inheritance as causal mechanisms, and there was no evidence of cryptic diversity. Mutation rates were found to differ significantly among mitochondrial lineages, and population genetic statistics revealed a signature of balancing selection that could be acting to maintain this diversity. The observed pattern is similar to that seen in lineages with inherited endosymbionts like Wolbachia infections, which merits further investigation. Although questions remain concerning the precise cause(s) of intraspecific mitochondrial diversity in pleurocerids, nuclear and/or genomic data, combined with anatomical and life history investigations in an integrative phylogenetic context, is the most promising avenue for resolving pleurocerid systematics.     

Mouse models of mitochondrial disease, oxidative stress, and senescence.

During the course of normal respiration, reactive oxygen species are produced which are particularly detrimental to mitochondrial function. This is shown by recent studies with a mouse that lacks the mitochondrial form of superoxide dismutase (Sod2). Tissues that are heavily dependent on mitochondrial function such as the brain and heart are most severely affected in the Sod2 mutant mouse. Recent work with a mouse mutant for the heart/muscle specific isoform of the mitochondrial adenine nuclear translocator (Ant1) demonstrates a potential link between mitochondrial oxidative stress and mitochondrial DNA mutations. These mutations can be detected by Long-extension PCR, a method for detecting a wide variety of mutations of the mitochondrial genome. Such mutations have also been observed in the mitochondrial genome with senescence regardless of the mean or maximal lifespan of the organism being studied. Mutations have been detected with age in Caenorhabditis elegans, mice, chimpanzees, and humans. This implies that a causal relationship may exist between mitochondrial reactive oxygen species production, and the senescence specific occurrence of mitochondrial DNA mutations.     

Ancient hybridization and mitochondrial capture between two species of chipmunks

Models that posit speciation in the face of gene flow are replacing classical views that hybridization is rare between animal species. We use a multilocus approach to examine the history of hybridization and gene flow between two species of chipmunks ( Tamias ruficaudus and T. amoenus ). Previous studies have shown that these species occupy different ecological niches and have distinct genital bone morphologies, yet appear to be incompletely isolated reproductively in multiple areas of sympatry. We compared data from four sequenced nuclear loci and from seven microsatellite loci to published cytochrome b sequences. Interspecific gene flow was primarily restricted to introgression of the T. ruficaudus mitochondrial genome into a sympatric subspecies of T. amoenus , T. a. canicaudus , with the four sequenced nuclear loci showing little to no interspecific allele sharing. Microsatellite data were consistent with high levels of differentiation between the species and also showed no current gene flow between broadly sympatric populations of T. a. canicaudus and T. ruficaudus . Coalescent analyses date the mtDNA introgression event from the mid-Pleistocene to late Pliocene. Overall, these data indicate that introgression has had a minimal impact on the nuclear genomes of T. amoenus and T. ruficaudus despite multiple independent hybridization events. Our findings challenge long-standing assumptions on patterns of reproductive isolation in chipmunks and suggest that there may be other examples of hybridization among the 23 species of Tamias that occur in western North America.     

Bridging the gap: western rock skinks (Trachylepis sulcata) have a short history in South Africa

Phylogeographic patterns in wide-ranging species in southern Africa remain largely unexplored, especially in areas north of South Africa. Here, we investigate population structuring, demographic history, and the colonization pattern of the western rock skink (Trachylepis sulcata), a rock-dwelling species with a range extending from southwestern South Africa into Angola. Using 1056 bp from the mitochondrial marker ND2 and > 2.5 kb from three nuclear genes (EXPH5, KIF24, RAG-1), we constructed allele networks, generated extended Bayesian skyline plots and performed population clustering analyses. Analyses of historical demographic patterns show an overall southward range expansion from Northern Namibia into Southern Namibia and South Africa, although we find contrasting genetic breaks across these geographic regions using nuclear and mitochondrial data. We suggest that mtDNA has introgressed across a nuclear break corresponding to the Knersvlakte region of South Africa, a previously proposed biogeographic barrier for rupicolous species. This pattern of mitochondrial variation contrasts sharply to that of other South African taxa previously investigated, which all show significant mtDNA differentiation across the Knersvlakte region. Additionally, while other taxa show divergences dating to the Pliocene, T. sulcata appears to be a recent arrival in southern Africa, having crossed this barrier and colonized South Africa in the mid-Pleistocene. The complex phylogeographic history of T. sulcata corroborates the intricate patterns of genetic variation found in South African taxa and provides novel insight into historical processes affecting species distributed across Namibia.     

The Influence of Pleistocene Climatic Changes and Ocean Currents on the Phylogeography of the Southern African Barnacle,Tetraclita serrata (Thoracica; Cirripedia)

The evolutionary effects of glacial periods are poorly understood for Southern Hemisphere marine intertidal species, particularly obligatory sessile organisms. We examined this by assessing the phylogeographic patterns of the southern African volcano barnacle, Tetraclita serrata, a dominant species on rocky intertidal shores. Restricted gene flow in some geographical areas was hypothesized based on oceanic circulation patterns and known biogeographic regions. Barnacle population genetic structure was investigated using the mitochondrial cytochrome oxidase subunit 1 (COI) region for 410 individuals sampled from 20 localities spanning the South African coast. The mtDNA data were augmented by generating nuclear internal transcribed spacer 1 (ITS1) sequences from a subset of samples. Phylogenetic and population genetic analyses of mitochondrial DNA data reveal two distinct clades with mostly sympatric distributions, whereas nuclear analyses reveal only a single lineage. Shallow, but significant structure (0.0041–0.0065, P<0.01) was detected for the mtDNA data set, with the south-west African region identified as harbouring the highest levels of genetic diversity. Gene flow analyses on the mtDNA data show that individuals sampled in south-western localities experience gene flow primarily in the direction of the Benguela Current, while south and eastern localities experience bi-directional gene flow, suggesting an influence of both the inshore currents and the offshore Agulhas Current in the larval distribution of T. serrata. The mtDNA haplotype network, Bayesian Skyline Plots, mismatch distributions and time since expansion indicate that T. serrata population numbers were not severely affected by the Last Glacial Maximum (LGM), unlike other southern African marine species. The processes resulting in the two morphologically cryptic mtDNA lineages may be the result of a recent historical allopatric event followed by secondary contact or could reflect selective pressures due to differing environmental conditions.     

Multilocus phylogeny and systematics of Iberian endemic Squalius (Actinopterygii,Leuciscidae)

Inferring the evolutionary history of a group of species can be challenging given the many factors involved. In recent years, the increased availability of sequences of multiple genes per species has spurred the development of new methodologies to analyse multilocus data sets. Two approaches that analyse such data are concatenated supermatrix and coalescent-based species-tree analyses. In this study, we used both of these methods to infer the phylogenetic relationships of Iberian species of the genus Squalius from one mitochondrial and six nuclear genes. We found mitonuclear discordance in the phylogenetic relationships of the group. According to the mitochondrial gene analysis, all species were recovered as monophyletic except S. pyrenaicus; besides, in the concatenated supermatrix analysis of the nuclear markers, this species resolved as polyphyletic with three divergent evolutionary lineages. The coalescent-based nuclear species-tree analysis rendered a well-resolved phylogeny compared with the supermatrix analysis, which was unable to discern between S. carolitertii, S. castellanus and one of the evolutionary lineages of S. pyrenaicus. This result is likely due to the better integration of population uncertainty in the coalescent approach. Furthermore, Bayesian multilocus species delimitation analyses based on a BPP approach strongly supported the distinct nuclear lineages as different species. Nevertheless, the supermatrix analysis was able to obtain well-supported relationships in the divergent lineages with low numbers of individuals. Our study highlights the usefulness of different analytical methodologies to obtain a more complete picture of the evolutionary history of taxa, especially when discordant patterns among genes are found.     

The on‐again,off‐again relationship between mitochondrial genomes and species boundaries

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co‐evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.     

Gene Expression Pattern in Transmitochondrial Cytoplasmic Hybrid Cells Harboring Type 2 Diabetes-Associated Mitochondrial DNA Haplogroups

Decreased mitochondrial function plays a pivotal role in the pathogenesis of type 2 diabetes mellitus (T2DM). Recently, it was reported that mitochondrial DNA (mtDNA) haplogroups confer genetic susceptibility to T2DM in Koreans and Japanese. Particularly, mtDNA haplogroup N9a is associated with a decreased risk of T2DM, whereas haplogroups D5 and F are associated with an increased risk. To examine functional consequences of these haplogroups without being confounded by the heterogeneous nuclear genomic backgrounds of different subjects, we constructed transmitochondrial cytoplasmic hybrid (cybrid) cells harboring each of the three haplogroups (N9a, D5, and F) in a background of a shared nuclear genome. We compared the functional consequences of the three haplogroups using cell-based assays and gene expression microarrays. Cell-based assays did not detect differences in mitochondrial functions among the haplogroups in terms of ATP generation, reactive oxygen species production, mitochondrial membrane potential, and cellular dehydrogenase activity. However, differential expression and clustering analyses of microarray data revealed that the three haplogroups exhibit a distinctive nuclear gene expression pattern that correlates with their susceptibility to T2DM. Pathway analysis of microarray data identified several differentially regulated metabolic pathways. Notably, compared to the T2DM-resistant haplogroup N9a, the T2DM-susceptible haplogroup F showed down-regulation of oxidative phosphorylation and up-regulation of glycolysis. These results suggest that variations in mtDNA can affect the expression of nuclear genes regulating mitochondrial functions or cellular energetics. Given that impaired mitochondrial function caused by T2DM-associated mtDNA haplogroups is compensated by the nuclear genome, we speculate that defective nuclear compensation, under certain circumstances, might lead to the development of T2DM.