Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei)

Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei) | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research article Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei) Digang Zeng, Chunling Yang, Qiangyong Li, Weilin Zhu, Xiuli Chen, and 8 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.2.18942/v3 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 02 Dec, 2020 Read the published version in BMC Genomics → Version 3 posted 4 You are reading this latest preprint version Show more versions Abstract Background: Ammonia is one of the most common toxicological environment factors affecting shrimp health. Although ammonia tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available. Results: In this study, we constructed a high-density genetic map of the Pacific white shrimp ( Litopenaeus vannamei ) using specific length amplified fragment sequencing (SLAF-seq). The constructed genetic map contained 17,338 polymorphic markers spanning 44 linkage groups, with a total distance of 6,360.12 centimorgans (cM) and an average distance of 0.37 cM. Using this genetic map, we identified a quantitative trait locus (QTL) that explained 7.41–8.46% of the phenotypic variance in L. vannamei survival time under acute ammonia stress. We then sequenced the transcriptomes of the most ammonia-tolerant and the most ammonia-sensitive individuals from each of four genetically distinct L. vannamei families. We found that 7546 genes were differentially expressed between the ammonia-tolerant and ammonia-sensitive individuals. Using QTL analysis and the transcriptomes, we identified one candidate gene (annotated as an ATP synthase g subunit) associated with ammonia tolerance. Conclusions: In this study, we constructed a high-density genetic map of L. vannamei and identified a QTL for ammonia tolerance. By combining QTL and transcriptome analyses, we identified a candidate gene associated with ammonia tolerance. Our work provides the basis for future genetic studies focused on molecular marker-assisted selective breeding. Epigenetics & Genomics Aquaculture and Mariculture genetic map QTL transcriptomic ammonia tolerance Litopenaeus vannamei Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Background The Pacific white shrimp ( Litopenaeus vannamei ) is the most widely cultivated and highest-yielding crustacean species in the world [1]. L. vannamei tolerates a wide range of salinities, grows rapidly, is highly disease resistant, and can be farmed at high densities [2]. However, high-density shrimp cultivation often leads to water quality deterioration [3]. The toxicological factors associated with poor quality water often negatively affect shrimp [4]. One of the most common toxicological factors affecting shrimp health is ammonia nitrogen (ammonia-N) [5]. In aquaculture water, ammonia-N is mainly found as non-ionic ammonia (NH 3 ) and ionic ammonia (NH 4 + ); these compounds are usually in dynamic equilibrium [6]. As NH 3 has no electric charge, it is highly fat-soluble and can easily penetrate organismal cell membranes, leading to toxic effects [7]. In aquatic organisms, NH 3 affects membrane stability, as well as physiology, biochemistry, and growth; shrimp exposed to NH 3 may exhibit dyspnea, lack of appetite, decreased disease resistance, and even death [8-11]. The maintenance of low aquatic ammonia-N concentrations is required for successful shrimp farming [12]. However, aquatic physical and chemical properties are complicated, and may be affected by various factors such as weather, the local environment, and the introduction of artificial feeds. Thus, new breeds of ammonia-tolerant shrimp may improve industrial production and reduce economic losses. Marker-assisted selection has proven to be a useful strategy for the development of new breeds with dramatically improved trait characteristics [13], and the first step towards developing a new shrimp breed is to identify genes or markers associated with the desired trait [14]. Several previous studies have focused on the genetic bases of ammonia tolerance in shrimp. For example, Lu et al. identified 12 single nucleotide polymorphisms (SNPs) associated with ammonia tolerance in L. vannamei using marker-trait correlation analyses [15]. At the same time, Lu et al. identified 202 proteins that were significantly differentially expressed between ammonia-tolerant and ammonia-sensitive L. vannamei families using a comparative proteome analysis based on isobaric tags for relative and absolute quantification (iTRAQ) [16]. In addition, Jie et al. identified several pathways and genes involved in ammonia tolerance in L. vannamei based on comparative transcriptomic and metabolomic analyses of ammonia-tolerant and ammonia-sensitive L. vannamei families [17]. Finally, several studies identified transcriptomic changes and differentially expressed genes in L. vannamei after ammonia stress [9, 18]. However, no studies have investigated the quantitative trait loci (QTL) associated with ammonia tolerance in shrimp. QTL analysis effectively identifies molecular markers or candidate genes associated with economically important traits in plants and animals [19]. QTL analyses usually require high-density genetic linkage maps. To date, genetic linkage maps have primarily been constructed using high-throughput sequencing technologies, such as restriction site-related DNA sequencing (RAD-seq), genotyping sequencing (GBS), and specific length amplified fragment sequencing (SLAF-seq) [20]. In particular, SLAF-seq efficiently identifies and genotypes large-scale SNPs [20]. SLAF-seq has been applied to many plant species, including spinach [21], sesame [22], walnut [23] soybean [24], cucumber [25], wax gourd [26], cauliflower [27], white jute [28], and maize [29]. SLAF-seq has also been successfully applied to L. vannamei [30]. Therefore, SLAF-seq was used in the current study to construct a high-density genetic map of L. vannamei . Furthermore, QTL analysis of ammonia tolerance in L. vannamei was performed. Transcriptomic differences between ammonia-tolerant and ammonia-sensitive individuals across several L. vannamei families were compared to identify potential candidate genes coferring ammonia tolerance within QTLs. Methods Preparation of the mapping family The L. vannamei used in experiments were obtained from the shrimp-breeding center at the Guangxi Academy of Fishery Sciences (Nanning, Guangxi, China). The L. vannamei family used for mapping was constructed using artificial insemination. In brief, a male shrimp from a family with a relatively high ammonia-tolerance (obtained via 10 consecutive generations of breeding) was mated with a female shrimp from a common family. The hatched offspring were reared for about 1 year. Then, a male and female shrimp were randomly selected from the year-old offspring and mated. The F1 progeny were used for mapping (LV-N). Measurement of ammonia tolerance A total of 284 shrimp (average body weight: 20.78 g) were randomly selected from the LV-N family. Selected shrimp were transferred to a 2 m × 4 m × 1 m indoor pool and allowed to acclimate for one week. Aquatic conditions during the acclimation and experimental periods were kept constant: temperature of 27.0 ± 0.5°C, pH of 8.1 ± 0.2, salinity of 30.2‰, and dissolved oxygen of 6–8 mg/L; culture water was kept aerated, and shrimp were fed formulated pellets (Zhengda Feed, China) daily at a ratio of 5% body weight. Following acclimation, an acute ammonia stress test was performed. The ammonia-N concentration used for the acute stress test was 345.94 mg/L, based on the results of a preliminary experiment. This was the concentration at which half of the experimental shrimp died in 72 hours under stress. The ammonia-N concentration of the water in the experimental pool was controlled by adding NH 4 Cl stock solution (prepared by dissolving analytically pure NH 4 Cl in filtered seawater). The concentration of ammonia-N in the water was measured daily using standard methods [31]. To keep the ammonia-N concentration constant, NH 4 Cl stock solution was added if the ammonia-N concentration was 345.94 mg/L. During the experiment, shrimp heath was observed every hour, and dead shrimp were removed immediately. Shrimp were considered dead when lying motionless on the bottom of the pool and not responding to external stimuli. Collected dead shrimp were immediately frozen in liquid nitrogen and stored at −20°C for DNA extraction. The survival time of each shrimp was used as a proxy for ammonia tolerance. The experiment ended when all shrimp had died. DNA extraction DNA was collected from the 284 F1 (LV-N) shrimp and the two parent shrimp. Marine animal genomic DNA extraction kits (Tiangen Biotech, China) were used to extract DNA from the tail muscle of each shrimp. DNA was quantified using a NanoDrop spectrophotometer and 1% agarose gel electrophoresis with a lambda DNA standard. SLAF library preparation and sequencing First, we predicted the digestion of the L. vannamei genome (https://www.ncbi.nlm.nih.gov/genome/?term=Vannamei) [32] using self-developed software. We digested the extracted genomic DNA of all LV-N shrimp using the endonucleases identified by the predictive software. Then, dual-index sequencing adaptors were ligated to the DNA fragments obtained by digestion with T4 ligase, and the fragments were amplified using polymerase chain reactions (PCRs). PCR products (314–414 bp including the adaptor sequences) were purified and re-amplified using PCR. SLAF sequencing was carried out on an Illumina HiSeq system, following the Illumina-recommended procedure. To assess the accuracy of library construction, the same library-construction and sequencing steps using the genome of Oryza sativa japonica as a control was performed. Library construction and sequencing were performed by Biomarker Technologies Corporation (Beijing, China). SLAF-seq data analysis and genotyping The raw sequencing reads were quality controlled by removing reads with a quality score <20. The remaining raw reads were grouped by individual based on the dual-index adaptor sequences. The dual-index adaptor and 5-bp end sequences were then trimmed to obtain clean reads. The clean reads were mapped to the L. vannamei genome (https://www.ncbi.nlm.nih.gov/genome/?term=Vannamei) [32] using BWA [33]. Reads mapped to the same position with >95% identity were considered the same SLAF. SNP-based polymorphic SLAF markers were identified by aligning reads from the same SLAF sequence. These polymorphic SLAF markers were then filtered by removing those with a parental sequencing depth less than 10-fold; those where the number of SNPs was >5; those where the proportion of genotypes covering offspring was <70%; and those with significant segregation distortion (chi-square test P<0.05). The remaining polymorphic SLAFs were classified into eight separate patterns: aa × bb, ab × cd, cc × ab, ab × cc, ef × eg, hk × hk, nn × np, and lm × ll. Because the mapping population used in this study was an F1 population, the polymorphic SLAF with the pattern aa × bb was removed, and the remaining polymorphic SLAFs were used for the construction of the genetic map. Genetic map construction and QTL analysis After coding the genotypes of the polymorphic SLAF markers, the genetic map was constructed using the single-chain clustering algorithm in HighMap [34], with the probability log threshold set to ≥5.0 and a maximum recombination rate of 0.4. The Kosanbi mapping function was used to convert percent recombination to genetic distance (cM). QTL analysis was conducted using the R/qtl software package [35]. The logarithm of odds (LOD) threshold was determined based on 1,000 permutations (P < 0.05). The phenotypic variance explained by the QTL was estimated using the formula 1 – 10 –2LOD/ n , where n was the sample size [36]. Transcriptome sequencing, candidate gene identification and quantitative real-time PCR (qRT-PCR) verification To identify differentially expressed genes (DEGs) between ammonia-tolerant and ammonia-sensitive L. vannamei , the transcriptomes of four L. vannamei families were sequenced: the mapping family (LV-N) and three other randomly chosen common families (LV-A, LV-C, and LV-F) with different genetic backgrounds. Our previous analysis indicated that the 24-h median lethal concentration of NH4Cl was 140.96 mg/L, 189.19 mg/L, 117.88 mg/L, and 137.26 mg/L for families LV-A, LV-C, LV-F and LV-N, respectively (Supplementary Material, Table S1). Two hundred shrimp from each family were randomly selected, and subjected to the acute ammonia stress test (345.94 mg/L ammonia-N), as described above. In each family, 20 shrimp with the longest survival times (i.e., the most ammonia tolerant) were collected, as were the 20 shrimp with the shortest survival times (i.e., the most ammonia sensitive). When collecting the ammonia-sensitive shrimp, specimens that were out of balance and lying on the bottom of the pool were judged to be dying, and were collected immediately, without waiting for death. The hepatopancreas of each shrimp was extracted, and hepatopancreases were pooled to form an ammonia-tolerant sample and an ammonia-sensitive sample per family. Total RNA was extracted from each pooled sample using TRIzol reagent (Invitrogen, USA), following the manufacturer's instructions. Residual genomic DNA was removed with DNase I. RNA purity (OD260 / 280), concentration, and absorption peak were measured using a NanoDrop 2000. RNA integrity was assessed using an RNA Nano 6000 Assay Kit with an Agilent Bioanalyzer 2100. The isolated mRNA was divided into 100–400 bp fragments using an RNA fragment reagent (Illumina, USA). cDNA libraries were then constructed using NEBNext Ultra RNA Library Prep Kits (Illumina, USA), following manufacturer’s recommendations, and sequenced on an Illumina HiSeq system (Illumina, USA). Library construction and sequencing were performed by Biomarker Technologies Corporation (Beijing, China). Raw sequencing reads were trimmed and filtered using in-house Perl scripts to remove adaptor sequences and low-quality reads; the Q20, Q30, GC-content, and sequence duplication levels of the clean data were calculated. Clean reads were then aligned to the L. vannamei genome (https://www.ncbi.nlm.nih.gov/genome/?term=Vannamei) [32] using Hisat2 2.1.0 (http://ccb.jhu.edu/software/hisat2/index.shtml) [37]. Matched reads were counted to determine gene expression levels using the fragments per kilobase of transcript per million mapped reads (FPKM) method [38]. DEGs were identified using edger [39]. unigenes were considered differentially expressed when the false discovery rate (FDR) was ≤ 0.01 and the fold change between groups was > 2. DEGs were functionally annotated against the following databases: Non-Redundant protein sequences (NR) (ftp://ftp.ncbi.nih.gov/blast/db/), Protein family (Pfam) (http://pfam.xfam.org/), Clusters of Orthologous Groups (http://www.ncbi.nlm.nih.gov/COG/), Swiss-Prot (http://www.uniprot.org/), Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/), and Gene Ontology (GO) (http://www.geneontology.org/). After obtaining DEGs, candidate genes among the DEGs were identified. We consider candidate genes associated with ammonia tolerance when (1) candidate genes located within the QTL interval; (2) candidate genes differentially expressed between the most ammonia-tolerant and the most ammonia-sensitive individuals in the mapping family (LV-N); (3) the regulation pattern (up- or down-regulated expression) of candidate genes between the most ammonia-tolerant and the most ammonia-sensitive individuals was consistent across the four families of Litopenaeus vannamei (LV-A, LV-C, LV-F, and LV-N). qRT-PCR was used to validate the RNA-seq results by quantifying the expression of the candidate gene (LOC113809108) in the ammonia-tolerant and ammonia-sensitive pooled samples from the four L. vannamei families (LV-A, LV-C, LV-F, and LV-N). RNA-seq and qRT-PCRs analyses were carried out using the same samples. qRT-PCRs were performed using SYBR Premix Ex TaqTM II kits (TaKaRa, Japan), according to the manufacturer's instructions. The primer sets used to detect LOC113809108 gene expression levels were designed using the Primer Premier software (version 5.0) [40] as follows: 5’-ACTTGGGTGCTGTAGCTCAA-3’ and 5’-CTCGACAGCAACCAGGGTAT-3’. L. vannamei 18S RNA was used as the internal reference gene; this gene was amplified using the primer sets as follows: 5’-GCCTGAGAAACGGCTACCACATC-3’ and 5’-GTAGTAGCGACGGGCGGTGTGT-3’ [41]. The qRT-PCR cycling program was as follows: preheating at 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. The qRT-PCR was carried out at 95 °C for 40 s, 95 °C for 5 s, and 62 °C for 30 s for 40 cycles. Three parallel qRT-PCRs were carried out for each sample. Relative gene expression levels were calculated using the 2-ΔΔCT method[42]. Results Phenotypic variation We developed an ammonia-tolerant shrimp family (designated LV-N) for mapping, and subjected 284 LV-N shrimp to an acute ammonia stress test. All shrimp died within 2–98 hours, with a mean survival time of 65 hours. Individual survival times were normally distributed and thus suitable for QTL detection. The accumulated mortality rate of shrimp is showed in Supplementary Material Fig S1. SLAF-seq and genotyping Based on the digestive enzyme prediction using the reference genome of L. vannamei , HaeIII and Hpy166II were used to digest the genomic DNA of the 284 LV-N shrimp for SLAF library construction. SLAF sequencing generated 439.77 gigabases (Gb) of data, consisting of 2,201 megabases (Mb) of 100-bp paired-end reads. Across all reads, the average Q30 was 95.81%, the average GC content was 40.60%, and the GC distribution was normal (Table 1). The rice ( Oryza sativa japonica ) genome was used as a control to estimate the validity of the library construction. For the rice library, 343.21 Mb of data (1.72 Mb paired-end reads) were generated, with a Q30 of 95.81% and a GC content of 40.96%. In L. vannamei , 57.83% of the paired-end reads mapped successfully to the genome, as compared to 91.43% of the paired-end reads in rice. In addition, enzymatic digestion efficiency was 87.75% for L. vannamei and 92.19% for rice (Supplementary Material, Table S2). These results indicated that SLAF library construction and sequencing were adequate. After filtering and clustering all reads, 807,505 SLAFs were identified. The average sequencing depth of these SLAFs was 42.8-fold for the male parent, 42.14-fold for the female parent, and 12.43-fold for the progeny (Table 1). Of the 807,505 high-quality SLAFs detected, 293,415 (36.34%) were polymorphic (Table 1). After further filtering, the remaining 115,973 SLAF markers were successfully classified into eight genotypic patterns: ab × cd, cc × ab, aa × bb, ab × cc, ef × eg, lm × ll, hk × hk, and nn × np. The most common pattern was aa × bb, followed by nn × np and lm × ll (Fig. 1). Because the mapped population was an F1 population, aa × bb were eliminated as a valid marker. Characteristics of the genetic map Linkage analysis labeled 17,338 SLAF markers on the genetic map: 11,512 on the male map, 10,293 on the female map, and 17,338 on the sex-average map (Fig. 2). Each map contained 44 linkage groups (LGs). The total distances on the male, female, and sex-average maps were 6,604.99 cM, 5,476.20 cM, and 6,360.12 cM, respectively. The mean distance between adjacent markers was 0.58 cM on the male map, 0.53 cM on the female map, and 0.37 cM on the sex-average map (Supplementary Material, Table S3, Table S4, and Table S5). The distribution of markers among LGs was not uniform: in the male map, LG31 contained the most markers (585), while LG26 contained the least (39); in the female map, LG36 contained the most markers (540), while LG27 contained the least (21); and in the sex-average map, LG31 contained the most markers (695), while LG26 contained the least (53). QTL mapping of ammonia-tolerance A QTL analysis of the ammonia-tolerance trait in the LV-N L. vannamei family was performed based on the genetic maps. The LOD threshold was 4.75 (1000 permutations, P 4.75 were considered effective QTLs. Using this criterion, we identified a QTL within LG19 for ammonia tolerance (Fig. 3). The phenotypic variation explained by this QTL was 7.41–8.46%, the LOD score was 4.75–5.45, and the confidence interval was 12.42–29.43 cM. Transcriptome sequencing, candidate gene identification and qRT-PCR verification The transcriptomes of the 20 most ammonia-tolerant and the 20 most ammonia-sensitive shrimp in each of four L. vannamei families (LV-N, LV-A, LV-C, and LV-F) with various genetic backgrounds were sequenced. Transcriptome sequencing generated 56.79 Gb of clean data. A total of 7546 DEGs were identified between the ammonia-tolerant and ammonia-sensitive shrimp across all four families: 1,869 in LV-A, 2,005 in LV-C, 1,875 in LV-F, and 1,797 in LV-N (Supplementary Material, Table S6). The numbers of DEG annotations recovered in the databases searched were similar across the four L. vannamei families. For instance, the COG terms mainly enriched in the DEGs from all four families were posttranslational modification, protein turnover, chaperones, and general function prediction only (Fig. 4); the GO terms primarily enriched in the DEGs from all four families were binding, catalytic activity, cellular process, metabolic process, cell, cell part, single-organism process and membrane functions (Fig. 5). By aligning the DEGs with the QTL region in LG19, 107 DEGs located in the QTL interval were identified. The expression levels and annotations of these DEGs are listed in Supplementary Material Table S7. Of these DEGs, only one gene (LOC113809108) met the criterion used to determine candidate genes associated with ammonia tolerance. This gene was annotated as an ATP synthase g subunit. LOC113809108 was located in the QTL interval, and was significantly upregulated in the most ammonia-tolerant shrimp compared to the most ammonia-sensitive shrimp from families LV-N and LV-C (Fig. 6). This gene was also upregulated in the most ammonia-tolerant shrimp from families LV-A and LV-F, but this difference in expression was not significant (Table 2). The qRT-PCR analysis showed that the patterns of LOC113809108 gene expression in ammonia-tolerant and ammonia-sensitive pooled samples from the families LV-A, LV-C, LV-F, and LV-N were similar to the patterns determined using RNA-seq: LOC113809108 gene expression was upregulated in the ammonia-tolerant shrimp as compared to the ammonia-sensitive shrimp across all four families (Fig. 6). Discussion This study was aimed at investigating the ammonia tolerance in L.vannamei by QTL analysis. A high-density genetic map of L. vannamei was constructed using SLAF-seq, and a QTL associated with ammonia tolerance was identified as well as a putative candidate gene associated with ammonia tolerance. The genome of L. vannamei is large (~2.45 Gb) [32]. Whole-genome deep resequencing is relatively costly for large genomes and is often not necessary for gene/QTL mapping [43, 44]. In recent years, several genetic linkage maps based on SNPs were constructed and QTL analyses were conducted in L. vannamei [30, 45-47]. Yu et al. constructed a high-density genetic map for L. vannamei and detected several QTLs for body weight and body length [30]. Yang et al. mapped the sex determination region in L. vannamei based on the data used for high-density linkage map construction [45]. Du et al.mapped a QTL for L. vannamei gender using a gene-based SNP linkage map [46]. In this study, a high-density genetic map of L. vannamei was constructed using SLAF-seq, which is an effective method for discovering large numbers of SNPs and to perform large-scale genotyping [20]. Compared to traditional methods of genetic map construction (e.g., random amplified polymorphic DNA (RAPD), amplified fragment length polymorphism (AFLP), and simple sequence repeat (SSR)), the SLAF-seq method has several advantages for large-scale SNP discovery and genotyping: high density, high throughput, high efficiency, and low cost [48]. Previous studies have developed genetic maps of L. vannamei using RAPD, AFLP, and SSR, but in these maps, the average distance between adjacent markers was 1–5 cM [49-52]. The average distance between adjacent markers in the SLAF-seq genetic maps of L. vannamei in this study was substantially shorter (0.34 cM). Notably, the average distance between adjacent markers found here was also less than in a previously reported SLAF-seq genetic map of L. vannamei (0.75 cM) [30], possibly because a larger sample size and a greater sequencing depth were used. However, the number of LGs in the genetic map of L. vannamei in this study (44) was consistent with the number of LGs in the previously reported genetic map of L. vannamei [30]. This indicated that L. vannamei had 44 chromosomes, which was first reported by CamposRamos [53]. As far as we know, this is the highest density genetic linkage map L. vannamei , and it is very useful for comparative analysis of genomic synteny, QTL mapping, positioning of candidate genes, and marker-assisted selection. As most animals cannot self-fertilize, it is difficult to develop common populations for genetic mapping (e.g., F2, recombinant inbred line (RIL), and nearlyisogenic line (NIL) populations). Therefore, an F1 population of L. vannamei was used to construct the genetic map, relying on a pseudo-testcross strategy. This strategy was based on the selection of single-dose markers present in one parent and absent in the other, and carried at a 1:1 ratio by the F1 offspring [54]. Therefore, gamete separation in each individual can be directly analyzed. The pseudo-testcross strategy has been widely used to construct animal F1 populations for genetic mapping [55-58]. In this study, an F1 population was developed using one ammonia-tolerant male parent (the result of 10 generations of selective breeding) and one female shrimp from a common family. Previous studies have suggested that the size of the mapped population might affect the accuracy of the genetic map and the QTL analysis, and have shown that genetic map accuracy increases with the size of the population used [59]. Specifically, populations of >200 individuals are considered sufficient for the construction of accurate genetic maps [59]. Thus, 284 randomly-selected individuals from the F1 population were used to construct the linkage map. However, the determination of shrimp survival time during the ammonia stress experiment depended on human observation, and thus may not have been perfectly accurate. To reduce the possible impacts of measurement inaccuracies on the QTL analysis, a relatively large F1 population was used. This larger population increased the accuracy of the QTL mapping, compensating for any instances of human error in survival time measurement. Ammonia stress is one of the biggest challenges facing shrimp aquaculture. Ammonia not only has a direct lethal effect on shrimp [7], but also inhibits the shrimp immune system and increases sensitivity to pathogens [60]. Breeding new varieties of ammonia-tolerant shrimp is therefore an important target of the shrimp aquaculture industry. Some researches on shrimp ammonia tolerance have been carried out in recent years. Lu et al. found 12 SNPs associated with ammonia tolerance in L. vannamei using marker-trait correlation analyses, and these SNPs were identified in the thrombospondin gene and X-box binding protein 1 gene [15]. At the same time, Lu et al. identified 202 differentially expressed proteins (DEPs) between ammonia-tolerant and ammonia-sensitive L. vannamei using a comparative proteome analysis based on iTRAQ technique, and 77.8% of the DEPs were reported mainly involving in immune defense and stress tolerant in crustacean species [16]. Here, a QTL for ammonia-tolerance, located on LG19 at 169.09–169.49 cM was identified, that explained 7.41–8.46% of the phenotypic variation in ammonia tolerance. To the best of our knowledge, this is the first QTL for ammonia-tolerance reported in shrimp. However, having one QTL that explains about 8% of the variance indeed is not enough to initiate genetic breeding. Therefore in future research, the cutoff of LOD score can be relaxed to get more QTLs, intervals, and candidate DEGs. And then, more candidate genes related to ammonia tolerance should be identified. Moreover, independent families should be used to verify whether the QTL is common across different populations. A QTL usually spans a large chromosomal region and may contain hundreds of genes. Therefore, in order to identify functional genes associated with ammonia tolerance, the transcriptomes of four L. vannamei families (LV-N, LV-A, LV-C, and LV-F) with different genetic backgrounds were sequenced. Combining QTL mapping and gene expression analysis, we identified a single DEG (LOC113809108), located in the QTL interval, that was annotated as an ATP synthase g subunit and was significantly upregulated in the ammonia-tolerant LV-N and LV-C shrimp. The ATP synthase g subunit is located in F0 portion of ATP synthase, which consists of a membrane-integrated portion (F0 complex) and a membrane-protruding portion (F1 complex) [61]. The ATP synthase g subunit is essential for the formation of the F0 complex [61]. In the previous studies [15,16], ATP synthase has not been identified as being related to shrimp ammonia tolerance. ATP synthase is a double-motor enzyme that is involved in ATP synthesis, ATP hydrolysis-dependent processes, and the regulation of the proton gradient across some membrane-dependent systems [62]. Several studies have shown that ammonia excretion in aquatic animals is associated with Na + /K + -ATPase, which is mainly located on the basolateral membrane of branchial cells; NH 4 + is excreted into the environment when K + is replaced by NH 4 + via the Na + /NH 4 + exchanger [63-65]. Indeed, a previous study suggested that in L. vannamei , high ammonia tolerance was mainly a result of improved ammonia excretion and detoxification, as well as an accelerated energy metabolism [17]. Therefore, we speculate that ATP synthesis might affect the ammonia tolerance of L. vannamei by regulating ATP synthesis and controlling cellular ammonia excretion. The results would provide useful information for further study of the molecular mechanisms of ammonia adaptive strategies in shrimps. Conclusions In this study, we constructed a high-density genetic map of L. vannamei and identified a QTL for ammonia tolerance. By combining QTL and transcriptome analyses, we identified a candidate gene associated with ammonia tolerance. The results help us better understand the molecular mechanism of ammonia tolerance in shrimp. Our work provides the basis for future genetic studies focused on molecular marker-assisted selective breeding. Abbreviations cM: centimorgans; SLAF-seq: specific length amplified fragment sequencing; QTL:quantitative trait locus; ammonia-N: ammonia nitrogen; SNPs: single nucleotide polymorphisms; PCRs: polymerase chain reactions; LOD: logarithm of odds; qRT-PCR: quantitative real-time PCR; FPKM: fragments per kilobase of transcript per million mapped reads; FDR: false discovery rate; DEGs: differentially expressed genes; Gb: gigabases; Mb: megabases; LGs: linkage groups Declarations Acknowledgments We thank the Guangxi Shrimp Breeding Engineering Technology Research Center for providing animal materials. Authors’ contributions YZ and QL conceived the study. DZ and CY wrote the manuscript and contributed in the bioinformatics analysis. QL, WZ, XC, and MP conducted the experiments and contributed for raw data analysis. XC, YL, HW, HL, and JL contributed in the bioinformatics analysis. All authors read and approved the final manuscript. Funding This work was supported by the scientific and technological innovation major project of guangxi (aa17204080-1), scientific and technological innovation major base of guangxi (aa17204088-1), guangxi scientific research and technology development project (ab16380189), the national modern agriculture industry technology system project (nycytxgxcxtd-14-01), and the national modern agroindustry technology research system (cars-48). Availability of data and materials All data generated or analysed during this study are included in this published article and its supplementary information files. Raw SLAF sequencing reads are deposited in NCBI database under the accession numbers PRJNA545592. Raw RNA-seq are deposited in NCBI database under the accession numbers SRR9822091, SRR9822090, SRR9822095, SRR9822094, SRR9822098, SRR9822099, SRR9822093, and SRR9822085. Ethics approval and consent to participate Not applicable Consent for publication Not applicable Competing interests The authors declare they have no competing interests. References Cruz P, Ibarra AM, Mejia-Ruiz H, Gaffney PM, Perez-Enriquez R: Genetic variability assessed by microsatellites in a breeding program of Pacific white shrimp ( Litopenaeus vannamei ). Mar Biotechnol (NY). 2004; 6(2):157-164. Javahery S, Noori A, Hoseinifar SH: Growth performance, immune response, and digestive enzyme activity in Pacific white shrimp, Penaeus vannamei Boone, 1931, fed dietary microbial lysozyme. Fish & shellfish immunology. 2019; 92:528-535. Alagappan KM, Deivasigamani B, Somasundaram ST, Kumaran S: Occurrence of Vibrio parahaemolyticus and its specific phages from shrimp ponds in east coast of India. Current microbiology. 2010; 61(4):235-240. Sun CB, Wang G, Chan SF: Effects of artificial infection of Litopenaeus vannamei by Micrococcus lysodeikticus and WSSV on the activity of immunity related enzymes. Fish & shellfish immunology. 2015; 46(2):778-786. Marazza D, Bornens P, Le Gal Y: Effect of ammonia on survival and adenylate energy charge in the shrimp Palaemonetes varians . Ecotoxicology and environmental safety. 1996; 34(2):103-108. Frias-Espericueta MG, Harfush-Melendez M, Paez-Osuna F: Effects of ammonia on mortality and feeding of postlarvae shrimp Litopenaeus vannamei . Bulletin of environmental contamination and toxicology. 2000; 65(1):98-103. Qiu L, Shi X, Yu S, Han Q, Diao X, Zhou H: Changes of Ammonia-Metabolizing Enzyme Activity and Gene Expression of Two Strains in Shrimp Litopenaeus vannamei Under Ammonia Stress. Frontiers in physiology. 2018; 9:211. Liang Z, Liu R, Zhao D, Wang L, Sun M, Wang M, Song L: Ammonia exposure induces oxidative stress, endoplasmic reticulum stress and apoptosis in hepatopancreas of pacific white shrimp ( Litopenaeus vannamei ). Fish & shellfish immunology. 2016; 54:523-528. Lu X, Kong J, Luan S, Dai P, Meng X, Cao B, Luo K: Transcriptome Analysis of the Hepatopancreas in the Pacific White Shrimp ( Litopenaeus vannamei ) under Acute Ammonia Stress. PloS one. 2016; 11(10):e0164396. Lu X, Luan S, Cao B, Meng X, Sui J, Dai P, Luo K, Shi X, Hao D, Han G et al : Estimation of genetic parameters and genotype-by-environment interactions related to acute ammonia stress in Pacific white shrimp ( Litopenaeus vannamei ) juveniles at two different salinity levels. PloS one. 2017; 12(3):e0173835. Duan Y, Zhang Y, Dong H, Wang Y, Zheng X, Zhang J: Effect of dietary Clostridium butyricum on growth, intestine health status and resistance to ammonia stress in Pacific white shrimp Litopenaeus vannamei . Fish & shellfish immunology. 2017; 65:25-33. Valencia-Castaneda G, Frias-Espericueta MG, Vanegas-Perez RC, Perez-Ramirez JA, Chavez-Sanchez MC, Paez-Osuna F: Acute Toxicity of Ammonia, Nitrite and Nitrate to Shrimp Litopenaeus vannamei Postlarvae in Low-Salinity Water. Bulletin of environmental contamination and toxicology. 2018; 101(2):229-234. Campbell NR, LaPatra SE, Overturf K, Towner R, Narum SR: Association mapping of disease resistance traits in rainbow trout using restriction site associated DNA sequencing. G3 (Bethesda). 2014; 4(12):2473-2481. Yu Y, Liu J, Li F, Zhang X, Zhang C, Xiang J: Gene set based association analyses for the WSSV resistance of Pacific white shrimp Litopenaeus vannamei . Scientific reports. 2017; 7:40549. Lu X, Kong J, Meng X, Cao B, Luo K, Dai P, Luan S: Identification of SNP markers associated with tolerance to ammonia toxicity by selective genotyping from de novo assembled transcriptome in Litopenaeus vannamei . Fish & shellfish immunology. 2018; 73:158-166. Lu X, Luan S, Dai P, Meng X, Cao B, Luo K, Kong J: iTRAQ-based comparative proteome analysis for molecular mechanism of defense against acute ammonia toxicity in Pacific White shrimp Litopenaeus vannamei . Fish & shellfish immunology. 2018; 74:52-61. Xiao J, Li QY, Tu JP, Chen XL, Chen XH, Liu QY, Liu H, Zhou XY, Zhao YZ, Wang HL: Stress response and tolerance mechanisms of ammonia exposure based on transcriptomics and metabolomics in Litopenaeus vannamei . Ecotoxicology and environmental safety. 2019; 180:491-500. Lu X, Luan S, Dai P, Luo K, Chen B, Cao B, Sun L, Yan Y, Kong J: Insights into the molecular basis of immunosuppression and increasing pathogen infection severity of ammonia toxicity by transcriptome analysis in pacific white shrimp Litopenaeus vannamei . Fish & shellfish immunology. 2019; 88:528-539. Zhu M, Zhao S: Candidate gene identification approach: progress and challenges. International journal of biological sciences. 2007; 3(7):420-427. Sun X, Liu D, Zhang X, Li W, Liu H, Hong W, Jiang C, Guan N, Ma C, Zeng H et al : SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PloS one. 2013; 8(3):e58700. Qian W, Fan G, Liu D, Zhang H, Wang X, Wu J, Xu Z: Construction of a high-density genetic map and the X/Y sex-determining gene mapping in spinach based on large-scale markers developed by specific-locus amplified fragment sequencing (SLAF-seq). BMC genomics. 2017; 18(1):276. Zhang Y, Wang L, Xin H, Li D, Ma C, Ding X, Hong W, Zhang X: Construction of a high-density genetic map for sesame based on large scale marker development by specific length amplified fragment (SLAF) sequencing. BMC plant biology. 2013; 13:141. Zhu Y, Yin Y, Yang K, Li J, Sang Y, Huang L, Fan S: Construction of a high-density genetic map using specific length amplified fragment markers and identification of a quantitative trait locus for anthracnose resistance in walnut ( Juglans regia L.). BMC genomics. 2015; 16:614. Li B, Tian L, Zhang J, Huang L, Han F, Yan S, Wang L, Zheng H, Sun J: Construction of a high-density genetic map based on large-scale markers developed by specific length amplified fragment sequencing (SLAF-seq) and its application to QTL analysis for isoflavone content in Glycine max . BMC genomics. 2014; 15:1086. Wei Q, Wang Y, Qin X, Zhang Y, Zhang Z, Wang J, Li J, Lou Q, Chen J: An SNP-based saturated genetic map and QTL analysis of fruit-related traits in cucumber using specific-length amplified fragment (SLAF) sequencing. BMC genomics. 2014; 15:1158. Jiang B, Liu W, Xie D, Peng Q, He X, Lin Y, Liang Z: High-density genetic map construction and gene mapping of pericarp color in wax gourd using specific-locus amplified fragment (SLAF) sequencing. BMC genomics. 2015; 16:1035. Zhao Z, Gu H, Sheng X, Yu H, Wang J, Huang L, Wang D: Genome-Wide Single-Nucleotide Polymorphisms Discovery and High-Density Genetic Map Construction in Cauliflower Using Specific-Locus Amplified Fragment Sequencing. Frontiers in plant science. 2016; 7:334. Tao A, Huang L, Wu G, Afshar RK, Qi J, Xu J, Fang P, Lin L, Zhang L, Lin P: High-density genetic map construction and QTLs identification for plant height in white jute ( Corchorus capsularis L.) using specific locus amplified fragment (SLAF) sequencing. BMC genomics. 2017; 18(1):355. Liu C, Zhou Q, Dong L, Wang H, Liu F, Weng J, Li X, Xie C: Genetic architecture of the maize kernel row number revealed by combining QTL mapping using a high-density genetic map and bulked segregant RNA sequencing. BMC genomics. 2016; 17(1):915. Yu Y, Zhang X, Yuan J, Li F, Chen X, Zhao Y, Huang L, Zheng H, Xiang J: Genome survey and high-density genetic map construction provide genomic and genetic resources for the Pacific White Shrimp Litopenaeus vannamei . Scientific reports. 2015; 5:15612. Gilcreas FW: Future of standard methods for the examination of water and wastewater. Health laboratory science. 1967; 4(3):137-141. Zhang X, Yuan J, Sun Y, Li S, Gao Y, Yu Y, Liu C, Wang Q, Lv X, Ma KY et al : Penaeid shrimp genome provides insights into benthic adaptation and frequent molting. Nature communications. 2019; 10(1):356. Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25(14):1754-1760. Liu D, Ma C, Hong W, Huang L, Liu M, Liu H, Zeng H, Deng D, Xin H, Song J et al : Construction and analysis of high-density linkage map using high-throughput sequencing data. PloS one. 2014; 9(6):e98855. Broman KW, Gatti DM, Simecek P, Furlotte NA, Prins P, Sen S, Yandell BS, Churchill GA: R/qtl2: Software for Mapping Quantitative Trait Loci with High-Dimensional Data and Multiparent Populations. Genetics. 2019; 211(2):495-502. Gardiner SE, Norelli JL, de Silva N, Fazio G, Peil A, Malnoy M, Horner M, Bowatte D, Carlisle C, Wiedow C et al : Putative resistance gene markers associated with quantitative trait loci for fire blight resistance in Malus 'Robusta 5' accessions. BMC genetics. 2012; 13:25. Kim D, Langmead B, Salzberg SL: HISAT: a fast spliced aligner with low memory requirements. Nature methods. 2015; 12(4):357-360. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature methods. 2008; 5(7):621-628. Robinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010; 26(1):139-140. Singh VK, Mangalam AK, Dwivedi S, Naik S: Primer premier: program for design of degenerate primers from a protein sequence. BioTechniques. 1998; 24(2):318-319. Zhang Q, Li F, Zhang X, Dong B, Zhang J, Xie Y, Xiang J: cDNA cloning, characterization and expression analysis of the antioxidant enzyme gene, catalase, of Chinese shrimp Fenneropenaeus chinensis . Fish & shellfish immunology. 2008; 24(5):584-591. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001; 25(4):402-408. Xu X, Chao J, Cheng X, Wang R, Sun B, Wang H, Luo S, Wu T, Li Y: Mapping of a Novel Race Specific Resistance Gene to Phytophthora Root Rot of Pepper ( Capsicum annuum ) Using Bulked Segregant Analysis Combined with Specific Length Amplified Fragment Sequencing Strategy. PloS one. 2016; 11(3):e0151401. Wang G, Chen B, Du H, Zhang F, Zhang H, Wang Y, He H, Geng S, Zhang X: Genetic mapping of anthocyanin accumulation-related genes in pepper fruits using a combination of SLAF-seq and BSA. PloS one. 2018; 13(9):e0204690. Yu Y, Zhang X, Yuan J, Wang Q, Li S, Huang H, Li F, Xiang J: Identification of Sex-determining Loci in Pacific White Shrimp Litopeneaus vannamei Using Linkage and Association Analysis. Marine Biotechnology. 2017; 19: 277–286. Du ZQ, Ciobanu DC, Onteru SK, Gorbach D, Mileham AJ, Jaramillo G, Rothschild MF: A gene-based SNP linkage map for pacific white shrimp, Litopenaeus vannamei . Anim Genet 2010, 41(3); 286-294. Du Z, Onteru S, Gorbach D, Rothschild M: A SNP Genetic Map for Pacific White Shrimp ( Litopenaeus vannamei ). Animal Industry Report. 2010; 656(1):35. Liang D, Chen M, Qi X, Xu Q, Zhou F, Chen X: QTL Mapping by SLAF-seq and Expression Analysis of Candidate Genes for Aphid Resistance in Cucumber. Frontiers in plant science. 2016; 7:1000. Garcia DK, Dhar AK, Alcivar-Warren A: Molecular analysis of a RAPD marker (B20) reveals two microsatellites and differential mRNA expression in Penaeus vannamei. Molecular marine biology and biotechnology. 1996; 5(1):71-83. Perez F, Ortiz J, Zhinaula M, Gonzabay C, Calderon J, Volckaert FA: Development of EST-SSR markers by data mining in three species of shrimp: Litopenaeus vannamei , Litopenaeus stylirostris , and Trachypenaeus birdy . Mar Biotechnol (NY). 2005; 7(5):554-569. Andriantahina F, Liu X, Huang H: Genetic map construction and quantitative trait locus (QTL) detection of growth-related traits in Litopenaeus vannamei for selective breeding applications. PloS one. 2013; 8(9):e75206. Zhang L, Yang C, Zhang Y, Li L, Zhang X, Zhang Q, Xiang J: A genetic linkage map of Pacific white shrimp ( Litopenaeus vannamei ): sex-linked microsatellite markers and high recombination rates. Genetica. 2007; 131(1):37-49. CamposRamos, R: Chromosome studies on the marine shrimps Penaeus vannamei and P-californiensis (Decapoda). Journal of Crustacean Biology. 1997; 17(4), 666–673. Grattapaglia D, Sederoff R: Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo-testcross: mapping strategy and RAPD markers. Genetics. 1994; 137(4):1121-1137. Grattapaglia D, Bertolucci FL, Sederoff RR: Genetic mapping of QTLs controlling vegetative propagation in Eucalyptus grandis and E. urophylla using a pseudo-testcross strategy and RAPD markers. TAG Theoretical and applied genetics Theoretische und angewandte Genetik. 1995; 90(7-8):933-947. Porceddu A, Albertini E, Barcaccia G, Falistocco E, Falcinelli M: Linkage mapping in apomictic and sexual Kentucky bluegrass ( Poa pratensis L.) genotypes using a two way pseudo-testcross strategy based on AFLP and SAMPL markers. TAG Theoretical and applied genetics Theoretische und angewandte Genetik. 2002; 104(2-3):273-280. Liao M, Zhang L, Yang G, Zhu M, Wang D, Wei Q, Zou G, Chen D: Development of silver carp ( Hypophthalmichthys molitrix ) and bighead carp ( Aristichthys nobilis ) genetic maps using microsatellite and AFLP markers and a pseudo-testcross strategy. Animal genetics. 2007; 38(4):364-370. Xie W, Zhang X, Cai H, Huang L, Peng Y, Ma X: Genetic maps of SSR and SRAP markers in diploid orchardgrass ( Dactylis glomerata L.) using the pseudo-testcross strategy. Genome. 2011; 54(3):212-221. Adésio FerreiraI MFdS, Luciano da Costa e SilvaI, Cosme Damião CruzI: Estimating the effects of population size and type on the accuracy of genetic maps. Genetics and Molecular Biology. 2006; 29(1):187-192. Jiang L, Feng J, Ying R, Yin F, Pei S, Lu J, Cao Y, Guo J, Li Z: Individual and combined effects of ammonia-N and sulfide on the immune function and intestinal microbiota of Pacific white shrimp Litopenaeus vannamei . Fish & shellfish immunology. 2019; 92:230-240. Fujikawa M, Sugawara K, Tanabe T, Yoshida M: Assembly of human mitochondrial ATP synthase through two separate intermediates, F1-c-ring and b-e-g complex. FEBS Letter. 2015; 589(19 Pt B):2707-2712. Hong S, Pedersen PL: ATP synthase and the actions of inhibitors utilized to study its roles in human health, disease, and other scientific areas. Microbiology and molecular biology reviews : MMBR. 2008; 72(4):590-641, Table of Contents. Wilkie MP: Mechanisms of ammonia excretion across fish gills. Comparative Biochemistry and Physiology Part A. 1997; 118A(1):39-50. Lucu C, Flik G: Na+-K+-ATPase and Na+/Ca2+ exchange activities in gills of hyperregulating Carcinus maenas . The American journal of physiology. 1999; 276(2):R490-499. Weihrauch D, Morris S, Towle DW: Ammonia excretion in aquatic and terrestrial crabs. The Journal of experimental biology. 2004; 207(Pt 26):4491-4504. Tables Table 1. Summary of the constructed genetic map of Litopenaeus vannamei. Map data Value Total bases 439.77 Gb Total reads 2,201.28 Mb Average Q30 95.81% Average GC 40.60% Enzyme digestion protocol HaeIII+Hpy166II Restriction fragment length 314-414 bp Percentage of reads matching the L. vannamei genome 57.83% Average enzymatic digestion efficiency 87.75% Predicted number of markers 339,517 Number of high-quality slafs 807,505 Number of polymorphic slafs 293,415 Number of SLAF markers on the map 17,338 Average depth in parents 208.90 × Average depth in offspring individual 38.55 × Number of linkage groups 44 Total distance of the map 6,360.12 cM Average distance of the map 0.37 cM Table 2. RNA-Seq analysis showing expression of gene LOC113809108 in the most ammonia-tolerant and the most ammonia-sensitive individuals across four families of Litopenaeus vannamei (LV-A, LV-C, LV-F, and LV-N). Family Gene LOC113809108 expression level Sensitive group (FPKM) Sensitive group (count) Tolerant group (FPKM) Tolerant group (count) FDR Log2-fold Change Regulated Difference LV-A 112.8410 3239 130.2663 2561 0.9589 0.2072 up Normal LV-C 63.2980 1651 148.3678 2444 0.0002 1.2289 up Significant LV-F 70.4520 1378 122.4527 2258 0.9402 0.7975 up Normal LV-N 114.9725 2659 329.0251 8025 0.0001 1.5169 up Significant Supplementary Files SupplementaryTableS1.docx Table S1. Median lethal concentration of NH4Cl for Litopenaeus vannamei families LV-A, LV-C, LV-F, and LV-N. SupplementaryTableS2.docx Table S2. Summary of the SLAF sequencing of Oryza sativa japonica. SupplementaryTableS3.docx Table S3. Basic information of the male map. SupplementaryTableS4.docx Table S4. Basic information of the female map. SupplementaryTableS5.docx Table S5. Basic information of the sex-average map SupplementaryTableS6.docx Table S6. Statistics of the number of differentially expressed genes (DEGs) between the most ammonia-sensitive and most ammonia-tolerant shrimp in the four experimental families. SupplementaryTableS7.xlsx Table S7. DEGs within the QTL interval and their expression levels and annotations figS1.tif Fig S1. The accumulated mortality rate of shrimp in LV-N family under an acute ammonia stress (ammonia-N concentration of 345.94 mg/L, temperature of 27.0 ± 0.5°C, pH of 8.1 ± 0.2, salinity of 30.2‰, and dissolved oxygen of 6–8 mg/L). Cite Share Download PDF Status: Published Journal Publication published 02 Dec, 2020 Read the published version in BMC Genomics → Version 3 posted Editorial decision: Accept 17 Nov, 2020 Editor assigned by journal 20 Sep, 2020 Submission checks completed at journal 19 Sep, 2020 Editor invited by journal 19 Sep, 2020 You are reading this latest preprint version Show more versions Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-9560","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research article","associatedPublications":[],"authors":[{"id":2806522,"identity":"dc252dd9-453e-4589-8109-67f9d7291ed0","order_by":0,"name":"Digang Zeng","email":"","orcid":"","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":false,"prefix":"","firstName":"Digang","middleName":"","lastName":"Zeng","suffix":""},{"id":2806523,"identity":"1cbb98d7-705a-4cef-ab82-e0c471b9e48b","order_by":1,"name":"Chunling Yang","email":"","orcid":"","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":false,"prefix":"","firstName":"Chunling","middleName":"","lastName":"Yang","suffix":""},{"id":2806524,"identity":"3d2c81ba-f9b9-445e-b0d2-589087aa44a7","order_by":2,"name":"Qiangyong Li","email":"","orcid":"","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":false,"prefix":"","firstName":"Qiangyong","middleName":"","lastName":"Li","suffix":""},{"id":2806525,"identity":"db58aeab-63b5-4e14-bdb5-8ac2337c5197","order_by":3,"name":"Weilin Zhu","email":"","orcid":"","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":false,"prefix":"","firstName":"Weilin","middleName":"","lastName":"Zhu","suffix":""},{"id":2806526,"identity":"25dca916-2d2a-42e0-af3b-af6862510d2a","order_by":4,"name":"Xiuli Chen","email":"","orcid":"","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":false,"prefix":"","firstName":"Xiuli","middleName":"","lastName":"Chen","suffix":""},{"id":2806527,"identity":"097ac984-0be8-4a9a-b685-1970fe075ff6","order_by":5,"name":"Min Peng","email":"","orcid":"","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":false,"prefix":"","firstName":"Min","middleName":"","lastName":"Peng","suffix":""},{"id":2806528,"identity":"70cb3d21-4ac0-4822-83bb-5b24ea3b4937","order_by":6,"name":"Xiaohan Chen","email":"","orcid":"","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":false,"prefix":"","firstName":"Xiaohan","middleName":"","lastName":"Chen","suffix":""},{"id":2806529,"identity":"0608678a-93ef-4cb1-bdd0-c5b63a939603","order_by":7,"name":"Yong Lin","email":"","orcid":"","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":false,"prefix":"","firstName":"Yong","middleName":"","lastName":"Lin","suffix":""},{"id":2806530,"identity":"07142c76-7801-410b-8d71-a6330d182eb3","order_by":8,"name":"Huanling Wang","email":"","orcid":"","institution":"Huazhong Agriculture University","correspondingAuthor":false,"prefix":"","firstName":"Huanling","middleName":"","lastName":"Wang","suffix":""},{"id":2806531,"identity":"4bb6fae0-8cf4-41f1-9033-7235911cf6b9","order_by":9,"name":"Hong Liu","email":"","orcid":"","institution":"Huazhong Agriculture University","correspondingAuthor":false,"prefix":"","firstName":"Hong","middleName":"","lastName":"Liu","suffix":""},{"id":2806532,"identity":"219ad624-060b-4800-98e6-cb04779e4a7f","order_by":10,"name":"Jingzhen Liang","email":"","orcid":"","institution":"GuangXi University","correspondingAuthor":false,"prefix":"","firstName":"Jingzhen","middleName":"","lastName":"Liang","suffix":""},{"id":2806533,"identity":"27804a1d-3305-4cd9-b4f0-e3ae7a57c0b2","order_by":11,"name":"Qingyun Liu","email":"","orcid":"","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":false,"prefix":"","firstName":"Qingyun","middleName":"","lastName":"Liu","suffix":""},{"id":2806534,"identity":"f99156b9-3589-4cae-9788-87fc5bd8e1b5","order_by":12,"name":"Yongzhen Zhao","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA/ElEQVRIiWNgGAWjYLACxgYgIQHEHxvA3MYDRGthnNkAphqI18LMC9bCwIBXi8Hx3sMvf+6wyZOf3fzwse0Omzrd9sNAW2psonFqOXMuzULyTFox45xjxsa5Z9IkzM4kArUcS8ttwKXlRo6ZgWHb4cRmiQQz6dy2wxJmB4BaGBsO49eS2PY/sU0i/Zu0JUjL+YcEtRg/ONh2ILFHIsdMmhGk5QYBWyTPnDFjbGxLTpwhkVNs2NuWJrntBtCWBDx+4TveY/zxZ5td4vwZ6Rsf/Gyz4Tc7n/7wwYcaG5xaFA4wsElgCifgUA4C8g0MzB/wyI+CUTAKRsEoYGAAAADhZ42j4fIpAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0002-9051-0383","institution":"Guangxi Academy of Fisheries Sciences","correspondingAuthor":true,"prefix":"","firstName":"Yongzhen","middleName":"","lastName":"Zhao","suffix":""}],"badges":[],"createdAt":"2019-12-12 12:55:17","currentVersionCode":3,"declarations":"","doi":"10.21203/rs.2.18942/v3","doiUrl":"https://doi.org/10.21203/rs.2.18942/v3","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s12864-020-07254-x","type":"published","date":"2020-12-02T15:00:46+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":2689471,"identity":"a35b7a4c-9629-414d-a150-778c999a761e","added_by":"auto","created_at":"2020-09-29 21:43:10","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":329449,"visible":true,"origin":"","legend":"Number of markers associated with each of the eight polymorphic specific length amplified fragment (SLAF) segregation patterns. ","description":"","filename":"fig1.png","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/fig1.png"},{"id":2689472,"identity":"44d83aec-258a-4f6d-8183-badb36f356ab","added_by":"auto","created_at":"2020-09-29 21:43:18","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":9775399,"visible":true,"origin":"","legend":"High-density linkage map of Litopenaeus vannamei showing genetic distances among specific length amplified fragment (SLAF) markers. Black bars represent SLAF markers. ","description":"","filename":"fig2.png","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/fig2.png"},{"id":2689475,"identity":"10a06d72-4ff3-4b76-b7f1-abbb8cf336b0","added_by":"auto","created_at":"2020-09-29 21:43:18","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":381374,"visible":true,"origin":"","legend":"Quantitative trait loci (QTL) for ammonia tolerance in Litopenaeus vannamei, showing the logarithm of odds (LOD) values of the linkage groups. The gray line indicates the LOD threshold (4.75; P = 0.05).","description":"","filename":"fig3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/fig3.jpg"},{"id":2689476,"identity":"02675521-a662-44e9-8bb5-a215cf146fe5","added_by":"auto","created_at":"2020-09-29 21:43:18","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":400258,"visible":true,"origin":"","legend":"Clusters of Orthologous Groups (COG) classifications of the putative functions of the differentially expressed genes between the most ammonia-tolerant and the most ammonia-sensitive individuals across four families of Litopenaeus vannamei (LV-A, LV-C, LV-F, and LV-N).","description":"","filename":"fig4.png","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/fig4.png"},{"id":2689478,"identity":"5f951dd5-2b5e-4753-8a50-5806e04b88e5","added_by":"auto","created_at":"2020-09-29 21:43:19","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":496834,"visible":true,"origin":"","legend":"Gene Ontology (GO) classifications of the putative functions of the differentially expressed genes between the most ammonia-tolerant and the most ammonia-sensitive individuals across four families of Litopenaeus vannamei (LV-A, LV-C, LV-F, and LV-N).","description":"","filename":"fig5.png","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/fig5.png"},{"id":2689480,"identity":"08c299ca-9628-4813-b52c-c76b778e5d19","added_by":"auto","created_at":"2020-09-29 21:43:19","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":363019,"visible":true,"origin":"","legend":"Expression of LOC113809108 gene from the transcriptomic analysis validated by qRT-PCR. Expression of LOC113809108 gene was detected in the most ammonia-tolerant and the most ammonia-sensitive individuals from four families of Litopenaeus vannamei (LV-A, LV-C, LV-F, and LV-N). Data were normalized to 18s rRNA as the reference and presented as a relative log2-fold change to validate the transcriptomic analysis results. Error bars show the standard deviation of three technical replicates.","description":"","filename":"fig6.png","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/fig6.png"},{"id":2689468,"identity":"3e3f59ba-9f67-4784-8121-1135987b217d","added_by":"auto","created_at":"2020-09-29 21:43:13","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":453169,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/b55bdb58-8e37-4be2-b812-4543e30420e0.pdf"},{"id":2689470,"identity":"147fcc0f-aae5-4cc9-a068-b79a38818f2f","added_by":"auto","created_at":"2020-09-29 21:43:09","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":12343,"visible":true,"origin":"","legend":"Table S1. Median lethal concentration of NH4Cl for Litopenaeus vannamei families LV-A, LV-C, LV-F, and LV-N.","description":"","filename":"SupplementaryTableS1.docx","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/SupplementaryTableS1.docx"},{"id":2689473,"identity":"7c201697-e45c-420e-82d0-c2f671c9dd0c","added_by":"auto","created_at":"2020-09-29 21:43:18","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":12113,"visible":true,"origin":"","legend":"Table S2. Summary of the SLAF sequencing of Oryza sativa japonica.","description":"","filename":"SupplementaryTableS2.docx","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/SupplementaryTableS2.docx"},{"id":2689474,"identity":"388988b7-27b5-413c-a097-bb4a5ed04fda","added_by":"auto","created_at":"2020-09-29 21:43:18","extension":"docx","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":16914,"visible":true,"origin":"","legend":"Table S3. Basic information of the male map.","description":"","filename":"SupplementaryTableS3.docx","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/SupplementaryTableS3.docx"},{"id":2689477,"identity":"840e5b62-79b5-4136-81d9-3e8074b268ee","added_by":"auto","created_at":"2020-09-29 21:43:18","extension":"docx","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":16590,"visible":true,"origin":"","legend":"Table S4. Basic information of the female map.","description":"","filename":"SupplementaryTableS4.docx","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/SupplementaryTableS4.docx"},{"id":2689479,"identity":"5b140fe7-4e0e-4c88-86c8-7ecdd35eaf6e","added_by":"auto","created_at":"2020-09-29 21:43:19","extension":"docx","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":16764,"visible":true,"origin":"","legend":"Table S5. Basic information of the sex-average map","description":"","filename":"SupplementaryTableS5.docx","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/SupplementaryTableS5.docx"},{"id":2689481,"identity":"88c2e73e-6bb0-4235-b52e-98b8688dc1da","added_by":"auto","created_at":"2020-09-29 21:43:19","extension":"docx","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":12039,"visible":true,"origin":"","legend":"Table S6. Statistics of the number of differentially expressed genes (DEGs) between the most ammonia-sensitive and most ammonia-tolerant shrimp in the four experimental families.","description":"","filename":"SupplementaryTableS6.docx","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/SupplementaryTableS6.docx"},{"id":2689482,"identity":"b6bb3f91-71a0-49ee-8615-db7d65971d07","added_by":"auto","created_at":"2020-09-29 21:43:19","extension":"xlsx","order_by":7,"title":"","display":"","copyAsset":false,"role":"supplement","size":51962,"visible":true,"origin":"","legend":"Table S7. DEGs within the QTL interval and their expression levels and annotations","description":"","filename":"SupplementaryTableS7.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/SupplementaryTableS7.xlsx"},{"id":2689483,"identity":"bb5be897-77b7-4b85-a1b9-51639de2cfbf","added_by":"auto","created_at":"2020-09-29 21:43:19","extension":"tif","order_by":8,"title":"","display":"","copyAsset":false,"role":"supplement","size":66928,"visible":true,"origin":"","legend":"Fig S1. The accumulated mortality rate of shrimp in LV-N family under an acute ammonia stress (ammonia-N concentration of 345.94 mg/L, temperature of 27.0 ± 0.5°C, pH of 8.1 ± 0.2, salinity of 30.2‰, and dissolved oxygen of 6–8 mg/L).","description":"","filename":"figS1.tif","url":"https://assets-eu.researchsquare.com/files/rs-9560/v3/figS1.tif"}],"financialInterests":"","formattedTitle":"Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei)","fulltext":[{"header":"Background","content":"\u003cp\u003eThe Pacific white shrimp (\u003cem\u003eLitopenaeus vannamei\u003c/em\u003e) is the most widely cultivated and highest-yielding crustacean species in the world [1]. \u003cem\u003eL. vannamei\u003c/em\u003e tolerates a wide range of salinities, grows rapidly, is highly disease resistant, and can be farmed at high densities [2]. However, high-density shrimp cultivation often leads to water quality deterioration [3]. The toxicological factors associated with poor quality water often negatively affect shrimp [4]. One of the most common toxicological factors affecting shrimp health is ammonia nitrogen (ammonia-N) [5]. In aquaculture water, ammonia-N is mainly found as non-ionic ammonia (NH\u003csub\u003e3\u003c/sub\u003e) and ionic ammonia (NH\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e+\u003c/sup\u003e); these compounds are usually in dynamic equilibrium [6]. As NH\u003csub\u003e3\u003c/sub\u003e has no electric charge, it is highly fat-soluble and can easily penetrate organismal cell membranes, leading to toxic effects [7]. In aquatic organisms, NH\u003csub\u003e3\u003c/sub\u003e affects membrane stability, as well as physiology, biochemistry, and growth; shrimp exposed to NH\u003csub\u003e3\u003c/sub\u003e may exhibit dyspnea, lack of appetite, decreased disease resistance, and even death [8-11].\u003c/p\u003e\n\u003cp\u003eThe maintenance of low aquatic ammonia-N concentrations is required for successful shrimp farming [12]. However, aquatic physical and chemical properties are complicated, and may be affected by various factors such as weather, the local environment, and the introduction of artificial feeds. Thus, new breeds of ammonia-tolerant shrimp may improve industrial production and reduce economic losses. Marker-assisted selection has proven to be a useful strategy for the development of new breeds with dramatically improved trait characteristics [13], and the first step towards developing a new shrimp breed is to identify genes or markers associated with the desired trait [14]. Several previous studies have focused on the genetic bases of ammonia tolerance in shrimp. For example, Lu et al. identified 12 single nucleotide polymorphisms (SNPs) associated with ammonia tolerance in \u003cem\u003eL. vannamei\u003c/em\u003e using marker-trait correlation analyses [15]. At the same time, Lu et al. identified 202 proteins that were significantly differentially expressed between ammonia-tolerant and ammonia-sensitive \u003cem\u003eL. vannamei\u003c/em\u003e families using a comparative proteome analysis based on isobaric tags for relative and absolute quantification (iTRAQ) [16]. In addition, Jie et al. identified several pathways and genes involved in ammonia tolerance in L. \u003cem\u003evannamei\u003c/em\u003e based on comparative transcriptomic and metabolomic analyses of ammonia-tolerant and ammonia-sensitive \u003cem\u003eL. vannamei\u003c/em\u003e families [17]. Finally, several studies identified transcriptomic changes and differentially expressed genes in \u003cem\u003eL. vannamei\u003c/em\u003e after ammonia stress [9, 18]. However, no studies have investigated the quantitative trait loci (QTL) associated with ammonia tolerance in shrimp.\u003c/p\u003e\n\u003cp\u003eQTL analysis effectively identifies molecular markers or candidate genes associated with economically important traits in plants and animals [19]. QTL analyses usually require high-density genetic linkage maps. To date, genetic linkage maps have primarily been constructed using high-throughput sequencing technologies, such as restriction site-related DNA sequencing (RAD-seq), genotyping sequencing (GBS), and specific length amplified fragment sequencing (SLAF-seq) [20]. In particular, SLAF-seq efficiently identifies and genotypes large-scale SNPs [20]. SLAF-seq has been applied to many plant species, including spinach [21], sesame [22], walnut [23] soybean [24], cucumber [25], wax gourd [26], cauliflower [27], white jute [28], and maize [29]. SLAF-seq has also been successfully applied to \u003cem\u003eL. vannamei\u003c/em\u003e [30].\u003c/p\u003e\n\u003cp\u003eTherefore, SLAF-seq was used in the current study to construct a high-density genetic map of \u003cem\u003eL. vannamei\u003c/em\u003e. Furthermore, QTL analysis of ammonia tolerance in \u003cem\u003eL. vannamei \u003c/em\u003ewas performed. Transcriptomic differences between ammonia-tolerant and ammonia-sensitive individuals across several \u003cem\u003eL.\u003c/em\u003e\u003cem\u003e vannamei\u003c/em\u003e families were compared to identify potential candidate genes coferring ammonia tolerance within QTLs.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003ePreparation of the mapping family\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe \u003cem\u003eL. vannamei\u003c/em\u003e used in experiments were obtained from the shrimp-breeding center at the Guangxi Academy of Fishery Sciences (Nanning, Guangxi, China). The \u003cem\u003eL. vannamei\u003c/em\u003e family used for mapping was constructed using artificial insemination. In brief, a male shrimp from a family with a relatively high ammonia-tolerance (obtained via 10 consecutive generations of breeding) was mated with a female shrimp from a common family. The hatched offspring were reared for about 1 year. Then, a male and female shrimp were randomly selected from the year-old offspring and mated. The F1 progeny were used for mapping (LV-N).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMeasurement of ammonia tolerance\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eA total of 284 shrimp (average body weight: 20.78 g) were randomly selected from the LV-N family. Selected shrimp were transferred to a 2 m \u0026times; 4 m \u0026times; 1 m indoor pool and allowed to acclimate for one week. Aquatic conditions during the acclimation and experimental periods were kept constant: temperature of 27.0 \u0026plusmn; 0.5\u0026deg;C, pH of 8.1 \u0026plusmn; 0.2, salinity of 30.2\u0026permil;, and dissolved oxygen of 6\u0026ndash;8 mg/L; culture water was kept aerated, and shrimp were fed formulated pellets (Zhengda Feed, China) daily at a ratio of 5% body weight. Following acclimation, an acute ammonia stress test was performed. The ammonia-N concentration used for the acute stress test was 345.94 mg/L, based on the results of a preliminary experiment. This was the concentration at which half of the experimental shrimp died in 72 hours under stress. The ammonia-N concentration of the water in the experimental pool was controlled by adding NH\u003csub\u003e4\u003c/sub\u003eCl stock solution (prepared by dissolving analytically pure NH\u003csub\u003e4\u003c/sub\u003eCl in filtered seawater). The concentration of ammonia-N in the water was measured daily using standard methods [31]. To keep the ammonia-N concentration constant, NH\u003csub\u003e4\u003c/sub\u003eCl stock solution was added if the ammonia-N concentration was \u0026lt;345.94 mg/L, and seawater was added if the ammonia-N concentration was \u0026gt;345.94 mg/L. During the experiment, shrimp heath was observed every hour, and dead shrimp were removed immediately. Shrimp were considered dead when lying motionless on the bottom of the pool and not responding to external stimuli. Collected dead shrimp were immediately frozen in liquid nitrogen and stored at \u0026minus;20\u0026deg;C for DNA extraction. The survival time of each shrimp was used as a proxy for ammonia tolerance. The experiment ended when all shrimp had died.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDNA extraction\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eDNA was collected from the 284 F1 (LV-N) shrimp and the two parent shrimp. Marine animal genomic DNA extraction kits (Tiangen Biotech, China) were used to extract DNA from the tail muscle of each shrimp. DNA was quantified using a NanoDrop spectrophotometer and 1% agarose gel electrophoresis with a lambda DNA standard.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003e\u0026nbsp;\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSLAF library preparation and sequencing\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFirst, we predicted the digestion of the \u003cem\u003eL. vannamei\u003c/em\u003e genome (https://www.ncbi.nlm.nih.gov/genome/?term=Vannamei) [32] using self-developed software. We digested the extracted genomic DNA of all LV-N shrimp using the endonucleases identified by the predictive software. Then, dual-index sequencing adaptors were ligated to the DNA fragments obtained by digestion with T4 ligase, and the fragments were amplified using polymerase chain reactions (PCRs). PCR products (314\u0026ndash;414 bp including the adaptor sequences) were purified and re-amplified using PCR. SLAF sequencing was carried out on an Illumina HiSeq system, following the Illumina-recommended procedure. To assess the accuracy of library construction, the same library-construction and sequencing steps using the genome of \u003cem\u003eOryza sativa\u003c/em\u003e \u003cem\u003ejaponica\u003c/em\u003e as a control was performed. Library construction and sequencing were performed by Biomarker Technologies Corporation (Beijing, China).\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSLAF-seq data analysis and genotyping\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe raw sequencing reads were quality controlled by removing reads with a quality score \u0026lt;20. The remaining raw reads were grouped by individual based on the dual-index adaptor sequences. The dual-index adaptor and 5-bp end sequences were then trimmed to obtain clean reads. The clean reads were mapped to the \u003cem\u003eL. \u003c/em\u003e\u003cem\u003evannamei\u003c/em\u003e genome (https://www.ncbi.nlm.nih.gov/genome/?term=Vannamei) [32] using BWA [33]. Reads mapped to the same position with \u0026gt;95% identity were considered the same SLAF. SNP-based polymorphic SLAF markers were identified by aligning reads from the same SLAF sequence. These polymorphic SLAF markers were then filtered by removing those with a parental sequencing depth less than 10-fold; those where the number of SNPs was \u0026gt;5; those where the proportion of genotypes covering offspring was \u0026lt;70%; and those with significant segregation distortion (chi-square test P\u0026lt;0.05). The remaining polymorphic SLAFs were classified into eight separate patterns: aa \u0026times; bb, ab \u0026times; cd, cc \u0026times; ab, ab \u0026times; cc, ef \u0026times; eg, hk \u0026times; hk, nn \u0026times; np, and lm \u0026times; ll. Because the mapping population used in this study was an F1 population, the polymorphic SLAF with the pattern aa \u0026times; bb was removed, and the remaining polymorphic SLAFs were used for the construction of the genetic map.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003e\u0026nbsp;\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGenetic map construction and QTL analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAfter coding the genotypes of the polymorphic SLAF markers, the genetic map was constructed using the single-chain clustering algorithm in HighMap [34], with the probability log threshold set to \u0026ge;5.0 and a maximum recombination rate of 0.4. The Kosanbi mapping function was used to convert percent recombination to genetic distance (cM). QTL analysis was conducted using the R/qtl software package [35]. The logarithm of odds (LOD) threshold was determined based on 1,000 permutations (P \u0026lt; 0.05). The phenotypic variance explained by the QTL was estimated using the formula 1\u0026thinsp;\u0026ndash;\u0026thinsp;10\u003csup\u003e\u0026ndash;2LOD/\u003cem\u003en\u003c/em\u003e\u003c/sup\u003e, where n was the sample size [36].\u003c/p\u003e\n\u003cp\u003e\u003cem\u003e\u0026nbsp;\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTranscriptome sequencing, candidate gene identification and quantitative real-time PCR (qRT-PCR) verification\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo identify differentially expressed genes (DEGs) between ammonia-tolerant and ammonia-sensitive \u003cem\u003eL. vannamei\u003c/em\u003e, the transcriptomes of four \u003cem\u003eL. vannamei\u003c/em\u003e families were sequenced: the mapping family (LV-N) and three other randomly chosen common families (LV-A, LV-C, and LV-F) with different genetic backgrounds. Our previous analysis indicated that the 24-h median lethal concentration of NH4Cl was 140.96 mg/L, 189.19 mg/L, 117.88 mg/L, and 137.26 mg/L for families LV-A, LV-C, LV-F and LV-N, respectively (Supplementary Material, Table S1). Two hundred shrimp from each family were randomly selected, and subjected to the acute ammonia stress test (345.94 mg/L ammonia-N), as described above. In each family, 20 shrimp with the longest survival times (i.e., the most ammonia tolerant) were collected, as were the 20 shrimp with the shortest survival times (i.e., the most ammonia sensitive). When collecting the ammonia-sensitive shrimp, specimens that were out of balance and lying on the bottom of the pool were judged to be dying, and were collected immediately, without waiting for death. The hepatopancreas of each shrimp was extracted, and hepatopancreases were pooled to form an ammonia-tolerant sample and an ammonia-sensitive sample per family.\u003c/p\u003e\n\u003cp\u003eTotal RNA was extracted from each pooled sample using TRIzol reagent (Invitrogen, USA), following the manufacturer's instructions. Residual genomic DNA was removed with DNase I. RNA purity (OD260 / 280), concentration, and absorption peak were measured using a NanoDrop 2000. RNA integrity was assessed using an RNA Nano 6000 Assay Kit with an Agilent Bioanalyzer 2100. The isolated mRNA was divided into 100\u0026ndash;400 bp fragments using an RNA fragment reagent (Illumina, USA). cDNA libraries were then constructed using NEBNext Ultra RNA Library Prep Kits (Illumina, USA), following manufacturer\u0026rsquo;s recommendations, and sequenced on an Illumina HiSeq system (Illumina, USA). Library construction and sequencing were performed by Biomarker Technologies Corporation (Beijing, China).\u003c/p\u003e\n\u003cp\u003eRaw sequencing reads were trimmed and filtered using in-house Perl scripts to remove adaptor sequences and low-quality reads; the Q20, Q30, GC-content, and sequence duplication levels of the clean data were calculated. Clean reads were then aligned to the \u003cem\u003eL. vannamei\u003c/em\u003e genome (https://www.ncbi.nlm.nih.gov/genome/?term=Vannamei) [32] using Hisat2 2.1.0 (http://ccb.jhu.edu/software/hisat2/index.shtml) [37]. Matched reads were counted to determine gene expression levels using the fragments per kilobase of transcript per million mapped reads (FPKM) method [38]. DEGs were identified using edger [39]. unigenes were considered differentially expressed when the false discovery rate (FDR) was \u0026le; 0.01 and the fold change between groups was \u0026gt; 2. DEGs were functionally annotated against the following databases: Non-Redundant protein sequences (NR) (ftp://ftp.ncbi.nih.gov/blast/db/), Protein family (Pfam) (http://pfam.xfam.org/), Clusters of Orthologous Groups (http://www.ncbi.nlm.nih.gov/COG/), Swiss-Prot (http://www.uniprot.org/), Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/), and Gene Ontology (GO) (http://www.geneontology.org/).\u003c/p\u003e\n\u003cp\u003eAfter obtaining DEGs, candidate genes among the DEGs were identified. We consider candidate genes associated with ammonia tolerance when (1) candidate genes located within the QTL interval; (2) candidate genes differentially expressed between the most ammonia-tolerant and the most ammonia-sensitive individuals in the mapping family (LV-N); (3) the regulation pattern (up- or down-regulated expression) of candidate genes between the most ammonia-tolerant and the most ammonia-sensitive individuals was consistent across the four families of Litopenaeus vannamei (LV-A, LV-C, LV-F, and LV-N).\u003c/p\u003e\n\u003cp\u003eqRT-PCR was used to validate the RNA-seq results by quantifying the expression of the candidate gene (LOC113809108) in the ammonia-tolerant and ammonia-sensitive pooled samples from the four \u003cem\u003eL. vannamei\u003c/em\u003e families (LV-A, LV-C, LV-F, and LV-N). RNA-seq and qRT-PCRs analyses were carried out using the same samples. qRT-PCRs were performed using SYBR Premix Ex TaqTM II kits (TaKaRa, Japan), according to the manufacturer's instructions. The primer sets used to detect LOC113809108 gene expression levels were designed using the Primer Premier software (version 5.0) [40] as follows: 5\u0026rsquo;-ACTTGGGTGCTGTAGCTCAA-3\u0026rsquo; and 5\u0026rsquo;-CTCGACAGCAACCAGGGTAT-3\u0026rsquo;. \u003cem\u003eL. vannamei\u003c/em\u003e 18S RNA was used as the internal reference gene; this gene was amplified using the primer sets as follows: 5\u0026rsquo;-GCCTGAGAAACGGCTACCACATC-3\u0026rsquo; and 5\u0026rsquo;-GTAGTAGCGACGGGCGGTGTGT-3\u0026rsquo; [41]. The qRT-PCR cycling program was as follows: preheating at 95\u0026deg;C for 30 s, followed by 40 cycles of 95\u0026deg;C for 5 s and 60\u0026deg;C for 30 s. The qRT-PCR was carried out at 95 \u0026deg;C for 40 s, 95 \u0026deg;C for 5 s, and 62 \u0026deg;C for 30 s for 40 cycles. Three parallel qRT-PCRs were carried out for each sample. Relative gene expression levels were calculated using the 2-\u0026Delta;\u0026Delta;CT method[42].\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003ePhenotypic variation \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe developed an ammonia-tolerant shrimp family (designated LV-N) for mapping, and subjected 284 LV-N shrimp to an acute ammonia stress test. All shrimp died within 2\u0026ndash;98 hours, with a mean survival time of 65 hours. Individual survival times were normally distributed and thus suitable for QTL detection. The accumulated mortality rate of shrimp is showed in Supplementary Material Fig S1.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSLAF-seq and genotyping\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBased on the digestive enzyme prediction using the reference genome of \u003cem\u003eL. vannamei\u003c/em\u003e, HaeIII and Hpy166II were used to digest the genomic DNA of the 284 LV-N shrimp for SLAF library construction. SLAF sequencing generated 439.77 gigabases (Gb) of data, consisting of 2,201 megabases (Mb) of 100-bp paired-end reads. Across all reads, the average Q30 was 95.81%, the average GC content was 40.60%, and the GC distribution was normal (Table 1). The rice (\u003cem\u003eOryza sativa japonica\u003c/em\u003e) genome was used as a control to estimate the validity of the library construction. For the rice library, 343.21 Mb of data (1.72 Mb paired-end reads) were generated, with a Q30 of 95.81% and a GC content of 40.96%. In \u003cem\u003eL. vannamei\u003c/em\u003e, 57.83% of the paired-end reads mapped successfully to the genome, as compared to 91.43% of the paired-end reads in rice. In addition, enzymatic digestion efficiency was 87.75% for \u003cem\u003eL. vannamei\u003c/em\u003e and 92.19% for rice (Supplementary Material, Table S2). These results indicated that SLAF library construction and sequencing were adequate.\u003c/p\u003e\n\u003cp\u003eAfter filtering and clustering all reads, 807,505 SLAFs were identified. The average sequencing depth of these SLAFs was 42.8-fold for the male parent, 42.14-fold for the female parent, and 12.43-fold for the progeny (Table 1). Of the 807,505 high-quality SLAFs detected, 293,415 (36.34%) were polymorphic (Table 1). After further filtering, the remaining 115,973 SLAF markers were successfully classified into eight genotypic patterns: ab \u0026times; cd, cc \u0026times; ab, aa \u0026times; bb, ab \u0026times; cc, ef \u0026times; eg, lm \u0026times; ll, hk \u0026times; hk, and nn \u0026times; np. The most common pattern was aa \u0026times; bb, followed by nn \u0026times; np and lm \u0026times; ll (Fig. 1). Because the mapped population was an F1 population, aa \u0026times; bb were eliminated as a valid marker.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCharacteristics of the genetic map\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLinkage analysis labeled 17,338 SLAF markers on the genetic map: 11,512 on the male map, 10,293 on the female map, and 17,338 on the sex-average map (Fig. 2). Each map contained 44 linkage groups (LGs). The total distances on the male, female, and sex-average maps were 6,604.99 cM, 5,476.20 cM, and 6,360.12 cM, respectively. The mean distance between adjacent markers was 0.58 cM on the male map, 0.53 cM on the female map, and 0.37 cM on the sex-average map (Supplementary Material, Table S3, Table S4, and Table S5). The distribution of markers among LGs was not uniform: in the male map, LG31 contained the most markers (585), while LG26 contained the least (39); in the female map, LG36 contained the most markers (540), while LG27 contained the least (21); and in the sex-average map, LG31 contained the most markers (695), while LG26 contained the least (53).\u003c/p\u003e\n\u003cp\u003e\u003cem\u003e\u0026nbsp;\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eQTL mapping of ammonia-tolerance\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eA QTL analysis of the ammonia-tolerance trait in the LV-N \u003cem\u003eL. vannamei\u003c/em\u003e family was performed based on the genetic maps. The LOD threshold was 4.75 (1000 permutations, P \u0026lt; 0.05). Thus, QTLs with LOD scores \u0026gt;4.75 were considered effective QTLs. Using this criterion, we identified a QTL within LG19 for ammonia tolerance (Fig. 3). The phenotypic variation explained by this QTL was 7.41\u0026ndash;8.46%, the LOD score was 4.75\u0026ndash;5.45, and the confidence interval was 12.42\u0026ndash;29.43 cM.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003e\u0026nbsp;\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTranscriptome sequencing, candidate gene identification and \u003c/strong\u003e\u003cstrong\u003eqRT-PCR verification\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe transcriptomes of the 20 most ammonia-tolerant and the 20 most ammonia-sensitive shrimp in each of four \u003cem\u003eL. vannamei \u003c/em\u003efamilies (LV-N, LV-A, LV-C, and LV-F) with various genetic backgrounds were sequenced. Transcriptome sequencing generated 56.79 Gb of clean data. A total of 7546 DEGs were identified between the ammonia-tolerant and ammonia-sensitive shrimp across all four families: 1,869 in LV-A, 2,005 in LV-C, 1,875 in LV-F, and 1,797 in LV-N (Supplementary Material, Table S6).\u003c/p\u003e\n\u003cp\u003eThe numbers of DEG annotations recovered in the databases searched were similar across the four \u003cem\u003eL. vannamei \u003c/em\u003efamilies. For instance, the COG terms mainly enriched in the DEGs from all four families were posttranslational modification, protein turnover, chaperones, and general function prediction only (Fig. 4); the GO terms primarily enriched in the DEGs from all four families were binding, catalytic activity, cellular process, metabolic process, cell, cell part, single-organism process and membrane functions (Fig. 5).\u003c/p\u003e\n\u003cp\u003eBy aligning the DEGs with the QTL region in LG19, 107 DEGs located in the QTL interval were identified. The expression levels and annotations of these DEGs are listed in Supplementary Material Table S7. Of these DEGs, only one gene (LOC113809108) met the criterion used to determine candidate genes associated with ammonia tolerance. This gene was annotated as an ATP synthase g subunit. LOC113809108 was located in the QTL interval, and was significantly upregulated in the most ammonia-tolerant shrimp compared to the most ammonia-sensitive shrimp from families LV-N and LV-C (Fig. 6). This gene was also upregulated in the most ammonia-tolerant shrimp from families LV-A and LV-F, but this difference in expression was not significant (Table 2).\u003c/p\u003e\n\u003cp\u003eThe qRT-PCR analysis showed that the patterns of LOC113809108 gene expression in ammonia-tolerant and ammonia-sensitive pooled samples from the families LV-A, LV-C, LV-F, and LV-N were similar to the patterns determined using RNA-seq: LOC113809108 gene expression was upregulated in the ammonia-tolerant shrimp as compared to the ammonia-sensitive shrimp across all four families (Fig. 6).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eThis study was aimed at investigating the ammonia tolerance in \u003cem\u003eL.vannamei\u003c/em\u003e by QTL analysis. A high-density genetic map of \u003cem\u003eL. vannamei\u003c/em\u003e was constructed using SLAF-seq, and a QTL associated with ammonia tolerance was identified as well as a putative candidate gene associated with ammonia tolerance.\u003c/p\u003e\n\u003cp\u003eThe genome of \u003cem\u003eL. vannamei\u003c/em\u003e is large (~2.45 Gb) [32]. Whole-genome deep resequencing is relatively costly for large genomes and is often not necessary for gene/QTL mapping [43, 44]. In recent years, several genetic linkage maps based on SNPs were constructed and QTL analyses were conducted in \u003cem\u003eL. vannamei\u003c/em\u003e [30, 45-47]. Yu et al. constructed a high-density genetic map for \u003cem\u003eL. vannamei\u003c/em\u003e and detected several QTLs for body weight and body length [30]. Yang et al. mapped the sex determination region in \u003cem\u003eL. vannamei\u003c/em\u003e based on the data used for high-density linkage map construction [45]. Du et al.mapped a QTL for \u003cem\u003eL. vannamei\u003c/em\u003e gender using a gene-based SNP linkage map [46]. In this study, a high-density genetic map of \u003cem\u003eL. vannamei\u003c/em\u003e was constructed using SLAF-seq, which is an effective method for discovering large numbers of SNPs and to perform large-scale genotyping [20]. Compared to traditional methods of genetic map construction (e.g., random amplified polymorphic DNA (RAPD), amplified fragment length polymorphism (AFLP), and simple sequence repeat (SSR)), the SLAF-seq method has several advantages for large-scale SNP discovery and genotyping: high density, high throughput, high efficiency, and low cost [48]. Previous studies have developed genetic maps of \u003cem\u003eL. vannamei\u003c/em\u003e using RAPD, AFLP, and SSR, but in these maps, the average distance between adjacent markers was 1\u0026ndash;5 cM [49-52]. The average distance between adjacent markers in the SLAF-seq genetic maps of \u003cem\u003eL. vannamei\u003c/em\u003e in this study was substantially shorter (0.34 cM). Notably, the average distance between adjacent markers found here was also less than in a previously reported SLAF-seq genetic map of \u003cem\u003eL. vannamei\u003c/em\u003e (0.75 cM) [30], possibly because a larger sample size and a greater sequencing depth were used. However, the number of LGs in the genetic map of \u003cem\u003eL. vannamei \u003c/em\u003ein this study (44) was consistent with the number of LGs in the previously reported genetic map of \u003cem\u003eL. vannamei\u003c/em\u003e [30]. This indicated that \u003cem\u003eL. vannamei \u003c/em\u003ehad 44 chromosomes, which was first reported by CamposRamos [53]. As far as we know, this is the highest density genetic linkage map \u003cem\u003eL. vannamei\u003c/em\u003e, and it is very useful for comparative analysis of genomic synteny, QTL mapping, positioning of candidate genes, and marker-assisted selection.\u003c/p\u003e\n\u003cp\u003eAs most animals cannot self-fertilize, it is difficult to develop common populations for genetic mapping (e.g., F2, recombinant inbred line (RIL), and nearlyisogenic line (NIL) populations). Therefore, an F1 population of \u003cem\u003eL. vannamei \u003c/em\u003ewas used to construct the genetic map, relying on a pseudo-testcross strategy. This strategy was based on the selection of single-dose markers present in one parent and absent in the other, and carried at a 1:1 ratio by the F1 offspring [54]. Therefore, gamete separation in each individual can be directly analyzed. The pseudo-testcross strategy has been widely used to construct animal F1 populations for genetic mapping [55-58]. In this study, an F1 population was developed using one ammonia-tolerant male parent (the result of 10 generations of selective breeding) and one female shrimp from a common family.\u003c/p\u003e\n\u003cp\u003ePrevious studies have suggested that the size of the mapped population might affect the accuracy of the genetic map and the QTL analysis, and have shown that genetic map accuracy increases with the size of the population used [59]. Specifically, populations of \u0026gt;200 individuals are considered sufficient for the construction of accurate genetic maps [59]. Thus, 284 randomly-selected individuals from the F1 population were used to construct the linkage map. However, the determination of shrimp survival time during the ammonia stress experiment depended on human observation, and thus may not have been perfectly accurate. To reduce the possible impacts of measurement inaccuracies on the QTL analysis, a relatively large F1 population was used. This larger population increased the accuracy of the QTL mapping, compensating for any instances of human error in survival time measurement.\u003c/p\u003e\n\u003cp\u003eAmmonia stress is one of the biggest challenges facing shrimp aquaculture. Ammonia not only has a direct lethal effect on shrimp [7], but also inhibits the shrimp immune system and increases sensitivity to pathogens [60]. Breeding new varieties of ammonia-tolerant shrimp is therefore an important target of the shrimp aquaculture industry. Some researches on shrimp ammonia tolerance have been carried out in recent years. Lu et al. found 12 SNPs associated with ammonia tolerance in \u003cem\u003eL. vannamei\u003c/em\u003e using marker-trait correlation analyses, and these SNPs were identified in the thrombospondin gene and X-box binding protein 1 gene [15]. At the same time, Lu et al. identified 202 differentially expressed proteins (DEPs) between ammonia-tolerant and ammonia-sensitive \u003cem\u003eL. vannamei\u003c/em\u003e using a comparative proteome analysis based on iTRAQ technique, and 77.8% of the DEPs were reported mainly involving in immune defense and stress tolerant in crustacean species [16]. Here, a QTL for ammonia-tolerance, located on LG19 at 169.09\u0026ndash;169.49 cM was identified, that explained 7.41\u0026ndash;8.46% of the phenotypic variation in ammonia tolerance. To the best of our knowledge, this is the first QTL for ammonia-tolerance reported in shrimp. However, having one QTL that explains about 8% of the variance indeed is not enough to initiate genetic breeding. Therefore in future research, the cutoff of LOD score can be relaxed to get more QTLs, intervals, and candidate DEGs. And then, more candidate genes related to ammonia tolerance should be identified. Moreover, independent families should be used to verify whether the QTL is common across different populations.\u003c/p\u003e\n\u003cp\u003eA QTL usually spans a large chromosomal region and may contain hundreds of genes. Therefore, in order to identify functional genes associated with ammonia tolerance, the transcriptomes of four \u003cem\u003eL. vannamei\u003c/em\u003e families (LV-N, LV-A, LV-C, and LV-F) with different genetic backgrounds were sequenced. Combining QTL mapping and gene expression analysis, we identified a single DEG (LOC113809108), located in the QTL interval, that was annotated as an ATP synthase g subunit and was significantly upregulated in the ammonia-tolerant LV-N and LV-C shrimp. The ATP synthase g subunit is located in F0 portion of ATP synthase, which consists of a membrane-integrated portion (F0 complex) and a membrane-protruding portion (F1 complex) [61]. The ATP synthase g subunit is essential for the formation of the F0 complex [61]. In the previous studies [15,16], ATP synthase has not been identified as being related to shrimp ammonia tolerance. ATP synthase is a double-motor enzyme that is involved in ATP synthesis, ATP hydrolysis-dependent processes, and the regulation of the proton gradient across some membrane-dependent systems [62]. Several studies have shown that ammonia excretion in aquatic animals is associated with Na\u003csup\u003e+\u003c/sup\u003e/K\u003csup\u003e+\u003c/sup\u003e-ATPase, which is mainly located on the basolateral membrane of branchial cells; NH\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e+\u003c/sup\u003e is excreted into the environment when K\u003csup\u003e+\u003c/sup\u003e is replaced by NH\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e+\u003c/sup\u003e via the Na\u003csup\u003e+\u003c/sup\u003e/NH\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e+\u003c/sup\u003e exchanger [63-65]. Indeed, a previous study suggested that in \u003cem\u003eL. vannamei\u003c/em\u003e, high ammonia tolerance was mainly a result of improved ammonia excretion and detoxification, as well as an accelerated energy metabolism [17]. Therefore, we speculate that ATP synthesis might affect the ammonia tolerance of \u003cem\u003eL. vannamei\u003c/em\u003e by regulating ATP synthesis and controlling cellular ammonia excretion. The results would provide useful information for further study of the molecular mechanisms of ammonia adaptive strategies in shrimps.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eIn this study, we constructed a high-density genetic map of L. vannamei and identified a QTL for ammonia tolerance. By combining QTL and transcriptome analyses, we identified a candidate gene associated with ammonia tolerance. The results help us better understand the molecular mechanism of ammonia tolerance in shrimp. Our work provides the basis for future genetic studies focused on molecular marker-assisted selective breeding.\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003ecM: centimorgans; SLAF-seq: specific length amplified fragment sequencing; QTL:quantitative trait locus; ammonia-N: ammonia nitrogen; SNPs: single nucleotide polymorphisms; PCRs: polymerase chain reactions; LOD: logarithm of odds; qRT-PCR: quantitative real-time PCR; FPKM: fragments per kilobase of transcript per million mapped reads; FDR: false discovery rate; DEGs: differentially expressed genes; Gb: gigabases; Mb: megabases; LGs: linkage groups\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe thank the Guangxi Shrimp Breeding Engineering Technology Research Center for providing animal materials.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026rsquo; contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eYZ and QL conceived the study. DZ and CY wrote the manuscript and contributed in the bioinformatics analysis. QL, WZ, XC, and MP conducted the experiments and contributed for raw data analysis. XC, YL, HW, HL, and JL contributed in the bioinformatics analysis. All authors read and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by the scientific and technological innovation major project of guangxi (aa17204080-1), scientific and technological innovation major base of guangxi (aa17204088-1), guangxi scientific research and technology development project (ab16380189), the national modern agriculture industry technology system project (nycytxgxcxtd-14-01), and the national modern agroindustry technology research system (cars-48).\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll data generated or analysed during this study are included in this published article and its supplementary information files. Raw SLAF sequencing reads are deposited in NCBI database under the accession numbers PRJNA545592. Raw RNA-seq are deposited in NCBI database under the accession numbers SRR9822091, SRR9822090, SRR9822095, SRR9822094, SRR9822098, SRR9822099, SRR9822093, and SRR9822085.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare they have no competing interests.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eCruz P, Ibarra AM, Mejia-Ruiz H, Gaffney PM, Perez-Enriquez R: Genetic variability assessed by microsatellites in a breeding program of Pacific white shrimp (\u003cem\u003eLitopenaeus vannamei\u003c/em\u003e). Mar Biotechnol (NY). 2004; 6(2):157-164.\u003c/li\u003e\n\u003cli\u003eJavahery S, Noori A, Hoseinifar SH: Growth performance, immune response, and digestive enzyme activity in Pacific white shrimp, \u003cem\u003ePenaeus vannamei\u003c/em\u003e Boone, 1931, fed dietary microbial lysozyme. Fish \u0026amp; shellfish immunology. 2019; 92:528-535.\u003c/li\u003e\n\u003cli\u003eAlagappan KM, Deivasigamani B, Somasundaram ST, Kumaran S: Occurrence of Vibrio parahaemolyticus and its specific phages from shrimp ponds in east coast of India. Current microbiology. 2010; 61(4):235-240.\u003c/li\u003e\n\u003cli\u003eSun CB, Wang G, Chan SF: Effects of artificial infection of \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e by Micrococcus lysodeikticus and WSSV on the activity of immunity related enzymes. Fish \u0026amp; shellfish immunology. 2015; 46(2):778-786.\u003c/li\u003e\n\u003cli\u003eMarazza D, Bornens P, Le Gal Y: Effect of ammonia on survival and adenylate energy charge in the shrimp \u003cem\u003ePalaemonetes varians\u003c/em\u003e. Ecotoxicology and environmental safety. 1996; 34(2):103-108.\u003c/li\u003e\n\u003cli\u003eFrias-Espericueta MG, Harfush-Melendez M, Paez-Osuna F: Effects of ammonia on mortality and feeding of postlarvae shrimp \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Bulletin of environmental contamination and toxicology. 2000; 65(1):98-103.\u003c/li\u003e\n\u003cli\u003eQiu L, Shi X, Yu S, Han Q, Diao X, Zhou H: Changes of Ammonia-Metabolizing Enzyme Activity and Gene Expression of Two Strains in Shrimp \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e Under Ammonia Stress. Frontiers in physiology. 2018; 9:211.\u003c/li\u003e\n\u003cli\u003eLiang Z, Liu R, Zhao D, Wang L, Sun M, Wang M, Song L: Ammonia exposure induces oxidative stress, endoplasmic reticulum stress and apoptosis in hepatopancreas of pacific white shrimp (\u003cem\u003eLitopenaeus vannamei\u003c/em\u003e). Fish \u0026amp; shellfish immunology. 2016; 54:523-528.\u003c/li\u003e\n\u003cli\u003eLu X, Kong J, Luan S, Dai P, Meng X, Cao B, Luo K: Transcriptome Analysis of the Hepatopancreas in the Pacific White Shrimp (\u003cem\u003eLitopenaeus vannamei\u003c/em\u003e) under Acute Ammonia Stress. PloS one. 2016; 11(10):e0164396.\u003c/li\u003e\n\u003cli\u003eLu X, Luan S, Cao B, Meng X, Sui J, Dai P, Luo K, Shi X, Hao D, Han G\u003cem\u003e et al\u003c/em\u003e: Estimation of genetic parameters and genotype-by-environment interactions related to acute ammonia stress in Pacific white shrimp (\u003cem\u003eLitopenaeus vannamei\u003c/em\u003e) juveniles at two different salinity levels. PloS one. 2017; 12(3):e0173835.\u003c/li\u003e\n\u003cli\u003eDuan Y, Zhang Y, Dong H, Wang Y, Zheng X, Zhang J: Effect of dietary Clostridium butyricum on growth, intestine health status and resistance to ammonia stress in Pacific white shrimp \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Fish \u0026amp; shellfish immunology. 2017; 65:25-33.\u003c/li\u003e\n\u003cli\u003eValencia-Castaneda G, Frias-Espericueta MG, Vanegas-Perez RC, Perez-Ramirez JA, Chavez-Sanchez MC, Paez-Osuna F: Acute Toxicity of Ammonia, Nitrite and Nitrate to Shrimp \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e Postlarvae in Low-Salinity Water. Bulletin of environmental contamination and toxicology. 2018; 101(2):229-234.\u003c/li\u003e\n\u003cli\u003eCampbell NR, LaPatra SE, Overturf K, Towner R, Narum SR: Association mapping of disease resistance traits in rainbow trout using restriction site associated DNA sequencing. G3 (Bethesda). 2014; 4(12):2473-2481.\u003c/li\u003e\n\u003cli\u003eYu Y, Liu J, Li F, Zhang X, Zhang C, Xiang J: Gene set based association analyses for the WSSV resistance of Pacific white shrimp \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Scientific reports. 2017; 7:40549.\u003c/li\u003e\n\u003cli\u003eLu X, Kong J, Meng X, Cao B, Luo K, Dai P, Luan S: Identification of SNP markers associated with tolerance to ammonia toxicity by selective genotyping from de novo assembled transcriptome in \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Fish \u0026amp; shellfish immunology. 2018; 73:158-166.\u003c/li\u003e\n\u003cli\u003eLu X, Luan S, Dai P, Meng X, Cao B, Luo K, Kong J: iTRAQ-based comparative proteome analysis for molecular mechanism of defense against acute ammonia toxicity in Pacific White shrimp \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Fish \u0026amp; shellfish immunology. 2018; 74:52-61.\u003c/li\u003e\n\u003cli\u003eXiao J, Li QY, Tu JP, Chen XL, Chen XH, Liu QY, Liu H, Zhou XY, Zhao YZ, Wang HL: Stress response and tolerance mechanisms of ammonia exposure based on transcriptomics and metabolomics in \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Ecotoxicology and environmental safety. 2019; 180:491-500.\u003c/li\u003e\n\u003cli\u003eLu X, Luan S, Dai P, Luo K, Chen B, Cao B, Sun L, Yan Y, Kong J: Insights into the molecular basis of immunosuppression and increasing pathogen infection severity of ammonia toxicity by transcriptome analysis in pacific white shrimp \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Fish \u0026amp; shellfish immunology. 2019; 88:528-539.\u003c/li\u003e\n\u003cli\u003eZhu M, Zhao S: Candidate gene identification approach: progress and challenges. International journal of biological sciences. 2007; 3(7):420-427.\u003c/li\u003e\n\u003cli\u003eSun X, Liu D, Zhang X, Li W, Liu H, Hong W, Jiang C, Guan N, Ma C, Zeng H\u003cem\u003e et al\u003c/em\u003e: SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PloS one. 2013; 8(3):e58700.\u003c/li\u003e\n\u003cli\u003eQian W, Fan G, Liu D, Zhang H, Wang X, Wu J, Xu Z: Construction of a high-density genetic map and the X/Y sex-determining gene mapping in spinach based on large-scale markers developed by specific-locus amplified fragment sequencing (SLAF-seq). BMC genomics. 2017; 18(1):276.\u003c/li\u003e\n\u003cli\u003eZhang Y, Wang L, Xin H, Li D, Ma C, Ding X, Hong W, Zhang X: Construction of a high-density genetic map for sesame based on large scale marker development by specific length amplified fragment (SLAF) sequencing. BMC plant biology. 2013; 13:141.\u003c/li\u003e\n\u003cli\u003eZhu Y, Yin Y, Yang K, Li J, Sang Y, Huang L, Fan S: Construction of a high-density genetic map using specific length amplified fragment markers and identification of a quantitative trait locus for anthracnose resistance in walnut (\u003cem\u003eJuglans regia\u003c/em\u003e L.). BMC genomics. 2015; 16:614.\u003c/li\u003e\n\u003cli\u003eLi B, Tian L, Zhang J, Huang L, Han F, Yan S, Wang L, Zheng H, Sun J: Construction of a high-density genetic map based on large-scale markers developed by specific length amplified fragment sequencing (SLAF-seq) and its application to QTL analysis for isoflavone content in \u003cem\u003eGlycine max\u003c/em\u003e. BMC genomics. 2014; 15:1086.\u003c/li\u003e\n\u003cli\u003eWei Q, Wang Y, Qin X, Zhang Y, Zhang Z, Wang J, Li J, Lou Q, Chen J: An SNP-based saturated genetic map and QTL analysis of fruit-related traits in cucumber using specific-length amplified fragment (SLAF) sequencing. BMC genomics. 2014; 15:1158.\u003c/li\u003e\n\u003cli\u003eJiang B, Liu W, Xie D, Peng Q, He X, Lin Y, Liang Z: High-density genetic map construction and gene mapping of pericarp color in wax gourd using specific-locus amplified fragment (SLAF) sequencing. BMC genomics. 2015; 16:1035.\u003c/li\u003e\n\u003cli\u003eZhao Z, Gu H, Sheng X, Yu H, Wang J, Huang L, Wang D: Genome-Wide Single-Nucleotide Polymorphisms Discovery and High-Density Genetic Map Construction in Cauliflower Using Specific-Locus Amplified Fragment Sequencing. Frontiers in plant science. 2016; 7:334.\u003c/li\u003e\n\u003cli\u003eTao A, Huang L, Wu G, Afshar RK, Qi J, Xu J, Fang P, Lin L, Zhang L, Lin P: High-density genetic map construction and QTLs identification for plant height in white jute (\u003cem\u003eCorchorus capsularis \u003c/em\u003eL.) using specific locus amplified fragment (SLAF) sequencing. BMC genomics. 2017; 18(1):355.\u003c/li\u003e\n\u003cli\u003eLiu C, Zhou Q, Dong L, Wang H, Liu F, Weng J, Li X, Xie C: Genetic architecture of the maize kernel row number revealed by combining QTL mapping using a high-density genetic map and bulked segregant RNA sequencing. BMC genomics. 2016; 17(1):915.\u003c/li\u003e\n\u003cli\u003eYu Y, Zhang X, Yuan J, Li F, Chen X, Zhao Y, Huang L, Zheng H, Xiang J: Genome survey and high-density genetic map construction provide genomic and genetic resources for the Pacific White Shrimp \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Scientific reports. 2015; 5:15612.\u003c/li\u003e\n\u003cli\u003eGilcreas FW: Future of standard methods for the examination of water and wastewater. Health laboratory science. 1967; 4(3):137-141.\u003c/li\u003e\n\u003cli\u003eZhang X, Yuan J, Sun Y, Li S, Gao Y, Yu Y, Liu C, Wang Q, Lv X, Ma KY\u003cem\u003e et al\u003c/em\u003e: Penaeid shrimp genome provides insights into benthic adaptation and frequent molting. Nature communications. 2019; 10(1):356.\u003c/li\u003e\n\u003cli\u003eLi H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25(14):1754-1760.\u003c/li\u003e\n\u003cli\u003eLiu D, Ma C, Hong W, Huang L, Liu M, Liu H, Zeng H, Deng D, Xin H, Song J\u003cem\u003e et al\u003c/em\u003e: Construction and analysis of high-density linkage map using high-throughput sequencing data. PloS one. 2014; 9(6):e98855.\u003c/li\u003e\n\u003cli\u003eBroman KW, Gatti DM, Simecek P, Furlotte NA, Prins P, Sen S, Yandell BS, Churchill GA: R/qtl2: Software for Mapping Quantitative Trait Loci with High-Dimensional Data and Multiparent Populations. Genetics. 2019; 211(2):495-502.\u003c/li\u003e\n\u003cli\u003eGardiner SE, Norelli JL, de Silva N, Fazio G, Peil A, Malnoy M, Horner M, Bowatte D, Carlisle C, Wiedow C\u003cem\u003e et al\u003c/em\u003e: Putative resistance gene markers associated with quantitative trait loci for fire blight resistance in Malus 'Robusta 5' accessions. BMC genetics. 2012; 13:25.\u003c/li\u003e\n\u003cli\u003eKim D, Langmead B, Salzberg SL: HISAT: a fast spliced aligner with low memory requirements. Nature methods. 2015; 12(4):357-360.\u003c/li\u003e\n\u003cli\u003eMortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature methods. 2008; 5(7):621-628.\u003c/li\u003e\n\u003cli\u003eRobinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010; 26(1):139-140.\u003c/li\u003e\n\u003cli\u003eSingh VK, Mangalam AK, Dwivedi S, Naik S: Primer premier: program for design of degenerate primers from a protein sequence. BioTechniques. 1998; 24(2):318-319.\u003c/li\u003e\n\u003cli\u003eZhang Q, Li F, Zhang X, Dong B, Zhang J, Xie Y, Xiang J: cDNA cloning, characterization and expression analysis of the antioxidant enzyme gene, catalase, of Chinese shrimp \u003cem\u003eFenneropenaeus chinensis\u003c/em\u003e. Fish \u0026amp; shellfish immunology. 2008; 24(5):584-591.\u003c/li\u003e\n\u003cli\u003eLivak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001; 25(4):402-408.\u003c/li\u003e\n\u003cli\u003eXu X, Chao J, Cheng X, Wang R, Sun B, Wang H, Luo S, Wu T, Li Y: Mapping of a Novel Race Specific Resistance Gene to Phytophthora Root Rot of Pepper (\u003cem\u003eCapsicum annuum\u003c/em\u003e) Using Bulked Segregant Analysis Combined with Specific Length Amplified Fragment Sequencing Strategy. PloS one. 2016; 11(3):e0151401.\u003c/li\u003e\n\u003cli\u003eWang G, Chen B, Du H, Zhang F, Zhang H, Wang Y, He H, Geng S, Zhang X: Genetic mapping of anthocyanin accumulation-related genes in pepper fruits using a combination of SLAF-seq and BSA. PloS one. 2018; 13(9):e0204690.\u003c/li\u003e\n\u003cli\u003eYu Y, Zhang X, Yuan J, Wang Q, Li S, Huang H, Li F, Xiang J: Identification of Sex-determining Loci in Pacific White Shrimp \u003cem\u003eLitopeneaus vannamei\u003c/em\u003e Using Linkage and Association Analysis. Marine Biotechnology. 2017; 19: 277\u0026ndash;286.\u003c/li\u003e\n\u003cli\u003eDu ZQ, Ciobanu DC, Onteru SK, Gorbach D, Mileham AJ, Jaramillo G, Rothschild MF: A gene-based SNP linkage map for pacific white shrimp, \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Anim Genet 2010, 41(3); 286-294.\u003c/li\u003e\n\u003cli\u003eDu Z, Onteru S, Gorbach D, Rothschild M: A SNP Genetic Map for Pacific White Shrimp (\u003cem\u003eLitopenaeus vannamei\u003c/em\u003e). Animal Industry Report. 2010; 656(1):35.\u003c/li\u003e\n\u003cli\u003eLiang D, Chen M, Qi X, Xu Q, Zhou F, Chen X: QTL Mapping by SLAF-seq and Expression Analysis of Candidate Genes for Aphid Resistance in Cucumber. Frontiers in plant science. 2016; 7:1000.\u003c/li\u003e\n\u003cli\u003eGarcia DK, Dhar AK, Alcivar-Warren A: Molecular analysis of a RAPD marker (B20) reveals two microsatellites and differential mRNA expression in Penaeus vannamei. Molecular marine biology and biotechnology. 1996; 5(1):71-83.\u003c/li\u003e\n\u003cli\u003ePerez F, Ortiz J, Zhinaula M, Gonzabay C, Calderon J, Volckaert FA: Development of EST-SSR markers by data mining in three species of shrimp: \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e, \u003cem\u003eLitopenaeus stylirostris\u003c/em\u003e, and \u003cem\u003eTrachypenaeus birdy\u003c/em\u003e. Mar Biotechnol (NY). 2005; 7(5):554-569.\u003c/li\u003e\n\u003cli\u003eAndriantahina F, Liu X, Huang H: Genetic map construction and quantitative trait locus (QTL) detection of growth-related traits in Litopenaeus vannamei for selective breeding applications. PloS one. 2013; 8(9):e75206.\u003c/li\u003e\n\u003cli\u003eZhang L, Yang C, Zhang Y, Li L, Zhang X, Zhang Q, Xiang J: A genetic linkage map of Pacific white shrimp (\u003cem\u003eLitopenaeus vannamei\u003c/em\u003e): sex-linked microsatellite markers and high recombination rates. Genetica. 2007; 131(1):37-49.\u003c/li\u003e\n\u003cli\u003eCamposRamos, R: Chromosome studies on the marine shrimps \u003cem\u003ePenaeus vannamei\u003c/em\u003e and P-californiensis (Decapoda). Journal of Crustacean Biology. 1997; 17(4), 666\u0026ndash;673.\u003c/li\u003e\n\u003cli\u003eGrattapaglia D, Sederoff R: Genetic linkage maps of \u003cem\u003eEucalyptus grandis\u003c/em\u003e and \u003cem\u003eEucalyptus urophylla\u003c/em\u003e using a pseudo-testcross: mapping strategy and RAPD markers. Genetics. 1994; 137(4):1121-1137.\u003c/li\u003e\n\u003cli\u003eGrattapaglia D, Bertolucci FL, Sederoff RR: Genetic mapping of QTLs controlling vegetative propagation in \u003cem\u003eEucalyptus grandis\u003c/em\u003e and \u003cem\u003eE. urophylla\u003c/em\u003e using a pseudo-testcross strategy and RAPD markers. TAG Theoretical and applied genetics Theoretische und angewandte Genetik. 1995; 90(7-8):933-947.\u003c/li\u003e\n\u003cli\u003ePorceddu A, Albertini E, Barcaccia G, Falistocco E, Falcinelli M: Linkage mapping in apomictic and sexual \u003cem\u003eKentucky bluegrass\u003c/em\u003e (\u003cem\u003ePoa pratensis\u003c/em\u003e L.) genotypes using a two way pseudo-testcross strategy based on AFLP and SAMPL markers. TAG Theoretical and applied genetics Theoretische und angewandte Genetik. 2002; 104(2-3):273-280.\u003c/li\u003e\n\u003cli\u003eLiao M, Zhang L, Yang G, Zhu M, Wang D, Wei Q, Zou G, Chen D: Development of silver carp (\u003cem\u003eHypophthalmichthys molitrix\u003c/em\u003e) and bighead carp (\u003cem\u003eAristichthys nobilis\u003c/em\u003e) genetic maps using microsatellite and AFLP markers and a pseudo-testcross strategy. Animal genetics. 2007; 38(4):364-370.\u003c/li\u003e\n\u003cli\u003eXie W, Zhang X, Cai H, Huang L, Peng Y, Ma X: Genetic maps of SSR and SRAP markers in diploid orchardgrass (\u003cem\u003eDactylis glomerata\u003c/em\u003e L.) using the pseudo-testcross strategy. Genome. 2011; 54(3):212-221.\u003c/li\u003e\n\u003cli\u003eAd\u0026eacute;sio FerreiraI MFdS, Luciano da Costa e SilvaI, Cosme Dami\u0026atilde;o CruzI: Estimating the effects of population size and type on the accuracy of genetic maps. Genetics and Molecular Biology. 2006; 29(1):187-192.\u003c/li\u003e\n\u003cli\u003eJiang L, Feng J, Ying R, Yin F, Pei S, Lu J, Cao Y, Guo J, Li Z: Individual and combined effects of ammonia-N and sulfide on the immune function and intestinal microbiota of Pacific white shrimp \u003cem\u003eLitopenaeus vannamei\u003c/em\u003e. Fish \u0026amp; shellfish immunology. 2019; 92:230-240.\u003c/li\u003e\n\u003cli\u003eFujikawa M, Sugawara K, Tanabe T, Yoshida M: Assembly of human mitochondrial ATP synthase through two separate intermediates, F1-c-ring and b-e-g complex. FEBS Letter. 2015; 589(19 Pt B):2707-2712.\u003c/li\u003e\n\u003cli\u003eHong S, Pedersen PL: ATP synthase and the actions of inhibitors utilized to study its roles in human health, disease, and other scientific areas. Microbiology and molecular biology reviews : MMBR. 2008; 72(4):590-641, Table of Contents.\u003c/li\u003e\n\u003cli\u003eWilkie MP: Mechanisms of ammonia excretion across fish gills. Comparative Biochemistry and Physiology Part A. 1997; 118A(1):39-50.\u003c/li\u003e\n\u003cli\u003eLucu C, Flik G: Na+-K+-ATPase and Na+/Ca2+ exchange activities in gills of hyperregulating \u003cem\u003eCarcinus maenas\u003c/em\u003e. The American journal of physiology. 1999; 276(2):R490-499.\u003c/li\u003e\n\u003cli\u003eWeihrauch D, Morris S, Towle DW: Ammonia excretion in aquatic and terrestrial crabs. The Journal of experimental biology. 2004; 207(Pt 26):4491-4504.\u003c/li\u003e\n\u003c/ol\u003e\n"},{"header":"Tables","content":"\u003cp\u003e\u003cstrong\u003eTable 1. Summary of the constructed genetic map of \u003cem\u003eLitopenaeus vannamei.\u003c/em\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003ctable border=\"1\" width=\"439\"\u003e\n\u003ctbody\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eMap data\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003eValue\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eTotal bases\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e439.77 Gb\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eTotal reads\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e2,201.28 Mb\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eAverage Q30\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e95.81%\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eAverage GC\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e40.60%\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eEnzyme digestion protocol\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003eHaeIII+Hpy166II\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eRestriction fragment length\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e314-414 bp\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003ePercentage of reads matching the \u003cem\u003eL. vannamei\u003c/em\u003e genome\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e57.83%\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eAverage enzymatic digestion efficiency\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e87.75%\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003ePredicted number of markers\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e339,517\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eNumber of high-quality slafs\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e807,505\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eNumber of polymorphic slafs\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e293,415\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eNumber of SLAF markers on the map\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e17,338\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eAverage depth in parents\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e208.90 \u0026times;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eAverage depth in offspring individual\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e38.55 \u0026times;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eNumber of linkage groups\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e44\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eTotal distance of the map\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e6,360.12 cM\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"311\"\u003e\n\u003cp\u003eAverage distance of the map\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"128\"\u003e\n\u003cp\u003e0.37 cM\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e\u003cbr /\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTable 2. \u003c/strong\u003e\u003cstrong\u003eRNA-Seq analysis showing expression of gene LOC113809108 in the most ammonia-tolerant and the most ammonia-sensitive \u003c/strong\u003e\u003cstrong\u003eindividuals \u003c/strong\u003e\u003cstrong\u003eacross four families \u003c/strong\u003e\u003cstrong\u003eof \u003c/strong\u003e\u003cstrong\u003e\u003cem\u003eLitopenaeus vannamei\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e (LV-A, LV-C, LV-F, and LV-N).\u003c/strong\u003e\u003c/p\u003e\n\u003ctable border=\"1\" width=\"521\"\u003e\n\u003ctbody\u003e\n\u003ctr\u003e\n\u003ctd rowspan=\"2\" width=\"49\"\u003e\n\u003cp\u003eFamily\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd colspan=\"8\" width=\"472\"\u003e\n\u003cp\u003eGene LOC113809108 expression level\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003eSensitive group (FPKM)\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003eSensitive group (count)\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003eTolerant group (FPKM)\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003eTolerant group (count)\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003eFDR\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003eLog2-fold Change\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003eRegulated\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"72\"\u003e\n\u003cp\u003eDifference\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"49\"\u003e\n\u003cp\u003eLV-A\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e112.8410\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003e3239\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e130.2663\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003e2561\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003e0.9589\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e0.2072\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003eup\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"72\"\u003e\n\u003cp\u003eNormal\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"49\"\u003e\n\u003cp\u003eLV-C\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e63.2980\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003e1651\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e148.3678\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003e2444\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003e0.0002\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e1.2289\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003eup\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"72\"\u003e\n\u003cp\u003eSignificant\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"49\"\u003e\n\u003cp\u003eLV-F\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e70.4520\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003e1378\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e122.4527\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003e2258\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003e0.9402\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e0.7975\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003eup\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"72\"\u003e\n\u003cp\u003eNormal\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"49\"\u003e\n\u003cp\u003eLV-N\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e114.9725\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003e2659\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e329.0251\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"56\"\u003e\n\u003cp\u003e8025\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003e0.0001\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"64\"\u003e\n\u003cp\u003e1.5169\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"48\"\u003e\n\u003cp\u003eup\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"72\"\u003e\n\u003cp\u003eSignificant\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"genetic map, QTL, transcriptomic, ammonia tolerance, Litopenaeus vannamei","lastPublishedDoi":"10.21203/rs.2.18942/v3","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.2.18942/v3","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground:\u003c/strong\u003e Ammonia is one of the most common toxicological environment factors affecting shrimp health. Although ammonia tolerance in shrimp\u003cem\u003e \u003c/em\u003eis closely related to successful industrial production, few genetic studies of this trait are available. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults:\u003c/strong\u003e In this study, we constructed a high-density genetic map of the Pacific white shrimp (\u003cem\u003eLitopenaeus vannamei\u003c/em\u003e) using specific length amplified fragment sequencing (SLAF-seq). The constructed genetic map contained 17,338 polymorphic markers spanning 44 linkage groups, with a total distance of 6,360.12 centimorgans (cM) and an average distance of 0.37 cM. Using this genetic map, we identified a quantitative trait locus (QTL) that explained 7.41–8.46% of the phenotypic variance in \u003cem\u003eL. vannamei\u003c/em\u003e survival time under acute ammonia stress. We then sequenced the transcriptomes of the most ammonia-tolerant and the most ammonia-sensitive individuals from each of four genetically distinct \u003cem\u003eL. vannamei\u003c/em\u003e families. We found that 7546 genes were differentially expressed between the ammonia-tolerant and ammonia-sensitive individuals. Using QTL analysis and the transcriptomes, we identified one candidate gene (annotated as an ATP synthase g subunit) associated with ammonia tolerance. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusions:\u003c/strong\u003e In this study, we constructed a high-density genetic map of \u003cem\u003eL. vannamei \u003c/em\u003eand identified a QTL for ammonia tolerance. By combining QTL and transcriptome analyses, we identified a candidate gene associated with ammonia tolerance. Our work provides the basis for future genetic studies focused on molecular marker-assisted selective breeding.\u003c/p\u003e","manuscriptTitle":"Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei)","msid":"","msnumber":"","nonDraftVersions":[{"code":3,"date":"2020-09-29 21:43:08","doi":"10.21203/rs.2.18942/v3","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Accept","date":"2020-11-18T00:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2020-09-20T12:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2020-09-19T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2020-09-19T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":2,"date":"2020-08-11 18:10:45","doi":"10.21203/rs.2.18942/v2","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Minor revision","date":"2020-09-16T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2020-09-11T12:00:00+00:00","index":1,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"reviewerAgreed","content":"","date":"2020-08-26T12:00:00+00:00","index":1,"fulltext":""},{"type":"editorAssigned","content":"","date":"2020-08-03T12:00:00+00:00","index":"","fulltext":""},{"type":"reviewersInvited","content":"","date":"2020-08-03T12:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2020-08-02T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2020-08-02T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":1,"date":"2019-12-16 20:30:30","doi":"10.21203/rs.2.18942/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revision","date":"2020-07-06T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2020-07-01T12:00:00+00:00","index":2,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"editorInvitedReview","content":"","date":"2020-02-04T12:00:00+00:00","index":1,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"reviewerAgreed","content":"","date":"2020-01-16T12:00:00+00:00","index":2,"fulltext":""},{"type":"reviewerAgreed","content":"","date":"2019-12-20T12:00:00+00:00","index":1,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2019-12-19T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2019-12-11T12:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2019-12-05T12:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2019-12-04T12:00:00+00:00","index":"","fulltext":""},{"type":"submitted","content":"","date":"2019-12-03T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"d69a608c-7d08-40ba-b5a3-9ffa84f41684","owner":[],"postedDate":"September 29th, 2020","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":278723,"name":"Epigenetics \u0026 Genomics"},{"id":278724,"name":"Aquaculture and Mariculture"}],"tags":[],"updatedAt":"2020-12-06T15:02:31+00:00","versionOfRecord":{"articleIdentity":"rs-9560","link":"https://doi.org/10.1186/s12864-020-07254-x","journal":{"identity":"bmc-genomics","isVorOnly":false,"title":"BMC Genomics"},"publishedOn":"2020-12-02 15:00:46","publishedOnDateReadable":"December 2nd, 2020"},"versionCreatedAt":"2020-09-29 21:43:08","video":"","vorDoi":"10.1186/s12864-020-07254-x","vorDoiUrl":"https://doi.org/10.1186/s12864-020-07254-x","workflowStages":[]},"version":"v3","identity":"rs-9560","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"identity":"rs-9560","version":["v3"]},"buildId":"_2-kVJe1T_tPrBINL-cwx","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}