The microevolution of the Galápagos marine iguana Amblyrhynchus cristatus assessed by nuclear and mitochondrial genetic analyses

Marine iguanas may have inhabited the Galápagos archipelago and its former, now sunken islands for more than 10 million years (Myr). It is therefore surprising that morphological and immunological data indicate little evolutionary divergence within the genus. We utilized mitochondrial DNA (mtDNA) sequence analyses and nuclear DNA fingerprinting to re-evaluate the level and pattern of genetic differentiation among 22 marine iguana populations from throughout the archipelago. Both genetic marker systems detect a low level of within-genus divergence, but they show contrasting levels of geographical subdivision among the populations. The mitochondrial gene pools of populations from different regions of the archipelago are isolated, and the mtDNA pattern appears to follow the sequence in which the islands were colonized by marine iguanas. Conversely, the nuclear DNA study indicates substantial interpopulational gene exchange, and the geographical distribution of the nuclear markers seems to be determined by isolation by distance among the populations. The natural history of marine iguanas suggests that the contrasting nuclear and mitochondrial DNA patterns result from an asymmetric migration behaviour of the two sexes, with higher (active and passive) interisland dispersal for males than females. Separate genetic analyses for the sexes appear to support this hypophesis. Based on these findings, a scenario is proposed that explains the marine iguanas' low genetic divergence, notwithstanding their long evolutionary history in the Galápagos archipelago.     

Mechanical barriers to introgressive hybridization revealed by mitochondrial introgression patterns in Ohomopterus ground beetle assemblages

To reveal the role of diverged body size and genital morphology in reproductive isolation among closely related species, we examined patterns of, and factors limiting, introgressive hybridization between sympatric Ohomopterus ground beetles in central Japan using mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequences. We sampled 17 local assemblages that consisted of two to five species and estimated levels of interspecific gene flow using the genetic distance, D(A), and maximum-likelihood estimates of gene flow. Sharing of haplotypes or haplotype lineages was detected between six of seven species that occurred in the study areas, indicating mitochondrial introgression. The intensity and direction of mitochondrial gene flow were variable among species pairs. To determine the factors affecting introgression patterns, we tested the relationships between interspecific D(A) and five independent variables: difference in body size, difference in genital size, phylogenetic relatedness (nuclear gene sequence divergence), habitat difference, and species richness of the assemblage. Body and genital size differences contributed significantly to preventing gene flow. Thus, mechanical isolation mechanisms reduce the chance of introgressive hybridization between closely related species. Our results highlight the role of morphological divergence in speciation and assemblage formation processes through mechanical isolation.     

On shortcomings of using mtDNA sequence divergence for the systematics of hares (genus Lepus): An example from cape hares

Despite the well-known fact that evolutionary patterns of single genes or sequences are not necessarily paralleled by organismic evolution, using mitochondrial sequence divergence for inferring phylogenetic relationships among hare species is not uncommon. Earlier studies reported interspecific introgression in the genus Lepus that may slur systematic conclusions drawn exclusively from mtDNA data. We examined pairwise divergence in partial mtHV-1 sequences and nuclear gene pools based on microsatellite and allozyme loci separately in a South and North African population of cape hares, Lepus capensis, and did not find significant correspondence on the individual level within either population. Also, the former population had a significantly higher level of mtHV-1 sequence divergence compared to the latter, but levels of individual nuclear gene pool divergence did not differ significantly between the two populations. Hence, the absence of correspondence between mtHV-1 sequence divergence and nuclear gene pool divergence among individuals within a population might be independent of the population-specific level of differentiation among individuals. Our results complement earlier findings in hares, and strongly recommend to include nuclear gene pool evidence for systematic inferences within this genus, even under the absence of mitochondrial introgression.     

Decoupling of molecular and morphological evolution in deep lineages of a meiobenthic harpacticoid copepod

Molecular and biochemical genetic analyses have revealed that many marine invertebrate taxa, including some well-studied and presumably cosmopolitan species, are actually complexes of sibling species. When morphological differences are slight and estimated divergence times are old, data suggest either unusually high rates of sequence evolution or long-term morphological stasis. Here, five gene regions (mitochondrial cytochrome oxidase subunit I and large-subunit ribosomal 16S rDNA and nuclear ITS1, 5.8S rDNA, and ITS2) were analyzed in four geographic samples of the meiobenthic harpacticoid copepod Cletocamptus deitersi. Molecular sequences revealed four extremely differentiated molecular lineages with unalignable nuclear intergenic spacers and mitochondrial uncorrected divergences reaching 25% (cytochrome oxidase) and 36% (16S rDNA). These levels of divergence are greater than those reported previously for congeneric species in diverse invertebrate taxa, including crustaceans. The nominally intraspecific divergence matches or exceeds the corresponding divergence from a known congener (Cletocamptus helobius). A molecular clock applied to the cytochrome oxidase subunit I data suggests that these lineages split in the Miocene, consistent with the fossil record of a North American Cletocamptus from the same period. Morphological differences among the major lineages are subtle but congruent with the patterns of genetic differentiation. Our conclusion, based on concordant patterns of variation in two mitochondrial and three nuclear gene regions, as well as morphological observations, is that C. deitersi in North America is composed of at least four separate species by the genealogical concordance, phylogenetic, and morphological-species criteria. Alternative explanations for the deep phylogenetic nodes and apparent morphological stasis, including high rates of sequence evolution, balancing selection, and genetic signatures of historical events, are considered unlikely.     

Pleistocene separation of mitochondrial lineages of Mytilus spp. mussels from Northern and Southern Hemispheres and strong genetic differentiation among southern populations

Smooth-shelled mussels, Mytilus spp., have an antitropical distribution. In the Northern Hemisphere, the M. edulis complex of species is composed of three genetically well delineated taxa: M. edulis, M. galloprovincialis and M. trossulus. In the Southern Hemisphere, morphological characters, allozymes and intron length polymorphisms suggest that Mytilus spp. populations from South America and Kerguelen Islands are related to M. edulis and those from Australasia to M. galloprovincialis. On the other hand, a phylogeny of the 16S rDNA mitochondrial locus demonstrates a clear distinctiveness of southern mussels and suggests that they are related to Mediterranean M. galloprovincialis. Here, we analysed the faster-evolving cytochrome oxidase subunit I locus. The divergence between haplotypes of populations from the two hemispheres was confirmed and was found to predate the divergence between haplotypes of northern M. edulis and M. galloprovincialis. In addition, strong genetic structure was detected among the southern samples, revealing three genetic entities that correspond to (1) South America and Kerguelen Island, (2) Tasmania, (3) New Zealand. Using the trans-Arctic interchange as a molecular clock calibration, we estimated the time since divergence of populations from the two hemispheres to be between 0.5 million years (MY) and 1.3 MY (average 0.84 MY). The contrasting patterns observed for the nuclear and the organelle genomes suggested two alternative, complex scenarios: two trans-equatorial migrations and the existence of differential barriers to mitochondrial and nuclear gene flow, or a single trans-equatorial migration and a view of the composition of the nuclear genome biased by taxonomic preconception.     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

Genetic structure of introduced populations: 120‐year‐old DNA footprint of historic introduction in an insular small mammal population

Wildlife populations have been introduced to new areas by people for centuries, but this human‐mediated movement can disrupt natural patterns of genetic structure by altering patterns of gene flow. Insular populations are particularly prone to these influences due to limited opportunities for natural dispersal onto islands. Consequently, understanding how genetic patterns develop in island populations is important, particularly given that islands are frequently havens for protected wildlife. We examined the evolutionary origins and extent of genetic structure within the introduced island population of red squirrels (Sciurus vulgaris) on the Channel Island of Jersey using mitochondrial DNA (mtDNA) control region sequence and nuclear microsatellite genotypes. Our findings reveal two different genetic origins and a genetic architecture reflective of the introductions 120 years ago. Genetic structure is marked within the maternally inherited mtDNA, indicating slow dispersal of female squirrels. However, nuclear markers detected only weak genetic structure, indicating substantially greater male dispersal. Data from both mitochondrial and nuclear markers support historic records that squirrels from England were introduced to the west of the island and those from mainland Europe to the east. Although some level of dispersal and introgression across the island between the two introductions is evident, there has not yet been sufficient gene flow to erase this historic genetic “footprint.” We also investigated if inbreeding has contributed to high observed levels of disease, but found no association. Genetic footprints of introductions can persist for considerable periods of time and beyond traditional timeframes of wildlife management.     

Genetic structure in a montane mayfly Baetis bicaudatus (Ephemeroptera: Baetidae), from the Rocky Mountains, Colorado

1. Populations of a number of sub‐tropical stream insect species have been found to show unexpected patterns of genetic variation, with more differences between samples from the same stream than between whole streams or between subcatchments. Many samples also showed deviations from Hardy–Weinberg proportions. It has been proposed that these patterns result from limited numbers of matings contributing to a given stream reach, because adults emerge throughout the year, and low levels of larval drift between reaches. These patterns may be less likely in a northern hemisphere montane species with synchronous emergence of adults and high levels of drift. We tested the hypothesis that patterns of genetic variation in a montane mayfly from the Rocky Mountains, Colorado, would reflect a pattern of ‘isolation by distance’ with samples from the same creek being more similar than samples from different creeks and that deviations from Hardy–Weinberg proportions would be minimal. 2. Based on allozyme variation, the hypothesis of minimal deviations from Hardy–Weinberg proportions was not supported and there was no evidence of isolation by distance. Nevertheless the levels of differentiation among samples from within the same stream were less than those reported for most subtropical species. 3. Results from analysis of a fragment of the cytochrome oxidase gene (subunit 1) revealed contrasting patterns. The levels of genetic differentiation were an order of magnitude higher between streams than among samples within streams. In addition, although there was no significant isolation by distance effect overall, a nested clade analysis provided evidence for restricted gene flow with isolation by distance for some clades. 4. We suggest that these contrasting results may reflect the differences in male and female dispersal patterns. While differentiation at nuclear gene markers (allozymes) give information about both male and female dispersal, mitochondrial DNA markers reflect only female dispersal. We suggest that in this species, female dispersal may be more restricted, perhaps mostly along stream channels, whereas male dispersal is more widespread. An alternative explanation for the different results is the different evolutionary rates of the mitochondrial and nuclear markers.     

Rates of molecular evolution and the fraction of nucleotide positions free to vary

Summary Selective constraints on DNA sequence change were incorporated into a model of DNA divergence by restricting substitutions to a subset of nucleotide positions. A simple model showed that both mutation rate and the fraction of nucleotide positions free to vary are strong determinants of DNA divergence over time.When divergence between two species approaches the fraction of positions free to vary, standard methods that correct for multiple mutations yield severe underestimates of the number of substitutions per site. A modified method appropriate for use with DNA sequence, restriction site, or thermal renaturation data is derived taking this fraction into account. The model also showed that the ratio of divergence in two gene classes (e.g., nuclear and mitochondrial) may vary widely over time even if the ratio of mutation rates remains constant.DNA sequence divergence data are used increasingly to detect differences in rates of molecular evolution. Often, variation in divergence rate is assumed to represent variation in mutation rate. The present model suggests that differing divergence rates among comparisons (either among gene classes or taxa) should be interpreted cautiously. Differences in the fraction of nucleotide positions free to vary can serve as an important alternative hypothesis to explain differences in DNA divergence rates.     

World phylogeography and male‐mediated gene flow in the sandbar shark,Carcharhinus plumbeus

The sandbar shark, Carcharhinus plumbeus, is a large, cosmopolitan, coastal species. Females are thought to show philopatry to nursery grounds while males potentially migrate long distances, creating an opportunity for male‐mediated gene flow that may lead to discordance in patterns revealed by mitochondrial DNA (mtDNA) and nuclear markers. While this dynamic has been investigated in elasmobranchs over small spatial scales, it has not been examined at a global level. We examined patterns of historical phylogeography and contemporary gene flow by genotyping 329 individuals from nine locations throughout the species’ range at eight nuclear microsatellite markers and sequencing the complete mtDNA control region. Pairwise comparisons often resulted in fixation indices and divergence estimates of greater magnitude using mtDNA sequence data than microsatellite data. In addition, multiple methods of estimation suggested fewer populations based on microsatellite loci than on mtDNA sequence data. Coalescent analyses suggest divergence and restricted migration among Hawaii, Taiwan, eastern and western Australia using mtDNA sequence data and no divergence and high migration rates, between Taiwan and both Australian sites using microsatellite data. Evidence of secondary contact was detected between several localities and appears to be discreet in time rather than continuous. Collectively, these data suggest complex spatial/temporal relationships between shark populations that may feature pulses of female dispersal and more continuous male‐mediated gene flow.     

HALDANE'S RULE AND SEX BIASSED GENE FLOW BETWEEN TWO HYBRIDIZING FLYCATCHER SPECIES (FICEDULA ALBICOLLIS AND F. HYPOLEUCA,AVES: MUSCICAPIDAE)

The collared flycatcher (Ficedula albicollis) and the pied flycatcher (F. hypoleuca) hybridize where their geographic ranges overlap. Restriction fragment comparison of 5% of the mitochondrial genome showed a sequence divergence of 10% between these flycatcher species. This degree of sequence divergence between a closely related pair of bird species is unusually high and contrasts with the low level of divergence between F. albicollis and F. hypoleuca in nuclear genes (Nei's D = 0.0006) revealed by enzyme electrophoresis. The low nuclear differentiation is explained by sex biassed gene flow and introgression in nuclear genes (via fertile male hybrids), while the high mitochondrial DNA sequence divergence is preserved by sterility of female hybrids, which prevents mitochondrial introgression. This pattern is in accordance with Haldane's rule and is supported by field data on hybrid fertility. The high mtDNA differentiation could be explained by transfer of mitochondrial DNA from a third species during a past period of hybridization.     

PREDICTING NUCLEAR GENE COALESCENCE FROM MITOCHONDRIAL DATA: THE THREE-TIMES RULE

Abstract.— Coalescence theory predicts when genetic drift at nuclear loci will result in fixation of sequence differences to produce monophyletic gene trees. However, the theory is difficult to apply to particular taxa because it hinges on genetically effective population size, which is generally unknown. Neutral theory also predicts that evolution of monophyly will be four times slower in nuclear than in mitochondrial genes primarily because genetic drift is slower at nuclear loci. Variation in mitochondrial DNA (mtDNA) within and between species has been studied extensively, but can these mtDNA data be used to predict coalescence in nuclear loci? Comparison of neutral theories of coalescence of mitochondrial and nuclear loci suggests a simple rule of thumb. The “three‐times rule” states that, on average, most nuclear loci will be monophyletic when the branch length leading to the mtDNA sequences of a species is three times longer than the average mtDNA sequence diversity observed within that species. A test using mitochondrial and nuclear intron data from seven species of whales and dolphins suggests general agreement with predictions of the three‐times rule. We define the coalescence ratio as the mitochondrial branch length for a species divided by intraspecific mtDNA diversity. We show that species with high coalescence ratios show nuclear monophyly, whereas species with low ratios have polyphyletic nuclear gene trees. As expected, species with intermediate coalescence ratios show a variety of patterns. Especially at very high or low coalescence ratios, the three‐times rule predicts nuclear gene patterns that can help detect the action of selection. The three‐times rule may be useful as an empirical benchmark for evaluating evolutionary processes occurring at multiple loci.     

Mitochondrial malate dehydrogenases in Brassica napus: altered protein patterns in different nuclear mitochondrial combinations

Two-dimensional analyses of mitochondrial proteins of Brassica napus revealed a set of differences in patterns of mitochondrial matrix proteins isolated from different nuclear backgrounds. One of these varying proteins was identified as mitochondrial malate dehydrogenase (mMDH; EC 1.1.1.37) by homology analyses of the partial amino acid sequence. Immunological detection identified additional mMDH subunits and detected different patterns of mMDH subunits in two distinct mitochondria types although they were isolated from plants with the same nuclear genotype. These differences are also reflected in isozym patterns, whereas Southern analyses showed no alteration in genome structure. Therefore mitochondria type-specific mMDH modifications are possible.     

POPULATION STRUCTURE OF THE BOTTLENOSE DOLPHIN (TURSIOPS TRUNCATUS) AS DETERMINED BY RESTRICTION ENDONUCLEASE ANALYSIS OF MITOCHONDRIAL DNA

Abstract: Restriction fragment length polymorphisms of mitochondrial DNA (mtDNA) were used to test for population subdivision in the bottlenose dolphin (Tursiops truncatus). Atlantic and Pacific dolphin mtDNA samples exhibited distinctly different haplotypes (approximately 2.4% sequence divergence), indicating a lack of gene exchange. Within the Atlantic Ocean, mtDNA samples from the Gulf of Mexico and the Atlantic Coast were also found to be distinct, with a sequence divergence of approximately 0.6%. The Atlantic Coast–Gulf of Mexico dichotomy is consistent with patterns of genetic variation from other marine and coastal organisms from this region, and supports the hypothesized role of bio-geographic events in promoting the divergence of these and other forms. Regional differentiation was identified along the Atlantic Coast, whereas low sequence divergences among haplotypes and consistent haplotype frequencies across populations suggested considerable gene exchange among Gulf of Mexico populations. A highly divergent haplotype found in two individuals from two localities in the Gulf of Mexico is best explained by dispersal from either a distinct offshore Gulf stock or an unsampled Atlantic Coast stock. Additional samples are required to test for the existence of a distinct offshore race and, if it exists, to identify its distribution and contribution to population structure.     

Comparative analysis of gender-associated complete mitochondrial genomes in marine mussels (Mytilus spp.)

Marine mussels of the genus Mytilus have an unusual mode of mitochondrial DNA (mtDNA) transmission termed doubly uniparental inheritance (DUI). Female mussels are homoplasmic for the F mitotype, which is inherited maternally, while males are usually heteroplasmic, carrying a mixture of the maternal F mitotype and the paternally inherited M genome. Two classes of M genomes have been observed: "standard" M genomes and "recently masculinized" M genomes. The latter are more similar to F genomes at the sequence level but are transmitted paternally like standard M genomes. In this study we report the complete sequences of two standard male M. edulis and one recently masculinized male M. trossulus mitochondrial genome. A comparative analysis, including the previously sequenced M. edulis F and M. galloprovincialis F and M mtDNAs, reveals that these genomes are identical in gene order, but highly divergent in nucleotide and amino acid sequence. The large amount (>20%) of nucleotide substitutions that fall in coding regions implies that there are several amino acid replacements between the F and M genomes, which likely have an impact on the structural and functional properties of the mitochondrial proteome. Correlation of the divergence rate of different protein-coding genes indicates that mtDNA-encoded proteins of the M genome are still under selective constraints, although less highly than genes of the F genome. The mosaic F/M control region of the masculinized F genome provides evidence for lineage-specific sequences that may be responsible for the different mode of transmission genetics. This analysis shows the value of comparative genomics to better understand the mechanisms of maintenance and segregation of mtDNA sequence variants in mytilid mussels.     

Range‐wide genetic homogeneity in the California sea mussel (Mytilus californianus): a comparison of allozymes,nuclear DNA markers,and mitochondrial DNA sequences

We tested for genetic differentiation among six populations of California sea mussels (Mytilus californianus) sampled across 4000 km of its geographical range by comparing patterns of variation at four independent types of genetic markers: allozymes, single‐copy nuclear DNA markers, and DNA sequences from the male and female mitochondrial genomes. Despite our extensive sampling and genotyping efforts, we detected no significant differences among localities and no signal of isolation by distance suggesting that M. californianus is genetically homogeneous throughout its range. This concordance differs from similar studies on other mytilids, especially in the role of postsettlement selection generating differences between exposed coastal and estuarine habitats. To assess if this homogeneity was due to M. californianus not inhabiting estuarine environments, we reviewed studies comparing allozymes with other classes of nuclear DNA markers. Although both types of markers gave broadly consistent results, there was a bias favouring studies in which allozymes were more divergent than DNA markers (nine to three) and a disproportionate number of these cases involved marine taxa (seven). Furthermore, allozymes were significantly more heterogeneous than DNA markers in three of the four studies that sampled coastal and estuarine habitats. We conclude that the genetic uniformity exhibited by M. californianus may result from a combination of extensive gene flow and the lack of exposure to strong selective gradients across its range.     

Rapid evolution of postzygotic reproductive isolation in stalk-eyed flies

We test the relative rates of evolution of pre- and postzygotic reproductive isolation using eight populations of the sexually dimorphic stalk-eyed flies Cyrtodiopsis dalmanni and C. whitei. Flies from these populations exhibit few morphological differences yet experience strong sexual selection on male eyestalks. To measure reproductive isolation we housed one male and three female flies from within and between these populations in replicate cages and then recorded mating behavior, sperm transfer, progeny production, and hybrid fertility. Using a phylogeny based on partial sequences of two mitochondrial genes, we found that premating isolation, postmating isolation prior to hybrid eclosion, and female hybrid sterility evolve gradually with respect to mitochondrial DNA sequence divergence. In contrast, male hybrid sterility evolves much more rapidly-at least twice as fast as any other form of reproductive isolation. Hybrid sterility, therefore, obeys Haldane's rule. Although some brood sex ratios were female biased, average brood sex ratio did not covary with genetic distance, as would be expected if hybrid inviability obeyed Haldane's rule. The likelihood that forces including sexual selection and intra- and intergenomic conflict may have contributed to these patterns is discussed.     

Differential patterns of male and female mtDNA exchange across the Atlantic Ocean in the blue mussel, Mytilus edulis

Comparisons among loci with differing modes of inheritance can reveal unexpected aspects of population history. We employ a multilocus approach to ask whether two types of independently assorting mitochondrial DNAs (maternally and paternally inherited: F- and M-mtDNA) and a nuclear locus (ITS) yield concordant estimates of gene flow and population divergence. The blue mussel, Mytilus edulis, is distributed on both North American and European coastlines and these populations are separated by the waters of the Atlantic Ocean. Gene flow across the Atlantic Ocean differs among loci, with F-mtDNA and ITS showing an imprint of some genetic interchange and M-mtDNA showing no evidence for gene flow. Gene flow of F-mtDNA and ITS causes trans-Atlantic population divergence times to be greatly underestimated for these loci, although a single trans-Atlantic population divergence time (1.2 MYA) can be accommodated by considering all three loci in combination in a coalescent framework. The apparent lack of gene flow for M-mtDNA is not readily explained by different dispersal capacities of male and female mussels. A genetic barrier to M-mtDNA exchange between North American and European mussel populations is likely to explain the observed pattern, perhaps associated with the double uniparental system of mitochondrial DNA inheritance.