The influence of mitonuclear genetic variation on personality in seed beetles

There is a growing awareness of the influence of mitochondrial genetic variation on life-history phenotypes, particularly via epistatic interactions with nuclear genes. Owing to their direct effect on traits such as metabolic and growth rates, mitonuclear interactions may also affect variation in behavioural types or personalities (i.e. behavioural variation that is consistent within individuals, but differs among individuals). However, this possibility is largely unexplored. We used mitonuclear introgression lines, where three mitochondrial genomes were introgressed into three nuclear genetic backgrounds, to disentangle genetic effects on behavioural variation in a seed beetle. We found within-individual consistency in a suite of activity-related behaviours, providing evidence for variation in personality. Composite measures of overall activity of individuals in behavioural assays were influenced by both nuclear genetic variation and by the interaction between nuclear and mitochondrial genomes. More importantly, the degree of expression of behavioural and life-history phenotypes was correlated and mitonuclear genetic variation affected expression of these concerted phenotypes. These results show that mitonuclear genetic variation affects both behavioural and life-history traits, and they provide novel insights into the maintenance of genetic variation in behaviour and personality.     

Integrative genomic phylogeography reveals signs of mitonuclear incompatibility in a natural hybrid goby population

Hybridization between divergent lineages generates new allelic combinations. One mechanism that can hinder the formation of hybrid populations is mitonuclear incompatibility, that is, dysfunctional interactions between proteins encoded in the nuclear and mitochondrial genomes (mitogenomes) of diverged lineages. Theoretically, selective pressure due to mitonuclear incompatibility can affect genotypes in a hybrid population in which nuclear genomes and mitogenomes from divergent lineages admix. To directly and thoroughly observe this key process, we de novo sequenced the 747‐Mb genome of the coastal goby, Chaenogobius annularis, and investigated its integrative genomic phylogeographics using RNA‐sequencing, RAD‐sequencing, genome resequencing, whole mitogenome sequencing, amplicon sequencing, and small RNA‐sequencing. Chaenogobius annularis populations have been geographically separated into Pacific Ocean (PO) and Sea of Japan (SJ) lineages by past isolation events around the Japanese archipelago. Despite the divergence history and potential mitonuclear incompatibility between these lineages, the mitogenomes of the PO and SJ lineages have coexisted for generations in a hybrid population on the Sanriku Coast. Our analyses revealed accumulation of nonsynonymous substitutions in the PO‐lineage mitogenomes, including two convergent substitutions, as well as signals of mitochondrial lineage‐specific selection on mitochondria‐related nuclear genes. Finally, our data implied that a microRNA gene was involved in resolving mitonuclear incompatibility. Our integrative genomic phylogeographic approach revealed that mitonuclear incompatibility can affect genome evolution in a natural hybrid population.     

Gut microbial diversity across a contact zone for California voles: Implications for lineage divergence of hosts and mitonuclear mismatch in the assembly of the mammalian gut microbiome

Gut microbial diversity is thought to reflect the co‐evolution of microbes and their hosts as well as current host‐specific attributes such as genetic background and environmental setting. To explore interactions among these parameters, we characterized variation in gut microbiome composition of California voles (Microtus californicus) across a contact zone between two recently diverged lineages of this species. Because this contact zone contains individuals with mismatched mitochondrial‐nuclear genomes (cybrids), it provides an important opportunity to explore how different components of the genotype contribute to gut microbial diversity. Analyses of bacterial 16S rRNA sequences and joint species distribution modelling revealed that host genotypes and genetic differentiation among host populations together explained more than 50% of microbial community variation across our sampling transect. The ranked importance (most to least) of factors contributing to gut microbial diversity in our study populations were: genome‐wide population differentiation, local environmental conditions, and host genotypes. However, differences in microbial communities among vole populations (β‐diversity) did not follow patterns of lineage divergence (i.e., phylosymbiosis). Instead, among‐population variation was best explained by the spatial distribution of hosts, as expected if the environment is a primary source of gut microbial diversity (i.e., dispersal limitation hypothesis). Across the contact zone, several bacterial taxa differed in relative abundance between the two parental lineages as well as among individuals with mismatched mitochondrial and nuclear genomes. Thus, genetic divergence among host lineages and mitonuclear genomic mismatches may also contribute to microbial diversity by altering interactions between host genomes and gut microbiota (i.e., hologenome speciation hypothesis).     

Widespread co‐occurrence of divergent mitochondrial haplotype lineages in a Central American species of poison frog (Oophaga pumilio)

Aim To analyse the phylogeographic structure of the strawberry poison frog, Oophaga pumilio (Dendrobatidae), across a large part of its range using a combination of mitochondrial and nuclear markers. Location Costa Rica and Panama. Methods Sequence analyses of a mitochondrial (cytochrome b) and a nuclear (RAG‐1) gene fragment as well as analyses of seven microsatellite loci were carried out on 269 individuals of O. pumilio sampled from 24 localities and on two individuals of O. vicentei. Results Two main mitochondrial haplotype lineages, corresponding to a northern (north Costa Rica) and a southern (south Costa Rica and eastern Panama) lineage, were identified. They differed by up to 7% uncorrected distance. We observed co‐occurrence of both lineages in seven populations in Costa Rica and Panama, indicating a pattern of extensive admixture. The two main mitochondrial lineages of O. pumilio roughly corresponded to a previously described phylogeographic pattern. Microsatellites indicate admixture spanning over a wide geographic area, but significant variation between the northern and southern groups was also found with microsatellite data. While microsatellite data reconstructed a separation south of an assumed Caribbean valley barrier, mitochondrial haplotypes of the ‘southern lineage’ shifted this barrier towards the north. Main conclusions Despite admixture, all three markers showed significant variation between the northern and southern groups. Phylogeographical breaks known from other anuran species in the study region could not be verified for O. pumilio. The unexpected clustering of the population from Escudo de Veraguas and the individuals of O. vincentei with the northern O. pumilio lineage indicates the need for a fundamental and careful taxonomic revision, including an interspecific phylogeography of the entire genus.     

Direct and indirect genetic effects of sex-specific mitonuclear epistasis on reproductive ageing

Mitochondria are involved in ageing and their function requires coordinated action of both mitochondrial and nuclear genes. Epistasis between the two genomes can influence lifespan but whether this also holds for reproductive senescence is unclear. Maternal inheritance of mitochondria predicts sex differences in the efficacy of selection on mitonuclear genotypes that should result in differences between females and males in mitochondrial genetic effects. Mitonuclear genotype of a focal individual may also indirectly affect trait expression in the mating partner. We tested these predictions in the seed beetle Callosobruchus maculatus, using introgression lines harbouring distinct mitonuclear genotypes. Our results reveal both direct and indirect sex-specific effects of mitonuclear epistasis on reproductive ageing. Females harbouring coadapted mitonuclear genotypes showed higher lifetime fecundity due to slower senescence relative to novel mitonuclear combinations. We found no evidence for mitonuclear coadaptation in males. Mitonuclear epistasis not only affected age-specific ejaculate weight, but also influenced male age-dependent indirect effects on traits expressed by their female partners (fecundity, egg size, longevity). These results demonstrate important consequences of sex-specific mitonuclear epistasis for both mating partners, consistent with a role for mitonuclear genetic constraints upon sex-specific adaptive evolution.     

GENETIC AND MORPHOLOGICAL ANALYSES OF THE SOUTHERN BULL KELP DURVILLAEA ANTARCTICA (PHAEOPHYCEAE: DURVILLAEALES) IN NEW ZEALAND REVEAL CRYPTIC SPECIES1

Many macroalgae exhibit considerable intraspecific morphological variation, but whether such variation reflects phenotypic plasticity or underlying genetic differences is often poorly understood. We quantified both morphological and genetic variation of 96 plants from seven field sites across eastern South Island, New Zealand, to assess genetic differences between morphotypes of the southern bull kelp Durvillaea antarctica (Cham.) Har. Consistent DNA sequence differentiation across mitochondrial, plastid, and nuclear loci was correlated with two broadly sympatric morphotypes: “cape” and “thonged.” These ecologically, morphologically, and genetically distinct bull‐kelp lineages were previously considered to be environmentally determined phenotypes with no underlying genetic basis. Interestingly, the sheltered “cape” lineage appears essentially genetically uniform across its South Island range, whereas the exposed “thonged” lineage exhibits marked phylogeographic structure across its range. Results suggest that D. antarctica in New Zealand comprises two reproductively isolated species.     

Genetic rescue of an inbred captive population of the critically endangered Puerto Rican crested toad (<Emphasis Type="Italic">Peltophryne lemur</Emphasis>) by mixing lineages

The Puerto Rican crested toad (Peltophryne lemur) is currently composed of a single wild population on the south coast of Puerto Rico and two captive populations founded by animals from the northern and southern coasts. The main factors contributing to its decline are habitat loss, inundation of breeding ponds during storms, and impacts of invasive species. Recovery efforts have been extensive, involving captive breeding and reintroductions, habitat restoration, construction of breeding ponds, and public education. To guide future conservation efforts, genetic variation and differentiation were assessed for the two captive colonies and the remaining wild population using the mitochondrial control region and six novel microsatellite loci. Only two moderately divergent mitochondrial haplotypes were found, with one fixed in each of the southern and northern lineages. Moderate genetic variation exists for microsatellite loci in all three groups. The captive southern population has not diverged substantially from the wild population at microsatellite loci (F _ST = 0.03), whereas there is little allelic overlap between the northern and southern lineages at five of six loci (F _ST > 0.3). Despite this differentiation, they are no more divergent than many populations of other amphibian species. As the northern breeding colony may not remain viable due to its small size and inbred nature, it is recommended that a third breeding colony be established in which northern and southern individuals are combined. This will preserve any northern adaptive traits that may exist, and provide animals for release in the event that the pure northern lineage becomes extirpated.     

Exploring the potential of nuclear and mitochondrial sequencing data generated through genome‐skimming for plant phylogenetics: A case study from a clade of neotropical lianas

Few botanical studies have explored the potential of nuclear ribosomal DNA (nrDNA) and mitochondrial DNA (mtDNA) data obtained through genome skimming for phylogeny reconstruction. Here, we analyzed the phylogenetic information included in the nrDNA and mtDNA of 44 species of the “Adenocalymma‐Neojobertia” clade (Bignoniaceae). To deal with intraindividual polymorphisms within the nrDNA, different coding schemes were explored through the analyses of four datasets: (i) “nrDNA contig,” with base call following the majority rule; (ii) “nrDNA ambiguous,” with ambiguous base calls; (iii) “nrDNA informative,” with ambiguities converted to multistate characters; and, (iv) “mitochondrial,” with 39 mitochondrial genes. Combined analyses using the nrDNA and mtDNA data and previously published “plastid” datasets were also conducted. Trees were obtained using Maximum Likelihood and Bayesian criteria. The congruence among genomes was assessed. The nrDNA datasets were shown to be highly polymorphic within individuals, while the “mitochondrial” dataset was the least informative, with 0.36% of informative bases within the ingroup. The topologies inferred using the nrDNA and mtDNA datasets were broadly congruent with the tree derived from the analyses of the “plastid” dataset. The topological differences recovered were generally poorly supported. The topology that resulted from the analyses of the “combined” dataset largely resembles the “plastid” tree. These results highlight limitations of nuclear ribosomal DNA and mitochondrial genes for phylogeny reconstruction obtained through genome skimming and the need to include more data from both genomes. The different topologies observed among genomes also highlight the importance of exploring data from various genomes in plant phylogenetics.     

Genetic Variation in Wild and Hatchery Stocks of Suminoe Oyster (Crassostrea ariakensis) Assessed by PCR-RFLP and Microsatellite Markers

Genetic variation in wild Asian populations and U.S. hatchery stocks of Crassostrea ariakensis was examined using polymerase chain reactions with restriction fragment length polymorphism (PCR-RFLP) analysis of both the mitochondrial COI gene and the nuclear internal transcribed spacer (ITS) 1 region and using 3 microsatellite markers. Hierarchical analysis of molecular variance and pairwise comparisons revealed significant differentiation (P < 0.05) between samples from the northern region, represented by collections from China and Japan, and 2 of 3 samples from southern China. PCR-RFLP patterns were identified that were diagnostic for the northern (N-type) and southern (S-type) groups. Microsatellite marker profiles were used to assign each oyster to one of the two northern or two southern populations. Results for more than 97% of the oysters were consistent with the PCR-RFLP patterns observed for each individual in that oysters with N-type patterns were assigned to one of the northern populations and those with S-type patterns to one of the southern populations. At one site of the Beihai (B) region in southern China a mix of individuals with either the N-type or S-type PCR-RFLP genotypes was found. No heterozygotes at the nuclear ITS-1 locus were found in the sample, possibly indicating reproductive isolation in sympatry. Microsatellite assignment test results of the B individuals were also consistent with identifications as either the N-type or S-type based on PCR-RFLP patterns. The parental population for one hatchery stock was this B sample, which initially was composed of almost equal numbers of northern and southern genetic types. After hatchery spawns, however, more than 97% of the progeny fell into the northern genetic group by PCR-RFLP and microsatellite assignment test analyses, indicating that the individuals with the southern genotype contributed little to the spawn, owing to gametic incompatibility, differential larval survival, or a difference in timing of sexual maturity. Overall, results suggested that oysters collected as C. ariakensis in this study, and likely in other studies as well, include two different sympatric species with some degree of reproductive isolation.     

Historic hybridization and persistence of a novel mito-nuclear combination in red-backed voles (genus Myodes)

Background  

The role of hybridization in generating diversity in animals is an active area of discovery and debate. We assess hybridization across a contact zone of northern (Myodes rutilus) and southern (M. gapperi) red-backed voles using variation in skeletal features and both mitochondrial and nuclear loci. This transect extends approximately 550 km along the North Pacific Coast of North America and encompasses 26 populations (n = 485). We establish the history, geographic extent and directionality of hybridization, determine whether hybridization is ongoing, and assess the evolutionary stability of novel genomic combinations.     

Identification, distribution and molecular evolution of the pacifastin gene family in Metazoa

Background

This study was motivated by the observation of unusual mitochondrial haplotype distributions and associated physiological differences between populations of the killifish Fundulus heteroclitus distributed along the Atlantic coast of North America. A distinct "northern" haplotype is fixed in all populations north of New Jersey, and does not appear south of New Jersey except in extreme upper-estuary fresh water habitats, and northern individuals are known to be more tolerant of hyposmotic conditions than southern individuals. Complete mitochondrial genomes were sequenced from individuals from northern coastal, southern coastal, and fresh water populations (and from out-groups). Comparative genomics approaches were used to test multiple evolutionary hypotheses proposed to explain among-population genome variation including directional selection and hybridization.

Results

Structure and organization of the Fundulus mitochondrial genome is typical of animals, yet subtle differences in substitution patterns exist among populations. No signals of directional selection or hybridization were detected. Mitochondrial genes evolve at variable rates, but all genes exhibit very low dN/dS ratios across all lineages, and the southern population harbors more synonymous polymorphism than other populations.

Conclusion

Evolution of mitochondrial genomes within Fundulus is primarily governed by interaction between strong purifying selection and demographic influences, including larger historical population size in the south. Though directional selection and hybridization hypotheses were not supported, adaptive processes may indirectly contribute to partitioning of variation between populations.     

Hybrid zones between two genetically differentiated forms of the pond frog Rana lessonae in southern Italy

Two distinct population groups of the pond frog Rana lessonae were detected in peninsular Italy and Sicily by multilocus electrophoresis: one group inhabits the peninsula down to northern Calabria, the second occurs in southern Calabria and on Sicily. Fixed alternative alleles distinguish the two groups at 5 of the 25 loci examined; marked allele frequency differences were observed at two additional loci. On average, the two groups differ by a Nei's standard genetic distance of 0.4. A wide hybrid zone (about 120 km) occurs between the two groups, with high genotypic diversity and absence of pure parental genotypes in central Calabria. Patterns of allozyme variation suggest that at least two distinct contact and hybridization events occurred, one in the Catanzaro, the other in the Crati-Sibari plains, about 70 km to the north. Geological evidence indicates that these areas correspond to two main marine-flooded grabens that would have repeatedly interrupted or reduced genetic exchange during Plio-Pleistocene times. The finding of a fixed difference at the Mdhp-1 locus between Sicilian and Calabrian R. lessonae witnesses their continuing differentiation following their last separation by definitive opening of the Strait of Messina, about 50 000 years ago. The wide hybrid zone, the diversity of genotypes and the agreement with Hardy-Weinberg expectations suggest complete hybrid fertility. Different patterns of introgression were observed at the various loci. The pattern of allelic variation at loci in R. lessonae is paralleled by the pattern of variation in lessonae genomes of the sympatric hemiclonal hybrid Rana esculenta, into which lessonae genomes are introduced each generation as a result of hybridogenesis.     

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.     

THE INTERFACE BETWEEN PHYLOGENETICS AND POPULATION GENETICS: INVESTIGATING GENE TREES,SPECIES TREES,AND POPULATION DYNAMICS IN THE PHYLLOPHAGA FRATERNA SPECIES GROUP

This study uses traditional and contemporary phylogenetic and population genetic analyses to assess the causes of discordance (i.e., lineage sorting and introgression) among mitochondrial and nuclear gene trees for a clade of eastern North American scarab beetles (fraterna species group, genus Phyllophaga). I estimated gene trees using individual and combined analysis of one mitochondrial and two nuclear loci in MrBayes , and inferred a species tree using a hierarchical coalescent approach based on all loci in the program Best . Because hybridization violates the assumptions of Best , I tested for introgression by comparing species monophyly between the mitochondrial and nuclear gene trees based on the prediction that cytoplasmic genomes introgress more readily than nuclear genomes. Haplotype exclusivity was identified using Bayesian tests of monophyly and the genealogical sorting index. I used the results of the phylogenetic analyses and monophyly tests to develop an explicit hypothesis of introgression that could be tested in the program IMa. Results from these analyses provided evidence for introgression across clades within the fraterna group. The tiered analytical approach used in this study demonstrated how the use of multiple methods can identify when assumptions are violated and methods are prone to yield misleading results.     

Mitonuclear coevolution as the genesis of speciation and the mitochondrial DNA barcode gap

Mitochondrial genes are widely used in taxonomy and systematics because high mutation rates lead to rapid sequence divergence and because such changes have long been assumed to be neutral with respect to function. In particular, the nucleotide sequence of the mitochondrial gene cytochrome c oxidase subunit 1 has been established as a highly effective DNA barcode for diagnosing the species boundaries of animals. Rarely considered in discussions of mitochondrial evolution in the context of systematics, speciation, or DNA barcodes, however, is the genomic architecture of the eukaryotes: Mitochondrial and nuclear genes must function in tight coordination to produce the complexes of the electron transport chain and enable cellular respiration. Coadaptation of these interacting gene products is essential for organism function. I extend the hypothesis that mitonuclear interactions are integral to the process of speciation. To maintain mitonuclear coadaptation, nuclear genes, which code for proteins in mitochondria that cofunction with the products of mitochondrial genes, must coevolve with rapidly changing mitochondrial genes. Mitonuclear coevolution in isolated populations leads to speciation because population‐specific mitonuclear coadaptations create between‐population mitonuclear incompatibilities and hence barriers to gene flow between populations. In addition, selection for adaptive divergence of products of mitochondrial genes, particularly in response to climate or altitude, can lead to rapid fixation of novel mitochondrial genotypes between populations and consequently to disruption in gene flow between populations as the initiating step in animal speciation. By this model, the defining characteristic of a metazoan species is a coadapted mitonuclear genotype that is incompatible with the coadapted mitochondrial and nuclear genotype of any other population.     

Mitonuclear linkage disequilibrium in human populations

There is extensive evidence from model systems that disrupting associations between co-adapted mitochondrial and nuclear genotypes can lead to deleterious and even lethal consequences. While it is tempting to extrapolate from these observations and make inferences about the human-health effects of altering mitonuclear associations, the importance of such associations may vary greatly among species, depending on population genetics, demographic history and other factors. Remarkably, despite the extensive study of human population genetics, the statistical associations between nuclear and mitochondrial alleles remain largely uninvestigated. We analysed published population genomic data to test for signatures of historical selection to maintain mitonuclear associations, particularly those involving nuclear genes that encode mitochondrial-localized proteins (N-mt genes). We found that significant mitonuclear linkage disequilibrium (LD) exists throughout the human genome, but these associations were generally weak, which is consistent with the paucity of population genetic structure in humans. Although mitonuclear LD varied among genomic regions (with especially high levels on the X chromosome), N-mt genes were statistically indistinguishable from background levels, suggesting that selection on mitonuclear epistasis has not preferentially maintained associations involving this set of loci at a species-wide level. We discuss these findings in the context of the ongoing debate over mitochondrial replacement therapy.     

Phylogeography of an Italian endemic salamander (genus Salamandrina): glacial refugia,postglacial expansions,and secondary contact

The Italian endemic genus Salamandrina has been historically regarded as monotypic but, recently, studies based on both mitochondrial and nuclear markers have indicated the existence of two distinct species of spectacled salamanders: Salamandrina perspicillata, in central and northern Italy, and Salamandrina terdigitata, in southern Italy. We analyzed nucleotide variation at mitochondrial and nuclear genes [cytochrome b, 12S and 16S rRNA, recombination activating gene (RAG 1)] in 223 individuals from 56 locations, aiming to investigate their genetic structure and recent evolutionary histories. Phylogenetic and phylogeographical analyses revealed the existence of three and two genetically distinct groups of populations in northern and southern salamander, respectively. Historical demographic analyses led to the inference of range expansion for both species in the late Pleistocene. During the last glacial stage, each salamander survived in a single refugium, namely the southern in Calabria and the northern in central Italy. At the end of this period, both lineages expanded northward and established secondary contact. Spatial distribution of RAG 1 haplotype variation revealed two differentiated population groups corresponding to the major mitochondrial (mt)DNA clades. Nuclear pattern of introgressive hybridization was more extensive than the highly limited introgression of mtDNA markers. From a conservation standpoint, southern Latium and Calabria proved to be the major genetic diversity reservoirs, thus deserving particular conservation efforts. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 903–922.     

RIBOSOMAL-DNA,MITOCHONDRIAL-DNA,CHROMOSOMAL, AND ALLOZYMIC STUDIES ON A CONTACT ZONE IN THE POCKET GOPHER,GEOMYS

We studied 75 individuals of the Plains pocket gopher, Geomys bursarius, from eastern New Mexico, where the subspecies major and knoxjonesi hybridize. Each individual was examined for chromosome number, ribosomal DNA, mitochondrial DNA, and three protein systems for which reference parental populations were fixed for alternative alleles. Twenty individuals were indistinguishable from parental major, 14 individuals were indistinguishable from parental knoxjonesi, and 41 individuals had genotypes composed of combinations of character states that distinguish the two parental types. The parental types appear to represent discrete genetic entities that have restricted introgression across a narrow hybrid zone (width approximately 3 km, using the 20/80 criterion). Parental types overlap in geographic distribution near the center of the zone, and changes in mitochondrial DNA and the five nuclear markers are concordant across the zone. It is probable that there is premating isolation between knoxjonesi males and major females. The frequencies of individuals with certain genotypic combinations within our sample imply differential reproductive success of certain genotypes. We propose that F₁'s and highly heterozygous males are sterile and that hybrid females are less fertile than parental females. These postmating factors, along with premating isolation for one of the reciprocal crosses, probably account for the restriction of gene flow across the contact zone. The structure of the zone can be explained by the “dynamic equilibrium” model.