Figures
Abstract
Background
Massive infection caused by oomycete fungus Saprolegnia parasitica is detrimental to freshwater fish. Recently, we showed that copper sulfate demonstrated good efficacy for controlling S. parasitica infection in grass carp. In this study, we investigated the mechanism of inhibition of S. parasitica growth by copper sulfate by analyzing the transcriptome of copper sulfate—treated S. parasitica. To examine the mechanism of copper sulfate inhibiting S. parasitica, we utilized RNA-seq technology to compare differential gene expression in S. parasitica treated with or without copper sulfate.
Results
The total mapped rates of the reads with the reference genome were 90.50% in the control group and 73.50% in the experimental group. In the control group, annotated splice junctions, partial novel splice junctions and complete novel splice junctions were about 83%, 3% and 14%, respectively. In the treatment group, the corresponding values were about 75%, 6% and 19%. Following copper sulfate treatment, a total 310 genes were markedly upregulated and 556 genes were markedly downregulated in S. parasitica. Material metabolism related GO terms including cofactor binding (33 genes), 1,3-beta-D-glucan synthase complex (4 genes), carboxylic acid metabolic process (40 genes) were the most significantly enriched. KEGG pathway analysis also determined that the metabolism-related biological pathways were significantly enriched, including the metabolic pathways (98 genes), biosynthesis of secondary metabolites pathways (42 genes), fatty acid metabolism (13 genes), phenylalanine metabolism (7 genes), starch and sucrose metabolism pathway (12 genes). The qRT-PCR results were largely consistent with the RNA-Seq results.
Citation: Hu K, Ma R-R, Cheng J-M, Ye X, Sun Q, Yuan H-L, et al. (2016) Analysis of Saprolegnia parasitica Transcriptome following Treatment with Copper Sulfate. PLoS ONE 11(2): e0147445. https://doi.org/10.1371/journal.pone.0147445
Editor: Jason E. Stajich, University of California-Riverside, UNITED STATES
Received: April 29, 2015; Accepted: January 3, 2016; Published: February 19, 2016
Copyright: © 2016 Hu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are available via Figshare (https://dx.doi.org/10.6084/m9.figshare.2062170.v1).
Funding: This study was supported by the 863 Program (Grant 2011AA10A216), the Special Fund for Agro-scientific Research in the Public Interest (Grant 201203085), the National Natural Resources Platform and the Shanghai University Knowledge Service Platform. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Saprolegnia parasitica, a destructive oomycete pathogen that infects fish and fish eggs, is detrimental to freshwater hatcheries [1, 2]. It often causes a severe disease, characterized by visible white or grey patches of filamentous mycelium on the freshwater fish or their eggs, which seriously affect their value, and can result in significant economic loss [3]. S. parasitica infection in aquaculture is prevalent, particularly under stress conditions such as rapid changing of temperatures or water quality [3]. The infection of S. parasitica can be controlled with effective substances. However, malachite green, an organic dye that is very efficient at controlling S. parasitica infection [4], was banned worldwide in 2002 [3] because of its carcinogenic and toxicological effects on animal health [5, 6]. Therefore, it is imperative to find an alternative.
Humic substances with higher molecular weights and aromaticity contain a high number of organic radicals that are efficient for reducing fungal growth. Humic substance development of internal oxidative stress could be the mechanism for the observed inhibition of S. parasitica growth [7]. Clotrimazole is another potent treatment agent of S. parasitica infection by inhibiting sterol 14 alpha-demethylase (CYP51)[8]. Meanwhile, Saprolmycin A-E are reported to have highly effective anti–S. parasitica selective activity as well [9]. Recently, our lab found copper sulfate could control S. parasitica infection in grass carp with fairly good efficacy, which suggested that copper sulfate might be used as a drug additive to control S. parasiticainfection in the aquaculture industry [10]. However, the mechanism of inhibition of S. parasitica growth by copper sulfate is unclear and it is necessary to study.
In this study, we performed transcriptome analysis to investigate the mechanism of copper sulfate inhibition of S. parasitica growth. We found that copper sulfate might inhibit protein synthesis, energy biogenesis, and metabolism in S. parasitica. Our findings may describe the mechanism of action of copper sulfate for controlling S. parasitica growth in water.
Materials and Methods
The experimental protocol was established, according to the ethical guidelines of the Helsinki Declaration and was approved by the Human Ethics Committee of Shanghai Ocean University, China.
Maintenance and treatment of S. parasitica
Saprolegnia parasitica (ATCC 200013) was obtained from American Type Culture Collection (Manassas, USA) and cultured on potato dextrose agar (Sinopharm Chemical Reagent Co., Ltd., Shanghai, China) for 4 days at 20°C. A mold colony, composed of the countless S. parasitica, was treated with copper sulfate (0.5 mg/L, Shanghai Guoyao Chemical Reagent Co. Shanghai, China) for 30 min. Additional S. parasitica were treated with sterile H2O as a control. Afterward, he S. parasitica were collected by centrifugation and stored in liquid nitrogen for further experiment.
RNA isolation, RNA sequencing (RNA-seq) library construction and sequencing
We extracted total RNA from samples using TRIzol reagent (Invitrogen, USA) according to the manufacturer’s instruction. We removed DNA contaminants after treatment with RNase-free DNase I (Takara Biotechnology, Dalian, China), and confirmed RNA integrity with a minimum RNA integrated number value of 8 by an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, USA). We isolated poly(A) mRNA using oligo-dT beads and then treated it with fragmentation buffer. We then transcribed the cleaved RNA fragments into first-strand complementary DNA (cDNA) using reverse transcriptase and random hexamer primers, followed by second-strand cDNA synthesis using DNA polymerase I and RNase H. The double-stranded cDNA was end-repaired using T4 DNA polymerase, Klenow fragment, and T4 polynucleotide kinase, followed by the addition of a single (A) base using Klenow 39–59 exo2 polymerase, and ligation with an adapter or index adapter using T4 quick DNA ligase. We selected adaptor-ligated fragments according to their size on agarose gel. The desired cDNA fragment size range was excised from the gel, and PCR was performed to selectively enrich and amplify the fragments. Following validation on an Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System (ABI, California, USA), the cDNA library was sequenced on a flow cell using high-throughput 101-bp pair-end mode on an Illumina HiSeq 2500 (Illumina, San Diego, USA) unit.
Sequencing, data processing, and quality control
We filtered low-quality reads and removed 3′ adapter sequences using Trim Galore. The obtained reads were cleaned using FastQC software (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), and we evaluated the content and quality of the nucleotide bases in the sequencing data. Next, we conducted a comparative analysis with the reference genome (ASM15154v2). For each sample, sequence alignment with the reference genome sequences was carried out using Tophat [11].
Identification of differential expression genes
To investigate the expression level of each transcript of the two different treatment groups, we used Cuffnorm, a program to estimate the expression level (relative abundance) of a specific transcript expressed using FRKM (Fragments per kilobase of transcript per million fragments mapped) [12]. Generally, when FPKM > 0.1, this indicates that a given transcript is expressed. The expression level of each transcript was transformed using base 2 log2(FPKM+1). Meanwhile, we used the Cuffdiff program to calculate the fold change of a transcript and to screen all differentially expressed genes (DEGs). Two-fold changes with p-values < 0.05 were considered significant.
GO functional annotation and enrichment analysis for differentially expressed genes
Initially, DEGs were annotated against the UniProt database (http://www.uniprot.org/). Next, we analyzed the functional annotation by gene ontology terms (GO; http://www.geneontology.org) with Blast 2GO (https://www.blast2go.com/) [13].
All DEGs were mapped to GO terms in the GO database, and gene numbers were calculated for every term, followed by an ultra-geometric test to find significantly enriched GO terms in DEGs compared to the transcriptome background. The formula used was as follows: Where N is the number of all genes with GO annotation; n is the number of DEGs in N; M is the number of all genes annotated to specific GO terms; m is the number of DEGs in M. The calculated p-value was subjected to Bonferroni correction. A corrected p-value < 0.05 was defined as a threshold, and then GO terms were considered significantly enriched in the DEGs.
DEGs pathway analysis
DEGs were annotated to the KEGG pathways database using the online KEGG Automatic Annotation Server (KAAS, available online: http://www.genome.jp/keg/kaas/). Enriched DEG pathways were identified using the same formula as that in GO analysis. In KEGG pathway analysis, N is the number of all genes with KEGG annotation, n is the number of DEGs in N, M is the number of all genes annotated to specific pathways, and m is the number of DEGs in M.
DEG verification using qRT-PCR
The expression levels of DEGs identified in the RNA-seq analysis were verified with quantitative RT-PCR (qRT-PCR). Primers were designed using the Primer 5.0 software(Premier company, Canada), and SpTub-b (S.parasitica Tub-b)was used as the reference gene [14, 15].
The reactions were performed in a 25 μl volume, composed of 2μl cDNA, 0.5μl of both the forward and reverse primer (10μM), 12.5μl SYBR Premix Ex Taq (2×) and 9.5μl RNase-free H2O. The thermal cycling program was 95°C for 30 s, followed by 40 cycles of 95°C for 5 s, 60°C for 30s, and 72°C for 30 s. Melting curve analysis was performed at the end of the qRT-PCR cycles to confirm the PCR specificity. Three replications were carried out and analysis of relative gene expression data using 2−△△CT Method [16].
Results
Illumina sequencing and quality assessment
We performed RNA-seq using the Illumina sequencing platform to examine the effect of copper sulfate on the S. parasitica transcriptome. After filtering and quality checks of the raw reads (47,614,574 and 51,163,206 for the control and treatment groups, respectively), there were approximately 44 million (44,676,046) and 45 million (45,325,202) trimmed reads with trim rates of 95.0% and 88.6% for control and the treatment samples, respectively. Meanwhile, the respective average length of reads was 94.5 and 93.8 bp (Table 1), indicating successful sequencing of the S. parasitica transcriptome. And trimmed reads were used for the subsequent analysis.
Comparative analysis with reference genome
The trimmed reads of S. parasitica transcriptome were compared with the reference genome sequence. The total mapped rates of the reads with the reference genome were 90.50% in the control group and 73.50% in the experimental group. Uniquely mapped reads were about 14 million (14,952,142) for the control and 29 million (29,202,690) for the experimental group, and accounted for 33.47% and 64.42% of total reads, respectively. Multiple mapped reads were about 25 million (25,476,517) for the control group and 4 million (4,116,175) of the experimental group, and make up 57.03% and 9.08% of total reads, respectively. Total unmapped reads accounted for 9.50% in the control group and 26.50% in the experimental group (Table 2).
Additionally, splice junctions were discovered using TopHat, which can identify splice variants of genes. In the control group, known splice junctions accounted for 83% of total splice junctions annotated in the reference genome sequence. Partial novel splice junctions were about 3% of total splice junctions, and complete novel splice junctions were about 14% of total splice junctions (Fig 1A). In the experimental group, known splice junctions, partial novel splice junctions and complete novel splice junctions were about 75%, 6% and 19%, respectively (Fig 1B). These results indicate that the transcriptome of S. parasiticatreated by copper sulfate contained more partial novel and complete novel splice junctions than the S. parasitica of control group.
Pie charts break down the known, partial novel and novel splice junctions identified in this study. Known splice junctions: both ends of splice junction belong to the annotated splice junction in the known genetic model; partial novel splice junctions: one end of splice junction belongs to the annotated splice junction in the known genetic model and the other end belongs to the new splice junction; complete novel splice junctions:both ends of splice junction belong to the new splice junction.
Analysis of DEGs
To identify DEGs in S. parasitica following copper sulfate treatment, we used Cuffdiff program to generate S. parasitica gene expression profiles (Figs 2 and 3). The program identified 310 genes were markedly upregulated and 556 genes were markedly down regulated in S. parasitica following copper sulfate treatment, which indicates that copper sulfate affects S. parasitica gene expression.
Volcanic plot of the degree of differences in the expression profile of S. parazitica after treatment with copper sulfate. X-axis, log2(fold change); Y-axis,-log10(Pvalue). Red, the significantly upregulated genes, green, the significantly downregulated genes. Each dot represents one gene.
Red, upregulated genes; green, downregulated genes. Each line represents one gene.
GO annotation of DEGs
To investigate the biological functions in which DEGs are involved in S. parasitica, we classified the DEGs according to GO classification following copper sulfate treatment in S. parasitica. This analysis identifies the main biological functions DEGs exercise. The GO functional enrichment analysis also involved cluster analysis of expression patterns. Thus, the expression patterns of DEGs annotated with a given GO term were easily obtained. All annotated genes were classified into three GO domains: biological processes, cellular component, and molecular function. Dissimilar expression profiles were obtained from the DEGs in the treated and control samples, revealing the obvious effect of copper sulfate on S. parasitica metabolism and physiology. The expression profiles of the three GO domains were as follows (each domain showed the top 10 terms).
Molecular function: 1,3-beta-D-glucan synthase activity (4 genes), transferase activity, transferring acyl groups other than amino-acyl groups (16 genes), coenzyme binding (22 genes), lipase activity (9 genes), ammonia-lyase activity (3 genes), amylase activity (3 genes), hydrolase activity, hydrolyzing O-glycosyl compounds (20 genes), hydrolase activity, acting on glycosyl bonds (21 genes), transferase activity, transferring acyl groups (25 genes), cofactor binding (33 genes) (Fig 4).
GO categories (y-axis) were grouped into three main ontologies: biological process, cellular component, and molecular function. The x-axis indicates the statistical significance of the enrichment. All annotated genes were classified into three GO domains: biological processes, cellular component, and molecular function. Red histogram represents the biological processes, the yellow histogram represents cellular component, and the blue histogram represents molecular function. Dissimilar expression profiles were obtained from DEGs in the treated and control samples, revealing the obvious effect of copper sulfate on S. parasitica metabolism and physiology.
Cellular component: organelle envelope lumen (1 gene), membrane part (65 genes), extracellular region (9 genes), endocytic vesicle (2 genes), intrinsic to membrane (62 genes), integral to membrane (62 genes), extracellular region part (4 genes), myosin complex (9 genes), actin cytoskeleton (11 genes), 1,3-beta-D-glucan synthase complex (4 genes) (Fig 4).
Biological processes: single-organism catabolic process (9 genes), small molecule catabolic process (9 genes), small molecule metabolic process (90 genes), organic acid catabolic process (9 genes), carboxylic acid catabolic process (9 genes), cellular amino acid metabolic process (32 genes), single-organism metabolic process (141 genes), organic acid metabolic process (40 genes), oxoacid metabolic process (40 genes), carboxylic acid metabolic process (40 genes) (Fig 4).
Cofactor binding (33 genes), 1,3-beta-D-glucan synthase complex (4 genes), carboxylic acid metabolic process (40 genes) were the most significant GO terms. Globally, the 1515 GO terms were annotated, and the 262 GO terms underwent dramatic expression changes.
Following copper sulfate treatment, genes related to molecular function, including SPRG_10490(ornithine-oxo-acid transaminase), SPRG_03427 (oxoglutarate dehydrogenase), SPRG_06771 (tyrosine aminotransferase), SPRG_15406 (trehalase), and biological process related genes such as SPRG_11730 (histidine ammonia-lyase), SPRG_05010 (tryptophanyl-tRNA synthetase), SPRG_01675 (homogentisate 1%2C2-dioxygenase), SPRG_11441 (uridine phosphorylase), SPRG_11440 (uridine phosphorylase) were repressed, while molecular function related genes like SPRG_07378 (threonine ammonia-lyase) and biological process related genes like SPRG_01336 (maleylacetoacetate isomerase), SPRG_00691 (imidazoleglycerol-phosphate dehydratase) were upregulated. Most of these genes are responsible for protein synthesis, suggesting that copper sulfate mostly affects protein synthesis in S. parasitica. The expression of genes related to cellular flow, signal transduction, and cellular transport was also altered in S. parasitica after treatment with copper sulfate (including downregulated genes such as SPRG_05020, SPRG_13728, SPRG_03925, SPRG_07166, SPRG_17366, SPRG_20259 and SPRG_00668, and upregulated genes, such as SPRG_19069, SPRG_06201 and SPRG_20593), which are closely linked with substance metabolism. Additionally, copper sulfate treatment specifically activated SPRG_08152 (hypothetical protein), SPRG_09968 (hypothetical protein) and SPRG_06695 (hypothetical protein), which only are only expressed in the copper sulfate treated samples. The functions of hypothetical proteins were not clear, but were closely related to the biological functions of organism.
Overall, the above analyses indicate that copper sulfate inhibits S. parasitica growth by affecting multiple biological functions, including protein synthesis and metabolism.
KEGG pathway analysis of differentially expressed genes
To further explore the DEG biological functions, we mapped DEGs to the KEGG database, and then enriched them to important pathways based on the whole transcriptome background. We mapped genes to 203 pathways. Many genes were found in multiple pathways; however, many genes were also restricted to a single pathway. There were 25 significantly enriched pathways, the most significantly enriched of which were metabolic pathways (Fig 5). The metabolism-related biological pathways included metabolic pathways, biosynthesis of secondary metabolites and biosynthesis of amino acids. A total of 98 genes were mapped to the metabolic pathways, 42 genes were mapped to the biosynthesis of secondary metabolites pathways, and 18 genes were mapped to the biosynthesis of amino acids. In addition, lipid metabolism-related biological pathways including fatty acid metabolism (13 genes), arachidonic acid metabolism (9 genes), and glycerolipid metabolism (6 genes), were also significantly enriched. Enriched amino acid metabolism-related biological pathways included phenylalanine metabolism (7 genes), tyrosine metabolism (7 genes), histidine metabolism (7 genes), valine, leucine and isoleucine biosynthesis (5 genes), and phenylalanine, tyrosine and tryptophan biosynthesis (5 genes).
Y-axis, KEGG pathway categories; x-axis, statistical significance of enrichment.
The other significantly enriched pathways included carbohydrate metabolism-related starch and sucrose metabolic pathway (12 genes), metabolism of cofactors and vitamins related ubiquinone and other terpenoid-quinone biosynthesis pathway (4 genes), biosynthesis of other secondary metabolites related phenylpropanoid biosynthesis (4 genes), novobiocin biosynthesis (2 genes), tropane, piperidine and pyridine alkaloid biosynthesis (2 genes), immune diseases related primary immunodeficiency (4 genes), endocrine system related PPAR signaling pathway (8 genes), cyanoamino acid metabolism (3 genes), xenobiotic biodegradation and metabolism related benzoate degradation (3 genes), styrene degradation (2 genes), nitrotoluene degradation (2 genes), membrane transport related ABC transporters (11 genes),digestive system related protein digestion and absorption (5 genes), cancers overview related chemical carcinogenesis (3 genes).
In the metabolism-related biological pathways, the expression of certain genes (such as SPRG_05378, SPRG_01336, SPRG_00691, SPRG_19437, SPRG_17560, SPRG_13611, SPRG_08603, SPRG_09764, and SPRG_01617) were upregulated, while SPRG_11730, SPRG_04491, SPRG_10490, SPRG_12259, SPRG_04186, SPRG_18377, SPRG_18063, SPRG_06771, SPRG_16261, and SPRG_02541 genes were downregulated after copper sulfate treatment. These findings were consistent with the GO enrichment analysis, which indicated that copper sulfate mostly prevented the growth of S. parasitica by affecting protein synthesis and metabolism.
In the lipid metabolism-related biological pathways, among the 13 genes mapped to fatty acid metabolism, 12 genes were downregulated (including SPRG_15243, SPRG_02257, SPRG_06557, and SPRG_01190), but SPRG_09764. In the arachidonic acid metabolism, three genes were upregulated in S. parasitica, (SPRG_12960, SPRG_16543 and SPRG_13907), six genes were downreglated, (SPRG_20516, SPRG_10696, SPRG_09110, SPRG_11344, SPRG_18867 and SPRG_11345) after copper sulfate treatment. These results showed that copper sulfate affected S. parasitica by influencing the energy biogenesis deregulation.
Overall, the results of the DEGs pathway analysis support the viewpoint that copper sulfate inhibits S. parasitica growth by affecting multiple biological functions, such as energy biogenesis, protein synthesis and metabolism.
Differential expression verification of differentially expressed genes
Nine genes with a clearly defined function were randomly selected from the differentially expressed genes identified by RNA-Seq to verify their expression using quantitative PCR. The primer sequences for all of the genes examined are listed in Table 3. The results verified that eight of the genes examined were consistent with the results of RNA-Seq (Fig 6). However, one gene (SPRG_01336), was not consistent with the RNA-Seq data. Overall, these results suggest that the RNA-Seq results are generally reliable; however, further studies are still needed to be carried out to verify the results from this study.
Negative values indicate that the gene expression of S. parasitica was downregulated following copper sulfate treatment; positive values indicate that the gene expression was upregulated.
Discussion
Transcriptome sequencing is a powerful technique for studying the mechanism of changes of biological characteristics of an organism, and has been used successfully in some species [17–19]. In our study, we examined the mechanism by which copper sulfate inhibits the growth of S. parasitica by examining changes in the transcriptome after treatment using the Illumina sequencing platform. The total mapped rates of reads were 90.50% in the S. parasitica transcriptome of control group as compared to the reference genome, which met the required quality of the sequencing data needed for follow-up studies. The total unmapped reads was 9.50% of the total reads in the controls; however, the unmapped reads was 26.50% in the experimental group. It was remarkable that unmapped reads rate of the transcriptome of S. parasitica following treatment with copper sulfate was higher. This discrepancy might be due to the fact that the existing annotation is incomplete, or because copper sulfate caused gene variation.
The transcriptome of S. parasitica treated by copper sulfate contained more partial novel and complete novel splice junctions than that of control group. These findings indicate that copper sulfate may change the alternative splicing of certain genes that affect the metabolism and activity of S. parasitica. Therefore, investigating splice junctions may help understand the mechanism by which copper sulfate inhibits S. parasitica.
In addition to differences in splicing, copper sulfate may also affect S. parasitica by directly altering the expression of genes involved in key cellular processes. In order to investigate this, we classified DGEs into 1515 GO terms, consisting of three domains: biological process, cellular component, and molecular function. Of these, 262 GO terms were found to have dramatic expression changes. We mapped the DGEs to 203 pathways, with 25 of these pathways significantly enriched. The most significantly enriched pathways in this analysis were the metabolism-related biological pathways, including metabolic pathways, biosynthesis of secondary metabolites and biosynthesis of amino acids, which are responsible for the main biological functions of S. parasitica. Our study supplements the previous studies by Jiang et al., who investigated virulence genes in S. parasitica by sequencing its genome [20], and by Torto-Alalibo et al., who investigated expressed sequence tags and disclosed S. parasitica putative virulence factors [21].
In GO annotation enrichment, most of DGEs were related to protein synthesis and metabolism. However, in KEGG pathway analysis, the results showed not only substance metabolism was affected, but the production of energy was also affected. Together, these results indicate that copper sulfate inhibit S. parasitica growth by affecting multiple biological functions, which differ from the mechanisms of other anti–S. parasitica reagents. Such reagents include humic substances with higher molecular weights and aromaticity [7], clotrimazole [8], and saprolmycin A-E [9]. These data indicate that copper sulfate represents a novel anti–S. parasitica treatment.
We found that copper sulfate treatment did not alter the expression of the host-targeting protein 1 (SpHtp1) gene, which is responsible for S. parasitica translocation into host cells, suggesting that copper sulfate may not interfere with S. parasitica–host interactions [22].
Most of the genes in S. parasitica were hypothetical proteins, and the functions were not clear in the existing data. In order to facilitate further research, nine of the DGEs with clear functions were chosen randomly from the RNA-Seq data for qPCR to verify the results of the RNA-Seq. The qPCR results were essentially consistent with the results of the analysis of the transcriptome. One of the genes analyzed had a discrepant result, however, but we still believe the quality of the transcriptome of S. parasitica meet the requirements for studies on functional genes.
Taken together, copper sulfate represents a novel anti–S. parasitica agent that likely functions of inhibiting S. parasitica by affecting energy biogenesis, protein synthesis and metabolism. Our findings provide the basis for further investigation of the potential application of copper sulfate for controlling S. parasitica growth in aquaculture in the future.
Author Contributions
Conceived and designed the experiments: KH RRM. Performed the experiments: JMC NL. Analyzed the data: XY WHF. Contributed reagents/materials/analysis tools: QS. Wrote the paper: HLY HRL. Checked experiments and the manuscript: XLY.
References
- 1. Meyer FP. Aquaculture disease and health management. J Anim Sci.1991; 69: 4201–4208. pmid:1778836
- 2. Rach JJ, Redman S., Bast D., and Gaikowski M. P.. Efficacy of hydrogen peroxide versus formalin treatments to control mortality associated with saprolegniasis on lake trout eggs. North American Journal of Aquaculture.2005; 67: 148–154.
- 3. van West P. Saprolegnia parasitica, an oomycete pathogen with a fishy appetite: new challenges for an old problem. Mycologist.2006; 20: 99–104.
- 4. Alderman DJ. Malachite green: a review. Journal of Fish Diseases.1985; 8:
- 5. Srivastava S, Sinha R, Roy D. Toxicological effects of malachite green. Aquat Toxicol.2004; 66: 319–329. pmid:15129773
- 6. Stammati A, Nebbia C, Angelis ID, Albo AG, Carletti M, Rebecchi C, et al. Effects of malachite green (MG) and its major metabolite, leucomalachite green (LMG), in two human cell lines. Toxicol In Vitro.2005; 19: 853–858. pmid:16061355
- 7. Meinelt T, Paul A, Phan TM, Zwirnmann E, Kruger A, Wienke A, et al. Reduction in vegetative growth of the water mold Saprolegnia parasitica (Coker) by humic substance of different qualities. Aquat Toxicol.2007; 83: 93–103. pmid:17475348
- 8. Warrilow AG, Hull CM, Rolley NJ, Parker JE, Nes WD, Smith SN, et al. Clotrimazole as a potent agent for treating the oomycete fish pathogen Saprolegnia parasitica through inhibition of sterol 14alpha-demethylase (CYP51). Appl Environ Microbiol.2014; 80: 6154–6166. pmid:25085484
- 9. Nakagawa K, Hara C, Tokuyama S, Takada K, Imamura N. Saprolmycins A-E, new angucycline antibiotics active against Saprolegnia parasitica. J Antibiot (Tokyo).2012; 65: 599–607.
- 10. Sun Q, Hu K., Yang X.L.. The efficacy of copper sulfate in controlling infection of Saprolegnia parasitica. Journal of the World Aquaculture Society.2014; 45: 220–225.
- 11. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics.2009; 25: 1105–1111. pmid:19289445
- 12. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc.2012; 7: 562–578. pmid:22383036
- 13. Conesa A, Gotz S. Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics.2008; 2008: 619832. pmid:18483572
- 14. Yan HZ, Liou RF. Selection of internal control genes for real-time quantitative RT-PCR assays in the oomycete plant pathogen Phytophthora parasitica. Fungal Genet Biol.2006; 43: 430–438. pmid:16531084
- 15. van West P, de Bruijn I, Minor KL, Phillips AJ, Robertson EJ, Wawra S, et al. The putative RxLR effector protein SpHtp1 from the fish pathogenic oomycete Saprolegnia parasitica is translocated into fish cells. FEMS Microbiol Lett.2010; 310: 127–137. pmid:20659163
- 16. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods.2001; 25: 402–408. pmid:11846609
- 17. Price DP, Nagarajan V, Churbanov A, Houde P, Milligan B, Drake LL, et al. The fat body transcriptomes of the yellow fever mosquito Aedes aegypti, pre- and post- blood meal. PLoS One.2011; 6: e22573. pmid:21818341
- 18. Peng Y, Gao X, Li R, Cao G. Transcriptome sequencing and De Novo analysis of Youngia japonica using the illumina platform. PLoS One.2014; 9: e90636. pmid:24595283
- 19. Li C, Weng S, Chen Y, Yu X, Lu L, Zhang H, et al. Analysis of Litopenaeus vannamei transcriptome using the next-generation DNA sequencing technique. PLoS One.2012; 7: e47442. pmid:23071809
- 20. Jiang RH, de Bruijn I, Haas BJ, Belmonte R, Lobach L, Christie J, et al. Distinctive expansion of potential virulence genes in the genome of the oomycete fish pathogen Saprolegnia parasitica. PLoS Genet.2013; 9: e1003272. pmid:23785293
- 21. Torto-Alalibo T, Tian M, Gajendran K, Waugh ME, van West P, Kamoun S. Expressed sequence tags from the oomycete fish pathogen Saprolegnia parasitica reveal putative virulence factors. BMC Microbiol.2005; 5: 46. pmid:16076392
- 22. Wawra S, Bain J, Durward E, de Bruijn I, Minor KL, Matena A, et al. Host-targeting protein 1 (SpHtp1) from the oomycete Saprolegnia parasitica translocates specifically into fish cells in a tyrosine-O-sulphate-dependent manner. Proc Natl Acad Sci U S A.2012; 109: 2096–2101. pmid:22308362