Search

The karyotype, which is the number and shape of chromosomes, is a fundamental characteristic of all eukaryotes. Karyotypic changes play an important role in many aspects of evolutionary processes, including speciation. In organisms with monocentric chromosomes, it was previously thought that chromosome number changes were mainly caused by centric fusions and fissions, whereas chromosome shape changes, that is, changes in arm numbers, were mainly due to pericentric inversions. However, recent genomic and cytogenetic studies have revealed examples of alternative cases, such as tandem fusions and centromere repositioning, found in the karyotypic changes within and between species. Here, we employed comparative genomic approaches to investigate whether centromere repositioning occurred during karyotype evolution in medaka fishes. In the medaka family (Adrianichthyidae), the three phylogenetic groups differed substantially in their karyotypes. The Oryzias latipes species group has larger numbers of chromosome arms than the other groups, with most chromosomes being metacentric. The O. javanicus species group has similar numbers of chromosomes to the O. latipes species group, but smaller arm numbers, with most chromosomes being acrocentric. The O. celebensis species group has fewer chromosomes than the other two groups and several large metacentric chromosomes that were likely formed by chromosomal fusions. By comparing the genome assemblies of O. latipes, O. javanicus, and O. celebensis, we found that repositioning of centromere‐associated repeats might be more common than simple pericentric inversion. Our results demonstrated that centromere repositioning may play a more important role in karyotype evolution than previously appreciated. [ABSTRACT FROM AUTHOR]

The genus Pristionchus (Kreis, 1932) consists of more than 30 soil nematode species that are often found in association with scarab beetles. Three major radiations have resulted in the “maupasi species group” in America, the “pacificus species group” in Asia, and the “lheritieri species group,” which contains species from Europe and Asia. Phylogenetic analysis indicates that a group of three species, including the gonochorists P. elegans and P. bucculentus and the hermaphrodite P. fissidentatus, is basal to the above-mentioned radiations. Two novel species are described here: Pristionchus paulseni sp. n. from Taiwan and P. yamagatae sp. n. from Japan by means of morphology, morphometrics and genome-wide transcriptome sequence analysis. Previous phylotranscriptomic analysis of the complete Pristionchus genus recognized P. paulseni sp. n. as the sister species of P. fissidentatus, and thus its importance for macro-evolutionary studies. Specifically, the gonochorist P. paulseni sp. n. and the hermaphrodite P. fissidentatus form a species pair that is the sister group to all other described Pristionchus species. P. paulseni sp. n. has two distinct mouth forms, supporting the notion that the mouth dimorphism is ancestral in the genus Pristionchus.

Karyotype, including the chromosome and arm numbers, is a fundamental genetic characteristic of all organisms and has long been used as a species-diagnostic character. Additionally, karyotype evolution plays an important role in divergent adaptation and speciation. Centric fusion and fission change chromosome numbers, whereas the intra-chromosomal movement of the centromere, such as pericentric inversion, changes arm numbers. A probabilistic model simultaneously incorporating both chromosome and arm numbers has not been established. Here, we built a probabilistic model of karyotype evolution based on the "karyograph", which treats karyotype evolution as a walk on the two-dimensional space representing the chromosome and arm numbers. This model enables analysis of the stationary distribution with a stable karyotype for any given parameter. After evaluating their performance using simulated data, we applied our model to two large taxonomic groups of fish, Eurypterygii and series Otophysi, to perform maximum likelihood estimation of the transition rates and reconstruct the evolutionary history of karyotypes. The two taxa significantly differed in the evolution of arm number. The inclusion of speciation and extinction rates demonstrated possibly high extinction rates in species with karyotypes other than the most typical karyotype in both groups. Finally, we made a model including polyploidization rates and applied it to a small plant group. Thus, the use of this probabilistic model can contribute to a better understanding of tempo and mode in karyotype evolution and its possible role in speciation and extinction. Author summary: Karyotype, including chromosome number and morphology, has been observed even before DNA was identified as the primary genetic material. Although chromosomal changes are thought to play an important role in speciation, we have not reached a consensus on how rapidly karyotype can evolve and whether a particular karyotype is favored. This can be attributed partly to the lack of good probabilistic models for karyotype evolution. This contrasts with DNA sequence evolution, for which many established probabilistic models are available. Such probabilistic models have contributed to the understanding of rates and driving forces of DNA sequence evolution. Here, we built a probabilistic model including both chromosome and arm numbers. Using this model, we could demonstrate differences in the tempo and mode in karyotype evolution between two fish taxonomic groups and possible roles of karyotype in speciation and extinction. The use of our model in diverse taxa will lead to a better understanding of the evolutionary trends and functional roles of karyotypes. [ABSTRACT FROM AUTHOR]

Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations. [ABSTRACT FROM AUTHOR]

Intragenomic conflict, the conflict of interest between different genomic regions within an individual, is proposed as a mechanism driving both the rapid evolution of heterochromatin‐related proteins and the establishment of intrinsic genomic incompatibility between species. Although molecular studies of laboratory model organisms have demonstrated the link between heterochromatin evolution and hybrid abnormalities, we know little about their link in natural systems. Previously, we showed that F1 hybrids between the Japan Sea stickleback and the Pacific Ocean stickleback show hybrid male sterility and found a region responsible for hybrid male sterility on the X chromosome, but did not identify any candidate genes. In this study, we first screened for genes rapidly evolving under positive selection during the speciation of Japanese sticklebacks to find genes possibly involved in intragenomic conflict. We found that the region responsible for hybrid male sterility contains a rapidly evolving gene encoding a heterochromatin‐binding protein TRIM24B. We conducted biochemical experiments and showed that the binding affinity of TRIM24B to a heterochromatin mark found at centromeres and transposons, histone H4 lysine 20 trimethylation (H4K20me3), is reduced in the Japan Sea stickleback. In addition, mRNA expression levels of Trim24b were different between the Japan Sea and the Pacific Ocean testes. Further expression analysis of genes possibly in the TRIM24B‐regulated pathway showed that some gypsy retrotransposons are overexpressed in the F1 hybrid testes. We, therefore, demonstrate that a heterochromatin‐binding protein can evolve rapidly under positive selection and functionally diverge during stickleback speciation. [ABSTRACT FROM AUTHOR]

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68–1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow. [ABSTRACT FROM AUTHOR]

Adaptation to different salinities can drive and maintain divergence between populations of aquatic organisms. Anadromous and stream ecotypes of threespine stickleback ( Gasterosteus aculeatus) are an excellent model to explore the genetic mechanisms underlying osmoregulation divergence. Using a parapatric pair of anadromous and stream stickleback ecotypes, we employed an integrated genomic approach to identify candidate genes important for adaptation to different salinity environments. Quantitative trait loci ( QTL) mapping of plasma sodium concentrations under a seawater challenge experiment identified a significant QTL on chromosome 16. To identify candidate genes within this QTL, we first conducted RNA-seq and microarray analysis on gill tissue to find ecotypic differences in gene expression that were associated with plasma Na+ levels. This resulted in the identification of ten candidate genes. Quantitative PCR analysis on gill tissue of additional Japanese stickleback populations revealed that the majority of the candidate genes showed parallel divergence in expression levels. Second, we conducted whole-genome sequencing and found five genes that are predicted to have functionally important amino acid substitutions. Finally, we conducted genome scan analysis and found that eight of these candidate genes were located in genomic islands of high differentiation, suggesting that they may be under divergent selection. The candidate genes included those involved in ATP synthesis and hormonal signalling, whose expression or amino acid changes may underlie the variation in salinity tolerance. Further functional molecular analysis of these genes will reveal the causative genetic and genomic changes underlying divergent adaptation. [ABSTRACT FROM AUTHOR]

Invasive species pose a major threat to biological diversity. Although introduced populations often experience population bottlenecks, some invasive species are thought to be originated from hybridization between multiple populations or species, which can contribute to the maintenance of high genetic diversity. Recent advances in genome sequencing enable us to trace the evolutionary history of invasive species even at whole-genome level and may help to identify the history of past hybridization that may be overlooked by traditional marker-based analysis. Here, we conducted whole-genome sequencing of eight threespine stickleback ( Gasterosteus aculeatus) individuals, four from a recently introduced crater lake population and four of the putative source population. We found that both populations have several small genomic regions with high genetic diversity, which resulted from introgression from a closely related species ( Gasterosteus nipponicus). The sizes of the regions were too small to be detected with traditional marker-based analysis or even some reduced-representation sequencing methods. Further amplicon sequencing revealed linkage disequilibrium around an introgression site, which suggests the possibility of selective sweep at the introgression site. Thus, interspecies introgression might predate introduction and increase genetic variation in the source population. Whole-genome sequencing of even a small number of individuals can therefore provide higher resolution inference of history of introduced populations. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Sex chromosomes turn over rapidly in some taxonomic groups, where closely related species have different sex chromosomes. Although there are many examples of sex chromosome turnover, we know little about the functional roles of sex chromosome turnover in phenotypic diversification and genomic evolution. The sympatric pair of Japanese threespine stickleback (Gasterosteus aculeatus) provides an excellent system to address these questions: the Japan Sea species has a neo-sex chromosome system resulting from a fusion between an ancestral Y chromosome and an autosome, while the sympatric Pacific Ocean species has a simple XY sex chromosome system. Furthermore, previous quantitative trait locus (QTL) mapping demonstrated that the Japan Sea neo-X chromosome contributes to phenotypic divergence and reproductive isolation between these sympatric species. To investigate the genomic basis for the accumulation of genes important for speciation on the neo-X chromosome, we conducted whole genome sequencing of males and females of both the Japan Sea and the Pacific Ocean species. No substantial degeneration has yet occurred on the neo-Y chromosome, but the nucleotide sequence of the neo-X and the neo-Y has started to diverge, particularly at regions near the fusion. The neo-sex chromosomes also harbor an excess of genes with sex-biased expression. Furthermore, genes on the neo-X chromosome showed higher non-synonymous substitution rates than autosomal genes in the Japan Sea lineage. Genomic regions of higher sequence divergence between species, genes with divergent expression between species, and QTL for inter-species phenotypic differences were found not only at the regions near the fusion site, but also at other regions along the neo-X chromosome. Neo-sex chromosomes can therefore accumulate substitutions causing species differences even in the absence of substantial neo-Y degeneration. [ABSTRACT FROM AUTHOR]

Background: Non-coding small RNAs, ranging from 20 to 30 nucleotides in length, mediate the regulation of gene expression and play important roles in many biological processes. One class of small RNAs, microRNAs (miRNAs), are highly conserved across taxa and mediate the regulation of the chromatin state and the post-transcriptional regulation of messenger RNA (mRNA). Another class of small RNAs is the Piwi-interacting RNAs, which play important roles in the silencing of transposons and other functional genes. Although the biological functions of the different small RNAs have been elucidated in several laboratory animals, little is known regarding naturally occurring variation in small RNA transcriptomes among closely related species. Results: We employed next-generation sequencing technology to compare the expression profiles of brain small RNAs between sympatric species of the Japanese threespine stickleback (Gasterosteus aculeatus). We identified several small RNAs that were differentially expressed between sympatric Pacific Ocean and Japan Sea sticklebacks. Potential targets of several small RNAs were identified as repetitive sequences. Female-biased miRNA expression from the old X chromosome was also observed, and it was attributed to the degeneration of the Y chromosome. Conclusions: Our results suggest that expression patterns of small RNA can differ between incipient species and may be a potential mechanism underlying differential mRNA expression and transposon activity. [ABSTRACT FROM AUTHOR]

The endemic cichlid fishes in Lake Victoria are a model system for speciation through adaptive radiation. Although the evolution of the sex-determination system may also play a role in speciation, little is known about the sex-determination system of Lake Victoria cichlids. To understand the evolution of the sex-determination system in these fish, we performed cytogenetic analysis in 11 cichlid species from Lake Victoria. B chromosomes, which are present in addition to standard chromosomes, were found at a high prevalence rate (85%) in these cichlids. In one species, B chromosomes were female-specific. Cross-breeding using females with and without the B chromosomes demonstrated that the presence of the B chromosomes leads to a female-biased sex ratio in this species. Although B chromosomes were believed to be selfish genetic elements with little effect on phenotype and to lack protein-coding genes, the present study provides evidence that B chromosomes have a functional effect on female sex determination. FISH analysis using a BAC clone containing B chromosome DNA suggested that the B chromosomes are derived from sex chromosomes. Determination of the nucleotide sequences of this clone (104.5 kb) revealed the presence of several protein-coding genes in the B chromosome, suggesting that B chromosomes have the potential to contain functional genes. Because some sex chromosomes in amphibians and arthropods are thought to be derived from B chromosomes, the B chromosomes in Lake Victoria cichlids may represent an evolutionary transition toward the generation of sex chromosomes. [ABSTRACT FROM AUTHOR]

Identifying the proximate and ultimate mechanisms of social behavior remains a major goal of behavioral biology. In particular, the complex social interactions mediating schooling behavior have long fascinated biologists, leading to theoretical and empirical investigations that have focused on schooling as a group-level phenomenon. However, methods to examine the behavior of individual fish within a school are needed in order to investigate the mechanisms that underlie both the performance and the evolution of schooling behavior. We have developed a technique to quantify the schooling behavior of an individual in standardized but easily manipulated social circumstances. Using our model school assay, we show that threespine sticklebacks (Gasterosteus aculeatus) from alternative habitats differ in behavior when tested in identical social circumstances. Not only do marine sticklebacks show increased association with the model school relative to freshwater benthic sticklebacks, they also display a greater degree of parallel swimming with the models. Taken together, these data indicate that marine sticklebacks exhibit a stronger tendency to school than benthic sticklebacks. We demonstrate that these population-level differences in schooling tendency are heritable and are shared by individuals within a population even when they have experienced mixed-population housing conditions. Finally, we begin to explore the stimuli that elicit schooling behavior in these populations. Our data suggest that the difference in schooling tendency between marine and benthic sticklebacks is accompanied by differential preferences for social vs. non-social and moving vs. stationary shelter options. Our study thus provides novel insights into the evolution of schooling behavior, as well as a new experimental approach to investigate the genetic and neural mechanisms that underlie this complex social behavior. [ABSTRACT FROM AUTHOR]

Nematodes such as Caenorhabditis elegans and Pristionchus pacificus are extremely successful model organisms for comparative biology. Several studies have shown that phenotypic novelty but also conserved processes are controlled by taxon-restricted genes. To trace back the evolution of such new or rapidly evolving genes, a robust phylogenomic framework is indispensable. Here, we present an improved version of the genome of Parapristionchus giblindavisi which is the only known member of the sister group of Pristionchus. Relative to the previous short-read assembly, the new genome is based on long reads and displays higher levels of contiguity, completeness, and correctness. Specifically, the number of contigs dropped from over 7,303 to 735 resulting in an N50 increase from 112 to 791 kb. We made use of the new genome to revisit the evolution of multiple gene families. This revealed Pristionchus-specific expansions of several environmentally responsive gene families and a Pristionchus-specific loss of the de novo purine biosynthesis pathway. Focusing on the evolution of sulfatases and sulfotransferases, which control the mouth form plasticity in P. pacificus, reveals differences in copy number and genomic configurations between the genera Pristionchus and Parapristionchus. Altogether, this demonstrates the utility of the P. giblindavisi genome to date and polarizes lineage-specific patterns. [ABSTRACT FROM AUTHOR]

Cichlid fishes in East Africa have undergone extensive adaptive radiation, which has led to spectacular diversity in their morphology and ecology. To date, genetic linkage maps have been constructed for several tilapias (riverine), Astatotilapia burtoni (Lake Tanganyika), and hybrid lines of Lake Malawi cichlids to facilitate genome-wide comparative analyses. In the present study, we constructed a genetic linkage map of the hybrid line of Lake Victoria cichlids, so that maps of cichlids from all the major areas of East Africa will be available. The genetic linkage map shown here is derived from the F 2 progeny of an interspecific cross between Haplochromis chilotes and Haplochromis sauvagei and is based on 184 microsatellite and two single-nucleotide polymorphism (SNP) markers. Most of the microsatellite markers used in the present study were originally designed for other genetic linkage maps, allowing us to directly compare each linkage group (LG) among different cichlid groups. We found 25 LGs, the total length of which was 1133.2 cM with an average marker spacing of about 6.09 cM. Our subsequent linkage mapping analysis identified two putative sex-determining loci in cichlids. Interestingly, one of these two loci is located on cichlid LG5, on which the female heterogametic ZW locus and several quantitative trait loci (QTLs) related to adaptive evolution have been reported in Lake Malawi cichlids. We also found that V1R1 and V1R2 , candidate genes for the fish pheromone receptor, are located very close to the recently detected sex-determining locus on cichlid LG5. The genetic linkage map study presented here may provide a valuable foundation for studying the chromosomal evolution of East African cichlids and the possible role of sex chromosomes in generating their genomic diversity. [ABSTRACT FROM AUTHOR]

Intrinsic postzygotic isolation typically appears as reduced viability or fertility of interspecific hybrids caused by genetic incompatibilities between diverged parental genomes. Dobzhansky-Muller interactions among individual genes, and chromosomal rearrangements causing problems with chromosome synapsis and recombination in meiosis, have both long been considered as major mechanisms behind intrinsic postzygotic isolation. Recent research has, however, suggested that the genetic basis of intrinsic postzygotic isolation can be more complex and involves, for example, overall divergence of the DNA sequence or epigenetic changes. Here, we review the mechanisms of intrinsic postzygotic isolation from genic, chromosomal, genomic, and epigenetic perspectives across diverse taxa. We provide empirical evidence for these mechanisms, discuss their importance in the speciation process, and highlight questions that remain unanswered., (Copyright © 2023 Cold Spring Harbor Laboratory Press; all rights reserved.)

Pristionchus pacificus Sommer, Carta, Kim, and Sternberg, 1996 is an important model organism in evolutionary biology that aims to integrate developmental biology and evo-devo with population genetics and ecology. Functional studies in P. pacificus are supported by a well-established phylogenetic framework of around 30 species of the genus Pristionchus that have been described in the last decade based on their entomophilic and necromenic association with scarab beetles. Biogeographically, East Asia has emerged as a hotspot of Pristionchus speciation and recent samplings have therefore focused on Islands and mainland settings in East Asia. Here, we describe in a series of three publications the results of our sampling efforts in Taiwan, Japan, and Hongkong in 2016 and 2017. We describe a total of nine new species that cover different phylogenetic species-complexes of the Pristionchus genus. In this first publication, we describe two new species, Pristionchus sikae sp. n. and Pristionchus kurosawai sp. n. that are closely related to P. pacificus . Together with five previously described species they form the " pacificus species-complex sensu stricto" that is characterized by all species forming viable, but sterile F1 hybrids indicating reproductive isolation. P. sikae sp. n. and P. kurosawai sp. n. have a gonochorist mode of reproduction and they are described using morphology, morphometrics, mating experiments, and genome-wide sequence analysis. We discuss the extreme diversification in the pacificus species-complex sensu stricto in East Asia and its potential power to study speciation processes., (© The Society of Nematologists 2018.)

The authors describe five new species of Pristionchus from Japan and Hongkong. Scarab beetle samplings in Hongkong identified P. hongkongensis sp. n. and P. neolucani sp. n., representing the first beetle-associated Pristionchus species from China. Surprisingly, samplings of millipedes in Japan revealed a previously unknown association of Pristionchus nematodes with these arthropods. Specifically, the authors found three previously known Pristionchus species, P. arcanus , P. entomophagus , and P. fukushimae on Japanese millipedes. In addition, the authors found three new Pristionchus species on millipedes, which are described as P. riukiariae sp. n., P. degawai sp. n., and P. laevicollis , sp. n., the latter of which was also found on stag beetles. These species are most closely related to P. maxplancki , P. japonicus , and P. quartusdecimus and belong to the pacificus species-complex. The authors describe all species based on morphology, morphometrics, and genome-wide sequence analysis. Mating experiments indicated that all species are reproductively isolated from each other and in contrast to the species of the " pacificus species-complex sensu stricto " they do not form F1 hybrids., The authors describe five new species of Pristionchus from Japan and Hongkong. Scarab beetle samplings in Hongkong identified P. hongkongensis sp. n. and P. neolucani sp. n., representing the first beetle-associated Pristionchus species from China. Surprisingly, samplings of millipedes in Japan revealed a previously unknown association of Pristionchus nematodes with these arthropods. Specifically, the authors found three previously known Pristionchus species, P. arcanus , P. entomophagus , and P. fukushimae on Japanese millipedes. In addition, the authors found three new Pristionchus species on millipedes, which are described as P. riukiariae sp. n., P. degawai sp. n., and P. laevicollis , sp. n., the latter of which was also found on stag beetles. These species are most closely related to P. maxplancki , P. japonicus , and P. quartusdecimus and belong to the pacificus species-complex. The authors describe all species based on morphology, morphometrics, and genome-wide sequence analysis. Mating experiments indicated that all species are reproductively isolated from each other and in contrast to the species of the “ pacificus species-complex sensu stricto ” they do not form F1 hybrids.

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., (Copyright © 2015 Kawajiri et al.)