A genetic linkage map for the tiger pufferfish, Takifugu rubripes

Feeding Strategy to Use Beef Tallow and Modify Farmed Tiger Puffer Fatty Acid Composition

A 12-week feeding experiment was conducted to evaluate the effects of replacing fish oil (FO) with beef tallow (BT) on the fatty acid composition of farmed tiger puffer (Takifugu rubripes). Two replacement strategies were used: a standard Graded Dietary Replacement of FO with BT (GDR strategy) and Alternate Feeding between FO- and BT-based Diets (AFD strategy). The positive and negative control diets were formulated with 6% FO (FO-C group) or BT (BT-C group) as the sole added lipid source. In the GDR strategy, three experimental diets were formulated, with 25, 50 and 75% of the added FO in the FO-C diet replaced with BT, named 25BT, 50BT and 75BT, respectively. In the AFD strategy, alternated feeding patterns between the FO-C and BT-C diet-namely, 1, 2 and 3 weeks with BT-C followed by 1 week feeding with FO-C (1BT-1FO, 2BT-1FO and 3BT-1FO, respectively)-were applied. Each diet or feeding strategy was assigned to triplicate tanks. The results showed that dietary BT inclusion reduced the contents of long-chain polyunsaturated fatty acids (LC-PUFA) in both the muscle and liver (edible tissues for this species) of the experimental fish, and the liver displayed a more drastic decrease than the muscle. The LC-PUFA content linearly decreased with the decreasing dietary FO levels in the GDR strategy. However, in the AFD strategy, a linear relationship was not observed between the LC-PUFA content and the FO feeding duration. The 3BT-1FO treatment resulted in higher LC-PUFA content than 2BT-1FO. When comparing the two strategies with the same final FO administration level-namely, 50BT vs. 1BT-1FO, and in particular, 75BT vs. 3BT-1FO-the AFD strategy resulted in higher LC-PUFA contents in both the muscle and liver than the GDR strategy. In conclusion, when FO was replaced with BT in the diets, alternate feeding between FO- and BT-based diets resulted in a higher LC-PUFA content than the standard direct replacement. Three weeks of feeding with BT-C followed by one week of feeding with FO-C appeared to be a good alternate feeding pattern. This study provided a promising strategy of FO-sparing in fish farming when the LC-PUFA contents were maintained as high as possible.

Comparative genomics analysis and genome assembly integration with the recombination landscape contribute to Takifugu bimaculatus assembly refinement

Takifugu genus has been brought to the fore in scientific and practical research due to its compact genome, explosive speciation progress and economic value. Here we updated the chromosome-level genome of Takifugu bimaculatus by an ultra-high-density linkage map, a classic and accurate way of chromosome assembly. The map constituted a robust assembly frame, with 92.2% (372.77 Mb) of the draft genome cumulatively placed. With intraspecies and interspecies comparative genomic analysis, we developed a criterion to quantify the differences between assemblies and established a novel way to integrate information from multiple assemblies. The integrated assembly rectified potential mis-assemblies, greatly improving the genome contiguity and correctness. Our results rendered profound information on the genetic recombination of T. bimaculatus and provided new insights into effective genome assembly. The consolidated assembly will be a contributory tool of T. bimaculatus and broadly across the Takifugu by providing a convincing reference for genomic research.

HSP70 gene expression responses to the temperature stress in pufferfish (Takifugu rubripes)

In this study, we isolated and characterized HSP70 cDNA from pufferfish (Takifugu rubripes). The 3053 bp full-length TrHSP70 sequence consisted of a 167 bp 5'-UTR (untranslated region), a 2535 bp open reading frame, and a 351 bp 3'-UTR. BLAST analysis revealed that the TrHSP70 shared high similarity with HSP70 sequences in other species. In our study, we set 3 experimental groups as H1 group (20 °C), H2 group (24 °C), and H3 group (28 °C) for checking the expression level of TrHSP70 in T. rubripes. Tissue-specific gene expression results showed that TrHSP70 had higher expression in the intestines than other tissues of the T. rubripes by RT-qPCR. In the experimental group, we found that the expression of TrHSP70 was upregulated in different tissues in the H3 group. The results show that TrHSP70 is a constitutively expressed gene, which plays an important role in maintaining normal physiological function and coping with stress.

Genome-Scale Phylogenetic and Population Genetic Studies Provide Insight Into Introgression and Adaptive Evolution of Takifugu Species in East Asia

As a typical marine adaptive radiation species, most Takifugu species are widely distributed in East Asian offshore, which have diversified morphological characteristics and different ecological habits. The phylogenetic relationship and population structure of the Takifugu species was complicated because of incomplete lineage sorting, widespread hybridization and introgression. Therefore, to systematically clarify the phylogenetic relationships of Takifugu genus, explore the introgression and natural hybridization between different Takifugu species, and detect the selective signatures in the adaptive evolution of diversified traits, whole-genome resequencing was used in 122 Takifugu samples from 10 species. Phylogenetic analysis showed solid sister-group relationships between Takifugu bimaculatus and Takifugu flavidus, Takifugu oblongus, and Takifugu niphobles, Takifugu rubripes, and Takifugu obscurus, Takifugu xanthoptreus, and Takifugu ocellatus. Further admixture analysis indicated the divergence of T. obscurus population and the bidirectional gene flow between T. bimaculatus and T. flavidus. Using species-specific homozygous genetic variance sites, we detected the asymmetric introgression between T. bimaculatus and T. flavidus at China East sea and southern Taiwan Strait. By genome-scale genetic diversity scanning, we detected two copies of syt1, zar1 and tgfbr1 related to the semilunar reproduction rhythm in T. niphobles, involved in memory formation, embryo maturation and female reproduction. Furthermore, we also found lots of T. niphobles specific mutations in CDS region of circadian rhythm related genes and endocrine hormone genes. For Takifugu species, our research provides reliable genetic resources and results for the phylogeny, introgression, hybridization and adaptive evolution, and could be used as a guide for the formulation of the protection and proliferation release policies.

Dietary taurine stimulates the hepatic biosynthesis of both bile acids and cholesterol in the marine teleost, tiger puffer (Takifugu rubripes)

Taurine (TAU) plays important roles in the metabolism of bile acids, cholesterol and lipids. However, little relevant information has been available in fish where TAU has been identified as a conditionally essential nutrient. The present study aimed to investigate the effects of dietary TAU on the metabolism of bile acids, cholesterol and lipids in tiger puffer, which is both an important aquaculture species and a good research model, having a unique lipid storage pattern. An 8-week feeding trial was conducted in a flow-through seawater system. Three experimental diets differed only in TAU level, that is, 1·7, 8·2 and 14·0 mg/kg. TAU supplementation increased the total bile acid content in liver but decreased the content in serum. TAU supplementation also increased the contents of total cholesterol and HDL-cholesterol in both liver and serum. The hepatic bile acid profile mainly includes taurocholic acid (94·48 %), taurochenodeoxycholic acid (4·17 %) and taurodeoxycholic acid (1·35 %), and the contents of all these conjugated bile acids were increased by dietary TAU. The hepatic lipidomics analysis showed that TAU tended to decrease the abundance of individual phospholipids and increase those of some individual TAG and ceramides. The hepatic mRNA expression study showed that TAU stimulated the biosynthesis of both bile acids and cholesterol, possibly via regulation of farnesoid X receptor and HDL metabolism. TAU also stimulated the hepatic expression of lipogenic genes. In conclusion, dietary TAU stimulated the hepatic biosynthesis of both bile acids and cholesterol and tended to regulate lipid metabolism in multiple ways.

Toll-like receptors (TLRs) respond to tributyltin chloride (TBT-Cl) exposure in the river pufferfish (Takifugu obscurus): Evidences for its toxic injury function

Tributyltin chloride (TBT-Cl) residual in water body had become a noticeable ecological problem for aquatic ecosystems. Toll-like receptors (TLRs) are an ancient family of pattern recognition receptors that play key roles in detecting nonself antigens and immune system activation. In this study, we explored the effect of TBT-Cl exposure on four TLRs expression in river pufferfish, Takifugu obscurus. The four T. obscurus Toll-like receptors (To-TLRs) contained different types of domains such as leucine-rich repeats (LRRs), leucine-rich repeats, typical subfamily (LRR_TYP) and other special domains. The To-TLRs mRNA transcripts expressed in all tissues, also To-TLR2 was investigated with higher level in kidney, as well as To-TLR3 in kidney, while To-TLR18 in liver and To-TLR22 in intestine. After the acute and chronic exposure of TBT-Cl, To-TLR2 and To-TLR3 mRNA transcripts were significantly down-regulated in gill. However, To-TLR18 and To-TLR22 were significantly up-regulated in gill and liver. Moreover, the histology and immunohistochemistry (IHC) results showed the different injury degrees of TBT-Cl in liver and gill and implied the cytoplasm reorganization after TBT-Cl stress and the function of immunoregulation for To-TLRs to TBT-Cl exposure. All the results indicated that To-TLRs might involve in sensing and mediating innate immune responses caused by TBT-Cl for keeping detoxification homeostasis.

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.

Improvement of the Pacific bluefin tuna (Thunnus orientalis) reference genome and development of male-specific DNA markers

The Pacific bluefin tuna, Thunnus orientalis, is a highly migratory species that is widely distributed in the North Pacific Ocean. Like other marine species, T. orientalis has no external sexual dimorphism; thus, identifying sex-specific variants from whole genome sequence data is a useful approach to develop an effective sex identification method. Here, we report an improved draft genome of T. orientalis and male-specific DNA markers. Combining PacBio long reads and Illumina short reads sufficiently improved genome assembly, with a 38-fold increase in scaffold contiguity (to 444 scaffolds) compared to the first published draft genome. Through analysing re-sequence data of 15 males and 16 females, 250 male-specific SNPs were identified from more than 30 million polymorphisms. All male-specific variants were male-heterozygous, suggesting that T. orientalis has a male heterogametic sex-determination system. The largest linkage disequilibrium block (3,174 bp on scaffold_064) contained 51 male-specific variants. PCR primers and a PCR-based sex identification assay were developed using these male-specific variants. The sex of 115 individuals (56 males and 59 females; sex was diagnosed by visual examination of the gonads) was identified with high accuracy using the assay. This easy, accurate, and practical technique facilitates the control of sex ratios in tuna farms. Furthermore, this method could be used to estimate the sex ratio and/or the sex-specific growth rate of natural populations.

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.

Exaggerated heterochiasmy in a fish with sex-linked male coloration polymorphisms

It is often stated that polymorphisms for mutations affecting fitness of males and females in opposite directions [sexually antagonistic (SA) polymorphisms] are the main selective force for the evolution of recombination suppression between sex chromosomes. However, empirical evidence to discriminate between different hypotheses is difficult to obtain. We report genetic mapping results in laboratory-raised families of the guppy (Poecilia reticulata), a sexually dimorphic fish with SA polymorphisms for male coloration genes, mostly on the sex chromosomes. Comparison of the genetic and physical maps shows that crossovers are distributed very differently in the two sexes (heterochiasmy); in male meiosis, they are restricted to the termini of all four chromosomes studied, including chromosome 12, which carries the sex-determining locus. Genome resequencing of male and female guppies from a population also indicates sex linkage of variants across almost the entire chromosome 12. More than 90% of the chromosome carrying the male-determining locus is therefore transmitted largely through the male lineage. A lack of heterochiasmy in a related fish species suggests that it originated recently in the lineage leading to the guppy. Our findings do not support the hypothesis that suppressed recombination evolved in response to the presence of SA polymorphisms. Instead, a low frequency of recombination on a chromosome that carries a male-determining locus and has not undergone genetic degeneration has probably facilitated the establishment of male-beneficial coloration polymorphisms.

Effect of tributyltin chloride (TBT-Cl) exposure on expression of HSP90β1 in the river pufferfish (Takifugu obscurus): Evidences for its immunologic function involving in exploring process

HSP90β1 (known as glyco-protein 96, GP96) is a vital endoplasmic reticulum (ER) depended chaperonin among the HSPs (heat shock proteins) family. Furthermore, it always processes and presents antigen of the tumor and keeps balance for the intracellular environment. In the present study, we explored the effect of tributyltin chloride (TBT-Cl) exposure on HSP90β1 expression in river pufferfish, Takifugu obscurus. The full length of To-HSP90β1 was gained with 2775 bp in length, with an ORF (open reading frame) encoding an 803 aa polypeptide. A phylogenetic tree was constructed and showed the close relationship to other fish species. The HSP90β1 mRNA transcript was expressed in all tissues investigated with higher level in the gill and liver. After the acute and chronic exposure of TBT-Cl, the To-HSP90β1 mRNA transcript significantly was up-regulated in gills. Moreover, the histology study indicated the different injury degree of TBT-Cl in liver and gill. Immunohistochemistry (IHC) staining results implied the cytoplasm reorganization after TBT-Cl stress and the function of immunoregulation for To-HSP90β1 to TBT-Cl exposure. All the results indicated that HSP90β1 may be involved in the resistance to the invasion of TBT-Cl for keeping autoimmune homeostasis.

The Guppy Sex Chromosome System and the Sexually Antagonistic Polymorphism Hypothesis for Y Chromosome Recombination Suppression

Sex chromosomes regularly evolve suppressed recombination, distinguishing them from other chromosomes, and the reason for this has been debated for many years. It is now clear that non-recombining sex-linked regions have arisen in different ways in different organisms. A major hypothesis is that a sex-determining gene arises on a chromosome and that sexually antagonistic (SA) selection (sometimes called intra-locus sexual conflict) acting at a linked gene has led to the evolution of recombination suppression in the region, to reduce the frequency of low fitness recombinant genotypes produced. The sex chromosome system of the guppy (Poecilia reticulata) is often cited as supporting this hypothesis because SA selection has been demonstrated to act on male coloration in natural populations of this fish, and probably contributes to maintaining polymorphisms for the genetic factors involved. I review classical genetic and new molecular genetic results from the guppy, and other fish, including approaches for identifying the genome regions carrying sex-determining loci, and suggest that the guppy may exemplify a recently proposed route to sex chromosome evolution.

Restriction site-associated DNA sequencing for SNP discovery and high-density genetic map construction in southern catfish (Silurus meridionalis)

Single-nucleotide polymorphism (SNP) markers and high-density genetic maps are important resources for marker-assisted selection, mapping of quantitative trait loci (QTLs) and genome structure analysis. Although linkage maps in certain catfish species have been obtained, high-density maps remain unavailable in the economically important southern catfish (Silurus meridionalis). Recently developed restriction site-associated DNA (RAD) markers have proved to be a promising tool for SNP detection and genetic map construction. The objective of the present study was to construct a high-density linkage map using SNPs generated by next-generation RAD sequencing in S. meridionalis for future genetic and genomic studies. An F1 population of 100 individuals was obtained by intraspecific crossing of two wild heterozygous individuals. In total, 77 634 putative high-quality bi-allelic SNPs between the parents were discovered by mapping the parents' paired-end RAD reads onto the reference contigs from both parents, of which 54.7% were transitions and 45.3% were transversions (transition/transversion ratio of 1.2). Finally, 26 714 high-quality RAD markers were grouped into 29 linkage groups by using de novo clustering methods (Stacks). Among these markers, 4514 were linked to the female genetic map, 23 718 to the male map and 6715 effective loci were linked to the integrated map spanning 5918.31 centimorgans (cM), with an average marker interval of 0.89 cM. High-resolution genetic maps are a useful tool for both marker-assisted breeding and various genome investigations in catfish, such as sequence assembly, gene localization, QTL detection and genome structure comparison. Hence, such a high-density linkage map will serve as a valuable resource for comparative genomics and fine-scale QTL mapping in catfish species.

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.

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.

First haploid genetic map based on microsatellite markers in Senegalese sole (Solea senegalensis, Kaup 1858)

The Senegalese sole (Solea senegalensis, Kaup 1858) is a flatfish species of great value for aquaculture. In this study, we develop the first linkage map in this species based on microsatellite markers characterized from genomic DNA libraries and EST databases of Senegalese sole and from other flatfish species. Three reference gynogenetic families were obtained by chromosome-manipulation techniques: two haploid gynogenetics, used to assign and order microsatellites to linkage groups and another diploid gynogenetic family, used for estimating marker-centromere distances. The consensus map consists of 129 microsatellites distributed in 27 linkage groups (LG), with an average density of 4.7 markers per LG and comprising 1,004 centimorgans (cM). Additionally, 15 markers remained unlinked. Through half-tetrad analysis, we were able to estimate the centromere distance for 81 markers belonging to 24 LG, representing an average of 3 markers per LG. Comparative mapping was performed between flatfish species LG and model fish species chromosomes (stickleback, Tetraodon, medaka, fugu and zebrafish). The usefulness of microsatellite markers and the genetic map as tools for comparative mapping and evolution studies is discussed.

Comparative mapping for bighead carp (Aristichthys nobilis) against model and non-model fishes provides insights into the genomic evolution of cyprinids

Comparative mapping provides an efficient method to connect genomes of non-model and model fishes. In this study, we used flanking sequences of the 659 microsatellites on a genetic map of bighead carp (Aristichthys nobilis) to comprehensively study syntenic relationships between bighead carp and nine model and non-model fishes. Of the five model and two food fishes with whole genome data, Cyprinus carpio showed the highest rate of positive BLAST hits (95.3 %) with bighead carp map, followed by Danio rerio (70.9 %), Oreochromis niloticus (21.7 %), Tetraodon nigroviridis (6.4 %), Gasterosteus aculeatus (5.2 %), Oryzias latipes (4.7 %) and Fugu rubripes (3.5 %). Chromosomal syntenic analyses showed that inversion was the basic chromosomal rearrangement during genomic evolution of cyprinids, and the extent of inversions and translocations was found to be positively correlated with evolutionary relationships among fishes studied. Among the five investigated cyprinids, linkage groups (LGs) of bighead carp, Hypophthalmichthys molitrix and Ctenopharyngodon idella exhibited a one-to-one relationship. Besides, LG 9 of bighead carp and homologous LGs of silver carp and grass carp all corresponded to the chromosomes 10 and 22 of zebrafish, suggesting that chromosomal fission may have occurred in the ancestor of zebrafish. On the other hand, LGs of bighead carp and common carp showed an approximate one-to-two relationship with extensive translocations, confirming the occurrence of a 4th whole genome duplication in common carp. This study provides insights into the understanding of genome evolution among cyprinids and would aid in transferring positional and functional information of genes from model fish like zebrafish to non-model fish like bighead carp.

Development and application of biological technologies in fish genetic breeding

Fish genetic breeding is a process that remolds heritable traits to obtain neotype and improved varieties. For the purpose of genetic improvement, researchers can select for desirable genetic traits, integrate a suite of traits from different donors, or alter the innate genetic traits of a species. These improved varieties have, in many cases, facilitated the development of the aquaculture industry by lowering costs and increasing both quality and yield. In this review, we present the pertinent literatures and summarize the biological bases and application of selection breeding technologies (containing traditional selective breeding, molecular marker-assisted breeding, genome-wide selective breeding and breeding by controlling single-sex groups), integration breeding technologies (containing cross breeding, nuclear transplantation, germline stem cells and germ cells transplantation, artificial gynogenesis, artificial androgenesis and polyploid breeding) and modification breeding technologies (represented by transgenic breeding) in fish genetic breeding. Additionally, we discuss the progress our laboratory has made in the field of chromosomal ploidy breeding of fish, including distant hybridization, gynogenesis, and androgenesis. Finally, we systematically summarize the research status and known problems associated with each technology.

Novel sex-determining genes in fish and sex chromosome evolution

Although the molecular mechanisms underlying many developmental events are conserved across vertebrate taxa, the lability at the top of the sex-determining (SD) cascade has been evident from the fact that four master SD genes have been identified: mammalian Sry; chicken DMRT1; medaka Dmy; and Xenopus laevis DM-W. This diversity is thought to be associated with the turnover of sex chromosomes, which is likely to be more frequent in fishes and other poikilotherms than in therian mammals and birds. Recently, four novel candidates for vertebrate SD genes were reported, all of them in fishes. These include amhy in the Patagonian pejerrey, Gsdf in Oryzias luzonensis, Amhr2 in fugu and sdY in rainbow trout. These studies provide a good opportunity to infer patterns from the seemingly chaotic picture of sex determination systems. Here, we review recent advances in our understanding of the master SD genes in fishes.

A RAD-based linkage map and comparative genomics in the gudgeons (genus Gnathopogon, Cyprinidae)

Background

The construction of linkage maps is a first step in exploring the genetic basis for adaptive phenotypic divergence in closely related species by quantitative trait locus (QTL) analysis. Linkage maps are also useful for comparative genomics in non-model organisms. Advances in genomics technologies make it more feasible than ever to study the genetics of adaptation in natural populations. Restriction-site associated DNA (RAD) sequencing in next-generation sequencers facilitates the development of many genetic markers and genotyping. We aimed to construct a linkage map of the gudgeons of the genus Gnathopogon (Cyprinidae) for comparative genomics with the zebrafish Danio rerio (a member of the same family as gudgeons) and for the future QTL analysis of the genetic architecture underlying adaptive phenotypic evolution of Gnathopogon.

Results

We constructed the first genetic linkage map of Gnathopogon using a 198 F₂ interspecific cross between two closely related species in Japan: river-dwelling Gnathopogon elongatus and lake-dwelling Gnathopogon caerulescens. Based on 1,622 RAD-tag markers, a linkage map spanning 1,390.9 cM with 25 linkage groups and an average marker interval of 0.87 cM was constructed. We also identified a region involving female-specific transmission ratio distortion (TRD). Synteny and collinearity were extensively conserved between Gnathopogon and zebrafish.

Conclusions

The dense SNP-based linkage map presented here provides a basis for future QTL analysis. It will also be useful for transferring genomic information from a “traditional” model fish species, zebrafish, to screen candidate genes underlying ecologically important traits of the gudgeons.

Construction of a high-density microsatellite genetic linkage map and mapping of sexual and growth-related traits in half-smooth tongue sole (Cynoglossus semilaevis)

High-density genetic linkage maps of half-smooth tongue sole were developed with 1007 microsatellite markers, two SCAR markers and an F1 family containing 94. The female map was composed of 828 markers in 21 linkage groups, covering a total of 1447.3 cM, with an average interval 1.83 cM between markers. The male map consisted of 794 markers in 21 linkage groups, spanning 1497.5 cM, with an average interval of 1.96 cM. The female and male maps had 812 and 785 unique positions, respectively. The genome length of half-smooth tongue sole was estimated to be 1527.7 cM for the females and 1582.1 cM for the males. Based on estimations of the map lengths, the female and male maps covered 94.74 and 94.65% of the genome, respectively. The consensus map was composed of 1007 microsatellite markers and two SCAR markers in 21 linkage groups, covering a total of 1624 cM with an average interval of 1.67 cM. Furthermore, 159 sex-linked SSR markers were identified. Five sex-linked microsatellite markers were confirmed in their association with sex in a large number of individuals selected from different families. These sex-linked markers were mapped on the female map LG1f with zero recombination. Two QTLs that were identified for body weight, designated as We-1 and We-2, accounted for 26.39% and 10.60% of the phenotypic variation. Two QTLs for body width, designated Wi-1 and Wi-2, were mapped in LG4f and accounted for 14.33% and 12.83% of the phenotypic variation, respectively. Seven sex-related loci were mapped in LG1f, LG14f and LG1m by CIM, accounting for 12.5–25.2% of the trait variation. The results should prove to be very useful for improving growth traits using molecular MAS.

Construction of high-density genetic linkage maps and mapping of growth-related quantitative trail loci in the Japanese flounder (Paralichthys olivaceus)

High-density genetic linkage maps were constructed for the Japanese flounder (Paralichthys olivaceus). A total of 1624 microsatellite markers were polymorphic in the reference family. Linkage analysis using JoinMap 4.0 resulted in the mapping of 1487 markers to 24 linkage groups, a result which was consistent with the 24 chromosomes seen in chromosome spreads. The female map was composed of 1257 markers, covering a total of 1663.8 cM with an average interval 1.35 cM between markers. The male map consisted of 1224 markers, spanning 1726.5 cM, with an average interval of 1.44 cM. The genome length in the Japanese flounder was estimated to be 1730.3 cM for the females and 1798.0 cM for the males, a coverage of 96.2% for the female and 96.0% for the male map. The mean recombination at common intervals throughout the genome revealed a slight difference between sexes, i.e. 1.07 times higher in the male than female. High-density genetic linkage maps are very useful for marker-assisted selection (MAS) programs for economically valuable traits in this species and for further evolutionary studies in flatfish and vertebrate species. Furthermore, four quantiative trait loci (QTL) associated with growth traits were mapped on the genetic map. One QTL was identified for body weight on LG 14 f, which explained 14.85% of the total variation of the body weight. Three QTL were identified for body width on LG14f and LG14m, accounting for 16.75%, 13.62% and 13.65% of the total variation in body width, respectively. The additive effects were evident as negative values. There were four QTL for growth traits clustered on LG14, which should prove to be very useful for improving growth traits using molecular MAS.

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.

Diverse systems of sex determination have evolved independently in the animal and plant kingdoms. In vertebrates, so far four master sex-determining (SD) genes, Sry, Dmrt1, Dmy, and Dm-W, have been identified. These genes code for transcription factors and are located on only one of the sex chromosomes surrounded by nonrecombining regions. It is hypothesized that these sex chromosomes evolved from a pair of homologous chromosomes that diverged after acquiring the SD gene. We investigated the SD locus in fugu by high-resolution genetic mapping and association mapping. We found that a SNP that changes an amino acid (His/Asp384) in the kinase domain of anti-Müllerian hormone receptor type II (Amhr2) is perfectly associated with phenotypic sex. A combination of the two alleles of the SNP (homozygous females and heterozygous males) is likely to be responsible for sex determination in fugu. While these alleles are conserved in two other species of Takifugu, they are absent in the freshwater pufferfish, Tetraodon. Furthermore, Fugu Amhr2 lies in a region that shows no evidence for recombination suppression between X and Y chromosomes. Thus, fugu sex chromosomes represent an unusual example of a pre-differentiated phase of sex chromosomes in vertebrates.

An Expressed Sequence Tag (EST)-enriched genetic map of turbot (Scophthalmus maximus): a useful framework for comparative genomics across model and farmed teleosts

Background

The turbot (Scophthalmus maximus) is a relevant species in European aquaculture. The small turbot genome provides a source for genomics strategies to use in order to understand the genetic basis of productive traits, particularly those related to sex, growth and pathogen resistance. Genetic maps represent essential genomic screening tools allowing to localize quantitative trait loci (QTL) and to identify candidate genes through comparative mapping. This information is the backbone to develop marker-assisted selection (MAS) programs in aquaculture. Expressed sequenced tag (EST) resources have largely increased in turbot, thus supplying numerous type I markers suitable for extending the previous linkage map, which was mostly based on anonymous loci. The aim of this study was to construct a higher-resolution turbot genetic map using EST-linked markers, which will turn out to be useful for comparative mapping studies.

Results

A consensus gene-enriched genetic map of the turbot was constructed using 463 SNP and microsatellite markers in nine reference families. This map contains 438 markers, 180 EST-linked, clustered at 24 linkage groups. Linkage and comparative genomics evidences suggested additional linkage group fusions toward the consolidation of turbot map according to karyotype information. The linkage map showed a total length of 1402.7 cM with low average intermarker distance (3.7 cM; ~2 Mb). A global 1.6:1 female-to-male recombination frequency (RF) ratio was observed, although largely variable among linkage groups and chromosome regions. Comparative sequence analysis revealed large macrosyntenic patterns against model teleost genomes, significant hits decreasing from stickleback (54%) to zebrafish (20%). Comparative mapping supported particular chromosome rearrangements within Acanthopterygii and aided to assign unallocated markers to specific turbot linkage groups.

Conclusions

The new gene-enriched high-resolution turbot map represents a useful genomic tool for QTL identification, positional cloning strategies, and future genome assembling. This map showed large synteny conservation against model teleost genomes. Comparative genomics and data mining from landmarks will provide straightforward access to candidate genes, which will be the basis for genetic breeding programs and evolutionary studies in this species.

A lineage specific enhancer drives Otx2 expression in teleost organizer tissues

In mouse Otx2 plays essential roles in anterior-posterior axis formation and head development in anterior visceral endoderm and anterior mesendoderm. The Otx2 expression in these sites is regulated by VE and CM enhancers at the 5' proximal to the translation start site, and we proposed that these enhancers would have been established in ancestral sarcoptergians after divergence from actinopterigians for the use of Otx2 as the head organizer gene (Kurokawa et al., 2010). This would make doubtful an earlier proposal of ours that a 1.1 kb fragment located at +14.4 to +15.5 kb 3' (3'En) of fugu Otx2a gene harbors enhancers phylogenetically and functionally homologous to mouse VE and CM enhancers (Kimura-Yoshida et al., 2007). In the present study, we demonstrate that fugu Otx2a is not expressed in the dorsal margin of blastoderm, shield and early anterior mesendoderm, and that the fugu Otx2a 3'En do not exhibit activities at these sites of fugu embryos. We conclude that the fugu Otx2a 3'En does not harbor an organizer enhancer, but encodes an enhancer for the expression in later anterior mesendodermal tissues. Instead, in fugu embryos Otx2b is expressed in the dorsal margin of blastoderm at blastula stage and shield at 50% epiboly, and this expression is directed by an enhancer, 5'En, located at -1000 to -800 bp, which is uniquely conserved among teleost Otx2b orthologues.

Correlation with larval body size of mRNA levels of growth hormone, growth hormone receptor I and insulin-like growth factor I in larval torafugu Takifugu rubripes

The full-length of insulin-like growth factor (IGF) complementary (c)DNAs encoded by igf-I and igf-II from torafugu pufferfish Takifugu rubripes were cloned in the present study. The deduced amino acid sequences of the two genes showed c. 80% identity each with those of Igf-I and Igf-II from other teleosts, respectively. Two growth hormone (GH) receptors, ghr1 and ghr2, were also cloned in silico using the T. rubripes Fugu genome database. The transcripts of T. rubripes igf-I were detected in slow muscle, heart, skin, gill, liver and intestine but not in fast muscle, spleen and testis of adult fish, whereas those of igf-II were found in all tissues examined. Subsequently, the accumulated messenger (m)RNA levels of igf-I and igf-II were investigated in an F(2) population derived from a male of an apparent fast-growing T. rubripes strain and a wild female T. rubripes together with those of other growth-related genes encoding Gh, Ghr1 and Ghr2, and with those of prolactin (Prl) and leptin (Lep) previously reported. The accumulated mRNA levels of igf-I, gh and ghr1 were significantly correlated to growth rate at larval stages in the population, but not for those of igf-II, prl, ghr2 and lep. Although it is unclear whether or not this phenotype is directly related to the heredity of the fast-growing strain, the findings suggest that the expression of igf-I, gh and ghr1 is involved in the regulation of growth rate at larval stages in T. rubripes.

Genetic diversity of microsatellite loci in hierarchically structured populations

Microsatellite loci are widely used for investigating patterns of genetic variation within and among populations. Those patterns are in turn determined by population sizes, migration rates, and mutation rates. We provide exact expressions for the first two moments of the allele frequency distribution in a stochastic model appropriate for studying microsatellite evolution with migration, mutation, and drift under the assumption that the range of allele sizes is bounded. Using these results, we study the behavior of several measures related to Wright's F(ST), including Slatkin's R(ST). Our analytical approximations for F(ST) and R(ST) show that familiar relationships between N(e)m and F(ST) or R(ST) hold when the migration and mutation rates are small. Using the exact expressions for F(ST) and R(ST), our numerical results show that, when the migration and mutation rates are large, these relationships no longer hold. Our numerical results also show that the diversity measures most closely related to F(ST) depend on mutation rates, mutational models (stepwise versus two-phase), migration rates, and population sizes. Surprisingly, R(ST) is relatively insensitive to the mutation rates and mutational models. The differing behaviors of R(ST) and F(ST) suggest that properties of the among-population distribution of allele frequencies may allow the roles of mutation and migration in producing patterns of diversity to be distinguished, a topic of continuing investigation.

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.

Effect of cadmium exposure on expression of antioxidant gene transcripts in the river pufferfish, Takifugu obscurus (Tetraodontiformes)

Cadmium (Cd) is a non-essential toxic heavy metal with the potential to induce oxidative stress. Cd toxicity and its capacity for accumulation in aquatic habitats have earned its recognition as a pollutant of immediate and widespread concern. To obtain a better understanding of oxidative stress-associated gene expression in different tissues, six antioxidant genes such as catalase (CAT), glutathione reductase (GR), glutathione peroxidase 1a (GPx1a), glutathione peroxidase 1b (GPx1b), Cu/Zn superoxide dismutase (Cu/Zn-SOD), and Mn superoxide dismutase (Mn-SOD) were cloned and fully sequenced in the river pufferfish, Takifugu obscurus. On tissue specific mRNA expression, the liver showed the highest expression when compared to other tissues, even though each antioxidant gene showed different modes of expression patterns in the examined tissues. Of the various antioxidant genes, GR was the most highly expressed in the liver, followed by CAT, GPx1, and Cu/Zn-SOD. For the time-course experiment, all the antioxidant genes were significantly induced over time except for Cu/Zn-SOD in the liver, and there was a 5-fold induction in hepatic GR, CAT, and Mn-SOD mRNA compared to the control. These findings indicate that the liver of T. obscurus has a robust antioxidant system. In addition, these results suggest that Cd exposure modulates the expression of antioxidant genes, and would indicate that the antioxidant genes would be a relevant biomarker of trace metal pollution such as Cd exposure in T. obscurus.

A consensus linkage map of the grass carp (Ctenopharyngodon idella) based on microsatellites and SNPs

Background

Grass carp (Ctenopharyngodon idella) belongs to the family Cyprinidae which includes more than 2000 fish species. It is one of the most important freshwater food fish species in world aquaculture. A linkage map is an essential framework for mapping traits of interest and is often the first step towards understanding genome evolution. The aim of this study is to construct a first generation genetic map of grass carp using microsatellites and SNPs to generate a new resource for mapping QTL for economically important traits and to conduct a comparative mapping analysis to shed new insights into the evolution of fish genomes.

Results

We constructed a first generation linkage map of grass carp with a mapping panel containing two F1 families including 192 progenies. Sixteen SNPs in genes and 263 microsatellite markers were mapped to twenty-four linkage groups (LGs). The number of LGs was corresponding to the haploid chromosome number of grass carp. The sex-specific map was 1149.4 and 888.8 cM long in females and males respectively whereas the sex-averaged map spanned 1176.1 cM. The average resolution of the map was 4.2 cM/locus. BLAST searches of sequences of mapped markers of grass carp against the whole genome sequence of zebrafish revealed substantial macrosynteny relationship and extensive colinearity of markers between grass carp and zebrafish.

Conclusions

The linkage map of grass carp presented here is the first linkage map of a food fish species based on co-dominant markers in the family Cyprinidae. This map provides a valuable resource for mapping phenotypic variations and serves as a reference to approach comparative genomics and understand the evolution of fish genomes and could be complementary to grass carp genome sequencing project.

Teleost B7 expressed on monocytes regulates T cell responses

In mammals, professional APCs induce adaptive immunity via the activation of T cells. During this process, B7 family molecules present upon APCs are known to play crucial roles in optimal T cell stimulation. In contrast, the confirmation of APCs in a nonmammalian vertebrate has yet to be achieved. To obtain further insights into the evolutionary origin of APCs, we have identified three members of the B7 family in the teleost Takifugu rubripes (fugu): B7-H1/DC, B7-H3, and B7-H4. The three fugu B7s were expressed on the surface of blood monocytes. The B7(+) monocytes, which are composed of at least two distinct populations, expressed the MHC class II component gene. The fugu B7 molecules bound to activated T cells, indicating that putative B7 receptors were expressed upon T cells. Fugu B7-H1/DC inhibited T cell proliferation concomitant with increasing levels of both IL-10 and IFN-gamma expression, whereas both B7-H3 and B7-H4 promoted T cell growth following IL-2 induction and the suppression of IL-10. These observations indicate that fugu B7s regulate T cell responses via receptors upon T cells. We suggest that fish B7(+) monocytes are APCs and that a costimulatory system has already developed in fish via the evolutionary process.

Genetic linkage map of the guppy, Poecilia reticulata, and quantitative trait loci analysis of male size and colour variation

We report construction of a genetic linkage map of the guppy genome using 790 single nucleotide polymorphism markers, integrated from six mapping crosses. The markers define 23 linkage groups (LGs), corresponding to the known haploid number of guppy chromosomes. The map, which spans a genetic length of 899 cM, includes 276 markers linked to expressed genes (expressed sequence tag), which have been used to derive broad syntenic relationships of guppy LGs with medaka chromosomes. This combined linkage map should facilitate the advancement of genetic studies for a wide variety of complex adaptive phenotypes relevant to natural and sexual selection in this species. We have used the linkage data to predict quantitative trait loci for a set of variable male traits including size and colour pattern. Contributing loci map to the sex LG for many of these traits.

A first-generation microsatellite-based genetic linkage map of the Siberian jay (Perisoreus infaustus): insights into avian genome evolution

BACKGROUND

Genomic resources for the majority of free-living vertebrates of ecological and evolutionary importance are scarce. Therefore, linkage maps with high-density genome coverage are needed for progress in genomics of wild species. The Siberian jay (Perisoreus infaustus; Corvidae) is a passerine bird which has been subject to lots of research in the areas of ecology and evolutionary biology. Knowledge of its genome structure and organization is required to advance our understanding of the genetic basis of ecologically important traits in this species, as well as to provide insights into avian genome evolution.

RESULTS

We describe the first genetic linkage map of Siberian jay constructed using 117 microsatellites and a mapping pedigree of 349 animals representing five families from a natural population breeding in western Finland from the years 1975 to 2006. Markers were resolved into nine autosomal and a Z-chromosome-specific linkage group, 10 markers remaining unlinked. The best-position map with the most likely positions of all significantly linked loci had a total sex-average size of 862.8 cM, with an average interval distance of 9.69 cM. The female map covered 988.4 cM, whereas the male map covered only 774 cM. The Z-chromosome linkage group comprised six markers, three pseudoautosomal and three sex-specific loci, and spanned 10.6 cM in females and 48.9 cM in males. Eighty-one of the mapped loci could be ordered on a framework map with odds of >1000:1 covering a total size of 809.6 cM in females and 694.2 cM in males. Significant sex specific distortions towards reduced male recombination rates were revealed in the entire best-position map as well as within two autosomal linkage groups. Comparative mapping between Siberian jay and chicken anchored 22 homologous loci on 6 different linkage groups corresponding to chicken chromosomes Gga1, 2, 3, 4, 5, and Z. Quite a few cases of intra-chromosomal rearrangements within the autosomes and three cases of inter-chromosomal rearrangement between the Siberian jay autosomal linkage groups (LG1, LG2 and LG3) and the chicken sex chromosome GgaZ were observed, suggesting a conserved synteny, but changes in marker order, within autosomes during about 100 million years of avian evolution.

CONCLUSION

The constructed linkage map represents a valuable resource for intraspecific genomics of Siberian jay, as well as for avian comparative genomic studies. Apart from providing novel insights into sex-specific recombination rates and patterns, the described maps - from a previously genomically uncharacterized superfamily (Corvidae) of passerine birds - provide new insights into avian genome evolution. In combination with high-resolution data on quantitative trait variability from the study population, they also provide a foundation for QTL-mapping studies.

Explosive speciation of Takifugu: another use of fugu as a model system for evolutionary biology

Although the fugu Takifugu rubripes has attracted attention as a model organism for genomic studies because of its compact genome, it is not generally appreciated that there are approximately 25 closely related species with limited distributions in the waters of East Asia. We performed molecular phylogenetic analyses and constructed a time tree using whole mitochondrial genome sequences from 15 Takifugu species together with 10 outgroups to examine patterns of diversification. The resultant time tree showed that the modern Takifugu species underwent explosive speciation during the Pliocene 1.8-5.3 Ma, which is comparable with that of the Malawi cichlids and tropheine cichlids in Lake Tanganyika. Considering their limited distributions and remarkable variations in coloration, morphology, and behavior, the results of the present study strongly suggest that Takifugu species are strong candidates as a model system for evolutionary studies of speciation mechanisms in marine environments where few such organisms are available.

Evolution of Otx paralogue usages in early patterning of the vertebrate head

To assess evolutional changes in the expression pattern of Otx paralogues, expression analyses were undertaken in fugu, bichir, skate and lamprey. Together with those in model vertebrates, the comparison suggested that a gnathostome ancestor would have utilized all of Otx1, Otx2 and Otx5 paralogues in organizer and anterior mesendoderm for head development. In this animal, Otx1 and Otx2 would have also functioned in specification of the anterior neuroectoderm at presomite stage and subsequent development of forebrain/midbrain at somite stage, while Otx5 expression would have already been specialized in epiphysis and eyes. Otx1 and Otx2 functions in anterior neuroectoderm and brain of the gnathostome ancestor would have been differentially maintained by Otx1 in a basal actinopterygian and by Otx2 in a basal sarcopterygian. Otx5 expression in head organizer and anterior mesendoderm seems to have been lost in the teleost lineage after divergence of bichir, and also from the amniotes after divergence of amphibians as independent events. Otx1 expression was lost from the organizer in the tetrapod lineage. In contrast, in a teleost ancestor prior to whole genome duplication, Otx1 and Otx2 would have both been expressed in the dorsal margin of blastoderm, embryonic shield, anterior mesendoderm, anterior neuroectoderm and forebrain/midbrain, at respective stages of head development. Subsequent whole genome duplication and the following genome changes would have caused different Otx paralogue usages in each teleost lineage. Lampreys also have three Otx paralogues; their sequences are highly diverged from gnathostome cognates, but their expression pattern is well related to those of skate Otx cognates.

A combined AFLP and microsatellite linkage map and pilot comparative genomic analysis of European sea bass Dicentrarchus labrax L

European sea bass (Dicentrarchus labrax L., Moronidae, Teleostei) sustains a regional fishery and is commonly farmed in the Mediterranean basin, but has not undergone much long-term genetic improvement. An updated genetic linkage map of the European sea bass was constructed using 190 microsatellites, 176 amplified fragment length polymorphisms and two single nucleotide polymorphisms. From the 45 new microsatellite markers (including 31 type I markers) reported in this study, 28 were mapped. A total of 368 markers were assembled into 35 linkage groups. Among these markers, 28 represented type I (coding) markers, including those located within the peptide Y, SOX10, PXN1, ERA and TCRB genes (linkage groups 1, 7, 16, 17 and 27 respectively). The sex-averaged map spanned 1373.1 centimorgans (cM) of the genome. The female map measured 1380.0 cM, whereas the male map measured 1046.9 cM, leading to a female-to-male (F:M) recombination rate ratio of 1.32:1. The intermarker spacing of the second-generation linkage map of the European sea bass was 3.67 cM, which is smaller than that of the first-generation linkage map (5.03 cM). Comparative mapping of microsatellite flanking regions was performed with five model teleosts and this revealed a high percentage (33.6%) of evolutionarily conserved regions with the three-spined stickleback.

SynBlast: assisting the analysis of conserved synteny information

Motivation

In the last years more than 20 vertebrate genomes have been sequenced, and the rate at which genomic DNA information becomes available is rapidly accelerating. Gene duplication and gene loss events inherently limit the accuracy of orthology detection based on sequence similarity alone. Fully automated methods for orthology annotation do exist but often fail to identify individual members in cases of large gene families, or to distinguish missing data from traceable gene losses. This situation can be improved in many cases by including conserved synteny information.

Results

Here we present the SynBlast pipeline that is designed to construct and evaluate local synteny information. SynBlast uses the genomic region around a focal reference gene to retrieve candidates for homologous regions from a collection of target genomes and ranks them in accord with the available evidence for homology. The pipeline is intended as a tool to aid high quality manual annotation in particular in those cases where automatic procedures fail. We demonstrate how SynBlast is applied to retrieving orthologous and paralogous clusters using the vertebrate Hox and ParaHox clusters as examples.

Software

The SynBlast package written in Perl is available under the GNU General Public License at .

A microsatellite genetic map of the turbot (Scophthalmus maximus)

A consensus microsatellite-based linkage map of the turbot (Scophthalmus maximus) was constructed from two unrelated families. The mapping panel was derived from a gynogenetic family of 96 haploid embryos and a biparental diploid family of 85 full-sib progeny with known linkage phase. A total of 242 microsatellites were mapped in 26 linkage groups, six markers remaining unlinked. The consensus map length was 1343.2 cM, with an average distance between markers of 6.5 +/- 0.5 cM. Similar length of female and male maps was evidenced. However, the mean recombination at common intervals throughout the genome revealed significant differences between sexes, approximately 1.6 times higher in the female than in the male. The comparison of turbot microsatellite flanking sequences against the Tetraodon nigroviridis genome revealed 55 significant matches, with a mean length of 102 bp and high sequence similarity (81-100%). The comparative mapping revealed significant syntenic regions among fish species. This study represents the first linkage map in the turbot, one of the most important flatfish in European aquaculture. This map will be suitable for QTL identification of productive traits in this species and for further evolutionary studies in fish and vertebrate species.

Description of a fugu CXC chemokine and two CXC receptor genes, and characterization of the effects of different stimulators on their expression

The primary structures of a CXC chemokine (CXCL8) and two CXC receptors (CXCR) have been characterized in fugu, Takifugu rubripes. Unlike mammalian and avian species, CXCL8 of teleosts including fugu lacks the ELR motif that appears to be important in ligand/receptor interactions on neutrophils. Genomic organization shows that fugu CXCL8 gene consists of four exons and three introns. As in other vertebrates, two CXCR genes isolated from fugu encode proteins CXCR1 and CXCR2 that possess characteristic seven transmembrane domains. Each receptor consists of two exons separated by an intron. Synteny analysis indicates that these two CXCRs were derived from whole genome duplication in teleosts, differing from mammalian CXCR1 and CXCR2. All of these genes are primarily expressed in the lymphoid tissues. Immune stimulation with PHA showed that the expression of both CXCL8 and CXCRs in PBL are upregulated even after only a short time period, but downregulated by LPS stimulation, implying that these genes are involved in the regulation of the immune response in fugu.

Development of silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) genetic maps using microsatellite and AFLP markers and a pseudo-testcross strategy

Silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) are two of the four most important pond-cultured fish species inhabiting the major river basins of China. In the present study, genetic maps of silver carp and bighead carp were constructed using microsatellite and AFLP markers and a two-way pseudo-testcross strategy. To create the maps, 60 individuals were obtained from a cross of a single bighead carp (female) and a single silver carp (male). The silver carp map consisted of 271 markers (48 microsatellites and 223 AFLPs) that were assembled into 27 linkage groups, of which 22 contained at least four markers. The total length of the silver carp map was 952.2 cM, covering 82.8% of the estimated genome size. The bighead carp map consisted of 153 markers (27 microsatellites and 126 AFLPs) which were organized into 30 linkage groups, of which 19 contained at least four markers. The total length of the bighead carp map was 852.0 cM, covering 70.5% of the estimated genome size. Eighteen microsatellite markers were common to both maps. These maps will contribute to discovery of genes and genetic regions controlling traits in the two species of carp.

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.