Y-specific amh allele, amhy, is the master sex-determining gene in Japanese flounder Paralichthys olivaceus.

Japanese flounder (Paralichthys olivaceus) is an important marine fish species of both fisheries and aquaculture in Northeast Asia. The commercial interest for all-female progenies due to several sex-related traits has prompted basic research on the mechanisms of sex determination in this species. By conducting a linkage analysis of the sex-determining locus, we initially identified 12 microsatellite markers linked to sex in 11 scaffolds, whose localization was restricted to a specific region of linkage group 9. Sequence analysis of this region identified 181 genes based on the UniProt database annotations. Among them, the amh gene was considered a potential candidate for sex determination because this gene is known to have taken over the role of sex determination in many teleosts. An in-depth sequence analysis of both the coding and non-coding regions of amh in XX and XY individuals detected nine SNPs linked with maleness. However, because these substitutions were synonymous, the upstream and downstream regions of amh were also investigated and a male-specific variant with deletions in the promoter region was detected. This truncated Y-specific amh variant was named amhy, and the amh shared by both sexes was named amhx. The association analysis using both females and males of the genotypic sex inferred by the presence/absence of amhy found complete association with phenotypic sex and genotype. Gene expression analysis in larvae derived from a single-pair progeny by quantitative real-time PCR detected amhy transcripts in the larval trunks between 20 and 100 days after hatching only in XY larvae. Localization of amhy by in situ hybridization was detected in presumptive Sertoli cells of XY gonads. Expression of amhx was almost undetectable in both XX and XY genotypes. Loss of Amh function by CRISPR-Cas9 induced male-to-female sex reversal, indicating that this gene was necessary for the masculinization of XY individuals. In conclusion, the complete linkage of amhy with males, its early expression in XY gonads before testicular differentiation, and the induction of sex reversal by loss-of-function mutation support the view that amhy is the sex-determining gene in this species.

Identification of Sex-Linked SNPs and Sex-Determining Regions in the Yellowtail Genome.

Unlike the conservation of sex-determining (SD) modes seen in most mammals and birds, teleost fishes exhibit a wide variety of SD systems and genes. Hence, the study of SD genes and sex chromosome turnover in fish is one of the most interesting topics in evolutionary biology. To increase resolution of the SD gene evolutionary trajectory in fish, identification of the SD gene in more fish species is necessary. In this study, we focused on the yellowtail, a species widely cultivated in Japan. It is a member of family Carangidae in which no heteromorphic sex chromosome has been observed, and no SD gene has been identified to date. By performing linkage analysis and BAC walking, we identified a genomic region and SNPs with complete linkage to yellowtail sex. Comparative genome analysis revealed the yellowtail SD region ancestral chromosome structure as medaka-fugu. Two inversions occurred in the yellowtail linage after it diverged from the yellowtail-medaka ancestor. An association study using wild yellowtails and the SNPs developed from BAC ends identified two SNPs that can reasonably distinguish the sexes. Therefore, these will be useful genetic markers for yellowtail breeding. Based on a comparative study, it was suggested that a PDZ domain containing the GIPC protein might be involved in yellowtail sex determination. The homomorphic sex chromosomes widely observed in the Carangidae suggest that this family could be a suitable marine fish model to investigate the early stages of sex chromosome evolution, for which our results provide a good starting point.

Construction of a high-coverage bacterial artificial chromosome library and comprehensive genetic linkage map of yellowtail Seriola quinqueradiata.

BACKGROUND: Japanese amberjack/yellowtail (Seriola quinqueradiata) is a commonly cultured marine fish in Japan. For cost effective fish production, a breeding program that increases commercially important traits is one of the major solutions. In selective breeding, information of genetic markers is useful and sufficient to identify individuals carrying advantageous traits but if the aim is to determine the genetic basis of the trait, large insert genomic DNA libraries are essential. In this study, toward prospective understanding of genetic basis of several economically important traits, we constructed a high-coverage bacterial artificial chromosome (BAC) library, obtained sequences from the BAC-end, and constructed comprehensive female and male linkage maps of yellowtail using Simple Sequence Repeat (SSR) markers developed from the BAC-end sequences and a yellowtail genomic library. RESULTS: The total insert length of the BAC library we constructed here was estimated to be approximately 11 Gb and hence 16-times larger than the yellowtail genome. Sequencing of the BAC-ends showed a low fraction of repetitive sequences comparable to that in Tetraodon and fugu. A total of 837 SSR markers developed here were distributed among 24 linkage groups spanning 1,026.70 and 1,057.83 cM with an average interval of 4.96 and 4.32 cM in female and male map respectively without any segregation distortion. Oxford grids suggested conserved synteny between yellowtail and stickleback. CONCLUSIONS: In addition to characteristics of yellowtail genome such as low repetitive sequences and conserved synteny with stickleback, our genomic and genetic resources constructed and revealed here will be powerful tools for the yellowtail breeding program and also for studies regarding the genetic basis of traits.

Construction of a radiation hybrid panel and the first yellowtail (Seriola quinqueradiata) radiation hybrid map using a nanofluidic dynamic array.

BACKGROUND: Yellowtail (Seriola quinqueradiata) are an economically important species in Japan. However, there are currently no methods for captive breeding and early rearing for yellowtail. Thus, the commercial cultivation of this species is reliant upon the capture of wild immature fish. Given this, there is a need to develop captive breeding techniques to reduce pressure on wild stocks and facilitate the sustainable development of yellowtail aquaculture. We constructed a whole genome radiation hybrid (RH) panel for yellowtail gene mapping and developed a framework physical map using a nanofluidic dynamic array to use SNPs (single nucleotide polymorphisms) in ESTs (expressed sequence tags) for the DNA-assisted breeding of yellowtail. RESULTS: Clonal RH cell lines were obtained after ionizing radiation; specifically, 78, 64, 129, 55, 42, and 53 clones were isolated after treatment with 3,000, 4,000, 5,000, 6,000, 8,000, or 10,000 rads, respectively. A total of 421 hybrid cell lines were obtained by fusion with mouse B78 cells. Ninety-four microsatellite markers used in the genetic linkage map were genotyped using the 421 hybrid cell lines. Based upon marker retention and genome coverage, we selected 93 hybrid cell lines to form an RH panel. Importantly, we performed the first genotyping of yellowtail markers in an RH panel using a nanofluidic dynamic array (Fluidigm, CA, USA). Then, 580 markers containing ESTs and SNPs were mapped in the first yellowtail RH map. CONCLUSIONS: We successfully developed a yellowtail RH panel to facilitate the localization of markers. Using this, a framework RH map was constructed with 580 markers. This high-density physical map will serve as a useful tool for the identification of genes related to important breeding traits using genetic structural information, such as conserved synteny. Moreover, in a comparison of 30 sequences in the RH group 1 (SQ1), yellowtail appeared to be evolutionarily closer to medaka and the green-spotted pufferfish than to zebrafish. We suggest that synteny analysis may be potentially useful as a tool to investigate chromosomal evolution by comparison with model fish.

Quantitative trait loci (QTL) associated with resistance to a monogenean parasite (Benedenia seriolae) in yellowtail (Seriola quinqueradiata) through genome wide analysis.

Benedenia infections caused by the monogenean fluke ectoparasite Benedenia seriolae seriously impact marine finfish aquaculture. Genetic variation has been inferred to play a significant role in determining the susceptibility to this parasitic disease. To evaluate the genetic basis of Benedenia disease resistance in yellowtail (Seriola quinqueradiata), a genome-wide and chromosome-wide linkage analyses were initiated using F1 yellowtail families (n = 90 per family) based on a high-density linkage map with 860 microsatellite and 142 single nucleotide polymorphism (SNP) markers. Two major quantitative trait loci (QTL) regions on linkage groups Squ2 (BDR-1) and Squ20 (BDR-2) were identified. These QTL regions explained 32.9-35.5% of the phenotypic variance. On the other hand, we investigated the relationship between QTL for susceptibility to B. seriolae and QTL for fish body size. The QTL related to growth was found on another linkage group (Squ7). As a result, this is the first genetic evidence that contributes to detailing phenotypic resistance to Benedenia disease, and the results will help resolve the mechanism of resistance to this important parasitic infection of yellowtail.

Linkage mapping of toll-like receptors (TLRs) in Japanese flounder, Paralichthys olivaceus.

Toll-like receptors (TLRs) are responsible for the recognition of specific pathogen-associated molecular patterns and consequently activate signal pathways leading to inflammatory and interferon responses. The region surrounding several TLRs was previously found to be associated with resistance to specific disease. Hence, we determined the location of 11 TLRs in Japanese flounder (Paralichthys olivaceus) using polymorphic microsatellite markers. TLR1 and TLR3 were located on linkage group (LG) 21 and 7, respectively. Membrane TLR5 and soluble TLR5 were mapped to LG22. TLR7 and TLR8 were mapped to LG3. TLR9 was found on LG1 and TLR14 and TLR21 were located on the same linkage group, LG10. TLR22 was found on LG8. Interestingly, TLR2 was mapped with the previously reported Poli9-8TUF microsatellite marker which is tightly associated with lymphocystis virus disease resistance. Therefore, TLR2 is a candidate gene for resistance to lymphocystis disease. These results imply that the location of a TLR associated with a particular disease may be valuable for the research on the relationship between host immune response and disease resistance.

A second generation genetic linkage map of Japanese flounder (Paralichthys olivaceus).

BACKGROUND: Japanese flounder (Paralichthys olivaceus) is one of the most economically important marine species in Northeast Asia. Information on genetic markers associated with quantitative trait loci (QTL) can be used in breeding programs to identify and select individuals carrying desired traits. Commercial production of Japanese flounder could be increased by developing disease-resistant fish and improving commercially important traits. Previous maps have been constructed with AFLP markers and a limited number of microsatellite markers. In this study, improved genetic linkage maps are presented. In contrast with previous studies, these maps were built mainly with a large number of codominant markers so they can potentially be used to analyze different families and populations. RESULTS: Sex-specific genetic linkage maps were constructed for the Japanese flounder including a total of 1,375 markers [1,268 microsatellites, 105 single nucleotide polymorphisms (SNPs) and two genes]; 1,167 markers are linked to the male map and 1,067 markers are linked to the female map. The lengths of the male and female maps are 1,147.7 cM and 833.8 cM, respectively. Based on estimations of map lengths, the female and male maps covered 79 and 82% of the genome, respectively. Recombination ratio in the new maps revealed F:M of 1:0.7. All linkage groups in the maps presented large differences in the location of sex-specific recombination hot-spots. CONCLUSIONS: The improved genetic linkage maps are very useful for QTL analyses and marker-assisted selection (MAS) breeding programs for economically important traits in Japanese flounder. In addition, SNP flanking sequences were blasted against Tetraodon nigroviridis (puffer fish) and Danio rerio (zebrafish), and synteny analysis has been carried out. The ability to detect synteny among species or genera based on homology analysis of SNP flanking sequences may provide opportunities to complement initial QTL experiments with candidate gene approaches from homologous chromosomal locations identified in related model organisms.

Hyper-expansion of large DNA segments in the genome of kuruma shrimp, Marsupenaeus japonicus.

BACKGROUND: Higher crustaceans (class Malacostraca) represent the most species-rich and morphologically diverse group of non-insect arthropods and many of its members are commercially important. Although the crustacean DNA sequence information is growing exponentially, little is known about the genome organization of Malacostraca. Here, we constructed a bacterial artificial chromosome (BAC) library and performed BAC-end sequencing to provide genomic information for kuruma shrimp (Marsupenaeus japonicus), one of the most widely cultured species among crustaceans, and found the presence of a redundant sequence in the BAC library. We examined the BAC clone that includes the redundant sequence to further analyze its length, copy number and location in the kuruma shrimp genome. RESULTS: Mj024A04 BAC clone, which includes one redundant sequence, contained 27 putative genes and seemed to display a normal genomic DNA structure. Notably, of the putative genes, 3 genes encode homologous proteins to the inhibitor of apoptosis protein and 7 genes encode homologous proteins to white spot syndrome virus, a virulent pathogen known to affect crustaceans. Colony hybridization and PCR analysis of 381 BAC clones showed that almost half of the BAC clones maintain DNA segments whose sequences are homologous to the representative BAC clone Mj024A04. The Mj024A04 partial sequence was detected multiple times in the kuruma shrimp nuclear genome with a calculated copy number of at least 100. Microsatellites based BAC genotyping clearly showed that Mj024A04 homologous sequences were cloned from at least 48 different chromosomal loci. The absence of micro-syntenic relationships with the available genomic sequences of Daphnia and Drosophila suggests the uniqueness of these fragments in kuruma shrimp from current arthropod genome sequences. CONCLUSIONS: Our results demonstrate that hyper-expansion of large DNA segments took place in the kuruma shrimp genome. Although we analyzed only a part of the duplicated DNA segments, our result suggested that it is difficult to analyze the shrimp genome following normal analytical methodology. Hence, it is necessary to avoid repetitive sequence (such as segmental duplications) when studying the other unique structures in the shrimp genome.