High-resolution genetic linkage mapping, high-temperature tolerance and growth-related quantitative trait locus (QTL) identification in Marsupenaeus japonicus

Penaeid Shrimp Chromosome Studies Entering the Post-Genomic Era

Chromosome studies provide the foundation for comprehending inheritance, variation, systematics, and evolution. Penaeid shrimps are a group of crustaceans with great economic importance. Basic cytogenetic information obtained from these shrimps can be used to study their genome structure, chromosome relationships, chromosome variation, polyploidy manipulation, and breeding. The study of shrimp chromosomes experienced significant growth in the 1990s and has been closely linked to the progress of genome research since the application of next-generation sequencing technology. To date, the genome sequences of five penaeid shrimp species have been published. The availability of these genomes has ushered the study of shrimp chromosomes into the post-genomic era. Currently, research on shrimp cytogenetics not only involves chromosome counting and karyotyping, but also extends to investigating submicroscopic changes; exploring genome structure and regulation during various cell divisions; and contributing to the understanding of mechanisms related to growth, sexual control, stress resistance, and genome evolution. In this article, we provide an overview of the progress made in chromosome research on penaeid shrimp. We emphasize the mutual promotion between studies on chromosome structure and genome research and highlight the impact of chromosome-level assembly on studies of genome structure and function. Additionally, we summarize the emerging trends in post-genomic-era shrimp chromosome research.

Construction of a high-density linkage map and detection of sex-specific markers in Penaeus japonicus

Penaeus japonicus is one of the most important farmed shrimp species in many countries. Sexual dimorphism is observed in P. japonicus, in which females grow faster and larger than males; therefore, a unisexual female culture of P. japonicus could improve the efficiency of productivity. However, the genetic mechanisms underlying sex determination in P. japonicus are unclear. In this study, we constructed a high-density genetic linkage map of P. japonicus using genotyping-by-sequencing (GBS) technology in a full-sib family. The final map was 3,481.98 cM in length and contained 29,757 single nucleotide polymorphisms (SNPs). These SNPs were distributed on 41 sex-averaged linkage groups, with an average inter-marker distance of 0.123 cM. One haplotype, harboring five sex-specific SNPs, was detected in linkage group 1 (LG1), and its corresponding confidence interval ranged from 211.840 to 212.592 cM. Therefore, this high-density genetic linkage map will be informative for genome assembly and marker-assisted breeding, and the sex-linked SNPs will be helpful for further studies on molecular mechanisms of sex determination and unisexual culture of P. japonicus in the future.

The Polymorphism of LvMMD2 and Its Association with Growth Traits in Litopenaeus vannamei

The Pacific white shrimp Litopenaeus vannamei is one of the major economic aquaculture species. The growth trait is considered as the most important trait in L. vannamei aquaculture. Identification of the genetic components underlying growth-related traits in L. vannamei could be useful for the selective breeding of growth trait. Our previous work identified several growth-related SNPs by genome-wide association study (GWAS). Based on the assembled genome, we identified a new candidate gene (LvMMD2) beside the associated marker. This gene encodes the progestin and AdipoQ receptor 10 (PAQR10) protein. We further investigate the polymorphisms of LvMMD2 and their association with body weight of L. vannamei. By resequencing the coding region of LvMMD2, a total of 8 SNPs were identified, including 6 synonymous mutations and 2 nonsynonymous mutations. Association analyses based on a population of 322 individuals revealed that several SNPs located in the coding region of LvMMD2 were significantly associated with the body weight, especially the nonsynonymous mutation named as MMD_5 contributed the most association to the trait and it could explain 10.5% of phenotypic variance. In addition, several genes involved in growth and development have been identified as LvMMD2-interacting genes. These findings strongly suggested that LvMMD2 might be an important gene regulating the shrimp growth. More importantly, the MMD_5 could be a promising candidate locus for marker-assisted selection (MAS) of the body weight in L. vannamei.

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.

Comparative physiology and aquaculture: Toward Omics-enabled improvement of aquatic animal health and sustainable production

Omics-technologies have revolutionized biomedical research over the past two decades, and are now poised to play a transformative role in aquaculture. This article serves as an introduction to a Virtual Special Issue of Comparative Biochemistry and Physiology - Part D: Genomics and Proteomics (CBPD), with the objective to showcase the state-of-the-science for Omics in aquaculture. In this editorial, we describe the role that Omics can play in aquaculture, and provide a synopsis for each of the Special Issue articles that use these technologies to improve aquaculture practices. Current genomic resources available for some aquaculture species are also described. The number of datasets is impressive for species such as Atlantic salmon and rainbow trout, totaling in the thousands (NCBI Gene Expression Omnibus and Sequence Read Archive). We present a conceptual framework that describes how Omics can be leveraged to understand complex responses of aquatic animals in culture for relevant physiological outcomes, such as fecundity, growth, and immunity. Lastly, knowledge gaps and new directions are identified to address current obstacles in aquaculture. Articles in this Special Issue on aquaculture in CBPD highlight the diversity and scope of Omics in aquaculture. As the technology becomes more cost-effective, it is anticipated that genomics, transcriptomics, proteomics, metabolomics and lipidomics will play increasingly important roles in stock diagnostics (e.g. genetics, health, performance). The timing is right, as global concerns are reaching critical levels over food availability/security and water restrictions for humankind.

A High-Density Genetic Linkage Map and QTL Mapping for Sex in Black Tiger Shrimp (Penaeus monodon)

The black tiger shrimp, Penaeus monodon, is important in both fishery and aquaculture and is the second-most widely cultured shrimp species in the world. However, the current strains cannot meet the market needs in various cultural environments, and the genome resources for P. monodon are still lacking. Restriction-site associated DNA sequencing (RADseq) has been widely used in genetic linkage map construction and in quantitative trait loci (QTL) mapping. We constructed a high-density genetic linkage map with RADseq in a full-sib family. This map contained 6524 single nucleotide polymorphisms (SNPs) and 2208 unique loci. The total length was 3275.4 cM, and the genetic distance was estimated to be 1.1 Mb/cM. The sex trait is a dichotomous phenotype, and the same interval was detected as a QTL using QTL mapping and genome-wide association analysis. The most significant locus explained 77.4% of the phenotype variance. The sex locus was speculated to be the same in this species based on the sequence alignments in Mozambique, India, and Hawaii populations. The constructed genetic linkage map provided a valuable resource for QTL mapping, genome assembly, and genome comparison for shrimp. The demonstrated common sex locus is a step closer to locating the underlying gene.

High-Density Genetic Linkage Maps Provide Novel Insights Into ZW/ZZ Sex Determination System and Growth Performance in Mud Crab (Scylla paramamosain)

Mud crab, Scylla paramamosain is one of the most important crustacean species in global aquaculture. To determine the genetic basis of sex and growth-related traits in S. paramamosain, a high-density genetic linkage map with 16,701 single nucleotide polymorphisms (SNPs) was constructed using SLAF-seq and a full-sib family. The consensus map has 49 linkage groups, spanning 5,996.66 cM with an average marker-interval of 0.81 cM. A total of 516 SNP markers, including 8 female-specific SNPs segregated in two quantitative trait loci (QTLs) for phenotypic sex were located on LG32. The presence of female-specific SNP markers only on female linkage map, their segregation patterns and lower female: male recombination rate strongly suggest the conformation of a ZW/ZZ sex determination system in S. paramamosain. The QTLs of most (90%) growth-related traits were found within a small interval (25.18–33.74 cM) on LG46, highlighting the potential involvement of LG46 in growth. Four markers on LG46 were significantly associated with 10–16 growth-related traits. BW was only associated with marker 3846. Based on the annotation of transcriptome data, 11 and 2 candidate genes were identified within the QTL regions of sex and growth-related traits, respectively. The newly constructed high-density genetic linkage map with sex-specific SNPs, and the identified QTLs of sex- and growth-related traits serve as a valuable genetic resource and solid foundation for marker-assisted selection and genetic improvement of crustaceans.

Genome survey, high-resolution genetic linkage map construction, growth-related quantitative trait locus (QTL) identification and gene location in Scylla paramamosain

Scylla paramamosain is one of the most economically important crabs in China. In this study, the first genome survey sequencing of this crab was performed, and the results revealed that the estimated genome size was 1.21 Gb with high heterozygosity (1.3%). Then, RAD technology was used to construct a high-resolution linkage map for this species. A total of 24,444 single nucleotide polymorphism (SNP) makers were grouped into 47 linkage groups. The total length of the linkage groups was 3087.53 cM with a markers interval of 0.92 cM. With the aid of transcriptome and genome scaffold data, 4,271 markers were linked to genes, including several important growth-related genes such as transforming growth factor-beta regulator I, immune related-gene C-type lectin and ecdysone pathway gene broad-complex-like protein. Further, 442 markers, representing 279 QTLs, associated with 24 traits were identified, and of these markers, 78 were linked to genes. Some interesting genes, such as dedicator of cytokinesis protein 3, tenascin-X and DNA helicase MCM8, were believed to have important relationship with specific traits and merit further exploration. The results of this study will accelerate the genetic improvement and genome sequencing analysis of the mud crab.

The State of "Omics" Research for Farmed Penaeids: Advances in Research and Impediments to Industry Utilization

Elucidating the underlying genetic drivers of production traits in agricultural and aquaculture species is critical to efforts to maximize farming efficiency. "Omics" based methods (i.e., transcriptomics, genomics, proteomics, and metabolomics) are increasingly being applied to gain unprecedented insight into the biology of many aquaculture species. While the culture of penaeid shrimp has increased markedly, the industry continues to be impeded in many regards by disease, reproductive dysfunction, and a poor understanding of production traits. Extensive effort has been, and continues to be, applied to develop critical genomic resources for many commercially important penaeids. However, the industry application of these genomic resources, and the translation of the knowledge derived from "omics" studies has not yet been completely realized. Integration between the multiple "omics" resources now available (i.e., genome assemblies, transcriptomes, linkage maps, optical maps, and proteomes) will prove critical to unlocking the full utility of these otherwise independently developed and isolated resources. Furthermore, emerging "omics" based techniques are now available to address longstanding issues with completing keystone genome assemblies (e.g., through long-read sequencing), and can provide cost-effective industrial scale genotyping tools (e.g., through low density SNP chips and genotype-by-sequencing) to undertake advanced selective breeding programs (i.e., genomic selection) and powerful genome-wide association studies. In particular, this review highlights the status, utility and suggested path forward for continued development, and improved use of "omics" resources in penaeid aquaculture.

Applications of genotyping by sequencing in aquaculture breeding and genetics

Selective breeding is increasingly recognized as a key component of sustainable production of aquaculture species. The uptake of genomic technology in aquaculture breeding has traditionally lagged behind terrestrial farmed animals. However, the rapid development and application of sequencing technologies has allowed aquaculture to narrow the gap, leading to substantial genomic resources for all major aquaculture species. While high-density single-nucleotide polymorphism (SNP) arrays for some species have been developed recently, direct genotyping by sequencing (GBS) techniques have underpinned many of the advances in aquaculture genetics and breeding to date. In particular, restriction-site associated DNA sequencing (RAD-Seq) and subsequent variations have been extensively applied to generate population-level SNP genotype data. These GBS techniques are not dependent on prior genomic information such as a reference genome assembly for the species of interest. As such, they have been widely utilized by researchers and companies focussing on nonmodel aquaculture species with relatively small research communities. Applications of RAD-Seq techniques have included generation of genetic linkage maps, performing genome-wide association studies, improvements of reference genome assemblies and, more recently, genomic selection for traits of interest to aquaculture like growth, sex determination or disease resistance. In this review, we briefly discuss the history of GBS, the nuances of the various GBS techniques, bioinformatics approaches and application of these techniques to various aquaculture species.

Large-scale SNP discovery and construction of a high-density genetic map of Colossoma macropomum through genotyping-by-sequencing

Colossoma macropomum, or tambaqui, is the largest native Characiform species found in the Amazon and Orinoco river basins, yet few resources for genetic studies and the genetic improvement of tambaqui exist. In this study, we identified a large number of single-nucleotide polymorphisms (SNPs) for tambaqui and constructed a high-resolution genetic linkage map from a full-sib family of 124 individuals and their parents using the genotyping by sequencing method. In all, 68,584 SNPs were initially identified using minimum minor allele frequency (MAF) of 5%. Filtering parameters were used to select high-quality markers for linkage analysis. We selected 7,734 SNPs for linkage mapping, resulting in 27 linkage groups with a minimum logarithm of odds (LOD) of 8 and maximum recombination fraction of 0.35. The final genetic map contains 7,192 successfully mapped markers that span a total of 2,811 cM, with an average marker interval of 0.39 cM. Comparative genomic analysis between tambaqui and zebrafish revealed variable levels of genomic conservation across the 27 linkage groups which allowed for functional SNP annotations. The large-scale SNP discovery obtained here, allowed us to build a high-density linkage map in tambaqui, which will be useful to enhance genetic studies that can be applied in breeding programs.

Development of a 690 K SNP array in catfish and its application for genetic mapping and validation of the reference genome sequence

Single nucleotide polymorphisms (SNPs) are capable of providing the highest level of genome coverage for genomic and genetic analysis because of their abundance and relatively even distribution in the genome. Such a capacity, however, cannot be achieved without an efficient genotyping platform such as SNP arrays. In this work, we developed a high-density SNP array with 690,662 unique SNPs (herein 690 K array) that were relatively evenly distributed across the entire genome, and covered 98.6% of the reference genome sequence. Here we also report linkage mapping using the 690 K array, which allowed mapping of over 250,000 SNPs on the linkage map, the highest marker density among all the constructed linkage maps. These markers were mapped to 29 linkage groups (LGs) with 30,591 unique marker positions. This linkage map anchored 1,602 scaffolds of the reference genome sequence to LGs, accounting for over 97% of the total genome assembly. A total of 1,007 previously unmapped scaffolds were placed to LGs, allowing validation and in few instances correction of the reference genome sequence assembly. This linkage map should serve as a valuable resource for various genetic and genomic analyses, especially for GWAS and QTL mapping for genes associated with economically important traits.