E cloned using pGEM-T easy vector system kit (Promega Corporation, Madison, WI, USA), transformed into chemically competent JM109 Escherichia coli (Promega Corporation), and plated onto Luria-Bertani (LB) agar with ampicillin, 5-bromo-4-chloro-3-indolyl–D-galactopyranoside (X-gal) and isopropyl -D-thiogalactopyranoside (IPTG). Selected white colonies were grown overnight in LB broth with ampicillin. The plasmids were extracted using the resin-based plasmid miniprep kit (Axygen Stattic web Biosciences, Union City, CA, USA). The plasmids were quantified by the NanoDrop ND-1000 spectrophotometer. Approximately 80?00 ng of plasmid DNA was used in BigDye Terminator v3.1 Cycle Sequencing Kit (Life Technologies Corporation) with 2 M T7 primers. Excess fluorescent nucleotides and salts were removed from the samples by ethanol precipitation. The dried samples were resuspended in Hi-Di Formamide (Life Technologies Corporation) before loading to the Prism 3130XL sequencer (Life Technologies Corporation). A total of 500 clones for each forward and reverse library were selected for sequencing. Sequence output was exported as text and edited manually to remove vector sequences using BioEdit Sequence Alignment Editor software version 7.0.9 [15]. BLAST searches were performed using the tBLASTx algorithm [16] and default search conditions. Proteins were considered significant when the E value was <1E-04. The annotated sequences were grouped based on Gene Ontology classification. The sequences were deposited in Genbank EST database and were assigned with accession numbers JZ575382 to JZ575617.Relative quantitative real-time PCR (qPCR)In order to validate the changes obtained in the SSH studies, nine genes were selected for the determination of mRNA expression using quantitative real-time PCR (qPCR). These include acyl-CoA desaturase (acd), argininosuccinate synthetase 1 (ass1), betaine-homocysteine Smethyltransferase 1 (bhmt1), ceruloplasmin (cp), carbamoyl phosphate synthetase III (cpsIII), fumarate hydratase (fh), ferritin light chain (ftl), glyceraldehyde-3-phosphate dehydrogenase (gapdh) and superoxide dismutase 1 (sod1). Prior to first strand cDNA synthesis, RNA from the liver of fish kept in fresh water, aestivated for 6 months in air or aroused for 1 day after 6 months of aestivation in air were treated separately with Deoxyribonuclease I (Qiagen Inc.) to remove any contaminating genomic DNA. First strand cDNA was synthesized from 1 g of total RNA using random hexamer jir.2013.0113 primer and the RevertAid first stand cDNA synthesis kit, following the manufacturer’s instruction (Thermo Fisher SIS3 custom synthesis Scientific Inc). mRNA expressions of the selected genes were quantified using a StepOnePlus Real-Time PCR System (Life Technologies Corporation). Each PCR reaction contained 5 l of 2x Fast SYBR Green Master Mix (Life Technologies Corporation), a certain aliquot of gene-specific primers (listed in Table 1) and 0.1? ng of cDNA in a total volume of 10 jir.2010.0097 l. Samples were run in triplicate. qPCR reactions were performed with the following cycling conditions: 95 for 20 s (1 cycle), followed by 40 cycles of 95 for 3 s and 60 of 30 s. Data was collected at each elongation step. Each run was followed by a melt curve analysis by increasing the temperature from 60 to 95 at 0.3 increment to confirm the presence of only a single PCR product. In addition, random PCR products were electrophoresed in a 1.8 agarose gel to verify that only one band was present. AllPLOS ONE | DOI:10.1371/journal.po.E cloned using pGEM-T easy vector system kit (Promega Corporation, Madison, WI, USA), transformed into chemically competent JM109 Escherichia coli (Promega Corporation), and plated onto Luria-Bertani (LB) agar with ampicillin, 5-bromo-4-chloro-3-indolyl–D-galactopyranoside (X-gal) and isopropyl -D-thiogalactopyranoside (IPTG). Selected white colonies were grown overnight in LB broth with ampicillin. The plasmids were extracted using the resin-based plasmid miniprep kit (Axygen Biosciences, Union City, CA, USA). The plasmids were quantified by the NanoDrop ND-1000 spectrophotometer. Approximately 80?00 ng of plasmid DNA was used in BigDye Terminator v3.1 Cycle Sequencing Kit (Life Technologies Corporation) with 2 M T7 primers. Excess fluorescent nucleotides and salts were removed from the samples by ethanol precipitation. The dried samples were resuspended in Hi-Di Formamide (Life Technologies Corporation) before loading to the Prism 3130XL sequencer (Life Technologies Corporation). A total of 500 clones for each forward and reverse library were selected for sequencing. Sequence output was exported as text and edited manually to remove vector sequences using BioEdit Sequence Alignment Editor software version 7.0.9 [15]. BLAST searches were performed using the tBLASTx algorithm [16] and default search conditions. Proteins were considered significant when the E value was <1E-04. The annotated sequences were grouped based on Gene Ontology classification. The sequences were deposited in Genbank EST database and were assigned with accession numbers JZ575382 to JZ575617.Relative quantitative real-time PCR (qPCR)In order to validate the changes obtained in the SSH studies, nine genes were selected for the determination of mRNA expression using quantitative real-time PCR (qPCR). These include acyl-CoA desaturase (acd), argininosuccinate synthetase 1 (ass1), betaine-homocysteine Smethyltransferase 1 (bhmt1), ceruloplasmin (cp), carbamoyl phosphate synthetase III (cpsIII), fumarate hydratase (fh), ferritin light chain (ftl), glyceraldehyde-3-phosphate dehydrogenase (gapdh) and superoxide dismutase 1 (sod1). Prior to first strand cDNA synthesis, RNA from the liver of fish kept in fresh water, aestivated for 6 months in air or aroused for 1 day after 6 months of aestivation in air were treated separately with Deoxyribonuclease I (Qiagen Inc.) to remove any contaminating genomic DNA. First strand cDNA was synthesized from 1 g of total RNA using random hexamer jir.2013.0113 primer and the RevertAid first stand cDNA synthesis kit, following the manufacturer’s instruction (Thermo Fisher Scientific Inc). mRNA expressions of the selected genes were quantified using a StepOnePlus Real-Time PCR System (Life Technologies Corporation). Each PCR reaction contained 5 l of 2x Fast SYBR Green Master Mix (Life Technologies Corporation), a certain aliquot of gene-specific primers (listed in Table 1) and 0.1? ng of cDNA in a total volume of 10 jir.2010.0097 l. Samples were run in triplicate. qPCR reactions were performed with the following cycling conditions: 95 for 20 s (1 cycle), followed by 40 cycles of 95 for 3 s and 60 of 30 s. Data was collected at each elongation step. Each run was followed by a melt curve analysis by increasing the temperature from 60 to 95 at 0.3 increment to confirm the presence of only a single PCR product. In addition, random PCR products were electrophoresed in a 1.8 agarose gel to verify that only one band was present. AllPLOS ONE | DOI:10.1371/journal.po.