Identification of SNP markers associated with tolerance to ammonia toxicity by selective genotyping from de novo assembled transcriptome in Litopenaeus vannamei

Transcriptomic Insights and the Development of Microsatellite Markers to Assess Genetic Diversity in the Broodstock Management of Litopenaeus stylirostris

The Pacific blue shrimp (Litopenaeus stylirostris) is a premium product in the international seafood market. However, intensified farming has increased disease incidence and reduced genetic diversity. In this study, we developed a transcriptome database for L. stylirostris and mined microsatellite markers to analyze their genetic diversity. Using the Illumina HiSeq 4000 platform, we identified 53,263 unigenes from muscle, hepatopancreas, the intestine, and lymphoid tissues. Microsatellite analysis identified 36,415 markers from 18,657 unigenes, predominantly dinucleotide repeats. Functional annotation highlighted key disease resistance pathways and enriched categories. The screening and PCR testing of 42 transcriptome-based and 58 literature-based markers identified 40 with successful amplification. The genotyping of 200 broodstock samples revealed that Na, Ho, He, PIC, and FIS values were 3, 0.54 ± 0.05, 0.43 ± 0.09, 0.41 ± 0.22, and 0.17 ± 0.27, respectively, indicating moderate genetic variability and significant inbreeding. Four universal microsatellite markers (CL1472.Contig13, CL517.Contig2, Unigene5692, and Unigene7147) were identified for precise diversity analysis in Pacific blue, Pacific white (Litopenaeus vannamei), and black tiger shrimps (Penaeus monodon). The transcriptome database supports the development of markers and functional gene analysis for selective breeding programs. Our findings underscore the need for an appropriate genetic management system to mitigate inbreeding depression, reduce disease susceptibility, and preserve genetic diversity in farmed shrimp populations.

A review of the toxic effects of ammonia on invertebrates in aquatic environments

Aquatic invertebrates are the organisms most susceptible to ammonia toxicity. However, the toxic effects of ammonia on invertebrates are still poorly understood. This study reviews the research progress in ammonia toxicology for the period from 1986 to 2023, focusing on the effects on invertebrates. Through examining the toxic effects of ammonia at different levels of organization (community, individual, tissue and physiology, and molecular) as well as the results from omics studies, we determined that the most significant effects were on the reproductive capacity of invertebrates and the growth of offspring, although different populations show variation in their tolerance to ammonia, and tissues have varied potential to respond to ammonia stress. A multicomponent analysis is an in-depth technique employed in toxicological studies, as it can be used to explore the enrichment pathways and functional genes expressed under ammonia stress. This study comprehensively discusses ammonia toxicity from multiple aspects in order to provide new insights into the toxic effects of ammonia on aquatic invertebrates.

First transcriptome profiling in gill and hepatopancrease tissues of Metapenaeus ensis in response to acute ammonia-N stress

The greasyback shrimp, Metapenaeus ensis, suffers from ammonia-N stress during intensive factory aquaculture. Optimizing ammonia-N stress tolerance has become an important issue in M. ensis breeding. The metabolic and adaptive mechanisms of ammonia-N toxicity in M. ensis have not been comprehensively understood yet. In this study, a large number of potential simple sequence repeats (SSRs) in the transcriptome of M. ensis were identified. Differentially expressed genes (DEGs) in the gill and hepatopancreas at 24 h post-challenges under high concentrations of ammonia-N treatment were detected. We obtained 20,108,851-27,681,918 clean reads from the control and high groups, assembled and clustered a total of 103,174 unigenes with an average of 876 bp and an N50 of 1189 bp. Comparative transcriptome analyses identified 2000 different expressed genes in the gill and 2010 different expressed genes in the hepatopancreas, a large number of which were related to immune function, oxidative stress, metabolic regulation, and apoptosis. The results suggest that M. ensis may counteract ammonia-N toxicity at the transcriptome level by increasing the expression of genes related to immune stress and detoxification metabolism, and that selected genes may serve as molecular indicators of ammonia-N. By exploring the genetic basis of M. ensis' ammonia-N stress adaptation, we constructed the genetic networks for ammonia-N adaptation. These findings will accelerate the understanding of M. ensis' ammonia-N adaptation, contribute to the research of future breeding, and promote the level of factory aquaculture of M. ensis.

New insights into the regulation mechanism of red claw crayfish (Cherax quadricarinatus) hepatopancreas under air exposure using transcriptome analysis

Red claw crayfish (Cherax quadricarinatus) is an important freshwater shrimp species worldwide with enormous economic value. Waterless transportation is an inherent feature of red claw crayfish transportation. However, the high mortality of red claw crayfish is a severe problem in the aquaculture of crayfish after waterless transportation. In this study, we investigated the responses of the hepatopancreas from the red claw crayfish undergoing air exposure stress and normal conditions on transcriptome levels. We used Illumina-based RNA sequencing (RNA-Seq) to perform a transcriptome analysis from the hepatopancreas of red claw crayfish challenged by air exposure. An average of 57,148,800 clean reads per library was obtained, and 33,567 unigenes could be predicted and classified according to their homology with matches in the National Center for Biotechnology Information (NCBI) non-redundant protein sequences (Nr), Gene Ontology (GO), a manually annotated and reviewed protein sequence database (Swiss-Prot), protein families (Pfam), Clusters of Orthologous Groups (COG) of proteins, and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. 690 and 3407 differentially expressed genes (DEGs) were identified between the two stress stages of the red claw crayfish. More DEGs were identified in 12 h, indicating that gene expressions were largely changed at 12 h. Some immune-related pathways and genes were identified according to KEGG and GO enrichment analysis. A total of 12 DEGs involved in immune response and trehalose mechanism were verified by quantitative real-time-polymerase chain reaction (qRT-PCR). The results indicated that the red claw crayfish might counteract the stress of air exposure at the transcriptomic level by increasing expression levels of antioxidant-, immune-, and trehalose metabolism-related genes. These transcriptome results from the hepatopancreas provide significant insights into the influence mechanism of air exposure to the trehalose mechanism and immune response in the red claw crayfish.

The effects of ammonia stress exposure on protein degradation, immune response, degradation of nitrogen-containing compounds and energy metabolism of Chinese mitten crab

BACKGROUND

The Chinese mitten crab is one of the most economically important crabs that are widely farmed in China. Ammonia, which is a main physiological challenge for crab culture, grows rapidly in the intensive culture system over time, but little information is available with Chinese mitten crab on the molecular mechanisms.

METHODS AND RESULTS

Therefore, to understand the mechanism of response to ammonia stress in Eriocheir japonica sinensis, comparative transcriptome analysis was used to identify the key genes and pathways in hepatopancreas challenged with ammonia stress (325.07 mg/L NH4Cl). By sequencing the transcriptome hepatopancreas of E. j. sinensis treated with ammonia, 366,472 unigenes were obtained and annotated into several public libraries for later analyses. Subsequently, 1775 differentially expressed genes (DEGs) were identified according to comparative transcriptome analysis, of which 307 were up-regulated and 1468 were down-regulated. According to the DEGs of GO and KEGG enrichment analyses, we focused on four aspects of significant enrichment in this study: protein degradation, immune response, degradation of nitrogen-containing compounds and energy metabolism. The genes involved in protein degradation and energy metabolism process showed a significant decrease which was consisting of overall biological activity of E. j. sinensis decreased. In addition, five genes involved in high concentration of ammonia were discovered and validated by qRT-PCR.

CONCLUSIONS

This study will help us understand the molecular mechanisms of E. j. sinensis under high ammonia exposure and provide valuable information to the future research of other crabs with ammonia exposure.

Identification of Candidate Genes Associated With Hypoxia Tolerance in Trachinotus blochii Using Bulked Segregant Analysis and RNA-Seq

Golden Pompano (Trachinotus blochii) has rapidly developed into the one of the main valuable fish species in Chinese marine aquaculture. Due to its rapid growth, active metabolism, and high oxygen consumption, hypoxia will increase its mortality and cause serious economic losses. We constructed two experimental groups of fish with different degrees of tolerance to hypoxia, used BSR-Seq analysis based on genome and genetic linkage groups to locate SNPs and genes that were related to the differences in hypoxia tolerance. The results showed that hypoxia tolerance SNPs of golden pompano may be jointly determined by multiple linkage groups, especially linkage groups 18 and 22. There were 768 and 348 candidate genes located in the candidate regions of the brain and liver, respectively. These genes were mainly involved in anaerobic energy metabolism, stress response, immune response, waste discharge, and cell death. The prostaglandin-endoperoxide synthase 2 (PTGS2) on LG8, which is involved in the metabolism of arachidonic acid, has a G/A nonsynonymous mutation at position 20641628, and the encoded amino acid was changed from hydrophobic aspartic acid to asparaginate. The specific pathway of the RIG-I-like receptor signaling pathway in the liver may mediate the metabolic system and the immune system, linking glucose metabolism with immune regulation. The death of the hypoxia-intolerant group may be due to the accumulation of lactic acid caused by the activation of anaerobic glycolysis during the early stage of hypoxia stress, and the activation of type I interferon was inhibited, which resulted in decreased immunity. Among the genes involved in the RIG-I-like receptor signaling pathway, the CYLD Lysine 63 Deubiquitinase (CYLD) located on LG16 had a G/T nonsynonymous mutation at position 13629651, and the encoded amino acid was changed from alanine acid to valine. The interferon induced with helicase C domain 1 (Ifih1) located on LG18 has a G/C nonsynonymous mutation at position 16153700, and the encoded hydrophilic glycine was changed to hydrophobic alanine. Our findings suggest these SNPs may assist in the molecular breeding of hypoxia-tolerant golden pompano, and speculate that the balance of glucose and lipid metabolism plays a key role in Trachinotus blochii under acute hypoxia.

SNP Development in Penaeus vannamei via Next-Generation Sequencing and DNA Pool Sequencing

Next-generation sequencing and pool sequencing have been widely used in SNP (single-nucleotide polymorphism) detection and population genetics research; however, there are few reports on SNPs related to the growth of Penaeus vannamei. The purpose of this study was to call SNPs from rapid-growing (RG) and slow-growing (SG) individuals’ transcriptomes and use DNA pool sequencing to assess the reliability of SNPs. Two parameters were applied to detect SNPs. One parameter was the p-values generated using Fisher’s exact test, which were used to calculate the significance of allele frequency differences between RG and SG. The other one was the AFI (minor allele frequency imbalance), which was defined to highlight the fold changes in MAF (minor allele frequency) values between RG and SG. There were 216,015 hypothetical SNPs, which were obtained based on the transcriptome data. Finally, 104 high-quality SNPs and 96,819 low-quality SNPs were predicted. Then, 18 high-quality SNPs and 17 low-quality SNPs were selected to assess the reliability of the detection process. Here, 72.22% (13/18) accuracy was achieved for high-quality SNPs, while only 52.94% (9/17) accuracy was achieved for low-quality SNPs. These SNPs enrich the data for population genetics studies of P. vannamei and may play a role in the development of SNP markers for future breeding studies.

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

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 6360.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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07254-x.

Single nucleotide polymorphism associated with disease resistance in Penaeus vannamei

Despite the considerable number of genetic markers published for Penaeus vannamei, the classification of these markers and their standardization in specific databases is still insufficient. As a consequence, access to these markers is difficult, hampering their application in genetic association studies. In this study, all previously described single nucleotide polymorphisms (SNPs) related to resistance for P. vannamei were revised, and 512 SNPs were identified and classified in detail. We observed that most of the SNPs occurred in the proteins including Toll like receptors 1 and 3, hemocyanin large and small subunits, and anti-lipopolysaccharide factors 1 and 2, allowing to propose to use them as targets in association studies involving resistance in P. vannamei. Additionally, the potential effects of the most frequent non-synonymous coding SNPs in the secondary structure of the main target proteins were evaluated using an in silico approach. These data can serve as the starting point for the development of new genetic and computational tools as well as for the design of new association studies that involve resistance in P. vannamei.

A SNP of the hemocyanin gene (LvHc) is a marker for high growth and ammonia-tolerance in Pacific white shrimp Litopenaeus vannamei

Expression levels of hemocyanin (LvHc), activating transcription factor 4 (LvAtf4), glutathione S-transferase (LvGst), caspase 2 (LvCasp2) and anti-lipopolysaccharide factor (LvAlf) were examined in the hepatopancreas of Pacific white shrimp Litopenaeus vannamei juveniles exposed to a lethal concentration of ammonia-N (32.15 mg/l). The expression levels of all transcripts except LvAlf were significantly greater (P < 0.05) in tolerant shrimp (Lv-AT; N = 30) that survived up to 72 h post treatment (hpt) than in susceptible shrimp (Lv-AS24 and Lv-AS72; N = 45 and 15), that died within 24 h or between 24 and 72 hpt, respectively. Subsequently, effects of non-lethal concentrations of ammonia-N (control, 10 and 20 mg/l) on the expression of LvHc in juvenile shrimp were examined. Compared to the control, expression levels of LvHc transcripts in hemocytes and the hepatopancreas of tested shrimp changed after exposure to ammonia-N. One SNP (C > T₅₄₅) was found in the LvHc₃₂₂ gene segment. Real-time PCR amplification of specific alleles (real-time PASA) was developed for detection of C > T₅₄₅ genotypes. Juveniles in the lethal exposure test that carried a C/T₅₄₅ genotype showed a greater average body weight and total length (8.46 ± 0.36 g and 10.05 ± 0.16 cm) than those with a C/C₅₄₅ genotype (7.48 ± 0.31 g and 9.60 ± 0.13 cm) (P < 0.05). Similar results were found in the second generation (G2) of a growth-improved stock (3 and 4 families of BIOTEC-G2-L1 and BIOTEC-G2-L2) and in commercially farmed shrimp (2 groups). Accordingly, expression levels and SNP of LvHc can serve as markers for selection high growth performance in ammonia-tolerant L. vannamei.

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

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

Molecular mapping of major QTL conferring resistance to orange wheat blossom midge (Sitodiplosis mosellana) in Chinese wheat varieties with selective populations

Two novel midge resistance QTL were mapped to a 4.9-Mb interval on chromosome arm 4AL based on the genetic maps constructed with SNP markers. Orange wheat blossom midge (OWBM) is a devastating insect pest affecting wheat production. In order to detect OWBM resistance genes and quantitative trait loci (QTL) for wheat breeding, two recombinant inbred line (RIL) populations were established and used for molecular mapping. A total of seven QTL were detected on chromosomes 2D, 4A, 4D and 7D, respectively, of which positive alleles were all from the resistant parents except for the QTL on 7D. Two stable QTL (QSm.hbau-4A.2-1 and QSm.hbau-4A.2-2) were detected in both populations with the LOD scores ranging from 5.58 to 29.22 under all three environments, and they explained a combined phenotypic variation of 24.4-44.8%. These two novel QTL were mapped to a 4.9-Mb physical interval. The single-nucleotide polymorphism (SNP) markers AX-109543456, AX-108942696 and AX-110928325 were closely linked to the QTL and could be used for marker-assisted selection (MAS) for OWBM resistance in wheat breeding programs.

Gill transcriptome alterations in Macrobrachium rosenbergii under copper exposure

Copper is one of common contaminants in estuaries and coastal zones, which may cause physiological dysfunction in aquatic organisms. However, molecular response triggered by Cu have remained largely unknown in freshwater prawn Macrobrachium rosenbergii. In the present study, we performed transcriptomic analysis to characterize molecular mechanisms of copper immunotoxicity in gills from M. rosenbergii. A large number of potential simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) loci in the transcriptome were identified. 19,417 and 8989 differentially expressed genes (DEGs) were obtained at 3 h and 48 h after exposure, respectively. Most of these DEGs were down-regulated implying that gene expressions were largely inhibited by Cu, which might lead to impairments of biological functions. Functional enrichment analysis of these DEGs revealed that immune, detoxification and apoptosis were the differentially regulated processes by Cu stress. 12 DGEs involved in immune response and heavy metal detoxification were discovered and validated by qRT-PCR. The results indicated that the M. rosenbergii might counteract the toxicity of Cu at the transcriptomic level by increasing expressions of immune- and heavy metal detoxification-related genes, and these selected genes could be used as molecular indicators for Cu stress. Our study firstly reported the stress response at transcriptional level in M. rosenbergii during Cu exposure. The genes and pathways identified here not only give us new insight into molecular mechanisms underlying Cu toxicity effects in prawn, but facilitate biomarker identification and stress-resistant breeding studies.

The Single-molecule long-read sequencing of Scylla paramamosain

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

Identification of Single Nucleotide Polymorphisms Related to the Resistance Against Acute Hepatopancreatic Necrosis Disease in the Pacific White Shrimp Litopenaeus vannamei by Target Sequencing Approach

Acute hepatopancreatic necrosis disease (AHPND) is a major bacterial disease in Pacific white shrimp Litopenaeus vannamei farming, which is caused by Vibrio parahaemolyticus. AHPND has led to a significant reduction of shrimp output since its outbreak. Selective breeding of disease-resistant broodstock is regarded as a key strategy in solving the disease problem. Understanding the relationship between genetic variance and AHPND resistance is the basis for marker-assisted selection in shrimp. The purpose of this study was to identify single nucleotide polymorphisms (SNPs) associated with the resistance against AHPND in L. vannamei. In this work, two independent populations were used for V. parahaemolyticus challenge and the resistant or susceptible shrimp were evaluated according to the survival time after Vibrio infection. The above two populations were genotyped separately by a SNP panel designed based on the target sequencing platform using a pooling strategy. The SNP panel contained 508 amplicons from DNA fragments distributed evenly along the genome and some immune-related genes of L. vannamei. By analyzing the allele frequency in the resistant and susceptible groups, 30 SNPs were found to be significantly associated with the resistance of the shrimp against V. parahaemolyticus infection (false discovery rate corrected at P < 0.05). Three SNPs were further validated by individual genotyping in all samples of population 1. Our study illustrated that target sequencing and pooling sequencing were effective in identifying the markers associated with economic traits, and the SNPs identified in this study could be used as molecular markers for breeding disease-resistant shrimp.

Identification of SNPs potentially related to immune responses and growth performance in Litopenaeus vannamei by RNA-seq analyses

Litopenaeus vannamei is one of the most important shrimp species for worldwide aquaculture. Despite this, little genomic information is available for this penaeid and other closely related taxonomic crustaceans. Consequently, genes, proteins and their respective polymorphisms are poorly known for these species. In this work, we used the RNA sequencing technology (RNA-seq) in L. vannamei shrimp evaluated for growth performance, and exposed to the White Spot Syndrome Virus (WSSV), in order to investigate the presence of Single Nucleotide Polymorphisms (SNPs) within genes related to innate immunity and growth, both features of great interest for aquaculture activity. We analyzed individuals with higher and lower growth rates; and infected (unhealthy) and non-infected (healthy), after exposure to WSSV. Approximately 7,000 SNPs were detected in the samples evaluated for growth, being 3,186 and 3,978 exclusive for individuals with higher and lower growth rates, respectively. In the animals exposed to WSSV we found about 16,300 unique SNPs, in which 9,338 were specific to non-infected shrimp, and 7,008 were exclusive to individuals infected with WSSV and symptomatic. In total, we describe 4,312 unigenes containing SNPs. About 60% of these unigenes returned GO blastX hits for Biological Process, Molecular Function and Cellular Component ontologies. We identified 512 KEGG unique KOs distributed among 275 pathways, elucidating the majority of metabolism roles related to high protein metabolism, growth and immunity. These polymorphisms are all located in coding regions, and certainly can be applied in further studies involving phenotype expression of complex traits, such as growth and immunity. Overall, the set of variants raised herein enriches the genomic databases available for shrimp, given that SNPs originated from nextgen are still rare for this relevant crustacean group, despite their huge potential of use in genomic selection approaches.