Genetic consequences of postglacial range expansion in two codistributed rodents (genus Dipodomys) depend on ecology and genetic locus

How does range expansion affect genetic diversity in species with different ecologies, and do different types of genetic markers lead to different conclusions? We addressed these questions by assessing the genetic consequences of postglacial range expansion using mitochondrial DNA (mtDNA) and nuclear restriction site‐associated DNA (RAD) sequencing in two congeneric and codistributed rodents with different ecological characteristics: the desert kangaroo rat (Dipodomys deserti), a sand specialist, and the Merriam's kangaroo rat (Dipodomys merriami), a substrate generalist. For each species, we compared genetic variation between populations that retained stable distributions throughout glacial periods and those inferred to have expanded since the last glacial maximum. Our results suggest that expanded populations of both species experienced a loss of private mtDNA haplotypes and differentiation among populations, as well as a loss of nuclear single‐nucleotide polymorphism (SNP) private alleles and polymorphic loci. However, only D. deserti experienced a loss of nucleotide diversity (both mtDNA and nuclear) and nuclear heterozygosity. For all indices of diversity and differentiation that showed reduced values in the expanded areas, D. deserti populations experienced a greater degree of loss than did D. merriami populations. Additionally, patterns of loss in genetic diversity in expanded populations were substantially less extreme (by two orders of magnitude in some cases) for nuclear SNPs in both species compared to that observed for mitochondrial data. Our results demonstrate that ecological characteristics may play a role in determining genetic variation associated with range expansions, yet mtDNA diversity loss is not necessarily accompanied by a matched magnitude of loss in nuclear diversity.     

Nuclear and mitochondrial DNA reveal contrasting evolutionary processes in populations of deer mice (Peromyscus maniculatus)

We investigated a major geographic break in mitochondrial DNA (mtDNA) haplotypes in deer mice, Peromyscus maniculatus, by analysing spatial variation in a 491‐bp fragment of the mtDNA control region from 455 samples distributed across a north–south transect of 2000 km in Western North America. To determine whether the mtDNA break was reflected in the nuclear genome, we then compared spatial variation in 13 nuclear microsatellites of 95 individuals surrounding the mtDNA break. Using a canonical correlation analysis we found that nuclear genomic variation was not correlated with mtDNA differentiation. The contrasting patterns of variation in mtDNA and nuclear DNA are consistent with a hypothesis of historic genetic drift that occurred in isolated refugia combined with recent gene flow between the formerly isolated refugial populations. A Mantel test of genetic vs. geographic distance revealed that recent gene flow between deer mouse populations has been high. We conclude that past vicariant events associated with Pleistocene climate changes together with recent gene flow have created the observed intra‐specific cytonuclear discordance in Western North America.     

Viability of cytochrome c genotypes depends on cytoplasmic backgrounds in Tigriopus californicus

Because of their extensive functional interaction, mitochondrial DNA (mtDNA) and nuclear genes may evolve to form coadapted complexes within reproductively isolated populations. As a consequence of coadaptation, the fitness of particular nuclear alleles may depend on mtDNA genotype. Among populations of the copepod Tigriopus californicus, there are high levels of amino acid substitutions in both the mtDNA genes encoding subunits of cytochrome c oxidase (COX) and the nuclear gene encoding cytochrome c (CYC), the substrate for COX. Because of the functional interaction between enzyme and substrate proteins, we hypothesized that the fitness of CYC genotypes would depend on mtDNA genotype. To test this hypothesis, segregation ratios for CYC and a second nuclear marker (histone H1) unrelated to mitochondrial function were scored in F2 progeny of several reciprocal interpopulation crosses. Genotypic ratios at the CYC locus (but not the H1 locus) differed between reciprocal crosses and differed from expected Mendelian ratios, suggesting that CYC genotypic fitnesses were strongly influenced by cytoplasmic (including mtDNA) background. However, in most cases the nature of the deviations from Mendelian ratios and differences between reciprocal crosses are not consistent with simple coevolution between CYC and mtDNA background. In a cross in which both newly hatched larvae and adults were sampled, only the adult sample showed deviations from Mendelian ratios, indicating that genotypic viabilities differed. In two of six crosses, large genotypic ratio differences for CYC were observed between the sexes. These results suggest that significant variation in nuclear-mtDNA coadaptation may exist between T. californicus populations and that the relative viability of specific cytonuclear allelic combinations is somehow affected by sex.     

Mitochondrial regulation of epigenetics and its role in human diseases

Most pathogenic mitochondrial DNA (mtDNA) mutations induce defects in mitochondrial oxidative phosphorylation (OXPHOS). However, phenotypic effects of these mutations show a large degree of variation depending on the tissue affected. These differences are difficult to reconcile with OXPHOS as the sole pathogenic factor suggesting that additional mechanisms contribute to lack of genotype and clinical phenotype correlationship. An increasing number of studies have identified a possible effect on the epigenetic landscape of the nuclear genome as a consequence of mitochondrial dysfunction. In particular, these studies demonstrate reversible or irreversible changes in genomic DNA methylation profiles of the nuclear genome. Here we review how mitochondria damage checkpoint (mitocheckpoint) induces epigenetic changes in the nucleus. Persistent pathogenic mutations in mtDNA may also lead to epigenetic changes causing genomic instability in the nuclear genome. We propose that “mitocheckpoint” mediated epigenetic and genetic changes may play key roles in phenotypic variation related to mitochondrial diseases or host of human diseases in which mitochondrial defect plays a primary role.     

Mito‐nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW–NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW–NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito‐nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within‐species discord. Male‐mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.     

Mitochondrial DNA under siege in avian phylogeography

Mitochondrial DNA (mtDNA) has been the workhorse of research in phylogeography for almost two decades. However, concerns with basing evolutionary interpretations on mtDNA results alone have been voiced since the inception of such studies. Recently, some authors have suggested that the potential problems with mtDNA are so great that inferences about population structure and species limits are unwarranted unless corroborated by other evidence, usually in the form of nuclear gene data. Here we review the relative merits of mitochondrial and nuclear phylogeographical studies, using birds as an exemplar class of organisms. A review of population demographic and genetic theory indicates that mitochondrial and nuclear phylogeographical results ought to concur for both geographically unstructured populations and for populations that have long histories of isolation. However, a relatively common occurrence will be shallow, but geographically structured mtDNA trees--without nuclear gene corroboration--for populations with relatively shorter periods of isolation. This is expected because of the longer coalescence times of nuclear genes (approximately four times that of mtDNA); such cases do not contradict the mtDNA inference of recent isolation and evolutionary divergence. Rather, the nuclear markers are more lagging indicators of changes in population structure. A review of the recent literature on birds reveals the existence of relatively few cases in which nuclear markers contradict mitochondrial markers in a fashion not consistent with coalescent theory. Preliminary information from nuclear genes suggests that mtDNA patterns will prove to be robust indicators of patterns of population history and species limits. At equilibrium, mitochondrial loci are generally a more sensitive indicator of population structure than are nuclear loci, and mitochondrial estimates of F(ST)-like statistics are generally expected to exceed nuclear ones. Hence, invoking behavioural or ecological explanations of such differences is not parsimonious. Nuclear genes will prove important for quantitative estimates of the depths of haplotype trees, rates of population growth and values of gene flow.     

Geography disentangles introgression from ancestral polymorphism in Lake Malawi cichlids

Phenotypically diverse Lake Malawi cichlids exhibit similar genomes. The extensive sharing of genetic polymorphism among forms has both intrigued and frustrated biologists trying to understand the nature of diversity in this and other rapidly evolving systems. Shared polymorphism might result from hybridization and/or the retention of ancestrally polymorphic alleles. To examine these alternatives, we used new genomic tools to characterize genetic differentiation in widespread, geographically structured populations of Labeotropheus fuelleborni and Metriaclima zebra. These phenotypically distinct species share mitochondrial DNA (mtDNA) haplotypes and show greater mtDNA differentiation among localities than between species. However, Bayesian analysis of nuclear single nucleotide polymorphism (SNP) data revealed two distinct genetic clusters corresponding perfectly to morphologically diagnosed L. fuelleborni and M. zebra. This result is a function of the resolving power of the multi‐locus dataset, not a conflict between nuclear and mitochondrial partitions. Locus‐by‐locus analysis showed that mtDNA differentiation between species (F_CT) was nearly identical to the median single‐locus SNP F_CT. Finally, we asked whether there is evidence for gene flow at sites of co‐occurrence. We used simulations to generate a null distribution for the level of differentiation between co‐occurring populations of L. fuelleborni and M. zebra expected if there was no hybridization. The null hypothesis was rejected for the SNP data; populations that co‐occur at rock reef sites were slightly more similar than expected by chance, suggesting recent gene flow. The coupling of numerous independent markers with extensive geographic sampling and simulations utilized here provides a framework for assessing the prevalence of gene flow in recently diverged species.     

Extensive sampling of polar bears (Ursus maritimus) in the Northwest Passage (Canadian Arctic Archipelago) reveals population differentiation across multiple spatial and temporal scales

As global warming accelerates the melting of Arctic sea ice, polar bears (Ursus maritimus) must adapt to a rapidly changing landscape. This process will necessarily alter the species distribution together with population dynamics and structure. Detailed knowledge of these changes is crucial to delineating conservation priorities. Here, we sampled 361 polar bears from across the center of the Canadian Arctic Archipelago spanning the Gulf of Boothia (GB) and M'Clintock Channel (MC). We use DNA microsatellites and mitochondrial control region sequences to quantify genetic differentiation, estimate gene flow, and infer population history. Two populations, roughly coincident with GB and MC, are significantly differentiated at both nuclear (F_ST = 0.01) and mitochondrial (Φ_ST = 0.47; F_ST = 0.29) loci, allowing Bayesian clustering analyses to assign individuals to either group. Our data imply that the causes of the mitochondrial and nuclear genetic patterns differ. Analysis of mtDNA reveals the matrilineal structure dates at least to the Holocene, and is common to individuals throughout the species’ range. These mtDNA differences probably reflect both genetic drift and historical colonization dynamics. In contrast, the differentiation inferred from microsatellites is only on the scale of hundreds of years, possibly reflecting contemporary impediments to gene flow. Taken together, our data suggest that gene flow is insufficient to homogenize the GB and MC populations and support the designation of GB and MC as separate polar bear conservation units. Our study also provide a striking example of how nuclear DNA and mtDNA capture different aspects of a species demographic history.     

Mitochondrial regulation of epigenetics and its role in human diseases

Most pathogenic mitochondrial DNA (mtDNA) mutations induce defects in mitochondrial oxidative phosphorylation (OXPHOS). However, phenotypic effects of these mutations show a large degree of variation depending on the tissue affected. These differences are difficult to reconcile with OXPHOS as the sole pathogenic factor suggesting that additional mechanisms contribute to lack of genotype and clinical phenotype correlationship. An increasing number of studies have identified a possible effect on the epigenetic landscape of the nuclear genome as a consequence of mitochondrial dysfunction. In particular, these studies demonstrate reversible or irreversible changes in genomic DNA methylation profiles of the nuclear genome. Here we review how mitochondria damage checkpoint (mitocheckpoint) induces epigenetic changes in the nucleus. Persistent pathogenic mutations in mtDNA may also lead to epigenetic changes causing genomic instability in the nuclear genome. We propose that “mitocheckpoint” mediated epigenetic and genetic changes may play key roles in phenotypic variation related to mitochondrial diseases or host of human diseases in which mitochondrial defect plays a primary role.     

Contrasting population structure from nuclear intron sequences and mtDNA of humpback whales

Powerful analyses of population structure require information from multiple genetic loci. To help develop a molecular toolbox for obtaining this information, we have designed universal oligonucleotide primers that span conserved intron-exon junctions in a wide variety of animal phyla. We test the utility of exon-primed, intron-crossing amplifications by analyzing the variability of actin intron sequences from humpback, blue, and bowhead whales and comparing the results with mitochondrial DNA (mtDNA) haplotype data. Humpback actin introns fall into two major clades that exist in different frequencies in different oceanic populations. It is surprising that Hawaii and California populations, which are very distinct in mtDNAs, are similar in actin intron alleles. This discrepancy between mtDNA and nuclear DNA results may be due either to differences in genetic drift in mitochondrial and nuclear genes or to preferential movement of males, which do not transmit mtDNA to offspring, between separate breeding grounds. Opposing mtDNA and nuclear DNA results can help clarify otherwise hidden patterns of structure in natural populations.      

Patterns of population subdivision, gene flow and genetic variability in the African wild dog (Lycaon pictus)

African wild dogs are large, highly mobile carnivores that are known to disperse over considerable distances and are rare throughout much of their geographical range. Consequently, genetic variation within and differentiation between geographically separated populations is predicted to be minimal. We determined the genetic diversity of mitochondrial DNA (mtDNA) control region sequences and microsatellite loci in seven populations of African wild dogs. Analysis of mtDNA nucleotide diversity suggests that, historically, wild dog populations have been small relative to other large carnivores. However, population declines due to recent habitat loss have not caused a dramatic reduction in genetic diversity. We found one historical and eight recent mtDNA genotypes in 280 individuals that defined two highly divergent clades. In contrast to a previous, more limited, mtDNA analysis, sequences from these clades are not geographically restricted to eastern or southern African populations. Rather, we found a large admixture zone spanning populations from Botswana, Zimbabwe and south-eastern Tanzania. Mitochondrial and microsatellite differentiation between populations was significant and unique mtDNA genotypes and alleles characterized the populations. However, gene flow estimates (Nm) based on microsatellite data were generally greater than one migrant per generation. In contrast, gene flow estimates based on the mtDNA control region were lower than expected given differences in the mode of inheritance of mitochondrial and nuclear markers which suggests a male bias in long-distance dispersal.     

Genetic and morphometric differentiation between island and mainland southern elephant seal populations

We compare genetic (both nuclear and mitochondrial) and morphometric measures between two putative populations of southern elephant seal (Mirounga leonina), and interpret the results in the context of data from mark-recapture and satellite-telemetric studies. One population is on the Argentine mainland, while the other is 2,400 km away on South Georgia island. We found pronounced differentiation at the mitochondrial DNA (mtDNA) control region that was distinct from the pattern of variation seen among island rookeries. Some morphometric characters and seven out of ten nuclear-DNA markers also showed differentiation between the island and mainland sites. Diversity at nuclear markers was high in both populations but mtDNA diversity was low in the mainland population, suggesting a founder event and little subsequent immigration of females. Morphological differences may suggest different selective environments at the two sites.     

Association between mitochondrial DNA and morphological evolution in Canada geese

Summary Variation within and between eight subspecies of Canada geese was assessed by restriction fragment analysis of mitochondrial DNA (mtDNA), electrophoresis of proteins encoded by nuclear DNA, and the morphometric analysis of skeletons. Estimates of mtDNA sequence divergence between Canada goose subspecies ranged from 0.04 to 2.54%. Pairwise comparisons of the three data matrices revealed that only mtDNA variation and body size are significantly correlated. Subspecies with northern breeding grounds are small-bodied and display small variations of one mtDNA clone, whereas those breeding further south are largebodied and show small differences in another mtDNA clone. Canada geese exhibit strong geographic differentiation with respect to mtDNA sequence, but weak structuring in protein-encoding nuclear DNA. This finding can be explained by a lower level of gene flow for the mitochondrial genome than for the nuclear genome, which in turn emanates from the maternal inheritance of mtDNA and male-biased dispersal in Canada geese. Despite male-mediated flow of nuclear genes, strong morphometric differentiation persists among Canada geese subspecies.     

Mitochondrial DNA variation is associated with measurable differences in life-history traits and mitochondrial metabolism in Drosophila simulans

Recent studies have used a variety of theoretical arguments to show that mitochondrial (mt) DNA rarely evolves as a strictly neutral marker and that selection operates on the mtDNA of many species. However, the vast majority of researchers are not convinced by these arguments because data linking mtDNA variation with phenotypic differences are limited. We investigated sequence variation in the three mtDNA and nine nuclear genes (including all isoforms) that encode the 12 subunits of cytochrome c oxidase of the electron transport chain in Drosophila. We then studied cytochrome c oxidase activity as a key aspect of mitochondrial bioenergetics and four life-history traits. In Drosophila simulans, sequence data from the three mtDNA encoded cytochrome c oxidase genes show that there are 76 synonymous and two nonsynonymous fixed differences among flies harboring siII compared with siIII mtDNA. In contrast, 13 nuclear encoded genes show no evidence of genetic subdivision associated with the mtDNA. Flies with siIII mtDNA had higher cytochrome c oxidase activity and were more starvation resistant. Flies harboring siII mtDNA had greater egg size and fecundity, and recovered faster from cold coma. These data are consistent with a causative role for mtDNA variation in these phenotypic differences, but we cannot completely rule out the involvement of nuclear genes. The results of this study have significant implications for the use of mtDNA as an assumed neutral marker and show that evolutionary shifts can involve changes in mtDNA despite the small number of genes encoded in the organelle genome.     

DNA extraction procedures meaningfully influence qPCR-based mtDNA copy number determination

Quantitative real time PCR (qPCR) is commonly used to determine cell mitochondrial DNA (mtDNA) copy number. This technique involves obtaining the ratio of an unknown variable (number of copies of an mtDNA gene) to a known parameter (number of copies of a nuclear DNA gene) within a genomic DNA sample. We considered the possibility that mtDNA:nuclear DNA (nDNA) ratio determinations could vary depending on the method of genomic DNA extraction used, and that these differences could substantively impact mtDNA copy number determination via qPCR. To test this we measured mtDNA:nDNA ratios in genomic DNA samples prepared using organic solvent (phenol–chloroform–isoamyl alcohol) extraction and two different silica-based column methods, and found mtDNA:nDNA ratio estimates were not uniform. We further evaluated whether different genomic DNA preparation methods could influence outcomes of experiments that use mtDNA:nDNA ratios as endpoints, and found the method of genomic DNA extraction can indeed alter experimental outcomes. We conclude genomic DNA sample preparation can meaningfully influence mtDNA copy number determination by qPCR.     

Effects of differentiation of embryonal carcinoma cells (P19) on mitochondrial DNA content in vitro

Summary The embryonal carcinoma cell line P19 is derived from mouse teratocarcinomas. These pluripotent cells can be induced to differentiate into a variety of cell types by exposure to various drugs. We used retinoic acid to induce embryonal carcinoma cells to differentiate into neuronlike cells. In this study, we show that changes occur in mitochondria during differentiation of embryonal carcinoma cells to neuronlike cells. We found that various morphologic parameters such as mitochondrial fractional area and mitochondrial size decrease as embryonal carcinoma cells differentiate into neuronlike cells. Similar changes were also observed in mitochondrial DNA content. Stereologic analysis of cell preparations provided a measure of mitochondrial fractional area per cell and mtDNA content was assessed by radiolabeled mtDNA probe. This study establishes that mitochondria are regulated as cells differentiate. This study was financially supported by the Medical Research Council of Canada.     

Deletion of mtDNA disrupts mitochondrial function and structure, but not biogenesis

The role of nuclear DNA (nDNA)-encoded proteins in the regulation of mitochondrial fission and fusion has been documented, yet the role of mitochondrial DNA (mtDNA) and encoded proteins in mitochondrial biogenesis remains unknown. Long-term treatment of a lymphoblastoid cell line Molt-4 with ethidium bromide generated mtDNA-deficient rho0 mutants. Depletion of mtDNA in rho0 cells produced functional and morphological changes in mitochondria without affecting the nuclear genome and encoded proteins. Indeed, the gene encoding subunit II of mitochondrial cytochrome c oxidase (COX II), a prototypical mitochondrial gene, was reduced in rho0 mutants blunting the activity of mitochondrial cytochrome coxidase. Yet, the amount of the nuclear beta-actin gene and the activity of citrate synthase, a mitochondrial matrix enzyme encoded by nDNA, remained unaffected in rho0 cells. Loss of mtDNA in rho0 cells was associated with significant distortion of mitochondrial structure, decreased electron density of the matrix and disorganized inner and outer membranes, resulting in the appearance of 'ghost-like' mitochondria. However, the number of mitochondria-like structures was not significantly different between mtDNA-deficient and parental cells. Thus, we conclude that cells lacking mtDNA still generate mitochondrial scaffolds, albeit with aberrant function.