HYBRIDIZATION OF BOMBINA BOMBINA AND B. VARIEGATA (ANURA, DISCOGLOSSIDAE) AT A SHARP ECOTONE IN WESTERN UKRAINE: COMPARISONS ACROSS TRANSECTS AND OVER TIME

Population genetic structure and hybrid zone analyses for species delimitation in the Japanese toad (Bufo japonicus)

Hybridization following secondary contact may produce different outcomes depending on the extent to which genetic diversity and reproductive barriers have accumulated during isolation. The Japanese toad, Bufo japonicus, is distributed on the main islands of Japan. In the present study, we applied multiplexed inter-simple sequence repeat genotyping by sequencing to achieve the fine-scale resolution of the genetic cluster in B. j. japonicus and B. j. formosus. We also elucidated hybridization patterns and gene flow degrees across contact zones between the clusters identified. Using SNP data, we found four genetic clusters in B. j. japonicus and B. j. formosus and three contact zones of the cluster pairs among these four clusters. The two oldest diverged lineages, B. j. japonicus and B. j. formosus, formed a narrow contact zone consistent with species distinctiveness. Therefore, we recommend that these two subspecies be elevated to the species level. In contrast, the less diverged pairs of two clusters in B. j. japonicus and B. j. formosus, respectively, admixed over a hundred kilometers, suggesting that they have not yet developed strong reproductive isolation and need to be treated as conspecifics. These results will contribute to resolving taxonomic confusion in Japanese toads.

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.

Evolutionary Relationships and Reproductive Isolating Mechanisms in the Rice Frog (Fejervarya limnocharis) Species Complex from Sri Lanka, Thailand, Taiwan and Japan, Inferred from mtDNA Gene Sequences, Allozymes, and Crossing Experiments

The rice frog (Fejervarya limnocharis) species complex is widely distributed, from India to Japan, and most prevalently in Southeast Asia. Conspicuous morphological variation has been reported for this species complex throughout its distribution range. In the present study, we used mtDNA gene sequence and allozyme analyses to infer evolutionary affinities within this species complex using eight populations (Sri Lanka; Bangkok and Ranong in Thailand; Taiwan; and Hiroshima, Okinawa, Ishigaki and Iriomote in Japan). We also conducted crossing experiments among four populations from Japan, Thailand, and Sri Lanka in order to find out more about the reproductive isolating mechanisms that might exist among the East, Southeast, and South Asian populations of this species complex. The crossing experiments revealed that the Sri Lanka population is reproductively isolated from the Hiroshima, Bangkok, and Ranong populations by complete hybrid inviability, and that the Bangkok population may be reproductively isolated from the Hiroshima population by partial hybrid inviability. Thus, it is not unreasonable to regard the Sri Lanka population as a species separated from F. limnocharis. The mtDNA and allozyme data showed that the Ranong population is most closely related to the Bangkok population in nuclear genome, but more similar to the Okinawa and Taiwan populations in mtDNA genome. The present, preliminary survey may raise questions about the species status of these particular populations and also about the nature of the biological species concept.