The genetic architecture of growth rate in juvenile Takifugu species

Host Hybridization Dominates over Cohabitation in Affecting Gut Microbiota of Intrageneric Hybrid Takifugu Pufferfish

Microbial symbionts are of great importance for macroscopic life, including fish, and both collectively comprise an integrated biological entity known as the holobiont. Yet little is known as to how the normal balance within the fish holobiont is maintained and how it responds to biotic and/or abiotic influences. Here, through amplicon profiling, the genealogical relationship between artificial F1 hybrid pufferfish with growth heterosis, produced from crossing female Takifugu obscurus with male Takifugu rubripes and its maternal halfsibling purebred, was well recapitulated by their gut microbial community similarities, indicating an evident parallelism between host phylogeny (hybridity) and microbiota relationships therein. Interestingly, modest yet significant fish growth promotion and gut microbiota alteration mediated by hybrid-purebred cohabitation were observed, in comparison with their respective monoculture cohorts that share common genetic makeups, implying a certain degree of environmental influences. Moreover, the underlying assemblage patterns of gut microbial communities were found associated with a trade-off between variable selection and dispersal limitation, which are plausibly driven by the augmented social interactions between hybrid and purebred cohabitants differing in behaviors. Results from this study not only can enrich, from a microbial perspective, the sophisticated understanding of complex and dynamic assemblage of the fish holobiont, but will also provide deeper insights into the ecophysiological factors imposed on the diversity-function relationships thereof. Our findings emphasize the intimate associations of gut microbiota in host genetics-environmental interactions and would have deeper practical implications for microbial contributions to optimize performance prediction and to improve the production of farmed fishes. IMPORTANCE Microbial symbionts are of great importance for macroscopic life, including fish, and yet little is known as to how the normal balance within the fish holobiont is maintained and how it responds to the biotic and/or abiotic influences. Through gut microbiota profiling, we show that host intrageneric hybridization and cohabitation can impose a strong disturbance upon pufferfish gut microbiota. Moreover, marked alterations in the composition and function of gut microbiota in both hybrid and purebred pufferfish cohabitants were observed, which are potentially correlated with different metabolic priorities and behaviors between host genealogy. These results can enrich, from a microbial perspective, the sophisticated understanding of the complex and dynamic assemblage of the fish holobiont and would have deeper practical implications for microbial contributions to optimize performance prediction and to improve farmed fish production.

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

SignificanceAlthough turnover of sex chromosomes is very common in many vertebrate lineages, the transition process is still elusive. We studied the sex-determining region (SDR) of 12 congeneric fish species of Takifugu fish that compose an important model for the study of genomics and sex determination. We found that while nine species retained their ancestral SDR, three species had acquired derived SDRs. Although the derived SDRs resided in three different chromosomes, they harbored a shared supergene flanked by two putative transposable elements. The results highlight the underestimated role of a mobile supergene in turnover of sex chromosomes in vertebrates.

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.

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.

Capacity for freshwater acclimation and differences in the transcription of ion transporter genes underlying different migratory life histories of Takifugu fish

The genus Takifugu is a group of approximately 20 species of puffer fishes living in a wide range of salinity environments around East Asian countries. This group presents a broad spectrum of evolutionary stages adapted to anadromy as a result of speciation that occurred a short time (2-5 million years) ago on an evolutionary timescale. This group thus can be considered as a model for studying the evolutionary mechanisms of anadromy. We firstly conducted a transfer experiment from seawater to low-salinity waters on five Takifugu species: two anadromous species T. obscurus and T. ocellatus, two euryhaline wanderer marine species T. rubripes and T. niphobles, and a strictly marine species T. snyderi, and confirmed that the capacity for acclimation to hypotonic environments was associated with their life history strategies. Next, transcriptomes of the gill and intestine of these species in hypotonic condition were compared to those under hypertonic condition for each species using RNA-Sequencing so as to determine possible candidate transporters playing an important role on freshwater adaptation. As this analysis suggested that cftr, encoding an important ion transporter for seawater acclimation in the gill, and ncc, encoding a transporter that is suggested to play important osmoregulatory roles in the intestine, are important candidates, their expression was validated by quantitative real-time PCR analysis. Expression of cftr was downregulated in the gills of the four euryhaline species under the hypotonic condition, but no change was detected in the gill of stenohaline T. snyderi, which may be one reason for the poor hypotonic acclimation capacity of T. snyderi. Expression of ncc was clearly upregulated in the intestines of the two anadromous species under the hypotonic condition, but not in other three species. Different ion transporter expression patterns between the five species indicate that the transcriptional regulation of cftr in the gill and ncc in the intestine may be important for the improvement of hypotonic acclimation capacity and evolution of anadromy in the Takifugu species.

D-mannose-specific Immunoglobulin M in Grass Puffer (Takifugu niphobles), a Nonhost Fish of a Monogenean Ectoparasite Heterobothrium okamotoi, Can Act as a Trigger for its Parasitism

Heterobothrium okamotoi, a monogenean gill parasite, exhibits high host specificity for the tiger puffer, Takifugu rubripes, and it has been experimentally verified that the parasite cannot colonize either closely related species such as the grass puffer Takifugu niphobles or distantly related fish such as the red seabream Pagrus major. Previously, we demonstrated in T. rubripes that immunoglobulin M (IgM) with d-mannose affinity induced deciliation of the oncomiracidia, the first step of parasitism, indicating that the parasite utilizes the molecule as a receptor for infection. In the present study, we purified mannose-specific IgM from 2 nonhost species, T. niphobles and P. major, by affinity and gel-filtration chromatography techniques and compared their deciliation-inducing activity against H. okamotoi oncomiracidia. The IgM of the former showed activity, whereas the latter had no effect, suggesting that in addition to d-mannose-binding ability, the crystallizable fragment domain of IgM, which is not part of the antigen-binding domain, plays an important role in host recognition by the oncomiracidia, such as direct binding to the parasites. It also suggests that the host specificity of H. okamotoi is relatively low upon initial recognition, and the specificity is established by exclusion in nonhosts during a later stage.

Genetic analysis of body weight in wild populations of medaka fish from different latitudes

The genetic bases of growth and body weight are of economic and scientific interest, and teleost fish models have proven useful in such investigations. The Oryzias latipes species complex (medaka) is an abundant freshwater fish in Japan and suitable for genetic studies. We compared two wild medaka stocks originating from different latitudes. The Maizuru population from higher latitudes weighed more than the Ginoza population. We investigated the genetic basis of body weight, using quantitative trait locus (QTL) analysis of the F2 offspring of these populations. We detected one statistically significant QTL for body weight on medaka chromosome 4 and identified 12 candidate genes that might be associated with body weight or growth. Nine of these 12 genes had at least one single nucleotide polymorphism that caused amino acid substitutions in protein-coding regions, and we estimated the effects of these substitutions. The present findings might contribute to the marker-assisted selection of economically important aquaculture species.

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.

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.

Identification of differentially expressed genes associated with differential body size in mandarin fish (Siniperca chuatsi)

Body size is an obvious and important characteristic of fish. Mandarin fish Siniperca chuatsi (Basilewsky) is one of the most valuable perciform species widely cultured in China. Individual differences in body size are common in mandarin fish and significantly influence the aquaculture production. However, little is currently known about its genetic control. In this study, digital gene expression profiling and transcriptome sequencing were performed in mandarin fish with differential body size at 30 and 180 days post-hatch (dph), respectively. Body weight, total length and body length of fish with big-size were significantly higher than those with small-size at both 30 and 180 dph (P < 0.05). 2171 and 2014 differentially expressed genes were identified between small-size and big-size fish at 30 and 180 dph, respectively. RT quantitative PCR (qPCR) analysis showed that the differential expression of 10 selected genes in mandarin fish that went through the same training procedure. The genes were involved in the growth hormone-insulin-like growth factor axis, cell proliferation and differentiation, appetite control, glucose metabolism, reproduction and sexual size dimorphism pathways. This study will help toward a comprehensive understanding of the complexity of regulation of body size in mandarin fish individuals and provide valuable information for future research.

A high-density genetic map and growth related QTL mapping in bighead carp (Hypophthalmichthys nobilis)

Growth related traits in fish are controlled by quantitative trait loci (QTL), but no QTL for growth have been detected in bighead carp (Hypophthalmichthys nobilis) due to the lack of high-density genetic map. In this study, an ultra-high density genetic map was constructed with 3,121 SNP markers by sequencing 117 individuals in a F₁ family using 2b-RAD technology. The total length of the map was 2341.27 cM, with an average marker interval of 0.75 cM. A high level of genomic synteny between our map and zebrafish was detected. Based on this genetic map, one genome-wide significant and 37 suggestive QTL for five growth-related traits were identified in 6 linkage groups (i.e. LG3, LG11, LG15, LG18, LG19, LG22). The phenotypic variance explained (PVE) by these QTL varied from 15.4% to 38.2%. Marker within the significant QTL region was surrounded by CRP1 and CRP2, which played an important role in muscle cell division. These high-density map and QTL information provided a solid base for QTL fine mapping and comparative genomics in bighead carp.

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.