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.

Parallel Evolution of Sex-Linked Genes across XX/XY and ZZ/ZW Sex Chromosome Systems in the Frog Glandirana rugosa.

Genetic sex-determination features male (XX/XY) or female heterogamety (ZZ/ZW). To identify similarities and differences in the molecular evolution of sex-linked genes between these systems, we directly compared the sex chromosome systems existing in the frog Glandirana rugosa. The heteromorphic X/Y and Z/W sex chromosomes were derived from chromosomes 7 (2n = 26). RNA-Seq, de novo assembly, and BLASTP analyses identified 766 sex-linked genes. These genes were classified into three different clusters (XW/YZ, XY/ZW, and XZ/YW) based on sequence identities between the chromosomes, probably reflecting each step of the sex chromosome evolutionary history. The nucleotide substitution per site was significantly higher in the Y- and Z-genes than in the X- and W- genes, indicating male-driven mutation. The ratio of nonsynonymous to synonymous nucleotide substitution rates was higher in the X- and W-genes than in the Y- and Z-genes, with a female bias. Allelic expression in gonad, brain, and muscle was significantly higher in the Y- and W-genes than in the X- and Z-genes, favoring heterogametic sex. The same set of sex-linked genes showed parallel evolution across the two distinct systems. In contrast, the unique genomic region of the sex chromosomes demonstrated a difference between the two systems, with even and extremely high expression ratios of W/Z and Y/X, respectively.

Genetic and morphological variation analyses of Glandirana rugosa with description of a new species (Anura, Ranidae).

Glandirana rugosa is known to include several geographic groups differing in sex chromosomes, and has been proven to be paraphyletic in mitochondrial phylogeny with respect to G. susurra. By analyzing genetic and morphological variation in a large number of individuals of Glandirana, we studied their taxonomic relationships. A mitochondrial DNA phylogeny, with the G. tientaiensis as outgroup, revealed two major lineages containing respectively (1) the East group of G. rugosa, G. susurra, and the Central and Southeast-Kyushu groups of G. rugosa; and (2) G. emeljanovi, and the North and West groups of G. rugosa. In contrast, in a nuclear DNA phylogeny based on SNP data, lineages of (1) G. susurra and East group, and (2) the remaining groups of G. rugosa and G. emeljanovi, were split, indicating a distinct status of the East group among G. rugosa. In adult morphology, there were only minor differences between the East group and the remaining groups of G. rugosa, but in larvae, the East group had significantly more sparse skin glands than the others. The exact type locality of G. rugosa is most probably in western Japan, not including the range of the East group. From these results, we describe the East group as a new species, G. reliquia, distinct from the remaining groups of G. rugosa. The new species with sexually homomorphic chromosomes is thought to represent a basic stock of Japanese Glandirana, which existed far before G. rugosa originated.

Identification of ancestral sex chromosomes in the frog Glandirana rugosa bearing XX-XY and ZZ-ZW sex-determining systems.

Sex chromosomes constantly exist in a dynamic state of evolution: rapid turnover and change of heterogametic sex during homomorphic state, and often stepping out to a heteromorphic state followed by chromosomal decaying. However, the forces driving these different trajectories of sex chromosome evolution are still unclear. The Japanese frog Glandirana rugosa is one taxon well suited to the study on these driving forces. The species has two different heteromorphic sex chromosome systems, XX-XY and ZZ-ZW, which are separated in different geographic populations. Both XX-XY and ZZ-ZW sex chromosomes are represented by chromosome 7 (2n = 26). Phylogenetically, these two systems arose via hybridization between two ancestral lineages of West Japan and East Japan populations, of which sex chromosomes are homomorphic in both sexes and to date have not yet been identified. Identification of the sex chromosomes will give us important insight into the mechanisms of sex chromosome evolution in this species. Here, we used a high-throughput genomic approach to identify the homomorphic XX-XY sex chromosomes in both ancestral populations. Sex-linked DNA markers of West Japan were aligned to chromosome 1, whereas those of East Japan were aligned to chromosome 3. These results reveal that at least two turnovers across three different sex chromosomes 1, 3 and 7 occurred during evolution of this species. This finding raises the possibility that cohabitation of the two different sex chromosomes from ancestral lineages induced turnover to another new one in their hybrids, involving transition of heterogametic sex and evolution from homomorphy to heteromorphy.

Draft genome sequence data of Indian rhinoceros, Rhinoceros unicornis.

The Indian rhinoceros (Rhinoceros unicornis) is a large herbivore found in northern India and southern Nepal. It is a critically endangered species, with an estimated population of approximately 3,600 in the wild. Genetic factors, such as the loss of genetic diversity and the accumulation of deleterious variations, are critical risk factors for the extinction of endangered species, such as the Indian rhinoceros. To support the conservation efforts of the Indian rhinoceros, we assembled its draft genome. The new genomic data will enable the study of functional genes associated with the ecological and physiological characteristics of Indian rhinoceros and help us establish more effective conservation measures. The muscles of an Indian rhinoceros that died from prostration at a zoo were collected, and the samples were stored at the National Institute for Environmental Studies (Tsukuba, Japan). Sequence data were obtained using an Illumina NovaSeq 6000 platform for short reads and an Oxford Nanopore Technologies PromethION for long reads. We generated approximately 235.2 Gbp of data. From these sequences, we assembled a 2,375,051,758 bp genome consisting of 7,615 contigs. The genome data are available from the National Center Biotechnology Information BioProject database under accession number BOSQ00000000.

Sex chromosome evolution from a heteromorphic to a homomorphic system by inter-population hybridization in a frog.

Sex chromosomes generally evolve from a homomorphic to heteromorphic state. Once a heteromorphic system is established, the sex chromosome system may remain stable for an extended period. Here, we show the opposite case of sex chromosome evolution from a heteromorphic to a homomorphic system in the Japanese frog Glandirana rugosa. One geographic group, Neo-ZW, has ZZ-ZW type heteromorphic sex chromosomes. We found that its western edge populations, which are geographically close to another West-Japan group with homomorphic sex chromosomes of XX-XY type, showed homozygous genotypes of sex-linked genes in both sexes. Karyologically, no heteromorphic sex chromosomes were identified. Sex-reversal experiments revealed that the males were heterogametic in sex determination. In addition, we identified another similar population around at the southwestern edge of the Neo-ZW group in the Kii Peninsula: the frogs had homomorphic sex chromosomes under male heterogamety, while shared mitochondrial haplotypes with the XY group, which is located in the east and bears heteromorphic sex chromosomes. In conclusion, our study revealed that the heteromorphic sex chromosome systems independently reversed back to or turned over to a homomorphic system around each of the western and southwestern edges of the Neo-ZW group through hybridization with the West-Japan group bearing homomorphic sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Comparative genomics of Glandirana rugosa using unsupervised AI reveals a high CG frequency.

The Japanese wrinkled frog (Glandirana rugosa) is unique in having both XX-XY and ZZ-ZW types of sex chromosomes within the species. The genome sequencing and comparative genomics with other frogs should be important to understand mechanisms of turnover of sex chromosomes within one species or during a short period. In this study, we analyzed the newly sequenced genome of G. rugosa using a batch-learning self-organizing map which is unsupervised artificial intelligence for oligonucleotide compositions. To clarify genome characteristics of G. rugosa, we compared its short oligonucleotide compositions in all 1-Mb genomic fragments with those of other six frog species (Pyxicephalus adspersus, Rhinella marina, Spea multiplicata, Leptobrachium leishanense, Xenopus laevis, and Xenopus tropicalis). In G. rugosa, we found an Mb-level large size of repeat sequences having a high identity with the W chromosome of the African bullfrog (P. adspersus). Our study concluded that G. rugosa has unique genome characteristics with a high CG frequency, and its genome is assumed to heterochromatinize a large size of genome via methylataion of CG.

Comparative analysis of sperm motility in liquid and seminal coagulum portions between Bornean orangutan (Pongo pygmaeus) and chimpanzee (Pan troglodytes).

Coagulum in the semen of some primates plays different roles, depending on the species. In the present study, we examined sperm motility in the coagulum and liquid portions of semen collected from captive individuals from two great ape species: two adult Bornean orangutans (Pongo pygmaeus) (n = 27) and three adult chimpanzees (Pan troglodytes) (n = 14). The results revealed that orangutan sperm remained motile for significantly longer in the coagulum than in the liquid portion (> 18 h). By contrast, chimpanzee sperm motility did not differ significantly over time between the two portions of the semen, although motility was slightly higher in the liquid portion than in the coagulum. The evolution of the seminal coagulum is thought to be related to postcopulatory sperm competition; however, functions of seminal coagulum have not been completely elucidated. Our data from the orangutan semen suggest that in this species, seminal coagulum may strengthen own-sperm survival. This report is the first to provide evidence for this distinctive function of the seminal coagulum. This unique property of orangutan seminal coagulum might be attributable to their reproductive traits, e.g., difficulty in predicting ovulation due to a lack of genital swelling during estrus. The orangutan is a Critically Endangered species, and captive breeding, including artificial insemination (AI), is expected. However, worldwide, only one case of orangutan AI has been successful. Our findings may contribute to an understanding of their basic semen characteristics and help improve the AI method.

Reconstruction of female heterogamety from admixture of XX-XY and ZZ-ZW sex-chromosome systems within a frog species.

Sex-determining mechanisms change repeatedly throughout evolution, and it is difficult to track this continual process. The Japanese soil-frog Glandirana rugosa is a remarkable evolutionary witness to the ongoing process of the evolution of sex-determining modes. The two geographic groups, designated XY and Neo-ZW, have homologous sex chromosomes, yet display opposite types of sex chromosomes, XX-XY and ZZ-ZW, respectively. These two groups are sympatric at the edges of their respective ranges in Central Japan. In this study, we discovered molecular evidence that the eastern part of the Neo-ZW group (Neo-ZW2 subgroup), which is found near the sympatric area, shares mitochondrial haplotypes with the XY group. By analysing single nucleotide polymorphism (SNP) loci, we have also discovered that the representative nuclear genome of the Neo-ZW2 subgroup shares allele clusters with both the XY group and another part of the Neo-ZW group (Neo-ZW1 subgroup), indicating a hybrid origin of the Neo-ZW2. Further analysis of sex-linked SNP loci revealed that the alleles on the W chromosomes of the Neo-ZW2 were derived mostly from X chromosomes, while alleles on the Z chromosomes originated from the Z chromosomes of the Neo-ZW1 subgroup and partly from the Y chromosomes of the XY group. Our study revealed that admixture of the two opposite sex-chromosome systems reconstructed a female heterogametic system by recycling the X chromosomes into new W chromosomes. This work offers an illustrative example of how de novo sex-chromosome systems can arise by recycling material from ancestral sex chromosomes.

The Distributions and Boundary of Two Distinct, Local Forms of Japanese Pond Frog, Pelophylax porosus brevipodus, Inferred From Sequences of Mitochondrial DNA.

The Nagoya Daruma pond frog Pelophylax porosus brevipodus is distributed in western Japan and is traditionally divided into two local forms: the Okayama form in the west and the Nagoya form in the east. These two forms are genetically differentiated, but have never been defined taxonomically because their distributions are unclear to date. To complete the distributions and identify the boundary of the two forms, we genetically investigated 16 populations including eight populations located within the unexamined area. We found that the distributional boundary is located within a small area of Hyogo Prefecture where haplotypes of mitochondrial cytochrome b (cytb) and D-loop region corresponding to the two forms co-existed. On the other hand, the polymorphic site of the nuclear gene SOX3 revealed introgression over the boundary into Okayama cytb clade. These results suggest that the two forms were geographically isolated from each other in the past, and secondarily contacted and then accepted one-way introgression. As a next step of the research, taxonomic approach is expected to define the two forms.

Spontaneous tyrosinase mutations identified in albinos of three wild frog species.

The present study reports spontaneous tyrosinase gene mutations identified in oculocutaneous albinos of three Japanese wild frog species, Pelophylax nigromaculatus, Glandirana rugosa and Fejervarya kawamurai. This represents the first molecular analyses of albinic phenotypes in frogs. Albinos of P. nigromaculatus collected from two different populations were found to suffer from frameshift mutations. These mutations were caused by the insertion of a thymine residue within each of exons 1 and 4, while albinos in a third population lacked three nucleotides encoding lysine in exon 1. Albinos from the former two P. nigromaculatus populations were also associated with splicing variants of mRNA that lacked either exons 2-4 or exon 4. In the other two frog species examined, missense mutations that resulted in amino acid substitutions from glycine to arginine and glycine to aspartic acid were identified in exons 1 and 3, respectively. The two glycines in F. kawamurai and G. rugosa, and the lysine deleted in one P. nigromaculatus albino, were highly conserved in vertebrates, which suggested that they were situated in regions of critical importance to tyrosinase function. In fact, the glycine of G. rugosa is located within a predicted copper-binding domain. The five mutations identified in the present study are candidates for causing the albinic phenotypes, and, if directly confirmed, they are all unique among vertebrates, which suggests that molecular analysis of albino frogs could contribute to research on albinos in humans and vertebrates by providing new information about tyrosinase structure and transcript processing.

Evolutionary Changes in Sensitivity to Hormonally Induced Gonadal Sex Reversal in a Frog Species.

The Japanese frog Glandirana rugosa is unique in that it shows geographic variation in sex chromosome differentiation and heterogametic sex determination. To elucidate the cause of interpopulation differences in gonadal sex differentiation, we investigated hormonally induced sex reversal and the expression patterns of genes associated with sex determination during early tadpole development. We found that sex reversal was easily induced in XX females and XY males of 2 forms (West-Japan and East-Japan) of G. rugosa with the ancestral homomorphic sex chromosomes under male heterogametic sex determination. During sex reversal, expression of CYP19 and/or FOXL2 was dependent on the phenotypic sex of the gonad. In contrast, sex reversal was not induced in ZW females of a population with a heteromorphic ZW sex chromosome system or in XX females or XY males in a population with a heteromorphic XY sex chromosome system. The latter 2 populations are evolutionarily derived forms. These results indicate an evolutionary direction for the gonadal sex differentiation mechanism. The original system was highly sensitive to sex hormones and allowed almost complete sex reversal. From this ancestral form, a new system evolved that was resistant to hormones and showed a change in the heterogametic sex and the sex chromosome differentiation mechanism.

Meiotic recombination counteracts male-biased mutation (male-driven evolution).

Meiotic recombination is believed to produce greater genetic variation despite the fact that deoxyribonucleic acid (DNA)-replication errors are a major source of mutations. In some vertebrates, mutation rates are higher in males than in females, which developed the theory of male-driven evolution (male-biased mutation). However, there is little molecular evidence regarding the relationships between meiotic recombination and male-biased mutation. Here we tested the theory using the frog Rana rugosa, which has both XX/XY- and ZZ/ZW-type sex-determining systems within the species. The male-to-female mutation-rate ratio (α) was calculated from homologous sequences on the X/Y or Z/W sex chromosomes, which supported male-driven evolution. Surprisingly, each α value was notably higher in the XX/XY-type group than in the ZZ/ZW-type group, although α should have similar values within a species. Interestingly, meiotic recombination between homologous chromosomes did not occur except at terminal regions in males of this species. Then, by subdividing α into two new factors, a replication-based male-to-female mutation-rate ratio (ß) and a meiotic recombination-based XX-to-XY/ZZ-to-ZW mutation-rate ratio (γ), we constructed a formula describing the relationship among a nucleotide-substitution rate and the two factors, ß and γ. Intriguingly, the ß- and γ-values were larger and smaller than 1, respectively, indicating that meiotic recombination might reduce male-biased mutations.

Genetic analysis of captive proboscis monkeys.

Information on the genetic relationships of captive founders is important for captive population management. In this study, we investigated DNA polymorphisms of four microsatellite loci and the mitochondrial control region sequence of five proboscis monkeys residing in a Japanese zoo as captive founders, to clarify their genetic relationship. We found that two of the five monkeys appeared to be genetically related. Furthermore, the haplotypes of the mitochondrial control region of the five monkeys were well differentiated from the haplotypes previously reported from wild populations from the northern area of Borneo, indicating a greater amount of genetic diversity in proboscis monkeys than previously reported.

Independent degeneration of W and Y sex chromosomes in frog Rana rugosa.

The frog Rana rugosa uniquely possesses two different sex-determining systems of XX/XY and ZZ/ZW, separately in the geographic populations. The sex chromosomes of both types share the same origin at chromosome 7, and the structural differences between X and Y or Z and W were evolved through two inversions. In order to ascertain the mechanisms of degeneration of W and Y chromosomes, we gynogenetically produced homozygous diploids WW and YY and examined their viability. Tadpoles from geographic group N (W(N)W(N)) containing three populations died of edema at an early developmental stage within 10 days after hatching, while tadpoles from the geographic group K (W(K)W(K)) that contained two populations died of underdeveloped growth at a much later stage, 40-50 days after fertilization. On the contrary, W(N)W(K) and W(K)W(N) hybrid embryos were viable, successfully passed the two lethal stages, and survived till the attainment of adulthood. The observed survival implies that the lethal genes of the W chromosomes are not shared by the two groups and thus demonstrates their independent degeneration histories between the local groups. In sharp contrast, a sex-linked gene of androgen receptor gene (AR) from the W chromosome was down-regulated in expression in both the groups, suggesting that inactivation of the W-AR allele preceded divergence of the two groups and appearance of the lethal genes. Besides, the YY embryos died of cardiac edema immediately after hatching. The symptom of lethality and the stage of developmental arrest differed from those for either of WW lethal embryos. We therefore conclude that the W and Y chromosomes involve no evolutionary common scenario for degeneration.

Noninvasive monitoring of reproductive activity based on fecal progestagen profiles and sexual behavior in Koalas, Phascolarctos cinereus.

Studies on the reproductive endocrinology of koalas have been performed mainly by using blood samples; however, in practice it is difficult to collect blood periodically because koalas are easily stressed. The purposes of the present study were to establish a noninvasive endocrine monitoring technique and to investigate the reproductive physiology of female koalas. Feces were collected from female northern and southern koalas, and progestagen was extracted from lyophilized fecal samples and determined by enzyme immunoassay. In nonpregnant northern and southern koalas, fecal progestagen markedly increased after copulation and remained high for 36.3 +/- 2.5 days and 38.9 +/- 1.4 days (luteal phase, mean +/- SEM), respectively. Mean (+/-SEM) progestagen levels (6.34 +/- 0.49 microg/g) during the luteal phase in northern koalas were significantly higher than in southern koalas (4.19 +/- 0.24 microg/g). Fecal progestagen in parturient northern koalas remained high for 36.2 +/- 1.9 days (gestation period, 34.1 +/- 0.3 days). In northern koalas, the mean levels and profiles of progestagen during pregnancy (6.44 +/- 0.37 microg/g) were consistent with those during nonpregnancy after copulation (6.34 +/- 0.49 microg/g). The duration of behavioral estrus in northern koalas was 13.5 +/- 0.9 days without copulation. In contrast, when estrous females mated, the estrous sign disappeared just after copulation. The mean (+/-SEM) length of the estrous cycle in northern koalas, as determined by behavioral estrus intervals, was 33.5 +/- 2.2 days without the luteal phase and 69.2 +/- 7.6 days with the luteal phase. Fecal progestagen analysis is a helpful and noninvasive tool to monitor ovulatory activity in northern and southern koalas and could help us to understand the reproductive activity of koalas by the combination approach with behavioral estrus.