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.

Migration Behavior of Anguilla celebesensis Silver Eels within their Tomini Bay Spawning Area.

The tropical Celebes eel, Anguilla celebesensis, has a short migration between its spawning and growth habitats. Its spawning areas were hypothesized to be in Tomini Bay and the Celebes Sea after collecting their small leptocephali. However, there is no information about the silver eel oceanic spawning migration behavior of A. celebesensis. To better understand their short-distance spawning migration behavior, four large female silver eels (Eel 1-4) were equipped with pop-up satellite archival tags (PSATs) and released near the mouth of the Poso River in Tomini Bay of Sulawesi Island on 22 February (Eel 1-3) and 11 March 2010 (Eel 4). All PSATs ascended in Tomini Bay and transmitted their data. Eel 3 and 4 provided clear records of consistent diel vertical migration (DVM: eight days-Eel 3, 13 days-Eel 4) with daytime dives to mean depths of 444.7 m (Eel 3) and 539.0 m (Eel 4), where mean temperatures were 9.1°C (Eel 3) and 7.7°C (Eel 4), and nighttime ascents to mean depths of 132.8 m (Eel 3) and 112.4 m (Eel 4), where mean temperatures were 20.6°C (Eel 3) and 23.4°C (Eel 4). Eel 3 and 4 started to dive to deeper water around nautical dawn and swam up to shallower water around sunset. During nighttime, both eels swam in deeper and colder water during nights with moonlight than during nights without moonlight, and there was a negative linear relationship between experienced water temperatures with the moon in the sky and the lunar age for the eels. The A. celebesensis daily rhythm of DVM behaviors was similar to spawning-migration DVM behaviors of other anguillid species. Essential life history characteristics of A. celebesensis appear to be a short migration between freshwater growth habitat and ocean spawning habitat, and high GSI values with advanced gonadal development in downstream-migrating silver eels.

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.

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.

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.

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.

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 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.

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.

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.

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.

Mitochondrial introgression by ancient admixture between two distant lacustrine fishes in Sulawesi Island.

Sulawesi, an island located in a biogeographical transition zone between Indomalaya and Australasia, is famous for its high levels of endemism. Ricefishes (family Adrianichthyidae) are an example of taxa that have uniquely diversified on this island. It was demonstrated that habitat fragmentation due to the Pliocene juxtaposition among tectonic subdivisions of this island was the primary factor that promoted their divergence; however, it is also equally probable that habitat fusions and resultant admixtures between phylogenetically distant species may have frequently occurred. Previous studies revealed that some individuals of Oryzias sarasinorum endemic to a tectonic lake in central Sulawesi have mitochondrial haplotypes that are similar to the haplotypes of O. eversi, which is a phylogenetically related but geologically distant (ca. 190 km apart) adrianichthyid endemic to a small fountain. In this study, we tested if this reflects ancient admixture of O. eversi and O. sarasinorum. Population genomic analyses of genome-wide single-nucleotide polymorphisms revealed that O. eversi and O. sarasinorum are substantially reproductively isolated from each other. Comparison of demographic models revealed that the models assuming ancient admixture from O. eversi to O. sarasinorum was more supported than the models assuming no admixture; this supported the idea that the O. eversi-like mitochondrial haplotype in O. sarasinorum was introgressed from O. eversi. This study is the first to demonstrate ancient admixture of lacustrine or pond organisms in Sulawesi beyond 100 km. The complex geological history of this island enabled such island-wide admixture of lacustrine organisms, which usually experience limited migration.

Genome editing reveals fitness effects of a gene for sexual dichromatism in Sulawesian fishes.

Sexual selection drives rapid phenotypic diversification of mating traits. However, we know little about the causative genes underlying divergence in sexually selected traits. Here, we investigate the genetic basis of male mating trait diversification in the medaka fishes (genus Oryzias) from Sulawesi, Indonesia. Using linkage mapping, transcriptome analysis, and genome editing, we identify csf1 as a causative gene for red pectoral fins that are unique to male Oryzias woworae. A cis-regulatory mutation enables androgen-induced expression of csf1 in male fins. csf1-knockout males have reduced red coloration and require longer for mating, suggesting that coloration can contribute to male reproductive success. Contrary to expectations, non-red males are more attractive to a predatory fish than are red males. Our results demonstrate that integrating genomics with genome editing enables us to identify causative genes underlying sexually selected traits and provides a new avenue for testing theories of sexual selection.

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.

Evidence for sympatric speciation in a Wallacean ancient lake.

Sympatric speciation has been demonstrated in few empirical case studies, despite intense searches, because of difficulties in testing the criteria for this mode of speciation. Here, we report a possible case of sympatric speciation in ricefishes of the genus Oryzias on Sulawesi, an island of Wallacea. Three species of Oryzias are known to be endemic to Lake Poso, an ancient tectonic lake in central Sulawesi. Phylogenetic analyses using RAD-seq-derived single nucleotide polymorphisms (SNPs) revealed that these species are monophyletic. We also found that the three species are morphologically distinguishable and clearly separated by population-structure analyses based on the SNPs, suggesting that they are reproductively isolated from each other. A mitochondrial DNA chronogram suggested that their speciation events occurred after formation of the tectonic lake, and existence of a historical allopatric phase was not supported by coalescent-based demographic inference. Demographic inference also suggested introgressive hybridization from an outgroup population. However, differential admixture among the sympatric species was not supported by any statistical tests. These results all concur with criteria necessary to demonstrate sympatric speciation. Ricefishes in this Wallacean lake provide a promising new model system for the study of sympatric speciation.

Convergent evolution of body color between sympatric freshwater fishes via different visual sensory evolution.

Although there are many examples of color evolution potentially driven by sensory drive, only few studies have examined whether distinct species inhabiting the same environments evolve similar body colors via shared sensory mechanisms. In this study, we tested whether two sympatric freshwater fish taxa, halfbeaks of the genus Nomorhamphus and ricefishes of the genus Oryzias in Sulawesi Island, converge in both body color and visual sensitivity. After reconstructing the phylogeny separately for Nomorhamphus and Oryzias using transcriptome-wide sequences, we demonstrated positive correlations of body redness between these two taxa across environments, even after phylogenetic corrections, which support convergent evolution. However, substantial differences were observed in the expression profiles of opsin genes in the eyes between Nomorhamphus and Oryzias. Particularly, the expression levels of the long wavelength-sensitive genes were negatively correlated between the taxa, indicating that they have different visual sensitivities despite living in similar light environments. Thus, the convergence of body colorations between these two freshwater fish taxa was not accompanied by convergence in opsin sensitivities. This system presents a case in which body color convergence can occur between sympatric species via different mechanisms.

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.

Coelacanth genomes reveal signatures for evolutionary transition from water to land.

Coelacanths are known as "living fossils," as they show remarkable morphological resemblance to the fossil record and belong to the most primitive lineage of living Sarcopterygii (lobe-finned fishes and tetrapods). Coelacanths may be key to elucidating the tempo and mode of evolution from fish to tetrapods. Here, we report the genome sequences of five coelacanths, including four Latimeria chalumnae individuals (three specimens from Tanzania and one from Comoros) and one L. menadoensis individual from Indonesia. These sequences cover two African breeding populations and two known extant coelacanth species. The genome is ∼2.74 Gbp and contains a high proportion (∼60%) of repetitive elements. The genetic diversity among the individuals was extremely low, suggesting a small population size and/or a slow rate of evolution. We found a substantial number of genes that encode olfactory and pheromone receptors with features characteristic of tetrapod receptors for the detection of airborne ligands. We also found that limb enhancers of bmp7 and gli3, both of which are essential for limb formation, are conserved between coelacanth and tetrapods, but not ray-finned fishes. We expect that some tetrapod-like genes may have existed early in the evolution of primitive Sarcopterygii and were later co-opted to adapt to terrestrial environments. These coelacanth genomes will provide a cornerstone for studies to elucidate how ancestral aquatic vertebrates evolved into terrestrial animals.