THE GENETIC STRUCTURE OF A GYNODIOECIOUS PLANT: NUCLEAR AND CYTOPLASMIC GENES

Sex expression in gynodioecious plants is often determined by an interaction between biparentally and maternally inherited genes. Their relative rates of gene flow should be considered when modeling the evolution of the sex ratio in structured populations. In order to understand patterns of gene flow in Silene vulgaris, a gynodioecious plant, genetic structure was estimated from biparentally inherited genetic markers (allozymes) and a maternally inherited marker (chloroplast DNA) using Wright's F_st. Based on data from 16 local populations, chloroplast DNA showed considerably more genetic structure than did allozymes (F_st values of 0.62 and 0.22, respectively). This suggests that the rate of gene flow is about three times greater for nuclear genes.     

Genetic differentiation and history of populations of the Italian treefrog Hyla intermedia: lack of concordance between mitochondrial and nuclear markers

The genetic differentiation among 33 populations of the Italian treefrog, Hyla intermedia (Anura: Hylidae), was investigated using both biparentally (23 allozyme loci) and maternally (partial mitochondrial cytochrome b gene) inherited markers. Two main population groups were evidenced by both markers, located north and south of the northern Apennines. However, the pattern of differentiation between these two groups was much less pronounced at allozymes than at mtDNA, leading to gene flow estimates that were 25 times lower at mitochondrial than at nuclear level. Also, the mtDNA divergence between the two groups was particularly marked for two cospecific lineages of anuran amphibians (the P-distance being on average 9.04%), while their average genetic distance at allozymes was comparatively low (D _NEI = 0.07). This contrasting pattern of nuclear versus mitochondrial genetic variation is discussed in the context of: (1) marker specific selection, (2) secondary contact and sex-biased gene flow and (3) ancestral polymorphism and colonization from north to south. Finally we emphasize how, for population genetic studies, the use of multiple markers having distinct evolutionary properties can help unravel the existence of more complex evolutionary histories.     

Conservation implications of complex population structure: lessons from the loggerhead turtle (Caretta caretta)

Complex population structure can result from either sex-biased gene flow or population overlap during migrations. Loggerhead turtles (Caretta caretta) have both traits, providing an instructive case history for wildlife management. Based on surveys of maternally inherited mtDNA, pelagic post-hatchlings show no population structure across the northern Atlantic (phi(ST) < 0.001, P = 0.919), subadults in coastal habitat show low structure among locations (phi(ST) = 0.01, P < 0.005), and nesting colonies along the southeastern coast of the United States have strong structure (phi(ST) = 0.42, P < 0.001). Thus the level of population structure increases through progressive life history stages. In contrast, a survey of biparentally inherited microsatellite DNA shows no significant population structure: R(ST) < 0.001; F(ST) = 0.002 (P > 0.05) across the same nesting colonies. These results indicate that loggerhead females home faithfully to their natal nesting colony, but males provide an avenue of gene flow between regional nesting colonies, probably via opportunistic mating in migratory corridors. As a result, all breeding populations in the southeastern United States have similar levels of microsatellite diversity (H(E) = 0.70-0.89), whereas mtDNA haplotype diversity varies dramatically (h = 0.00-0.66). Under a conventional interpretation of the nuclear DNA data, the entire southeastern United States would be regarded as a single management unit, yet the mtDNA data indicate multiple isolated populations. This complex population structure mandates a different management strategy at each life stage. Perturbations to pelagic juveniles will have a diffuse impact on Atlantic nesting colonies, mortality of subadults will have a more focused impact on nearby breeding populations, and disturbances to adults will have pinpoint impact on corresponding breeding populations. These findings demonstrate that surveys of multiple life stages are desirable to resolve management units in migratory marine species.     

Sex-biased gene flow and colonization in the Formosan lesser horseshoe bat: inference from nuclear and mitochondrial markers

Sex-biased behaviours are expected to play an important role in partitioning genetic variance in animal populations. Comparing genetic structure at markers with different modes of inheritance provides a means of detecting these behaviours and their consequences for population genetic structure. In colonially breeding mammals, the common combination of female philopatry and male vagility can promote contrasting patterns of genetic differentiation between the sexes, both via their effects on recurrent gene flow and on colonization. We examined sex differences in gene flow and structure by comparing maternally inherited mitochondrial DNA (mtDNA) and biparentally inherited autosomal loci in the Formosan lesser horseshoe bat Rhinolophus monoceros . We found that genetic partitioning was higher at mtDNA than autosomal markers in both sexes, indicative of female-biased philopatry and male-biased dispersal. Across Taiwan, isolation-by-distance was detected for all sex/marker combinations but was steeper for mtDNA than for nuclear markers. We suggest that isolation-by-distance shown from mtDNA at large scales is likely to reflect the stepwise founding of new breeding colonies by females during colonization. In contrast, no isolation-by-distance was found at smaller distances of up to 100 km, indicating that gene flow and/or recent shared ancestry homogenises genetic structure among nearby sites. Our results highlight the value of an indirect genetic approach to understanding sex-biased behaviours and their consequences in a little-studied species.     

GENETIC STRUCTURE AND MALE-MEDIATED GENE FLOW IN THE GHOST BAT (MACRODERMA GIGAS)

The Australian ghost bat is a large, opportunistic carnivorous species that has undergone a marked range contraction toward more mesic, tropical sites over the past century. Comparison of mitochondrial DNA (mtDNA) control region sequences and six nuclear microsatellite loci in 217 ghost bats from nine populations across subtropical and tropical Australia revealed strong population subdivision (mtDNA φ_ST = 0.80; microsatellites UR_ST = 0.337). Low-latitude (tropical) populations had higher heterozygosity and less marked phylogeographic structure and lower subdivision among sites within regions (within Northern Territory [NT] and within North Queensland [NQ]) than did populations at higher latitudes (subtropical sites; central Queensland [CQ]), although sampling of geographically proximal breeding sites is unavoidably restricted for the latter. Gene flow among populations within each of the northern regions appears to be male biased in that the difference in population subdivision for mtDNA and microsatellites (NT φ_ST = 0.39, UR_ST = 0.02; NQ φ_ST = 0.60, UR_ST = ?0.03) is greater than expected from differences in the effective population size of haploid versus diploid loci. The high level of population subdivision across the range of the ghost bat contrasts with evidence for high gene flow in other chiropteran species and may be due to narrow physiological tolerances and consequent limited availability of roosts for ghost bats, particularly across the subtropical and relatively arid regions. This observation is consistent with the hypothesis that the contraction of the species' range is associated with late Holocene climate change. The extreme isolation among higher-latitude populations may predispose them to additional local extinctions if the processes responsible for the range contraction continue to operate.     

Range-wide genetic population structure of common pochard (Aythya ferina): a potentially important vector of highly pathogenic avian influenza viruses

An understanding of the distribution and spatial structure of the natural vectors of zoonothic pathogens is of interest for effective disease control and prevention. Here, we investigate the range-wide population genetic structure of common pochard (Aythya ferina), a long-distance migratory duck and potential vector of highly pathogenic avian influenza. We collected several hundred samples from breeding and wintering grounds across Eurasia including some H5N1-positive individuals and generated partial sequences of the mitochondrial control region and multilocus microsatellite genotypes. Genetic differentiation among breeding populations was significant for both marker types but higher for maternally inherited mtDNA than for biparentally inherited nuclear markers. There was only weak genetic divergence between ducks sampled in Europe and East Asia, and genetic differentiation between populations was not generally associated with geographical distance. No evidence of genetic substructure was detected for ducks sampled on the European wintering grounds. Our results suggest limited breeding-site fidelity, especially in females, but extensive population admixture on the wintering grounds. The specific role of pochards as natural vectors of zoonotic pathogens and in particular H5N1 remains to be clarified but our results point to wintering grounds as potential hotspots for disease transmission.     

Identification of conservation units in the European Mergus merganser based on nuclear and mitochondrial DNA markers

The conservation status of small breeding areas of the Goosander (Mergus merganser merganser) in Central Europe is unclear. Geographic isolation of these areas suggests restricted gene flow to and from large North-European populations. On the other hand, migrating Goosanders from northern Europe join the Central European breeding population for wintering. To evaluate the conservation status of the small breeding areas we assessed the genetic structure of M. merganser populations in Europe by examining two nuclear marker systems (microsatellites and Single Nucleotide Polymorphisms, SNP) and mitochondrial (mtDNA) control region sequence variation for Goosanders in 11 sampling areas representing three of five distinct breeding areas and two subspecies (M. m. merganser and M. m. americanus). Overall population differentiation estimates including both subspecies were high, both based on mtDNA () and nuclear markers (θ _ST = 0.219; 95% CI 0.088–0.398, SNP and microsatellites combined). Within Europe, mtDNA revealed a strong overall () and significant pairwise population differentiation between almost all comparisons. In contrast, both nuclear marker systems combined revealed only a small overall genetic differentiation (θ _ST = 0.022; 95% CI 0.003–0.041). The strong genetic differentiation based on female-inherited mtDNA but not on biparentally inherited nuclear markers can be explained by sex-biased dispersal and strong female philopatry. Therefore, small breeding areas in Europe are endangered despite large male-mediated gene-flow, because when these populations decline, only males—but due to strong philopatry not females—can be efficiently supplemented by migration from the large North European populations. We therefore propose to manage the small breeding areas independently and to strengthen conservation efforts for this species in Central Europe.     

Multiscale spatial genetic structure within and between populations of wild cherry trees in nuclear genotypes and chloroplast haplotypes

Spatial genetic structure (SGS) of plants mainly depends on the effective population size and gene dispersal. Maternally inherited loci are expected to have higher genetic differentiation between populations and more intensive SGS within populations than biparentally inherited loci because of smaller effective population sizes and fewer opportunities of gene dispersal in the maternally inherited loci. We investigated biparentally inherited nuclear genotypes and maternally inherited chloroplast haplotypes of microsatellites in 17 tree populations of three wild cherry species under different conditions of tree distribution and seed dispersal. As expected, interpopulation genetic differentiation was 6–9 times higher in chloroplast haplotypes than in nuclear genotypes. This difference indicated that pollen flow 4–7 times exceeded seed flow between populations. However, no difference between nuclear and chloroplast loci was detected in within‐population SGS intensity due to their substantial variation among the populations. The SGS intensity tended to increase as trees became more aggregated, suggesting that tree aggregation biased pollen and seed dispersal distances toward shorter. The loss of effective seed dispersers, Asian black bears, did not affect the SGS intensity probably because of mitigation of the bear loss by other vertebrate dispersers and too few tree generations after the bear loss to alter SGS. The findings suggest that SGS is more variable in smaller spatial scales due to various ecological factors in local populations.     

High gene flow on a continental scale in the polyandrous Kentish plover Charadrius alexandrinus

Gene flow promotes genetic homogeneity of species in time and space. Gene flow can be modulated by sex‐biased dispersal that links population genetics to mating systems. We investigated the phylogeography of the widely distributed Kentish plover Charadrius alexandrinus. This small shorebird has a large breeding range spanning from Western Europe to Japan and exhibits an unusually flexible mating system with high female breeding dispersal. We analysed genetic structure and gene flow using a 427‐bp fragment of the mitochondrial (mtDNA) control region, 21 autosomal microsatellite markers and a Z microsatellite marker in 397 unrelated individuals from 21 locations. We found no structure or isolation‐by‐distance over the continental range. However, island populations had low genetic diversity and were moderately differentiated from mainland locations. Genetic differentiation based on autosomal markers was positively correlated with distance between mainland and each island. Comparisons of uniparentally and biparentally inherited markers were consistent with female‐biased gene flow. Maternally inherited mtDNA was less structured, whereas the Z‐chromosomal marker was more structured than autosomal microsatellites. Adult males were more related than females within genetic clusters. Taken together, our results suggest a prominent role for polyandrous females in maintaining genetic homogeneity across large geographic distances.     

Genomic scans reveal multiple mito‐nuclear incompatibilities in population crosses of the copepod Tigriopus californicus

The evolution of intrinsic postzygotic isolation can be explained by the accumulation of Dobzhansky‐Muller incompatibilities (DMI). Asymmetries in the levels of hybrid inviability and hybrid sterility are commonly observed between reciprocal crosses, a pattern that can result from the involvement of uniparentally inherited factors. The mitochondrial genome is one such factor that appears to participate in DMI in some crosses but the frequency of its involvement versus biparentally inherited factors is unclear. Here we assess the relative importance of incompatibilities between nuclear factors (nuclear‐nuclear) versus those between mitochondrial and nuclear factors (mito‐nuclear) in a species that lacks sex chromosomes. We used a Pool‐seq approach to survey three crosses among genetically divergent populations of the copepod, Tigriopus californicus, for regions of the genome that are affected by hybrid inviability. Results from reciprocal crosses suggest that mito‐nuclear incompatibilities are more common than nuclear‐nuclear incompatibilities overall. These results suggest that in the presence of very high levels of nucleotide divergence between mtDNA haplotypes, mito‐nuclear incompatibilities can be important for the evolution of intrinsic postzygotic isolation. This is particularly interesting considering this species lacks sex chromosomes, which have been shown to harbor a particularly high number of nuclear‐nuclear DMI in several other species.     

Relative information content of polymorphic microsatellites and mitochondrial DNA for inferring dispersal and population genetic structure in the olive sea snake, Aipysurus laevis

Polymorphic microsatellites are widely considered more powerful for resolving population structure than mitochondrial DNA (mtDNA) markers, particularly for recently diverged lineages or geographically proximate populations. Weaker population subdivision for biparentally inherited nuclear markers than maternally inherited mtDNA may signal male-biased dispersal but can also be attributed to marker-specific evolutionary characteristics and sampling properties. We discriminated between these competing explanations with a population genetic study on olive sea snakes, Aipysurus laevis. A previous mtDNA study revealed strong regional population structure for A. laevis around northern Australia, where Pleistocene sea-level fluctuations have influenced the genetic signatures of shallow-water marine species. Divergences among phylogroups dated to the Late Pleistocene, suggesting recent range expansions by previously isolated matrilines. Fine-scale population structure within regions was, however, poorly resolved for mtDNA. In order to improve estimates of fine-scale genetic divergence and to compare population structure between nuclear and mtDNA, 354 olive sea snakes (previously sequenced for mtDNA) were genotyped for five microsatellite loci. F statistics and Bayesian multilocus genotype clustering analyses found similar regional population structure as mtDNA and, after standardizing microsatellite F statistics for high heterozygosities, regional divergence estimates were quantitatively congruent between marker classes. Over small spatial scales, however, microsatellites recovered almost no genetic structure and standardized F statistics were orders of magnitude smaller than for mtDNA. Three tests for male-biased dispersal were not significant, suggesting that recent demographic expansions to the typically large population sizes of A. laevis have prevented microsatellites from reaching mutation-drift equilibrium and local populations may still be diverging.     

Patterns of hybridization in a multispecies hybrid zone in the Ranunculus cantoniensis complex (Ranunculaceae)

Simultaneous hybridizations among three or more parental taxa occur in nature, but these have rarely been analysed. The present study investigates the natural hybridization of diploid members of the Ranunculus cantoniensis complex over a microgeographical area in south‐western China. Sequence information from maternally inherited plastid DNA (trnQ‐rps16 and rpL32‐trnL) and biparentally inherited nuclear ribosomal DNA was used to identify hybrids, and these were further confirmed by a dramatic reduction in pollen viability. The populations from the contact zone contained individuals that had more than two nuclear ribosomal internal transcribed spacer (nrITS) types; in addition, a network analysis revealed the presence of hard conflicts between the nuclear and plastid DNA data. These results prove that R. trigonus and R. silerifolius var. dolicanthus have relatively little reproductive isolation from each other, leading to a high gene flow rate in contact populations. The four investigated species can hybridize with each other in this contact zone, forming seven hybrid types. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 174 , 227–239.     

Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China

How coniferous trees in northern China changed their distribution ranges in response to Quaternary climatic oscillations remains largely unknown. Here we report a study of the phylogeography of Pinus tabulaeformis, an endemic and dominant species of coniferous forest in northern China. We examined sequence variation of maternally inherited, seed‐dispersed mitochondrial DNA (mtDNA) (nad5 intron 1 and nad4/3–4) and paternally inherited, pollen‐ and seed‐dispersed chloroplast DNA (cpDNA) (rpl16 and trnS‐trnG) within and among 30 natural populations across the entire range of the species. Six mitotypes and five chlorotypes were recovered among 291 trees surveyed. Population divergence was high for mtDNA variation (G_ST = 0.738, N_ST = 0.771) indicating low levels of seed‐based gene flow and significant phylogeographical structure (N_ST > G_ST, P < 0.05). The spatial distribution of mitotypes suggests that five distinct population groups exist in the species: one in the west comprising seven populations, a second with a north–central distribution comprising 15 populations, a third with a southern and easterly distribution comprising five populations, a fourth comprising one central and one western population, and a fifth comprising a single population located in the north‐central part of the species’ range. Each group apart from the fourth group is characterized by a distinct mitotype, with other mitotypes, if present, occurring at low frequency. It is suggested, therefore, that most members of each group apart from Group 4 are derived from ancestors that occupied different isolated refugia in a previous period of range fragmentation of the species, possibly at the time of the Last Glacial Maximum. Possible locations for these refugia are suggested. A comparison of mitotype diversity between northern and southern subgroups within the north‐central group of populations (Group 2) showed much greater uniformity in the northern part of the range both within and between populations. This could indicate a northward migration of the species from a southern refugium in this region during the postglacial period, although alternative explanations cannot be ruled out. Two chlorotypes were distributed across the geographical range of the species, resulting in lower levels of among‐population chlorotype variation. The geographical pattern of variation for all five chlorotypes provided some indication of the species surviving past glaciations in more than one refugium, although differentiation was much less marked, presumably due to the greater dispersal of cpDNA via pollen.     

Phylogeography and recolonization of the Swiss Alps by the Valais shrew (Sorex antinorii), inferred with autosomal and sex‐specific markers

Using one male‐inherited, one female‐inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii) in the Swiss Alps. Bayesian analysis on autosomal microsatellites suggests a clear genetic differentiation between two groups of populations. This geographically based structure is consistent with two separate postglacial recolonization routes of the species into Switzerland from Italian refugia after the last Pleistocene glaciations. Sex‐specific markers also confirm genetic structuring among western and eastern areas, since very few haplotypes for either Y chromosome or mtDNA genome are shared between the two regions. Overall, these results suggest that two already well‐differentiated genetic lineages colonized the Swiss Alps and came into secondary contact in the Rhône Valley. Low level of admixture between the two lineages is likely explained by the mountainous landscape structure of lateral valleys orthogonal to the main Rhône valley.     

Natural hybridization of Yezo and Sakhalin spruce in central Hokkaido,revealed by DNA markers with contrasting modes of inheritance

Yezo spruce (Picea jezoensis var. jezoensis) and Sakhalin spruce (Picea glehnii) occur across Hokkaido and co‐occur in some forest habitats. This leads to the potential for natural hybridization between these two species, which has been shown to occur at low frequencies. The purpose of this study was to identify these hybrids and their possible mating patterns, using various Pinaceae DNA markers with different modes of inheritance. The markers used were maternally inherited mitochondrial DNA (mtDNA), paternally inherited chloroplast DNA (cpDNA) and biparentally inherited nuclear microsatellites (nSSRs). Seven putative natural hybrids, four artificially‐crossed F₁ hybrids, four parent plants from each species, and two artificially‐backcrossed hybrids of putative natural hybrids and their parents were analyzed using the diagnostic DNA markers developed in this study. We found Yezo spruce and Sakhalin spruce to be distinct (J and G types, respectively), and the modes of inheritance held true for the two species, as was previously reported to be the case in Pinaceae. Four of the seven putative natural hybrids harbored J‐type cpDNA, G‐type mtDNA and J/G‐type nSSRs, indicating that natural F₁ hybrids are likely to arise from a G (female) × J (male) crossing. One natural hybrid harbored G‐type cpDNA, J‐type mtDNA and J/G‐type nSSRs, which implies that hybrids produced by J (female) × G (male) crossings occur at low frequencies. The two remaining hybrids harbored J‐type cpDNA and mtDNA with either J/G or J/J‐type nSSRs, suggesting that they may be F₂ hybrids resulting from backcrossing between an F₁ hybrid and a Yezo spruce.     

Male-biased dispersal promotes large scale gene flow in a subterranean army ant, Dorylus (Typhlopone) fulvus

Sex-biased dispersal is a widespread phenomenon in the animal kingdom, which strongly influences gene flow and population structure. Particularly army ants, important key-stone predators in tropical ecosystems, are prone to population fragmentation and isolation due to their extraordinary mating system: queens are permanently wingless, propagate via colony fission, and only the males disperse in mating flights. Here we report on sex-biased dispersal and the genetic population structure of an African subterranean army ant, Dorylus (Typhlopone) fulvus. Using maternally inherited mtDNA markers and bi-parentally inherited nuclear microsatellites we found strong geographical structuring of mtDNA haplotypes, whereas the nuclear genetic population structure was less pronounced. Strong mtDNA (Φ _ST = 0.85), but significantly lower nuclear (F _ST = 0.23) genetic differentiation translated into a more than an order of magnitude larger male migration rate compared to that of queens, reflecting the low motility of queens and strong, promiscuous dispersal by males. Thus, the well flying D. fulvus males appear to be the sex to promote large scale gene flow, and D. fulvus is indeed a species in which sex specific dispersal patterns and the mating system profoundly affect the population structure and phylogeography.     

A re-evaluation of management units based on gene flow of a rare waterbird in the Americas

The maintenance of gene flow in species that have experienced population contractions and are geographically fragmented is important to the maintenance of genetic variation and evolutionary potential; thus, gene flow is also important to conservation and management of these species. For example, the Reddish Egret (Egretta rufescens) has recovered after severe population reductions during the 19th and 20th centuries, but population numbers remain below historical levels. In this study, we characterized gene flow among management units of the Reddish Egret by using ten nuclear microsatellite markers and part of the mitochondrial (mtDNA) control region from 176 nestlings captured at eight localities in Mexico (Baja California, Chiapas, Tamaulipas, and Yucatan), the USA (Texas, Louisiana, and Florida), and the Bahamas. We found evidence of population structure and that males disperse more often and across longer distances compared with females, which is congruent with previous banding and telemetry data. The maternally inherited mtDNA and biparentally inherited microsatellite data supported slightly different MU models; however, when interpreted together, a four MU model that considered population structure and geographic proximity was most optimal. Namely, MU 1 (Baja California); MU 2 (Chiapas); MU 3 (Yucatan, Tamaulipas, Texas, and Louisiana); and MU 4 (Florida and the Bahamas). Regions outside our sampled localities (e.g., the Greater Antilles and South America) require additional sampling to fully understand gene flow and movement of individuals across the species’ entire range. However, the four MUs we have defined group nesting localities into genetically similar subpopulations, which can guide future management plans.     

DNA‐based approaches for evaluating historical demography in terrestrial vertebrates

Contemporary DNA sequences can provide information about the historical demography of a species. However, different molecular markers are informative under different circumstances. In particular, mitochondrial (mt)DNA is uniparentally inherited and haploid in most vertebrates and thus has a smaller effective population size than diploid, biparentally inherited nuclear (n)DNA. Here, we review the characteristics of mtDNA and nDNA in the context of historical demography. In particular, we address how their contrasting rates of evolution and sex‐biased dispersal can lead to different demographic inferences. We do so in the context of an extensive review of the vertebrate literature that describes the use of mtDNA and nDNA sequence data in demographic reconstruction. We discuss the effects of coalescence, effective population size, substitution rates, and sex‐biased dispersal on informative timeframes and expected patterns of genetic differentiation. We argue that mtDNA variationin species with male‐biased dispersal can imply deviations from neutrality that do not reflect actual population expansion or selection. By contrast, mtDNA can be more informative when coalescence has occurred within the recent past, which appears to be the case with many vertebrates. We also compare the application and interpretation of demographic and neutrality test statistics in historical demography studies. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 112 , 367–386.     

Differential introgression causes genealogical discordance in host races of Acrocercops transecta (Insecta: Lepidoptera)

Recently diverged populations often exhibit incomplete reproductive isolation, with a low level of gene flow continuing between populations. Previous studies have shown that, even under a low level of gene flow, genetic divergence between populations can proceed at the loci governing local adaptation and reproductive isolation but not at other neutral loci. A leaf‐mining moth, Acrocercops transecta, consists of Juglans‐ and Lyonia‐associated host races. The two host races differ in host preferences of ovipositing females and in larval adaptation to host plants but mate readily in the laboratory, producing fertile hybrids. The Juglans and Lyonia races are often sympatric in the wild, implying that gene introgression could occur in nature between the two host races. We tested this hypothesis by combining phylogenetic analyses with coalescent simulations, focusing on mitochondrial genes (COI and ND5) and the nuclear Tpi, Per and Ldh genes located on the Z‐chromosome. The mitochondrial genes clearly distinguished the Lyonia race from the Juglnas race, whereas the Tpi, Per and Ldh genealogies did not reflect the two host races. Coalescent simulations indicated gene flow at the three Z‐linked genes in both directions, whereas there was no introgression in the mitochondrial genes. The lack of introgression in mitochondrial genes suggests that female host preference is the primary force leading to the bifurcation of maternally inherited loci. Thus, the results show that a low level of gene flow coupled with the inflexible female host preference differentiates histories of divergence between maternally and biparentally inherited genes in this host race system.     

Population structure of nuclear and mitochondrial DNA variation among humpback whales in the North Pacific

The population structure of variation in a nuclear actin intron and the control region of mitochondrial DNA is described for humpback whales from eight regions in the North Pacific Ocean: central California, Baja Peninsula, nearshore Mexico (Bahia Banderas), offshore Mexico (Socorro Island), southeastern Alaska, central Alaska (Prince Williams Sound), Hawaii and Japan (Ogasawara Islands). Primary mtDNA haplotypes and intron alleles were identified using selected restriction fragment length polymorphisms of target sequences amplified by the polymerase chain reaction (PCR–RFLP). There was little evidence of heterogeneity in the frequencies of mtDNA haplotypes or actin intron alleles due to the year or sex composition of the sample. However, frequencies of four mtDNA haplotypes showed marked regional differences in their distributions (Φ_ST = 0.277; P < 0.001; n = 205 individuals) while the two alleles showed significant, but less marked, regional differences (Φ_ST = 0.033; P < 0.013; n = 400 chromosomes). An hierarchical analysis of variance in frequencies of haplotypes and alleles supported the grouping of six regions into a central and eastern stock with further partitioning of variance among regions within stocks for haplotypes but not for alleles. Based on available genetic and demographic evidence, the southeastern Alaska and central California feeding grounds were selected for additional analyses of nuclear differentiation using allelic variation at four microsatellite loci. All four loci showed significant differences in allele frequencies (overall F_ST = 0.043; P < 0.001; average n = 139 chromosomes per locus), indicating at least partial reproductive isolation between the two regions as well as the segregation of mtDNA lineages. Although the two feeding grounds were not panmictic for nuclear or mitochondrial loci, estimates of long-term migration rates suggested that male-mediated gene flow was several-fold greater than female gene flow. These results include and extend the range and sample size of previously published work, providing additional evidence for the significance of genetic management units within oceanic populations of humpback whales.