A standard panel of microsatellites for Asian seabass (Lates calcarifer)

Mapping of a major QTL for increased robustness and detection of genome assembly errors in Asian seabass (Lates calcarifer)

BACKGROUND

For Asian seabass (Lates calcarifer, Bloch 1790) cultured at sea cages various aquatic pathogens, complex environmental and stress factors are considered as leading causes of disease, causing tens of millions of dollars of annual economic losses. Over the years, we conducted farm-based challenges by exposing Asian seabass juveniles to complex natural environmental conditions. In one of these challenges, we collected a total of 1,250 fish classified as either 'sensitive' or 'robust' individuals during the 28-day observation period.

RESULTS

We constructed a high-resolution linkage map with 3,089 SNPs for Asian seabass using the double digest Restriction-site Associated DNA (ddRAD) technology and a performed a search for Quantitative Trait Loci (QTL) associated with robustness. The search detected a major genome-wide significant QTL for increased robustness in pathogen-infected marine environment on linkage group 11 (ASB_LG11; 88.9 cM to 93.6 cM) with phenotypic variation explained of 81.0%. The QTL was positioned within a > 800 kb genomic region located at the tip of chromosome ASB_LG11 with two Single Nucleotide Polymorphism markers, R1-38468 and R1-61252, located near to the two ends of the QTL. When the R1-61252 marker was validated experimentally in a different mass cross population, it showed a statistically significant association with increased robustness. The majority of thirty-six potential candidate genes located within the QTL have known functions related to innate immunity, stress response or disease. By utilizing this ddRAD-based map, we detected five mis-assemblies corresponding to four chromosomes, namely ASB_LG8, ASB_LG9, ASB_LG15 and ASB_LG20, in the current Asian seabass reference genome assembly.

CONCLUSION

According to our knowledge, the QTL associated with increased robustness is the first such finding from a tropical fish species. Depending on further validation in other stocks and populations, it might be potentially useful for selecting robust Asian seabass lines in selection programs.

Toward Genome-Based Selection in Asian Seabass: What Can We Learn From Other Food Fishes and Farm Animals?

Due to the steadily increasing need for seafood and the plateauing output of fisheries, more fish need to be produced by aquaculture production. In parallel with the improvement of farming methods, elite food fish lines with superior traits for production must be generated by selection programs that utilize cutting-edge tools of genomics. The purpose of this review is to provide a historical overview and status report of a selection program performed on a catadromous predator, the Asian seabass (Lates calcarifer, Bloch 1790) that can change its sex during its lifetime. We describe the practices of wet lab, farm and lab in detail by focusing onto the foundations and achievements of the program. In addition to the approaches used for selection, our review also provides an inventory of genetic/genomic platforms and technologies developed to (i) provide current and future support for the selection process; and (ii) improve our understanding of the biology of the species. Approaches used for the improvement of terrestrial farm animals are used as examples and references, as those processes are far ahead of the ones used in aquaculture and thus they might help those working on fish to select the best possible options and avoid potential pitfalls.

BAC-pool sequencing and analysis confirms growth-associated QTLs in the Asian seabass genome

The Asian seabass is an important marine food fish that has been cultured for several decades in Asia Pacific. However, the lack of a high quality reference genome has hampered efforts to improve its selective breeding. A 3D BAC pool set generated in this study was screened using 22 SSR markers located on linkage group 2 which contains a growth-related QTL region. Seventy-two clones corresponding to 22 FPC contigs were sequenced by Illumina MiSeq technology. We co-assembled the MiSeq-derived scaffolds from each FPC contig with error-corrected PacBio reads, resulting in 187 sequences covering 9.7 Mb. Eleven genes annotated within this region were found to be potentially associated with growth and their tissue-specific expression was investigated. Correlation analysis demonstrated that SNPs in ctsb, skp1 and ppp2ca can be potentially used as markers for selecting fast-growing fingerlings. Conserved syntenies between seabass LG2 and five other teleosts were identified. This study i) provided a 10 Mb targeted genome assembly; ii) demonstrated NGS of BAC pools as a potential approach for mining candidates underlying QTLs of this species; iii) detected eleven genes potentially responsible for growth in the QTL region; and iv) identified useful SNP markers for selective breeding programs of Asian seabass.

Fine mapping QTL for resistance to VNN disease using a high-density linkage map in Asian seabass

Asian seabass has suffered from viral nervous necrosis (VNN) disease. Our previous study has mapped quantitative trait loci (QTL) for resistance to VNN disease. To fine map these QTL and identify causative genes, we identified 6425 single nucleotide polymorphisms (SNPs) from 85 dead and 94 surviving individuals. Combined with 155 microsatellites, we constructed a genetic map consisting of 24 linkage groups (LGs) containing 3000 markers, with an average interval of 1.27 cM. We mapped one significant and three suggestive QTL with phenotypic variation explained (PVE) of 8.3 to 11.0%, two significant and two suggestive QTL with PVE of 7.8 to 10.9%, for resistance in three LGs and survival time in four LGs, respectively. Further analysis one QTL with the largest effect identified protocadherin alpha-C 2-like (Pcdhac2) as the possible candidate gene. Association study in 43 families with 1127 individuals revealed a 6 bp insertion-deletion was significantly associated with disease resistance. qRT-PCR showed the expression of Pcdhac2 was significantly induced in the brain, muscle and skin after nervous necrosis virus (NNV) infection. Our results could facilitate marker-assisted selection (MAS) for resistance to NNV in Asian seabass and set up the basis for functional analysis of the potential causative gene for resistance.

Construction of a high-density linkage map and fine mapping of QTL for growth in Asian seabass

A high-density genetic map is essential for comparative genomic studies and fine mapping of QTL, and can also facilitate genome sequence assembly. Here, a high density genetic map of Asian seabass was constructed with 3321 SNPs generated by sequencing 144 individuals in a F₂ family. The length of the map was 1577.67 cM with an average marker interval of 0.52 cM. A high level of genomic synteny among Asian seabass, European seabass, Nile tilapia and stickleback was detected. Using this map, one genome-wide significant and five suggestive QTL for growth traits were detected in six linkage groups (i.e. LG4, LG5, LG11, LG13, LG14 and LG15). These QTL explained 10.5–16.0% of phenotypic variance. A candidate gene, ACOX1 within the significant QTL on LG5 was identified. The gene was differentially expressed between fast- and slow-growing Asian seabass. The high-density SNP-based map provides an important tool for fine mapping QTL in molecular breeding and comparative genome analysis.

Microsatellite multiplex assay for the analysis of Atlantic sturgeon populations

We have developed a multiplex assay covering 16 microsatellite loci, amplified in four polymerase chain reaction (PCR) assays, and loaded on the ABI DNA Analyzer in two separate panels. The assay was tested on 603 individuals originating from wild populations and hatchery stocks of Atlantic sturgeon. The assay was also tested on 12 individuals of European sturgeon and appeared to be almost equally useful. The multiplex assay designed in this study can be successfully applied in studies requiring high genetic resolution, such as relatedness analysis, selective breeding programs, and stock identification of Atlantic sturgeon.

Tracing Asian seabass individuals to single fish farms using microsatellites

Traceability through physical labels is well established, but it is not highly reliable as physical labels can be easily changed or lost. Application of DNA markers to the traceability of food plays an increasingly important role for consumer protection and confidence building. In this study, we tested the efficiency of 16 polymorphic microsatellites and their combinations for tracing 368 fish to four populations where they originated. Using the maximum likelihood and Bayesian methods, three most efficient microsatellites were required to assign over 95% of fish to the correct populations. Selection of markers based on the assignment score estimated with the software WHICHLOCI was most effective in choosing markers for individual assignment, followed by the selection based on the allele number of individual markers. By combining rapid DNA extraction, and high-throughput genotyping of selected microsatellites, it is possible to conduct routine genetic traceability with high accuracy in Asian seabass.

Molecular parentage analysis is essential in breeding Asian seabass

In aquaculture species, maintaining pedigree information and genetic variation in each generation is essential, but very difficult. In this study, we used nine microsatellites to genotype 2,520 offspring from four independent full-factorial crosses (10 males ×10 females) of Asian seabass to reconstruct pedigree and monitor the change of genetic variations. In all four crosses, over 96.8% of the offspring could be assigned to their parents, indicating the high power of the nine microsatellites for parentage assignment. This study revealed several interesting results: (1). In all four crosses, the contribution of parents to offspring was significantly uneven, and some dominant breeding fishes (i.e. brooders) were found; (2). In two mass crosses where the brooders were carefully checked for reproductive status, a majority (≥90%) of brooders contributed to offspring, whereas in another two crosses, where the brooders were randomly picked without checking reproductive status, only a few brooders (40.0–45.0%) produced offspring; (3). Females had more problems in successful spawning compared to males; and (4). In the two crosses where a few brooders produced offspring, there was a substantial loss in allelic (24.1–34.3%) and gene (20.5–25.7%) diversities in offspring, while in the other two crosses, the majority of allelic (96.8–97.0%) and gene diversities (94.8–97.1%) were maintained. These observations suggest that a routine molecular parentage analysis is required to maintain both allelic and gene diversity in breeding Asian seabass.

Cryptocaryon irritans infection induces the acute phase response in Lates calcarifer: a transcriptomic perspective

Cryptocaryoniasis (also known as marine white spot disease) is mediated by Cryptocaryon irritans. This obligate ectoparasitic protozoan infects virtually all marine teleosts, which includes Lates calcarifer, a highly valuable aquaculture species. Little is known about L. calcarifer-C. irritans interactions. This study was undertaken to gain an informative snapshot of the L. calcarifer transcriptomic response over the course of C. irritans infection. An in-house fabricated cDNA microarray slides containing 3872 probes from L. calcarifer liver and spleen cDNA libraries were used as a tool to investigate the response of L. calcarifer to C. irritans infection. Juvenile fish were infected with parasites for four days, and total RNA was extracted from liver tissue, which was harvested daily. We compared the transcriptomes of C. irritans-infected liver to uninfected liver over an infection period of four days; the comparison was used to identify the genes with altered expression levels in response to C. irritans infection. The greatest number of infection-modulated genes was recorded at 2 and 3 days post-infection. These genes were mainly associated with the immune response and were associated in particular with the acute phase response. Acute phase proteins such as hepcidin, C-type lectin and serum amyloid A are among the highly modulated genes. Our results indicate that an induced acute phase response in L. calcarifer toward C. irritans infection is similar to the responses observed in bacterial infections of teleosts. This response demonstrates the importance of first line defenses in teleost innate immune responses against ectoparasite infection.

Evidence for female-biased dispersal in the protandrous hermaphroditic Asian Seabass, Lates calcarifer

Movement of individuals influences individual reproductive success, fitness, genetic diversity and relationships among individuals within populations and gene exchange among populations. Competition between males or females for mating opportunities and/or local resources predicts a female bias in taxa with monogamous mating systems and a male-biased dispersal in polygynous species. In birds and mammals, the patterns of dispersal between sexes are well explored, while dispersal patterns in protandrous hermaphroditic fish species have not been studied. We collected 549 adult individuals of Asian seabass (Lates calcarifer) from four locations in the South China Sea. To assess the difference in patterns of dispersal between sexes, we genotyped all individuals with 18 microsatellites. Significant genetic differentiation was detected among and within sampling locations. The parameters of population structure (F(ST)), relatedness (r) and the mean assignment index (mAIC), in combination with data on tagging-recapture, supplied strong evidences for female-biased dispersal in the Asian seabass. This result contradicts our initial hypothesis of no sex difference in dispersal. We suggest that inbreeding avoidance of females, female mate choice under the condition of low mate competition among males, and male resource competition create a female-biased dispersal. The bigger body size of females may be a cause of the female-biased movement. Studies of dispersal using data from DNA markers and tagging-recapture in hermaphroditic fish species could enhance our understanding of patterns of dispersal in fish.

Analysis of the Asian seabass transcriptome based on expressed sequence tags

Analysis of transcriptomes is of great importance in genomic studies. Asian seabass is an important fish species. A number of genomic tools in it were developed, while large expressed sequence tag (EST) data are lacking. We sequenced ESTs from nine normalized cDNA libraries and obtained 11 431 high-quality ESTs. We retrieved 8524 ESTs from dbEST database and analyzed all 19 975 ESTs using bioinformatics tools. After clustering, we obtained 8837 unique sequences (2838 contigs and 5999 singletons). The average contig length was 574 bp. Annotation of these unique sequences revealed that 48.9% of them showed significant homology to RNA sequences in GenBank. Functional classification of the unique ESTs identified a broad range of genes involved in different functions. We identified 6114 putative single-nucleotide polymorphisms and 634 microsatellites in ESTs. We discovered different temporal and spatial expression patterns of some immune-related genes in the Asian seabass after challenging with a pathogen Vibrio harveyi. The unique EST sequences are being used in developing a cDNA microarray to examine global gene expression and will also facilitate future whole-genome sequence assembly and annotation of Asian seabass and comparative genomics.

Identification of triploid individuals and clonal lines in Carassius auratus complex using microsatellites

The Carassius auratus complex in natural populations includes diploid triploid and polyploidy individuals. Diploid individuals belong to the species Carassius auratus whereas triploid and polyploidy individuals are from the subspecies Carassius auratus gibelio. Triploid individuals are all female and reproduce clonally by gynogenesis. Therefore the Carassius auratus complex is an ideal system for studying evolution of unisexual reproduction. Identification of triploid individuals and clonal lines is the first step towards understanding of the evolution of unisexual clonal lines. We examined the ability of 10 microsatellites in identifying triploid individuals in 94 individuals from Japan and China. In 40 confirmed triploid individuals and eight confirmed diploid individuals, all triploid and diploid individuals can be identified by genotyping 10 microsatellite, and four triploid clonal lines were identified. Using the 10 microsatellites we genotyped 46 adult individuals (40 females and six males) from a natural population in China and found that all six males were diploid whereas the majority of females (36 of 40) were triploid and three triploid clonal lines were detected. In 18 diploid individuals from China, all individuals showed different genotypes, suggesting there is no diploid clonal line in diploid crucian carp. A phylogenetic analysis of 94 individuals from China and Japan showed that triploid individuals and clonal lines have originated recurrently.