Identification of the causative gene of a transparent phenotype of juvenile red sea bream Pagrus major.

Deformities in cultured fish species may be genetic, and identifying causative genes is essential to expand production and maintain farmed animal welfare. We previously reported a genetic deformity in juvenile red sea bream, designated a transparent phenotype. To identify its causative gene, we conducted genome-wide linkage analysis and identified two single nucleotide polymorphisms (SNP) located on LG23 directly linked to the transparent phenotype. The scaffold on which the two SNPs were located contained two candidate genes, duox and duoxa, which are related to thyroid hormone synthesis. Four missense mutations were found in duox and one in duoxa, with that in duoxa showing perfect association with the transparent phenotype. The mutation of duoxa was suggested to affect the transmembrane structure and thyroid-related traits, including an enlarged thyroid gland and immature erythrocytes, and lower thyroxine (T4) concentrations were observed in the transparent phenotype. The transparent phenotype was rescued by T4 immersion. Loss-of-function of duoxa by CRISPR-Cas9 induced the transparent phenotype in zebrafish. Evidence suggests that the transparent phenotype of juvenile red sea bream is caused by the missense mutation of duoxa and that this mutation disrupts thyroid hormone synthesis. The newly identified missense mutation will contribute to effective selective breeding of red sea bream to purge the causative gene of the undesirable phenotype and improve seed production of red sea bream as well as provide basic information of the mechanisms of thyroid hormones and its related diseases in fish and humans.

Identification, Characterization, and Mapping of a Novel SNP Associated with Body Color Transparency in Juvenile Red Sea Bream (Pagrus major).

We previously reported a body color deformity in juvenile red sea bream, which shows transparency in the juvenile stage because of delayed chromatophore development compared with normal individuals, and this finding suggested a genetic cause based on parentage assessments. To conduct marker-assisted selection to eliminate broodstock inheriting the causative gene, developing DNA markers associated with the phenotype was needed. We first conducted SNP mining based on AFLP analysis using bulked-DNA from normal and transparent individuals. One SNP was identified from a transparent-specific AFLP fragment, which significantly associated with transparent individuals. Two alleles (A/G) were observed in this locus, and the genotype G/G was dominantly observed in the transparent groups (97.1%) collected from several production lots produced from different broodstock populations. A few normal individuals inherited the G/G genotype (5.0%), but the A/A and A/G genotypes were dominantly observed in the normal groups. The homologs region of the SNP was searched using a medaka genome database, and intron 12 of the Nell2a gene (located on chromosome 6 of the medaka genome) was highly matched. We also mapped the red sea bream Nell2a gene on the previously developed linkage maps, and this gene was mapped on a male linkage group, LG4-M. The newly found SNP was useful in eliminating broodstock possessing the causative gene of the body color transparency observed in juvenile stage of red sea bream.

Identification of Quantitative Trait Loci for Resistance to RSIVD in Red Sea Bream (Pagrus major).

Red sea bream iridoviral disease (RSIVD) is a major viral disease in red sea bream farming in Japan. Previously, we identified one candidate male individual of red sea bream that was significantly associated with convalescent individuals after RSIVD. The purpose of this study is to identify the quantitative trait loci (QTL) linked to the RSIVD-resistant trait for future marker-assisted selection (MAS). Two test families were developed using the candidate male in 2014 (Fam-2014) and 2015 (Fam-2015). These test families were challenged with RSIV, and phenotypes were evaluated. Then, de novo genome sequences of red sea bream were obtained through next-generation sequencing, and microsatellite markers were searched and selected for linkage map construction. One immune-related gene, MHC class IIß, was also used for linkage map construction. Of the microsatellite markers searched, 148 and 197 were mapped on 23 and 27 linkage groups in the female and male linkage maps, respectively, covering approximately 65% of genomes in both sexes. One QTL linked to an RSIVD-resistant trait was found in linkage group 2 of the candidate male in Fam-2014, and the phenotypic variance of the QTL was 31.1%. The QTL was closely linked to MHC class IIß. Moreover, the QTL observed in Fam-2014 was also significantly linked to an RSIVD-resistant trait in the candidate male of Fam-2015. Our results suggest that the RSIVD-resistant trait in the candidate male was controlled by one major QTL closely linked to the MHC class IIß gene and could be useful for MAS of red sea bream.

DNA barcoding and morphological analyses revealed validity of Diadema clarki Ikeda, 1939 (Echinodermata, Echinoidea, Diadematidae).

A long-spined sea urchin Diadema-sp reported from Japanese waters was genetically distinct from all known Diadema species, but it remained undescribed. Extensive field surveys in Japan with molecular identification performed in the present study determined five phenotypes (I to V) in Diadema-sp according to the presence and/or shape of a white streak and blue iridophore lines in the naked space of the interambulacral area. All phenotypes were distinct from Diadema setosum (Leske, 1778) and Diadema savignyi (Audouin, 1829), of which a major type (I) corresponded to Diadema clarki Ikeda, 1939 that was questioned and synonymized with Diadema setosum by Mortensen (1940). The holotype of Diadema clarki has not been found, but three unlabeled dried tests of Diadema were found among Ikeda's original collection held in the Kitakyushu Museum of Natural History and Human History, Fukuoka, Japan. A short mtDNA COI fragment (ca. 350bp) was amplified from one of the tests, and the nucleotide sequence determined (275bp) was nearly identical with that of Diadema-sp. Arrangements of the primary tubercles on the coronal plates in Diadema-sp and the museum specimen also conformed with Diadema clarki, indicating that Diadema-sp is identical to Diadema clarki and a valid species. Narrow latitudinal distribution (31°N to 35°N) of Diadema clarki in Japan was observed, where it co-existed with abundant Diadema setosum and rare Diadema savignyi. No Diadema clarki was found in the southern islands in Japan, such as Satsunan Islands to Ryukyu Islands and Ogasawara Island, where Diadema setosum and Diadema savignyi were commonly observed.

Population genetic structure of the messmate pipefish Corythoichthys haematopterus in the northwest pacific: evidence for a cryptic species.

The population genetic structure of the messmate pipefish, Corythoichthys haematopterus, in the northwest Pacific was investigated based on the partial mitochondrial DNA cytochrome b (589 bp) and 16S rRNA (528 bp) region sequences of 108 individuals collected from six sites along the coast of the Japanese archipelago and one site on Mactan Island, the Philippines. A total of 60 and 28 haplotypes were obtained from the cytochrome b and 16S rRNA regions, respectively. Two genetically distinct lineages were detected: lineage A and B, which are separated by mean pairwise genetic distances of 23.3 and 14.1% in the partial cytochrome b and 16S rRNA genes, respectively. Such a huge genetic divergence between lineages, which is comparable to or even higher than the interspecific level, and the difference in their geographical distributions and habitat preferences suggests that they are distinct species, although there is no marked difference in their morphology. Haplotype network and gene and nucleotide diversity statistics indicate that the two lineages have different biogeographic histories: lineage A experienced rapid population expansion after a population bottleneck whereas lineage B has a long evolutionary history in a large stable population. In contrast, the levels of genetic variation among populations are relatively low in both lineages, probably because of frequent gene flow among populations resulting from the dispersal of pelagic larvae by the Kuroshio Current. These results indicate that past climatic events and contemporary oceanographic features have played a major role in establishing the population genetic structure of C. haematopterus.

DNA markers indicate low genetic diversity and high genetic divergence in the landlocked freshwater goby, Rhinogobius sp. YB, in the Ryukyu Archipelago, Japan.

Genetic diversity and genetic divergence were investigated in the landlocked goby Rhinogobius sp. YB by analysis of seven microsatellite DNA loci and the mtDNA control region sequence, and were compared with those of the closely related amphidromous species Rhinogobius sp. DA. Samples of Rhinogobius sp. YB and Rhinogobius sp. DA were collected from seven and four rivers, respectively. All pairwise Fst tests based on microsatellite DNA showed significant genetic differences, except for one pair of populations of Rhinogobius sp. DA (P<0.00064, alpha=78). The average Nei's genetic distance was 0.616 in Rhinogobius sp. YB and 0.394 in Rhinogobius sp. DA. Forty-two haplotypes were detected in both species, and almost all Rhinogobius sp. YB populations included different haplotypes. The means of allelic richness, Ho, and He in Rhinogobius sp. YB (2.057, 0.149, and 0.156, respectively) were significantly lower than in Rhinogobius sp. DA (4.868, 0.366, and 0.403, respectively; P<0.05). The high genetic divergence and low genetic diversity in Rhinogobius sp. YB may have resulted from repeated colonizations of rivers by different founders. Efforts to conserve genetic resources should take these evolutionarily significant units (ESU) of Rhinogobius sp. YB into account. The genetic markers used in this study provide simple and highly informative indicators for Rhinogobius sp. YB population management.

Isolation and characterization of 13 microsatellite markers for the viviparous surfperch Ditrema temmincki (Embiotocidae) and cross-species amplification.

Thirteen microsatellite loci were isolated from a size-selected genomic library of the surfperch (Ditrema temmincki Bleeker). All loci displayed a high degree of length polymorphism, as observed in the total number of alleles per locus (two to 23) and a high degree of estimated heterozygosity, ranging from 0.080 to 0.893. The primers developed for D. temmincki were also tested for their ability to amplify homologous sequences in D. viride and Neoditrema ransonetii. Distinct differences were observed among three species of surfperches, in both genetic variability and the frequency distribution of the alleles.