Hybridization of Bombina bombina and B. variegata (Anura, Discoglossidae) at a sharp ecotone in western Ukraine: comparisons across transects and over time

Abstract Bombina bombina and B. variegata are two anciently diverged toad taxa that have adapted to different breeding habitats yet hybridize freely in zones of overlap where their parapatric distributions meet. Here, we report on a joint genetic and ecological analysis of a hybrid zone in the vicinity of Stryi in western Ukraine. We used five unlinked allozyme loci, two nuclear single nucleotide polymorphisms and a mitochondrial DNA haplotype as genetic markers. Parallel allele frequency clines with a sharp central step occur across a sharp ecotone, where transitions in aquatic habitat, elevation, and terrestrial vegetation coincide. The width of the hybrid zone, estimated as the inverse of the maximum gradient in allele frequency, is 2.3 km. This is the smallest of four estimates derived from different clinal transects across Europe. We argue that the narrow cline near Stryi is mainly due to a combination of habitat distribution and habitat preference. Adult toads show a preference for either ponds ( B. bombina ) or puddles ( B. variegata ), which is known to affect the distribution of genotypes within the hybrid zones. At Stryi, it should cause a reduction of the dispersal rate across the ecotone and thus narrow the cline. A detailed comparison of all five intensively studied Bombina transects lends support to the hypothesis that habitat distribution plus habitat preference can jointly affect the structure of hybrid zones and, ultimately, the resulting barriers to gene flow between differentiated gene pools. This study also represents a resampling of an area that was last studied more than 70 years ago. Our allele-frequency clines largely coincide with those that were described then on the basis of morphological variation. However, we found asymmetrical introgression of B. variegata genes into B. bombina territory along the bank of a river.     

Mitochondrial DNA variation in the hybridizing fire-bellied toads, Bombina bombina and B. variegata

Using five restriction enzymes, geographical variation of mitochondrial DNA (mtDNA) in Bombina bombina and B. variegata was studied in samples from 20 locations. Each restriction enzyme produced a species-specific fragment pattern. B. bombina haplotypes A and B were closely related to each other. In contrast, haplotypes A and B of B. variegata formed two distinct lineages. A very distinctive haplotype (C) was found in the Carpathian Mountains, whereas two other haplotypes, D and E (differing by a single AvaI site), were present in western Europe and the Balkans, respectively. Populations polymorphic for haplotypes D and E occurred in the central Balkans where the haplotypes could replace each other clinally. mtDNA sequence divergence between B. bombina and B. variegata was estimated as 6.0-8.1% and 4.7-5.2% between type C and types D/E of B. variegata. The latter divergence is contrary to allozyme and morphological data that place the western and Carpathian B. v. variegata together (Nei's D = 0.07) and separate them from the Balkan subspecies B. v. scabra (Nei's D = 0.18). Broad interspecific correlation among morphology, allozymes and mtDNA types in European fire-bellied toads argues that, despite continuous hybridization (interrupted perhaps during Pleistocene glacial maxima), little or no mtDNA introgression between the species has occurred outside the narrow hybrid zones that separate these parapatric species.     

Comparative Analysis of Mitochondrial Genomes in <Emphasis Type="Italic">Bombina</Emphasis> (Anura; Bombinatoridae)

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.     

Mapping Bombina mitochondrial genomes: the conundrum of Carpathian Bombina variegata (Anura: Discoglossidae)

Complete mitochondrial genomes (mtDNA) of five individuals representing two haplotypes of Bombina bombina and three of Bombina variegata were compared using restriction site maps. Phylogenetic analyses reveal three ancient mitochondrial lineages: (1) two very similar haplotypes A and B of B. bombina ; (2) almost identical haplotypes D and E of B. variegata ; and (3) haplotype C of B. variegata . Haplotype C is as different from haplotypes D/E as from A/B. These data are strikingly discordant with relationships based on morphology and allozymes. Haplotypes C and D/E represent a pre-Pleistocene mitochondrial divergence within B. variegata , nearly coincident with speciation between B. variegata and B. bombina . Geographical partitioning of the two divergent B. variegata mitochondrial lineages indicates repeated localization of the lineages in separate glacial refugia during the Pleistocene. That nuclear genes do not show a similar divergence, but rather indicate relatively free genetic exchange between populations with divergent mtDNAs, suggests that males dispersed much more widely than females during expansions from glacial refugia. Comparison of Bombina mtDNA maps with a restriction site map of Xenopus laevis mtDNA revealed 16 homologous sites; 12 of these may be nearly invariant across primitive anuran mtDNAs. Two distinct regions of heteroplasmy, representing two regions with variable numbers of sequence repeats [length variable (LV) regions], were characterized. Comparison with the Xenopus map places LV1, present in all five haplotypes, near the 5'-end of the control region, and LV2 present only in B. variegata , near the 3'-end. Although phylogenetic analyses did not group the two major B. variegata lineages together, presence of LV2 in both lineages supports placement of both within B. variegata .     

Phylogeography of the fire-bellied toads Bombina: independent Pleistocene histories inferred from mitochondrial genomes

The fire-bellied toads Bombina bombina and Bombina variegata, interbreed in a long, narrow zone maintained by a balance between selection and dispersal. Hybridization takes place between local, genetically differentiated groups. To quantify divergence between these groups and reconstruct their history and demography, we analysed nucleotide variation at the mitochondrial cytochrome b gene (1096 bp) in 364 individuals from 156 sites representing the entire range of both species. Three distinct clades with high sequence divergence (K2P = 8-11%) were distinguished. One clade grouped B. bombina haplotypes; the two other clades grouped B. variegata haplotypes. One B. variegata clade included only Carpathian individuals; the other represented B. variegata from the southwestern parts of its distribution: Southern and Western Europe (Balkano-Western lineage), Apennines, and the Rhodope Mountains. Differentiation between the Carpathian and Balkano-Western lineages, K2P approximately 8%, approached interspecific divergence. Deep divergence among European Bombina lineages suggests their preglacial origin, and implies long and largely independent evolutionary histories of the species. Multiple glacial refugia were identified in the lowlands adjoining the Black Sea, in the Carpathians, in the Balkans, and in the Apennines. The results of the nested clade and demographic analyses suggest drastic reductions of population sizes during the last glacial period, and significant demographic growth related to postglacial colonization. Inferred history, supported by fossil evidence, demonstrates that Bombina ranges underwent repeated contractions and expansions. Geographical concordance between morphology, allozymes, and mtDNA shows that previous episodes of interspecific hybridization have left no detectable mtDNA introgression. Either the admixed populations went extinct, or selection against hybrids hindered mtDNA gene flow in ancient hybrid zones.     

Inferences of selection and migration in the Danish house mouse hybrid zone

We analysed the patterns of allele frequency change for ten diagnostic autosomal allozyme loci in the hybrid zone between the house mouse subspecies Mus musculus domesticus and M. m. musculus in central Jutland. After determining the general orientation of the clines of allele frequencies, we analysed the cline shapes along the direction of maximum gradient. Eight of the ten clines are best described by steep central steps with coincident positions and an average width of 8.9 km (support limits 7.6–12.4) flanked by tails of introgression, indicating the existence of a barrier to gene flow and only weak selection on the loci studied. We derived estimates of migration from linkage disequilibrium in the centre of the zone, and by applying isolation by distance methods to microsatellite data from some of these populations. These give concordant estimates of σ =  0.5–0.8 km generation     . The barrier to gene flow is of the order of 20 km (support limits 14–28), and could be explained by selection of a few per cent at 43–120 underdominant loci that reduces the mean fitness in the central populations to 0.45. Some of the clines appear symmetrical, whereas others are strongly asymmetrical, and two loci appear to have escaped the central barrier to gene flow, reflecting the differential action of selection on different parts of the genome. Asymmetry is always in the direction of more introgression into musculus , indicating either a general progression of domesticus into the musculus territory, possibly mediated by differential behaviour, or past movement of the hybrid zone in the opposite direction, impeded by potential geographical barriers to migration in domesticus territory.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 593–616.     

Inferences from a rapidly moving hybrid zone

Anartia fatima and Anartia amathea (Lepidoptera: Nymphalidae) are sister taxa whose ranges abut in a narrow hybrid zone in eastern Panama. At the center of the zone, hybrids are abundant, although deviations from Hardy-Weinberg and linkage disequilibria are strong, due in part to assortative mating. We measured differences across the zone in four wing color-pattern characters, three allozyme loci, and mitochondrial haplotype. Wing pattern, allozyme, and mitochondrial clines were coincident (i.e., had the same positions) and concordant (i.e.. all markers had similar cline shapes, about 28 km wide). Repeated samples demonstrated that the hybrid zone has been moving eastwards at an average rate of 2.5 km/year over the past 20 years, accompanied by an equivalent movement of the mtDNA cline. No introgression of mtDNA haplotypes were found in the 'wake" of the moving cline, as might be expected for a neutral marker. The concordance of morphological and mtDNA clines between 1994 and 2000, in spite of hybrid zone movement, suggests strong epistasis between the mitochondrial genome and nuclear loci. Cline movement is achieved mainly by pure fatima immigrating into amathea populations; hybrids had little effect, and were presumably outcompeted by fitter pure fatima genotypes. This movement can be explained if random dispersal of 7-19 km x gen(1/2) is coupled with a competitive advantage to A. fatima genomes of 2-5%. Hybrid zone motion is equivalent to Phase III of Wright's shifting balance. Hybrid zone movement has rarely been considered likely in the past, but our results show that it may be more important in biogeography and evolution than generally realized.     

Reproductive biology of the fire-bellied toads, Bombina bombina and B. variegata (Anura: Discoglossidae): egg size, clutch size and larval period length differences

Clutch size, egg size, larval size and duration of larval period were estimated for samples of Bombina bombina and B. variegata from several populations. Bombina variegata lays fewer eggs per batch (mean 17.4) than B. bombina (mean 32.5). The volume of B. variegata eggs is 2.4 times greater and the dry weight 2.3 times greater than that of B. bombina. As a consequence, the volume of B. variegata hatchlings and their dry weight are respectively 1.9 times and 2.5 times greater than that of B. bombina. Analysis of variance showed that variation of egg size in B. variegata is significantly greater than in B. bombina. When reared under laboratory conditions freshly metamorphosed individuals of both species showed no difference in body size or weight. Although growth rates are equal, tadpoles of B. variegata reached metamorphosis 14 days earlier than B. bombina tadpoles because of differences in size of hatchlings. Differences in reproductive parameters found for the two European Bombina species can be interpreted as adaptations to contrasting environments in which they breed.     

Evolution of Bombina bombina and Bombina variegata (Anura: Discoglossidae) in the Carpathian Basin: a history of repeated mt-DNA introgression across species

The structure and geographic location of hybrid zones change through time. Current patterns result from present and historical population-environment interactions that act on each of the hybridizing taxa. This is particularly evident for species involved in complex hybrid zones, such as that formed by the toad species Bombina bombina and Bombina variegata (Anura: Discoglossidae), which interact along extensive areas in Central Europe. We used data on external morphology and partial sequences of the cytochrome oxidase I (cox1) and nicotinamide adenine dinucleotid dehydrogenase subunit 4 (nad4) mitochondrial DNA (mt-DNA) genes to analyze the current patterns of genetic structure shown by both species of Bombina along their contact zone in Hungary. Phylogenetic, phylogeographic, and historical demography analyses were applied to 1.5kb mt-DNA obtained from 119 individuals representing 24 populations from Hungary and additional specimens from Slovakia, Albania, and Bosnia-Herzegovina. We use these data to infer the evolutionary history of the isolated populations of B. variegata in Hungary and to discriminate between competing biogeographic scenarios accounting for the historical interactions between species in this region. Results from the inferred phylogenetic branching pattern and sequence divergence among species and populations support the following: (i) recent population expansion has occurred in Hungarian populations of B. bombina, which are genetically very homogeneous; (ii) the Hungarian populations of B. variegata correspond to two distinct mitochondrial lineages (Carpathian and Alpine, respectively); average maximum-likelihood-corrected sequence divergence between these lineages is 8.96% for cox1 and 10.85% for nad4; (iii) mt-DNA divergence among the three isolated western populations of B. variegata from Transdanubia is low, with four closely related haplotypes, which suggests that the isolation between these populations is the result of a recent process, possibly mediated by the invasion of B. bombina; and (iv) we have detected discordances between morphology and mt-DNA data in the Transdanubia region (Bakony Mountains, Mecsek Mountains, Orség area), suggesting mt-DNA introgression across species in this regions. These results are discussed with reference to previous biogeographic hypotheses.     

Comparative analysis of the genetic structure of hybrid populations of European toads of the genus Bombina (Oken, 1816) from Ukraine and Croatia

The comparative analysis of hybrid populations of toads Bombina bombina and B. variegata from Ukraine and Croatia revealed the essential differences in their genetic structure. The population composition and the direction of genetic introgression are determined initially by the ecological factors (the locality type, the country relief). The clear predominance of the red-bellied toad Bombina bombina was found in two cases: in the lowland populations and in the big samples. The second parental species, B. variegata, predominates in the highlands and in the small samples, respectively.     

Mitochondrial DNA in the hybrid zone between Mus musculus musculus and Mus musculus domesticus: a comparison of two transects

We studied mtDNA introgression across the contact zone between Mus musculus musculus and M. m. domesticus in two independent transects in the Czech Republic and Bavaria, Germany. A total of 1270 mice from 98 localities in the Czech transect and 456 mice from 41 localities in the Bavarian transect were examined for presence or absence of a Bam HI restriction site in the mt-Nd1 gene. Using this simple mtDNA marker, variants that belonged to the M. m. domesticus lineage (presence of restriction site) could be unequivocally distinguished from those belonging to the M. m. musculus lineage (absence of restriction site). The extent of introgression of mtDNA, three autosomal allozymes and the X chromosome was compared. The introgression of X markers was more limited than was that of the allozymes and mtDNA. In the Czech transect, the centre for the mtDNA cline was shifted about 3.6 km to the west relative to the X chromosome cline, with asymmetric introgression from M. m. musculus to M. m. domesticus . Interestingly, in the Bavarian transect, the centre of the mtDNA cline was shifted about 10.9 km to the east relative to the X chromosome cline, with asymmetric introgression from M. m. domesticus to M. m. musculus, opposite in direction to that observed in the Czech transect.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 363–378.     

Quantitative immunological studies of the albumins of several species of fire bellied toads, genus Bombina.

1. Rabbit antisera against purified serum albumin of Bombina bombina were used to study relationships between B. bombina, B. variegata, and B. orientalis. 2. Quantitative micro-complement fixation tests indicated the albumins of B. bombina from central Poland and Bulgaria were indistinguishable. The albumins from several populations of B. variegata differed very slightly from that of B. bombina. The albumin of B. orientalis was quite distinct from that of B. bombina. 3. Using albumin as a molecular clock, we estimated B. bombina and B. variegata diverged within the last million years, whereas the B. orientalis lineage diverged roughly 10-12 mil yr ago.     

A CONTACT ZONE WITH NONCOINCIDENT CLINES FOR SEX-SPECIFIC MARKERS IN THE FIELD VOLE (MICROTUS AGRESTIS)

A field vole (Microtus agrestis) population characterized by unique mitochondrial DNA (mtDNA) and Y chromosome markers occurs in southwest Sweden. A contact zone between this Lund (Lu) population and a standard (St) population was examined with two sex-specific genetic markers. The field vole mtDNA and Y chromosome clines display a remarkable lack of coincidence, rarely observed in contact zones. The cline width for both markers is about 50 km, but the two clines are displaced from each other: the mtDNA cline is found in the central part of the study area, whereas the cline for the Y chromosome is located in the eastern part. Thus, the absolute width of the Lu-St zone spans about 95–110 km. As a result of the cline shift, all male hybrids carry the Lu-Y chromosome and St-type mtDNA. The other possible male hybrid class is lacking. The distinct noncoincidence of the mtDNA and Y chromosome clines is most likely explained by selection against male hybrids with the St-Y chromosome and Lu-mtDNA. It is possible that incompatibilities between the maternal genome of Lu-type animals and the paternal genome of St-type individuals exist. However, alternative explanations based on neutrality or selective advantage cannot be totally dismissed.     

Morphological,genetic and behavioural analyses of a hybrid zone between the ground beetles Carabus lewisianus and C. albrechti (Coleoptera,Carabidae): asymmetrical introgression caused by movement of the zone?

We analysed the putative hybrid zone between the ground beetles Carabus lewisianus Breuning and C. albrechti Morawitz from the Kanto and Tanzawa Mountains in central Honshu, Japan, using morphological, genetic, and behavioural data. Canonical discriminant analysis of three external and five genital morphological characters revealed an apparent morphological gap, suggesting restricted sympatry in the contact zone. RFLP and molecular phylogenetic analyses of the mitochondrial gene ND5 revealed that haplotypes originating from C. lewisianus introgressed extensively into C. albrechti . Fitting tanh curves to morphological (external and genital) and genetic characters showed differential cline widths and centres, suggesting stronger selection for genital morphology with narrower cline widths, and extensive introgression of the nuclear genes responsible for external characters relative to those controlling genital morphology into the range of C. albrechti . Linkage disequilibria between the morphological characters did not differ from zero, suggesting sufficient recombination occurred in the zone but the effect of small sample sizes was not negligible. Mate choice trials demonstrated that males showed a weak but not significant preference for conspecific mates. We interpret these patterns as evidence for (1) partial reproductive isolation presumably by postmating processes; (2) selection against intermediate genital morphologies, and (3) independent responses of the characters to different evolutionary forces that may cause movement of genital clines and introgression of the mitochondrial gene and external morphology. The movement of the boundaries may have resulted in the particular species distribution patterns in the study area.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 86 , 79–94.     

Steep,coincident, and concordant clines in mitochondrial and nuclear‐encoded genes in a hybrid zone between subspecies of Atlantic killifish,Fundulus heteroclitus

Steep genetic clines resulting from recent secondary contact between previously isolated taxa can either gradually erode over time or be stabilized by factors such as ecological selection or selection against hybrids. We used patterns of variation in 30 nuclear and two mitochondrial SNPs to examine the factors that could be involved in stabilizing clines across a hybrid zone between two subspecies of the Atlantic killifish, Fundulus heteroclitus. Increased heterozygote deficit and cytonuclear disequilibrium in populations near the center of the mtDNA cline suggest that some form of reproductive isolation such as assortative mating or selection against hybrids may be acting in this hybrid zone. However, only a small number of loci exhibited these signatures, suggesting locus‐specific, rather than genomewide, factors. Fourteen of the 32 loci surveyed had cline widths inconsistent with neutral expectations, with two SNPs in the mitochondrial genome exhibiting the steepest clines. Seven of the 12 putatively non‐neutral nuclear clines were for SNPs in genes related to oxidative metabolism. Among these putatively non‐neutral nuclear clines, SNPs in two nuclear‐encoded mitochondrial genes (SLC25A3 and HDDC2), as well as SNPs in the myoglobin, 40S ribosomal protein S17, and actin‐binding LIM protein genes, had clines that were coincident and concordant with the mitochondrial clines. When hybrid index was calculated using this subset of loci, the frequency distribution of hybrid indices for a population located at the mtDNA cline center was non‐unimodal, suggesting selection against advanced‐generation hybrids, possibly due to effects on processes involved in oxidative metabolism.     

Phylogeography and subspecies taxonomy of dunlins (Calidris alpina) in western Palearctic analysed by DNA microsatellites and amplified fragment length polymorphism markers

In a circumpolar wader, the dunlin ( Calidris alpina ), there are 11 named subspecies, but only five mitochondrial DNA (mtDNA) lineages have been found. In the present study, we investigated the genetic structure of dunlins in western Palearctic (from East Greenland to Taimyr peninsula) using DNA microsatellites and amplified fragment length polymorphism (AFLP) markers that may detect more recent differentiation than mtDNA. In this region, we consider four described subspecies; alpina , schinzii , arctica and centralis , together comprising two mtDNA lineages. We analyse seven polymorphic microsatellite loci and 91 AFLP markers in 287 and 152 unrelated individuals, respectively, originating from 17 populations. Neither microsatellites nor AFLPs reveal distinct groups that correspond to currently recognized subspecies. There is a clear pattern of isolation by distance in microsatellites. Our results do not contradict the former mtDNA results that there are two phylogenetic lineages (approximately corresponding to schinzii and centralis ) that have met and formed a cline ( alpina ). We find no difference between schinzii and arctica (East Greenland). We conclude that, given the lack of distinct groups and the gradual changes in microsatellite allele frequencies, these markers provide little genetic support for the dunlin subspecies taxonomy in the western Palearctic.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 92 , 713–726.     

Falling apart and merging: diversification of slender salamanders (Plethodontidae: Batrachoseps) in the American West

The plethodontid genus Batrachoseps , the slender salamanders, is the most diverse clade of salamanders in western North America, but it has posed taxonomic difficulties because it contains many morphologically cryptic species. A segment of the mitochondrial DNA gene cytochrome b was studied for 278 individuals densely sampled from throughout the range of all 18 described species and several undescribed species. Phylogenetic analyses of the mtDNA data identify six major clades, one corresponding to the subgenus Plethopsis and five within a monophyletic subgenus Batrachoseps. All major clades and most species within these clades display strong phylogeographic structuring. Comparisons of mtDNA and allozyme data show that several allozymically cohesive groups are not monophyletic with respect to mtDNA. We suggest that this phenomenon results from fragmentation of populations, divergence in allopatry, and then recontact and gradual merging of units caused predominantly by male-mediated gene flow. The mtDNA offers evidence that populations were once more isolated than they are now, while the patterns of allozyme variation reflect recent and current interactions among populations. The complex patterns of morphological, allozymic and mtDNA variation associated with the constantly changing geological landscape give insight into the nature of processes responsible for species formation in Batrachoseps .  © 2002 The Linnean Society of London. Biological Journal of the Linnean Society , 2002, 76 , 361–391.     

Phylogenetic analyses of mtDNA sequences reveal three cryptic lineages in the widespread neotropical frog Leptodactylus fuscus (Schneider, 1799) (Anura, Leptodactylidae)

Leptodactylus fuscus is a neotropical frog ranging from Panamá to Argentina, to the east of the Andes mountains, and also inhabiting Margarita, Trinidad, and the Tobago islands. We performed phylogenetic analyses of 12S rRNA, 16S rRNA, tRNA-Leu, and ND1 mitochondrial (mt) DNA sequences from specimens collected across the geographic distribution of L. fuscus to examine two alternative hypotheses: (i) L. fuscus is a single, widely distributed species, or (ii) L. fuscus is a species complex. We tested statistically for geographic association and partitioning of genetic variation among mtDNA clades. The mtDNA data supported the hypothesis of several cryptic species within L. fuscus. Unlinked mtDNA and nuclear markers supported independently the distinctness of a 'northern' phylogenetic unit. In addition, the mtDNA data divided the southern populations into two clades that showed no sister relationship to each other, consistent with high differentiation and lack of gene flow among southern populations as suggested by allozyme data. Concordance between mtDNA and allozyme patterns suggests that cryptic speciation has occurred in L. fuscus without morphological or call differentiation. This study illustrates a case in which lineage splitting during the speciation process took place without divergence in reproductive isolation mechanisms (e.g. advertisement call in frogs), contrary to expectations predicted using a biological species framework.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society, 2006, 87 , 325–341. No claim to original US government works.     

Markedly discordant mitochondrial DNA and allozyme phylogenies of tube-nosed fruit bats, Nyctimene, at the Australian–Oriental biogeographical interface

Australo-Papuan tube-nosed fruit bats of the genus Nyctimene reach their most westerly distribution on the island archipelagos of eastern Indonesia. A recent morphological examination indicates three species occur on Moluccan islands. Both allozyme electrophoresis and mitochondrial DNA (mtDNA) sequencing reveal there are only two species present, Nyctimene cephalotes and Nyctimene albiventer , but there is considerable disparity between the results obtained from the two genetic data sets. Allozyme data indicate N. albiventer occurs on Wokam, which sits on the Australian continental shelf and was joined to Australo-Papua during the last glacial maximum, and N. cephalotes on the other Moluccan islands, all of which are off the continental shelf. Divergence of these two species is dated at approximately 1.5 Mya. By contrast, the mtDNA gene tree shows two deep clades, one containing all specimens from Wokam and Yamdena, and the second all the specimens from the other islands. This especially marked incongruence between the two genetic data sets is ascribed to either a single introgression event of N. albiventer mtDNA into N. cephalotes on Yamdena at around 100 000 ya or lineage sorting of very old mtDNA lineages that coalesce a rather long time before the speciation event. These results highlight that caution should be exercised when relying on mtDNA as barcodes in species taxonomy. The continental-associated N. albiventer on Wokam has higher levels of allozyme heterozygosity and mtDNA nucleotide diversity than the N. cephalotes populations occurring on the more remote islands, indicating that this colonizing species has experienced bottlenecks and/or a low effective population size since speciation.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 93 , 589–602.     

Evolution of mosaic hybrid zone between invasive and endemic species of Cyprinidae through space and time

Hybridization phenomena can be considered in different ways. Hybridization studies can be used to elucidate some aspects of speciation and adaptation. However, hybridization may be a threat to endemic species. Here, we studied bidirectional introgression between two sympatric species of Cyprinidae: Chondrostoma toxostoma toxostoma and C. nasus nasus . Analysis of morphology, allozymes and mtDNA sequences revealed that the hybridization between these two cyprinids takes different forms in the same river. The mosaic hybrid zone (so called because of the absence of a simple cline) appeared to be partitioned due to the proportions of the two species along a spatial scale. The proportion of each hybrid group in this zone was unstable over a 5-year period. We propose various protection/conservation states for the C. t. toxostoma populations based on the level of introgressive hybridization and the genetic structure of both species in 'pure' populations.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 85 , 135–155.