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.

Genome-wide DNA Methylation and Gene Expression Patterns of Androgenetic Haploid Tiger Pufferfish (Takifugu Rubripes) Provide Insights into Haploid Syndrome

Androgenesis is an important chromosome set manipulation technique used in sex control in aquaculture. Haploid embryos exhibit haploid syndrome and show abnormalities and even die during early embryonic development. In this study, we used whole genome bisulfite sequencing (WGBS) to investigate the genome-wide DNA methylation profiles in haploid females (1n-X) and males (1n-Y), and diploid females (2n-XX) and males (2n-XY) of tiger pufferfish (Takifugu rubripes), an economically important fish in China. A total of 96.32 Gb clean data was produced. Analysis of differentially methylated regions (DMRs) showed that haploids had more hyper-methylated regions than diploids, which may be related to abnormal development and early embryonic death in haploids. There were 7,838 hyper-methylated differentially methylated genes (DMGs) and 4,755 hypo-methylated DMGs in haploid vs. diploid comparisons in both females and males. These DMGs were mainly related to genomic stability maintenance and cell cycle regulation. slf1, actr8, gas2, and pbrm1 genes were detected to validate the methylation sequencing. After combining the methylation data with the corresponding transcriptome data, we identified several genes, including guca2a, myoc, fezf2, rprml, telo2, s100a1, and marveld1, which exhibited differential expression levels modulated by DNA methylation. In conclusion, our study revealed different methylation and expression profiles between haploid and diploid T. rubripes for the first time. Several DMGs were identified between different ploidy levels, which may be related to haploid syndrome formation. The results expand the understanding of the effects of ploidy on the early development of teleosts and provide knowledge about target genes and networks to improve the survival rate of haploids.

Identification of sex differentiation-related microRNA and long non-coding RNA in Takifugu rubripes gonads

Although sex determination and differentiation are key developmental processes in animals, the involvement of non-coding RNA in the regulation of this process is still not clarified. The tiger pufferfish (Takifugu rubripes) is one of the most economically important marine cultured species in Asia, but analyses of miRNA and long non-coding RNA (lncRNA) at early sex differentiation stages have not been conducted yet. In our study, high-throughput sequencing technology was used to sequence transcriptome libraries from undifferentiated gonads of T. rubripes. In total, 231 (107 conserved, and 124 novel) miRNAs were obtained, while 2774 (523 conserved, and 2251 novel) lncRNAs were identified. Of these, several miRNAs and lncRNAs were predicted to be the regulators of the expression of sex-related genes (including fru-miR-15b/foxl2, novel-167, novel-318, and novel-538/dmrt1, novel-548/amh, lnc_000338, lnc_000690, lnc_000370, XLOC_021951, and XR_965485.1/gsdf). Analysis of differentially expressed miRNAs and lncRNAs showed that three mature miRNAs up-regulated and five mature miRNAs were down-regulated in male gonads compared to female gonads, while 79 lncRNAs were up-regulated and 51 were down-regulated. These findings could highlight a group of interesting miRNAs and lncRNAs for future studies and may reveal new insights into the function of miRNAs and lncRNAs in sex determination and differentiation.

Correction to: Precocious maturation in male tiger pufferfish Takifugu rubripes: genetics and endocrinology

A correction to this paper has been published: https://doi.org/10.1007/s12562-021-01505-w

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.

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.