l-fucoside localization in the gills of the genus Takifugu and its possible implication in the parasitism of Heterobothrium okamotoi (Monogenea: Diclidophoridae).

BACKGROUND: The monogenean parasite Heterobothrium okamotoi only parasitizes the gills of Takifugu rubripes. In this study, we hypothesized that the carbohydrates contribute to high host specificity of H. okamotoi. METHODS: T. rubripes, T. niphobles, T. snyderi, and T. pardalis were used for UEA I staining of the gills and an in vivo challenge test against H. okamotoi. To examine the effect of l-fucose, an in vitro detachment test was conducted using the host's gills. Additionally, fucosylated proteins were isolated from the membrane proteins of T. niphobles gills. RESULTS: The location of l-fucoside and the infection dynamics in four species were correlated to some extent; H. okamotoi detached relatively quickly from T. niphobles possessing l-fucoside both on the surface of the gills and in certain types of cells, including mucus cells, but detached slowly from T. snyderi possessing l-fucoside in only certain types of cells, including mucus cells. Under the conditions examined, H. okamotoi exhibited minimal detachment from T. rubripes and T. pardalis, and l-fucoside was not detected. The significantly higher detachment rate of H. okamotoi from the host's gills incubated in l-fucose-containing medium compared with the controls suggests that l-fucose in the non-host gills induced detachment of H. okamotoi. Four fucosylated proteins, including mucin5AC-like, were identified as potential factors for the detachment of H. okamotoi. CONCLUSIONS: Fucosylated proteins covering the surface of non-host gills might contribute to H. okamotoi detachment. GENERAL SIGNIFICANCE: This research shows the possible involvement of oligosaccharides in the host specificity of monogenean parasites.

Repeated translocation of a supergene underlying rapid sex chromosome turnover in Takifugu pufferfish.

Recent studies have revealed a surprising diversity of sex chromosomes in vertebrates. However, the detailed mechanism of their turnover is still elusive. To understand this process, it is necessary to compare closely related species in terms of sex-determining genes and the chromosomes harboring them. Here, we explored the genus Takifugu, in which one strong candidate sex-determining gene, Amhr2, has been identified. To trace the processes involved in transitions in the sex-determination system in this genus, we studied 12 species and found that while the Amhr2 locus likely determines sex in the majority of Takifugu species, three species have acquired sex-determining loci at different chromosomal locations. Nevertheless, the generation of genome assemblies for the three species revealed that they share a portion of the male-specific supergene that contains a candidate sex-determining gene, GsdfY, along with genes that potentially play a role in male fitness. The shared supergene spans ∼100 kb and is flanked by two duplicated regions characterized by CACTA transposable elements. These results suggest that the shared supergene has taken over the role of sex-determining locus from Amhr2 in lineages leading to the three species, and repeated translocations of the supergene underlie the turnover of sex chromosomes in these lineages. These findings highlight the underestimated role of a mobile supergene in the turnover of sex chromosomes in vertebrates.

Reproductive Ability of Hybrids between Japanese Pond Turtle (Mauremys japonica) and Reeves' Pond Turtle (Mauremys reevesii).

Hybridization induced by human activities, such as crossbreeding between invasive and native species, can adversely affect the natural biodiversity of an ecosystem. In Japan, the endemic turtle species Mauremys japonica is known to hybridize with the alien species Mauremys reevesii, and putative hybrids have been encountered in the wild. If M. japonica × M. reevesii hybrids can readily crossbreed with M. japonica, the hybridization with M. reevesii could lead to the extinction of pure M. japonica populations. However, information on the reproductive ability of M. japonica × M. reevesii hybrids is limited. In this study, we collected wild-caught hybrids from across western Japan to assess their reproductive ability. We investigated the nesting season timing, clutch size, embryonic development, hatching success, and sperm viability. The results showed that female hybrids nested during the same months as the parental species and had similar clutch sizes and hatching success. No embryonic development abnormalities were detected, and viable sperm were observed in all hybrid male semen samples. In conclusion, the fertility of M. japonica × M. reevesii hybrids appears to be similar to the fertilities of the parental species, posing a potential challenge for M. japonica conservation.

Genomic prediction for testes weight of the tiger pufferfish, Takifugu rubripes, using medium to low density SNPs.

Aquaculture production is expected to increase with the help of genomic selection (GS). The possibility of performing GS using only a small number of SNPs has been examined in order to reduce genotyping costs; however, the practicality of this approach is still unclear. Here, we tested whether the effects of reducing the number of SNPs impaired the prediction accuracy of GS for standard length, body weight, and testes weight in the tiger pufferfish (Takifugu rubripes). High values for predictive ability (0.563-0.606) were obtained with 4000 SNPs for all traits under a genomic best linear unbiased predictor (GBLUP) model. These values were still within an acceptable range with 1200 SNPs (0.554-0.588). However, predictive abilities and prediction accuracies deteriorated using less than 1200 SNPs largely due to the reduced power in accurately estimating the genetic relationship among individuals; family structure could still be resolved with as few as 400 SNPs. This suggests that the SNPs informative for estimation of genetic relatedness among individuals differ from those for inference of family structure, and that non-random SNP selection based on the effects on family structure (e.g., site-FST, principal components, or random forest) is unlikely to increase the prediction accuracy for these traits.

Genetic Dissection of a Precocious Phenotype in Male Tiger Pufferfish (Takifugu rubripes) using Genotyping by Random Amplicon Sequencing, Direct (GRAS-Di).

The novel non-targeted PCR-based genotyping system, namely Genotyping by Random Amplicon Sequencing, Direct (GRAS-Di), is characterized by the simplicity in library construction and robustness against DNA degradation and is expected to facilitate advancements in genetics, in both basic and applied sciences. In this study, we tested the utility of GRAS-Di for genetic analysis in a cultured population of the tiger pufferfish Takifugu rubripes. The genetic analyses included family structure analysis, genetic map construction, and quantitative trait locus (QTL) analysis for the male precocious phenotype using a population consisting of four full-sib families derived from a genetically precocious line. An average of 4.7 million raw reads were obtained from 198 fish. Trimmed reads were mapped onto a Fugu reference genome for genotyping, and 21,938 putative single-nucleotide polymorphisms (SNPs) were obtained. These 22 K SNPs accurately resolved the sibship and parent-offspring pairs. A fine-scale linkage map (total size: 1,949 cM; average interval: 1.75 cM) was constructed from 1,423 effective SNPs, for which the allele inheritance patterns were known. QTL analysis detected a significant locus for testes weight on Chr_14 and three suggestive loci on Chr_1, Chr_8, and Chr_19. The significant QTL was shared by body length and body weight. The effect of each QTL was small (phenotypic variation explained, PVE: 3.1-5.9%), suggesting that the precociousness seen in the cultured pufferfish is polygenic. Taken together, these results indicate that GRAS-Di is a practical genotyping tool for aquaculture species and applicable for molecular breeding programs, such as marker-assisted selection and genomic selection.

Integrative genomic phylogeography reveals signs of mitonuclear incompatibility in a natural hybrid goby population.

Hybridization between divergent lineages generates new allelic combinations. One mechanism that can hinder the formation of hybrid populations is mitonuclear incompatibility, that is, dysfunctional interactions between proteins encoded in the nuclear and mitochondrial genomes (mitogenomes) of diverged lineages. Theoretically, selective pressure due to mitonuclear incompatibility can affect genotypes in a hybrid population in which nuclear genomes and mitogenomes from divergent lineages admix. To directly and thoroughly observe this key process, we de novo sequenced the 747-Mb genome of the coastal goby, Chaenogobius annularis, and investigated its integrative genomic phylogeographics using RNA-sequencing, RAD-sequencing, genome resequencing, whole mitogenome sequencing, amplicon sequencing, and small RNA-sequencing. Chaenogobius annularis populations have been geographically separated into Pacific Ocean (PO) and Sea of Japan (SJ) lineages by past isolation events around the Japanese archipelago. Despite the divergence history and potential mitonuclear incompatibility between these lineages, the mitogenomes of the PO and SJ lineages have coexisted for generations in a hybrid population on the Sanriku Coast. Our analyses revealed accumulation of nonsynonymous substitutions in the PO-lineage mitogenomes, including two convergent substitutions, as well as signals of mitochondrial lineage-specific selection on mitochondria-related nuclear genes. Finally, our data implied that a microRNA gene was involved in resolving mitonuclear incompatibility. Our integrative genomic phylogeographic approach revealed that mitonuclear incompatibility can affect genome evolution in a natural hybrid population.

Genomic selection for heterobothriosis resistance concurrent with body size in the tiger pufferfish, Takifugu rubripes.

Parasite resistance traits in aquaculture species often have moderate heritability, indicating the potential for genetic improvements by selective breeding. However, parasite resistance is often synonymous with an undesirable negative correlation with body size. In this study, we first tested the feasibility of genomic selection (GS) on resistance to heterobothriosis, caused by the monogenean parasite Heterobothrium okamotoi, which leads to huge economic losses in aquaculture of the tiger pufferfish Takifugu rubripes. Then, using a simulation study, we tested the possibility of simultaneous improvement of parasite resistance, assessed by parasite counts on host fish (HC), and standard length (SL). Each trait showed moderate heritability (square-root transformed HC: h2 = 0.308 ± 0.123, S.E.; SL: h2 = 0.405 ± 0.131). The predictive abilities of genomic prediction among 12 models, including genomic Best Linear Unbiased Predictor (GBLUP), Bayesian regressions, and machine learning procedures, were also moderate for both transformed HC (0.248‒0.344) and SL (0.340‒0.481). These results confirmed the feasibility of GS for this trait. Although an undesirable genetic correlation was suggested between transformed HC and SL (rg = 0.228), the simulation study suggested the desired gains index can help achieve simultaneous genetic improvements in both traits.

The Genetic Basis of Scale-Loss Phenotype in the Rapid Radiation of Takifugu Fishes.

Rapid radiation associated with phenotypic divergence and convergence provides an opportunity to study the genetic mechanisms of evolution. Here we investigate the genus Takifugu that has undergone explosive radiation relatively recently and contains a subset of closely-related species with a scale-loss phenotype. By using observations during development and genetic mapping approaches, we show that the scale-loss phenotype of two Takifugu species, T. pardalis Temminck & Schlegel and T. snyderi Abe, is largely controlled by an overlapping genomic segment (QTL). A search for candidate genes underlying the scale-loss phenotype revealed that the QTL region contains no known genes responsible for the evolution of scale-loss phenotype in other fishes. These results suggest that the genes used for the scale-loss phenotypes in the two Takifugu are likely the same, but the genes used for the similar phenotype in Takifugu and distantly related fishes are not the same. Meanwhile, Fgfrl1, a gene predicted to function in a pathway known to regulate bone/scale development was identified in the QTL region. Since Fgfr1a1, another memebr of the Fgf signaling pathway, has been implicated in scale loss/scale shape in fish distantly related to Takifugu, our results suggest that the convergence of the scale-loss phenotype may be constrained by signaling modules with conserved roles in scale development.

A highly flexible and repeatable genotyping method for aquaculture studies based on target amplicon sequencing using next-generation sequencing technology.

Studies using genome-wide single nucleotide polymorphisms (SNPs) have become commonplace in genetics and genomics, due to advances in high-throughput sequencing technologies. Since the numbers of required SNPs and samples vary depending on each research goal, genotyping technologies with high flexibility in the number of SNPs/samples and high repeatability have been intensively investigated. For example, the ultrahigh-multiplexed amplicon sequencing, Ion AmpliSeq, has been used as a high-throughput genotyping method mainly for diagnostic purposes. Here, we designed a custom panel targeting 3,187 genome-wide SNPs of fugu, Takifugu rubripes, and applied it for genotyping farmed fugu to test its feasibility in aquaculture studies. We sequenced two libraries consisting of different pools of individuals (n = 326 each) on the Illumina MiSeq sequencer. Consequently, over 99% target regions (3,178 SNPs) were amplified and 2,655 SNPs were available after filtering steps. Strong correlation was observed in the mean depth of coverage of each SNP between duplicate runs (r = 0.993). Genetic analysis using these genotype data successfully detected the known population structure and the sex determining locus of fugu. These results show the method is superior in repeatability and flexibility, and suits genetic studies including molecular breeding, such as marker assisted and genomic selection.

A SNP in a Steroidogenic Enzyme Is Associated with Phenotypic Sex in Seriola Fishes.

Vertebrate sex development consists largely of two processes: "sex determination," the initial bifurcation of sexual identity, and "sex differentiation," which subsequently facilitates maleness or femaleness according to the sex determination signal. Steroid hormones promote multiple types of sexual dimorphism in eutherian mammals and avians [1-3], in which they are indispensable for proper sex differentiation. By contrast, in many poikilothermic vertebrates, steroid hormones have been proposed to be key players in sex determination as well as sex differentiation [4-8]. This hypothesis was introduced more than 50 years ago but has never been rigorously tested due to difficulties in discriminating the roles of steroids in sex determination and differentiation. We found that a missense SNP in the gene encoding the steroidogenic enzyme 17ß-hydroxysteroid dehydrogenase 1 (Hsd17b1) was perfectly associated with ZZ/ZW sex determination in Seriola fishes. Biochemical analyses revealed that a glutamate residue present specifically in Z-type HSD17B1 attenuated interconversion between 17-keto and 17ß-hydroxy steroids relative to the allelic product from the W chromosome, which harbors glycine at that position, by disrupting the hydrogen bond network between the steroid and the enzyme's catalytic residues. Hsd17b1 mRNA is constitutively expressed in undifferentiated and differentiating gonads of both genotypic sexes, whereas W-type mRNA is expressed only in genotypic females. Meanwhile, Cyp19a1 is predominantly expressed in differentiating ovary. We conclude that the combination of Hsd17b1 alleles determines sex by modulating endogenous estrogen levels in Seriola species. These findings strongly support the long-standing hypothesis on steroids in sex determination.

Random PCR-based genotyping by sequencing technology GRAS-Di (genotyping by random amplicon sequencing, direct) reveals genetic structure of mangrove fishes.

While various technologies for high-throughput genotyping have been developed for ecological studies, simple methods tolerant to low-quality DNA samples are still limited. In this study, we tested the availability of a random PCR-based genotyping-by-sequencing technology, genotyping by random amplicon sequencing, direct (GRAS-Di). We focused on population genetic analysis of estuarine mangrove fishes, including two resident species, the Amboina cardinalfish (Fibramia amboinensis, Bleeker, 1853) and the Duncker's river garfish (Zenarchopterus dunckeri, Mohr, 1926), and a marine migrant, the blacktail snapper (Lutjanus fulvus, Forster, 1801). Collections were from the Ryukyu Islands, southern Japan. PCR amplicons derived from ~130 individuals were pooled and sequenced in a single lane on a HiSeq2500 platform, and an average of three million reads was obtained per individual. Consensus contigs were assembled for each species and used for genotyping of single nucleotide polymorphisms by mapping trimmed reads onto the contigs. After quality filtering steps, 4,000-9,000 putative single nucleotide polymorphisms were detected for each species. Although DNA fragmentation can diminish genotyping performance when analysed on next-generation sequencing technology, the effect was small. Genetic differentiation and a clear pattern of isolation-by-distance was observed in F. amboinensis and Z. dunckeri by means of principal component analysis, FST and the admixture analysis. By contrast, L. fulvus comprised a genetically homogeneous population with directional recent gene flow. These genetic differentiation patterns reflect patterns of estuary use through life history. These results showed the power of GRAS-Di for fine-grained genetic analysis using field samples, including mangrove fishes.

Assessment of genetic diversity in Coho salmon (Oncorhynchus kisutch) populations with no family records using ddRAD-seq.

OBJECTIVE: Selective breeding for desirable traits is becoming popular in aquaculture. In Miyagi prefecture, Japan, a selectively bred population of Coho salmon (Oncorhynchus kisutch) has been established with the original, randomly breeding population maintained separately. Since they have been bred without family records, the genetic diversity within these populations remains unknown. In this study, we estimated the genetic diversity and key quantitative genetic parameters such as heritability and genomic breeding value for body size traits by means of genomic best linear unbiased prediction to assess the genetic health of these populations. RESULTS: Ninety-nine and 83 females from the selective and random groups, respectively, were genotyped at 2350 putative SNPs by means of double digest restriction associated DNA sequencing. The genetic diversity in the selectively bred group was low, as were the estimated heritability and prediction accuracy for length and weight (h2 = 0.26-0.28; accuracy = 0.34), compared to the randomly bred group (h2 = 0.50-0.60; accuracy = 0.51-0.54). Although the tested sample size was small, these results suggest that further selection is difficult for the selectively bred population, while there is some potential for the randomly bred group, especially with the aid of genomic information.

Identification of the sex-determining locus in grass puffer (Takifugu niphobles) provides evidence for sex-chromosome turnover in a subset of Takifugu species.

There is increasing evidence for frequent turnover in sex chromosomes in vertebrates. Yet experimental systems suitable for tracing the detailed process of turnover are rare. In theory, homologous turnover is possible if the new sex-determining locus is established on the existing sex-chromosome. However, there is no empirical evidence for such an event. The genus Takifugu includes fugu (Takifugu rubripes) and its two closely-related species whose sex is most likely determined by a SNP at the Amhr2 locus. In these species, males are heterozygous, with G and C alleles at the SNP site, while females are homozygous for the C allele. To determine if a shift in the sex-determining locus occurred in another member of this genus, we used genetic mapping to characterize the sex-chromosome systems of Takifugu niphobles. We found that the G allele of Amhr2 is absent in T. niphobles. Nevertheless, our initial mapping suggests a linkage between the phenotypic sex and the chromosome 19, which harbors the Amhr2 locus. Subsequent high-resolution analysis using a sex-reversed fish demonstrated that the sex-determining locus maps to the proximal end of chromosome 19, far from the Amhr2 locus. Thus, it is likely that homologous turnover involving these species has occurred. The data also showed that there is a male-specific reduction of recombination around the sex-determining locus. Nevertheless, no evidence for sex-chromosome differentiation was detected: the reduced recombination depended on phenotypic sex rather than genotypic sex; no X- or Y-specific maker was obtained; the YY individual was viable. Furthermore, fine-scale mapping narrowed down the new sex-determining locus to the interval corresponding to approximately 300-kb of sequence in the fugu genome. Thus, T. niphobles is determined to have a young and small sex-determining region that is suitable for studying an early phase of sex-chromosome evolution and the mechanisms underlying turnover of sex chromosome.

Mucosal IgM Antibody with d-Mannose Affinity in Fugu Takifugu rubripes Is Utilized by a Monogenean Parasite Heterobothrium okamotoi for Host Recognition.

How parasites recognize their definitive hosts is a mystery; however, parasitism is reportedly initiated by recognition of certain molecules on host surfaces. Fish ectoparasites make initial contact with their hosts at body surfaces, such as skin and gills, which are covered with mucosa that are similar to those of mammalian guts. Fish are among the most primitive vertebrates with immune systems that are equivalent to those in mammals, and they produce and secrete IgM into mucus. In this study, we showed that the monogenean parasite Heterobothrium okamotoi utilizes IgM to recognize its host, fugu Takifugu rubripes Oncomiracidia are infective larvae of H. okamotoi that shed their cilia and metamorphose into juveniles when exposed to purified d-mannose-binding fractions from fugu mucus. Using liquid chromatography-tandem mass spectrometry analysis, proteins contained in the fraction were identified as d-mannose-specific IgM with two d-mannose-binding lectins. However, although deciliation was significantly induced by IgM and was inhibited by d-mannose or a specific Ab against fugu IgM, other lectins had no effect, and IgM without d-mannose affinity induced deciliation to a limited degree. Subsequent immunofluorescent staining experiments showed that fugu d-mannose-specific IgM binds ciliated epidermal cells of oncomiracidium. These observations suggest that deciliation is triggered by binding of fugu IgM to cell surface Ags via Ag binding sites. Moreover, concentrations of d-mannose-binding IgM in gill mucus were sufficient to induce deciliation in vitro, indicating that H. okamotoi parasites initially use host Abs to colonize host gills.

Genetic Basis Underlying Behavioral Correlation Between Fugu Takifugu rubripes and a Closely Related Species, Takifugu niphobles.

Correlated suits of behaviors (behavioral syndrome) are commonly observed in both inter- and intraspecific studies. In order to understand the genetic basis of such a correlation between species, we compared ten behaviors classified into five categories (acclimation, feeding, normal swimming, reaction to a novel object and activity in a novel environment) between two pufferfish species, Takifugu rubripes and T. niphobles. The two species showed consistent differences in nine behaviors with a significant correlation among behaviors. Quantitative trait locus (QTL) analysis using second generation hybrids revealed that different sets of small effect QTL are associated with the observed interspecific behavioral disparity. This indicates that correlations in temperament traits between them are governed by many genes with small effects, and each behavior has been selected to form particular combination patterns. One of the QTL showing small pleiotropic effect includes the Drd4 gene known for its association with behavioral traits in some animal taxa including mammals.

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.

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 F(2) progenies of their F(1) 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.

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.

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.

Cortisol response and immune-related effects of Atlantic salmon (Salmo salar Linnaeus) subjected to short- and long-term stress.

It is generally considered that stress causes decreased immune function in fish. In this study we examined in Atlantic salmon (Salmo salar Linnaeus) the effects of both short- (a single 15s out of water) and long-term (4 weeks of daily handling 15s out of water) stress on plasma cortisol (free and total) and glucose levels, expression of interleukin-1beta (IL-1beta) and survival of head kidney (HK) macrophages under culture with Aeromonas salmonicida. In the short-term study, samples were collected prior to the application of the stressor, and at 1, 3, 6, 12 and 24h post stress. Free and total plasma cortisol levels and the percentage of free cortisol increased significantly in the stressed group at 1 and 3h post stress. Plasma glucose levels were significantly higher than those of control fish at 1, 3 and 6h post stress. Constitutive expression of IL-1beta in macrophages isolated from head kidneys in stressed fish was significantly higher at 1 and 3h post stress. However, lipopolysaccharide (LPS) stimulated expression of IL-1beta in HK macrophages, exhibited significantly higher fold increases in unstressed fish compared to stressed fish. In the long-term study, with the exception of an increase in plasma glucose levels at 1 week, there were no significant differences in stress parameters between groups. There was a significantly higher constitutive IL-1beta expression in macrophages isolated from stressed fish over the first 2 weeks. At weeks 1, 2 and 3 the magnitude of IL-1beta response of isolated HK macrophages to LPS stimulation was reduced in >90% of the stressed fish. At 4 weeks there was no significant difference in inducible IL-1beta expression between the groups. Macrophages isolated from stressed fish also showed significantly decreased survival when exposed to A. salmonicida. This study shows a clear pattern from repeated handling stress, whereby effects on immune cells begin with increased constitutive expression of IL-1beta, followed by decreased stimulation of leucocytes by extracellular antigen, and finally decreased leukocyte survival when exposed to A. salmonicida. The implications of these changes in the immune system will be discussed with respect to the use of classical indicators of stress to predict possible effects on the immune system of fish.