Highly Complex Mitochondrial Dna Genealogy in an Endemic Japanese Subterranean Breeding Brown FrogRana Tagoi(Amphibia, Anura, Ranidae)

Multiple Transitions between Y Chromosome and Autosome in Tago's Brown Frog Species Complex

Sex chromosome turnover is the transition between sex chromosomes and autosomes. Although many cases have been reported in poikilothermic vertebrates, their evolutionary causes and genetic mechanisms remain unclear. In this study, we report multiple transitions between the Y chromosome and autosome in the Japanese Tago's brown frog complex. Using chromosome banding and molecular analyses (sex-linked and autosomal single nucleotide polymorphisms, SNPs, from the nuclear genome), we investigated the frogs of geographic populations ranging from northern to southern Japan of two species, Rana tagoi and Rana sakuraii (2n = 26). Particularly, the Chiba populations of East Japan and Akita populations of North Japan in R. tagoi have been, for the first time, investigated here. As a result, we identified three different sex chromosomes, namely chromosomes 3, 7, and 13, in the populations of the two species. Furthermore, we found that the transition between the Y chromosome (chromosome 7) and autosome was repeated through hybridization between two or three different populations belonging to the two species, followed by restricted chromosome introgression. These dynamic sex chromosome turnovers represent the first such findings in vertebrates and imply that speciation associated with inter- or intraspecific hybridization plays an important role in sex chromosome turnover in frogs.

Nuclear phylogenies and genomics of a contact zone establish the species rank of Podarcis lusitanicus (Squamata, Lacertidae)

Unravelling when divergent lineages constitute distinct species can be challenging, particularly in complex scenarios combining cryptic diversity and phylogenetic discordances between different types of molecular markers. Combining a phylogenetic approach with the study of contact zones can help to overcome such difficulties. The Podarcis hispanicus species complex has proven to be prosperous in independent evolutionary units, sometimes associated with cryptic diversity. Previous studies have revealed that one of the species of this complex, P. guadarramae, comprises two deeply divergent yet morphologically indistinguishable evolutionary units, currently regarded as subspecies (P. g. guadarramae and P. g. lusitanicus). In this study we used molecular data to address the systematics of the two lineages of Podarcis guadarramae and the closely related P. bocagei. Firstly, we reconstructed the species tree of these three and two additional taxa based on 30 nuclear loci using the multispecies coalescent with and without gene flow. Secondly, we used SNPs obtained from RADseq data to analyze the population structure across the distribution limits P. g. lusitanicus and P. g. guadarramae, and for comparison, a contact zone between P. bocagei and P. g. lusitanicus. Nuclear phylogenetic relationships between these three taxa are clearly difficult to determine due to the influence of gene flow, but our results give little support to the monophyly of P. guadarramae, potentially due to a nearly simultaneous divergence between them. Genetic structure and geographic cline analysis revealed that the two lineages of P. guadarramae replace each other abruptly across the sampled region and that gene flow is geographically restricted, implying the existence of strong reproductive isolation. Podarcis bocagei and P. g. lusitanicus show a similar degree of genetic differentiation and reproductive isolation, with very low levels of admixture in syntopy. These results support that all three forms are equally differentiated and reproductively isolated. In consequence, we conclude that the two former subspecies of Podarcis guadarramae constitute valid, yet cryptic species, that should be referred to as P. lusitanicus and P. guadarramae.

Evolution of Sex Chromosome Heteromorphy in Geographic Populations of the Japanese Tago's Brown Frog Complex

The sex chromosomes of most anuran amphibians are characterized by homomorphy in both sexes, and evolution to heteromorphy rarely occurs at the species or geographic population level. Here, we report sex chromosome heteromorphy in geographic populations of the Japanese Tago's brown frog complex (2n = 26), comprising Rana sakuraii and R. tagoi. The sex chromosomes of R. sakuraii from the populations in western Japan were homomorphic in both sexes, whereas chromosome 7 from the populations in eastern Japan were heteromorphic in males. Chromosome 7 of R. tagoi, which is distributed close to R. sakuraii in eastern Japan, was highly similar in morphology to the Y chromosome of R. sakuraii. Based on this and on mitochondrial gene sequence analysis, we hypothesize that in the R. sakuraii populations from eastern Japan the XY heteromorphic sex chromosome system was established by the introduction of chromosome 7 from R. tagoi via interspecies hybridization. In contrast, chromosome 13 of R. tagoi from the 2 large islands in western Japan, Shikoku and Kyushu, showed a heteromorphic pattern of constitutive heterochromatin distribution in males, while this pattern was homomorphic in females. Our study reveals that sex chromosome heteromorphy evolved independently at the geographic lineage level in this species complex.