Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

A traits-based approach to assess aquaculture's contributions to food, climate change, and biodiversity goals

Aquaculture has the potential to support a sustainable and equitable food system in line with the United Nations Sustainable Development Goals (SDG) on food security, climate change, and biodiversity (FCB). Biological diversity amongst aquaculture organisms can drive diverse contributions to such goals. Existing studies have assessed the performance of a limited number of taxa in the general context of improving aquaculture production, but few explicitly consider the biological attributes of farmed aquatic taxa at the FCB nexus. Through a systematic literature review, we identify key traits associated with FCB and evaluate the potential of aquaculture to contribute to FCB goals using a fuzzy logic model. The majority of identified traits are associated with food security, and two-thirds of traits linked with food security are also associated with climate change or biodiversity, revealing potential co-benefits of optimizing a single trait. Correlations between FCB indices further suggest that challenges and opportunities in aquaculture are intertwined across FCB goals, but low mean FCB scores suggest that the focus of aquaculture research and development on food production is insufficient to address food security, much less climate or biodiversity issues. As expected, production-maximizing traits (absolute fecundity, the von Bertalanffy growth function coefficient K, macronutrient density, maximum size, and trophic level as a proxy for feed efficiency) highly influence a species' FCB potential, but so do species preferences for environmental conditions (tolerance to phosphates, nitrates, and pH levels, as well as latitudinal and geographic ranges). Many highly farmed species that are typically associated with food security, especially finfish, score poorly for food, climate, and biodiversity potential. Algae and mollusc species tend to perform well across FCB indices, revealing the importance of non-fish species in achieving FCB goals and potential synergies in integrated multi-trophic aquaculture systems. Overall, this study provides decision-makers with a biologically informed assessment of desirable aquaculture traits and species while illuminating possible strategies to increase support for FCB goals. Our findings can be used as a foundation for studying the socio-economic opportunities and barriers for aquaculture transitions to develop equitable pathways toward FCB-positive aquaculture across nuanced regional contexts.

Genome-wide QTL and eQTL mapping reveal genes associated with growth rate trait of the Pacific white shrimp (Litopenaeus vannamei)

BACKGROUND

Growth rate is a crucial economic trait for farmed animals, but the genetic regulation of this trait is largely unknown in non-model organisms such as shrimp.

RESULTS

In this study, we performed genome-wide phenotypic quantitative trait loci (QTL) and expression quantitative trait loci (eQTL) mapping analyses to identify genes affecting the growth rate of Pacific white shrimp (Litopenaeus vannamei), which is the most commercially-farmed crustacean worldwide. We used RNA-sequencing of 268 individuals in a mapping population, and subsequently validated our findings through gene silencing and shrimp growth experiments. We constructed a high-density genetic linkage map comprising 5533 markers spanning 44 linkage groups, with a total distance of 6205.75 cM and an average marker interval of 1.12 cM. Our analyses identified 11 QTLs significantly correlated with growth rate, and 117,525 eQTLs. By integrating QTL and eQTL data, we identified a gene (metalloreductase STEAP4) highly associated with shrimp growth rate. RNA interference (RNAi) analysis and growth experiments confirmed that STEAP4 was significantly correlated with growth rate in L. vannamei.

CONCLUSIONS

Our results indicate that the comprehensive analysis of QTL and eQTL can effectively identify genes involved in complex animal traits. This is important for marker-assisted selection (MAS) of animals. Our work contributes to the development of shrimp breeding and available genetic resources.

FishSNP: a high quality cross-species SNP database of fishes

The progress of aquaculture heavily depends on the efficient utilization of diverse genetic resources to enhance production efficiency and maximize profitability. Single nucleotide polymorphisms (SNPs) have been widely used in the study of aquaculture genomics, genetics, and breeding research since they are the most prevalent molecular markers on the genome. Currently, a large number of SNP markers from cultured fish species are scattered in individual studies, making querying complicated and data reuse problematic. We compiled relevant SNP data from literature and public databases to create a fish SNP database, FishSNP ( http://bioinfo.ihb.ac.cn/fishsnp ), and also used a unified analysis pipeline to process raw data that the author of the literature did not perform SNP calling on to obtain SNPs with high reliability. This database presently contains 45,690,243 (45 million) nonredundant SNP data for 13 fish species, with 30,288,958 (30 million) of those being high-quality SNPs. The main function of FishSNP is to search, browse, annotate and download SNPs, which provide researchers various and comprehensive associated information.

Unravelling the genomic secrets of bacterial fish pathogens: a roadmap to aquaculture sustainability

In the field of aquaculture, bacterial pathogens pose significant challenges to fish health and production. Advancements in genomic technologies have revolutionized our understanding of bacterial fish pathogens and their interactions with their host species. This review explores the application of genomic approaches in the identification, classification, and characterization of bacterial fish pathogens. Through an extensive analysis of the literature, we have compiled valuable data on 79 bacterial fish pathogens spanning 13 different phyla, encompassing their whole genome sequences. By leveraging high-throughput sequencing techniques, researchers have gained valuable insights into the genomic makeup of these pathogens, enabling a deeper understanding of their virulence factors and mechanisms of host interaction. Furthermore, genomic approaches have facilitated the discovery of potential vaccine and drug targets, opening up new avenues for the development of effective interventions against fish pathogens. Additionally, the utilization of genomics in fish disease resistance and control in aquaculture has shown promising results, enabling the identification of genetic markers associated with disease resistance traits. This review highlights the significant contributions of genomics to the field of fish pathogen research and underscores its potential for improving disease management strategies and enhancing the sustainability of aquaculture practices.

Versatile generation of precise gene edits in bovines using SEGCPN

BACKGROUND

Gene knockout and knock-in have been widely performed in large farm animals based on genome editing systems. However, many types of precise gene editing, including targeted deletion, gene tagging, and large gene fragment replacement, remain a challenge in large farm animals.

RESULTS

Here, we established versatile self-excising gene-targeting technology in combination with programmable nucleases (SEGCPN) to efficiently generate various types of precise gene editing in bovine. First, we used this versatile method to successfully generate bovine embryos with point mutations and 11-bp deletions at the MSTN locus. Second, we successfully generated bulls with EGFP labeling at the SRY locus. Finally, we successfully generated humanized cows in which the endogenous 18-kb α-casein gene was replaced with a 2.6-kb human α-lactalbumin gene.

CONCLUSIONS

In summary, our new SEGCPN method offers unlimited possibilities for various types of precise gene editing in large animals for application both in agriculture and disease models.

Reference genomes of channel catfish and blue catfish reveal multiple pericentric chromosome inversions

BACKGROUND

Channel catfish and blue catfish are the most important aquacultured species in the USA. The species do not readily intermate naturally but F₁ hybrids can be produced through artificial spawning. F₁ hybrids produced by mating channel catfish female with blue catfish male exhibit heterosis and provide an ideal system to study reproductive isolation and hybrid vigor. The purpose of the study was to generate high-quality chromosome level reference genome sequences and to determine their genomic similarities and differences.

RESULTS

We present high-quality reference genome sequences for both channel catfish and blue catfish, containing only 67 and 139 total gaps, respectively. We also report three pericentric chromosome inversions between the two genomes, as evidenced by long reads across the inversion junctions from distinct individuals, genetic linkage mapping, and PCR amplicons across the inversion junctions. Recombination rates within the inversional segments, detected as double crossovers, are extremely low among backcross progenies (progenies of channel catfish female × F₁ hybrid male), suggesting that the pericentric inversions interrupt postzygotic recombination or survival of recombinants. Identification of channel catfish- and blue catfish-specific genes, along with expansions of immunoglobulin genes and centromeric Xba elements, provides insights into genomic hallmarks of these species.

CONCLUSIONS

We generated high-quality reference genome sequences for both blue catfish and channel catfish and identified major chromosomal inversions on chromosomes 6, 11, and 24. These perimetric inversions were validated by additional sequencing analysis, genetic linkage mapping, and PCR analysis across the inversion junctions. The reference genome sequences, as well as the contrasted chromosomal architecture should provide guidance for the interspecific breeding programs.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Proteomic Applications in Aquatic Environment Studies

Genome determines the unique individualities of organisms; however, proteins play significant roles in the generation of the colorful life forms below water. Aquatic systems are usually complex and multifaceted and can take on unique modifications and adaptations to environmental changes by altering proteins at the cellular level. Proteomics is an essential strategy for exploring aquatic ecosystems due to the diverse involvement of proteins, proteoforms, and their complexity in basic and advanced cellular functions. Proteomics can expedite the analysis of molecular mechanisms underlying biological processes in an aquatic environment. Previous proteomic studies on aquatic environments have mainly focused on pollution assessments, ecotoxicology, their role in the food industry, and extraction and identification of natural products. Aquatic protein biomarkers have been comprehensively reported and are currently extensively applied in the pharmaceutical and medical industries. Cellular- and molecular-level responses of organisms can be used as indicators of environmental changes and stresses. Conversely, environmental changes are expedient in predicting aquatic health and productivity, which are crucial for ecosystem management and conservation. Recent advances in proteomics have contributed to the development of sustainable aquaculture, seafood safety, and high aquatic food production. Proteomic approaches have expanded to other aspects of the aquatic environment, such as protein fingerprinting for species identification. In this review, we encapsulated current proteomic applications and evaluated the potential strengths, weaknesses, opportunities, and threats of proteomics for future aquatic environmental studies. The review identifies both pros and cons of aquatic proteomics and projects potential challenges and recommendations. We postulate that proteomics is an emerging, powerful, and integrated omics approach for aquatic environmental studies.

Development of a high-throughput SNP array for sea cucumber (Apostichopus japonicus) and its application in genomic selection with MCP regularized deep neural networks

High-throughput single nucleotide polymorphism (SNP) genotyping assays are powerful tools for genetic studies and genomic breeding applications for many species. Though large numbers of SNPs have been identified in sea cucumber (Apostichopus japonicus), but, as yet, no high-throughput genotyping platform is available for this species. In this study, we designed and developed a high-throughput 24 K SNP genotyping array named HaishenSNP24K for A. japonicus, based on the multi-objective-local optimization (MOLO) algorithm and HD-Marker genotyping method. The SNP array exhibited a relatively high genotyping call rate (> 96%), genotyping accuracy (>95%) and exhibited highly polymorphic in sea cucumber populations. In addition, we also assessed its application in genomic selection (GS). Deep neural networks (DNN) that can capture the complicated interactions of genes have been proposed as a promising tool in GS for SNP-based genomic prediction of complex traits in animal breeding. To overcome the problem of over-fitting when using the HaishenSNP24K array as high-dimensional DNN input, we developed minmax concave penalty (MCP) regularization for sparse deep neural networks (DNN-MCP) that finds an optimal sparse structure of a DNN by minimizing the square error subject to the non-convex penalty MCP on the parameters (weights and biases). Compared to two linear models, namely RR-GBLUP and Bayes B, and the nonlinear model DNN, DNN-MCP has greatly improved the genomic prediction ability for three quantitative traits (e.g., wet weight, dry weight and survival time) in the sea cucumber population. To the best of our knowledge, this is the first work to develop a high-throughput SNP array for A. japonicus and a new model DNN-MCP for genomic prediction of complex traits in GS. The present results provide evidence that supports the HaishenSNP24K array with DNN-MCP will be valuable for genetic studies and molecular breeding in A. japonicus.

Pentaploidization Enriches the Genetic Diversity of Wheat by Enhancing the Recombination of AB Genomes

Allohexaploidization and continuous introgression play a key role in the origin and evolution of bread wheat. The genetic bottleneck of bread wheat resulting from limited germplasms involved in the origin and modern breeding may be compensated by gene flow from tetraploid wheat through introgressive hybridization. The inter-ploidy hybridization between hexaploid and tetraploid wheat generates pentaploid hybrids first, which absorbed genetic variations both from hexaploid and tetraploid wheat and have great potential for re-evolution and improvement in bread wheat. Therefore, understanding the effects of the pentaploid hybrid is of apparent significance in our understanding of the historic introgression and in informing breeding. In the current study, two sets of F2 populations of synthetic pentaploid wheat (SPW1 and SPW2) and synthetic hexaploid wheat (SHW1 and SHW2) were created to analyze differences in recombination frequency (RF) of AB genomes and distorted segregation of polymorphic SNP markers through SNP genotyping. Results suggested that (1) the recombination of AB genomes in the SPW populations was about 3- to 4-fold higher than that in the SHW populations, resulting from the significantly (P < 0.01) increased RF between adjacent and linked SNP loci, especially the variations that occurred in a pericentromeric region which would further enrich genetic diversity; (2) the crosses of hexaploid × tetraploid wheat could be an efficient way to produce pentaploid derivatives than the crosses of tetraploid × hexaploid wheat according to the higher germination rate found in the former crosses; (3) the high proportion of distorted segregation loci that skewed in favor of the female parent genotype/allele in the SPW populations might associate with the fitness and survival of the offspring. Based on the presented data, we propose that pentaploid hybrids should increasingly be used in wheat breeding. In addition, the contribution of gene flow from tetraploid wheat to bread wheat mediated by pentaploid introgressive hybridization also was discussed in the re-evolution of bread wheat.

Integrated Analyses of DNA Methylation and Gene Expression of Rainbow Trout Muscle under Variable Ploidy and Muscle Atrophy Conditions

Rainbow trout, Oncorhynchus mykiss, is an important cool, freshwater aquaculture species used as a model for biological research. However, its genome reference has not been annotated for epigenetic markers affecting various biological processes, including muscle growth/atrophy. Increased energetic demands during gonadogenesis/reproduction provoke muscle atrophy in rainbow trout. We described DNA methylation and its associated gene expression in atrophying muscle by comparing gravid, diploid females to sterile, triploid females. Methyl Mini-seq and RNA-Seq were simultaneously used to characterize genome-wide DNA methylation and its association with gene expression in rainbow trout muscle. Genome-wide enrichment in the number of CpGs, accompanied by depleted methylation levels, was noticed around the gene transcription start site (TSS). Hypermethylation of CpG sites within ±1 kb on both sides of TSS (promoter and gene body) was weakly/moderately associated with reduced gene expression. Conversely, hypermethylation of the CpG sites in downstream regions of the gene body +2 to +10 kb was weakly associated with increased gene expression. Unlike mammalian genomes, rainbow trout gene promotors are poor in CpG islands, at <1% compared to 60%. No signs of genome-wide, differentially methylated (DM) CpGs were observed due to the polyploidy effect; only 1206 CpGs (0.03%) were differentially methylated, and these were primarily associated with muscle atrophy. Twenty-eight genes exhibited differential gene expression consistent with methylation levels of 31 DM CpGs. These 31 DM CpGs represent potential epigenetic markers of muscle atrophy in rainbow trout. The DM CpG-harboring genes are involved in apoptosis, epigenetic regulation, autophagy, collagen metabolism, cell membrane functions, and Homeobox proteins. Our study also identified genes explaining higher water content and modulated glycolysis previously shown as characteristic biochemical signs of rainbow trout muscle atrophy associated with sexual maturation. This study characterized DNA methylation in the rainbow trout genome and its correlation with gene expression. This work also identified novel epigenetic markers associated with muscle atrophy in fish/lower vertebrates.

A High-Quality Genome Assembly of Striped Catfish (Pangasianodon hypophthalmus) Based on Highly Accurate Long-Read HiFi Sequencing Data

The HiFi sequencing technology yields highly accurate long-read data with accuracies greater than 99.9% that can be used to improve results for complex applications such as genome assembly. Our study presents a high-quality chromosome-scale genome assembly of striped catfish (Pangasianodon hypophthalmus), a commercially important species cultured mainly in Vietnam, integrating HiFi reads and Hi-C data. A 788.4 Mb genome containing 381 scaffolds with an N50 length of 21.8 Mb has been obtained from HiFi reads. These scaffolds have been further ordered and clustered into 30 chromosome groups, ranging from 1.4 to 57.6 Mb, based on Hi-C data. The present updated assembly has a contig N50 of 14.7 Mb, representing a 245-fold and 4.2-fold improvement over the previous Illumina and Illumina-Nanopore-Hi-C based version, respectively. In addition, the proportion of repeat elements and BUSCO genes identified in our genome is remarkably higher than in the two previously released striped catfish genomes. These results highlight the power of using HiFi reads to assemble the highly repetitive regions and to improve the quality of genome assembly. The updated, high-quality genome assembled in this work will provide a valuable genomic resource for future population genetics, conservation biology and selective breeding studies of striped catfish.

Gene Editing of the Catfish Gonadotropin-Releasing Hormone Gene and Hormone Therapy to Control the Reproduction in Channel Catfish, Ictalurus punctatus

Transcription activator-like effector nuclease (TALEN) plasmids targeting the channel catfish gonadotropin-releasing hormone (cfGnRH) gene were delivered into fertilized eggs with double electroporation to sterilize channel catfish (Ictalurus punctatus). Targeted cfGnRH fish were sequenced and base deletion, substitution, and insertion were detected. The gene mutagenesis was achieved in 52.9% of P1 fish. P1 mutants (individuals with human-induced sequence changes at the cfGnRH locus) had lower spawning rates (20.0−50.0%) when there was no hormone therapy compared to the control pairs (66.7%) as well as having lower average egg hatch rates (2.0% versus 32.3−74.3%) except for one cfGnRH mutated female that had a 66.0% hatch rate. After low fertility was observed in 2016, application of luteinizing hormone-releasing hormone analog (LHRHa) hormone therapy resulted in good spawning and hatch rates for mutants in 2017, which were not significantly different from the controls (p > 0.05). No exogenous DNA fragments were detected in the genome of mutant P1 fish, indicating no integration of the plasmids. No obvious effects on other economically important traits were observed after the knockout of the reproductive gene in the P1 fish. Growth rates, survival, and appearance between mutant and control individuals were not different. While complete knock-out of reproductive output was not achieved, as these were mosaic P1 brood stock, gene editing of channel catfish for the reproductive confinement of gene-engineered, domestic, and invasive fish to prevent gene flow into the natural environment appears promising.

New Insight on Vitality Differences for the Penaeid Shrimp, Fenneropenaeus chinensis, in Low Salinity Environment Through Transcriptomics

Excessive rainfall changes salinity in shrimp farming ponds in short period and exerts low salinity stress on the outdoor breeding shrimp under global warming. Fenneropenaeus chinensis can have different performance on vitality in low salinity environments. To reveal mechanisms of vitality difference in shrimp living in low saline environments. This study based on the normal and moribund F. chinensis in 10 ppt salinity environment using high-throughput sequencing identifies 1,429 differentially expressed genes (DEGs), 586 of which are upregulated, while 843 of which are downregulated in the normal group (FCN10) as compared to the moribund group (FCM10). Meanwhile, another transcriptomic analysis is conducted on the normal and moribund shrimp from 25 ppt (FCN25 vs. FCM25) salinity environment as the control, in which 1,311 DEGs (upregulated: 327 genes, downregulated: 984 genes) are identified. In this study, intersective pathways, GO (Gene Ontology) categories and DEGs from the two groups of comparative transcriptome are investigated. The two intersective pathways (Metabolism of xenobiotics by cytochrome P450, Pentose, and glucuronate interconversions) significantly enriched by DEGs are related to detoxification. In these two pathways, there is one vitality regulation-related gene (VRRG), the Dhdh (dihydrodiol dehydrogenase), which is upregulated in both the groups of FCN10 and FCN25 as compared to the groups of FCM10 and FCM25, respectively. Similarly, in the 25 top intersective GO categories, four VRRGs are revealed. Three of them are upregulated (Itgbl, kielin/chordin-like protein, Slc2a8, solute carrier family 2, facilitated glucose transporter member 8-like protein and Cyp3a30, cytochrome P450 3A30-like protein); one of them is downregulated (Slc6a9, sodium-dependent nutrient amino acid transporter 1-like protein isoform X2). These GO categories are related to transmembrane transporter activity of substance, enzyme inhibitor activity, monooxygenase activity. RT-qPCR analysis further verifies the VRRGs. The study gives new insight into understanding the vitality differences for F. chinensis, in low salinity environment. The pathways and DEGs in response to low salinity stress in modulating the vitality of F. chinensis that could serve as tools in future genetic studies and molecular breeding.

Integrated application of multi-omics strategies provides insights into the environmental hypoxia response in Pelteobagrus vachelli muscle

Increasing pressures on aquatic ecosystems because of pollutants, nutrient enrichment, and global warming have severely depleted oxygen concentrations. This sudden and significant lack of oxygen has resulted in persistent increases in fish mortality rates. Revealing the molecular mechanism of fish hypoxia adaptation will help researchers to find markers for hypoxia induced by environmental stress. Here, we used a multiomics approach to identify several hypoxia-associated miRNAs, mRNAs, proteins, and metabolites involved in diverse biological pathways in the muscles of Pelteobagrus vachelli. Our findings revealed significant hypoxia-associated changes in muscles over 4 h of hypoxia exposure and discrete tissue-specific patterns. We have previously reported that P. vachelli livers exhibit increased anaerobic glycolysis, heme synthesis, erythropoiesis, and inhibit apoptosis when exposed to hypoxia for 4 h. However, the opposite was observed in muscles. According to our comprehensive analysis, fishes show an acute response to hypoxia, including activation of catabolic pathways to generate more energy, reduction of biosynthesis to decrease energy consumption, and shifting from aerobic to anaerobic metabolic contributions. Also, we found that hypoxia induced muscle dysfunction by impairing mitochondrial function, activating inflammasomes, and apoptosis. The hypoxia-induced mitochondrial dysfunction enhanced oxidative stress, apoptosis, and further triggered interleukin-1β production via inflammasome activation. In turn, interleukin-1β further impaired mitochondrial function or apoptosis by suppressing downstream mitochondrial biosynthesis–related proteins, thus resulting in a vicious cycle of inflammasome activation and mitochondrial dysfunction. Our findings contribute meaningful insights into the molecular mechanisms of hypoxia, and the methods and study design can be utilized across different fish species.

•

First multiomics analysis of mRNA, miRNA, protein, and metabolite in fishes.

•

Liver and muscle were tissue-specific induced by hypoxia.

•

About 70 genes and 16 miRNAs related to hypoxia adaptation were detected.

•

Hypoxia affects muscle function by mediating energy metabolism via HIF pathway.

Genomic Selection in Aquaculture Species

To date, genomic prediction has been conducted in about 20 aquaculture species, with a preference for intra-family genomic selection (GS). For every trait under GS, the increase in accuracy obtained by genomic estimated breeding values instead of classical pedigree-based estimation of breeding values is very important in aquaculture species ranging from 15% to 89% for growth traits, and from 0% to 567% for disease resistance. Although the implementation of GS in aquaculture is of little additional investment in breeding programs already implementing sib testing on pedigree, the deployment of GS remains sparse, but could be boosted by adaptation of cost-effective imputation from low-density panels. Moreover, GS could help to anticipate the effect of climate change by improving sustainability-related traits such as production yield (e.g., carcass or fillet yields), feed efficiency or disease resistance, and by improving resistance to environmental variation (tolerance to temperature or salinity variation). This chapter synthesized the literature in applications of GS in finfish, crustaceans and molluscs aquaculture in the present and future breeding programs.

Zika Virus Overview: Transmission, Origin, Pathogenesis, Animal Model and Diagnosis

Zika virus (ZIKV) was first discovered in 1947 in Uganda. ZIKV did not entice much attention until Brazil hosted the 2016 Summer Olympics Game, where ZIKV attracted a global audience. ZIKV is a flavivirus that can be transmitted chiefly through the biting of the mosquito or sexually or by breastfeeding at a lower scale. As time passed, the recent discovery of how the ZIKV causes congenital neurodevelopmental defects, including microcephaly, makes us reevaluate the importance of ZIKV interaction with centrosome organization because centrosome plays an important role in cell division. When the ZIKV disrupts centrosome organization and mitotic abnormalities, this will alter neural progenitor differentiation. Altering the neural progenitor differentiation will lead to cell cycle arrest, increase apoptosis, and inhibit the neural progenitor cell differentiation, as this can lead to abnormalities in neural cell development resulting in microcephaly. Understanding the importance of ZIKV infection throughout the years, this review article gives an overview of the history, transmission routes, pathogenesis, animal models, and diagnosis.

Cross-species Exon Capture and Whole Exome Sequencing: Application, Utility and Challenges for Genomic Resource Development in Non-model Species

Comprehending the genetic architecture of complex traits has many applications in evolution, ecology, conservation biology and plant and animal production systems. Underlying research questions in these fields are diverse species that often have limited genetic information available. In aquaculture, for example, genetic progress has been slow in many species due to a lack in such genetic information. In this study, zebrafish (as a well-studied model species) was used in cross-species transfer to develop genomic resources and identify candidate genes underling growth differentials in dusky kob. Dusky kob is a Sciaenid finfish and an emerging aquaculture species. The zebrafish All Exon Predesigned Probe-set capture protocol was used to enrich fractionated DNA samples from kob, classified as either large or small, before massive parallel sequencing on the Ion Torrent platform. Although vast quantities of sequence data were generated, only about 30% of contigs could be identified as zebrafish homologues. There were numerous species-specific sequences and inconsistent coverage of sequencing products across samples, likely due to non-specific binding of the probe-set as a result of the evolutionary divergence between zebrafish and kob. Nonetheless, more than 55,000 SNPs could be reliably identified and genotyped to the individual level. Using SNP genotypic divergence estimates, between large and small cohorts, a number of candidate genes associated with growth was also identified for future investigation. These findings contribute to the growing body of evidence demonstrating the utility of a cross-species capture approach in the development of important genomic resources for understanding traits of interest in species without reference genomes.

Genome-wide association analysis of adaptation to oxygen stress in Nile tilapia (Oreochromis niloticus)

Background

Tilapia is one of the most abundant species in aquaculture. Hypoxia is known to depress growth rate, but the genetic mechanism by which this occurs is unknown. In this study, two groups consisting of 3140 fish that were raised in either aerated (normoxia) or non-aerated pond (nocturnal hypoxia). During grow out, fish were sampled five times to determine individual body weight (BW) gains. We applied a genome-wide association study to identify SNPs and genes associated with the hypoxic and normoxic environments in the 16th generation of a Genetically Improved Farmed Tilapia population.

Results

In the hypoxic environment, 36 SNPs associated with at least one of the five body weight measurements (BW1 till BW5), of which six, located between 19.48 Mb and 21.04 Mb on Linkage group (LG) 8, were significant for body weight in the early growth stage (BW1 to BW2). Further significant associations were found for BW in the later growth stage (BW3 to BW5), located on LG1 and LG8. Analysis of genes within the candidate genomic region suggested that MAPK and VEGF signalling were significantly involved in the later growth stage under the hypoxic environment. Well-known hypoxia-regulated genes such as igf1rb, rora, efna3 and aurk were also associated with growth in the later stage in the hypoxic environment. Conversely, 13 linkage groups containing 29 unique significant and suggestive SNPs were found across the whole growth period under the normoxic environment. A meta-analysis showed that 33 SNPs were significantly associated with BW across the two environments, indicating a shared effect independent of hypoxic or normoxic environment. Functional pathways were involved in nervous system development and organ growth in the early stage, and oocyte maturation in the later stage.

Conclusions

There are clear genotype-growth associations in both normoxic and hypoxic environments, although genome architecture involved changed over the growing period, indicating a transition in metabolism along the way. The involvement of pathways important in hypoxia especially at the later growth stage indicates a genotype-by-environment interaction, in which MAPK and VEGF signalling are important components.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07486-5.

A de novo Full-Length mRNA Transcriptome Generated From Hybrid-Corrected PacBio Long-Reads Improves the Transcript Annotation and Identifies Thousands of Novel Splice Variants in Atlantic Salmon

Atlantic salmon (Salmo salar) is a major species produced in world aquaculture and an important vertebrate model organism for studying the process of rediploidization following whole genome duplication events (Ss4R, 80 mya). The current Salmo salar transcriptome is largely generated from genome sequence based in silico predictions supported by ESTs and short-read sequencing data. However, recent progress in long-read sequencing technologies now allows for full-length transcript sequencing from single RNA-molecules. This study provides a de novo full-length mRNA transcriptome from liver, head-kidney and gill materials. A pipeline was developed based on Iso-seq sequencing of long-reads on the PacBio platform (HQ reads) followed by error-correction of the HQ reads by short-reads from the Illumina platform. The pipeline successfully processed more than 1.5 million long-reads and more than 900 million short-reads into error-corrected HQ reads. A surprisingly high percentage (32%) represented expressed interspersed repeats, while the remaining were processed into 71 461 full-length mRNAs from 23 071 loci. Each transcript was supported by several single-molecule long-read sequences and at least three short-reads, assuring a high sequence accuracy. On average, each gene was represented by three isoforms. Comparisons to the current Atlantic salmon transcripts in the RefSeq database showed that the long-read transcriptome validated 25% of all known transcripts, while the remaining full-length transcripts were novel isoforms, but few were transcripts from novel genes. A comparison to the current genome assembly indicates that the long-read transcriptome may aid in improving transcript annotation as well as provide long-read linkage information useful for improving the genome assembly. More than 80% of transcripts were assigned GO terms and thousands of transcripts were from genes or splice-variants expressed in an organ-specific manner demonstrating that hybrid error-corrected long-read transcriptomes may be applied to study genes and splice-variants expressed in certain organs or conditions (e.g., challenge materials). In conclusion, this is the single largest contribution of full-length mRNAs in Atlantic salmon. The results will be of great value to salmon genomics research, and the pipeline outlined may be applied to generate additional de novo transcriptomes in Atlantic Salmon or applied for similar projects in other species.

The Chinese mitten crab genome provides insights into adaptive plasticity and developmental regulation

The infraorder Brachyura (true or short-tailed crabs) represents a successful group of marine invertebrates yet with limited genomic resources. Here we report a chromosome-anchored reference genome and transcriptomes of the Chinese mitten crab Eriocheir sinensis, a catadromous crab and invasive species with wide environmental tolerance, strong osmoregulatory capacity and high fertility. We show the expansion of specific gene families in the crab, including F-ATPase, which enhances our knowledge on the adaptive plasticity of this successful invasive species. Our analysis of spatio-temporal transcriptomes and the genome of E. sinensis and other decapods shows that brachyurization development is associated with down-regulation of Hox genes at the megalopa stage when tail shortening occurs. A better understanding of the molecular mechanism regulating sexual development is achieved by integrated analysis of multiple omics. These genomic resources significantly expand the gene repertoire of Brachyura, and provide insights into the biology of this group, and Crustacea in general.

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.

Identification of Growth-Associated Genes by Genome-Wide Association Study and Their Potential Application in the Breeding of Pacific White Shrimp (Litopenaeus vannamei)

The Pacific white shrimp (Litopenaeus vannamei) is the most widely cultured shrimp in the world. A great attention has been paid to improve its body weight (BW) at harvest through genetic selection for decades. Genome-wide association study (GWAS) is a tool to dissect the genetic basis of the traits. In this study, a GWAS approach was conducted to find genes related to BW through genotyping 94,113 single nucleotide polymorphisms (SNPs) in 200 individuals from a breeding population. Four BW-related SNPs located in LG19 and LG39 were identified. Through further candidate gene association analysis, the SNPs in two candidate genes, deoxycytidylate deaminase and non-receptor protein tyrosine kinase, were found to be related with the body weight of the shrimp. Marker-assisted best linear unbiased prediction (MA-BLUP) based on the SNPs in these two genes was used to estimate the breeding values, and the result showed that the highest prediction accuracy of MA-BLUP was increased by 9.4% than traditional BLUP. These results will provide useful information for the marker-assisted breeding in L. vannamei.

Direct and pleiotropic effects of the Masou Salmon Delta-5 Desaturase transgene in F1 channel catfish (Ictalurus punctatus)

Channel catfish (Ictalurus punctatus) is the primary culture species in the US along with its hybrid made with male blue catfish, I. furcatus. In an effort to improve the nutritional value of channel catfish, the masou salmon Δ5-desaturase like gene (D5D) driven by the common carp beta-actin promoter (βactin) was inserted into channel catfish. The objectives of this study were to determine the effectiveness of βactin-D5D for improving n-3 fatty acid production in F1 transgenic channel catfish, as well as examine pleiotropic effects on growth, proximate analysis, disease resistance, and other performance traits. Transgenic F1 channel catfish showed a 33% increase in the relative proportion of n-3 fatty acids coupled with a 15% decrease in n-6 fatty acids and a 17% decrease in n-9 fatty acids when compared to non-transgenic full-siblings (P < 0.01, P < 0.01, P < 0.01). However, while the relative proportion of n-3 fatty acids was achieved, the total amount of fatty acids in the transgenic fish decreased resulting in a reduction of all fatty acids. Insertion of the βactin-D5D transgene into channel catfish also had large effects on the body composition, and growth of channel catfish. Transgenic channel catfish grew faster, were more disease resistant, had higher protein and moisture percentage, but lower fat percentage than full-sib controls. There were sex effects as performance changes were more dramatic and significant in males. The βactin-D5D transgenic channel catfish were also more uniform in their fatty acid composition, growth and other traits.

Differential expression analysis and identification of sex-related genes by gonad transcriptome sequencing in estradiol-treated and non-treated Ussuri catfish Pseudobagrus ussuriensis

The Ussuri catfish (Pseudobagrus ussuriensis) has an XX/XY sex determination system but its sex determination gene(s) remain unknown. To better understand the molecular sex determination mechanism, transcriptome analysis was conducted to obtain sex-related gene expression profiles. Transcriptome analyses were made of male and female developing/differentiating gonads by high-throughput RNA sequencing, including gonads from fish given an estradiol-induced sex reversal treatment. A total of 81,569 unigenes were assembled and 39,904 were significantly matched to known unique proteins by comparison with public databases. Twenty specifically expressed and 142 differentially expressed sex-related genes were extracted from annotated data by comparing the treatment groups. These genes are involved in spermatogenesis (e.g., Dnali1, nectin3, klhl10, mybl1, Katnal1, Eno4, Mns1, Spag6, Tsga10, Septin7), oogenesis (e.g., Lagr5, Fmn2, Npm2, zar1, Fbxo5, Fbxo43, Prdx4, Nrip1, Lfng, Atrip), gonadal development/differentiation (e.g., Cxcr4b, Hmgb2, Cftr, Ch25h, brip1, Prdm9, Tdrd1, Star, dmrt1, Tut4, Hsd17b12a, gdf9, dnd, arf1, Spata22), and estradiol response (e.g., Mmp14, Lhcgr, vtg1, vtg2, esr2b, Piwil1, Aifm1, Hsf1, gdf9). Dmrt1 and gdf9 may play an essential role in sex determination in P. ussuriensis. The expression patterns of six random genes were validated by quantitative real-time PCR, which confirmed the reliability and accuracy of the RNA-seq results. These data provide a valuable resource for future studies of gene expression and for understanding the molecular mechanism of sex determination/differentiation and gonadal development/differentiation (including hormone-induced sexual reversal) in Ussuri catfish. This has the potential to assist in producing monosex Ussuri catfish to increase aquacultural productivity.

Genetic Parameters and Genome-Wide Association Studies of Quality Traits Characterised Using Imaging Technologies in Rainbow Trout, Oncorhynchus mykiss

One of the top priorities of the aquaculture industry is the genetic improvement of economically important traits in fish, such as those related to processing and quality. However, the accuracy of genetic evaluations has been hindered by a lack of data on such traits from a sufficiently large population of animals. The objectives of this study were thus threefold: (i) to estimate genetic parameters of growth-, yield-, and quality-related traits in rainbow trout (Oncorhynchus mykiss) using three different phenotyping technologies [invasive and non-invasive: microwave-based, digital image analysis, and magnetic resonance imaging (MRI)], (ii) to detect quantitative trait loci (QTLs) associated with these traits, and (iii) to identify candidate genes present within these QTL regions. Our study collected data from 1,379 fish on growth, yield-related traits (body weight, condition coefficient, head yield, carcass yield, headless gutted carcass yield), and quality-related traits (total fat, percentage of fat in subcutaneous adipose tissue, percentage of fat in flesh, flesh colour); genotypic data were then obtained for all fish using the 57K SNP Axiom® Trout Genotyping array. Heritability estimates for most of the 14 traits examined were moderate to strong, varying from 0.12 to 0.67. Most traits were clearly polygenic, but our genome-wide association studies (GWASs) identified two genomic regions on chromosome 8 that explained up to 10% of the genetic variance (cumulative effects of two QTLs) for several traits (weight, condition coefficient, subcutaneous and total fat content, carcass and headless gutted carcass yields). For flesh colour traits, six QTLs explained 1-4% of the genetic variance. Within these regions, we identified several genes (htr1, gnpat, ephx1, bcmo1, and cyp2x) that have been implicated in adipogenesis or carotenoid metabolism, and thus represent good candidates for further functional validation. Finally, of the three techniques used for phenotyping, MRI demonstrated particular promise for measurements of fat content and distribution, while the digital image analysis-based approach was very useful in quantifying colour-related traits. This work provides new insights that may aid the development of commercial breeding programmes in rainbow trout, specifically with regard to the genetic improvement of yield and flesh-quality traits as well as the use of invasive and/or non-invasive technologies to predict such traits.

Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper (Epinephelus tukula)

Simple Summary

The potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits.

Abstract

The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F₁ juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.

An ultra-high density SNP-based linkage map for enhancing the pikeperch (Sander lucioperca) genome assembly to chromosome-scale

Pikeperch (Sander lucioperca) is a fish species with growing economic significance in the aquaculture industry. However, successful positioning of pikeperch in large-scale aquaculture requires advances in our understanding of its genome organization. In this study, an ultra-high density linkage map for pikeperch comprising 24 linkage groups and 1,023,625 single nucleotide polymorphisms markers was constructed after genotyping whole-genome sequencing data from 11 broodstock and 363 progeny, belonging to 6 full-sib families. The sex-specific linkage maps spanned a total of 2985.16 cM in females and 2540.47 cM in males with an average inter-marker distance of 0.0030 and 0.0026 cM, respectively. The sex-averaged map spanned a total of 2725.53 cM with an average inter-marker distance of 0.0028 cM. Furthermore, the sex-averaged map was used for improving the contiguity and accuracy of the current pikeperch genome assembly. Based on 723,360 markers, 706 contigs were anchored and oriented into 24 pseudomolecules, covering a total of 896.48 Mb and accounting for 99.47% of the assembled genome size. The overall contiguity of the assembly improved with a scaffold N50 length of 41.06 Mb. Finally, an updated annotation of protein-coding genes and repetitive elements of the enhanced genome assembly is provided at NCBI.

Effects of gonadotropin-releasing hormone analog (GnRHa) immunization on the gonadal transcriptome and proteome of tilapia (Oreochromis niloticus)

Gonadotropin releasing hormone (GnRH) plays an important role in the regulation of vertebrate reproduction. Studies have shown that immunization against GnRHa can induce sexually sterile tilapia. To explore the mechanism behind this, in this study, RNA-seq and data-independent acquisition (DIA) techniques were used to study the transcriptome and proteome of the gonad of tilapia immunized with GnRHa. 644 differentially expressed genes (80 upregulated and 564 downregulated) and 1150 differentially expressed proteins (351 upregulated and 799 downregulated) were identified. There were 209 genes with consistent differential expression patterns in the transcriptomic and proteomic analyses, of which 9 were upregulated and 200 downregulated, indicating that the gonad gene expression was inhibited by GnRHa immunization. The downregulated genes were particularly involved in the functions of single-organism process, binding, cellular process, metabolic process and catalytic activity, and associated with the pathways including ECM-receptor interaction, focal adhesion, cardiac muscle contraction and oxidative phosphorylation. The expression of six differentially expressed genes involved in the GnRH signaling pathway was all downregulated. In addition, several important functional genes related to gonadal development after GnRHa immunization were screened. This study confirmed the expression of corresponding genes was affected by GnRHa on the gonad development in tilapia at the molecular level, and laid a foundation for elucidating the mechanism of GnRHa immunization.

Development and Evaluation of a High-Throughput Single-Nucleotide Polymorphism Array for Large Yellow Croaker (Larimichthys crocea)

High-density single-nucleotide polymorphism (SNP) genotyping array is an essential tool for genetic analyses of animals and plants. Large yellow croaker (Larimichthys crocea) is one of the most commercially important marine fish species in China. Although plenty of SNPs have been identified in large yellow croaker, no high-throughput genotyping array is available. In this study, a high-throughput SNP array named NingXin-I with 600K SNPs was developed and evaluated. A set of 82 large yellow croakers were collected from different locations of China and re-sequenced. A total of 9.34M SNPs were identified by mapping sequence reads to the large yellow croaker reference genome. About 1.98M candidate SNPs were selected for further analyses by using criteria such as SNP quality score and conversion performance in the final array. Finally, 579.5K SNPs evenly distributed across the large yellow croaker genome with an average spacing of 1.19 kb were proceeded to array production. The performance of NingXin-I array was evaluated in 96 large yellow croaker individuals from five populations, and 83.38% SNPs on the array were polymorphic sites. A further test of the NingXin-I array in five closely related species in Sciaenidae identified 26.68–56.23% polymorphic SNP rate across species. A phylogenetic tree inferred by using the genotype data generated by NingXin-I confirmed the phylogenetic distance of the species in Sciaenidae. The performance of NingXin-I in large yellow croaker and the other species in Sciaenidae suggested high accuracy and broad application. The NingXin-I array should be valuable for quantitative genetic studies, such as genome-wide association studies (GWASs), high-density linkage map construction, haplotype analysis, and genome-based selection.

Full-length transcriptome analysis provides new insights into the innate immune system of Marsupenaeus japonicus

As invertebrates, shrimp are generally thought to solely rely on their innate immune system to combat invading pathogens. Recently, an increasing number of studies have revealed that the innate immune response of invertebrates exhibits diversity and specificity based on their diverse immune molecules. Herein, a full-length transcriptome analysis of several immune-related tissues (hepatopancreas, gill, hemocytes, stomach and intestine) in the kuruma shrimp (Marsupenaeus japonicus) was conducted to identify immune-related molecules with a focus on transcript variations. In total, 11,222 nonredundant full-length transcripts with an N50 length of 5174 were obtained, and most of these transcripts (94.84%) were successfully annotated. In addition, a total of 147 long noncoding RNAs (lncRNAs) were also predicted. Importantly, transcript variants of several vital immune-related genes were observed, including twenty-five alpha-2-macroglobulins (α2-Ms), ten Toll-like receptors (TLRs), six C-type lectins (CTLs), five M-type lectins (MTLs) and three Down syndrome cell adhesion molecules (Dscams). Furthermore, 509 nonredundant full-length transcripts were predicted to be generated from alternative splicing (AS) events, which contribute to the diversity of immune molecules. Overall, our study provides valuable data on the full-length transcripts of M. japonicus, which will facilitate the exploration of immune molecules in this species. Moreover, numerous transcript variants of immune molecules detected in this study provide clues for further investigating the diversity and specificity of the innate immune response in shrimp.

Twelve quick steps for genome assembly and annotation in the classroom

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

Construction of a High-Density Genetic Map and Identification of Quantitative Trait Loci for Nitrite Tolerance in the Pacific White Shrimp (Litopenaeus vannamei)

Nitrite is a major environmental toxin in aquaculture systems that disrupts multiple physiological functions in aquatic animals. Although nitrite tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available. In this study, we constructed a high-density genetic map of Litopenaeus vannamei with 17,242 single nucleotide polymorphism markers spanning 6,828.06 centimorgans (cM), with an average distance of 0.4 cM between adjacent markers on 44 linkage groups (LGs). Using this genetic map, we identified two markers associated with nitrite tolerance. We then sequenced the transcriptomes of the most nitrite-tolerant and nitrite-sensitive individuals from each of four genetically distinct L. vannamei families (LV-I–4). We found 2,002, 1,983, 1,954, and 1,867 differentially expressed genes in families LV-1, LV-2, LV-3, and LV-4, respectively. By integrating QTL and transcriptomics analyses, we identified a candidate gene associated with nitrite tolerance. This gene was annotated as solute carrier family 26 member 6 (SLC26A6). RNA interference (RNAi) analysis demonstrated that SLC26A6 was critical for nitrite tolerance in L. vannamei. The present study increases our understanding of the molecular mechanisms underlying nitrite tolerance in shrimp and provides a basis for molecular-marker-assisted shrimp breeding.

Identification of a Growth-Associated Single Nucleotide Polymorphism (SNP) in Cyclin C of the Giant Tiger Shrimp Penaeus monodon

The full-length cDNA of cyclin C of the giant tiger shrimp Penaeus monodon (PmCyC) was isolated by RACE-PCR. It was 1443 bp in length containing an open reading frame (ORF) of 804 bp and 267 deduced amino acids. Tissue distribution analysis indicated that PmCyC was more abundantly expressed in ovaries and testes than other tissues of female and male juveniles (P < 0.05). A pair of primers was designed, and an amplification product of 403 bp containing an intron of 123 bp was obtained. Polymorphism of amplified PmCyC gene segments of the 5th (3-month-old G5, N = 30) and 7th (5-month-old G7, N = 18) generations of domesticated juveniles was analyzed. Four conserved SNPs (T>C₁₃₄, T>C₁₈₈, G>A₃₇₉, and T>C₃₈₂) were found within the examined sequences. A TaqMan genotyping assay was developed for detection of a T>C₁₃₄ SNP. Association analysis indicated that this SNP displayed significant association with body weight (P < 4.2e-10) and total length (P < 2e-09) of the examined G7 P. monodon (N = 419) with an allele substitution effect of 5.02 ± 0.78 g and 1.41 ± 0.19 cm, respectively. Juveniles with C/C₁₃₄ (22.80 ± 2.51 g and 12.97 ± 0.53 cm, N = 19) and T/C₁₃₄ (20.41 ± 0.93 g and 12.77 ± 0.21 cm, N = 129) genotypes exhibited a significantly greater average body weight and total length than those with a T/T₁₃₄ genotype (14.72 ± 0.53 g and 11.37 ± 0.13 cm, N = 271) (P < 0.05).

Antibiotic-induced alterations and repopulation dynamics of yellowtail kingfish microbiota

BACKGROUND

The use of antibiotics in aquaculture is a common infection treatment and is increasing in some sectors and jurisdictions. While antibiotic treatment can negatively shift gut bacterial communities, recovery and examination of these communities in fish of commercial importance is not well documented. Examining the impacts of antibiotics on farmed fish microbiota is fundamental for improving our understanding and management of healthy farmed fish. This work assessed yellowtail kingfish (Seriola lalandi) skin and gut bacterial communities after an oral antibiotic combination therapy in poor performing fish that displayed signs of enteritis over an 18-day period. In an attempt to promote improved bacterial re-establishment after antibiotic treatment, faecal microbiota transplantation (FMT) was also administered via gavage or in the surrounding seawater, and its affect was evaluated over 15 days post-delivery.

RESULTS

Antibiotic treatment greatly perturbed the global gut bacterial communities of poor-performing fish - an effect that lasted for up to 18 days post treatment. This perturbation was marked by a significant decrease in species diversity and evenness, as well as a concomitant increase in particular taxa like an uncultured Mycoplasmataceae sp., which persisted and dominated antibiotic-treated fish for the entire 18-day period. The skin-associated bacterial communities were also perturbed by the antibiotic treatment, notably within the first 3 days; however, this was unlike the gut, as skin microbiota appeared to shift towards a more 'normal' (though disparate) state after 5 days post antibiotic treatment. FMT was only able to modulate the impacts of antibiotics in some individuals for a short time period, as the magnitude of change varied substantially between individuals. Some fish maintained certain transplanted gut taxa (i.e. present in the FMT inoculum; namely various Aliivibrio related ASVs) at Day 2 post FMT, although these were lost by Day 8 post FMT.

CONCLUSION

As we observed notable, prolonged perturbations induced by antibiotics on the gut bacterial assemblages, further work is required to better understand the processes/dynamics of their re-establishment following antibiotic exposure. In this regard, procedures like FMT represent a novel approach for promoting improved microbial recovery, although their efficacy and the factors that support their success requires further investigation.

High accuracy of pooled DNA genotyping by 2b-RAD sequencing in the Pacific white shrimp, Litopenaeus vannamei

Using pooled DNA genotyping to estimate the proportional contributions from multiple families in a pooled sample is of particular interest for selective breeding in aquaculture. We compared different pooled libraries with separate 2b-RAD sequencing of Litopenaeus vannamei individuals to assess the effect of different population structures (different numbers of individuals and families) on pooled DNA sequencing, the accuracy of parent sequencing of the DNA pools and the effect of SNP numbers on pooled DNA sequencing. We demonstrated that small pooled DNA genotyping of up to 53 individuals by 2b-RAD sequencing could provide a highly accurate assessment of population allele frequencies. The accuracy increased as the number of individuals and families increased. The allele frequencies of the parents from each pool were highly correlated with those of the pools or the corresponding individuals in the pool. We chose 500-28,000 SNPs to test the effect of SNP number on the accuracy of pooled sequencing, and no linear relationship was found between them. When the SNP number was fixed, increasing the number of individuals in the mixed pool resulted in higher accuracy of each pooled genotyping. Our data confirmed that pooled DNA genotyping by 2b-RAD sequencing could achieve higher accuracy than that of individual-based genotyping. The results will provide important information for shrimp breeding programs.

Comparative Analysis of the Transcriptome and Distribution of Putative SNPs in Two Rainbow Trout (Oncorhynchus mykiss) Breeding Strains by Using Next-Generation Sequencing

Selective breeding can significantly improve the establishment of sustainable and profitable aquaculture fish farming. For rainbow trout (Oncorhynchus mykiss), one of the main aquaculture coldwater species in Europe, a variety of selected hatchery strains are commercially available. In this study, we investigated the genetic variation between the local Born strain, selected for survival, and the commercially available Silver Steelhead strain, selected for growth. We sequenced the transcriptome of six tissues (gills, head kidney, heart, liver, spleen, and white muscle) from eight healthy individuals per strain, using RNA-seq technology to identify strain-specific gene-expression patterns and single nucleotide polymorphisms (SNPs). In total, 1760 annotated genes were differentially expressed across all tissues. Pathway analysis assigned them to different gene networks. We also identified a set of SNPs, which are heterozygous for one of the two breeding strains: 1229 of which represent polymorphisms over all tissues and individuals. Our data indicate a strong genetic differentiation between Born and Silver Steelhead trout, despite the relatively short time of evolutionary separation of the two breeding strains. The results most likely reflect their specifically adapted genotypes and might contribute to the understanding of differences regarding their robustness toward high stress and pathogenic challenge described in former studies.

An osmolality/salinity-responsive enhancer 1 (OSRE1) in intron 1 promotes salinity induction of tilapia glutamine synthetase

Euryhaline tilapia (Oreochromis mossambicus) are fish that tolerate a wide salinity range from fresh water to > 3× seawater. Even though the physiological effector mechanisms of osmoregulation that maintain plasma homeostasis in fresh water and seawater fish are well known, the corresponding molecular mechanisms that control switching between hyper- (fresh water) and hypo-osmoregulation (seawater) remain mostly elusive. In this study we show that hyperosmotic induction of glutamine synthetase represents a prominent part of this switch. Proteomics analysis of the O. mossambicus OmB cell line revealed that glutamine synthetase is transcriptionally regulated by hyperosmolality. Therefore, the 5' regulatory sequence of O. mossambicus glutamine synthetase was investigated. Using an enhancer trapping assay, we discovered a novel osmosensitive mechanism by which intron 1 positively mediates glutamine synthetase transcription. Intron 1 includes a single, functional copy of an osmoresponsive element, osmolality/salinity-responsive enhancer 1 (OSRE1). Unlike for conventional enhancers, the hyperosmotic induction of glutamine synthetase by intron 1 is position dependent. But irrespective of intron 1 position, OSRE1 deletion from intron 1 abolishes hyperosmotic enhancer activity. These findings indicate that proper intron 1 positioning and the presence of an OSRE1 in intron 1 are required for precise enhancement of hyperosmotic glutamine synthetase expression.

Harnessing genomics to fast-track genetic improvement in aquaculture

Aquaculture is the fastest-growing farmed food sector and will soon become the primary source of fish and shellfish for human diets. In contrast to crop and livestock production, aquaculture production is derived from numerous, exceptionally diverse species that are typically in the early stages of domestication. Genetic improvement of production traits via well-designed, managed breeding programmes has great potential to help meet the rising seafood demand driven by human population growth. Supported by continuous advances in sequencing and bioinformatics, genomics is increasingly being applied across the broad range of aquaculture species and at all stages of the domestication process to optimize selective breeding. In the future, combining genomic selection with biotechnological innovations, such as genome editing and surrogate broodstock technologies, may further expedite genetic improvement in aquaculture.

Single-Nucleotide Polymorphisms (SNP) Mining and Their Effect on the Tridimensional Protein Structure Prediction in a Set of Immunity-Related Expressed Sequence Tags (EST) in Atlantic Salmon (Salmo salar)

Single-nucleotide polymorphisms (SNPs) are single genetic code variations considered one of the most common forms of nucleotide modifications. Such SNPs can be located in genes associated to immune response and, therefore, they may have direct implications over the phenotype of susceptibility to infections affecting the productive sector. In this study, a set of immune-related genes (cc motif chemokine 19 precursor [ccl19], integrin β2 (itβ2, also named cd18), glutathione transferase omega-1 [gsto-1], heat shock 70 KDa protein [hsp70], major histocompatibility complex class I [mhc-I]) were analyzed to identify SNPs by data mining. These genes were chosen based on their previously reported expression on infectious pancreatic necrosis virus (IPNV)-infected Atlantic salmon phenotype. The available EST sequences for these genes were obtained from the Unigene database. Twenty-eight SNPs were found in the genes evaluated and identified most of them as transition base changes. The effect of the SNPs located on the 5’-untranslated region (UTR) or 3’-UTR upon transcription factor binding sites and alternative splicing regulatory motifs was assessed and ranked with a low-medium predicted FASTSNP score risk. Synonymous SNPs were found on itβ2 (c.2275G > A), gsto-1 (c.558G > A), and hsp70 (c.1950C > T) with low FASTSNP predicted score risk. The difference in the relative synonymous codon usage (RSCU) value between the variant codons and the wild-type codon (ΔRSCU) showed one negative (hsp70 c.1950C > T) and two positive ΔRSCU values (itβ2 c.2275G > A; gsto-1 c.558G > A), suggesting that these synonymous SNPs (sSNPs) may be associated to differences in the local rate of elongation. Nonsynonymous SNPs (nsSNPs) in the gsto-1 translatable gene region were ranked, using SIFT and POLYPHEN web-tools, with the second highest (c.205A > G; c484T > C) and the highest (c.499T > C; c.769A > C) predicted score risk possible. Using homology modeling to predict the effect of these nonsynonymous SNPs, the most relevant nucleotide changes for gsto-1 were observed for the nsSNPs c.205A > G, c484T > C, and c.769A > C. Molecular dynamics was assessed to analyze if these GSTO-1 variants have significant differences in their conformational dynamics, suggesting these SNPs could have allosteric effects modulating its catalysis. Altogether, these results suggest that candidate SNPs identified may play a crucial potential role in the immune response of Atlantic salmon.

Improved genomic resources for the black tiger prawn (Penaeus monodon)

World production of farmed crustaceans was 7.8 million tons in 2016. While only making up approximately 10% of world aquaculture production, crustaceans are generally high-value species and can earn significant export income for producing countries. Viet Nam is a major seafood producing country earning USD 7.3 billion in 2016 in export income with shrimp as a major commodity. However, there is a general lack of genomic resources available for shrimp species, which is challenging to obtain due to the need to deal with large repetitive genomes, which characterize many decapod crustaceans. The first tiger prawn (P. monodon) genome assembly was assembled in 2016 using the standard Illumina PCR-based pair-end reads and a computationally-efficient but relatively suboptimal assembler, SOAPdenovo v2. As a result, the current P. monodon draft genome is highly fragmented (> 2 million scaffolds with N₅₀ length of <1000 bp), exhibiting only moderate genome completeness (< 35% BUSCO complete single-copy genes). We sought to improve upon the recently published P. monodon genome assembly and completeness by generating Illumina PCR-free pair-end sequencing reads to eliminate genomic gaps associated with PCR-bias and performing de novo assembly using the updated MaSurCA de novo assembler. Furthermore, we scaffolded the assembly with low coverage Nanopore long reads and several recently published deep Illumina transcriptome paired-end sequencing data, producing a final genome assembly of 1.6 Gbp (1,211,364 scaffolds; N₅₀ length of 1982 bp) with an Arthropod BUSCO completeness of 96.8%. Compared to the previously published P. monodon genome assembly from China (NCBI Accession Code: NIUS01), this represents an almost 20% increase in the overall BUSCO genome completeness that now consists of more than 90% of Arthropod BUSCO single-copy genes. The revised P. monodon genome assembly (NCBI Accession Code: VIGR01) will be a valuable resource to support ongoing functional genomics and molecular-based breeding studies in Vietnam.

Fine Mapping of the High-pH Tolerance and Growth Trait-Related Quantitative Trait Loci (QTLs) and Identification of the Candidate Genes in Pacific White Shrimp (Litopenaeus vannamei)

High-pH tolerance and growth are important traits for the shrimp culture industry in areas with saline-alkali water. In the present study, an F1 full-sib family of Pacific white shrimp (Litopenaeus vannamei) was generated with a new "semidirectional cross" method, and double-digest restriction site-associated DNA sequencing (ddRAD-Seq) technology was applied to genotype the 2 parents and 148 progenies. A total of 3567 high-quality markers were constructed for the genetic linkage map, and the total map length was 4161.555 centimorgans (cM), showing 48 linkage groups (LGs) with an average interlocus length of 1.167 cM. With a constrained logarithm of odds (LOD) score ≥ 2.50, 12 high-pH tolerance and 2 growth (body weight) QTLs were located. L. vannamei genomic scaffolds were used to assist with the detection of 21 stress- and 5 growth-related scaffold genes. According to the high-pH transcriptome data of our previous study, 6 candidate high-pH response genes were discovered, and 5 of these 6 genes were consistently expressed with the high-pH transcriptome data, validating the locations of the high-pH tolerance trait-related QTLs in this study. This paper is the first report of fine-mapping high-pH tolerance and growth (body weight) trait QTLs in one L. vannamei genetic map. Our results will further benefit marker-assisted selection work and might be useful for promoting genomic research on the shrimp L. vannamei.

Applications of Next-Generation Sequencing Technologies and Computational Tools in Molecular Evolution and Aquatic Animals Conservation Studies: A Short Review

Aquatic ecosystems that form major biodiversity hotspots are critically threatened due to environmental and anthropogenic stressors. We believe that, in this genomic era, computational methods can be applied to promote aquatic biodiversity conservation by addressing questions related to the evolutionary history of aquatic organisms at the molecular level. However, huge amounts of genomics data generated can only be discerned through the use of bioinformatics. Here, we examine the applications of next-generation sequencing technologies and bioinformatics tools to study the molecular evolution of aquatic animals and discuss the current challenges and future perspectives of using bioinformatics toward aquatic animal conservation efforts.

Transcriptome-wide analysis of immune responses in Eriocheir sinensis hemocytes after challenge with different microbial derivatives

Knowledge about how Eriocheir sinensis interacts with microorganisms in its ambient environment is still lacking. Using RNA-Seq, we determined the most conserved genes and pathways compared with other animals and detected highly-induced immune genes in E. sinensis hemocytes post-in vivo challenge with different microbial derivatives. In total, 33.2 million high-quality reads were generated and assembled into 177,679 contigs. Completeness assessment and functional annotation were performed. Lipopolysaccharide, peptidoglycan, and β-1, 3-glucan stimulation induced 373, 173, and 108 differentially expressed (DE) transcripts, respectively. GO terms such as 'G-protein-coupled receptor binding', 'negative regulation of mitogen-activated protein kinase activity', and 'positive regulation of blood circulation' were enriched in the DE transcripts. Quantitative real-time PCR validated the data for selected genes. Our data contribute to understanding the immune defense mechanism in E. sinensis and the development of the innate immune system, thereby providing insights into disease control and prevention in aquaculture.

Transcriptome and Gene Coexpression Network Analyses of Two Wild Populations Provides Insight into the High-Salinity Adaptation Mechanisms of Crassostrea ariakensis

Crassostrea ariakensis naturally distributes in the intertidal and estuary region with relative low salinity ranging from 10 to 25‰. To understand the adaptive capacity of oysters to salinity stress, we conducted transcriptome analysis to investigate the metabolic pathways of salinity stress effectors in oysters from two different geographical sites, namely at salinities of 16, 23, and 30‰. We completed transcriptome sequencing of 18 samples and a total of 52,392 unigenes were obtained after assembly. Differentially expressed gene (DEG) analysis and weighted gene correlation network analysis (WGCNA) were performed using RNA-Seq transcriptomic data from eye-spot larvae at different salinities and from different populations. The results showed that at moderately high salinities (23 and 30‰), genes related to osmotic agents, oxidation-reduction processes, and related regulatory networks of complex transcriptional regulation and signal transduction pathways dominated to counteract the salinity stress. Moreover, there were adaptive differences in salinity response mechanisms, especially at high salinity, in oyster larvae from different populations. These results provide a framework for understanding the interactions of multiple pathways at the system level and for elucidating the complex cellular processes involved in responding to osmotic stress and maintaining growth. Furthermore, the results facilitate further research into the biological processes underlying physiological adaptations to hypertonic stress in marine invertebrates and provide a molecular basis for our subsequent search for high salinity-tolerant populations.

The Single-molecule long-read sequencing of Scylla paramamosain

Scylla paramamosain is an important aquaculture crab, which has great economical and nutritional value. To the best of our knowledge, few full-length crab transcriptomes are available. In this study, a library composed of 12 different tissues including gill, hepatopancreas, muscle, cerebral ganglion, eyestalk, thoracic ganglia, intestine, heart, testis, ovary, sperm reservoir, and hemocyte was constructed and sequenced using Pacific Biosciences single-molecule real-time (SMRT) long-read sequencing technology. A total of 284803 full-length non-chimeric reads were obtained, from which 79005 high-quality unique transcripts were obtained after error correction and sequence clustering and redundant. Additionally, a total of 52544 transcripts were annotated against protein database (NCBI nonredundant, Swiss-Prot, KOG, and KEGG database). A total of 23644 long non-coding RNAs (lncRNAs) and 131561 simple sequence repeats (SSRs) were identified. Meanwhile, the isoforms of many genes were also identified in this study. Our study provides a rich set of full-length cDNA sequences for S. paramamosain, which will greatly facilitate S. paramamosain research.