Repositioning of centromere-associated repeats during karyotype evolution in Oryzias fishes.

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.

Mechanisms of Intrinsic Postzygotic Isolation: From Traditional Genic and Chromosomal Views to Genomic and Epigenetic Perspectives.

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.

Do sex-linked male meiotic drivers contribute to intrinsic hybrid incompatibilities? Recent empirical studies from flies and rodents.

Intrinsic hybrid incompatibility is one of the important isolating barriers between species. In organisms with sex chromosomes, intrinsic hybrid incompatibility often follows two rules: Haldane's rule and large-X effects. One explanation for these two rules is that sex chromosomes are hotspots for meiotic drivers that can cause intrinsic hybrid incompatibility between geographically isolated populations. Although this hypothesis seems plausible and several empirical data are consistent with it, we are still unsure whether such mechanisms occur in nature, particularly with respect to speciation with gene flow. Here, we review empirical studies that have investigated the roles of meiotic drive in sex-chromosome evolution and speciation and propose future studies necessary for testing this hypothesis.

Chromosome fusions repatterned recombination rate and facilitated reproductive isolation during Pristionchus nematode speciation.

Large-scale genome-structural evolution is common in various organisms. Recent developments in speciation genomics revealed the importance of inversions, whereas the role of other genome-structural rearrangements, including chromosome fusions, have not been well characterized. We study genomic divergence and reproductive isolation of closely related nematodes: the androdioecious (hermaphroditic) model Pristionchus pacificus and its dioecious sister species Pristionchus exspectatus. A chromosome-level genome assembly of P. exspectatus using single-molecule and Hi-C sequencing revealed a chromosome-wide rearrangement relative to P. pacificus. Strikingly, genomic characterization and cytogenetic studies including outgroup species Pristionchus occultus indicated two independent fusions involving the same chromosome, ChrIR, between these related species. Genetic linkage analysis indicated that these fusions altered the chromosome-wide pattern of recombination, resulting in large low-recombination regions that probably facilitated the coevolution between some of the ~14.8% of genes across the entire genomes. Quantitative trait locus analyses for hybrid sterility in all three sexes revealed that major quantitative trait loci mapped to the fused chromosome ChrIR. While abnormal chromosome segregations of the fused chromosome partially explain hybrid female sterility, hybrid-specific recombination that breaks linkage of genes in the low-recombination region was associated with hybrid male sterility. Thus, recent chromosome fusions repatterned recombination rate and drove reproductive isolation during Pristionchus speciation.

The improved genome of the nematode Parapristionchus giblindavisi provides insights into lineage-specific gene family evolution.

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.

Tempo and mode in karyotype evolution revealed by a probabilistic model incorporating both chromosome number and morphology.

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.

Accumulation of Deleterious Mutations in Landlocked Threespine Stickleback Populations.

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.

Two new Species of Pristionchus (Nematoda: Diplogastridae) include the Gonochoristic Sister Species of P. fissidentatus.

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

Functional divergence of a heterochromatin-binding protein during stickleback speciation.

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.

Two New Species of Pristionchus (Nematoda: Diplogastridae) from Taiwan and the Definition of the pacificus Species-Complex Sensu Stricto.

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.

Phylotranscriptomics of Pristionchus Nematodes Reveals Parallel Gene Loss in Six Hermaphroditic Lineages.

Mutation and recombination are main drivers of phenotypic diversity, but the ability to create new allelic combinations is strongly dependent on the mode of reproduction. While most animals are dioecious (i.e., separated male and female sexes), in a number of evolutionary lineages females have gained the ability to self-fertilize [1, 2], with drastic consequences on effective recombination rate, genetic diversity, and the efficacy of selection [3]. In the genus Caenorhabditis, such hermaphroditic or androdioecious lineages, including C. briggsae and C. tropicalis, display a genome shrinkage relative to their dioecious sister species C. nigoni and C. brenneri, respectively [4, 5]. However, common consequences of reproductive modes on nematode genomes remain unknown, because most taxa contain single or few androdioecious species. One exception is the genus Pristionchus, with seven androdioecious species. Pristionchus worms are found in association with scarab beetles in worldwide samplings, resulting in deep taxon sampling and currently 39 culturable and available species. Here, we use phylotranscriptomics of all 39 Pristionchus species to provide a robust phylogeny based on an alignment of more than 2,000 orthologous clusters, which indicates that the seven androdioecious species represent six independent lineages. We show that gene loss is more prevalent in all hermaphroditic lineages than in dioecious relatives and that the majority of lost genes evolved recently in the Pristionchus genus. Further, we provide evidence that genes with male-biased expression are preferentially lost in hermaphroditic lineages. This supports a contribution of adaptive gene loss to shaping nematode genomes following the evolution of hermaphroditism.

The genomic landscape at a late stage of stickleback speciation: High genomic divergence interspersed by small localized regions of introgression.

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.

Samplings of Millipedes in Japan and Scarab Beetles in Hong Kong result in five new Species of Pristionchus (Nematoda: Diplogastridae).

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.

Contribution of gene flow to the evolution of recombination suppression in sex chromosomes.

Polymorphism of alleles that benefit one sex but harm the other (sexually antagonistic alleles) generates selective pressures for reduced recombination between themselves and sex-determination loci. Such polymorphism can be maintained within a population when selection coefficients are sufficiently balanced between males and females. However, if regulatory mutations restrict gene expression only to one sex, these alleles become neutral in the other sex and easily fixed within a population, removing the selective pressures for recombination suppression in sex chromosomes. When there is spatial variation in selection regimes, however, alleles that are deleterious in one sex and neutral in the other can be maintained in other neighboring populations and gene flow may continuously supply deleterious alleles. We hypothesized that this maintenance of genetic variation may promote the establishment of recombination suppression in sex chromosomes even in cases where selection is limited to one sex. Using individual-based simulations, we show that spatial variation in male-limited selection and gene flow can promote the establishment of Y-autosome fusions, a special case of recombination suppression in sex chromosomes. This can be explained by the fact that fused Y-chromosomes that capture alleles that are beneficial for local males have a higher mean fitness compared to unfused Y chromosomes in the presence of deleterious gene flow. We also simulated the case of sex-concordant selection and found that gene flow of alleles that are deleterious in both sexes did not substantially increase the establishment rates of Y-autosome fusions across the parameter space examined. This can be accounted for by the fact that foreign alleles that are deleterious in both sexes can be efficiently removed from the population compared to alleles that are neutral in females. These results indicate that how gene flow affects the establishment rates of Y-autosome fusions depends largely on selection regimes. Spatial variation in sex-specific selection and gene flow should be appreciated as a factor affecting sex chromosome evolution.

Different contributions of local- and distant-regulatory changes to transcriptome divergence between stickleback ecotypes.

Differential gene expression can play an important role in phenotypic evolution and divergent adaptation. Although differential gene expression can be caused by both local- and distant-regulatory changes, we know little about their relative contribution to transcriptome evolution in natural populations. Here, we conducted expression quantitative trait loci (eQTL) analysis to investigate the genetic architecture underlying transcriptome divergence between marine and stream ecotypes of threespine sticklebacks (Gasterosteus aculeatus). We identified both local and distant eQTLs, some of which constitute hotspots, regions with a disproportionate number of significant eQTLs relative to the genomic background. The majority of local eQTLs including those in the hotspots caused expression changes consistent with the direction of transcriptomic divergence between ecotypes. Genome scan analysis showed that many local eQTLs overlapped with genomic regions of high differentiation. In contrast, nearly half of the distant eQTLs including those in the hotspots caused opposite expression changes, and few overlapped with regions of high differentiation, indicating that distant eQTLs may act as a constraint of transcriptome evolution. Finally, a comparison between two salinity conditions revealed that nearly half of eQTL hotspots were environment specific, suggesting that analysis of genetic architecture in multiple conditions is essential for predicting response to selection.

Genetic basis for variation in salinity tolerance between stickleback ecotypes.

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.

Accumulation of Deleterious Mutations on the Neo-Y Chromosome of Japan Sea Stickleback (Gasterosteus nipponicus).

Degeneration of Y chromosomes is a common evolutionary path of XY sex chromosome systems. Recent genomic studies in flies and plants have revealed that even young neo-sex chromosomes with the age of a few million years show signs of Y degeneration, such as the accumulation of nonsense and frameshift mutations. However, it remains unclear whether neo-Y chromosomes also show rapid degeneration in fishes, which often have homomorphic sex chromosomes. Here, we investigated whether a neo-Y chromosome of Japan Sea stickleback (Gasterosteus nipponicus), which was formed by a Y-autosome fusion within the last 2 million years, accumulates deleterious mutations. Our previous genomic analyses did not detect excess nonsense and frameshift mutations on the Japan Sea stickleback neo-Y. In the present study, we found that the nonrecombining region of the neo-Y near the fusion end has accumulated nonsynonymous mutations altering amino acids of evolutionarily highly conserved residues. Enrichment of gene ontology terms related to protein phosphorylation and cellular protein modification process was found in the genes with potentially deleterious mutations on the neo-Y. These results suggest that the neo-Y of the Japan Sea stickleback has already accumulated mutations that may impair protein functions.

Whole-genome sequencing reveals small genomic regions of introgression in an introduced crater lake population of threespine stickleback.

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.

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.

A microsatellite-based genetic linkage map and putative sex-determining genomic regions in Lake Victoria cichlids.

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 F2 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.2cM with an average marker spacing of about 6.09cM. 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.