Gsdf is not indispensable for male differentiation in the medaka species Oryzias hubbsi.

Sex determination mechanisms differ widely among vertebrates, particularly in fish species, where diverse sex chromosomes and sex-determining genes have evolved. However, the sex-differentiation pathways activated by these sex-determining genes appear to be conserved. Gonadal soma-derived growth factor (Gsdf) is one of the genes conserved across teleost fish, especially in medaka fishes of the genus Oryzias, and is implicated in testis differentiation and germ cell proliferation. However, its role in sex differentiation remains unclear. In this study, we investigated Gsdf function in Oryzias hubbsi, a species with a ZW sex-determination system. We confirmed its male-dominant expression, as in other species. However, histological analyses revealed no male-to-female sex reversal in Gsdf-knockout fish, contrary to findings in other medaka species. Genetic sex determination remained intact without Gsdf function, indicating a Gsdf-independent sex-differentiation pathway in O. hubbsi. Instead, Gsdf loss led to germ cell overproliferation in both sexes and accelerated onset of meiosis in testes, suggesting a role in germ cell proliferation. Notably, the feminizing effect of germ cells observed in O. latipes was absent, suggesting diverse germ cell-somatic cell relationships in Oryzias gonad development. Our study highlights species-specific variations in the molecular pathways governing sex determination and differentiation, emphasizing the need for further exploration to elucidate the complexities of sexual development.

Mesozoic origin and 'out-of-India' radiation of ricefishes (Adrianichthyidae).

The Indian subcontinent has an origin geologically different from Eurasia, but many terrestrial animal and plant species on it have congeneric or sister species in other parts of Asia, especially in the Southeast. This faunal and floral similarity between India and Southeast Asia is explained by either of the two biogeographic scenarios, 'into-India' or 'out-of-India'. Phylogenies based on complete mitochondrial genomes and five nuclear genes were undertaken for ricefishes (Adrianichthyidae) to examine which of these two biogeographic scenarios fits better. We found that Oryzias setnai, the only adrianichthyid distributed in and endemic to the Western Ghats, a mountain range running parallel to the western coast of the Indian subcontinent, is sister to all other adrianichthyids from eastern India and Southeast-East Asia. Divergence time estimates and ancestral area reconstructions reveal that this western Indian species diverged in the late Mesozoic during the northward drift of the Indian subcontinent. These findings indicate that adrianichthyids dispersed eastward 'out-of-India' after the collision of the Indian subcontinent with Eurasia, and subsequently diversified in Southeast-East Asia. A review of geographic distributions of 'out-of-India' taxa reveals that they may have largely fuelled or modified the biodiversity of Eurasia.

Origin of Boundary Populations in Medaka (Oryzias latipes Species Complex).

The Japanese wild population of the medaka fish (Oryzias latipes species complex) comprises two genetically distinct groups, the Northern and the Southern Populations, with boundary populations having a unique genotype. It is thought that the boundary populations have been formed through introgressive hybridization between the two groups, because they are fixed with the Northern alleles at two allozymic loci, with the Southern alleles at two other loci, and have a unique allele at one locus. In this study, we examined the genetic population structure of the boundary populations using genome-wide single nucleotide polymorphism (SNP) data. Most SNPs of the Toyooka population, a typical boundary population, were shared with the Northern Population, some were shared with the Southern Population, and the remaining SNPs were unique to this population, suggesting that the boundary populations originated and diverged from the Northern Population. Further analyses of different populations using SNPs at eight genomic loci indicated that the boundary populations at different locations share similar genomic constitutions, and can be genetically distinguished from typical Northern Populations by unique SNPs. In addition, the boundary populations in the Maruyama River Basin had Northern mitochondrial DNA (mtDNA), while others, from the Fukuda and Kishida River Basins and from the Kumihama Bay area, had Southern mtDNA. These findings suggested that the boundary populations originated from the Northern Population, and then their genomes diverged as a result of geographical isolation, followed by mtDNA introgression from the Southern Population that occurred independently in some populations.

Diversity and Convergence of Sex-Determination Mechanisms in Teleost Fish.

Sexual reproduction is prevalent across diverse taxa. However, sex-determination mechanisms are so diverse that even closely related species often differ in sex-determination systems. Teleost fish is a taxonomic group with frequent turnovers of sex-determining mechanisms and thus provides us with great opportunities to investigate the molecular and evolutionary mechanisms underlying the turnover of sex-determining systems. Here, we compile recent studies on the diversity of sex-determination mechanisms in fish. We demonstrate that genes in the TGF-ß signaling pathway are frequently used for master sex-determining (MSD) genes. MSD genes arise via two main mechanisms, duplication-and-transposition and allelic mutations, with a few exceptions. We also demonstrate that temperature influences sex determination in many fish species, even those with sex chromosomes, with higher temperatures inducing differentiation into males in most cases. Finally, we review theoretical models for the turnover of sex-determining mechanisms and discuss what questions remain elusive.

Sub-functionalization of duplicated genes in the evolution of nine-spined stickleback hatching enzyme.

Gene duplication is the primary source of novel genes, and is followed by non-, sub-, or neo-functionalization. In this study, we compared the egg envelope digestion mechanism of hatching enzymes between three-spined stickleback and nine-spined stickleback species, and found that the function of the hatching enzymes of nine-spined sticklebacks was uniquely derived by gene duplication, followed by sub-functionalization. The hatching enzyme of euteleosts consists of two metalloproteases, high choriolytic enzyme (HCE), and low choriolytic enzyme (LCE). LCE, especially, has an important role in solubilizing egg envelope protein by cleaving two specific sites. Three-spined stickleback had a single copy of the LCE gene, like other euteleosts. However, nine-spined stickleback had two types of LCE genes, α-type and ß-type, suggesting that a duplication of the LCE gene occurred during the evolution of sticklebacks. The α-type and ß-type each cleaved one of the two sites. Therefore, in the nine-spined stickleback, the function of the ancestral LCE was driven by a single copy gene, which was partitioned into two functions separately driven by two duplicated genes, and egg envelope was solubilized by the cooperative action of the two LCEs, α-type and ß-type. Herein, we provide a molecular mechanism for an evolutionary adaptation driven by gene duplication and sub-functionalization.

Molecular cloning and characterization of the repetitive DNA sequences that comprise the constitutive heterochromatin of the W chromosomes of medaka fishes.

Among the medaka fishes of the genus Oryzias, most species have homomorphic sex chromosomes, while some species, such as Oryzias hubbsi and Oryzias javanicus, have heteromorphic ZW sex chromosomes. In this study, a novel family of repetitive sequence was molecularly cloned from O. hubbsi and characterized by chromosome in situ and filter hybridization, respectively. This repetitive element, which we designated as a BstNI family element, localized at heterochromatin regions on the W chromosome, as well as on two pairs of autosomes. Homologous sequences to this element were found only in O. javanicus, which is a sister species of O. hubbsi, suggesting that this repeated element originated in the common ancestor of these two species. However, the intensity of the hybridization signals was lower in O. javanicus than in O. hubbsi, and the chromosomal location of this element in O. javanicus was confined to heterochromatin regions on one pair of autosomes. Thus, we hypothesize that this repetitive element was extensively amplified in the O. hubbsi lineage, especially on its W chromosome, after the separation of the O. javanicus lineage. In addition, we also found the W chromosomal location of the 18S-28S ribosomal RNA genes in both O. hubbsi and O. javanicus. Our previous studies showed no linkage homology of the sex chromosomes in these species, indicating that the RNA genes were shared between W chromosomes of different origins. This situation may be explained by a translocation of the sex-determining region with the ribosomal RNA genes in either species or an independent accumulation of the RNA genes as a convergent process during W chromosome degeneration.

Establishment of cytochrome P450 1a gene-knockout Javanese medaka, Oryzias javanicus, which distinguishes toxicity modes of the polycyclic aromatic hydrocarbons, pyrene and phenanthrene.

Cytochrome P450 1a (Cyp1a) is an important enzyme for metabolism of organic pollutants. To understand its reaction to polycyclic aromatic hydrocarbons (PAHs), we knocked out this gene in a marine model fish, Javanese medaka, Oryzias javanicus, using the CRISPR/Cas 9 system. A homozygous mutant (KO) strain with a four-base deletion was established using an environmental DNA (eDNA)-based genotyping technique. Subsequently, KO, heterozygous mutant (HT), and wild-type (WT) fish were exposed to model pollutants, pyrene and phenanthrene, and survivorship and swimming behavior were analyzed. Compared to WT, KO fish were more sensitive to pyrene, suggesting that Cyp1a transforms pyrene into less toxic metabolites. Conversely, WT fish were sensitive to phenanthrene, suggesting that metabolites transformed by Cyp1a are more toxic than the original compound. HT fish showed intermediate results. Thus, comparative use of KO and WT fish can distinguish modes of pollutant toxicity, providing a deeper understanding of fish catabolism of environmental pollutants.

Effects of synthetic sex steroid hormone exposures on gonadal sex differentiation and dynamics of a male-related gene, Gonadal soma-derived factor (Gsdf) and an estrogen up-regulated gene, Choriogenine-H (ChgH) gene expression in the euryhaline Javafish medaka, Oryzias javanicus, based on genetic sexes.

To clarify the basal aspects of sexual development in Javafish medaka, Oryzias javanicus (ZZ/ZW), a model marine species for ecotoxicity testing, we examined the details of gonadal sex differentiation and exogenous sex hormone-dependent sex reversals using genetic sex-linked DNA markers. Sex differences in germ cell numbers were observed at 5 days post hatching (dph), in which there was a significant increase in the germ cells of ZW. In ZW, diplotene oocytes and the ovarian cavity appeared at approximately 10, and 30 dph, respectively. In ZZ, spermatogonial proliferation was observed at approximately 20 dph. A ZZ-dominant expression of Gonadal soma-derived factor (Gsdf) mRNA was detected before hatching. The exposure of embryos to 17α-ethinylestradiol (EE2; 0.1, 1, 10 ng/mL) did not cause sex reversals in most cases. However, EE2 exposures led to significant Choriogenin-H (ChgH) mRNA expression, an estrogen up-regulated gene, in all fry; these exposures did not suppress Gsdf expression in ZZ fry. The exposure of embryos to 17α-methyltestosterone (MT; 0.1, 1, 10 ng/mL) caused sex reversals but only at low frequencies in ZW and ZZ fish. Although the 10 ng/mL MT exposure was accompanied by induction of significant Gsdf expression in ZW fry, induction of ChgH expression was also observed in several fry. Together, the present study indicates for the first time that male-dominant sexual dimorphic expression of Gsdf precedes the first morphological sex difference, i.e., the sex difference in germ cell number, and results strongly suggest that exogenous sex hormone-dependent sex reversal is not induced easily in O. javanicus.

Genome Sequence of the Euryhaline Javafish Medaka, Oryzias javanicus: A Small Aquarium Fish Model for Studies on Adaptation to Salinity.

The genus Oryzias consists of 35 medaka-fish species each exhibiting various ecological, morphological and physiological peculiarities and adaptations. Beyond of being a comprehensive phylogenetic group for studying intra-genus evolution of several traits like sex determination, behavior, morphology or adaptation through comparative genomic approaches, all medaka species share many advantages of experimental model organisms including small size and short generation time, transparent embryos and genome editing tools for reverse and forward genetic studies. The Java medaka, Oryzias javanicus, is one of the two species of medaka perfectly adapted for living in brackish/sea-waters. Being an important component of the mangrove ecosystem, O. javanicus is also used as a valuable marine test-fish for ecotoxicology studies. Here, we sequenced and assembled the whole genome of O. javanicus, and anticipate this resource will be catalytic for a wide range of comparative genomic, phylogenetic and functional studies. Complementary sequencing approaches including long-read technology and data integration with a genetic map allowed the final assembly of 908 Mbp of the O. javanicus genome. Further analyses estimate that the O. javanicus genome contains 33% of repeat sequences and has a heterozygosity of 0.96%. The achieved draft assembly contains 525 scaffolds with a total length of 809.7 Mbp, a N50 of 6,3 Mbp and a L50 of 37 scaffolds. We identified 21454 predicted transcripts for a total transcriptome size of 57, 146, 583 bps. We provide here a high-quality chromosome scale draft genome assembly of the euryhaline Javafish medaka (321 scaffolds anchored on 24 chromosomes (representing 97.7% of the total bases)), and give emphasis on the evolutionary adaptation to salinity.

Construction of High-Resolution RAD-Seq Based Linkage Map, Anchoring Reference Genome, and QTL Mapping of the Sex Chromosome in the Marine Medaka Oryzias melastigma.

Medaka (Oryzias sp.) is an important fish species in ecotoxicology and considered as a model species due to its biological features including small body size and short generation time. Since Japanese medaka Oryzias latipes is a freshwater species with access to an excellent genome resource, the marine medaka Oryzias melastigma is also applicable for the marine ecotoxicology. In genome era, a high-density genetic linkage map is a very useful resource in genomic research, providing a means for comparative genomic analysis and verification of de novo genome assembly. In this study, we developed a high-density genetic linkage map for O. melastigma using restriction-site associated DNA sequencing (RAD-seq). The genetic map consisted of 24 linkage groups with 2,481 single nucleotide polymorphism (SNP) markers. The total map length was 1,784 cM with an average marker space of 0.72 cM. The genetic map was integrated with the reference-assisted chromosome assembly (RACA) of O. melastigma, which anchored 90.7% of the assembled sequence onto the linkage map. The values of complete Benchmarking Universal Single-Copy Orthologs were similar to RACA assembly but N50 (23.74 Mb; total genome length 779.4 Mb; gap 5.29%) increased to 29.99 Mb (total genome length 778.7 Mb; gap 5.2%). Using MapQTL analysis with SNP markers, we identified a major quantitative trait locus for sex traits on the Om10. The integration of the genetic map with the reference genome of marine medaka will serve as a good resource for studies in molecular toxicology, genomics, CRISPR/Cas9, and epigenetics.

Seasonal regulation of the lncRNA LDAIR modulates self-protective behaviours during the breeding season.

To cope with seasonal environmental changes, animals adapt their physiology and behaviour in response to photoperiod. However, the molecular mechanisms underlying these adaptive changes are not completely understood. Here, using genome-wide expression analysis, we show that an uncharacterized long noncoding RNA (lncRNA), LDAIR, is strongly regulated by photoperiod in Japanese medaka fish (Oryzias latipes). Numerous transcripts and signalling pathways are activated during the transition from short- to long-day conditions; however, LDAIR is one of the first genes to be induced and its expression shows a robust daily rhythm under long-day conditions. Transcriptome analysis of LDAIR knockout fish reveals that the LDAIR locus regulates a gene neighbourhood, including corticotropin releasing hormone receptor 2, which is involved in the stress response. Behavioural analysis of LDAIR knockout fish demonstrates that LDAIR affects self-protective behaviours under long-day conditions. Therefore, we propose that photoperiodic regulation of corticotropin releasing hormone receptor 2 by LDAIR modulates adaptive behaviours to seasonal environmental changes.

A key metabolic gene for recurrent freshwater colonization and radiation in fishes.

Colonization of new ecological niches has triggered large adaptive radiations. Although some lineages have made use of such opportunities, not all do so. The factors causing this variation among lineages are largely unknown. Here, we show that deficiency in docosahexaenoic acid (DHA), an essential ω-3 fatty acid, can constrain freshwater colonization by marine fishes. Our genomic analyses revealed multiple independent duplications of the fatty acid desaturase gene Fads2 in stickleback lineages that subsequently colonized and radiated in freshwater habitats, but not in close relatives that failed to colonize. Transgenic manipulation of Fads2 in marine stickleback increased their ability to synthesize DHA and survive on DHA-deficient diets. Multiple freshwater ray-finned fishes also show a convergent increase in Fads2 copies, indicating its key role in freshwater colonization.

Evolution of different Y chromosomes in two medaka species, Oryzias dancena and O. latipes.

Although the sex-determining gene DMY has been identified on the Y chromosome in the medaka (Oryzias latipes), this gene is absent in most Oryzias species, suggesting that closely related species have different sex-determining genes. Here, we investigated the sex-determination mechanism in O. dancena, which does not possess the DMY gene. Since heteromorphic sex chromosomes have not been reported in this species, a progeny test of sex-reversed individuals produced by hormone treatment was performed. Sex-reversed males yielded all-female progeny, indicating that O. dancena has an XX/XY sex-determination system. To uncover the cryptic sex chromosomes, sex-linked DNA markers were screened using expressed sequence tags (ESTs) established in O. latipes. Linkage analysis of isolated sex-linked ESTs showed a conserved synteny between the sex chromosomes in O. dancena and an autosome in O. latipes. Fluorescence in situ hybridization (FISH) analysis of these markers confirmed that sex chromosomes of these species are not homologous. These findings strongly suggest an independent origin of sex chromosomes in O. dancena and O. latipes. Further analysis of the sex-determining region in O. dancena should provide crucial insights into the evolution of sex-determination mechanisms in vertebrates.

Evidence for different origins of sex chromosomes in closely related Oryzias fishes: substitution of the master sex-determining gene.

The medaka Oryzias latipes and its two sister species, O. curvinotus and O. luzonensis, possess an XX-XY sex-determination system. The medaka sex-determining gene DMY has been identified on the orthologous Y chromosome [O. latipes linkage group 1 (LG1)] of O. curvinotus. However, DMY has not been discovered in other Oryzias species. These results and molecular phylogeny suggest that DMY was generated recently [approximately 10 million years ago (MYA)] by gene duplication of DMRT1 in a common ancestor of O. latipes and O. curvinotus. We identified seven sex-linked markers from O. luzonensis (sister species of O. curvinotus) and constructed a sex-linkage map. Surprisingly, all seven sex-linked markers were located on an autosomal linkage group (LG12) of O. latipes. As suggested by the phylogenetic tree, the sex chromosomes of O. luzonensis should be "younger" than those of O. latipes. In the lineage leading to O. luzonensis after separation from O. curvinotus approximately 5 MYA, a novel sex-determining gene may have arisen and substituted for DMY. Oryzias species should provide a useful model for evolution of the master sex-determining gene and differentiation of sex chromosomes from autosomes.

CRISPR/Cas9-based heritable targeted mutagenesis in Thermobia domestica: A genetic tool in an apterygote development model of wing evolution.

Despite previous developmental studies on basally branching wingless insects and crustaceans, the evolutionary origin of insect wings remains controversial. Knowledge regarding genetic regulation of tissues hypothesized to have given rise to wings would help to elucidate how ancestral development changed to allow the evolution of true wings. However, genetic tools available for basally branching wingless species are limited. The firebrat Thermobia domestica is an apterygote species, phylogenetically related to winged insects. T. domestica presents a suitable morphology to investigate the origin of wings, as it forms the tergal paranotum, from which wings are hypothesized to have originated. Here we report the first successful CRISPR/Cas9-based germline genome editing in T. domestica. We provide a technological platform to understand the development of tissues hypothesized to have given rise to wings in an insect with a pre-wing evolution body plan.

pnp4a Is the Causal Gene of the Medaka Iridophore Mutant guanineless.

See-through medaka lines are suitable for observing internal organs throughout life. They were bred by crossing multiple color mutants. However, some of the causal genes for these mutants have not been identified. The medaka has four pigment cell types: black melanophores, yellow xanthophores, white leucophores, and silvery iridophores. The causal genes of melanophore, xanthophore, and leucophore mutants have been elucidated, but the causal gene for the iridophore mutant remains unknown. Here, we describe the iridophore mutant, guanineless (gu), which exhibits a strong reduction in visible iridophores throughout its larval to adult stages. The gu locus was previously mapped to chromosome 5, but was located near the telomeric region, making it difficult to integrate into the chromosome. We sought the causal gene of gu using synteny analysis with the zebrafish genome and found a strong candidate, purine nucleoside phosphorylase 4a (pnp4a). Gene targeting and complementation testing showed that pnp4a is the causal gene of gu This result will allow the establishment of inbred medaka strains or other useful strains with see-through phenotypes without major disruption in the genetic background of each strain.

Evolution of Shh endoderm enhancers during morphological transition from ventral lungs to dorsal gas bladder.

Shh signalling plays a crucial role for endoderm development. A Shh endoderm enhancer, MACS1, is well conserved across terrestrial animals with lungs. Here, we first show that eliminating mouse MACS1 causes severe defects in laryngeal development, indicating that MACS1-directed Shh signalling is indispensable for respiratory organogenesis. Extensive phylogenetic analyses revealed that MACS1 emerged prior to the divergence of cartilaginous and bony fishes, and even euteleost fishes have a MACS1 orthologue. Meanwhile, ray-finned fishes evolved a novel conserved non-coding sequence in the neighbouring region. Transgenic assays showed that MACS1 drives reporter expression ventrally in laryngeal epithelium. This activity has been lost in the euteleost lineage, and instead, the conserved non-coding sequence of euteleosts acquired an enhancer activity to elicit dorsal epithelial expression in the posterior pharynx and oesophagus. These results implicate that evolution of these two enhancers is relevant to the morphological transition from ventral lungs to dorsal gas bladder.

Arginine vasotocin neuronal development and its projection in the adult brain of the medaka.

The neurohypophysial peptide arginine vasotocin (AVT) and its mammalian ortholog arginine vasopressin function in a wide range of physiological and behavioral events. Here, we generated a new line of transgenic medaka (Oryzias latipes), which allowed us to monitor AVT neurons by enhanced green fluorescent protein (EGFP) and demonstrate AVT neuronal development in the embryo and the projection of AVT neurons in the adult brain of avt-egfp transgenic medaka. The onset of AVT expression manifested at 2 days postfertilization (dpf) as a pair of signals in the telencephalon of the brain. The telencephalic AVT neurons migrated and converged on the preoptic area (POA) by 4dpf. At the same stage, another onset of AVT expression manifested in the central optic tectum (OT), and they migrated to the ventral part of the hypothalamus (VH) by 6dpf. In the adult brain, the AVT somata with EGFP signals existed in the gigantocellular POA (gPOA), magnocellular POA (mPOA), and parvocellular POA (pPOA) and in the VH. Whereas the major projection of AVT fibers was found from the pPOA and VH to the posterior pituitary, it was also found that AVT neurons in the three POAs send their fibers into wide regions of the brain such as the telencephalon, mesencephalon and diencephalon. This study suggests that the avt-egfp transgenic medaka is a useful model to explore AVT neuronal development and function.

Turnover of Sex Chromosomes in Celebensis Group Medaka Fishes.

Sex chromosomes and the sex-determining (SD) gene are variable in vertebrates. In particular, medaka fishes in the genus Oryzias show an extremely large diversity in sex chromosomes and the SD gene, providing a good model to study the evolutionary process by which they turnover. Here, we investigated the sex determination system and sex chromosomes in six celebensis group species. Our sex-linkage analysis demonstrated that all species had an XX-XY sex determination system, and that the Oryzias marmoratus and O. profundicola sex chromosomes were homologous to O. latipes linkage group (LG) 10, while those of the other four species, O. celebensis, O. matanensis, O. wolasi, and O. woworae, were homologous to O. latipes LG 24. The phylogenetic relationship suggested a turnover of the sex chromosomes from O. latipes LG 24 to LG 10 within this group. Six sex-linkage maps showed that the former two and the latter four species shared a common SD locus, respectively, suggesting that the LG 24 acquired the SD function in a common ancestor of the celebensis group, and that the LG 10 SD function appeared in a common ancestor of O. marmoratus and O. profundicola after the divergence of O. matanensis. Additionally, fine mapping and association analysis in the former two species revealed that Sox3 on the Y chromosome is a prime candidate for the SD gene, and that the Y-specific 430-bp insertion might be involved in its SD function.

Genetic structure of korean wild populations of the Medaka Oryzias latipes inferred from allozymic variation.

Previous allozymic studies have revealed that Korean wild populations of Oryzias latipes have differentiated regionally, and are composed of two distinct groups, the East Korean Population and the China-West Korean Population. Recently, mitochondrial DNA (mtDNA) sequencing and restriction fragment length polymorphism (RFLP) analyses have confirmed these two groups, and shown that the distribution ranges of the two groups overlap in western Korea. In order to describe the detailed distributions of the two groups and the gene flow between them, genotypes of 13 allozymic loci were determined in 444 specimens from 96 localities in Korea. The two major groups were supported by remarkable allele frequency differences at six diagnostic loci: ACP*, AMY*, CK-A*, LDH-A*, PGM* and TF*. Individuals with the typical "eastern" genotype were mainly distributed in eastern and southern areas. In contrast, fish with the "western" genotype were predominant in the western area, and were further divided into two subgroups (the Han River and Geum River Subpopulations) by unique alleles at the ADH* locus. In the western coast, two distinct (eastern and western) genotypes were distributed in a mosaic fashion. This distribution pattern was identical to those from mtDNA analyses. Although the distribution patterns of the alleles at three loci (GPI-A*, LDH-C* and SOD*) showed introgressive conditions between the two groups, each population was nearly fixed as either the eastern or western genotype at all six diagnostic loci despite the proximity among samples. Therefore, it is suggested that some reproductive isolation mechanisms exist between the two groups in natural habitats.