Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Physiol., 01 August 2025

Sec. Aquatic Physiology

Volume 16 - 2025 | https://doi.org/10.3389/fphys.2025.1642696

Epizootic shell disease induces systemic transcriptomic shifts in Homarus americanus, characterized by increased shell degradation and impaired energy metabolism across tissues

  • 1Department of Microbiology, College of Bio-convergence, Dankook University, Cheonan, Republic of Korea
  • 2Lewyt College of Veterinary Medicine, Long Island University, Brookville, NY, United States
  • 3Division of Marine Resources, New York State Department of Environmental Conservation, East Setauket, NY, United States
  • 4Smart Animal Bio Institute, Dankook University, Cheonan, Republic of Korea
  • 5Center for Bio Medical Engineering Core Facility, Dankook University, Cheonan, Republic of Korea
  • 6Department of Cosmedical and Materials, Dankook University, Cheonan, Republic of Korea

Epizootic shell disease (ESD) is characterized by shell erosion, pitting, and melanization in the American lobster (Homarus americanus) and is associated with a polymicrobial infection. The disease is multifactorial, with several contributing factors such as rising water temperatures and environmental pollution, which may facilitate bacterial invasion and increase host susceptibility. In a previous study, we found that the microbiome composition of the carapace in lobsters with ESD differed from that of healthy individuals, with ESD-associated bacteria enriched in the green gland and testis. However, the effects of bacterial infection on internal organs have not been clearly identified. In this study, we investigated the effects of ESD on four major tissues of the lobster (testis, intestine, hepatopancreas, and green gland) using transcriptomic analysis. A total of 564 genes were differentially expressed in the testis, 105 in the intestine, 333 in the hepatopancreas, and 112 in the green gland. The expression of the anti-lipopolysaccharide factor gene was increased in all tissues, indicating a systemic immune response to bacterial infection. Notably, chitinase genes involved in chitin degradation were upregulated, while the acetyl-coenzyme A transporter 1-like gene related to energy metabolism was significantly downregulated in the testis. In the intestine, expression of phosphoenolpyruvate carboxykinase cytosolic [GTP] and cytochrome P450 genes, which are involved in gluconeogenesis and xenobiotic metabolism, respectively, was reduced. The hepatopancreas showed decreased expression of hemocyanin genes, which play key roles in oxygen transport and immune defense in crustaceans. The green gland exhibited reduced expression of heat shock proteins involved in the cellular stress response, organic cation transporter proteins that mediate the excretion of organic cations, and UDP-xylose and UDP-N-acetylglucosamine transporters required for glycosylation and chitin biosynthesis. Together, these transcriptional changes suggest that ESD may compromise physiological functions such as immune defense, energy metabolism, and stress response, while promoting chitin degradation and cuticle remodeling in response to shell infection. This study revealed tissue-specific transcriptomic responses to ESD in the American lobster, providing a foundation for elucidating the molecular mechanisms underlying disease progression.

1 Introduction

The American lobster (Homarus americanus) is distributed throughout the coastal waters of the western North Atlantic. Long Island Sound, a large estuary between Connecticut and Long Island, was one of the largest lobster habitats in the United States until the late 1990s. Over the past 2 decades, however, the region has experienced a substantial decline in lobster populations (Long Island Sound Study, 2024). This decline has been associated with the increasing prevalence of epizootic shell disease (Castro et al., 2012). First reported in 1996, disease prevalence remains moderate to high, particularly in the eastern Long Island Sound (Castro and Somers, 2012; Tanaka et al., 2017).

Epizootic shell disease (ESD) is characterized by erosion, pitting, and melanization of the shell, resulting from bacterial infection (Smolowitz et al., 2005). Elevated water temperature is a key contributing factor, as high temperatures and other stressors can disrupt lobsters’ homeostasis, molting, and resistance to infection, increasing their susceptibility to ESD (Dove et al., 2005; Ishaq et al., 2022; Lorenzon et al., 2007). The genus Aquimarina has been consistently identified as the primary causative pathogen (Ooi et al., 2020; Quinn et al., 2017; Quinn et al., 2012), although various rod-shaped bacteria colonize shell lesions (Davies et al., 2014). In addition to shell infection, carapace bacteria have been detected in internal organs (Schaubeck et al., 2023), indicating the possibility of systemic infection. In lobsters with ESD, the expression of the ecdysteroid receptor and cytochrome P450 is increased in the hepatopancreas, while the expression of arginine kinase is downregulated in the muscle (Tarrant et al., 2010; Tarrant et al., 2012). Additionally, female lobsters with ESD show elevated levels of macroglobulin, a protease inhibitor related to the innate immune system, in the ovary (Tarrant et al., 2010). Male lobsters with ESD exhibit testicular abnormalities, resulting in over 50% of spermatozoa being nonviable (Taylor et al., 2021). These findings suggest that ESD is associated with systemic infection, which may impair immunity, metabolism, and reproductive function.

Transcriptome analysis provides a powerful approach for unraveling the complex biological alterations and pathogenic mechanisms underlying disease, enabling the identification of key molecular pathways. Under simulating ocean warming conditions, lobsters exposed to elevated temperatures after the larval stage showed up to a 7.1-fold increase in the expression of metabolism-related transcripts, suggesting that failure to meet rising energy demands could lead to increased mortality from disease and starvation (Harrington et al., 2020). In addition, reduced expression of chitin-related and pseudohemocyanin genes in the hepatopancreas has been linked to the progression of impoundment shell disease (Argenta et al., 2024). However, research on transcriptomic changes across internal tissues during ESD remains limited.

We have previously revealed that the carapace microbiota of lobsters with ESD harbors a high abundance of Aquimarina, Halocynthiibacter, and Tenacibaculum (Schaubeck et al., 2023). These carapace bacteria have also been detected in internal organs, including the green gland, hepatopancreas, intestine, and testis, and appear to be more common in lobsters with ESD compared to healthy lobsters. Therefore, it is hypothesized that ESD-associated bacteria translocate from the carapace to internal organs, where they disrupt local tissue homeostasis, leading to systemic dysfunction and disease progression. To test this hypothesis, this study aimed to identify tissue-specific molecular responses in lobsters affected by ESD.

2 Materials and methods

2.1 Lobster sampling

American lobsters were collected in August 2020 from eastern Long Island Sound and approximately 50 miles south of Montauk, in collaboration with local lobstermen and the New York State Department of Environmental Conservation. Carapace lengths ranged from 24.0 to 28.5 cm, and body weights ranged from 540.5 to 874.6 g (Supplementary Table S1). To prevent cross-contamination, lobsters were placed in a cooler on ice and sorted into two groups based on the presence (ESD, n = 37) or absence (HTH, n = 33) of shell lesions, then immediately transported to the laboratory. Following ice chilling for 30 min to minimize mobility and stress, the lobsters were dissected to collect various tissues, including the testis, intestine, hepatopancreas (HP), and green gland (GG). The collected tissues were gently washed with phosphate-buffered saline (PBS, pH 7.4) and stored in RNAlater stabilization solution (Qiagen) at −80°C until RNA extraction.

2.2 RNA sequencing

Total RNA was extracted from 10 mg of homogenized tissues using the RNeasy® Plus Mini Kit (Qiagen) following the manufacturer’s instructions. RNA quantity and purity were evaluated with a NanoDrop OneC microvolume UV-Vis spectrophotometer (Thermo Scientific). Purity was assessed based on the A260/A280 ratio, with values between 2.0 and 2.2, indicating acceptable RNA quality. The RNA samples were stored at −80°C until use in RNA sequencing (RNA-seq). The frozen RNA samples were shipped on dry ice to Genewiz (South Plainfield, NJ) for RNA-seq. The RNA-seq workflow comprised Poly A selection-based mRNA enrichment, mRNA fragmentation, and ramdom priming with subsequent first- and second-strand complementary DNA (cDNA) synthesis. End-repair 5′phosphorylation and adenine nucleotide (dA)-tailing were then performed. Finally, adaptor ligation, polymerase chain reaction (PCR) enrichment, and sequencing using the Illumina HiSeq 2,500 platform with 2 × 150 bp paired-end reads were performed. Quality controls and basic statistical methods were performed using FastQC (https://www.bioinformatics.babraham.ac.uk) to remove low quality reads as follows: Reads containing more than 10% skipped bases (labeled as ‘N'), sequencing reads containing more than 40% of bases with a quality score of less than 27, and average quality score (<27). Quality distributions of nucleotide, GC content, PCR duplicate characteristics, and k-mer sequencing data frequency were calculated.

2.3 Transcript read alignment

Sequencing reads were aligned to the GMGI_Hamer_2.0 H. americanus genome using STAR software (v.2.7.10a) with default parameters (accession: GCA_018991925.1, https://www.ncbi.nlm.nih.gov/data-hub/genome/GCF_018991925.1/). Then, gene expression was quantified using RSEM (v.1.3.3) software with default parameters and annotation from H. americanus Ensembl genes.

2.4 Analysis process for differentially expressed genes

Uniquely mapped read pairs were used for differentially expressed genes (DEGs) analysis. Gene expression counts were generated using RSEM (v.1.3.3) and normalized to transcripts per million (TPM). A total of 70 samples were initially used for RNA-seq. To ensure consistency, samples with a mapping rate below 65% and those identified as outliers were excluded. As a result, 53 samples were retained for downstream analysis. The final dataset consisted of the following tissues: testis (HTH = 6; ESD = 6), intestine (HTH = 6; ESD = 8), HP (HTH = 5; ESD = 5), and GG (HTH = 9; ESD = 8). All four tissue types used for RNA-seq were collected from the same individual lobsters (Supplementary Table S2). Gene expression differences between the two groups were considered statistically significant if |Log2 (fold change)| ≥ 1 and p-value <0.05 using the DESeq2 (v.1.34.0) in R package (v.4.4.1) (Love et al., 2014). A volcano plot was generated using the EnhancedVolcano package (Blighe et al., 2021), and a heatmap was generated using Pheatmap (v.1.0.12).

2.5 Gene Ontology function enrichment analysis

Gene Ontology (GO) term enrichment analysis was conducted using g:profiler (Raudvere et al., 2019). All genes with known statistical domain coverage were included, and multiple testing correction was applied using the Benjamini–Hochberg method with a significance threshold of 0.05. Gene annotations were obtained from Ensembl metazoan, which provides genomic data for non-vertebrate metazoan species including crustaceans (Howe et al., 2020). Pathway visualization was carried out using ggplot2 (v.3.5.2) in R.

2.6 Kyoto Encyclopedia of Genes and Genomes analysis

Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed to understand various biological functions and pathways in disease states. KEGG enrichment analysis of DEGs by tissue was performed using the clusterProfiler (v.4.12.6) in R (Kanehisa and Goto, 2000).

2.7 Selection of representative genes for expression analysis

To explore gene expression patterns, we visualized 16 DEGs (four from each organ) related to gluconeogenesis, innate immune system, and molting based on GO and KEGG pathway analysis. Gene expression values for each sample were converted to log2 (TPM +0.1) and visualized as violin plots using GraphPad Prism 10.

2.8 Validation of quantitative real-time polymerase chain reaction (qRT-PCR)

First-strand cDNA was synthesized from RNA samples using the Verso cDNA Synthesis Kit (Thermo Scientific) with a 3:1 blend of random hexamers and anchored oligo-dT. The synthesis conditions were as follows: 42°C for 30 min, 50°C for 30 min, 95°C for 2 min, followed by holding at 4°C. The synthesized cDNAs were stored at −80°C until further use. cDNA amplification was performed using gene-specific primers targeting acidic mammalian chitinase-like (LOC121854211), probable chitinase 2 (LOC121878393), chitinase-3-like protein 1 (LOC121860919), and acetyl-coenzyme A transporter 1-like (LOC121868265). To verify the functionality of the primers used in qRT-PCR, preliminary amplification was carried out using the T100 thermal cycler (Bio-Rad) at the Bio-medical engineering Core Facility of Dankook University. Following validation, qRT-PCR was performed using the CFX96 Duet Real-Time PCR System (Bio-Rad). All reactions were performed using 1 µL of cDNA, 10 pmol of primers and 2X Taq Pro Universal SYBR qRT-PCR Master Mix in a final volume of 20 µL. The PCR reaction started with a heating step at 95°C for 3 min, followed by 40 cycles of denaturation at 95°C for 10 s, annealing at the optimum temperature for 10 s, and extension at 72°C for 15 s. A melt curve analysis was performed from 62°C to 95°C, with a ramp rate of 0.5°C per minute. Gene expression was normalized to the expression of glycerol aldehyde-3-phosphate dehydrogenase (GAPDH), a housekeeping gene. A positive qRT-PCR reaction was detected based on the accumulation of fluorescent signals, and the cycle threshold (CT; the number of cycles required for the fluorescent signal to cross a fixed threshold) is inversely proportional to the amount of target nucleic acid in the sample. Relative gene expression was calculated using the 2−ΔΔCt value of each sample, and qRT-PCR was performed using five technical replicates per group for each cDNA sample. The primer sequences used are listed in Supplementary Table S1.

2.9 Statistical analysis

The pairwise Pearson correlation coefficient was computed using log-transformed TPM values for all expressed genes. The correlation coefficients were calculated using the base R function and visualized as a correlation matrix using the ggcorrplot (v.0.1.4.1) in R. Principal component analysis (PCA) was performed to assess sample distribution patterns using the built-in prcomp () function in R, and the results were visualized with ggplot2. Permutational multivariate analysis of variance (PERMANOVA) was conducted to statistically evaluate group differences by tissue type and disease status. TPM values were log2-transformed, followed by a Hellinger transformation, and then used to calculate a Euclidean distance matrix. PERMANOVA was performed with 10,000 permutations using the adonis function from the vegan (v.2.6–10) in R. TPM distributions across all samples were visualized using log10-transformed values, and density plots were generated with ggplot2 in R. The statistical significance of KEGG pathway enrichment was determined based on the Benjamini–Hochberg adjusted p-value (FDR <0.05). To validate RNA-seq results, qRT-PCR was performed, and relative gene expression levels (fold changes) were calculated using the 2−ΔΔCt method. Statistical comparisons between the ESD and HTH groups were conducted on ΔCt values using a Student’s t-test. Significance levels were set at p* <0.05 and p** <0.01.

3 Results

3.1 Clinical symptom and visual characteristics of epizootic shell disease

American lobsters were collected from eastern Long Island Sound in August 2020. Lobsters with ESD exhibited carapace lesions characterized by pitting, erosion, and melanization (Figure 1A). White triangles indicate pits, pink triangles indicate erosion of the carapace, and blue triangles denote both pits and melanin deposition.

Figure 1
A set of images analyzing lobster biology and gene expression: A) three images of lobster with color-coded indicators on different body sections.B) A pairwise correlation plot showing values with a color gradient from blue (0.00) to red (1.00), highlighting areas labeled tissue specificity.C) A heatmap representing Z-scores for various conditions, using a blue to red gradient.D) A PCA plot with labeled clusters representing different tissues and conditions.E) A density plot showing distribution of Log10(TPM) for various tissues and conditions, with color-coded legend.

Figure 1. Sample collection and RNA sequencing quality assessment. (A) American lobsters with epizootic shell disease (ESD). White triangles indicate pits; pink triangles indicate erosion; and blue triangles denote both pits and melanin deposition. (B) Pairwise correlation analysis showing tissue-specificity. (C) Heatmap of differentially expressed genes across tissues. Each column represents a gene, and each row indicates a sample. T: Testis; I: Intestine; HP: Hepatopancreas; GG: Green gland. (D) Principal component analysis showing tissue-specific transcriptional profiles (PERMANOVA, P < 0.0001). (E) Density plot showing consistent expression distributions across tissue samples. HTH: Healthy lobsters; ESD: Lobsters with ESD.

3.2 Transcriptome sequencing analysis

The Illumina HiSeq 2,500 instrument was used, and the read length was 150 bp. The number of uniquely mapped reads averaged 18,007,700 out of 24,950,208 total reads, with a mapping rate of 72.25% against the H. americanus reference genome (Supplementary Table S2). We performed pairwise correlation analysis based on the total gene expression values across samples to determine the reproducibility of technical replicates and differences in gene expression (Figure 1B). Consequently, robust correlations emerged across multiple tissues, suggesting tissue-dependent responses, with the testis exhibiting pronounced variability among lobsters. We visualized the DEGs from all organs using a heatmap and found that the gene expression differed by tissues and groups (Figure 1C). Notably, the testis showed clusters of genes that were strongly upregulated (red) or downregulated (blue) compared to other tissues, highlighting significant transcriptional changes compared to other tissues. PCA was performed to measure the distance between samples and confirmed significant separation by tissue (Figure 1D; PERMANOVA, p < 0.0001). Density plots showed overlapping curves across samples, indicating that the RNA-seq data were consistent in quality, reproducible across replicates, and comparable across tissue types (Figure 1E).

3.3 Differentially expressed genes analysis

We analyzed gene expression across tissues using DESeq2 to identify DEGs in the ESD group relative to the HTH group. In the testis, 513 genes were upregulated and 51 were downregulated in ESD compared to HTH. In the intestine, 51 genes were upregulated and 54 were downregulated. The HP exhibited 107 upregulated and 226 downregulated genes, while the GG showed 42 upregulated and 70 downregulated genes (Figure 2A). Among these DEGs, the anti-lipopolysaccharide factor gene (LOC121856986) was consistently upregulated across all tissues examined (Figure 2B). The testis, HP, and GG shared elevated expression of the soma ferritin-like (LOC121859498) and two anti-lipopolysaccharide factor-like genes (LOC121860186 and LOC121860188). Additionally, the testis, intestine, and HP shared increased expression of the lipase three gene (LOC121863816). The HP and GG demonstrated a consistent decrease in the expression of nine genes, including heat shock 70 kDa protein (LOC121855754, LOC121855757, and LOC121870242), la-related protein 6 (LOC121859532), cytochrome P450 2L1-like (LOC121876874), D-β hydroxybutyrate dehydrogenase (LOC121866996 and LOC121878121), estradiol 17-β-dehydrogenase 2 (LOC121878174), and an uncharacterized gene (LOC121880200).

Figure 2
Bar chart and Venn diagrams display differentially expressed genes (DEGs) across tissues. In the bar chart,

Figure 2. Differentially expressed genes (DEGs) across tissues in lobsters with epizootic shell disease (ESD). (A) Number of DEGs identified in each tissue. HP: Hepatopancreas; GG: Green gland. (B) Venn diagrams showing the number of shared upregulated and downregulated DEGs among tissues. The anti-lipopolysaccharide factor gene was consistently upregulated across all tissues in lobsters with ESD.

3.4 Gene Ontology of differentially expressed genes

All DEGs identified in the ESD group compared to the HTH group were analyzed for functional enrichment based on the GO terms. In the testis, the heatmap revealed distinct clustering of samples by health status (HTH vs ESD), although some variability among individuals was observed (Figure 3A). Upregulated pathways were associated with the GO terms chitin binding (GO:0008061), hydrolase activity, hydrolyzing O-glycosyl compounds (GO:0004553), proteolysis (GO:0006508), and oxidoreductase activity (GO:0016491). The chitin binding pathway was enriched with 13 DEGs, and six genes were predicted to encode chitinase (Supplementary Table S3). In contrast, the pathways for acetyl-CoA transmembrane transporter activity (GO:0008521) and DNA damage sensor activity (GO:0140612) were downregulated. In the intestine (Figure 3B), the most enriched upregulated pathways included serine hydrolase activity (GO:0017171) and ATPase-coupled transmembrane transporter activity (GO:0042626), while the most significantly downregulated pathways were phosphoenolpyruvate carboxykinase activity (GO:0004611) and monosaccharide biosynthetic process (GO:0046364). Additionally, genes involved in metabolism and redox balance were downregulated (Supplementary Table S4). In the HP (Figure 3C), the upregulated pathways included monooxygenase activity (GO:0004497), heme binding (GO:0020037), and chitin binding (GO:0008061), while the downregulated pathways were oxidoreductase activity (GO:0016491) and transmembrane transport (GO:0055085). Notably, ten downregulated genes associated with oxidoreductase activity (GO:0016491) were identified as hemocyanin A and B. The pathway related to ATP-dependent protein folding chaperone primarily contained six genes belonging to the heat shock protein 70 kDa (Supplementary Table S5). In the GG (Figure 3D), the spermine synthase activity (GO:0016768) pathway was upregulated, while transmembrane transport (GO:0055085) and ATP-dependent protein folding chaperone (GO:0140662) pathways were downregulated. Consistent with the pattern observed in the HP, heat shock protein showed decreased expression (Supplementary Table S6).

Figure 3
Clustered dot plots with heat maps depict gene set enrichment analysis for four panels (A, B, C, D), showing downregulation and upregulation of various gene functions. Dot size indicates the number of genes, and color denotes P-value significance. Heat maps display Z-scores, with blue to red indicating score range. Panel annotations include gene activities like oxidoreductase, iron ion binding, and transmembrane transport.

Figure 3. Gene ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in lobsters with epizootic shell disease (ESD) compared to healthy individuals. (A) Testis, (B) Intestine, (C) Hepatopancreas, (D) Green gland. Each dot represents a significantly enriched GO term. Dot size indicates the number of associated DEGs, and color represents statistical significance (adjusted p-value). HTH: Healthy lobsters; ESD: Lobsters with ESD. The heatmap illustrates the expression patterns of DEGs. Expression values were normalized using Z-score. Red and blue indicate upregulation and downregulation, respectively.

3.5 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis

KEGG pathway analysis was performed to identify metabolic alterations associated with ESD across lobster tissues. In the testis, four KEGG pathways were significantly enriched in the ESD group. The most significantly increased pathways (Figure 4A) included: (1) lysosome (hame04142), (2) phagosome (hame04145), (3) amino sugar and nucleotide sugar metabolism (hame00520), and (4) autophagy - animal (hame04140). These enrichments may indicate activation of cellular degradation pathways and immune defense, along with increased metabolic activity to meet elevated energy demands. In the HP, the steroid biosynthesis pathway (hame00100) was significantly increased, which may reflect increased synthesis of molting- and reproduction-related hormones. In the GG, the porphyrin metabolism pathway (hame00860) was significantly upregulated, suggesting a potential increase in demand for iron metabolism and detoxification in response to bacterial infection. In contrast, downregulated pathways were primarily observed in the HP and intestine (Figure 4B). In the HP, tyrosine metabolism (hame00350) and pentose and glucuronate interconversions (hame00040) were significantly decreased, suggesting dysregulation of amino acid metabolism and detoxification pathways. In the intestine, downregulated genes were more prominent, with seven KEGG pathways enriched: (1) citrate cycle (TCA cycle) (hame00020), (2) pyruvate metabolism (hame00620), (3) glycolysis/gluconeogenesis (hame00010), (4) FoxO signaling pathway (hame04068), (5) cysteine and methionine metabolism (hame00270), (6) glycine, serine and threonine metabolism (hame00260), and (7) lysine degradation (hame00310). These findings suggest that ESD may be associated with metabolic suppression in the intestine, potentially leading to reduced energy production and impaired tissue function. Detailed results of the KEGG enrichment analysis are provided in Supplementary Table S7.

Figure 4
Bar charts titled A and B display metabolic pathways with varying -Log10(P) values. In A,

Figure 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially expressed genes (DEGs) in lobsters with ESD. (A) Enriched KEGG pathways associated with upregulated genes. (B) Enriched KEGG pathways associated with downregulated genes. The y-axis represents pathway names, and the x-axis indicates statistical significance (p-values).

3.6 Expression analysis of selected target genes

We selected 16 genes involved in glucose biosynthesis (energy metabolism), immune function, and the chitin synthesis pathway to further examine differential expression between the ESD and HTH groups. In the testis of the ESD group, three chitinase-related genes—acidic mammalian chitinase-like (LOC121854211), chitinase-3-like protein 1 (LOC121860919), and probable chitinase 2 (LOC121878393) — were significantly increased, while acetyl-CoA transporter-1 (LOC121868265) was decreased (Figure 5A). In the intestine, all three genes from the phosphoenolpyruvate carboxykinase cytosolic [GTP] (LOC121879005, LOC121879006, and LOC121879007), as well as cytochrome P450 6k1-like (LOC121865918), were decreased in the ESD group (Figure 5B). In the HP, UDP-xylose and UDP-N-acetylglucosamine transporter-like (LOC121853811), heat shock cognate 70 kDa protein (LOC121862092), hemocyanin A chain-like (LOC 121877422), and hemocyanin B chain-like (LOC121863701) were decreased in the ESD group (Figure 5C). Similarly, in the GG, several genes showed decreased expression, including heat shock protein 70 kDa protein-like (LOC121855757), heat shock 70 kDa protein cognate 4-like (LOC121870242), UDP-xylose and UDP-N-acetylglucosamine transporter-like (LOC121869218), and organic cation transporter protein-like (LOC121863009) (Figure 5D). These results demonstrate tissue-specific molecular responses to ESD, characterized by enhanced chitinase activity in the testis and suppressed metabolic and stress-related pathways in the intestine, HP, and GG.

Figure 5
Four panels (A, B, C, D) display volcano plots and violin plots of gene expression data. Each panel highlights significant gene expressions with red and blue colors. Volcano plots show the log2 of fold change versus the negative log10 of the p-value, indicating differential expression. Violin plots illustrate expression levels for specific genes between HTH and ESD conditions. Panels A and C emphasize chitinase and transporter genes, while panels B and D focus on unidentified genes and heat shock proteins. Significant genes are marked by squares.

Figure 5. Distribution of differentially expressed genes (DEGs) and expression of key genes across tissues. (A) Testis, (B) Intestine, (C) Hepatopancreas, and (D) Green gland. Volcano plots display DEGs, with the y-axis representing -log10 (p-value) and the x-axis showing log2 (fold change). Upregulated genes are shown in red, downregulated genes in blue, and non-significant genes in gray. Violin plots illustrate the expression levels of selected key genes. The y-axis indicates log2 (TPM +1), and the x-axis indicates sample groups. The horizontal line within each violin plot represents the median expression level. *LOC121879005, *LOC121879006, and *LOC121879007: annotated as phosphoenolpyruvate carboxykinase, cytosolic [GTP]-like enzymes involved in gluconeogenesis. UDP-xylose transporter-like: annotated as UDP-xylose and UDP-N-acetylglucosamine transporter involved in glycosylation pathways.

3.7 Validation of quantitative reverse transcription polymerase chain reaction (qRT-PCR)

We performed qRT-PCR to confirm the expression levels of three chitinase genes and the acetyl-coenzyme A transporter 1-like gene in the testis of the ESD and HTH groups (Figure 6A). These genes were selected based on enrichment on gene expression and functional relevance to chitin metabolism and immune response pathways. qPCR validation was performed using the same testis samples as those used for RNA-seq. Expression of the chitinase genes, including acidic mammalian chitinase, chitinase-2, and chitinase-3, increased by 3.09-fold (p = 0.09), 2.66-fold (p = 0.28), and 2.99-fold (p = 0.11), respectively, in the ESD group compared to the HTH group. Although these changes were not statistically significant, the qRT-PCR results confirmed the gene expression patterns observed in the RNA-seq data (Figure 6B). In contrast, expression of the acetyl-coenzyme A transporter 1-like gene was significantly decreased by 0.59-fold (p < 0.01) in the ESD group relative to the HTH group. These findings suggest that ESD may involve dysregulation of genes related to chitin metabolism and energy transport.

Figure 6
Bar charts comparing gene expression in two conditions, HTH and ESD. A: Expression measured by 2^(-ΔΔCt) for four genes. Notable difference in Acetyl-CoA transporter-1, with lower expression in ESD, marked by double asterisks. B: Expression measured by TPM for the same genes. Significant increase in ESD for all genes except Acetyl-CoA transporter-1, where lower expression is marked by an asterisk.

Figure 6. Validation of expression patterns for selected differentially expressed genes (DEGs) in the testis. (A) Relative expression levels measured by qRT-PCR, presented as fold change using the 2−ΔΔCt method. Statistical analyses were performed on ΔCt values using an unpaired Student’s t-test. (B) Transcript abundance from RNA-seq, expressed as transcripts per million (TPM). Statistical significance was assessed using the Wald test. The genes analyzed include acidic mammalian chitinase-like (LOC121854211), probable chitinase 2 (LOC121878393), chitinase-3-like (LOC121860919), and acetyl-CoA transporter-1 (LOC121868265). Data represent mean values, and error bars indicate standard deviation (SD) based on biological replicates (n = 5). *p < 0.05; **p < 0.01.

4 Discussion

Our study is the first to characterize gene expression changes in the internal organs of the American lobster with ESD using high-throughput sequencing technology. Among 1,114 DEGs, the anti-lipopolysaccharide factor-like gene was upregulated in all tested organs of lobsters with ESD. Anti-lipopolysaccharide factors (ALF) such as ALFHa-1 and ALFHa-2 have been identified in H. americanus, and notably, the expression of the ALFHa-1 gene has been shown to increase in the gills following exposure to Vibrio fluvialis (Beale et al., 2008). ALF function as key components of the innate immune system by binding to lipopolysaccharides (LPS) in the outer membrane of Gram-negative bacteria, thereby neutralizing endotoxins and increasing permeability of bacterial membrane (Chaby, 2004). Therefore, the upregulation of ALF-like gene expression in the testis, intestine, HP, and GG indicates an innate immune response to Gram-negative bacterial invasion and supports the presence of systemic infection associated with ESD (Schaubeck et al., 2023).

We identified the chitin binding pathway as the most enriched differentially expressed gene pathway in the testis, with chitinase genes significantly upregulated within this pathway. Chitin is a homopolymeric carbohydrate composed of repeating N-acetyl-β-D-glucosamine (GlcNAc) units (Rowley and Coates, 2023). This compound is primarily found in the shells of crustaceans and the cell walls of fungi, providing the rigidity and mechanical strength needed to support the structure and shape of the organism (Zhang et al., 2021). Crustaceans such as shrimp and crabs contain the highest levels of chitin, which is highly associated with their growth and development (Zhang et al., 2021). Chitinase hydrolyzes chitin at intermediate regions to produce GlcNAc oligomers and is known to be induced in shrimp shells before molting (Rocha et al., 2012). Therefore, the observed upregulation of chitinase in lobsters with ESD may represent a molting-like response induced by shell infection. Although the upregulation of chitinase expression in the testis has not been widely reported, it is plausible that systemic immune and stress responses induced by ESD may also affect reproductive tissues. Chitinases are known to be involved not only in chitin degradation but also in immune defense and tissue remodeling across various tissues (Niu et al., 2018; Liu et al., 2019; Muthukrishnan et al., 2019; Lee et al., 2011). Therefore, the increased expression of chitinase observed in the testis may reflect a localized immune activation or structural remodeling of damaged tissue. Further investigation is required to clarify the functional role of chitinase in testicular physiology, particularly under disease conditions. In contrast, we observed significant downregulation of the acetyl-CoA transmembrane transport pathway in the testis. Acetyl-CoA is directly involved in major cellular processes, including energy metabolism, mitosis, and autophagy (Pietrocola et al., 2015), and also serves as a key substrate for chitin biosynthesis (Abo Elsoud and Kady, 2019). Additionally, we identified downregulation of genes encoding the transport protein responsible for the translocation of UDP-N-acetylglucosamine, a key precursor in the chitin synthesis pathway. Given that molting and reproduction are tightly coupled processes in crustaceans (Waddy et al., 2017), disruptions in chitin biosynthesis within the testis may reduce the systemic availability of chitin precursors, thereby impairing exoskeletal regeneration during infection.

We found that the expression of phosphoenolpyruvate carboxykinase cytosolic [GTP]-like gene was downregulated in the intestine of lobsters with ESD. The intestine, as the first line of defense, is one of the main sites of microbial and parasitic infections and serves as a critical organ that requires an efficient and specific immune response to eliminate invading pathogens (Du et al., 2016; Su et al., 2024). Phosphoenolpyruvate carboxykinase catalyzes the conversion of oxaloacetate to phosphoenolpyruvate within the gluconeogenic pathway, thereby playing a key role in maintaining glucose homeostasis (Jane et al., 2024; Sudhakaran et al., 2022; Xing et al., 2025; Yu et al., 2021). For example, increased activity of phosphoenolpyruvate carboxykinase enhances glucose production, which can exacerbate diabetes, whereas enzyme deficiency impairs gluconeogenesis, leading to hypoglycemia (Yu et al., 2021). In addition, phosphoenolpyruvate carboxykinase has been implicated in cold tolerance in lobsters (Farhadi et al., 2023). Thus, the downregulation of phosphoenolpyruvate carboxykinase in lobsters with ESD may indicate impaired glucose production and reduced capacity for stress adaptation.

The HP, also known as the digestive gland, is a major organ involved in energy metabolism, digestion, and nutrient storage, and is known to play an important role in various physiological processes such as stress adaptation and immune response (Farhadi et al., 2023; Ren et al., 2019; Roszer, 2014; Xin et al., 2024; Zhou et al., 2017). Therefore, the gene expression profile of the HP is considered an ideal indicator of stress, immune, and metabolic status (Xin et al., 2024). We found that the expression of the heat shock protein (HSP) 70 kDa genes was significantly downregulated in the HP of lobsters with ESD. HSPs are important for the maintenance of protein homeostasis in cells (Junprung et al., 2021) and play a role in cellular processes during and after exposure to oxidative stress caused by harmful environmental and microbial agents (Kumar et al., 2022). Thus, a decrease in these genes may reflect dysregulation of the cellular stress response. We also found that hemocyanin gene expression was significantly downregulated in the HP of lobsters with ESD. Hemocyanin is a copper-containing protein responsible for oxygen binding, transport, and storage in the blood of invertebrates, and performs several additional functions, including phenoloxidase-like activity (melanization, cuticle hardening), hormone transport, ecdysis, and blood clot formation (Coates and Costa-Paiva, 2020; Kuballa and Elizur, 2008). Hemocyanin also functions as a non-specific immune protein and is an essential component of innate immunity in aquatic invertebrates (Coates and Talbot, 2018; Ji et al., 2024). Together, the significant downregulation of hemocyanin gene expression in the HP may reflect impaired immune function and reduced shell integrity, possibly increasing susceptibility to bacterial infection.

The GG, also known as the antennal gland, is an excretory system responsible for maintaining water and mineral homeostasis, osmotic stabilization, ion regulation, and nitrogen excretion (Buranajitpirom et al., 2010; Khodabandeh et al., 2005; Kruangkum et al., 2024; Ross et al., 2019; Tsai and Lin, 2014). It has recently been shown to act as a natural entry point for viruses and bacteria and has also been implicated in the excretion of metabolites such as molting hormone (De Gryse, 2021; De Gryse et al., 2020). Supporting this role as an entry point, carapace bacteria, including Aquimarina, Halocynthiibacter, and Tenacibaculum, have been more frequently detected in the GG than in other organs, particularly in lobsters with ESD (Schaubeck et al., 2023). In the present study, we observed significant downregulation of genes encoding HSPs in the GG of lobsters with ESD, consistent with findings in the HP. In addition, genes encoding organic cation transporter proteins were also significantly downregulated in the GG. In humans, organic cation transporter one is primarily expressed in the intestines and liver, while transporter two is highly expressed in the kidneys (Motohashi et al., 2002; Nies et al., 2009; Zhang et al., 1997). Notably, the downregulation of these transporters has been associated with various cancers (Heise et al., 2012; Visentin et al., 2018). Therefore, the reduced expression of these transporter genes in the GG of lobsters suggests a potential disruption of excretory and detoxification functions, which may contribute to disease progression. Chitin, a major component of the exoskeleton, is composed of N-acetylglucosamine (Rowley and Coates, 2023; Zhang et al., 2021), which requires specific transporters for intracellular delivery. In this study, we identified downregulation of genes encoding UDP-xylose and UDP-N-acetylglucosamine transporters in the GG of lobsters with ESD. This downregulation may impair the transport of N-acetylglucosamine, thereby compromising chitin synthesis and exoskeletal integrity.

ESD primarily manifests as a localized infection of the shell, and lobsters may be able to overcome the disease by shedding the infected shell through molting. However, in our study, we observed a significant reduction in the expression of the estradiol 17β-dehydrogenase 2-like gene in both the HP and GG of lobsters with ESD. This gene is known to play a role in the metabolic regulation of molting in crustaceans. Previous studies have reported that 17β-estradiol is found in the HP and molting gland of crustaceans (Fu et al., 2022; Liu et al., 2021) and is involved in regulating molting-related processes (Fu et al., 2022; Head et al., 2019). In addition, RNA interference of estradiol 17β-dehydrogenase in the HP of Scylla paramamosain resulted in a significant decrease in shell hardness (Fu et al., 2022), further supporting the gene’s role in maintaining exoskeletal integrity. Together, these findings suggest that reduced expression of 17E2DH may impair molting, preventing the lobster from shedding the infected shell. This molting failure or delay could exacerbate the infection and potentially lead to systemic disease progression (Schaubeck et al., 2023). Therefore, the downregulation of the estradiol 17β-dehydrogenase 2-like gene serves as an important molecular indicator of impaired molting and disease severity in lobsters with ESD.

Our transcriptomic analysis revealed gene expression changes across multiple tissues in lobsters with ESD, highlighting enhanced chitin degradation, suppressed energy metabolism, and reduced cellular stress responses and immune defense. These results provide valuable insight into systemic host responses and the molecular mechanisms underlying disease progression. By uncovering ESD-related genes and tissue-specific expression patterns, this study enables the identification of potential diagnostic markers for the early detection of disease onset. These markers can also support the selection of disease-resistant stocks, contributing to more effective population management and improved long-term health of lobster populations. However, this study is limited by the incomplete gene annotation of H. americanus. Among the analyzed DEGs, some are uncharacterized genes with unknown functions, making it difficult to fully interpret their biological significance. Despite this limitation, our analysis identified key functional pathways involved in ESD and provides new evidence that ESD may be driven not only by localized shell infection but also by systemic dysfunction across internal organs. Continued efforts to improve genome annotation and functional characterization of genes in this species will be essential for bridging knowledge gaps in ESD pathogenesis and for developing effective strategies for disease management.

Data availability statement

BioSample metadata have been deposited in the NCBI BioSample database (http://www.ncbi.nlm.nih.gov/biosample/) under accession numbers SAMN48892171-SAMN48892223.

Ethics statement

Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because this study involved invertebrate species (lobsters), which are not subject to institutional animal care and use regulations. The animals were collected as part of routine fishery practices by local fishermen and staff of the Division of Marine Resources within the New York State Department of Environmental Conservation.

Author contributions

MK: Validation, Visualization, Data curation, Writing – review and editing, Writing – original draft. DC: Writing – review and editing. VC: Resources, Writing – review and editing. KH: Writing – review and editing. SM: Writing – review and editing. SJ: Writing – review and editing, Supervision, Writing – original draft, Resources, Conceptualization, Funding acquisition.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

The authors gratefully acknowledge the Center for Bio-Medical Engineering Core Facility at Dankook University for providing critical reagents and access to PCR equipment (NFEC-2019-09-257785).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2025.1642696/full#supplementary-material

References

Abo Elsoud M. M., Kady E. M. E. I. (2019). Current trends in fungal biosynthesis of chitin and chitosan. Bull. Natl. Res. Centre 43 (1), 1–12. doi:10.1186/s42269-019-0105-y

CrossRef Full Text | Google Scholar

Argenta N., Rampaul M., Clark K. F. (2024). American lobster (Homarus americanus) hepatopancreas transcriptome reveals the significance of chitin-related genes during impoundment shell disease. Facets 9, 1–10. doi:10.1139/facets-2024-0049

CrossRef Full Text | Google Scholar

Beale K. M., Towle D. W., Jayasundara N., Smith C. M., Shields J. D., Small H. J., et al. (2008). Anti-lipopolysaccharide factors in the American lobster homarus Americanus: molecular characterization and transcriptional response to Vibrio fluvialis challenge. Comp. Biochem. Physiol. Part D. Genomics Proteomics 3 (4), 263–269. doi:10.1016/j.cbd.2008.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Blighe K., Rana S., Lewis M. (2021). Publication-ready volcano plots with enhanced colouring and labeling. Available online at: https://github.com/kevinblighe/EnhancedVolcano.

Google Scholar

Buranajitpirom D., Asuvapongpatana S., Weerachatyanukul S., Wongprasert K., Namwong W., Poltana P., et al. (2010). Adaptation of the Black tiger shrimp, penaeus Monodon, to different salinities through an excretory function of the antennal gland. Cell Tissue Res. 340 (3), 481–489. doi:10.1007/s00441-010-0971-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Castro K. M., Cobb J. S., Gomez-Chiarri M., Tlusty M. (2012). Epizootic shell disease in American lobsters Homarus americanus in southern new England: past, present and future. Dis. Aquat. Organ. 100 (2), 149–158. doi:10.3354/dao02507

PubMed Abstract | CrossRef Full Text | Google Scholar

Castro K. M., Somers B. A. (2012). Observations of epizootic shell disease in American lobsters, Homarus americanus in southern new England. J. Shellfish Res. 31 (2), 423–430. doi:10.2983/035.031.0202

CrossRef Full Text | Google Scholar

Chaby R. (2004). Lipopolysaccharide-binding molecules: transporters, blockers and sensors. Cell Mol. Life Sci. 61 (14), 1697–1713. doi:10.1007/s00018-004-4020-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Coates C. J., Costa-Paiva E. M. (2020). Multifunctional roles of hemocyanins. Subcell. Biochem. 94, 233–250. doi:10.1007/978-3-030-41769-7_9

PubMed Abstract | CrossRef Full Text | Google Scholar

Coates C. J., Talbot J. (2018). Hemocyanin-derived phenoloxidase reaction products display anti-infective properties. Dev. Comp. Immunol. 86, 47–51. doi:10.1016/j.dci.2018.04.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Davies C. E., Whitten M. M. A., Kim A., Wootton E. C., Maffeis T. G. G., Tlusty M., et al. (2014). A comparison of the structure of American (Homarus americanus) and european (homarus Gammarus) lobster cuticle with particular reference to shell disease susceptibility. J. Invertebr. Pathol. 117, 33–41. doi:10.1016/j.jip.2014.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

De Gryse G. (2021). The discovery of the shrimp nephrocomplex (revision of the antennal gland) and its role during WSSV-entry and the moulting process. Merelbeke, Belgium: Ghent University, Faculty of Veterinary Medicine.

Google Scholar

De Gryse G. M. A., Khuong T. V., Descamps B., Van den Broeck W., Vanhove C., Cornillie P., et al. (2020). The shrimp nephrocomplex serves as a major portal of pathogen entry and is involved in the molting process. Proc. Natl. Acad. Sci. U. S. A. 117 (45), 28374–28383. doi:10.1073/pnas.2013518117

PubMed Abstract | CrossRef Full Text | Google Scholar

Dove A. D. M., Allam B., Powers J. J., Sokolowski M. S. (2005). A prolonged thermal stress experiment on the American lobster. J. Shellfish Res. 24 (3), 761–765. doi:10.2983/0730-8000(2005)24[761:APTSEO]2.0.CO;2

CrossRef Full Text | Google Scholar

Du Z., Jin Y., Ren D. (2016). In-depth comparative transcriptome analysis of intestines of red swamp crayfish, Procambarus clarkii, infected with WSSV. Sci. Rep. 6, 26780. doi:10.1038/srep26780

PubMed Abstract | CrossRef Full Text | Google Scholar

Farhadi A., Tang S., Huang M., Yu Q., Xu C., Li E. (2023). Identification of key immune and stress related genes and pathways by comparative analysis of the gene expression profile under multiple environmental stressors in Pacific white shrimp (litopenaeus Vannamei). Fish. Shellfish Immunol. 135, 108695. doi:10.1016/j.fsi.2023.108695

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu Y. Y., Liu X., Liu L., Fang W., Wang C. L. (2022). Identification and functional analysis of the estradiol 17β-dehydrogenase gene on the shell hardness of Scylla paramamosain during the molting cycle. Aquaculture 553, 738113. doi:10.1016/j.aquaculture.2022.738113

CrossRef Full Text | Google Scholar

Harrington A. M., Clark K. F., Hamlin H. J. (2020). Expected ocean warming conditions significantly alter the transcriptome of developing postlarval American lobsters (homarus Americanus): implications for energetic trade-offs. Comp. Biochem. Physiology D-Genomics Proteomics 36, 100716. doi:10.1016/j.cbd.2020.100716

PubMed Abstract | CrossRef Full Text | Google Scholar

Head T. B., Mykles D. L., Tomanek L. (2019). Proteomic analysis of the Crustacean molting gland (Y-organ) over the course of the molt cycle. Comp. Biochem. Physiol. Part D. Genomics Proteomics 29, 193–210. doi:10.1016/j.cbd.2018.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Heise M., Lautem A., Knapstein J., Schattenberg J. M., Hoppe-Lotichius M., Foltys D., et al. (2012). Downregulation of organic cation transporters OCT1 (SLC22A1) and OCT3 (SLC22A3) in human hepatocellular carcinoma and their prognostic significance. BMC Cancer 12, 109. doi:10.1186/1471-2407-12-109

PubMed Abstract | CrossRef Full Text | Google Scholar

Howe K. L., Contreras-Moreira B., De Silva N., Maslen G., Akanni W., Allen J., et al. (2020). Ensembl genomes 2020-enabling non-vertebrate genomic research. Nucleic Acids Res. 48 (D1), D689–D695. doi:10.1093/nar/gkz996

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishaq S. L., Turner S. M., Lee G., Tudor M. S., MacRae J. D., Hamlin H., et al. (2022). Warmer water temperature and epizootic shell disease reduces diversity but increases cultivability of bacteria on the shells of American lobster (Homarus americanus). bioRxiv. 2022–11. doi:10.1101/2022.11.11.512360

CrossRef Full Text | Google Scholar

Jane A., Rasher D. B., Waller J., Annis E., Frederich M. (2024). Rearing condition influences gene expression in postlarval American lobster (Homarus americanus). PLoS One 19 (7), e0307169. doi:10.1371/journal.pone.0307169

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji R., Guan L., Hu Z., Cheng Y., Cai M., Zhao G., et al. (2024). A comprehensive review on hemocyanin from marine products: structure, functions, its implications for the food industry and beyond. Int. J. Biol. Macromol. 269 (Pt 1), 132041. doi:10.1016/j.ijbiomac.2024.132041

PubMed Abstract | CrossRef Full Text | Google Scholar

Junprung W., Supungul P., Tassanakajon A. (2021). Structure, gene expression, and putative functions of crustacean heat shock proteins in innate immunity. Dev. Comp. Immunol. 115, 103875. doi:10.1016/j.dci.2020.103875

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa M., Goto S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28 (1), 27–30. doi:10.1093/nar/28.1.27

PubMed Abstract | CrossRef Full Text | Google Scholar

Khodabandeh S., Charmantier G., Charmantier-Daures M. (2005). Ultrastructural studies and Na+,K+-ATPase immunolocalization in the antennal urinary glands of the lobster homarus Gammarus (crustacea, Decapoda). J. Histochem Cytochem 53 (10), 1203–1214. doi:10.1369/jhc.4A6540.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Kruangkum T., Jaiboon K., Pakawanit P., Saetan J., Pudgerd A., Wannapaiboon S., et al. (2024). Anatomical and molecular insights into the antennal gland of the giant freshwater prawn Macrobrachium rosenbergii. Cell Tissue Res. 397 (2), 125–146. doi:10.1007/s00441-024-03898-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuballa A. V., Elizur A. (2008). Differential expression profiling of components associated with exoskeletal hardening in crustaceans. BMC genomics 9, 575–14. doi:10.1186/1471-2164-9-575

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar V., Roy S., Behera B. K., Das B. K. (2022). Heat shock proteins (hsps) in cellular homeostasis: a promising tool for health management in crustacean aquaculture. Life 12 (11), 1777. doi:10.3390/life12111777

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee C. G., Da Silva C. A., Dela C. C. S., Ahangari F., Ma B., Kang M. J., et al. (2011). Role of chitin and chitinase/chitinase-like proteins in inflammation, tissue remodeling, and injury. Annu. Rev. physiology 17 (1), 479–501. doi:10.1146/annurev-physiol-012110-142250

CrossRef Full Text | Google Scholar

Liu L., Fu Y., Xiao L., Liu X., Fang W., Wang C. (2021). iTRAQ-based quantitative proteomic analysis of the hepatopancreas in Scylla paramamosain during the molting cycle. Comp. Biochem. Physiol. Part D. Genomics Proteomics. 40, 100870. doi:10.1016/j.cbd.2021.100870

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu X., Zhang J., Zhu K. Y. (2019). “Chitin in arthropods: biosynthesis, modification, and metabolism,” in Targeting chitin-containing organisms, 169–207. doi:10.1007/978-981-13-7318-3_9

CrossRef Full Text | Google Scholar

Long Island Sound Study (2024). Lobster Abundance. Available online at: https://longislandsoundstudy.net/ecosystem-target-indicators/lobster-abundance/.

Google Scholar

Lorenzon S., Giulianini P. G., Martinis M., Ferrero E. A. (2007). Stress effect of different temperatures and air exposure during transport on physiological profiles in the American lobster Homarus americanus. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 147 (1), 94–102. doi:10.1016/j.cbpa.2006.11.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Love M. I., Huber W., Anders S. (2014). Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome Biol. 15 (12), 550. doi:10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Motohashi H., Sakurai Y., Saito H., Masuda S., Urakami Y., Goto M., et al. (2002). Gene expression levels and immunolocalization of organic ion transporters in the human kidney. J. Am. Soc. Nephrol. 13 (4), 866–874. doi:10.1681/ASN.V134866

PubMed Abstract | CrossRef Full Text | Google Scholar

Muthukrishnan S., Merzendorfer H., Arakane Y., Yang Q. (2019). “Chitin organizing and modifying enzymes and proteins involved in remodeling of the insect cuticle,” in Targeting chitin-containing organisms, 83–114. doi:10.1007/978-981-13-7318-3_5

CrossRef Full Text | Google Scholar

Nies A. T., Koepsell H., Winter S., Burk O., Klein K., Kerb R., et al. (2009). Expression of organic cation transporters OCT1 (SLC22A1) and OCT3 (SLC22A3) is affected by genetic factors and cholestasis in human liver. Hepatology 50 (4), 1227–1240. doi:10.1002/hep.23103

PubMed Abstract | CrossRef Full Text | Google Scholar

Niu S., Yang L., Zuo H., Zheng J., Weng S., He J., et al. (2018). A chitinase from Pacific white shrimp litopenaeus Vannamei involved in immune regulation. Dev. and Comp. Immunol. 85, 161–169. doi:10.1016/j.dci.2018.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ooi M. C., Goulden E. F., Trotter A. J., Smith G. G., Bridle A. R. (2020). Aquimarina sp. associated with a cuticular disease of cultured larval palinurid and scyllarid lobsters. Front. Microbiol. 11, 573588. doi:10.3389/fmicb.2020.573588

PubMed Abstract | CrossRef Full Text | Google Scholar

Pietrocola F., Galluzzi L., Bravo-San Pedro J. M., Madeo F., Kroemer G. (2015). Acetyl coenzyme A: a central metabolite and second messenger. Cell Metab. 21 (6), 805–821. doi:10.1016/j.cmet.2015.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinn R. A., Hazra S., Smolowitz R., Chistoserdov A. Y. (2017). Real-time PCR assay for Aquimarina macrocephali subsp. homaria and its distribution in shell disease lesions of Homarus americanus, milne-edwards, 1837, and environmental samples. J. Microbiol. Methods 139, 61–67. doi:10.1016/j.mimet.2017.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinn R. A., Metzler A., Smolowitz R. M., Tlusty M., Chistoserdov A. Y. (2012). Exposures of Homarus americanus shell to three bacteria isolated from naturally occurring epizootic shell disease lesions. J. Shellfish Res. 31 (2), 485–493. doi:10.2983/035.031.0208

CrossRef Full Text | Google Scholar

Raudvere U., Kolberg L., Kuzmin I., Arak T., Adler P., Peterson H., et al. (2019). g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 47 (W1), W191–W198. doi:10.1093/nar/gkz369

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren X., Liu P., Li J. (2019). Comparative transcriptomic analysis of marsupenaeus Japonicus hepatopancreas in response to Vibrio parahaemolyticus and white spot syndrome virus. Fish. Shellfish Immunol. 87, 755–764. doi:10.1016/j.fsi.2019.02.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Rocha J., Garcia-Carreño F. L., Muhlia-Almazán A., Peregrino-Uriarte A. B., Yépiz-Plascencia G., Córdova-Murueta J. H. (2012). Cuticular chitin synthase and chitinase mRNA of whiteleg shrimp litopenaeus Vannamei during the molting cycle. Aquaculture 330, 111–115. doi:10.1016/j.aquaculture.2011.12.024

CrossRef Full Text | Google Scholar

Ross E. P., Behringer D. C., Munoz A., Diaz D., Bojko J. (2019). A histological atlas for the palinuridae (crustacea: decapoda: achelata): a guide to parasite discovery and spotting the abnormal in spiny lobsters. J. Invertebr. Pathol. 163, 21–33. doi:10.1016/j.jip.2019.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Roszer T. (2014). The invertebrate midintestinal gland (“Hepatopancreas”) is an evolutionary forerunner in the integration of immunity and metabolism. Cell Tissue Res. 358 (3), 685–695. doi:10.1007/s00441-014-1985-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Rowley A. F., Coates C. J. (2023). Shell disease syndromes of decapod crustaceans. Environ. Microbiol. 25 (5), 931–947. doi:10.1111/1462-2920.16344

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaubeck A., Cao D., Cavaleri V., Mun S., Jeon S. J. (2023). Carapace microbiota in American lobsters (Homarus americanus) associated with epizootic shell disease and the green gland. Front. Microbiol. 14, 1093312. doi:10.3389/fmicb.2023.1093312

PubMed Abstract | CrossRef Full Text | Google Scholar

Smolowitz R., Chistoserdov A., Hsu A. (2005). “Epizootic shell disease in the American lobster,” in State of lobster science: Lobster shell disease workshop (Boston, MA: New England Aquarium), 2–11.

Google Scholar

Su X., Zhu X., Liang Z., Bao Z., Zhang J., Guo J., et al. (2024). Biochemical, histological and transcriptional response of intestines in litopenaeus Vannamei under chronic zinc exposure. Chemosphere 354, 141646. doi:10.1016/j.chemosphere.2024.141646

PubMed Abstract | CrossRef Full Text | Google Scholar

Sudhakaran G., Rajesh R., Guru A., Haridevamuthu B., Murugan R., Bhuvanesh N., et al. (2022). Deacetylated nimbin analog N2 fortifies alloxan-induced pancreatic beta-cell damage in insulin-resistant zebrafish larvae by upregulating phosphoenolpyruvate carboxykinase (PEPCK) and insulin levels. Toxicol. Appl. Pharmacol. 454, 116229. doi:10.1016/j.taap.2022.116229

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka K. R., Belknap S. L., Homola J. J., Chen Y. (2017). A statistical model for monitoring shell disease in inshore lobster fisheries: a case study in Long Island Sound. PLoS One 12 (2), e0172123. doi:10.1371/journal.pone.0172123

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarrant A. M., Franks D. G., Verslycke T. (2012). Gene expression in American lobster (Homarus americanus) with epizootic shell disease. J. Shellfish Res. 31 (2), 505–513. doi:10.2983/035.031.0210

CrossRef Full Text | Google Scholar

Tarrant A. M., Stegeman J. J., Verslycke T. (2010). Altered gene expression associated with epizootic shell disease in the American lobster, Homarus americanus. Fish. Shellfish Immunol. 29 (6), 1003–1009. doi:10.1016/j.fsi.2010.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Taylor L. J., Clark K. F., Daoud D., van den Heuvel M. R., Greenwood S. J. (2021). Exposure of American lobster (Homarus americanus) to the pesticide chlorpyrifos results in changes in gene expression. Comp. Biochem. Physiol. Part D. Genomics Proteomics 40, 100918. doi:10.1016/j.cbd.2021.100918

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsai J. R., Lin H. C. (2014). Functional anatomy and ion regulatory mechanisms of the antennal gland in a semi-terrestrial crab, Ocypode stimpsoni. Biol. Open 3 (6), 409–417. doi:10.1242/bio.20147336

PubMed Abstract | CrossRef Full Text | Google Scholar

Visentin M., Torozi A., Gai Z., Hausler S., Li C., Hiller C., et al. (2018). Fluorocholine transport mediated by the organic cation transporter 2 (OCT2, SLC22A2): implication for imaging of kidney tumors. Drug Metab. Dispos. 46 (8), 1129–1136. doi:10.1124/dmd.118.081091

PubMed Abstract | CrossRef Full Text | Google Scholar

Waddy S. L., Feindel N., Hamilton-Gibson N., Aiken D. E., Merritt V., Leavitt N. (2017). Reproductive cycles and mating capacity in Male American lobsters (Homarus americanus). Fish. Res. 1 (186), 358–366. doi:10.1016/j.fishres.2016.09.034

CrossRef Full Text | Google Scholar

Xin Z. Z., Zhang X. T., Zhou M., Chen J. Y., Zhu Z. Q., Zhang J. Y. (2024). Differential molecular responses of hemolymph and hepatopancreas of swimming crab, Portunus trituberculatus, infected with ameson portunus (microsporidia). Fish. Shellfish Immunol. 145, 109324. doi:10.1016/j.fsi.2023.109324

PubMed Abstract | CrossRef Full Text | Google Scholar

Xing Y. F., Feng Y. F., Tian J., Zhang Y., Han C., Liu X. J., et al. (2025). Nitrite and sulfide stress affects physiological and metabolic functions of gills in crab Eriocheir sinensis. Aquaculture 594, 741437. doi:10.1016/j.aquaculture.2024.741437

CrossRef Full Text | Google Scholar

Yu S., Meng S., Xiang M., Ma H. (2021). Phosphoenolpyruvate carboxykinase in cell metabolism: roles and mechanisms beyond gluconeogenesis. Mol. Metab. 53, 101257. doi:10.1016/j.molmet.2021.101257

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang L., Dresser M. J., Gray A. T., Yost S. C., Terashita S., Giacomini K. M. (1997). Cloning and functional expression of a human liver organic cation transporter. Mol. Pharmacol. 51 (6), 913–921. doi:10.1124/mol.51.6.913

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang X., Yuan J., Li F., Xiang J. (2021). Chitin synthesis and degradation in crustaceans: a genomic view and application. Mar. Drugs 19 (3), 153. doi:10.3390/md19030153

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou M., Abbas M. N., Kausar S., Jiang C. X., Dai L. S. (2017). Transcriptome profiling of red swamp crayfish (Procambarus clarkii) hepatopancreas in response to lipopolysaccharide (LPS) infection. Fish. Shellfish Immunol. 71, 423–433. doi:10.1016/j.fsi.2017.10.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: American lobster, epizootic shell disease, transcriptome analysis, antilipopolysaccharide factor, chitinase, acetyl-coenzyme A transporter, heat shock protein, polymicrobial infection

Citation: Kim M, Cao D, Cavaleri V, Han K, Mun S and Jeon SJ (2025) Epizootic shell disease induces systemic transcriptomic shifts in Homarus americanus, characterized by increased shell degradation and impaired energy metabolism across tissues. Front. Physiol. 16:1642696. doi: 10.3389/fphys.2025.1642696

Received: 09 June 2025; Accepted: 14 July 2025;
Published: 01 August 2025.

Edited by:

Yafei Duan, South China Sea Fisheries Research Institute, China

Reviewed by:

Ngoc Tuan Tran, Shantou University, China
Andressa Teles, Centro de Investigación Biológica del Noroeste (CIBNOR), Mexico

Copyright © 2025 Kim, Cao, Cavaleri, Han, Mun and Jeon. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Soo Jin Jeon, U29vamluLmplb25AbGl1LmVkdQ==; Seyoung Mun, bXVuc2V5b3VuZ0BnbWFpbC5jb20=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.