A Meta-Transcriptomics Survey Reveals Changes in the Microbiota of the Chinese Mitten Crab Eriocheir sinensis Infected with Hepatopancreatic Necrosis Disease

Infection of the freshwater Chinese mitten crab Eriocheir sinensis with hepatopancreatic necrosis disease (HPND) has been a major problem in the crab-cultivated Chinese Province of Jiangsu since 2015. To explore the etiology of HPND, meta-transcriptomic libraries of the hepatopancreata from crabs with and without HPND were constructed. Comparison analyses showed that there were no statistically significant differences in viral and microsporidial communities in the hepatopancreata of diseased and healthy crabs. Bacteroidetes, Proteobacteria, and Firmicutes were the most dominant bacterial phyla in the hepatopancreata of healthy crabs, with a combined prevalence of 93%. However, a decrease in bacterial diversity and a striking shift in the microbial composition were found in the hepatopancreata of crabs infected with HPND. Tenericutes was the most prevalent bacterial phylum in diseased crabs (31.82%), whereas its prevalence was low in healthy crabs (0.02%). By contrast, the prevalence of Bacteroidetes was significantly lower in crabs with HPND (3.49%) than in crabs without HPND (41.04%). We also found that the prevalence of Actinobacteria was higher in crabs with HPND (16.70%) than in crabs without the disease (4.03%). The major bacterial family within the Tenericutes phylum in crabs with HPND was detected by polymerase chain reaction and determined to be Mycoplasmataceae. In conclusion, there were striking changes in the microbiota of diseased and healthy crabs. Specifically, the prevalence of bacteria belonging to Tenericutes and Actinobacteria phyla increased, whereas the prevalence of bacteria belonging to the Bacteroidetes phylum decreased in crabs with HPND, clearly pointing to an association with HPND.

Introduction

The freshwater Chinese mitten crab E. sinensis is an economically important crustacean cultured in the Chinese provinces of Jiangsu, Anhui, Hubei, and Liaoning (Shen et al., 2014). With the rapid increase of the E. sinensis aquacultural industry, numerous diseases have recently evolved, thus resulting in huge economic losses (Shen et al., 2015). For example, hepatopancreatic necrosis disease (HPND, “shuibiezi” in Chinese), which affects E. sinensis, has been a major problem in the crab-cultivated area of Jiangsu Province since 2015. Crabs with HPND exhibit multiple clinical symptoms, including hepatopancreata that are lighter in color (turning from golden yellow and light yellow to gray–white) and soft shells that are darker in color than usual, as well as muscle atrophy and edema (Ding et al., 2016). Furthermore, the stomachs and intestines are empty in many diseased crabs. The majority of the diseased crabs do not die immediately; instead, they continue to molt, albeit at a later stage of breeding. However, these crabs are of little/no economic value, and they eventually die. A previous study reported that microsporidia were detected in the hepatopancreata of crabs infected with HPND (Ding et al., 2016). In addition, hepatopancreatic injury caused by environmental toxicants is believed to be one of the main causes of HPND. However, the etiology of HPND is unknown.

The hepatopancreas of crustaceans functions in both digestion and absorption (Wang et al., 2004b). A dense microbial colonization has been observed in the hepatopancreas of several isopods (Wood and Griffiths, 1988; Zimmer et al., 2001; Wang et al., 2004b) and the shrimp Neocaridina denticulate (Cheung et al., 2015). In addition to nutrition and digestion, the immanent microbiota play crucial roles in the immune response of their animal hosts (Cheung et al., 2015), and an imbalance of host–microbiota homeostasis has been shown to be responsible for certain illnesses (Chen et al., 2015). A previous study has reported symbiotic bacteria in crabs to belong to the phyla Bacteroidetes, Proteobacteria, Firmicutes, and Tenericutes (Li et al., 2007; Givens et al., 2013; Chen et al., 2015; Zhang et al., 2016). Of these symbiotic bacteria, Tenericutes are unique, because they lack a cell wall, are small in size, and possess a reduced genome (Ryan and Ray, 2004). Furthermore, bacteria of the genus Candidatus Hepatoplasma, which belong to the Tenericutes phylum, were identified in symbionts colonizing the midgut glands of terrestrial isopods, and a positive correlation was identified between the survival of hosts and the ingestion of low-quality food (Wang et al., 2004a, 2007; Fraune and Zimmer, 2008; Leclercq et al., 2014). In addition, bacteria of the phylum Tenericutes have been shown to associate with several crustacean diseases, such as mycoplasma spp., which was isolated from the moribund prawn Peneaux monodon (Ghadersohi and Owens, 1999); spiroplasma, which was linked to tremor disease in the Chinese mitten crab (Wang et al., 2003); and Acholeplasma sp., which was implicated in clearwater disease of the mud crab Scylla serrata (Chen et al., 2011).

In this study, meta-transcriptomic libraries were constructed from the hepatopancreata of crabs with and without HPND, and viral, microsporidial, and bacterial communities were compared between the two groups. We report herein on statistically significant changes in the microbiota of diseased and healthy crabs. Our data provide valuable insights on the etiology of HPND.

Materials and Methods

Sample Collection

In June 2015, two crabs with clinical symptoms typical of HPND and two healthy crabs were obtained from an aquatic breeding pond in Yandu District, Yancheng City, Jiangsu Province, China for meta-transcriptomic sequencing. In June 2016, more than 100 crabs exhibiting clinical symptoms typical of HPND were obtained from several aquatic breeding ponds in Jiangdu District, Yangzhou City and Yandu District, Yancheng City, Jiangsu Province, China.

cDNA Library Construction and Illumina Sequencing

Total RNA from the hepatopancreata of two crabs infected with HPND, as well as two healthy crabs, was isolated using the RNeasy^® Plus Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's protocol. To remove genomic DNA, RNA was treated with RNase-free DNase (Qiagen, Germany) according to the manufacturer's protocol. The quality and quantity of the total RNA was estimated with a NanoDrop 2000 Spectrophotometer (Thermo Scientific, USA) and an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). RNA integrity was assessed by electrophoresis on a 1% agarose gel. In order to ensure that there was no contamination of the genomic DNA, 1 μg RNA was used as PCR template, the expression of the β-actin gene of Eriocheir sinensis (accession number: HM053699.1) was selected as a control, using the primer pair EsACTIN-F (GCATCCACGAGACCACTTACA) and EsACTIN-R (CTCCTGCTTGCTGATCCACATC). The PCR program was as follows: denaturation at 94°C for 3 min; 35 cycles of 94°C for 30 s, 58°C for 30 s, and 72°C for 1 min; and 72°C for 10 min. Ribosomal RNA (rRNA) was depleted using the Ribo-Zero™ Magnetic Kit (Epicenter, Charlotte, NC, USA). The cDNA library was prepared using the TruSeq™ RNA Sample Prep Kit (Illumina, San Diego, CA, USA). The quality of the cDNA library was assessed using an Agilent 2100 Bioanalyzer. Sequencing was carried out on a HiSeq 2500 Ultra-High-Throughput Sequencer using the Mid Output Kit (both from Illumina), and 150-bp pair-end reads were obtained for each run.

Meta-Transcriptomic Data Analysis

The analysis of the entire meta-transcriptome was carried out as follows: the quality of each raw read (Phred score) was evaluated using the FastQC toolkit (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), and the adaptor was eliminated with SeqPrep software (https://github.com/jstjohn/SeqPrep). Thereafter, low-quality bases (Phred score < 20) were trimmed, and reads shorter than 50 bp were discarded with Sickle software (https://github.com/najoshi/sickle). rRNA reads were discarded after alignment to SILVA SSU (16S/18S) and SILVA LSU (23S/28S) databases with SortMeRNA software (http://bioinfo.lifl.fr/RNA/sortmerna/). The resulting high-quality reads were then used in the subsequent assembly. The mega-transcriptome was de novo assembled with Trinity software (http://trinityrnaseq.github.io/; version trinityrnaseq-r2013-02-25) using default parameters as previously described (Grabherr et al., 2011). ORFs were predicted using TransGeneScan software (http://sourceforge.net/projects/transgenescan/). Non-redundant gene catalogs were constructed with an identity of 95% and a coverage of 90% using CD-HIT software (http://www.bioinformatics.org/cd-hit/). To evaluate the expression level of each transcript, the FPKM (fragments per kilobase of exon per million fragments mapped) value of each transcript was obtained with RSEM (RNASeq by expectation maximization) software (http://deweylab.biostat.wisc.edu/rsem/). Differentially expressed genes were identified using edgeR (empirical analysis of digital gene expression data in R) software (http://www.bioconductor.org/packages/release/bioc/html/edgeR.html; Reiner et al., 2003; Robinson and Smyth, 2007). For this analysis, the filtering threshold was set to an FDR (false discovery rate) < 0.5 and a |log2FC| > 1.

All genes were characterized with BALSTX comparisons against the integrated NCBI NR database with an expectation value of 1e-5 (BLAST Version 2.2.28+, http://blast.ncbi.nlm.nih.gov/Blast.cgi). Species information was obtained from the respective taxonomy annotation NR database. The species abundance in each specimen was determined by calculating the FPKM value of each transcript in the respective species. The taxonomic abundance in each specimen was calculated at the domain, kingdom, phylum, class, order, family, and genus levels. Thus, the abundance profiles at the corresponding taxonomic levels were built. Differences in the abundance profile at each level between two groups were identified by Welch's t-test using STAMP software (http://kiwi.cs.dal.ca/Software/STAMP).

The meta-transcriptomic sequences were deposited into the NCBI Sequence Read Archive (SRA) under accession numbers SRR4308592, SRR4308642, SRR4330895, and SRR4333242.

Phylogenetic Analysis

A maximum likelihood tree was constructed from the alignment of 12 different 16S rRNA sequences obtained from the GenBank database with MEGA (Molecular Evolutionary Genetics Analysis, version 5) software (Tamura et al., 2007). The names and accession numbers of the related bacterial 16S genes were as follows: Candidatus Hepatoplasma crinochetorum (Accession number: AY500249.1), Mycoplasma pneumonia (AB680604.1), Mycoplasma ovipneumoniae (NR_025989.1), Ureaplasma parvum (NR_074762.1), Ureaplasma cati (NR_115604.1), Spiroplasma clarkii (NR_104750.1), Spiroplasma eriocheiris (NR_125517.1), Mesoplasma florum (AB681225.1), Entomoplasma ellychniae (NR_104951.1), Entomoplasma luminosum (AY155670.1), and Mesoplasma grammopterae (AY174170.1).

DNA Extraction and Polymerase Chain Reaction

The 16S rRNA gene from C. Hepatoplasma was assembled using rRNA raw reads and verified by PCR (polymerase chain reaction) using the primer pair Es.myco.F01 (5′-AGGGTTTGATTATGGCTCAGGA-3′) and Es.myco.R01 (5′-ACAAGACCAGAGAACGTATTCACC-3′). Based on the sequence of the 16S rRNA gene from C. Hepatoplasma, a primer pair (Es.myco.F02: 5′-ACTCCTACGGGAGGCAGCAG-3′ and Es.myco.R02: 5′-GCGGCTGCTGGCACATAGTT-3′) was designed to detect this species in the hepatopancreata of crabs infected with HPND. DNA was isolated using the Genomic DNA Extraction Kit according to the manufacturer's protocol (Tiangen, Beijing, China).

Results

Reads and Assembly Statistics

Four meta-transcriptomic libraries were constructed to identify differences in organism communities, especially the bacterial and viral communities from the hepatopancreata of crabs with and without HPND. A total of 195,152,066 reads were generated from four samples, and 162,963 transcripts were assembled (Tables S1, S2). The expression levels of entire transcripts were evaluated (Table S3). The distribution of transcripts to general functional categories was assessed on the basis of best BLAST matches to the COG database, GO database, and KEGG database (Tables S4– S6).

Taxonomic Data and Comparison of Crabs with and without HPND

Sequences were classified into 60 phyla, which included 16 bacterial phyla. Ten and 11 phyla were found to exist in the hepatopancreata of two crabs with HPND, while 12 and 16 phyla were found to exist in the hepatopancreata of two crabs without HPND. This analysis showed that the diversity of the organism communities in the hepatopancreata of healthy crabs was greater than that in diseased crabs (Table S7). Similar results were found at class, order, family, and genus levels. For example, 572 and 929 genera were identified in crabs with and without HPND, respectively, whereas 544 genera were identified in both groups (Table S8).

Among the 16 bacterial phyla, we identified four dominant phyla, namely, Bacteroidetes (41.04%), Firmicutes (26.43%), Proteobacteria (26.03%), and Actinobacteria (4.03%) in the hepatopancreata of crabs without HPND. We also identified five dominant phyla, namely, Tenericutes (31.82%), Firmicutes (21.38%), Proteobacteria (25.45%), Actinobacteria (16.70%), and Bacteroidetes (3.49%) in the hepatopancreata of crabs with HPND. The percentage of bacteria belonging to the Tenericutes phylum was higher in diseased crabs (31.82%) than in healthy crabs (0.02%; P < 0.01; Figure 1 and Tables S9, S10). By contrast, the percentage of bacteria belonging to the Bacteroidetes phylum was lower (3.49%) in diseased crabs than in healthy crabs (41.04%). There was a significant decrease in bacteria belonging to the Bacteroidetes phylum in Bacteroidia and Flavobacteria classes in crabs with HPND (Figure 2). In addition, the percentage of bacteria belonging to the Actinobacteria phylum was higher in crabs with HPND (16.70%) than in crabs without the disease (4.03%).

FIGURE 1

Figure 1. Relative abundance of bacteria at the phylum level in crabs with and without HPND based on meta-transcriptomic data. HPND #01, HPND #02: two samples from crabs with HPND; Healthy #01, Healthy #02: two samples from crabs without HPND.

FIGURE 2

Figure 2. Relative abundance of bacteria from Bacteroidetes at the class level in crabs with and without HPND. HPND #01, HPND #02: two samples from crabs with HPND; Healthy #01, Healthy #02: two samples from crabs without HPND.

Eleven families of viruses were identified from the analysis of the meta-transcriptome data. The expression levels of these viruses were low in the hepatopancreata of crabs with and without HPND. Five families of virus were present in only one of the crab samples. The remaining six families of virus, which were present in both groups, were compared, and there were no significant differences between groups (Figure 3). We also identified a new nodavirus, whose expression was higher in diseased crabs than in healthy crabs.

FIGURE 3

Figure 3. Relative abundance of viruses at the family level in crabs with and without HPND. HPND #01, HPND #02: two samples from crabs with HPND; Healthy #01, Healthy #02: two samples from crabs without HPND.

Five families of microsporidia, including three dominant families, were identified from the analysis of the meta-transcriptomic data. These microsporidia were found in the hepatopancreata of crabs with and without HPND, and there were no significant differences between groups (Figure 4).

FIGURE 4

Figure 4. Relative abundance of Microsporidia at the family level in crabs with and without HPND. HPND #01, HPND #02: two samples from crabs with HPND; Healthy #01, Healthy #02: two samples from crabs without HPND.

Mycoplasmataceae Were Abundant in Crabs with HPND

The different bacteria belonging to the Tenericutes phylum identified in the hepatopancreata of crabs with HPND was further studied by analyzing the meta-transcriptomic data. The majority of the Tenericutes bacteria (>93.8%) belonged to the Mycoplasmataceae family, with 72.9% belonging to the genus C. Hepatoplasma and 20.9% belonging to the genus Mycoplasma (Figure 5). In addition, the three genera Acholeplasma (3.05%), Spiroplasma (2.56%), and Entomoplasma (0.19%) were also found in Tenericutes bacterial communities in the hepatopancreata of diseased crabs.

FIGURE 5

Figure 5. Bacterial composition from Tenericutes at the genus level in crabs with and without HPND. (A) Comparison of absolute sequence counts for the major genera of Tenericutes in diseased and healthy crabs. (B) Comparison of the relative abundance of the major genera of Tenericutes in iseased and healthy crabs. HPND #01, HPND #02: two samples from crabs with HPND; Healthy #01, Healthy #02: two samples from crabs without HPND.

At the gene expression level, 173 protein-coding genes were found to belong to the Mycoplasmataceae family. Compared with healthy crabs, all genes were highly expressed in the hepatopancreata of diseased crabs (Table S11). A BLAST search revealed the majority of these genes (134 genes) to be highly homologous with those from Candidatus Hepatoplasma crinochetorum, a species in the genus C. Hepatoplasma, family Mycoplasmataceae. The remaining genes (39 genes) were homologous with genes from other species belonging to the Mycoplasmataceae family. For example, an alignment showed all 134 genes to be 33–93% homologous with their corresponding genes in Candidatus Hepatoplasma crinochetorum. The gene with the highest homology (93%) was rpoB. These results indicate that these genes are from a microorganism belonging to the genus C. Hepatoplasma. For the purpose of this study, this microorganism was designated as “C. Hepatoplasma sp.”

Phylogenetic Tree Analysis

The 16S rRNA gene is the most common gene; thus, it was used to identify and to characterize the taxonomic status of the new species described in this study. The partial sequence (1369 bp) of the 16S rRNA gene from C. Hepatoplasma was cloned from DNA isolated from the hepatopancreata of crabs with HPND. Sequence alignment showed that the 16S rRNA was 87% homologous with that of Candidatus Hepatoplasma crinochetorum. Phylogenic analysis indicated that C. Hepatoplasma was a close relative of Candidatus Hepatoplasma crinochetorum and clustered in one clade (Figure 6). These results confirmed that this microorganism classified under the Mycoplasmataceae family, which was abundant in crabs with HPND, was a new member of the genus C. Hepatoplasma.

FIGURE 6

Figure 6. Phylogenetic tree based on the 16S rRNA gene sequence showing the relationship between Candidatus Hepatopancreas and Tenericutes.

Polymerase Chain Reaction Detection

DNA isolated from more than 100 crabs with clinical symptoms typical of HPND was assayed by PCR using a primer pair specific for the 16S RNA gene of Candidatus Hepatopancreas. More than 98% of crabs were positive for HPND (Figure S1). Healthy crabs were also screened by PCR, and a weak/no positive band was found.

Discussion

Meta-transcriptomic analysis provides insights into disease etiology by cataloging and comparing sequences from suspected organisms (Rosales and Thurber, 2015). This approach is powerful, because it simultaneously evaluates both viral and bacterial communities, and even suspected pathogens from the Eukaryota domain. Considering that RNA viruses can associate with certain diseases in E. sinensis (Zhang et al., 2004; Zhang and Bonami, 2007, 2012; Shen et al., 2015), and microsporidia were found in crabs infected with HPND (Ding et al., 2016), we analyzed the meta-transcriptomes of hepatopancreata from crabs with and without HPND to identify potential pathogens from viruses, bacteria, and eukaryotes.

We studied the differences between the virome and the microbiome in the hepatopancreata of crabs with and without HPND. Although there was no significant difference in the virome of both groups, we identified a new nodavirus, which was more abundant in the hepatopancreata of diseased crabs than in healthy crabs. The bacterial community in the hepatopancreata of crabs without HPND was more species-rich than that in the hepatopancreata of crabs with HPND. Interestingly, similar findings were reported for the human oral and gut microbiomes during health and disease (Jorth et al., 2014; Ling et al., 2014). We identified four dominant phyla, namely, Proteobacteria, Bacteroidetes, Firmicutes, and Actinobacteria in the hepatopancreata of healthy crabs. A previous study has reported that 90–95% of the bacterial phylotypes in the intestine of E. sinensis were Proteobacteria and Bacteroidetes, and that Bacteroidetes were common to all crab types (Li et al., 2007). For example, the phyla Proteobacteria, Bacteroidetes, and Firmicutes were most abundant in the digestive system of the black tiger shrimp P. monodon (Shakibazadeh et al., 2009; Chaiyapechara et al., 2012). Similarly, bacteria belonging to the Bacteroidetes phylum were more plentiful in the hepatopancreas, whereas those of the Firmicutes phylum were more plentiful in the foregut and intestine of the shrimp N. denticulate (Cheung et al., 2015). In this study, Bacteroidetes, Proteobacteria, and Firmicutes were the most dominant bacterial phyla in the hepatopancreata of healthy crabs, with a combined prevalence of 93%. However, there was a significant difference in the prevalence of the major bacterial phyla in the hepatopancreata of crabs with and without HPND. Tenericutes was the most prevalent bacterial phylum in crabs with HPND crabs (31.82%), whereas its prevalence was low in healthy crabs (0.02%). We also found a significantly low proportion of Bacteroidetes in diseased crabs (3.49%) compared to healthy crabs (41.04%). Interestingly, microbial diversity is lower in obese individuals than in lean individuals, with obese individuals having a lower proportion of Bacteroidetes and a higher proportion of Actinobacteria (Turnbaugh et al., 2009). There was also a negative relationship between the abundance of Bacteroidetes and the presence of fat deposits (Guo et al., 2008). The higher proportion of Bacteroidetes in the hepatopancreas may be related to the higher cellulolytic activity in crabs (Cheung et al., 2015). In this study, a similar shift in the microbial composition was found in the hepatopancreata of crabs with HPND compared with healthy crabs. The evolutionarily-stable commensal microbes positively and negatively affect host health, and when this balance is disrupted, symbiotic microbes can induce disease (Breznak and Brune, 1994; Douglas, 1998; Fraune and Zimmer, 2008; Jorth et al., 2014). The lower proportion of Bacteroidetes and the higher proportion of Actinobacteria indicate that there were abnormalities in the metabolic breakdown of fat and cellulose in the hepatopancreata of crabs with HPND.

The majority of Tenericutes bacteria (>93.8%) belonged to the Mycoplasmataceae family, with 72.9% belonging to C. Hepatoplasma and 20.9% belonging to Mycoplasma. Mycoplasmataceae bacteria were found in the symbiotic bacterial community in the gill and the gut of E. sinensis. An alignment has shown that C. Hepatoplasma shares more than 95% sequence identity with uncultured Mycoplasmataceae detected in the mud crab Scylla paramamosain (Li et al., 2012), more than 94% sequence identity with uncultured Mycoplasmataceae bacteria in the carapace, gut, and hemolymph of the Atlantic Blue Crab Callinectes Sapidus (Givens et al., 2013), and more than 94% sequence identity with uncultured Mycoplasmataceae symbionts in the midgut of the isopod Ligia oceanic (Fraune and Zimmer, 2008). Similar results were reported for the CMC Mollicutes group I that was found in the midgut bacterial community of E. sinensis (Chen et al., 2015), suggesting that C. Hepatoplasma and the CMC Mollicutes group I might belong to the same species of Mycoplasmataceae bacteria. Mycoplasmataceae bacteria were not detected in the guts of crabs from the Yangtze River Estuary and Chongming Islands (Li et al., 2007). Tenericutes were the dominant phyla in the gut of the Chinese mitten crab (Chen et al., 2015). Mycoplasmataceae are one of the most dominant families in the gut; however, they are hardly detected in the gills or water (Zhang et al., 2016). These reports suggest that Mycoplasmataceae bacteria in E. sinensis are variable, and that Tenericutes are not indispensable to the crab. A previous study has shown that Candidatus Hepatoplasma crinochetorum was beneficial to its isopod host under low-nutrient conditions (Fraune and Zimmer, 2008). In this study, the abundance of Candidatus Hepatopancreas in the hepatopancreata of crabs with HPND was indicative of abnormalities in the absorption of nutrients.

We found that the expression of carboxylesterase family genes was significantly up-regulated (Table S12). Carboxylesterase is the major enzyme that breaks down the insecticide pyrethroid, which is widely used in Chinese mitten crab cell cultures, although it is known to harm the hepatopancreas (Wheelock et al., 2006). Further studies are needed to determine how pyrethroid injures the heptopancreas, thus leading to abnormalities in metabolism, nutrient absorption, microbial dysbiosis, and HPND in crabs.

In addition to Mycoplasmataceae, the genera Acholeplasma and Spiroplasma, also from Mullicutes, were abundant in the hepatopancreata of crabs with HPND. Acholeplasma and Spiroplasma bacteria are responsible for crab disease (Wang W. et al., 2004; Chen et al., 2011). An alignment of sequences showed that Spiroplasma and Spiroplasma eriocheir, which is believed to cause crab tremor disease, are two different species. Further studies are needed to determine whether Acholeplasma and Spiroplasma associate with HPND infection in crabs.

Microsporidia can also associate with HPND (Ding et al., 2016). Analysis of meta-transcriptome data indicate that the microsporidial community included three major families, namely, Enterocytozoonidae, Nosematidae, and Unikaryonidae. These microsporidia were found in the hepatopancreata of crabs with and without HPND, although there were no significant differences between groups. Regardless, these data indicate that microsporidia are unlikely to cause HPND.

In conclusion, there were statistically significant changes in the microbiota of HPND-infected crabs. Specifically, we found an increased prevalence of bacteria belonging to the Tenericutes and Actinobacteria phyla and a decreased prevalence of bacteria belonging to the Bacteroidetes phylum, thus reflecting abnormalities in metabolism and the absorption of nutrients in HPND-infected crabs and indicating a connection with HPND.

Author Contributions

HS performed the experiment work and statistical analysis, produced the tables and figures, and wrote the paper. YZ, KS, YM, and TD assist the experiment. AS read and made improvements to the manuscript.

Funding

This study was supported by the Aquatic Three Update Project of Jiangsu Province (Y2016-35) and the Central Public-Interest Scientific Institution Basal Research Fund, Freshwater Fisheries Research Center, CAFS (2017JBFM01).

Conflict of Interest Statement

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.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2017.00732/full#supplementary-material

Figure S1. PCR detection of the 16S rRNA gene of Candidatus Hepatopancreas in the hepatopancreata of different crabs. 01–26: crabs with HPND; ctr01, ctr02: crabs without HPND; M, 500 bp DNA Ladder.

Table S1. Statistics for sequencing data.

Table S2. Statistics of assembly result.

Table S3. Gene expression level in different samples.

Table S4. The distribution of transcripts to general functional categories was assessed on the basis of best BLAST matches to the COG database.

Table S5. The distribution of transcripts to general functional categories was assessed on the basis of best BLAST matches to the GO database.

Table S6. The distribution of transcripts to general functional categories was assessed on the basis of best BLAST matches to the KEGG database.

Table S7. Abundance of phylum in different samples.

Table S8. Abundance of genus in different samples.

Table S9. T-test between HPND group and Healthy grup in phylum level.

Table S10. T-test between HPND group and Healthy grup in genus level.

Table S11. 173 protein-coding genes were annotated belong to the Mycoplasmataceae family.

Table S12. Expression level of five carboxylesterase genes in crabs with and without HPND.

References

Breznak, J., and Brune, A. (1994). Role of microorganisms in the digestion of lignocellulose by termites. Annu. Rev. Entomol. 39, 453–487. doi: 10.1146/annurev.en.39.010194.002321

CrossRef Full Text | Google Scholar

Chaiyapechara, S., Rungrassamee, W., Suriyachay, I., Kuncharin, Y., Klanchui, A., Karoonuthaisiri, N., et al. (2012). Bacterial community associated with the intestinal tract of P. monodon in commercial farms. Microb. Ecol. 63, 938–953. doi: 10.1007/s00248-011-9936-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J. G., Lou, D., and Yang, J. F. (2011). Isolation and Identification of Acholeplasma sp. from the Mud Crab, Scylla serrata. Evid. Based Complement. Alternat. Med. 2011:209406. doi: 10.1155/2011/209406

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Di, P., Wang, H., Li, B., Pan, Y., Yan, S., et al. (2015). Bacterial community associated with the intestinal tract of Chinese mitten crab (Eriocheir sinensis) farmed in Lake Tai, China. PLoS ONE 10:e0123990. doi: 10.1371/journal.pone.0123990

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheung, M. K., Yip, H. Y., Nong, W., Law, P. T., Chu, K. H., Kwan, H. S., et al. (2015). Rapid change of microbiota diversity in the gut but not the Hepatopancreas during gonadal development of the new shrimp model Neocaridina denticulata. Mar. Biotechnol. 17, 811–819. doi: 10.1007/s10126-015-9662-8

CrossRef Full Text | Google Scholar

Ding, Z., Meng, Q., Liu, H., Yuan, S., Zhang, F., Sun, M., et al. (2016). First case of hepatopancreatic necrosis disease in pond-reared Chinese mitten crab, Eriocheir sinensis, associated with microsporidian. J. Fish Dis. 39, 1043–1051. doi: 10.1111/jfd.12437

PubMed Abstract | CrossRef Full Text | Google Scholar

Douglas, A. E. (1998). Nutritional interactions in insect-microbial symbioses: aphids and their symbiotic bacteria Buchnera. Annu. Rev. Entomol. 43, 17–37. doi: 10.1146/annurev.ento.43.1.17

PubMed Abstract | CrossRef Full Text | Google Scholar

Fraune, S., and Zimmer, M. (2008). Host-specificity of environmentally transmitted Mycoplasma-like isopod symbionts. Environ. Microbiol. 10, 2497–2504. doi: 10.1111/j.1462-2920.2008.01672.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghadersohi, A., and Owens, L. (1999). Isolation, characterisation and DNA analysis of Mycoplasma spp. from moribund prawns Penaeus monodon cultured in Australia. Dis. Aquat. Organ. 35, 53–61. doi: 10.3354/dao035053

CrossRef Full Text | Google Scholar

Givens, C. E., Burnett, K. G., Burnett, L. E., and Hollibaugh, M. B. T. (2013). Microbial communities of the carapace, gut, and hemolymph of the Atlantic blue crab, Callinectes sapidus. 160, 2841–2851. doi: 10.1007/s00227-013-2275-8

CrossRef Full Text | Google Scholar

Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Xia, X., Tang, R., Zhou, J., Zhao, H., and Wang, K. (2008). Development of a real-time PCR method for Firmicutes and Bacteroidetes in faeces and its application to quantify intestinal population of obese and lean pigs. Lett. Appl. Microbiol. 47, 367–373. doi: 10.1111/j.1472-765X.2008.02408.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jorth, P., Turner, K. H., Gumus, P., Nizam, N., Buduneli, N., and Whiteley, M. (2014). Metatranscriptomics of the human oral microbiome during health and disease. MBio 5, e01012–e01014. doi: 10.1128/mBio.01012-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Leclercq, S., Dittmer, J., Bouchon, D., and Cordaux, R. (2014). Phylogenomics of “Candidatus Hepatoplasma crinochetorum,” a lineage of mollicutes associated with noninsect arthropods. Genome Biol. Evol. 6, 407–415. doi: 10.1093/gbe/evu020

CrossRef Full Text | Google Scholar

Li, K., Guan, W., Wei, G., Liu, B., Xu, J., Zhao, L., et al. (2007). Phylogenetic analysis of intestinal bacteria in the Chinese mitten crab (Eriocheir sinensis). J. Appl. Microbiol. 103, 675–682. doi: 10.1111/j.1365-2672.2007.03295.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Sun, L., Wu, H., Hu, Z., Liu, W., Li, Y., et al. (2012). The intestinal microbial diversity in mud crab (Scylla paramamosain) as determined by PCR-DGGE and clone library analysis. J. Appl. Microbiol. 113, 1341–1351. doi: 10.1111/jam.12008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ling, Z., Liu, X., Jia, X., Cheng, Y., Luo, Y., Yuan, L., et al. (2014). Impacts of infection with different toxigenic Clostridium difficile strains on faecal microbiota in children. Sci. Rep. 4:7485. doi: 10.1038/srep07485

PubMed Abstract | CrossRef Full Text | Google Scholar

Reiner, A., Yekutieli, D., and Benjamini, Y. (2003). Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics 19, 368–375. doi: 10.1093/bioinformatics/btf877

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, M. D., and Smyth, G. K. (2007). Moderated statistical tests for assessing differences in tag abundance. Bioinformatics 23, 2881–2887. doi: 10.1093/bioinformatics/btm453

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosales, S. M., and Thurber, R. V. (2015). Brain meta-transcriptomics from harbor seals to infer the role of the microbiome and virome in a stranding event. PLoS ONE 10:e0143944. doi: 10.1371/journal.pone.0143944

CrossRef Full Text | Google Scholar

Ryan, K., and Ray, C. (2004). Sherris Medical Microbiology, 4th Edn. New York, NY: McGraw Hill.

Shakibazadeh, S., Saad, C., Christianus, A., Kamarudin, M., Sijam, K., Shamsudin, M., et al. (2009). Bacteria flora associated with different body parts of hatchery reared juvenile Penaeus monodon, tanks water and sediment. Ann. Microbiol. 59, 425–430. doi: 10.1007/BF03175126

CrossRef Full Text | Google Scholar

Shen, H., Hu, Y., and Zhou, X. (2014). Sex-lethal gene of the Chinese mitten crab Eriocheir sinensis: cDNA cloning, induction by eyestalk ablation, and expression of two splice variants in males and females. Dev. Genes Evol. 224, 97–105. doi: 10.1007/s00427-014-0467-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, H., Ma, Y., and Hu, Y. (2015). Near-full-length genome sequence of a novel reovirus from the Chinese mitten crab, Eriocheir sinensis. Genome Announc. 3, e00447–e00415. doi: 10.1128/genomeA.00447-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, K., Dudley, J., Nei, M., and Kumar, S. (2007). MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol. Biol. Evol. 24, 1596–1599. doi: 10.1093/molbev/msm092

PubMed Abstract | CrossRef Full Text | Google Scholar

Turnbaugh, P. J., Hamady, M., Yatsunenko, T., Cantarel, B. L., Duncan, A., Ley, R. E., et al. (2009). A core gut microbiome in obese and lean twins. Nature 457, 480–484. doi: 10.1038/nature07540

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Rong, L., Gu, W., Du, K., and Chen, J. (2003). Study on experimental infections of Spiroplasma from the Chinese mitten crab in crayfish, mice and embryonated chickens. Res. Microbiol. 154, 677–680. doi: 10.1016/j.resmic.2003.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Wen, B., Gasparich, G. E., Zhu, N., Rong, L., Chen, J., et al. (2004). A spiroplasma associated with tremor disease in the Chinese mitten crab (Eriocheir sinensis). Microbiology 150, 3035–3040. doi: 10.1099/mic.0.26664-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Brune, A., and Zimmer, M. (2007). Bacterial symbionts in the hepatopancreas of isopods: diversity and environmental transmission. FEMS Microbiol. Ecol. 61, 141–152. doi: 10.1111/j.1574-6941.2007.00329.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Stingl, U., Anton-Erxleben, F., Geisler, S., Brune, A., and Zimmer, M. (2004a). “Candidatus Hepatoplasma crinochetorum,” a new, stalk-forming lineage of Mollicutes colonizing the midgut glands of a terrestrial isopod. Appl. Environ. Microbiol. 70, 6166–6172. doi: 10.1128/AEM.70.10.6166-6172.2004

CrossRef Full Text | Google Scholar

Wang, Y., Stingl, U., Anton-Erxleben, F., Zimmer, M., and Brune, A. (2004b). ‘Candidatus Hepatincola porcellionum’ gen. nov., sp. nov., a new, stalk-forming lineage of Rickettsiales colonizing the midgut glands of a terrestrial isopod. Arch. Microbiol. 181, 299–304. doi: 10.1007/s00203-004-0655-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wheelock, C. E., Miller, J. L., Miller, M. J., Phillips, B. M., Huntley, S. A., Gee, S. J., et al. (2006). Use of carboxylesterase activity to remove pyrethroid-associated toxicity to Ceriodaphnia dubia and Hyalella azteca in toxicity identification evaluations. Environ. Toxicol. Chem. 25, 973–984. doi: 10.1897/05-334R.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wood, S., and Griffiths, B. (1988). Bacteria associated with the hepatopancreas of the woodlice Oniscus asellus and Porcellio scaber (Crustacea, Isopoda). Pedobiologia 31, 89–94.

Google Scholar

Zhang, M., Sun, Y., Chen, L., Cai, C., Qiao, F., Du, Z., et al. (2016). Symbiotic bacteria in gills and guts of Chinese mitten crab (Eriocheir sinensis) differ from the free-living bacteria in water. PLoS ONE 11:e0148135. doi: 10.1371/journal.pone.0148135

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., and Bonami, J. R. (2007). A roni-like virus associated with mortalities of the freshwater crab, Eriocheir sinensis Milne edwards, cultured in China, exhibiting ‘sighs disease’ and black gill syndrome. J. Fish Dis. 30, 181–186. doi: 10.1111/j.1365-2761.2007.00796.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., and Bonami, J. R. (2012). Isolation and partial characterization of a new reovirus in the Chinese mitten crab, Eriocheir sinensis H Milne edwards. J. Fish Dis. 35, 733–739. doi: 10.1111/j.1365-2761.2012.01398.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Shi, Z., Zhang, J., and Bonami, J. R. (2004). Purification and characterization of a new reovirus from the Chinese mitten crab, Eriocheir sinensis. J. Fish Dis. 27, 687–692. doi: 10.1111/j.1365-2761.2004.00587.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zimmer, M., Danko, J., Pennings, S., Danford, A., Ziegler, A., Uglow, R., et al. (2001). Hepatopancreatic endosymbionts in coastal isopods (Crustacea: Isopoda), and their contribution to digestion. Mar. Biol. 138, 955–963. doi: 10.1007/s002270000519

CrossRef Full Text | Google Scholar