Origin and intra-island diversification of Sulawesi endemic Adrianichthyidae

The Medaka Approach to Evolutionary Social Neuroscience

Previously, the integration of comparative biological and neuroscientific approaches has led to significant advancements in social neuroscience. This review elucidates the potential and future directions of evolutionary social neuroscience research utilizing medaka fishes (the family Adrianichthyidae) including Japanese medaka (Oryzias latipes). We focus on medaka's social cognitive capabilities and mate choice behavior, particularly emphasizing mate preference using visual cues. Medaka fishes are also advantageous due to their abundant genetic resources, extensive genomic information, and the relative ease of laboratory breeding and genetic manipulation. Here we present some research examples of both the conventional neuroscience approach and evolutionary approach involving medaka fishes and other species. We also discuss the prospects of uncovering the molecular and cellular mechanisms underlying the diversity of visual mate preference among species. Especially, we introduce that the single-cell transcriptome technology, particularly in conjunction with 'Adaptive Circuitry Census', is an innovative tool that bridges comparative biological methods and neuroscientific approaches. Evolutionary social neuroscience research using medaka has the potential to unveil fundamental principles in neuroscience and elucidate the mechanisms responsible for generating diversity in mating strategies.

Rapid body colouration changes in Oryzias celebensis as a social signal influenced by environmental background

Rapid body colouration changes in some animals, such as chameleons and octopuses, serve dual functions: camouflage and intraspecific communication. It has been hypothesized that these colouration changes originally evolved to provide camouflage and subsequently were co-opted as social signals; however, experimental model systems that are suitable for studying such evolutionary processes are limited. Here, we investigated the relationship between rapid colouration changes of the blackened markings and aggressive behaviours in male Oryzias celebensis, an Indonesian medaka fish, under triadic relationships (two males and one female) or three males conditions with two different environmental backgrounds. In an algae-covered tank, mimicking the common laboratory rearing conditions, males with blackened markings exhibited more frequent attacks towards different conspecific individuals compared with non-blackened males and females. The blackened males were seldom attacked by non-blackened males and females. By contrast, neither aggressive behaviours nor black colouration changes were observed in the transparent background condition with a brighter environment. These indicated that the blackened markings in O. celebensis serve as a social signal depending on the environmental backgrounds. Considering that such colouration changes for camouflage are widely conserved among teleost fishes, the traits are likely to be co-opted for displaying social signals in O. celebensis.

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.

Gene Duplication of Androgen Receptor As An Evolutionary Driving Force Underlying the Diversity of Sexual Characteristics in Teleost Fishes

Sexual dimorphism allows species to meet their fitness optima based on the physiological availability of each sex. Although intralocus sexual conflict appears to be a genetic constraint for the evolution of sex-specific traits, sex-linked genes and the regulation of sex steroid hormones contribute to resolving this conflict by allowing sex-specific developments. Androgens and their receptor, androgen receptor (Ar), regulate male-biased phenotypes. In teleost fish, ar ohnologs have emerged as a result of teleost-specific whole genome duplication (TSGD). Recent studies have highlighted the evolutionary differentiation of ar ohnologs responsible for the development of sexual characteristics, which sheds light on the need for comparative studies on androgen regulation among different species. In this review, we discuss the importance of ar signaling as a regulator of male-specific traits in teleost species because teleost species are suitable experimental models for comparative studies owing to their great diversity in male-biased morphological and physiological traits. To date, both in vivo and in vitro studies on teleost ar ohnologs have shown a substantial influence of ars as a regulator of male-specific reproductive traits such as fin elongation, courtship behavior, and nuptial coloration. In addition to these sexual characteristics, ar substantially influences immunity, inducing a sex-biased immune response. This review aims to provide a comprehensive understanding of the current state of teleost ar studies and emphasizes the potential of teleost fishes, given their availability, to find molecular evidence about what gives rise to the spectacular diversity among fish species.

A landmark-free analysis of the pelvic girdle in Sulawesi ricefishes (Adrianichthyidae): How 2D and 3D geometric morphometrics can complement each other in the analysis of a complex structure

Geometric morphometrics (GM) enable the quantification of morphological variation on various scales. Recent technical advances allow analyzing complex three-dimensional shapes also in cases where landmark-based approaches are not appropriate. Pelvic girdle bones (basipterygia) of Sulawesi ricefishes are 3D structures that challenge traditional morphometrics. We hypothesize that the pelvic girdle of ricefishes experienced sex-biased selection pressures in species where females provide brood care by carrying fertilized eggs supported by elongated pelvic fins ("pelvic brooding"). We test this by comparing pelvic bone shapes of both sexes in species exhibiting pelvic brooding and the more common reproductive strategy "transfer brooding," by using landmark-free 2D and 3D GM, as well as qualitative shape descriptions. Both landmark-free approaches revealed significant interspecific pelvic bone variation in the lateral process, medial facing side of the pelvic bone, and overall external and internal wing shape. Within pelvic brooders, the three analyzed species are clearly distinct, while pelvic bones of the genus Adrianichthys are more similar to transfer brooding Oryzias. Female pelvic brooding Oryzias exhibit prominent, medially pointing tips extending from the internal wing and basipterygial plate that are reduced or absent in conspecific males, Adrianichthys and transfer brooding Oryzias, supporting our hypothesis that selection pressures affecting pelvic girdle shape are sex-biased in Sulawesi ricefishes. Furthermore, both sexes of pelvic brooding Oryzias have overall larger pelvic bones than other investigated ricefishes. Based on these differences, we characterized two reproductive strategy- and sex-dependent pelvic girdle types for Sulawesi ricefishes. Morphological differences between the investigated pelvic brooding genera Adrianichthys and Oryzias provide additional evidence for two independent origins of pelvic brooding. Overall, our findings add to a better understanding on traits related to pelvic brooding in ricefishes and provide a basis for upcoming studies on pelvic girdle function and morphology.

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.

Species Delimitation, Phylogenomics, and Biogeography of Sulawesi Flying Lizards: A Diversification History Complicated by Ancient Hybridization, Cryptic Species, and Arrested Speciation

The biota of Sulawesi is noted for its high degree of endemism and for its substantial levels of in situ biological diversification. While the island's long period of isolation and dynamic tectonic history have been implicated as drivers of the regional diversification, this has rarely been tested in the context of an explicit geological framework. Here, we provide a tectonically informed biogeographical framework that we use to explore the diversification history of Sulawesi flying lizards (the Draco lineatus Group), a radiation that is endemic to Sulawesi and its surrounding islands. We employ a framework for inferring cryptic speciation that involves phylogeographic and genetic clustering analyses as a means of identifying potential species followed by population demographic assessment of divergence-timing and rates of bi-directional migration as means of confirming lineage independence (and thus species status). Using this approach, phylogenetic and population genetic analyses of mitochondrial sequence data obtained for 613 samples, a 50-SNP data set for 370 samples, and a 1249-locus exon-capture data set for 106 samples indicate that the current taxonomy substantially understates the true number of Sulawesi Draco species, that both cryptic and arrested speciations have taken place, and that ancient hybridization confounds phylogenetic analyses that do not explicitly account for reticulation. The Draco lineatus Group appears to comprise 15 species-9 on Sulawesi proper and 6 on peripheral islands. The common ancestor of this group colonized Sulawesi ~11 Ma when proto-Sulawesi was likely composed of two ancestral islands, and began to radiate ~6 Ma as new islands formed and were colonized via overwater dispersal. The enlargement and amalgamation of many of these proto-islands into modern Sulawesi, especially during the past 3 Ma, set in motion dynamic species interactions as once-isolated lineages came into secondary contact, some of which resulted in lineage merger, and others surviving to the present. [Genomics; Indonesia; introgression; mitochondria; phylogenetics; phylogeography; population genetics; reptiles.].

Potential Contribution of Ancient Introgression to the Evolution of a Derived Reproductive Strategy in Ricefishes

Transitions from no parental care to extensive care are costly and involve major changes in life history, behavior, and morphology. Nevertheless, in Sulawesi ricefishes, pelvic brooding evolved from transfer brooding in two distantly related lineages within the genera Adrianichthys and Oryzias, respectively. Females of pelvic brooding species carry their eggs attached to their belly until the fry hatches. Despite their phylogenetic distance, both pelvic brooding lineages share a set of external morphological traits. A recent study found no direct gene flow between pelvic brooding lineages, suggesting independent evolution of the derived reproductive strategy. Convergent evolution can, however, also rely on repeated sorting of preexisting variation of an admixed ancestral population, especially when subjected to similar external selection pressures. We thus used a multispecies coalescent model and D-statistics to identify gene-tree-species-tree incongruencies, to evaluate the evolution of pelvic brooding with respect to interspecific gene flow not only between pelvic brooding lineages but also between pelvic brooding lineages and other Sulawesi ricefish lineages. We found a general network-like evolution in Sulawesi ricefishes, and as previously reported, we detected no gene flow between the pelvic brooding lineages. Instead, we found hybridization between the ancestor of pelvic brooding Oryzias and the common ancestor of the Oryzias species from the Lake Poso area. We further detected signs of introgression within the confidence interval of a quantitative trait locus associated with pelvic brooding in O. eversi. Our results hint toward a contribution of ancient standing genetic variation to the evolution of pelvic brooding in Oryzias.

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.

Evolutionary differentiation of androgen receptor is responsible for sexual characteristic development in a teleost fish

Teleost fishes exhibit complex sexual characteristics in response to androgens, such as fin enlargement and courtship display. However, the molecular mechanisms underlying their evolutionary acquisition remain largely unknown. To address this question, we analyse medaka (Oryzias latipes) mutants deficient in teleost-specific androgen receptor ohnologs (ara and arb). We discovered that neither ar ohnolog was required for spermatogenesis, whilst they appear to be functionally redundant for the courtship display in males. However, both were required for reproductive success: ara for tooth enlargement and the reproductive behaviour eliciting female receptivity, arb for male-specific fin morphogenesis and sexual motivation. We further showed that differences between the two ar ohnologs in their transcription, cellular localisation of their encoded proteins, and their downstream genetic programmes could be responsible for the phenotypic diversity between the ara and arb mutants. These findings suggest that the ar ohnologs have diverged in two ways: first, through the loss of their roles in spermatogenesis and second, through gene duplication followed by functional differentiation that has likely resolved the pleiotropic roles derived from their ancestral gene. Thus, our results provide insights into how genome duplication impacts the massive diversification of sexual characteristics in the teleost lineage.

Cryptic genetic divergence of the red dwarf rasbora, Microrasbora rubescens, in and around Inle Lake: implications for the origin of endemicity in the ancient lake in Myanmar

Inle Lake, an ancient lake located in the Shan Plateau of Myanmar, is a biogeographically attractive region with high fish endemism. Some endemic species inhabit the lake as well as the surrounding areas. The genetic and ecological relationships between populations in the lake and surrounding areas provide important insights into the process underlying ichthyofaunal formation in Inle Lake. In this study, the authors focused on red dwarf rasbora Microrasbora rubescens, an endemic genus and species in this region, and estimated its population structure and evolutionary scenario based on genome-wide polymorphism, mtDNA and geometric morphometric analyses using samples from Inle Lake and three areas surrounding the lake. The results showed that M. rubescens comprises at least three genetically divergent lineages (Inle, Heho and Hopong) with distinct geographic structures consistent with nuclear and mtDNA data. In contrast, there was no clear regional differentiation in morphology. The divergence time estimation based on mtDNA suggests that the Hopong lineage diverged at 2.7 Ma and the Inle and Heho lineages diverged at 1.9 Ma - consistent with the nuclear DNA results. The deep divergence observed in the endemic species supports the ancient history of ichthyofaunal development in this region. The distinct regional differentiation and morphological conservatism of this species might have been shaped by niche conservatism in stagnant water environments that limit dispersal and morphological diversification. Future comprehensive genetic and morphological analyses and comparisons for other native species should reveal the geographic and ecological processes that shaped the ichthyofauna in this region.

Population genetic characteristics of Hainan medaka with whole-genome resequencing

The DMY gene is deleted in all males of the Sanya population (SY-medaka) of the Hainan medaka, Oryzias curvinotus, as recently reported by us. However, due to limited knowledge regarding their population genetic background, it is difficult to explore the possible evolutionary pathway. Herein, we resequenced the whole genome of four populations, including SY-medaka. A total of 56 mitogenomes and 32,826,105 SNPs were identified. We found that the genetic differentiation is highest between SY-medaka and the other populations. The results of the population history of the O. curvinotus suggest that the SY-medaka has been in a bottleneck period recently. Further analysis shows that SY-medaka are the most strongly affected by environmental selection. Moreover, we screened some potential genomic regions, and the genes contained in these regions may explain the potential mechanism of the selection process of the SY-medaka. In conclusion, our study can provide new clues for the adaptation process of medaka in the new environment of Sanya.

Molecular phylogeny and phylogeography of ricefishes (Teleostei: Adrianichthyidae: Oryzias) in Sri Lanka

Ricefishes of the genus Oryzias occur commonly in the fresh and brackish waters in coastal lowlands ranging from India across Southeast Asia and on to Japan. Among the three species of Oryzias recorded from peninsular India, two widespread species, O. carnaticus and O. dancena, have previously been reported from Sri Lanka based on museum specimens derived from a few scattered localities. However, members of the genus are widespread in the coastal lowlands of Sri Lanka, a continental island separated from India by the shallow Palk Strait. Although recent molecular phylogenies of Adrianichthyidae represent near-complete taxon representation, they lack samples from Sri Lanka. Here, based on sampling at 13 locations representative of the entire geographic and climatic regions of the island's coastal lowlands, we investigate for the first time the molecular phylogenetic relationships and phylogeography of Sri Lankan Oryzias based on one nuclear and two mitochondrial markers. Sri Lankan Oryzias comprise two distinct non-sister lineages within the javanicus species group. One of these is represented by samples exclusively from the northern parts of the island; it is recognized as O. dancena. This lineage is recovered as the sister group to the remaining species in the javanicus group. The second lineage represents a species that is widespread across the island's coastal lowlands. It is recovered as the sister group of O. javanicus and is identified as O. cf. carnaticus. Ancestral-range estimates suggest two independent colonizations of Indian subcontinent and Sri Lanka by widespread ancestral species of Oryzias during two discrete temporal windows: late Miocene and Plio-Pleistocene. No phylogeographic structure is apparent in Sri Lankan Oryzias, suggesting that there are no strong barriers to gene flow and dispersal along the coastal floodplains, as is the case also for other generalist freshwater fishes in the island.

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

The genetic basis of a novel reproductive strategy in Sulawesi ricefishes: How modularity and a low number of loci shape pelvic brooding

The evolution of complex phenotypes like reproductive strategies is challenging to understand, as they often depend on multiple adaptations that only jointly result in a specific functionality. Sulawesi ricefishes (Adrianichthyidae) evolved a reproductive strategy termed as pelvic brooding. In contrast to the more common transfer brooding, female pelvic brooders carry an egg bundle connected to their body for weeks until the fry hatches. To examine the genetic architecture of pelvic brooding, we crossed the pelvic brooding Oryzias eversi and the transfer brooding Oryzias nigrimas (species divergence time: ∼3.6 my). We hypothesize, that a low number of loci and modularity have facilitated the rapid evolution of pelvic brooding. Traits associated to pelvic brooding, like rib length, pelvic fin length, and morphology of the genital papilla, were correlated in the parental species but correlations were reduced or lost in their F1 and F2 hybrids. Using the Castle-Wright estimator, we found that generally few loci underlie the studied traits. Further, both parental species showed modularity in their body plans. In conclusion, morphological traits related to pelvic brooding were based on a few loci and the mid-body region likely could evolve independently from the remaining body parts. Both factors presumably facilitated the evolution of pelvic brooding.

Inflammation and convergent placenta gene co-option contributed to a novel reproductive tissue

The evolution of pregnancy exposes parental tissues to new, potentially stressful conditions, which can trigger inflammation.¹ Inflammation is costly^2,3 and can induce embryo rejection, which constrains the evolution of pregnancy.¹ In contrast, inflammation can also promote morphological innovation at the maternal-embryonic interface as exemplified by co-option of pro-inflammatory signaling for eutherian embryo implantation.^1,4,5 Given its dual function, inflammation could be a key process explaining how innovations such as pregnancy and placentation evolved many times convergently. Pelvic brooding ricefishes evolved a novel “plug” tissue,^6,7 which forms inside the female gonoduct after spawning, anchors egg-attaching filaments, and enables pelvic brooders to carry eggs externally until hatching.^6,8 Compared to pregnancy, i.e., internal bearing of embryos, external bearing should alleviate constraints on inflammation in the reproductive tract. We thus hypothesized that an ancestral inflammation triggered by the retention of attaching filaments gave rise to pathways orchestrating plug formation. In line with our hypothesis, histological sections of the developing plug revealed signs of gonoduct injuries by egg-attaching filaments in the pelvic brooding ricefish Oryzias eversi. Tissue-specific transcriptomes showed that inflammatory signaling dominates the plug transcriptome and inflammation-induced genes controlling vital processes for plug development such as tissue growth and angiogenesis were overexpressed in the plug. Finally, mammalian placenta genes were enriched in the plug transcriptome, indicating convergent gene co-option for building, attaching, and sustaining a transient tissue in the female reproductive tract. This study highlights the role of gene co-option and suggests that recruiting inflammatory signaling into physiological processes provides a fast-track to evolutionary innovation.

•

Pelvic brooding induces tissue-specific changes in gene expression

•

Inflammatory signaling characterizes transcriptome of the egg-anchoring plug

•

Similar to embryo implantation, the plug likely evolved from an inflammatory response

•

Mammalian placenta genes were independently co-opted into the plug

Female preference for color-enhanced males: a test of the sensory bias model in medaka, a drab fish

Sexual selection research has long focused on the evolution of female mate preferences. Most of the models that have been developed posit that mate preferences evolve in a mating context. In contrast, the sensory bias model proposes that mate choice preferences arise in a non-mating context, as a by-product of natural selection acting on a female’s perceptual system. Recent research has shown that many species of fishes, from across a large clade including poeciliids, goodeids, and medaka, have a bias for long wavelength (LW) colors (yellow, orange, red) in a non-mating context. Even species that do not have LW-colored ornaments, apparently because they have been lost secondarily, retain this latent bias for LW colors. Here, we predicted that female Oryzias latipes (Japanese medaka), a drab species with a latent preference for LW colors, would show a mate choice preference for males with an artificial secondary sexual trait—a colored stripe added to their flank. We confirmed that females were more responsive to red and orange objects in a non-mating context than to other colors. We also showed that females were less resistant towards males with an LW-colored stripe than to those enhanced with a non-LW stripe and that, for many females, responses towards specific LW colors were consistent across these non-mating and mating contexts. Therefore, our results provide support for the sensory bias model by providing a link between a sensory bias in a non-mating context and a mate choice preference in a drab species like medaka.

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.

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.

Oryzias curvinotus in Sanya Does Not Contain the Male Sex-Determining Gene dmy

Hainan medaka (Oryzias curvinotus) is distributed in the coastal waters of the South China Sea and is able to adapt to a wide range of salinities. In this study, we characterized O. curvinotus in Sanya River (SY-medaka), which lacks dmy (a male sex-determining gene in O. latipes and O. curvinotus). In a comparison of SY-medaka and Gaoqiao medaka (GQ-medaka), the morphological difference between the two populations does not reach the subspecies level and they can be considered two geographic populations of O. curvinotus. A mitochondrial cytochrome oxidase subunit I (CoI) sequence alignment showed that the sequence identities between SY-medaka and other geographic populations of O. curvinotus are as high as 95%. A phylogenetic analysis of the mitochondrial genome also indicated that SY-medaka belongs to O. curvinotus. Molecular marker-based genetic sex assays and whole genome re-sequencing showed that SY-medaka does not contain dmy. Further, in RNA-Seq analyses of the testis and ovaries of sexually mature SY-medaka, dmy expression was not detected. We speculate that high temperatures resulted in the loss of dmy in SY-medaka during evolution, or the lineage has another sex-determining gene. This study provides a valuable dataset for elucidating the mechanism underlying sex determination in Oryzias genus and advances research on functional genomics or reproduction biology in O. curvinotus.

Complex sexually dimorphic traits shape the parallel evolution of a novel reproductive strategy in Sulawesi ricefishes (Adrianichthyidae)

Background

Pelvic brooding is a form of uni-parental care, and likely evolved in parallel in two lineages of Sulawesi ricefishes. Contrary to all other ricefishes, females of pelvic brooding species do not deposit eggs at a substrate (transfer brooding), but carry them until the fry hatches. We assume that modifications reducing the costs of egg carrying are beneficial for pelvic brooding females, but likely disadvantageous in conspecific males, which might be resolved by the evolution of sexual dimorphism via sexual antagonistic selection. Thus we hypothesize that the evolution of pelvic brooding gave rise to female-specific skeletal adaptations that are shared by both pelvic brooding lineages, but are absent in conspecific males and transfer brooding species. To tackle this, we combine 3D-imaging and morphometrics to analyze skeletal adaptations to pelvic brooding.

Results

The morphology of skeletal traits correlated with sex and brooding strategy across seven ricefish species. Pelvic brooding females have short ribs caudal of the pelvic girdle forming a ventral concavity and clearly elongated and thickened pelvic fins compared to both sexes of transfer brooding species. The ventral concavity limits the body cavity volume in female pelvic brooders. Thus body volumes are smaller compared to males in pelvic brooding species, a pattern sharply contrasted by transfer brooding species.

Conclusions

We showed in a comparative framework that highly similar, sexually dimorphic traits evolved in parallel in both lineages of pelvic brooding ricefish species. Key traits, present in all pelvic brooding females, were absent or much less pronounced in conspecific males and both sexes of transfer brooding species, indicating that they are non-beneficial or even maladaptive for ricefishes not providing extended care. We assume that the combination of ventral concavity and robust, elongated fins reduces drag of brooding females and provides protection and stability to the egg cluster. Thus ricefishes are one of the rare examples where environmental factors rather than sexual selection shaped the evolution of sexually dimorphic skeletal adaptations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01791-z.

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.

Latitudinal variation in sexual dimorphism in a freshwater fish group

Tropical animals are characterized by showy ornaments and conspicuous body colours as compared with their temperate relatives. Some recent studies have hypothesized that sexual selection pressures are stronger in the tropics than in the temperate zone. Although negative correlations between latitude and the degree of sexual dimorphism would support this hypothesis, phylogeny should be taken into account in such comparative studies. Comparisons of the degree of sexual dimorphism in body size and fin lengths among species of the Adrianichthyidae, a freshwater fish family having a wide geographical range throughout Southeast and East Asia, revealed that lower latitude species are sexually more dimorphic in all characters than higher latitude species. Phylogenetic generalized least squares analyses using a mitochondrial DNA phylogeny demonstrated that the negative correlations between latitude and the degree of sexual dimorphism become non-significant when phylogeny is considered, but that the variance in the degree of sexual dimorphism is explained not only by phylogeny but also almost equally by latitude. Ancestral state reconstruction indicated that sexual dimorphisms have evolved independently even within major clades. These findings are consistent with the view that tropical species are exposed to stronger sexual selection pressures than temperate species. We discuss possible causes of the latitudinal variation in sexual selection pressure.

Miocene geologic dynamics of the Australian Sahul Shelf determined the biogeographic patterns of freshwater planorbid snails (Miratestinae) in the Indo-Australian Archipelago

The complex geological and climatic processes that have shaped the Indo-Australian Archipelago since the Cenozoic likely also gave rise to its species-rich biota. Strictly freshwater organisms might be particularly suitable for understanding the influence of these abiotic factors on their biogeography in such a insular setting as their distribution may reflect past abiotic events at large and small geographical scales. We here investigate the historical biogeography of the Miratestinae, a subfamily of Planorbidae. These freshwater gastropods are widely distributed in the eastern IAA from Australia, New Guinea, the Moluccas, and Sulawesi to the Philippines. The first comprehensive molecular phylogeny of the Miratestinae was inferred based on two mitochondrial and two nuclear genetic markers using maximum likelihood and Bayesian inference. Four species delimitation methods were applied to identify molecular operational taxonomic units (MOTUs). Divergence times were inferred using an uncorrelated lognormal relaxed-clock model by applying a taxon- and marker-specific substitution rate. Ancestral geographic ranges were estimated based on the dated phylogeny using BioGeoBEARS. The species delimitation revealed a total of 23 MOTUs, 16 of which might represent species new to science. The BioGeoBEARS analyses suggest an Australian origin for the Miratestinae at c. 22 Ma and identified jump dispersal to be the main process of colonization. The first colonization events from Australia to the IAA occurred in the Middle-Late Miocene (12-13 Ma), whereas intra-island diversification took mainly place since the Late Miocene-Pliocene. Colonization and diversification events remarkably coincide with major geologic events that shaped the geography of the region. The increasing availability of landmasses along the Sahul Shelf likely promoted stepping-stone dispersal to New Guinea, Sulawesi and the Philippines as early as the islands emerged. Major geological and climatic events such as the amalgamation of the island Sulawesi, the regional aridification in Australia or the uplift of massive mountain ranges in New Guinea likely played a considerable role for intra-island diversification.

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.

The untapped potential of medaka and its wild relatives

The medaka is a fish that has served as a model organism for over a century, yet there is still much to learn about its life in the wild. Here we summarize the current knowledge, highlight recent progress and outline remaining gaps in our understanding of the natural history of medaka. It has also become clear over time that rather than being a single species, medaka comprises an entire species complex, so disentangling the species boundaries is an important goal for future research. Moreover, medaka and other ricefishes exhibit striking functional diversity, little of which has been investigated to date. As such, there are opportunities to use the resources developed for medaka to study other ricefishes, and to learn more about medaka itself in an evolutionary context.

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.

Hyperosmotic tolerance of adult fish and early embryos are determined by discrete, single loci in the genus Oryzias

The acquisition of environmental osmolality tolerance traits in individuals and gametes is an important event in the evolution and diversification of organisms. Although teleost fish exhibit considerable intra- and interspecific variation in salinity tolerance, the genetic mechanisms underlying this trait remain unclear. Oryzias celebensis survives in sea and fresh water during both the embryonic and adult stages, whereas its close relative Oryzias woworae cannot survive in sea water at either stage. A linkage analysis using backcross progeny identified a single locus responsible for adult hyperosmotic tolerance on a fused chromosome that corresponds to O. latipes linkage groups (LGs) 6 and 23. Conversely, O. woworae eggs fertilised with O. celebensis sperm died in sea water at the cleavage stages, whereas O. celebensis eggs fertilised with O. woworae sperm developed normally, demonstrating that maternal factor(s) from O. celebensis are responsible for hyperosmotic tolerance during early development. A further linkage analysis using backcrossed females revealed a discrete single locus relating to the maternal hyperosmotic tolerance factor in a fused chromosomal region homologous to O. latipes LGs 17 and 19. These results indicate that a maternal factor governs embryonic hyperosmotic tolerance and maps to a locus distinct from that associated with adult hyperosmotic tolerance.

Genetic Architecture of the Variation in Male-Specific Ossified Processes on the Anal Fins of Japanese Medaka

Traits involved in reproduction evolve rapidly and show great diversity among closely related species. However, the genetic mechanisms that underlie the diversification of courtship traits are mostly unknown. Japanese medaka fishes (Oryzias latipes) use anal fins to attract females and to grasp females during courtship; the males have longer anal fins with male-specific ossified papillary processes on the fin rays. However, anal fin morphology varies between populations: the southern populations tend to have longer anal fins and more processes than the northern populations. In the present study, we conducted quantitative trait locus (QTL) mapping to investigate the genetic architecture underlying the variation in the number of papillary processes of Japanese medaka fish and compared the QTL with previously identified QTL controlling anal fin length. First, we found that only a few QTL were shared between anal fin length and papillary process number. Second, we found that the numbers of papillary processes on different fin rays often were controlled by different QTL. Finally, we produced another independent cross and found that some QTL were repeatable between the two crosses, whereas others were specific to only one cross. These results suggest that variation in the number of papillary processes is polygenic and controlled by QTL that are distinct from those controlling anal fin length. Thus, different courtship traits in Japanese medaka share a small number of QTL and have the potential for independent evolution.

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.

Mitochondrial and nuclear phylogenetic trees and divergence time estimations of Sulawesi endemic Adrianichthyidae

This data article is related to the research article entitled “Origin and intra-island diversification of Sulawesi endemic Adrianichthyidae” by Mokodongan and Yamahira [1]. In this data article, we present phylogenetic trees of Sulawesi adrianichthyids separately reconstructed using mitochondrial (cytochrome b: cyt b and NADH dehydrogenase subunit 2: ND2) and nuclear (tyrosinase) sequences. We also present Bayesian chronograms of Sulawesi adrianichthyids separately estimated using a substitution rate for cyt b and for ND2.