Wataru Kai

A trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, Takifugu rubripes (fugu).

Heterogametic sex chromosomes have evolved independently in various lineages of vertebrates. Such sex chromosome pairs often contain nonrecombining regions, with one of the chromosomes harboring a master sex-determining (SD) gene. It is hypothesized that these sex chromosomes evolved from a pair of autosomes that diverged after acquiring the SD gene. By linkage and association mapping of the SD locus in fugu (Takifugu rubripes), we show that a SNP (C/G) in the anti-Müllerian hormone receptor type II (Amhr2) gene is the only polymorphism associated with phenotypic sex. This SNP changes an amino acid (His/Asp384) in the kinase domain. While females are homozygous (His/His384), males are heterozygous. Sex in fugu is most likely determined by a combination of the two alleles of Amhr2. Consistent with this model, the medaka hotei mutant carrying a substitution in the kinase domain of Amhr2 causes a female phenotype. The association of the Amhr2 SNP with phenotypic sex is conserved in two other species of Takifugu but not in Tetraodon. The fugu SD locus shows no sign of recombination suppression between X and Y chromosomes. Thus, fugu sex chromosomes represent an unusual example of proto–sex chromosomes. Such undifferentiated X-Y chromosomes may be more common in vertebrates than previously thought.

Integration of the Genetic Map and Genome Assembly of Fugu Facilitates Insights into Distinct Features of Genome Evolution in Teleosts and Mammals

Abstract The compact genome of fugu (Takifugu rubripes) has been used widely as a reference genome for understanding the evolution of vertebrate genomes. However, the fragmented nature of the fugu genome assembly has restricted its use for comparisons of genome architecture in vertebrates. To extend the contiguity of the assembly to the chromosomal level, we have generated a comprehensive genetic map of fugu and anchored the scaffolds of the assembly to the 22 chromosomes of fugu. The map consists of 1,220 microsatellite markers that provide anchor points to 697 scaffolds covering 86% of the genome assembly (http://www.fugu-sg.org/). The integrated genome map revealed a higher recombination rate in fugu compared with other vertebrates and a wide variation in the recombination rate between sexes and across chromosomes of fugu. We used the extended assembly to explore recent rearrangement events in the lineages of fugu, Tetraodon, and medaka and compared them with rearrangements in three mammalian (human, mouse, and opossum) lineages. Between the two pufferfishes, fugu has experienced fewer chromosomal rearrangements than Tetraodon. The gene order is more highly conserved in the three teleosts than in mammals largely due to a lower rate of interchromosomal rearrangements in the teleosts. These results provide new insights into the distinct patterns of genome evolution between teleosts and mammals. The consolidated genome map and the genetic map of fugu are valuable resources for comparative genomics of vertebrates and for elucidating the genetic basis of the phenotypic diversity of ∼25 species of Takifugu that evolved within the last 5 My.

The Sex-Determining Locus in the Tiger Pufferfish, Takifugu rubripes

The tiger pufferfish (fugu), Takifugu rubripes, is a model fish that has had its genome entirely sequenced. By performing genomewide linkage analyses, we show that the sex of fugu is determined by a single chromosomal region on linkage group 19 in an XX–XY system.

THE GENETIC ARCHITECTURE OF GROWTH RATE IN JUVENILE TAKIFUGU SPECIES

Closely related species have often evolved dramatic differences in body size. Takifugu rubripes (fugu) is a large marine pufferfish whose genome has been sequenced, whereas T. niphobles is the smallest species among Takifugu. We show that, unsurprisingly, the juvenile growth rate of T. rubripes is higher than that of T. niphobles in a laboratory setting. We produced F2 progenies of their F1 hybrids and found one quantitative trait locus (QTL) significantly associated with variation in juvenile body size. This QTL region (3.5 Mb) contains no known genes directly related to growth phenotype (such as IGFs) except Fgf21, which inhibits growth hormone signaling in mouse. The QTL in Takifugu spp. is distinct from the region previously known to control body size variations in stickleback or tilapia. Our results suggest that in the fish tested herein, genomic regions underlying body size evolution might have different genetic origins. They also suggest that many diverse traits in Takifugu spp. are amenable to genetic mapping.

Second generation physical and linkage maps of yellowtail (Seriola quinqueradiata) and comparison of synteny with four model fish

BackgroundPhysical and linkage maps are important aids for the assembly of genome sequences, comparative analyses of synteny, and to search for candidate genes by quantitative trait locus analysis. Yellowtail, Seriola quinqueradiata, is an economically important species in Japanese aquaculture, and genetic information will be useful for DNA-assisted breeding. We report the construction of a second generation radiation hybrid map, its synteny analysis, and a second generation linkage map containing SNPs (single nucleotide polymorphisms) in yellowtail.ResultsApproximately 1.4 million reads were obtained from transcriptome sequence analysis derived from 11 tissues of one individual. To identify SNPs, cDNA libraries were generated from a pool of 500 whole juveniles, and the gills and kidneys of 100 adults. 9,356 putative SNPs were detected in 6,025 contigs, with a minor allele frequency ≥25%. The linkage and radiation hybrid maps were constructed based on these contig sequences. 2,081 markers, including 601 SNPs markers, were mapped onto the linkage map, and 1,532 markers were mapped in the radiation hybrid map.ConclusionsThe second generation linkage and physical maps were constructed using 6,025 contigs having SNP markers. These maps will aid the de novo assembly of sequencing reads, linkage studies and the identification of candidate genes related to important traits. The comparison of marker contigs in the radiation hybrid map indicated that yellowtail is evolutionarily closer to medaka than to green-spotted pufferfish, three-spined stickleback or zebrafish. The synteny analysis may aid studies of chromosomal evolution in yellowtail compared with model fish.

Genomic regions of pufferfishes responsible for host specificity of a monogenean parasite, Heterobothrium okamotoi.

The genetic mechanisms underlying host specificity of parasitic infections are largely unknown. After hatching, the larvae of the monogenean parasite, Heterobothrium okamotoi, attach to the gill filaments of hosts and the post-larval worms develop there by consuming nutrients from the host. The susceptibility to H. okamotoi infection differs markedly among fish species. While this parasite can grow on tiger pufferfish (also called fugu), Takifugu rubripes, it appears to be rejected by a close congener, grass pufferfish, Takifugu niphobles, after initial attachment to the gills. To determine the genetic architecture of the pufferfish responsible for this host specificity, we performed genome-wide quantitative trait loci analysis. We raised second generation (F2) hybrids of the two pufferfish species and experimentally infected them with the monogenean in vivo. To assess possible differences in host mechanisms between early and later periods of infection, we sampled fish three h and 21days after exposure. Genome scanning of fish from the 3h infection trial revealed suggestive quantitative trait loci on linkage groups 2 and 14, which affected the number of parasites on the gill. However, analysis of fish 21days p.i. detected a significant quantitative trait locus on linkage group 9 and three other suggestive quantitative trait loci on linkage groups 7, 18 and 22. These results indicated the polygenic nature of the host mechanisms involved in the infection/rejection of H. okamotoi. Moreover the analyses suggested that host factors may play a more important role during the growth period of the parasite than during initial host recognition at the time of attachment. Within the 95% confidence interval of the linkage group 9 quantitative trait locus in the fugu genome, there were 214 annotated protein-coding genes, including immunity-related genes such as Irak4, Muc2 and Muc5ac.

Complete Genome Sequences of Edwardsiella tarda-Lytic Bacteriophages KF-1 and IW-1.

ABSTRACT We report the complete genome sequences of two Edwardsiella tarda-lytic bacteriophages isolated from flounder kidney (KF-1) and seawater (IW-1). These newly sequenced phage genomes provide a novel resource for future studies on phage-host interaction mechanisms and various applications of the phages for control of edwardsiellosis in aquaculture.

Complete Genome Sequence of a Novel Myovirus Which Infects Atypical Strains of Edwardsiella tarda

ABSTRACT We present the genome sequence of a novel Edwardsiella tarda-lytic bacteriophage, MSW-3, which specifically infects atypical E. tarda strains. The morphological and genomic features of MSW-3 suggest that this phage is a new member of the dwarf myoviruses, which have been much less studied than other groups of myoviruses.

Complete Genome Sequence of the Edwardsiella ictaluri-Specific Bacteriophage PEi21, Isolated from River Water in Japan.

ABSTRACT We present the complete genome sequence for a novel Edwardsiella ictaluri-specific bacteriophage, PEi21, isolated from river water in Japan. An initial comparative genome analysis revealed that the phage was closely related to the previously reported Edwardsiella tarda phage MSW-3 isolated from a red sea bream farm in Japan.

Effective de novo assembly of fish genome using haploid larvae

Recent improvements in next-generation sequencing technology have made it possible to do whole genome sequencing, on even non-model eukaryote species with no available reference genomes. However, de novo assembly of diploid genomes is still a big challenge because of allelic variation. The aim of this study was to determine the feasibility of utilizing the genome of haploid fish larvae for de novo assembly of whole-genome sequences. We compared the efficiency of assembly using the haploid genome of yellowtail (Seriola quinqueradiata) with that using the diploid genome obtained from the dam. De novo assembly from the haploid and the diploid sequence reads (100 million reads per each datasets) generated by the Ion Proton sequencer (200 bp) was done under two different assembly algorithms, namely overlap-layout-consensus (OLC) and de Bruijn graph (DBG). This revealed that the assembly of the haploid genome significantly reduced (approximately 22% for OLC, 9% for DBG) the total number of contigs (with longer average and N50 contig lengths) when compared to the diploid genome assembly. The haploid assembly also improved the quality of the scaffolds by reducing the number of regions with unassigned nucleotides (Ns) (total length of Ns; 45,331,916 bp for haploids and 67,724,360 bp for diploids) in OLC-based assemblies. It appears clear that the haploid genome assembly is better because the allelic variation in the diploid genome disrupts the extension of contigs during the assembly process. Our results indicate that utilizing the genome of haploid larvae leads to a significant improvement in the de novo assembly process, thus providing a novel strategy for the construction of reference genomes from non-model diploid organisms such as fish.