- Research
- Open access
- Published:
Genome-wide identification of ATG genes and their expression profiles under biotic and abiotic stresses in Fenneropenaeus chinensis
BMC Genomics volume 25, Article number: 625 (2024)
Abstract
Background
Autophagy is a conserved catabolic process in eukaryotes that contributes to cell survival in response to multiple stresses and is important for organism fitness. Extensive research has shown that autophagy plays a pivotal role in both viral infection and replication processes. Despite the increasing research dedicated to autophagy, investigations into shrimp autophagy are relatively scarce.
Results
Based on three different methods, a total of 20 members of the ATGs were identified from F. chinensis, all of which contained an autophagy domain. These genes were divided into 18 subfamilies based on their different C-terminal domains, and were found to be located on 16 chromosomes. Quantitative real-time PCR (qRT-PCR) results showed that ATG genes were extensively distributed in all the tested tissues, with the highest expression levels were detected in muscle and eyestalk. To clarify the comprehensive roles of ATG genes upon biotic and abiotic stresses, we examined their expression patterns. The expression levels of multiple ATGs showed an initial increase followed by a decrease, with the highest expression levels observed at 6 h and/or 24 h after WSSV injection. The expression levels of three genes (ATG1, ATG3, and ATG4B) gradually increased until 60 h after injection. Under low-salt conditions, 12 ATG genes were significantly induced, and their transcription abundance peaked at 96 h after treatment.
Conclusions
These results suggested that ATG genes may have significant roles in responding to various environmental stressors. Overall, this study provides a thorough characterization and expression analysis of ATG genes in F. chinensis, laying a strong foundation for further functional studies and promising potential in innate immunity.
Background
Fenneropenaeus chinensis, commonly known as Chinese shrimp, is a highly valued species in China’s aquaculture industry. Its distribution is mainly concentrated in the Yellow Sea and Bohai Sea of China, as well as the western and southern coast of the Korean Peninsula [1]. Due to the advancement of aquaculture techniques, F. chinensis emerged as the most important shrimp species for cultivation in China during the 1990s [2]. F. chinensis and other crustaceans often encounter changes in environmental factors, such as extreme temperatures, salinity, hypoxia and abnormal acid-base levels. Environmental factors can disrupt homeostatic equilibrium and cause the fluctuations in the neuroendocrine, physiological, and behavioral status of aquatic animals, negatively impacting their health [3]. Saline-alkaline water aquaculture is becoming a promising solution to accommodate the growing needs of the aquaculture industry. Therefore, it is important to investigate the physiological changes and molecular responses in aquatic animals adapting to this environment. F. chinensis is highly sensitive to changes in salinity levels, particularly in low-salt environments, which is considered a major stressor. It is essential to understand the effects of low salinity on F. chinensis to develop effective mitigation strategies in F. chinensis aquaculture. The sensitivity to low salinity is rooted in the osmoregulatory mechanisms employed by these organisms to maintain internal salt and water balance. The balance is crucial for their metabolic processes and impacts the reproductive and developmental processes of shrimps. In addition to abiotic stress, the significant environmental variations were accompanied by severe bacterial and viral infections, resulting in an epizootic breakout that caused great losses to the aquaculture industry [4]. White spot syndrome virus (WSSV) is a major pathogen in shrimp aquaculture. It is a double-stranded DNA virus that causes white spot disease [5,6,7]. The mortality rate of F. chinensis infected with WSSV can reach 100% within a week. WSSV first emerged in the late 20th century, causing huge economic losses. Despite extensive research, no effective solution to the WSSV problem has been conducted, and it remains one of the most fatal pathogens in F. chinensis aquaculture [8].
Autophagy is a conserved intracellular degradation system found in eukaryotes. In autophagy, an isolation membrane emerges suddenly in the cytoplasm, which then expands and transforms into double membrane-bound structure called an autophagosome [9,10,11]. During this process, a portion of cytoplasm, including proteins and organelles, is sequestered into the autophagosome. The autophagosome then fuses with a lysosome (or a vacuole, in yeast and plants). The inner membrane, referred to as an autophagic body in yeast, is exposed to lysosomal hydrolases and degraded along with its contents [9,10,11]. The primary function of autophagy is to maintain cellular homeostasis by recycling intracellular materials, which promotes stress resistance and longevity [11]. Autophagy has emerged as a crucial cellular process with implications for health and immunity in all eukaryotic organisms. Research has shown that a range of autophagy-related genes play important roles in both biotic and abiotic stresses in plants and animals, participating in sugar metabolism, antimicrobial peptide regulation, and other pathways [12, 13]. In yeast, there are over 36 autophagy-related genes (ATG) have been identified. Most of these genes have corresponding homologous genes in mammals [14, 15]. Drought and salt stress cause ion stress which induces oxidative damage in plant cells. Autophagy repairs this damage under the activation of various biological factors [16]. Silencing OsATG2 and OsATG7 inhibits autophagy and reduces wheat salt tolerance [17]. Overexpression of MdATG18a enhances the adaptability of apples to drought stress [18]. In animal studies, it has been found that the ATG5-ATG12 conjugate functions in the mouse innate antiviral immune system, enhancing their bactericidal activity [19]. For the banana shrimp, interference with either ATG3 or ATG6 during WSSV challenge resulted in a decrease in autophagic levels [20]. Therefore, autophagy-related research has recently gained increasing attention from the scientific community, but it is relatively uncommon in crustaceans. Given the growing importance of shrimp aquaculture and the risks posed by viral infections and environmental stress, it is crucial to comprehend the autophagy-related mechanisms in this economically vital species. However, the precise details of the F. chinensis autophagy-related gene family, its structure, and its specific roles in the response to WSSV infection and low-salt stress have not been comprehensively elucidated.
The objective of this study is to identify and characterize the ATG gene family in F. chinensis and investigate their expression patterns in response to biotic and abiotic stresses. The study aims to contribute to the advancement of shrimp aquaculture and the management of viral diseases and environmental factors that pose a threat to this industry.
Results
Identification and characterization of ATG genes in F. chinensis
In the F. chinensis, a total of 20 ATG genes were identified and named according to the rules for ATG gene nomenclature (Table 1). These ATG genes (designated as FcATGs) were classified into18 subfamilies, including FcATG1-3, FcATG4B, FcATG4D, FcATG5-8, FcATG8B, FcATG9-10, FcATG12-14, FcVPS15, FcATG16 -18, and FcATG101. The predicted molecular weight of F. chinensis ATG proteins were in the range of 13.63 kDa to 241.9 kDa and the deduced isoelectric points ranged from 4.57 to 9.69. The FcATG12 gene had the smallest molecular weight of 13.63 kDa, while the FcATG2 gene had the largest molecular weight of 241.9 kDa. The instability index of three ATG proteins, namely FcATG12, FcATG14, and FcATG 101, was less than 40, indicating their stability. On the other hand, the remaining 17 proteins were predicted to be unstable. All ATG proteins had GRAVY values less than 0, indicating their hydrophilic nature. The open reading frames (ORFs) of ATG genes ranged from 360 to 6,546 bp, with predicted protein lengths ranging from 119 to 2181 amino acids (aa). Subcellular location prediction revealed that 12 ATG proteins are located in the nucleus, while the remaining proteins are mainly present in the Golgi apparatus, cytoplasm, and mitochondrion. The cDNA sequences of these ATG genes have been submitted to the GenBank database, and their characteristics are summarized in Table 1.
Phylogenetic tree construction
To explore the evolutionary relationships of FcATG genes, we constructed a phylogenetic tree using ATG protein sequences from fifteen species, including F. chinensis, Mus musculus, Penaeus monodon, Litopenaeus vannamei, Macrobrachium nipponense, Bombyx mori, and Drosophila melanogaster et al. (Fig. 1). The phylogenetic analysis revealed that the maximum likelihood (ML) tree for ATGs split into 18 subfamilies, with the subfamilies in F. chinensis clustering with their corresponding counterparts from other species as expected. The FcATG1, FcATG7, FcATG13, FcATG14, FcATG16, FcATG18, and FcATG101 subfamilies clustered first, followed by the FcATG2, FcATG3, FcATG4, FcATG6, FcATG10, FcATG12, and FcATG15 subfamilies, as well as the FcATG5, FcATG8, FcATG9, and FcATG12 subfamily. These results are consistent with the classification described in previous research and provide insights into the evolution of FcATGs orthologous genes in different species.
Chromosomal distribution and structure analysis
Chromosomal distribution analysis revealed that 20 ATG genes are distributed across 16 shrimp chromosomes, from NC-061821.1 to NC-061859.1. However, the distribution of ATG genes was uneven on each chromosome. For instance, FcATG3 and FcATG14, FcATG8 and FcATG17, as well as FcATG5 and FcATG4B are all located on the same chromosome, while the remaining chromosomes each hosted only a single ATG gene (Fig. 2). The non-random distribution of ATG genes in the shrimp genome highlights the importance of their genomic arrangement in autophagy-related processes.
Based on the structural analysis of the F. chinensis ATG genes, the number of exons among ATG members varies from 1 to 29. Each ATG member comprises both UTR and CDS regions, with coding regions of similar lengths for members within the same subfamily, such as FcATG4B, FcATG4D, FcATG8, and FcATG8B (Fig. 3a). Protein domain analysis reveals that FcATG6, FcVPS15, and FcATG16 each possess two domains: APG6_N and APG6 for FcATG6, VPS15 and WD40 for FcVPS15, and FcATG16 and WD40 for FcATG16. In contrast, other proteins have only one domain. Different subfamilies exhibit distinct domain compositions. However, within the same subfamily, the domain structure is generally conserved. For example, both FcATG4B and FcATG4D contain the Peptidase_C54 domain, while FcATG8 and FcATG8B contain the Ubl_ATG8 domain (see Fig. 3b). This structural analysis offers insight into the diversity and conservation of gene architecture within the FcATG family and sheds light on potential functional implications.
Expression profiles of ATG genes in F. chinensis
The relative expression levels of FcATG were measured in various tissues, including the eyestalk, gill, heart, hepatopancreas, intestine, muscle, and stomach, using qRT-PCR with 18 S rDNA as the internal control in nine untreated shrimps. As shown in Fig. 4, FcATG1, FcATG4B, FcATG7, and FcATG12 genes exhibited higher expression levels in the eyestalk compared to the gills, heart, intestine, hepatopancreas, and stomach. The other 14 FcATG genes had the highest level of expression in muscles.
Expression profiles of ATG genes after WSSV infection
This experiment utilized qRT-PCR analysis to study the expression levels of the 20 identified members of the FcATG gene family. The cDNA samples from the hepatopancreas were collected at 3 h, 6 h, 12 h, 24 h, 36 h, 48 h, and 60 h after WSSV injection to investigate the expression levels of each member during WSSV infection. Figure 5 shows that FcATG1, FcATG3, and FcATG4B had a significant upregulation trend, with the highest expression at 60 h. FcATG7 and FcATG10 exhibited a noticeable downregulation trend. Most ATG genes showed an up-regulation followed by a down-regulation trend, with FcATG2, FcATG8, and FcATG8B peaked at 3 h. FcATG5, FcATG6, FcATG9A, FcATG13, FcATG14, FcVPS15, FcATG17, and FcATG101 reached their peak expression at 6 h. FcATG4D, FcATG12, FcATG16, and FcATG18 reached their highest expression at 24 h.
Expression profiles of ATG genes after low-salt stress
The expression of FcATG genes in hepatopancreas tissues after low-salt stress was analyzed using qRT-PCR (Fig. 6). The results showed an up-regulation trend for FcATG1-3, FcATG4B, FcATG4D, FcATG5, FcATG9, FcATG12-13 FcVPS15, FcATG16, and FcATG8, peaking at 72–96 h. FcATG10 exhibited a notable upregulation, reaching peak expression at 24 h. Meanwhile, FcATG6- FcATG8, FcATG8B, FcATG14, FcATG17, FcATG101 initially showed a downregulation followed by an upregulation trend. FcATG17, and FcATG101 peaked at 72 h, while the rest reached their highest expression at 96 h.
Protein interaction analysis
To explore the intricate interaction networks of ATGs in F. chinensis, putative interaction networks were constructed using the online STRING server based on orthologous genes in D. melanogaster (Fig. 7). The results revealed that within the ATG gene family, 20 ATG proteins displayed mutual interactions. Among these, proteins such as Pi3K59F, Tango5 and CG42554 exhibited high connectivity to multiple ATG genes, indicating a close predicted interaction with ATG proteins in F. chinensis.
Discussion
Autophagy is a highly conserved catabolic pathway that is involved in the cellular degradation of long-lived proteins or dysfunctional cellular components through lysosomes action. It is considered a pro-survival mechanism and for maintaining homeostasis [21]. Meanwhile, autophagy is a sensitive process that underlies cell response to almost every stressful condition affecting cellular homeostasis. Its role is complex and likely depends on the cell’s genetic background and environmental cues [22, 23]. Recent research has shown that autophagy is crucial for maintaining homeostasis at both the cellular and organismal levels [24,25,26]. While the study of ATG genes has been characterized in yeasts and mammals for decades, there has been a lack of systematic identification and analysis of the ATG gene family in F. chinensis. Additionally, research on ATGs under biotic and abiotic stresses of F. chinensis has not been conducted. This study identified and characterized the ATG genes of F. chinensis, including their phylogenetic analysis, protein structures, and physicochemical properties of ATGs. In addition, the expression profiles of ATG genes were analyzed after WSSV infection, low salt stress, and in healthy tissues. These analyses provide insights into the involvement of ATGs in response to toxicological and environmental stresses.
This study reports the identification of 20 ATG family members in the F. chinensis genome, three of which have instability coefficients below 40, indicating relative stability of the gene family. This finding is consistent with the results reported by Liu et al. [27] in their investigation of the physical and chemical properties of the ATG gene in Punica granatum. Meanwhile, the 20 ATG gene family members of F. chinensis are distributed unevenly across 16 chromosomes. This distribution pattern is similar to that found in Arabidopsis [28], suggesting that even within the same family, members can have different chromosomal positions. The ATG genes are relatively conserved from yeast to humans, but the detailed evolutionary history remains unclear. Analyzing phylogenetic relationships can provide insight into the evolution of the ATG genes. A phylogenetic tree was constructed based on the amino acid sequences of ATG homologous genes from 15 different species. The resulting tree revealed that 20 shrimp ATG genes were grouped into 18 subfamilies. Moreover, the evolutionary relationship between L. vannamei and the ATG genes of the Chinese shrimp ATG family was found to be similar, suggesting a close association. As per Zhu’s findings, the ATG gene of Procambarus clarkii and Homarus americanus are closely related [26]. In general, proteins within the same subfamily on a branch exhibit similar structures and functional domains, thereby indicating a tendency towards functional conservation.
Previous studies have profiled the expression patterns of ATG family members in crustacean [20, 29,30,31]. However, such research has not been conducted in F. chinensis. To better understand the characteristics and functions of ATG genes in F. chinensis, we used qRT-PCR to profile the expression patterns of ATG genes in different tissues. ATGs are expressed ubiquitously in various F. chinensis tissues with varying levels, indicating their divergent functions in the organism. Sixteen ATG genes have the highest expression levels in muscles, while four ATG genes have the highest expression levels in the eye stalk, suggesting that ATG genes may be involved in the growth and development process. These finding are consistent with the research results in P. monodon, L. vannamei and Ctenopharyngodon idellus [20, 31, 32]. Intriguingly, the expression level of FcATG16 was notably higher in the gill compared to other ATG family members. The gill is a multipurpose organ that provides gas exchange and osmotic regulation. This suggests that FcATG16 may play an important role in the physiological function of the gill [20, 30]. Besides, ATG3 is highly expressed in hepatopancreas, which is an important immune organ, suggesting ATG3 may be involved in immune responses, lipid metabolism, and detoxification in F. chinensis, as previous studies have demonstrated in other species such as Pelteobagrus fulvidraco, Macrobrachium rosenbergii, and some mammals [33,34,35,36]. Taken together, ATG genes may participate in various biological processes, with some members specifically related to the immune system of F. chinensis.
Studies have shown that the expression of ATGs can be affected by various stimuli, such as environmental stress and pathogen invasion, and these genes play important roles in mediating autophagy, participating in embryonic development, host resistance to viral and bacterial infection, and immunity process [37, 38]. Recently, there have been several efforts to study the role of the ATG gene family in pathogen invasion. The results have shown that different ATG members are stimulated after virus infection in various species. For instance, in olive flounder Paralichthys olivaceus, the mRNA level of ATG6 significantly increased after viral hemorrhagic septicemia virus (VHSV) infection [39]. In Procambarus clarkii, ATG14 expression was initially upregulated upon WSSV infection and then stabilized [26]. Due to the limited information on ATG in crustaceans, only the expression patterns of individual members in response to virus invasion were detected. Therefore, we conducted WSSV infection experiments and measured the expression levels of all FcATG members in F. chinensis. The results showed that in the hepatopancreas, 12 FcATG genes exhibited an initial up-regulation followed by a down-regulation trend. Among them, ATG6, ATG9A, and ATG14 genes underwent significant changes, with peak expression levels occurring at either 6–24 h post-treatment. This is consistent with the changes in the ATG6 gene observed in P. vannamei’s responds to WSSV infection [20], suggesting that FcATG6 and FcATG14 play an crucial role in F. chinensis’s response to biotic stress. Unlike other ATG family members, the expression levels of three genes (FcATG1, FcATG3, and FcATG4B) gradually increased until 60 h after injection. This indicated that these genes can be significantly induced by WSSV infection. It is believed that autophagy-related genes may be regulated thus cellular autophagy be affected in the early stage of viral infection. In the late stage of infection, host cells downregulate the expression of ATG to resist the virus and prevent its proliferation [40]. In contrast, the expression trend of FcATG7 and FcATG10 was continuously down-regulated after WSSV infection. This suggests that these genes may have species-specific or environment-dependent effects on the response of F. chinensis to WSSV infection. On the other hand, the expression trends of FcATG8 and FcATG101 did not change significantly. It is speculated that they may not be the main genes involved in coping with WSSV in F. chinensis, and further research is need to confirm this hypothesis. In summary, the results of F. chinensis suggest that ATGs play a crucial rule in the immune response against WSSV, emphasizing their significance in immunity.
Further, autophagy is essential for maintaining the intracellular stability under normal and stress conditions and it ensures quality control within cells. However, autophagy induced by ion concentrations, pH, and other stressors, particularly when the stress level is not fatal, constitutes a strategy to adapt and cope with stress. This can promote cell survival by maintaining a sufficient amino acid pool and cell energy level [41, 42]. For example, L. vannamei and Macrobrachium japonicum are able to regulate autophagy in response to abiotic stresses such as low temperature and hypoxia [31, 43]. In this study, 12 FcATG genes were significantly induced, and their transcript abundance peaked at 96 h after treatment under low-salt conditions. The expression levels of FcATG5 and FcATG12 were significantly upregulated, suggesting that these FcATG genes are activated to promote cell survival by regulating autophagy when challenged by low-salt conditions. Similarly, the expression levels of ATG5 and ATG12 were also significantly upregulated under pH or carbonate alkalinity stress [44], indicating that the FcATG5 and FcATG12 genes may play an crucial role in the adaptation of F. chinensis to abiotic stress. Previous studies have shown that F. chinensis increases body energy consumption and enhances sugar metabolism activity in response to abiotic stress, including the acute changes in salinity [45]. Autophagy is proposed to participate in the process of sugar metabolism by regulating glucose uptake, key enzymes of glycolysis, and mitochondria [46]. Taken together, these results demonstrate that abiotic stress alters the expression of ATG genes in organisms and plays an important role in their response to such stress.
In addition, the connectivity of ATG genes in F. chinensis was firstly analyzed via PPI network analysis to further investigate the regulatory mechanisms of these ATG genes. The ATG genes exhibited high connectivity with genes related to the PI3K pathway, which has now been demonstrated to be involved in numerous important physiological processes and immune responses in aquatic animals. PI3Ks can regulate NK cytotoxic activity and inhibit NADPH oxidase [47]. Following a Vibrio infection, the expression level of PI3K in Haliotis disversicolor was found to be significantly elevated in various tissues [48]. Similarly, Scophthalmus maximus employed PI3K signaling pathway to regulate ion transport in response to low-salt stress [49]. According to our previous experiments, the induced upregulation of PI3K-related genes responded to high alkali and elevated pH stresses by activating the PI3K pathway in F. chinensis. Therefore, we hypothesize that the ATG genes’s response to WSSV injection and low-salt stress may also be mediated through the PI3K pathway in F. chinensis. This study will enhance our understanding of the molecular mechanisms underlying how FcATGs respond to various environmental stressors.
Conclusions
In summary, we identified 20 ATG genes in F. chinensis, and characterized their structure, domains, phylogenetic tree, chromosomal distribution, and bioinformatics information. The qPCR examination results indicated that ATG genes played an essential role in WSSV infection. Furthermore, we analyzed their expression profiles based on generated from low-salt stress, providing potential functional information on the response of ATG genes to abiotic stressors in F. chinensis. This study enhances our comprehension of the molecular basis of the ATG gene family’s response to toxicological and environmental stresses in F. chinensis. It provides important clues for future research on their functions.
Materials and methods
Experimental Shrimp
In the present study, F. chinensis (weight 15.4 ± 2.1 g, length 11.1 ± 0.87 cm) were obtained from the Changyi Haifeng Aquaculture Co., Ltd., Shandong, China. The shrimps were cultured in cement breeding ponds, with continuous aeration, at a temperature of 25 ℃, salinity of 30, and a pH of 8.4. During the acclimation period, the shrimps were fed thrice daily at 00:00, 09:00, and17:00 for two weeks. After acclimation, nine healthy and energetic F. chinensis were randomly selected for the study. Seven tissues, including eye stalks, gills, heart, hepatopancreas, muscles, stomach, and intestines, were collected and stored in liquid nitrogen for later use.
Identification of ATG gene family in F. chinensis
The ATG gene sequences of humans, mouse, zebrafish and arthropods were obtained from the UniProt (https://www.uniprot.org/) and NCBI (http://www.ncbi.nlm.nih.gov/) databases, and were used as the query sequences to search against the whole genome databases of F. chinensis using the TBLASTN alignment tool (E-value = 1e− 5). The amino acid sequences were predicted and translated from the open reading frame (ORF) using ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/). The results were validated by BLASTP against the NCBI non-redundant protein (NR) database. The domain architectures of the ATG in F. chinensis were examined using SMART to confirm all of the family genes. The functional interaction network of ATG genes in F. chinensis was investigated by utilizing the STRING online server (https://cn.string-db.org/) based on orthologous genes identified in D. melanogaster. The resulting network was then visualized using Cytoscape v 3.7.
Phylogenetic analysis
The amino acid sequences of ATG subunits from F. chinensis and several representative animals were obtained from NCBI and aligned using the ClustalW program. Phylogenetic trees were constructed using the Maximum Likelihood (ML) method with a bootstrap value of 1000 in MEGA 11. The resulting tree was annotated using iTOL (https://itol.embl.de/login.cgi).
Bioinformatics analysis
To analyse the characteristics of ATG, we predicted the amino acids, molecular weight (MW), theoretical isoelectric points (pI), instability index, and grand average of hydropathicity (GRAVY) of the ATG genes using ExPASy ProtParam (https://web.expasy.org/protparam/). We used SignalP 5.0 [50] to predict signal peptides and WOLF PSORT [51] to predict subcellular localization. NCBI-CDD (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi) was used to predict the conserved domain of ATG genes. The structures of ATG genes were visualized using TBtools.
The WSSV infection experiment
Preparation of WSSV stock solution: the carapace and hepatopancreas were removed from WSSV - infected shrimp before shearing in an ice - bath environment, and 1 g tissue was collected and mixed with 10 ml phosphate - buffered saline (PBS). The mixture was homogenized at 10,000 rpm for 5 s and then centrifuged at 8000 rpm for 5 min at 4 °C. The supernatant fluid was filtered through a 0.45 μm microporous membrane as a virus stock solution [52]. 103-diluted-stock-solution (using pre-chilled PBS) was used in this experiment to infect healthy shrimps. Shrimps sampled from the shrimp ponds were verified to be WSSV-free after being tested for related pathogens before experimentation.
After acclimation, a total of 180 healthy F. chinensis were selected for the WSSV infection experiment. The shrimps were randomly divided into two groups, with three replicates in each group and 30 shrimps in each replicate. The second abdominal segment muscle of the shrimps was injected with PBS or WSSV suspensions. The experimental group received 10 µL of WSSV virus suspension, while the control group was injected with 10 µL sterile phosphate (PBS). The injected shrimps were raised separately in aquariums and all showed symptoms of WSSV infection. We randomly selected hepatopancreas tissues from three healthy and complete shrimps were randomly selected at 0, 3, 6, 12, 24, 36, 48, and 60 h after injection and frozen immediately in liquid nitrogen and then stored at − 80 ℃.
The low salinity stress experiment
Ninety healthy and energetic shrimps were randomly selected and divided into three groups, each containing 30 shrimps. The stress treatment group was subjected to low-salt stress conditions with a salinity of 15. Salinity was corrected every 6 h using with light brine during the experiment. The hepatopancreas tissues were collected after exposure to low-salinity stress at 0, 3, 6, 12, 24, 48, 72, and 96 h. Three parallel samples were taken at each time point, frozen in liquid nitrogen and then stored at − 80 ℃ for later use.
Total RNA extraction, and cDNA synthesis
The total RNA of tissues was extracted using TransZol Up (TransGen Biotech, China) following the manufacturer’s instructions. RNA concentration and integrity were measured using an Ultra-trace UV spectrophotometer (Thermo, USA) and 1.5% agarose gel electrophoresis (AGE). According to the manufacturer’s instructions, the first-strand cDNA synthesis was performed with HiScript III RT SuperMix for qPCR (+ gDNA wiper) (Vazyme, China).
Quantitative real-time PCR (qRT-PCR) analysis
The expressions of ATGs in different tissues of healthy F. chinensis, as well as their temporal expression patterns in the hepatopancreas after WSSV infection and low-salt stress were detected by qRT - PCR. Primers were designed for each gene using Primer5 software, avoiding regions with hairpin structures as identified by Mfold at 60 °C [53]. The 18 S rRNA was selected as the reference gene. Primer specificity and efficiency were checked for all qPCR conditions. The primer sequences are listed in Table 2. The qRT-PCR assay was conducted using ChamQ SYBR Color qPCR master mix (Vazyme, China) in a final volume of 10 µl on the FAST 7500 Real-time system (FAST, USA). The reaction mixture contained ChamQ SYBR Color qPCR Master Mix (Low ROX Premixed) (Vazyme, China), forward and reverse primers (final concentration 100 nM) and 4 µl of diluted cDNA. The PCR reactions were initiated with a denaturation step at 95 °C for 30 s, followed by 40 cycles of two-step amplification as follows: denaturation at 95 °C for 10 s, annealing at 60 °C for 30 s. Fluorescence data were acquired during the final step. Gene-specific amplification was confirmed by a single peak in the melting curve analysis. Three technical repeats were performed for each biological repeat. The data was analysed using FAST 7500 software, and the relative expression ratio was calculated using 2−ΔΔCT method [54]. Statistical analysis was performed using SPSS 19.0. A One-Way ANOVA was exerted to compare the differences between groups. Finally, GraphPad Prism was used for scientific graphing.
Data availability
The sequence information of Fenneropenaeus chinensis ATG family genes were collected from Fenneropenaeus chinensis genome to Breeding Database (https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_019202785.1/), and the ATG protein sequences of Mus musculus, Litopenaeus vannamei, Macrobrachium nipponense, Bombyx mori, and Drosophila melanogaster et al. were downloaded from the NCBI website (https://www.ncbi.nlm.nih.gov/). All data used during the current study are included in this published article and its supplementary information files or available from the corresponding author on reasonable request.
References
Wang MZ, Kong J, Meng XH, Luan S, Luo K, Sui J, Chen BL, Cao JW, Shi XL, Zhang Q. Evaluation of genetic parameters for growth and cold tolerance traits in Fenneropenaeus chinensis juveniles. PLoS ONE. 2017;12(8):e0183801.
Wang Q, Ren XY, Liu P, Li JT, Lv JJ, Wang JJ, Zhang HE, Wei W, Zhou YX, He YY. Improved genome assembly of Chinese shrimp (Fenneropenaeus chinensis) suggests adaptation to the environment during evolution and domestication. Mol Ecol Resour. 2021;22(1):334–44.
He YY, Wang Q, Li J, Li ZX. Comparative proteomic profiling in Chinese shrimp Fenneropenaeus chinensis under low pH stress. Fish Shellfish Immunol. 2021;120:526–35.
Le Moullac G, Haffner P. Environmental factors affecting immune responses in Crustacea. Aquaculture. 2000;191(1):121–31.
Huang J, Cai SL, Song XL, Wang CM, Yu J, Yang CH. Study on artificial infection for penaeus chinensis by the pathogen of the explosive epidemic disease of shrimp. Mar Fisheries Reseach. 1995;16(1):51–8.
Lo CF, HO CH, peng SE, Chen CH, Chen H, Lin Y, Chang CF, Fu K, Su MS, Wang CH, et al. White spot syndrome baculovirus (WSBV) detected in cultured and captured shrimp, crabs and other arthropods. Dis Aquat Organ. 1996;27(3):215–25.
Van Hulten MCW, Witteveldt J, Peters S, Kloosterboer N, Tarchini R, Fiers M, Sandbrink H, Lankhorst RK, Vlak JM. The white spot syndrome virus DNA genome sequence. Virology. 2001;286(1):7–22.
Li XP, Luan S, Luo K, Cao BX, Chen BL, Kong J, Meng XH. Comparative transcriptomic analysis of Chinese shrimp Fenneropenaeus chinensis infected with white spot syndrome virus. Aquaculture Rep. 2022;22:100986.
Klionsky DJ, Ohsumi Y. Vacuolar import of proteins and organelles from the cytoplasm. Annual Rev Cell Dev Biology. 1999;15(1):1–32.
Mizushima N, Yoshimori T, Ohsumi Y. The role of ATG proteins in autophagosome formation. Annu Rev Cell Dev Biol. 2011;27:107–32.
Noda NN, Fujioka Y. Atg1 family kinases in autophagy initiation. Cell Mol Life Sci. 2015;72(16):3083–96.
He YY, Li ZX, Zhang HE, Hu S, Wang QY, Li J. Genome-wide identification of Chinese shrimp (Fenneropenaeus chinensis) microRNA responsive to low pH stress by deep sequencing. Cell Stress Chaperones. 2019;24(4):689–95.
Yang W, Liu C, Xu QS, Qu C, Sun J, Huang S, Kong N, Lv XJ, Liu ZQ, Wang LL et al. Beclin-1 is involved in the regulation of antimicrobial peptides expression in Chinese mitten crab Eriocheir sinensis. Fish & shellfish immunology. 2019;89:207–216.
Garcia Maurino S, Alcaide A, Dominguez C. Pharmacological control of autophagy: therapeutic perspectives in inflammatory bowel disease and colorectal cancer. Curr Pharm Design. 2012;18(26):3853–73.
Qiang L, Sample A, Shea CR, Soltani K, Macleod KF, He Y. Autophagy gene ATG7 regulates ultraviolet radiation-induced inflammation and skin tumorigenesis. Autophagy. 2017;13(12):2086–103.
Chen H, Dong JL, Wang T. Autophagy in plant abiotic stress management. Int J Mol Sci. 2021;22(8):954–63.
Yu JY, Wang YJ, Jiao JL, Wang HZ. Silencing of ATG2 and ATG7 promotes programmed cell death in wheat via inhibition of autophagy under salt stress. Ecotoxicol Environ Saf. 2021;225:112761.
Sun X, Wang P, Jia X, Huo L, Che RM, Ma F. Improvement of drought tolerance by overexpressing MdATG18a is mediated by modified antioxidant system and activated autophagy in transgenic apple. Plant Biotechnol. 2018;16(2):545–57.
Fumihiko T, Kouji K, Atsushi M, Nao J, Kenji O. The non-canonical role of Atg family members as suppressors of innate antiviral immune signaling. Autophagy. 2008;4(1):67–9.
Chu JX, Meng FJ, Zhang GC. Cloning of ATG6 in Litopenaeus vannamei and research on its expression under WSSV infection. J Anhui Agricultural Sci. 2018;46(12):107–11.
J.Klionsky RPK. An overview of autophagy: morphology, mechanism, and regulation. Mary Ann Liebert. 2014;20(3):460–73.
Lerner E, Kimchi. The paradox of autophagy and its implication in cancer etiology and therapy. Apoptosis: Int J Program cell Death. 2009;14(4):376–91.
Galati S, Boni C, Gerra MC, Lazzaretti M, Buschini A. Autophagy: a player in response to oxidative stress and DNA damage. Oxidative Med Cell Longev. 2019;2019:5692958.
Brittany E, Nimrod M, Chao M. Impaired autophagy and defective mitochondrial function:converging paths on the road to motor neuron degeneration. Front Cell Neurosci. 2016;10(44):00044.
Jacob JA, Mohammad JM, Salmani, Jiang Z, Feng L, Song J, Jia X, Chen B. Autophagy: an overview and its roles in cancer and obesity. Clin Chim Acta. 2017;468:85–9.
Zhu MR, Zhan M, Xi CJ, Gong J, Shen HS. Molecular characterization and expression of the autophagy-related gene Atg14 in WSSV-infected Procambarus clarkii. Fish Shellfish Immunol. 2022;125:200–11.
Liu C, Zhao YJ, Zhao X, Dong J, Yuan Z. Identification of ATG gene family of Punica granatum and analysis on their expression pattern under abiotic stress. J Plant Resour Environ. 2022;31(5):37–49.
Zhang HZ, Zhu L, Zuo ZY, Liu WH, Xu J. Prediction of the function of autophagy-related genes (ATGs) in development and abiotic stress based on expression profiling in Arabidopsi. Genomics Appl Biology. 2020;39(6):2671–82.
Sun SM, Chen YX, Zheng C, Zhao QQ, Xue C. Cloning and expression analysis of autophagy genes ATG13 and ATG101 in Macrobrachium nipponense under hypoxic stress. J Fisheries China. 2021;45(6):846–61.
Yin WJ, Liu F, Meng FJ, Zhang GC, Wang LY, Sun JS. Cloning of ATG3 in Litopenaeus vannamei and research on its expression under WSSV infection. J Tianjin Normal University(Natural Sci Edition). 2021;41(2):51–5.
Liu W. Mechanism of miR-7562 mediated ATG5-ATG12 conjugation system in response to Vibrio harveyi stress in Penaeus monodon. Shanghai, China: Shanghai Ocean University; 2018.
He L, Ruan JM, Liu Y, Fu JP, Liu L, Wei LL. Cloning of beclin1, an autophagy gene, and it’s expression under microcystinlr stress in grass carp (Ctenopharygodon Idella). Acta Hydrobiol Sin. 2019;43(3):479–85.
Wei CC, Luo Z, Song YF, Pan YX, Wu K, You WJ. Identification of autophagy related genes LC3 and ATG4 from yellow catfish Pelteobagrus fulvidraco and their transcriptional responses to waterborne and dietborne zinc exposure. Chemosphere. 2018;175:228–38.
Suwansa Ard S, Kankuan W, Thongbuakaew T, Saetan J, Kornthong N, Kruangkum T, Khornchatri K, Cummins SF, Isidoro C, Sobhon P. Transcriptomic analysis of the autophagy machinery in crustaceans. BMC Genomics. 2016;17:1–14.
He H, Dang YJ, Dai FY, Guo ZK, Wu JX, She XY, Pei Y, Chen YJ, Ling WH, Wu CQ, et al. Post-translational modifications of three members of the human MAP1LC3 family and detection of a novel type of modification for MAP1LC3B. J Biol Chem. 2003;278(31):29278–87.
Wu JX, Dang YJ, Su W, Liu C, Ma HJ, Shan YX, Pei Y, Wan B, Guo JH, Yu L. Molecular cloning and characterization of rat LC3A and LC3B–two novel markers of autophagosome. Biochem Biophys Res Commun. 2006;339(1):437–42.
Wang JR. Beclin 1 bridges autophagy, apoptosis and differentiation. Autophagy. 2008;4(7):947–8.
Tang H, Da L, Mao Y, Li Y, Li D, Xu ZH, Li F, Wang YF, Tiollais P, Li T, et al. Hepatitis B virus X protein sensitizes cells to starvation-induced autophagy via up‐regulation of beclin 1 expression. Hepatology. 2009;49(1):60–71.
Kong HJ, Moon J, Nam B, Kim Y, Kim W, Lee J, Kim K, Kim B, Yeo S, Lee CH, et al. Molecular characterization of the autophagy-related gene Beclin-1 from the olive flounder (Paralichthys olivaceus). Fish Shellfish Immunol. 2011;31(2):189–95.
Liu W, Li FF, Sun HJ, Wang YF, Yu GH, Wang FL, Qian YM. Tobacco mosaic virus infection on tobacco plants induces autophagy. Acta Phytopathologica Sinica. 2016;46(6):42–9.
Kroemer G, Mariño G, Levine B. Autophagy and the integrated stress response. Mol Cell. 2010;40(2):280–93.
Mariño G, Santano MN, Baehrecke EH, Kroemer G. Self-consumption: the interplay of autophagy and apoptosis. Nat Rev Mol Cell Biol. 2014;15(2):81–94.
Dong WN, Liao MQ, Zhuang XQ, Huang L, Liu C, Wang FF, Yin XL, Liu Y, Liang QJ, Wang WN. MYC drives autophagy to adapt to stress in Penaeus vannamei. Fish Shellfish Immunol. 2022;126:187–96.
Wei W, He YY, Li ZX, Zhou YX, Li J. Cloning of ATG5 gene of Fenneropenaeus chinensis and expression analysis under pH and carbonate alkalinity stress. Progress Fish Sci. 2022;43(3):84–94.
He YY, Li ZX, Zhang HE, Hu S, Wang QY, Li J. Genome-wide identification of Chinese shrimp (Fenneropenaeus chinensis) microRNA responsive to low pH stress by deep sequencing. Cell Stress Caperones. 2019;24(4):689–95.
Yang ZG, Zhang XD, Tang SL, Li HY. Research progress of autophagy on regulation of energy metabolism in tumor cells. Chin J Bases Clin Gen Surg. 2015;22(12):1525–9.
Dustin ML, Cooper JA. The immunological synapse and the actin cytoskeleton: molecular hardware for T cell signaling. Nat Immunol. 2000;1(1):23–9.
Sun Y. Transcriptome analysis of genes and signaling pathways involved in immunity under thermal /hypoxic stresses in Haliotis diversicolor. Fujian, China: JiMei University; 2018.
Cui W, Ma A, Wang X. Response of the PI3K-AKT signalling pathway to low salinity and the effect of its inhibition mediated by wortmannin on ion channels in turbot Scophthalmus maximus. Aquac Res. 2020(4).
Armenteros JJA, Tsirigos KD, Sønderby CK, Petersen TN, Winther O, Brunak S, Heijne Gv, Nielsen H. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat Biotechnol. 2019;37(4):420–3.
Horton P, Keun Joon Park, Obayashi T, Naoya Fujita, Harada H, Collier CJA, Nakai K. Wolf Psort protein localization predictor. Nucleic Acids Res. 2007;35(Web Server Issue):585–7.
Huang J, Song XL, Zhang YJ. The components of an inorganic physiological buffer for Penaeus chinensis. Methods cell Science: Official J Soc Vitro Biology. 1999;21(4):225–30.
Michael Z. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003;31(13):3406–15.
Schmittgen TD, Livak KJ, Schmittgen TD. Livak KJAnalyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.
Funding
This work was supported by grants from the National Natural Science Foundation of China [Grant No. 42376111], National Key R & D Program of China [No. 2022YFD2400104-03], China Agriculture Research System of MOF and MARA [No. CARS-48], Central Public-interest Scientific Institution Basal Research Fund, CAFS [No. 2023TD50] and ‘First Class Fishery Discipline’ program in Shandong Province, China.
Author information
Authors and Affiliations
Contributions
GCH conducted analysis and wrote the draft of manuscript. GCH, LYL, WQ, WJJ and TCJ performed data analysis. LZX and HYY conceived and designed the experiments, reviewed, and edited the writing of the manuscript. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The Guidelines for the Care and Use of Laboratory Animals in China were followed when carrying out all the experiments. The Institutional Animal Care and Use Committee (IACUC) of the Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences (Qingdao, China), approved the study, the study is reported in accordance with ARRlVE guidelines.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Guan, C., Li, Y., Wang, Q. et al. Genome-wide identification of ATG genes and their expression profiles under biotic and abiotic stresses in Fenneropenaeus chinensis. BMC Genomics 25, 625 (2024). https://doi.org/10.1186/s12864-024-10529-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-024-10529-2