Morphological and Molecular Functional Evidence of the Pharyngeal Sac in the Digestive Tract of Silver Pomfret, Pampus argenteus

The pharyngeal sac is a comparatively rare organ in the digestive tract among teleost fishes. However, our understanding of this remarkable organ in the silver pomfret (Pampus argenteus) is limited. In the present study, we examined the various morphological and histological characteristics of the pharyngeal sac using histochemical techniques and electron microscopy. The pharyngeal sac showed unique characteristics such as well-developed muscular walls, weakly keratinized epithelium, numerous goblet cells, and needle-like processes on the papillae. The porous cavity of the papillae contained numerous adipocytes and was tightly enveloped by type I collagen fibers. These structures might provide mechanical protection and excellent biomechanical properties for grinding and shredding prey. A comparison of gene expression levels between the pharyngeal sac and esophagus using RNA-seq showed that phenotype-associated genes (epithelial genes and muscle genes) were upregulated, whereas genes related to nutrient digestion and absorption were downregulated in the pharyngeal sac. These results support the role of the pharyngeal sac in shredding and predigesting food. Overall, these findings provide a clearer understanding of the pharyngeal sac morphology and explain the morphological adaptations of the digestive tract for feeding on gelatinous prey. To our knowledge, this is the first report on pharyngeal sac gene expression in P. argenteus.


Introduction
The digestive systems of vertebrates have the classic function of transforming complex, non-absorbable food into simple, absorbable forms [1]. According to their distinct feeding habits, diets, and environments, fishes exhibit specific morphological characteristics of the digestive tract and specialized functions necessary to obtain food and survive [2][3][4][5]. Notably, most stromateoidei fishes (Teleostei: Perciformes) are well characterized by a pharyngeal sac in the digestive tract and have an unusual diet of medusae during part or all of their lives [6,7]. The pharyngeal sac is a saccular outgrowth located on the proximal portion of the esophagus, with internally muscular walls covered by numerous projections commonly referred to as papillae [6,[8][9][10]. Early research on the pharyngeal sac focused on its morphology and the morphological characteristics that can be used as taxonomic characters [11][12][13]. The pharyngeal sac was initially described as an extra-esophageal growth or pharyngeal origin [7,8]. However, recent research analyzing its skeletal and muscular components has suggested that this organ originated from both the gill arches and the esophagus [14]. In addition, previous studies have examined the morphology, histology, and digestive enzyme activity of the pharyngeal sac [15][16][17]. However, to our knowledge, no detailed studies about the gross anatomy and histology of the pharyngeal sac and papillae are available, and related molecular research on the pharyngeal sac has not been reported to date.

Gross Anatomy of the Pharyngeal Sac
It was observed that the digestive tract of P. argenteus consists of the oropharyngeal cavity, pharyngeal sac, esophagus, stomach, pyloric caeca, and long intestine, with a total length of 45.12 ± 0.27 cm ( Figure 1A). The pharyngeal sac was located at the vertebral column on the dorsal side of the body, opening and connecting to the last gill arch, followed by a short esophagus. The pharyngeal sac was divided into two bilateral compartments by a sagittal sulcus (muscle ridge), with a mean size 6.3 ± 0.27 cm and a mean length 2.3 ± 0.5 cm. The pharyngeal sac was brown-black, with a muscular external appearance ( Figure 1B,C). The inner muscular wall of the pharyngeal sac was coated with sticky mucous secretions and showed numerous, hard papillae with needle-like secondary processes. The longest papillae were located in the center of the sac; from the center outward, the number of papillae increased, but their height gradually reduced ( Figure 1D). The papilla had an elongated axis and a bottom with a 5-8 radial cone-shaped stellate base, which was anchored in the muscle layer to support the papilla ( Figure 1E).

Analysis of Key Morphological Characteristics of the Pharyngeal Sac
The pharyngeal sac consisted of four layers-mucosa, submucosa, muscularis, and serosa-from internal to external (Figure 2A). The serosa was composed of connective tissue, melanocytes, and blood vessels ( Figure 2B,C). Next, we focused on the distinctive morphological features of the pharyngeal sac.

Analysis of Key Morphological Characteristics of the Pharyngeal Sac
The pharyngeal sac consisted of four layers-mucosa, submucosa, muscularis, and serosa-from internal to external (Figure 2A). The serosa was composed of connective tissue, melanocytes, and blood vessels ( Figure 2B,C). Next, we focused on the distinctive morphological features of the pharyngeal sac. and a bottom with a 5-8 radial cone-shaped stellate base, which was anchored in the muscle layer to support the papilla ( Figure 1E).

Analysis of Key Morphological Characteristics of the Pharyngeal Sac
The pharyngeal sac consisted of four layers-mucosa, submucosa, muscularis, and serosa-from internal to external (Figure 2A). The serosa was composed of connective tissue, melanocytes, and blood vessels ( Figure 2B,C). Next, we focused on the distinctive morphological features of the pharyngeal sac.

Muscularis
Histologically, the muscularis of the pharyngeal sac was well developed and composed of striated muscle. There was no obvious boundary between the two muscle layers, namely, the outer longitudinal layer and the inner circular layer ( Figure 3A,B). TEM showed a small amount of intercellular matrix and loose connective tissue between muscle fiber bundles ( Figure 3C,D). The sarcomere length was about 3.5 µm, with the A band in the bright area and the I band and the Z band in the dark area ( Figure 3E). In addition, a muscular crest was observed in the middle of the pharyngeal sac ( Figure 4A). Histologically, the muscular crest mainly consisted of striated muscle, collagen fibers, and blood vessels. The short papillae and stellate bases were visible in the muscular crest. (Figure 4B). With Masson's trichrome staining, the muscle fiber bundles of the muscular crest were stained red, and the collagen fibers were stained blue ( Figure 4C,D). Histologically, the muscularis of the pharyngeal sac was well developed and composed of striated muscle. There was no obvious boundary between the two muscle layers, namely, the outer longitudinal layer and the inner circular layer ( Figure 3A,B). TEM showed a small amount of intercellular matrix and loose connective tissue between muscle fiber bundles ( Figure 3C,D). The sarcomere length was about 3.5 μm, with the A band in the bright area and the I band and the Z band in the dark area ( Figure 3E). In addition, a muscular crest was observed in the middle of the pharyngeal sac ( Figure 4A). Histologically, the muscular crest mainly consisted of striated muscle, collagen fibers, and blood vessels. The short papillae and stellate bases were visible in the muscular crest. (Figure 4B). With Masson's trichrome staining, the muscle fiber bundles of the muscular crest were stained red, and the collagen fibers were stained blue ( Figure 4C,D).

Papilla
The papillae were located in the mucosa, submucosa, and muscularis layers. Through the full view of the papilla, three well-defined regions were observed: the elongated axis, secondary processes, and stellate base. For the first region, the elongated axis of the papilla was composed of the epithelium layer, collagen fiber layer, and adipose tissue layer from outside to inside. The epithelial layer contained stratified squamous cells at the top of papillae, interspersed with goblet cells. The goblet cells showed mucus-filled cytoplasm and were mainly found in the basal region of the papilla ( Figure 5A). Under the SEM, oval or circular columnar epithelial cells appeared embossed and contained densely packed, stubby microridges on the surface of the papilla. In some regions, the epithelial cells were disjointed by the opening of mucous cells ( Figure 5B,C). TEM showed the presence of microvilli on the surface of epithelial cells, whereas goblet cells were observed to be filled with mucus granules, vesicles, and electron-dense bodies ( Figure 5D). In addition, the collagen fiber layer was composed of type I collagen fibers arranged in parallel to form a close-knit envelope, separating the epithelial layer from the adipose tissue layer ( Figure 5H-K). A longitudinal section of the papilla showed the collagen fiber layer and the internal reticular porous structure ( Figure 6C). Moreover, many fat cells, blood vessels, and loose connective tissues were observed in the internal structure of the papillae ( Figure 5A). Oil red O staining showed that lipid droplets were distributed abundantly in the porous cavity of the papilla ( Figure 6D,E). For the second region, a cross section of the papilla showed eight secondary processes symmetrically arranged around the papilla, the bottom was connected to the papilla, and the tip was solid ( Figure 6A,B). SEM showed that the tips of the secondary processes penetrated through the mucosal epithelium and were exposed to the pharyngeal sac cavity with their non-cellular and smooth surface ( Figure 5E,F). For the third region, the stellate bases were wrapped by epithelial cells and anchored in a staggered manner in the submucosa and muscularis layers, with spaces between them ( Figure 5A). The SEM results showed the internal structure of the stellate base ( Figure 6G,H). TEM and Masson's trichrome staining results showed that it was composed of mineralized chitin arranged in parallel and collagen fibers

Papilla
The papillae were located in the mucosa, submucosa, and muscularis layers. Through the full view of the papilla, three well-defined regions were observed: the elongated axis, secondary processes, and stellate base. For the first region, the elongated axis of the papilla was composed of the epithelium layer, collagen fiber layer, and adipose tissue layer from outside to inside. The epithelial layer contained stratified squamous cells at the top of papillae, interspersed with goblet cells. The goblet cells showed mucus-filled cytoplasm and were mainly found in the basal region of the papilla ( Figure 5A). Under the SEM, oval or circular columnar epithelial cells appeared embossed and contained densely packed, stubby microridges on the surface of the papilla. In some regions, the epithelial cells were disjointed by the opening of mucous cells ( Figure 5B,C). TEM showed the presence of microvilli on the surface of epithelial cells, whereas goblet cells were observed to be filled with mucus granules, vesicles, and electron-dense bodies ( Figure 5D). In addition, the collagen fiber layer was composed of type I collagen fibers arranged in parallel to form a close-knit envelope, separating the epithelial layer from the adipose tissue layer ( Figure 5H-K). A longitudinal section of the papilla showed the collagen fiber layer and the internal reticular porous structure ( Figure 6C). Moreover, many fat cells, blood vessels, and loose connective tissues were observed in the internal structure of the papillae ( Figure 5A). Oil red O staining showed that lipid droplets were distributed abundantly in the porous cavity of the papilla ( Figure 6D,E). For the second region, a cross section of the papilla showed eight secondary processes symmetrically arranged around the papilla, the bottom was connected to the papilla, and the tip was solid ( Figure 6A,B). SEM showed that the tips of the secondary processes penetrated through the mucosal epithelium and were exposed to the pharyngeal sac cavity with their non-cellular and smooth surface ( Figure 5E,F). For the third region, the stellate bases were wrapped by epithelial cells and anchored in a staggered manner in the submucosa and muscularis layers, with spaces between them ( Figure 5A). The SEM results showed the internal structure of the stellate base ( Figure 6G,H). TEM and Masson's trichrome staining results showed that it was composed of mineralized chitin arranged in parallel and collagen fibers ( Figure 6I). On the basis of these findings, an illustration was created to depict the structure of the papillae (Figure 7). ( Figure 6I). On the basis of these findings, an illustration was created to depict the structure of the papillae (Figure 7).

De Novo Assembly and Functional Annotation
A total of six libraries were constructed from the pharyngeal sac and the esophagus. The clean reads of each sample ranged from 19.89 million to 23.78 million, Q30 was ≥92%, and GC content was ≥46% (Table 1). The mapping rate for the large yellow croaker (Larimichthys crocea) genome was approximately 71-78% (Table S1). An optimized assembly of the reads generated 41,683 unigenes and 105,485 transcripts. The unigenes' length distribution of minimum length, mean length, and N50 length were 301 bp, 1898 bp, and 2910 bp, respectively (Table S2). In addition, the BUSCO analysis results yielded 978 BUSCOs, including 926 (94.6%) complete BUSCOs, 31 (3.16%) fragmented BUSCOs, and 21 (2.14%) missing BUSCOs ( Figure S1A). The Pearson correlation analysis results among samples showed that R 2 values ranged from 0.68 to 0.75, indicating the difference between the two groups was small ( Figure S1B). Raw read data were uploaded to the NCBI Sequence Read Archive (SRA) under the accession number PRJNA889559.
Among the 41,683 unigenes, 33,654 (80.73%) were annotated in the Nt database, and only 10,032 (24.06%) were annotated in the KOG database (Table 1). These transcripts were also tested against the NR database, and the species distribution showed that Lates calcarifer accounted for 19.3% of all homologous genes, followed by Seriola dumeri (12.2%) ( Figure S1C). The Venn diagrams regarding annotation in five databases displayed that 8485 unigenes were annotated in all databases, while 9361 unigenes were uniquely anno-tated in the Nt database ( Figure S1D). All these results indicated that the sequences were appropriate for further analysis.
All unigenes were classified into three major functional categories. There were 42 significant GO categories, including biological process (25 subclasses), cellular component (5 subclasses), and molecular function (12 subclasses). The GO category of biological processes showed a high percentage of biological processes. In terms of biological processes, "cellular process", "metabolic process, "regulation of biological process", "biological regulation", and "response to stimulus" were highly enriched. The cellular component category contained "cellular anatomical entity" and "intracellular and protein containing complex" subclasses. With regard to molecular function, the top subclasses were "binding" and "catalytic activity" (Figure S2). In addition, in the functional classification of the KEGG pathways, there were 291 different pathways classified into five main categories, involving metabolism, environmental information processing, cellular processes, genetic information processing, and organismal systems. Among them, metabolism contained 127 sublevel pathways, including "carbohydrate metabolism", "lipid metabolism", "amino acid metabolism", etc. Other representative pathways enriched many unigenes such as "signal transduction" (2509 unigenes), "immune system" (1220 unigenes), and "endocrine system" (1122 unigenes) ( Figure S3). In addition, 10,032 (24 %) unigenes were annotated in KOG and grouped into 26 classifications. The largest cluster was "signal transduction mechanisms" (1842 unigenes), followed by "general function prediction only" (1586 unigenes) and "posttranslational modification", "protein turnover, chaperones" (999 unigenes) ( Figure S4). The top transcription factors in the pharyngeal sac were compiled and listed in Table S3.

DEGs Potentially Correlated with Structure and Function of the Pharyngeal Sac
The DEGs between pharyngeal sac and esophagus were analyzed to identify tissuespecific gene expression patterns. The heatmap results showed that there was a large number of genes with different expression patterns between the pharyngeal sac and the esophagus ( Figure 8A). The padj value < 0.05 and |log2Fold Change| > 1 were used as the threshold for detecting significant DEGs. A total of 8538 unigenes were identified as DEGs, including 4311 upregulated unigenes and 4227 downregulated unigenes ( Figure 8B). DEGs were grouped into four categories associated with enzyme genes (e.g., PGA, AMSase, SLC6A18, PRSS2, and CRTB), muscle-related genes (e.g., MYH, MYLPF, COL1A, and COL10A1), immune-related genes (e.g., MRC, TOLLIP, and IGHM), and metamorphosis-related genes (e.g., WNT5, TBX4, and PAX9) (Table S4).

Validation of Gene Expression Levels
To validate the expression levels of the DEGs obtained from the comparative RNA-seq analysis, 26 DEGs from various categories were selected for confirmation. Figure 11 showed that that the RT-qPCR data were basically consistent with the RNA-seq analysis results, supporting the reliability and accuracy of the RNA-seq analysis ( Figure 11).

Validation of Gene Expression Levels
To validate the expression levels of the DEGs obtained from the comp analysis, 26 DEGs from various categories were selected for confirmation. F that that the RT-qPCR data were basically consistent with the RNA-seq supporting the reliability and accuracy of the RNA-seq analysis ( Figure 11

Morphology of the Pharyngeal Sac
In the present study, we provided a clear description of the morphol Figure 11. Validation of RNA-seq data by RT-qPCR.

Morphology of the Pharyngeal Sac
In the present study, we provided a clear description of the morphological and histological features of the pharyngeal sac in P. argenteus. We found that the pharyngeal sac is composed of four tissue layers-the mucosa, submucosa, muscularis, and serosa-which is consistent with the previous descriptions of the digestive tract in vertebrates [25,31,32]. This study examined certain prominent features of the pharyngeal sac. First, the muscularis of the pharyngeal sac was particularly well developed, implying that food could be controlled to pass through the sac. Another prominent feature is the papillae lining the lateral walls; these papillae structures have been referred to as esophageal or pharyngeal sac "teeth" in previous studies [7,11,13], and were considered to be composed of homogeneous calcified tissue, pulp tissue, or bone tissue [11,33]. Contrary to these previous findings, we found that the papillae were composed of an epithelium layer, collagen fiber layer, and adipose tissue layer. The epithelium layer may protect the papilla from mechanical and chemical abrasions, as the surface of the papilla was most extensively exposed to friction during feeding. Moreover, the papillae had an internal honeycomb structure and secondary processes, which may make them lightweight and facilitate their ability to pierce food. Another key feature of the papillae was the stellate base anchored in the different muscle layers to improve space utilization. Although direct observation of the feeding behavior and mechanics was not feasible in the present study, the presence of features such as well-developed muscular walls; weakly keratinized epithelium cells; numerous goblet cells; and needle-like processes on the papillae confirmed the biomechanical role of papillae in the grinding and shredding of prey. Moreover, these findings suggest that the presence of the pharyngeal sac compensates for the small mouth and underdeveloped teeth in P. argenteus [34]. Similar papillae with finger-like or polypoid morphology have been reported in the masked angelfish (Genicanthus personatus). However, histological findings suggested that they contained no blood vessels or capillaries to absorb nutrients, which is inconsistent with our observations in P. argenteus [35]. Likewise, studies on the yellowtail fusilier (Caesio cuning) reported that forked papillae were located in the center of the esophagus and pointed backward to keep food moving down the esophagus [36]. Moreover, keratinized oropharynx papillae have been reported in turtles feeding on jellyfish; these papillae aid in shredding and provide mechanical protection from dangerous prey such as jellyfish [37]. The existence of these distinctive structure in different species, including P. argenteus, signifies adaptations related to food habits, dietary specialization, and effective manipulation of food. An important direction for future research on this specialized digestive tract is the underlying evolutionary mechanism.
The presence of papillae has a physiological significance likely associated with prey selection. In the natural environment, P. argenteus feed on gelatinous organisms such as medusae, ctenophores, and siphonophores [38], indicating that this fish possesses some immunity or protection mechanism against toxins. Previous studies have reported that the pharyngeal sac has a protective mucus-secreting system against the toxins of medusae [9]. In the present study, a novel finding was that the porous cavity of the papillae was filled with abundant adipocytes and enveloped by tightly arranged collagen fibers. However, few studies have been reported about the abundant adipocytes found in the digestive tract. It has been reported that large adipocytes typically show high metabolic activity and secrete higher quantities of immune cell chemoattractants [39]. Most stromateoidei fish have been reported to have an intradermal canal plexus situated under the epidermis, which is distributed all over the body and filled with a viscid, oily substance [9]; a recent report indicates that this system potentially protects the fish against jellyfish stings [14]. Inspired by this report, we speculate that the large amount of lipids accumulated in the porous cavity of the pharyngeal sac papillae is associated with feeding on and detoxification of jellyfish. Further research is required to identify the role of adipocytes in the papillae.

Candidate Genes Determined from Correlation between Structure and Function
To date, comparative transcriptome analyses have been performed in teleost species under normal physiological conditions. For example, RNA-seq analysis of three major organs (stomach, intestine, and rectum) was performed in S. quinqueradiata to understand digestion and nutrient absorption at the molecular level [29]. Moreover, the functional and transcriptional compartmentalization of the digestive tract in zebrafish (Danio rerio) has also been studied [40]. In the present study, we compared the gene expression levels between the pharyngeal sac and esophagus in P. argenteus by RNA-seq to identify tissue-specific gene expression patterns. The following mainly discusses the genes determined on the basis of correlation between structure and function.
Based on the RNA-seq data, we identified the candidate genes associated with complex phenotypes, such as epithelial genes and muscle growth genes. As keratin is a key part of the epithelial cytoskeleton, the keratin gene family (KRT) is important for the mechanical stability and integrity of epithelial cells and tissues. Fibrous proteins of this type form the structural framework of cells and constitute the tough skin, hair, and nails [41]. In the present study, the differentially expressed genes KRT8 and KRT13 were significantly enriched in the pharyngeal sac. Previous research has shown that KRT8 is expressed in the epithelial cells of the ectoderm and periderm in human fetal skin [41,42]. Moreover, KRT13 is expressed in the oral mucosa, esophagus, and forestomach, as the characteristic acidic keratin is produced in the superbase cells of non-cornified stratified epithelia [43]. These results suggested KRT8 and KRT13 play a role in forming the stratified epithelia of the pharyngeal sac. Similarly, muscle growth-related genes (e.g., MYL, MYH, MYLK, and COL1A) were significantly enriched in the pharyngeal sac. MYH and MYL have been reported as typical regulators of muscle contraction and are expressed in striated muscles [44,45]. To some degree, the high expression levels of muscle growth-related genes may play an important role in the muscle growth and motility of the pharyngeal sac.
As expected, our result showed that genes encoding proteolytic digestive enzymes (e.g., PGA, PRSS2, and AMCase) and lipid metabolism-related genes (e.g., APOE, LRP1, and FASN) were enriched significantly. Previous studies reported that PGA is involved in breaking down long-chain peptides during the earliest stage of protein digestion in the stomach of S. quinqueradiata [29]. PRSS plays a key role in hydrolyzing proteins and activating other digestive zymogenes [46]. AMCase exhibited chitinase activity in the esophagus, which was associated with the digestion of chitinous substances from crustaceans [47]. However, our results showed that the above enzyme-encoding genes were downregulated in the pharyngeal sac, indicating that the esophagus has a stronger digestive function than the pharyngeal sac. This finding is consistent with our previous finding that the pharyngeal sac exhibits the lowest digestive enzyme activity [17]. Moreover, the functions of the digestive tract include not only digestion and nutrient absorption but also immune defense against pathogens [48]. The intracellular TIR domain of Toll-like receptors (TLRs) interacts with adaptor proteins and triggers a signaling cascade that ultimately induces an immune response in the body [48]. TOLLIP encodes a key negative regulator of TLR-mediated innate immune responses [49]. These immune genes (e.g., MRC, TOLLIP, TLR3, and IGHM) were highly expressed in the pharyngeal sac. Moreover, the KEGG pathway related to human diseases was also significantly enriched. Although the role of these genes in the pharyngeal sac remains unclear, the present findings indicate that the pharyngeal sac, as the first segment of the digestive tract, has its own defense mechanism to protect against pathogens and toxins. Further research is needed to determine whether the pharyngeal sac has an immunological role in the digestive tract.
In conclusion, our findings provide a clearer understanding of the morphological and histological features of the pharyngeal sac in P. argenteus, confirming that its structure has an important biomechanical function in the grinding and shredding of prey. Moreover, to the best of our knowledge, this is the first report on gene expression in the pharyngeal sac. Candidate genes and KEGG pathways were associated with phenotype (epithelial and muscle genes), digestion, nutrient absorption, and immune function. Overall, the study of the pharyngeal sac is still a developing field of research with many knowledge gaps. Our results are expected to serve as a reference for future research aimed at studying the pharyngeal sac.

Experimental Fish and Sampling
Thirty adult individuals of P. argenteus (mean body weight 35.53 ± 0.32 g, mean body length 14.65 ± 2.16 cm) were obtained from Xiangshan Harbor Aquaculture and Larva Company (Ningbo, China). Healthy fish were cultivated in breeding ponds (temperature = 22.0−24.0 • C, dissolved oxygen = 9.85 mg/L, salinity = 33%, pH = 8.56). The fish were fed three times a day with commercial feed (Larvae Love, #6, Tokyo, Japan) equivalent to 2-3% of their body weight and were starved for 8 h before the experiment. Then, fish were euthanized using 75 mg/L buffered MS-222 (Finquel MS-222, Sigma Inc., Marlborough, MA, USA), and the pharyngeal sac and the esophagus were dissected, incised longitudinally, and spread out (n = 3). The samples were washed with 0.1 M phosphatebuffered saline and divided into four parts (Solarbio, Beijing, China): one part was fixed in 4% paraformaldehyde (Beyotime, Shanghai, China) to be used for paraffin sections; one part was embedded in optimal cutting temperature (OCT) compound to be used for frozen sections; a third part was placed in 2.5% glutaraldehyde solution (Innochem, Beijing, China) for electron microscopy observations; and for the fourth part, samples from three fishes were pooled as one sample and then rapidly preserved at −80 • C for RNA-seq and quantitative real-time PCR (RT-qPCR).

Light Miroscopy
The gross anatomy of the digestive tract was examined using a stereomicroscope (Leica, Germany), and photographs were taken in situ. Histological experiments were performed as previously described [50]. In brief, the samples were fixed in 4% paraformaldehyde for more than 48 h, dehydrated in a graded series of ethanol, and then cleared in xylene. Then, they were embedded in paraffin and sectioned into 5 µm thick sections using a table microtome (Leica, Wetzlar, Germany); these sections were stained with hematoxylin and eosin (HE), Masson's trichrome, Alizarin red, and Sirius red. Samples embedded in OCT were cut into 10 µm thick sections, which were stained with oil red O and stored at 4 • C. All sections were observed and imaged using a positive fluorescence microscope (Nikon, Tokyo, Japan), except for sections with Sirius red staining, which were observed using polarized light microscopy (Nikon, Tokyo, Japan).

Electron Microscopy
For SEM, the fixed pharyngeal sac and papillae were washed in 0.1 M sodium phosphate buffer (Real-Times, Beijing, China); dehydrated in a graded series of tert-butanol (Aladdin, Shanghai, China); and then dried using a critical point dryer (Nikon, Tokyo, Japan) for 24 h, mounted on carbon tape, and coated with gold palladium (approximately 30 nm thick). The samples were observed under a Hitachi S-3400N (Hitachi, Tokyo, Japan) operating at 25 kV.
For TEM, small fragments of the pharyngeal sac and papillae were washed in 0.1 M sodium phosphate buffer and post-fixed in 1% osmium tetroxide for 2 h at 25 • C. The samples were rinsed thoroughly with the same buffer (Real-Times, Beijing, China), dehydrated in an ascending series of acetone (Aladdin, Shanghai, China), and embedded in araldite resin. Then, semithin (60-90 nm thick) sections were cut using an ultramicrotome (Leica EM UC7, Wetzlar, Germany) and mounted on copper grids. Finally, the sections were stained with uranyl acetate and lead citrate, and images were captured under a Hitachi S-7650 (Hitachi, Tokyo, Japan) operating at 200 kV.

RNA Isolation, Library Preparation and Sequencing
Total RNA was extracted in 1.5 mL centrifuge tubes (Jet Biofil, Guangzhou, China) by using FastPure ® Cell/Tissue Total RNA Isolation Kit V2 (Vazyme, Nanjing, China) following the manufacturer's instructions. RNA integrity was assessed on 1% agarose gels and an RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA). The mRNA was purified from the total RNA by using poly-T oligoattached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in First-Strand Synthesis Reaction Buffer (5X). First-strand cDNA was synthesized using random hexamer primer and M-MLV Reverse Transcriptase, and then RNase was used to degrade the RNA. Second-strand cDNA synthesis was subsequently performed using DNA polymerase I and dNTP. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3 ends of DNA fragments, adaptors with a hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 370~420 bp in length, the library fragments were purified with an AMPure XP system (Beckman Coulter, Beverly, MA, USA). Finally, the Illumina NovaSeq 6000 platform (Beijing, China) was used for sequencing, and the end reads of 150 bp pairing were generated.

De Novo Assembly, Unigene Annotation and Functional Classification
Clean reads were filtered from raw data by removing adapter and low-quality reads. BUSCO analysis was used to evaluate the splicing quality. Then, Pearson correlation analysis was conducted to visualize the relative differences in the six samples. Gene function was annotated based on seven databases: the NCBI non-redundant protein sequences (Nr), NCBI non-redundant nucleotide sequences (Nt), Protein family (Pfam), Clusters of Orthologous Groups of proteins (COG/KOG), a manually annotated and reviewed protein sequence database (Swiss-Prot), KEGG Ortholog database (KO), and Gene Ontology (GO). The GOseq (1.10.0) and KOBAS (v2.0.12) software were used for GO function enrichment analysis and KEGG pathway enrichment analysis. The unigene expression was calculated and normalized with the reads per kb per million reads (RPKM) method (Mortazavi, Williams, Mccue, Schaeffer, & Wold, 2008). The genes with |log 2 Fold Change| ≥ 1 and padj ≤ 0.05 were identified as differentially expressed genes by using the DESeq2 R package (1.20.0). The GO function and KEGG pathway enrichment analysis of differentially expressed genes (DEGs) with padj ≤ 0.05 were defined as enriched significantly.

Validation of Transcriptome Data by RT-qPCR
To verify the reliability of the transcriptome results, 26 DEGs were selected for validation by RT-qPCR. The Primer 5.0 software was used to design the primers that are shown in Table 2. The system for RT-qPCR contains 5 µL SYBR Green I Master Mix (TransGen Biotech, Beijing, China), 0.2 µL forward and reverse primers, respectively, 1 µL cDNA template, and dd H 2 O added to 10 µL. Each DEG was repeated four times. RT-qPCR was performed on a q225 fluorescence quantitative qPCR machine (Kubo Technology, Beijing, China) with one cycle at 94 • C for 30 s, followed by 40 cycles at 94 • C for 5 s, 59 • C for 15 s, and 72 • C for 10 s. The relative expression levels were normalized to the endogenous control gene β-actin [51], and expression ratios were calculated by using the 2 −∆∆Ct method [52].