Evidence for hybridization-driven heteroplasmy maintained across generations in a ricefish endemic to a Wallacean ancient lake.

Heteroplasmy, the presence of multiple mitochondrial DNA (mtDNA) haplotypes within cells of an individual, is caused by mutation or paternal leakage. However, heteroplasmy is usually resolved to homoplasmy within a few generations because of germ-line bottlenecks; therefore, instances of heteroplasmy are limited in nature. Here, we report heteroplasmy in the ricefish species Oryzias matanensis, endemic to Lake Matano, an ancient lake in Sulawesi Island, in which one individual was known to have many heterozygous sites in the mitochondrial NADH dehydrogenase subunit 2 (ND2) gene. In this study, we cloned the ND2 gene for some additional individuals with heterozygous sites and demonstrated that they are truly heteroplasmic. Phylogenetic analysis revealed that the extra haplotype within the heteroplasmic O. matanensis individuals clustered with haplotypes of O. marmoratus, a congeneric species inhabiting adjacent lakes. This indicated that the heteroplasmy originated from paternal leakage due to interspecific hybridization. The extra haplotype was unique and contained two non-synonymous substitutions. These findings demonstrate that this hybridization-driven heteroplasmy was maintained across generations for a long time to the extent that the extra mitochondria evolved within the new host.

Distribution and Population Genetic Structure of the Hau Giang Medaka, Oryzias haugiangensis, Along the East Coast of the Indochinese Peninsula.

The east coast of the Indochinese Peninsula is a well-known transition zone from subtropical to tropical systems, yet only a small number of studies have been conducted on the biogeography and phylogeography of aquatic organisms in this region. The Hau Giang medaka, Oryzias haugiangensis, was originally described from the Mekong Delta in southern Vietnam, and later reported also from southeastern Thailand, west of the Mekong Delta region. However, the species' full geographic range and population genetic structures remain unknown. Field surveys showed a widespread distribution of this species along the east coast of the Indochinese Peninsula, as far as northern Vietnam. A mitochondrial gene phylogeny and population genetic structure analysis using genome-wide single nucleotide polymorphisms revealed that the populations of O. haugiangensis are highly structuralized along the east coast of Vietnam, with the southernmost Mekong Delta population clearly separated from three populations north of central Vietnam. Further field collections are necessary to determine the boundary between the southern and northern populations, and the presence or absence of a hybrid zone.

A Cryptic Sex-Linked Locus Revealed by the Elimination of a Master Sex-Determining Locus in Medaka Fish.

AbstractSex chromosomes rapidly turn over in several taxonomic groups. Sex chromosome turnover is generally thought to start with the appearance of a new sex-determining gene on an autosome while an old sex-determining gene still exists, followed by the fixation of the new one. However, we do not know how prevalent the transient state is, where multiple sex-determining loci coexist within natural populations. Here, we removed a Y chromosome with a master male-determining gene DMY from medaka fish using high temperature-induced sex-reversed males. After four generations, the genomic characteristics of a sex chromosome were found on one chromosome, which was an autosome in the original population. Thus, the elimination of a master sex-determining locus can reveal a cryptic locus with a possible sex-determining effect, which can be the seed for sex chromosome turnover. Our results suggest that populations that seem to have a single-locus XY system may have other chromosomal regions with sex-determining effects. In conclusion, the coexistence of multiple sex-determining genes in a natural population may be more prevalent than previously thought. Experimental elimination of a master sex-determining locus may serve as a promising method for finding a locus that can be a protosex chromosome.

Multiple colonizations and hybridization of a freshwater fish group on a satellite island of Sulawesi.

Repeated colonizations and resultant hybridization may increase lineage diversity on an island if introgression occurs only in a portion of the indigenous island lineage. Therefore, to precisely understand how island biodiversity was shaped, it is essential to reconstruct the history of secondary colonization and resultant hybridization both in time and space. In this study, we reconstructed the history of multiple colonizations of the Oryzias woworae species group, a freshwater fish group of the family Adrianichthyidae, from Sulawesi Island to its southeast satellite island, Muna Island. Phylogenetic and species tree analyses using genome-wide single-nucleotide polymorphisms revealed that all local populations on Muna Island were monophyletic, but that there were several genetically distinct lineages within the island. Population structure and phylogenetic network analyses demonstrated that colonization of this island occurred more than once, and that secondary colonization and resultant introgressive hybridization occurred only in one local population on the island. The spatially heterogeneous introgression induced by the multiple colonizations were also supported by differential admixture analyses. In addition, the differential admixture analyses detected reverse colonization from Muna Island to the Sulawesi mainland. Coalescence-based demographic inference estimated that these mutual colonizations occurred during the middle to late Quaternary period, during which sea level repeatedly declined; this indicates that the colonizations occurred via land bridges. We conclude that these mutual colonizations between Muna Island and the Sulawesi mainland, and the resultant spatially heterogeneous introgression shaped the current biodiversity of this species group in this area.

Ghost introgression in ricefishes of the genus Adrianichthys in an ancient Wallacean lake.

Because speciation might have been promoted by ancient introgression from an extinct lineage, it is important to detect the existence of 'ghost introgression' in focal taxa and examine its contribution to their diversification. In this study, we examined possible ghost introgression and its contributions to the diversification of ricefishes of the genus Adrianichthys in Lake Poso, an ancient lake on Sulawesi Island, in which some extinctions are known to have occurred. Population-genomic analysis revealed that two extant Adrianichthys species, A. oophorus and A. poptae are reproductively isolated from each other. Comparisons of demographic models demonstrated that introgression from a ghost population, which diverged from the common ancestor of A. oophorus and A. poptae, is essential for reconstructing the demographic history of Adrianichthys. The best model estimated that the divergence of the ghost population greatly predated the divergence between A. oophorus and A. poptae, and that the ghost population secondarily contacted the two extant species within Lake Poso more recently. Genome scans and simulations detected a greatly divergent locus, which cannot be explained without ghost introgression. This locus was also completely segregated between A. oophorus and A. poptae. These findings suggest that variants that came from a ghost population have contributed to the divergence between A. oophorus and A. poptae, but the large time-lag between their divergence and ghost introgression indicates that the contribution of introgression may be restricted.

Genetic basis for the evolution of pelvic-fin brooding, a new mode of reproduction, in a Sulawesian fish.

Modes of reproduction in animals are diverse, with different modes having evolved independently in multiple lineages across a variety of taxa. However, an understanding of the genomic change driving the transition between different modes of reproduction is limited. Several ricefishes (Adrianichthyidae) on the island of Sulawesi have a unique mode of reproduction called "pelvic-fin brooding," wherein females carry externally fertilized eggs until hatching using their pelvic fins. Phylogenomic analysis demonstrated pelvic-fin brooders to have evolved at least twice in two distant clades of the Adrianichthyidae. We investigated the genetic architecture of the evolution of this unique mode of reproduction. Morphological analyses and laboratory observations revealed that females of pelvic-fin brooders have longer pelvic fins and a deeper abdominal concavity, and that they can carry an egg clutch for longer than nonbrooding adrianichthyids, suggesting that these traits play important roles in this reproductive mode. Quantitative trait locus mapping using a cross between a pelvic-fin brooder Oryzias eversi and a nonbrooding O. dopingdopingensis reveals different traits involved in pelvic-fin brooding to be controlled by different loci on different chromosomes. Genomic analyses of admixture detected no signatures of introgression between two lineages with pelvic-fin brooders, indicating that introgression is unlikely to be responsible for repeated evolution of pelvic-fin brooding. These findings suggest that multiple independent mutations may have contributed to the convergent evolution of this novel mode of reproduction.

Deeply divergent freshwater fish species within a single river system in central Sulawesi.

Sulawesi is a biodiversity hotspot for ricefishes (Adrianichthyidae), with many species endemic to the central part of this island in single ancient lakes or lake systems. Frequent vicariance by lake fragmentation since the Pliocene may be largely responsible for diversification in this family. In this study, we demonstrate that not only lacustrine species but also riverine species in this area are also deeply divergent even within a single river system. A mitochondrial phylogeny revealed that a ricefish population newly discovered from Cerekang River is sister to Oryzias dopingdopingensis Mandagi, Mokodongan, Tanaka, & Yamahira, another riverine species endemic to Doping-doping River. However, the Cerekang Oryzias was genetically isolated from O. dopingdopingensis, despite that Cerekang River and Doping-doping River share a connection across estuarine waters. This separation was supported by phylogenomic trees and population genetic structure analyses based on genome-wide single nucleotide polymorphisms. Coalescent-based demographic inference demonstrated that the ancestral population of these two riverine ricefishes had experienced a substantial population decrease and subsequently diverged into two sub-populations. Because the Cerekang Oryzias was also morphologically distinguished from O. dopingdopingensis, we described it as a new species, O. landangiensis. We infer that O. landangiensis and O. dopingdopingensis are of lake-origin and are relic species which were left in these rivers after the lake disappeared, and that they have lost their dispersal ability when inhabiting the ancient lake. The lost dispersal ability possibly contributed to the formation of the biodiversity hotspot for this fish group on this island.

Diversity of sex chromosomes in Sulawesian medaka fishes.

Recent genetic and genomic studies have revealed tremendous diversity in sex chromosomes across diverse taxa. Closely related species with different sex chromosomes provide us excellent opportunities to investigate the driving forces and the consequences of sex chromosome turnover. In the present study, we investigated the diversity of sex chromosomes of 13 Oryzias species from Sulawesi, Indonesia, which diversified during the last 4.86 million years. Using pooled sequencing, we found sex chromosomes in nine species that all had XY systems, with a species being possibly modified by multiple loci. Seven species (O. woworae, O. asinua, O. wolasi, O. matanensis, O. celebensis, O. hadiatyae, and O. dopingdopingensis) share linkage group (LG) 24 as sex chromosomes; however, they differed in the length and magnitude of sequence divergence between the X and Y chromosomes. The sex chromosome of O. eversi was LG4, which has not been reported as a sex chromosome in any other medaka species. In O. sarasinorum, LG16 and LG22 are associated with sex. Although LG16 was found to be sex-linked in another medaka species previously examined, the sex-determining regions did not overlap. No significant signatures for sex chromosomes were identified in the other four species (O. marmoratus, O. nigrimas, O. nebulosus, and O. orthognathus). Frequent turnovers and the great diversity of the sex chromosomes will make Sulawesian medaka species a model system for investigating the driving forces and consequences of sex chromosome turnover.

Ancient Admixture in Freshwater Halfbeaks of the Genus Nomorhamphus in Southeast Sulawesi.

Freshwater halfbeaks of the genus Nomorhamphus (Zenarchopteridae) uniquely diversified on Sulawesi Island, where tectonic movements have been very active since the Pliocene. Most species of this genus have quite limited distributions, which indicates that geographic isolations have contributed to their diversification. In this study, we demonstrated that secondary contacts and resultant admixtures between long-isolated species/populations may have also been important. We found that the mitochondrial phylogeny of a group of Nomorhamphus in Southeast Sulawesi was discordant with the nuclear phylogeny. Most notably, individuals in the upper and lower streams of the Moramo River, a small river in this region, clustered with each other in the mitochondrial phylogeny but not in the nuclear phylogeny; in the latter, the lower-stream individuals formed a clade with individuals in the Anduna River, a different river with no present water connection to the Moramo River. Phylogenetic network and population structure analyses using genomic data obtained from RNA-seq revealed that the lower-stream Moramo population admixed with the upper-stream Moramo lineage in ancient times. These findings indicate that the observed mito-nuclear discordance was caused by mitochondrial introgression and not incomplete lineage sorting. The phylogenetic network also revealed several other admixtures between ancient lineages. Repeated admixtures were also evidenced by topological incongruence in population trees estimated using the RNA-seq data. We propose that activities of many fault systems dissecting Southeast Sulawesi caused repeated secondary contact.

Resource partitioning is not coupled with assortative mating in sympatrically divergent ricefish in a Wallacean ancient lake.

Sympatric speciation is considered to be difficult without the coupling between ecological traits that allow resource partitioning and reproductive traits that allow assortative mating. Such "magic traits" are known to be involved in most of the compelling examples of sympatric speciation. In this study, we report a possible case of sympatric speciation without magic traits. Three species of ricefish (genus Oryzias) are suggested to have diverged sympatrically within Lake Poso, an ancient lake in Sulawesi. An analysis of genome-wide single-nucleotide polymorphisms showed that these three species are reproductively isolated from each other throughout the lake. Stable isotope analyses revealed that the three species use different food resources, which reflect differences in their feeding morphologies (gill rakers and digestive tracts) and feeding sites. Field and laboratory observations showed that O. nebulosus and O. orthognathus share a mating habitat of cobbles, where they scatter fertilized eggs, whereas this site is never used by O. nigrimas, indicating that assortative mating is partly achieved by spatial isolation. The small, less-adhesive eggs of O. nebulosus and O. orthognathus probably reflect their adaptation to spawning on cobble beaches. Laboratory mating experiments showed strong prezygotic isolation between O. nebulosus and O. orthognathus, which is achieved by strong species recognition presumably by both sexes based on species-specific mating dances and nuptial coloration. In summary, the assortative mating of O. nebulosus and O. orthognathus is probably not coupled to resource partitioning. We discussed how sympatric speciation among these species might have been achieved even without magic traits.

Species divergence and repeated ancient hybridization in a Sulawesian lake system.

An increasing volume of empirical studies demonstrated that hybridization between distant lineages may have promoted speciation in various taxa. However, the timing, extent and direction of introgressive hybridization remain unknown in many cases. Here, we report a possible case in which repeated hybridization promoted divergence of Oryzias ricefishes (Adrianichthyidae) on Sulawesi, an island of Wallacea. Four Oryzias species are endemic to the Malili Lake system in central Sulawesi, which is composed of five tectonic lakes; of these, one lake is inhabited by two species. Morphological and population genomic analyses of genome-wide single-nucleotide polymorphisms revealed that these two sympatric species are phylogenetically sister to but substantially reproductively isolated from each other. Analyses of admixture and comparison of demographic models revealed that the two sympatric species experienced several substantial introgressions from outgroup populations that probably occurred soon after they had secondary contact with each other in the lake. However, the ratio of migrants from the outgroups was estimated to be different between the two species, which is consistent with the hypothesis that these introgressions aided their divergence or prevented them from forming a hybrid swarm. Repeated lake fragmentations and fusions may have promoted diversification of this freshwater fish species complex that is endemic to this ancient lake system.

Diversity of lateral line patterns and neuromast numbers in the genus Oryzias.

A remarkable diversity of lateral line patterns exists in adult teleost fishes, the basis of which is largely unknown. By analysing the lateral line patterns and organ numbers in 29 Oryzias species and strains we report a rapid diversification of the lateral line system within this genus. We show a strong dependence of lateral line elaboration (number of neuromasts per cluster, number of parallel lateral lines) on adult species body size irrespective of phylogenetic relationships. In addition, we report that the degree of elaboration of the anterior lateral line, posterior lateral line and caudal neuromast clusters is tightly linked within species, arguing for a globally coordinated mechanism controlling lateral line organ numbers and patterns. We provide evidence for a polygenic control over neuromast numbers and positioning in the genus Oryzias. Our data also indicate that the diversity in lateral lines can arise as a result of differences in patterning both during embryonic development and post-embryonically, where simpler embryonic patterns generate less complex adult patterns and organ numbers, arguing for a linkage between the two processes.

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.

Variation in Mating Behaviors Between a Tropical and a Temperate Species of Medaka Fishes.

In medaka fishes, the family Adrianichthyidae, tropical species are known to be sexually more dimorphic than temperate species. If this reflects that tropical species are exposed to stronger sexual-selection pressures than temperate species, mating behaviors may also differ between tropical and temperate species. Our mating experiments revealed that males of Oryzias woworae, a tropical species of the family, perform "chasing" another male more frequently than males of O. sakaizumii, a temperate congener, and that male-male "combats" of O. woworae tended to be followed by chasing compared with combats of O. sakaizumii males, indicating that O. woworae males are more aggressive in male-male interactions than O. sakaizumii males. Males of O. woworae also performed "approaching" a female and "mating dance" more frequently than O. sakaizumii males, indicating that O. woworae males are also more active in courting females. Males of O. sakaizumii often omitted "mating dance" in their mating sequences, supporting this view. Moreover, O. woworae females tended to reject male "wrapping", an attempt for fertilization, more frequently than O. sakaizumii females, suggesting that O. woworae females are choosier in mating than O. sakaizumii females. These findings are concordant with the view that O. woworae is exposed to stronger sexual-selection pressures than O. sakaizumii.

Phylogenomics reveals habitat-associated body shape divergence in Oryzias woworae species group (Teleostei: Adrianichthyidae).

The Oryzias woworae species group, composed of O. asinua, O. wolasi, and O. woworae, is widely distributed in southeastern Sulawesi, an island in the Indo-Australian Archipelago. Deep-elongated body shape divergence is evident among these three species to the extent that it is used as a species-diagnostic character. These fishes inhabit a variety of habitats, ranging from upper streams to ponds, suggesting that the body shape divergence among the three species may reflect adaptation to local environments. First, our geometric morphometrics among eight local populations of this species group revealed that the three species cannot be separated by body shape and that riverine populations had more elongated bodies and longer caudal parts than lacustrine populations. Second, their phylogenetic relationships did not support the presence of three species; phylogenies using mitochondrial DNA and genomic data obtained from RNA-Seq revealed that the eight populations could not be sorted into three different clades representing three described species. Third, phylogenetic corrections of body shape variations and ancestral state reconstruction of body shapes demonstrated that body shape divergence between riverine and lacustrine populations persisted even if the phylogenies were considered and that body shape evolved rapidly irrespective of phylogeny. Sexual dimorphism in body shape was also evident, but the degree of dimorphism did not significantly differ between riverine and lacustrine populations after phylogenetic corrections, suggesting that sexual selection may not substantially contribute to geographical variations in body shape. Overall, these results indicate that the deep-elongated body shape divergence of the O. woworae species group evolved locally in response to habitat environments, such as water currents, and that a thorough taxonomic reexamination of the O. woworae species group may be necessary.

Origin and intra-island diversification of Sulawesi endemic Adrianichthyidae.

Although the family Adrianichthyidae is broadly distributed throughout East and Southeast Asia, 19 endemic species, over half of the family, are distributed in Sulawesi, which is an island in Wallacea. However, it remains unclear how this Adrianichthyidae biodiversity hotspot was shaped. In this study, we reconstructed molecular phylogenies for the Sulawesi adrianichthyids and estimated the divergence times of major lineages to infer the detailed history of their origin and subsequent intra-island diversification. The mitochondrial and nuclear phylogenies revealed that Sulawesi adrianichthyids are monophyletic, which indicates that they diverged from a single common ancestor. Species in the earliest branching lineages are currently distributed in the central and southeastern parts of the island, indicating that the common ancestor colonized Sula Spur, which is a large promontory that projects from the Australian continental margin, from Asia by oversea dispersal c.a. 20Mya. The first diversification event on Sulawesi, the split of the genus Adrianichthys, occurred c.a. 16Mya, and resulted in the nesting of Adrianichthys within Oryzias. Strong geographic structure was evident in the phylogeny; many species in the lineages branching off early are riverine and widely distributed in the southeastern and southwestern arms of Sulawesi, which suggests that oversea dispersal between tectonic subdivisions of this island during the late Miocene (7-5Mya) contributed to the distributions and diversification of the early branching lineages. In contrast, most species in the lineages branched off later are endemic to a single tectonic lake or lake system in the central Sulawesi, suggesting that habitat fragmentation due to the Pliocene collisions (c.a. 4Mya) among the tectonic subdivisions was the primary factor for diversification of the late branching, lacustrine lineages. Adrianichthys and some Oryzias in a certain late branching lineage are sympatric in Lake Poso, which indicates multiple colonizations of these distinct lineages into this tectonic lake. Thus, the diversification of Sulawesi adrianichthyids largely reflects the complex geological history of this island.

Discrepancy in the degree of population differentiation between color-morph frequencies and neutral genetic loci in the damselfly Ischnura senegalensis in Okinawa Island, Japan.

Evaluation of relative contribution of natural selection and stochastic processes to population differentiation has been of great interest in evolutionary biology. In a damselfly, Ischnura senegalensis, females show color dimorphism (gynochrome vs. androchrome), and color-morph frequencies are known to greatly vary among local populations within Okinawa Island, a small island of Ryukyu Archipelago, Japan. In this study, to examine the effects of natural selection and stochastic processes on the within-island variation in color-morph frequencies, we compared the degree of population differentiation at the color-morph locus with that at a mitochondrial DNA region and ten nuclear microsatellite loci. F ST values at the neutral loci were close to zero, indicating presence of sufficient gene flow (dispersal of adult individuals) between the local populations. In contrast, F ST values at the color-morph locus were significantly different from zero. These results suggest that variation in female color-morph frequencies observed among local populations in Okinawa Island has been caused by divergent selection acting on the phenotype and/or genes tightly linked with the color locus.

Ontogenetic stage-specific quantitative trait loci contribute to divergence in developmental trajectories of sexually dimorphic fins between medaka populations.

Sexual dimorphism can evolve when males and females differ in phenotypic optima. Genetic constraints can, however, limit the evolution of sexual dimorphism. One possible constraint is derived from alleles expressed in both sexes. Because males and females share most of their genome, shared alleles with different fitness effects between sexes are faced with intralocus sexual conflict. Another potential constraint is derived from genetic correlations between developmental stages. Sexually dimorphic traits are often favoured at adult stages, but selected against as juvenile, so developmental decoupling of traits between ontogenetic stages may be necessary for the evolution of sexual dimorphism in adults. Resolving intralocus conflicts between sexes and ages is therefore a key to the evolution of age-specific expression of sexual dimorphism. We investigated the genetic architecture of divergence in the ontogeny of sexual dimorphism between two populations of the Japanese medaka (Oryzias latipes) that differ in the magnitude of dimorphism in anal and dorsal fin length. Quantitative trait loci (QTL) mapping revealed that few QTL had consistent effects throughout ontogenetic stages and the majority of QTL change the sizes and directions of effects on fin growth rates during ontogeny. We also found that most QTL were sex-specific, suggesting that intralocus sexual conflict is almost resolved. Our results indicate that sex- and age-specific QTL enable the populations to achieve optimal developmental trajectories of sexually dimorphic traits in response to complex natural and sexual selection.

Female mate preference for longer fins in medaka.

Medaka, Oryzias latipes complex, display sexual dimorphisms in anal- and dorsal-fin lengths that suggest that females may prefer males with longer fins. However, female preference for longer anal and/or dorsal fins has not yet been described for the medaka. One reason that previous studies have not investigated this relationship may be because variations in male fin lengths within a single population are too small to experimentally detect female preference. In this study, we artificially crossed individuals from two wild populations (Aomori and Okinawa) that differed in male anal- and dorsal-fin lengths to increase phenotypic variation. We then tested female mate preference using these hybrid males. The results of the mating experiments and stepwise multiple regression analyses indicate that anal- and/or dorsal-fin lengths of the males contributed to female preference (i.e., males with longer anal and/or dorsal fins were less likely to be rejected by females). Variation in male standard length did not affect female preference. The evolution of female preference for longer fins in the medaka species complex may be explained by the "sexy son" hypothesis or the direct benefit hypothesis.