- Research article
- Open Access
A gene catalogue for post-diapause development of an anhydrobiotic arthropod Artemia franciscana
© Chen et al; licensee BioMed Central Ltd. 2009
- Received: 28 October 2008
- Accepted: 27 January 2009
- Published: 27 January 2009
Diapause is a reversible state of developmental suspension and found among diverse taxa, from plants to animals, including marsupials and some other mammals. Although previous work has accumulated ample data, the molecular mechanism underlying diapause and reactivation from it remain elusive.
Using Artemia franciscana, a model organism to study the development of post-diapause embryos in Arthropod, we sequenced random clones up to a total of 28,039 ESTs from four cDNA libraries made from dehydrated cysts and three time points after rehydration/reactivation, which were assembled into 8,018 unigene clusters. We identified 324 differentially-expressed genes (DEGs, P < 0.05) based on pairwise comparisons of the four cDNA libraries. We identified a group of genes that are involved in an anti-water-deficit system, including proteases, protease inhibitors, heat shock proteins, and several novel members of the late embryogenesis abundant (LEA) protein family. In addition, we classified most of the up-regulated genes after cyst reactivation into metabolism, biosynthesis, transcription, and translation, and this result is consistent with the rapid development of the embryo. Some of the specific expressions of DEGs were confirmed experimentally based on quantitative real-time PCR.
We found that the first 5-hour period after rehydration is most important for embryonic reactivation of Artemia. As the total number of expressed genes increases significantly, the majority of DEGs were also identified in this period, including a group of water-deficient-induced genes. A group of genes with similar functions have been described in plant seeds; for instance, one of the novel LEA members shares ~70% amino-acid identity with an Arabidopsis EM (embryonic abundant) protein, the closest animal relative to plant LEA families identified thus far. Our findings also suggested that not only nutrition, but also mRNAs are produced and stored during cyst formation to support rapid development after reactivation.
- cDNA Library
- Late Embryogenesis Abundant
- Late Embryogenesis Abundant Protein
- Small Heat Shock Protein
- Late Embryogenesis Abundant Gene
Diapause, also known as discontinuous development, is a reversible state of developmental suspension often promoted by seasonal environmental adversity cues. It is a strategy that ensures successful species survival through timing the post-diapause development for favorable environmental conditions. This protective mechanism is widely distributed among different taxa, including plants, insects, vertebrates, and even some mammals and it occurs at different developmental stages, such as embryonic, larval, pupal, or adult stages, often varying from species to species.
Artemia is a group of small ancient crustaceans living in saline waters. Induced by signals from forthcoming seasonal adversities (not by the adversities themselves), Artemia produces shelled embryos (cysts) that suspend development and metabolism at the gastrula stage, which contain around 4,000 cells [1, 2]. These cysts are in a state of diapause, a physiological dormancy with specific releasing conditions. The primary important feature of diapause is that once initiated, it can only be released by certain stimuli; until then, the cysts wouldn't continue to develop even putting them into favorable environments. This feature is essential in distinguishing diapause as a different phenomenon from other forms of dormancy such as hibernation and quiescence . In Artemia, the diapause could be broken by dehydration (often complete dehydration is required) and is followed by a state of quiescence, from which the cysts can begin direct development, should conditions change to become more favorable. The diapause-broken was often considered as a mechanical process instead of physiological transition, because a) the time of the dehydration process doesn't matter to the hatching rate; b) multiple-round of dehydration- rehydration is sometime necessary to diapause-break and c) little physiological and biological changes were observed . These dehydrated cysts are able to survive years even decades, with little signs of metabolism and energy consumption, but remain viable [5–7]. Once reactivated in favorable environments, they resume development and give rise to free-swimming larva within 24 hours. The developing of embryonic development of Artemia is generally divided into three stages, namely pre-diapause, diapause and post-diapause .
Physiologically and biologically, the encysted embryos resemble plant seeds in many aspects. For example, both contain embryos that develop to certain stages and undergo complete dehydration without any obvious loss of viability even stored for years. Recent studies also revealed that certain molecules such as heat shock proteins are significantly enriched in both cysts and plant seeds.
As an initial step to address the molecular mechanisms of diapause, we first of all constructed a cDNA library using dehydrated cysts that were in a state of quiescence and likely to contain genes important to the maintenance of diapause, since no physiological changes observed during the transition. Then we analyzed its gene expression profiles and compared with that of cDNA libraries that were derived from three time points of subsequent development after rehydration. Expression sequence tags (ESTs) were used to get quick access to gene sequences and expression information of the developing Artemia. We obtained ~30,000 high-quality ESTs from the four cDNA libraries in approximately equal numbers from each, and clustered them into 8,018 unigenes after sequence assembly and annotation. We also validated a few gene expression profiles using quantitative real-time PCR. All EST sequences were deposited into NCBI's dbEST under accessions from ES492186 to ES529129.
Library construction and sequence assembly
Summary of cDNA libraries
Annotation of Artemia unigenes
Expression profiles of Artemia unigenes in four cDNA libraries
We obtained 1,698 unigenes from AfD0 that derived from dehydrated cysts where metabolism and development are nearly suspended. We obtained significantly more unigenes from the other three cDNA libraries: 3,146, 2,787, and 3,006 from AfR5, AfR10, and AfR15, respectively (Table 1). We also carried out pair-wise comparisons for expression profiles among the libraries. For example, we identified 362 unigenes expressed in AfD0 and AfR5, 560 in AfR5 and AfR10, and 558 in AfR10 and AfR15. Only 194 unigenes were found universally expressed and nearly 50% unigenes were found unique to each library. Expression abundances of the unigenes correlate to their distributions in four libraries. For example, universally presented genes were highly expressed, with a mean cluster size of 46.7; the cluster sizes of those unique unigenes strongly biased toward very low values (around one EST per unigene). Genes expressed in two libraries showed in-between abundances with a mean cluster size of 21.2.
Down-regulated and up-regulated DEGs after reactivation
We identified 142 down-regulated genes after reactivation and annotated 77 of them. We grouped cysts-enriched transcripts into several categories based on their putative biological functions, such as stored in cysts and utilized after reactivation or involved in anti-desiccation and diapause. We also look into ribosomal proteins as they are stored in cysts and used after reactivation to facilitate quick physiological and biological transition from the quiescence to the active states , and identified 37 ribosomal protein coding genes as down-regulated DEGs. In addition, we noticed a cathepsin-like cysteine proteinase (CLAP) that is stored in cysts but activated after rehydration and involved in yolk utilization in Artemia embryos . However, other studies shown that proteinases were involved in desiccation-tolerance and the accumulation of proteases was often correlated with increased amount of aggregated/denatured proteins induced by water-loss . As proteases are often multi-functional, our efforts to assign them into function-related groups became difficult sometimes. For example, we identified another down-regulated HAD-like protease (haloacid dehalogenase-like hydrolase) but categorized it as water-loss induced genes.
Water-deficit induces the expression of protease inhibitors and chaperones and serve as counteracting mechanisms to prevent protein degradation. Among down-regulated DEGs, we identified a cysteine protease inhibitor (Cystatin B) that binds tightly to and inhibit papain, cathepsin B and lysosomal cathepsin L  and several chaperone genes, including late embryogenesis abundant proteins (LEAs), a small heat shock protein p26, and HSP70. As classic chaperone, p26 was previous found abundantly expressed in encysted embryos of Artemia and it was suggested to protect proteins from irreversible denaturation in an energy-independent manner  and cells against oxidative damage . In stressed cells, HSP70 and p26 move to the nucleus upon stresses and played a role in stabilization of nuclear matrix proteins .
In a previous study, Qiu Z. et al identified ~85 genes enriched in the diapause-destined embryos of Artemia franciscana by using a subtractive cDNA library . We found only nine down-regulated DEGs that are shred between the two datasets, including those of p26, superoxide dismutase, and several ribosomal proteins (Additional file 2). Since only less than 300 cDNA clones were sequenced in the study and our sampling depth is also rather shallow, we did not anticipate much overlapping data.
Through our annotation pipeline, we were able to annotate 56 of 137 up-regulated DEGs in four cDNA libraries. From them, we identified several functionally related groups, including those involved in (1) metabolism and ATP generation, (2) cell growth and cell division, (3) transcription and translation, and (4) folding of newly translated proteins. The increasing expression of these genes was consistent with the rapid development after reactivation of Artemia cysts, which usually gives rise to free-swimming larva within 24 hours.
Validation of EST expressions by qRT-PCR
qRT-PCR validation of selected DEGs
Cystatin B inhibitor
cathepsin L-like protease
In this report, we described 28,039 Artemia franciscana ESTs and their 8,018 unigenes, representing the largest sequence resource in the public databases for this organism. We confirmed differential expression for p26, a small heat shock protein that was abundantly expressed in encysted embryos and that serves as a multifunctional molecular chaperone. Comparing to a public available dataset that derived from diapause-destined Artemia embryos, we noticed the critical role of p26 in both entering and maintaining dipause. We also identified two novel late embryo abundant (LEA) genes homologous to the plant LEA proteins, and together with HSP70, proteases, and protease inhibitors, these down-regulated genes after cysts reactivation provided the ability of desiccation tolerance for the encysted embryos. A similar set of genes have also been reported in plant seeds .
Among down-regulated DEGs, we identified 37 ribosomal protein coding genes; this is consistent with a previous finding that mRNA activity for ribosomal proteins is stored in the cytoplasm of dormant cysts and associated with polysomes . Other protein coding genes may also be stored in encysted embryos. For example, a cathepsin like cysteine proteinase that was involved in yolk utilization in Artemia embryos  was enriched in dehydrated embryos. Several groups of functionally related genes were up-regulated up to 20 folds after reactivation, which include genes that are involved in energy generation, translation/transcription, metabolism, growth, cell division, and differentiation. The activation of these genes allows Artemia embryos to convert the stored yolk palette into energy and the result is consistent with their rapid development after reactivation.
We obtained dehydrated Artemia franciscana cysts that were in the quiescence state from Dr. Jianhai Xiang's laboratory at the Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China. We confirmed the taxonomic identity based on mitochondrial cox1 sequence from single cysts, a technique known as DNA Bar-coding. We evaluated the hatching rate of the cysts following a previously-described procedure  to optimize hatching rates. The first free-swing larva usually appear at the 10th hour of hatching, 15% of the cysts hatch after 15 hours, and more than 90% of the cysts develop into free-swing larva at the 20th hour. We started our experiments from dehydrated cysts and collected samples from the rehydrated cysts after 5-, 10-, and 15-hour hatching for the cDNA library construction.
cDNA library construction and quality estimation
We extracted total RNAs from the samples using Trizol agent (Invitrogen), and isolated polyA mRNAs using PolyATtract mRNA isolation system (Promega). To obtain a broad coverage of Artemia transcripts, we size-fractioned double-stranded cDNAs before cloning. We constructed cDNA libraries using the directional pBluescript® II XR vector (Stratagene), exploiting the EcoRI and XhoI restriction sites, according to the manufacturer's instruction. The cDNA libraries were not normalized. To assess the quality of cDNA libraries, we performed colony PCR on 96 randomly picked clones to determine the average insert size and percentage of clones without inserts and sequenced 384 randomly picked clones from each cDNA library to determine the ratio of host and vector sequences (empty clones), the sequence length after masking vector sequences, and the ratio of unique sequences (contigs + singlets)/reads.
EST sequencing, assembly, and annotation
We acquired 5' ESTs from ~10,000 colonies of each library, using MegaBase® 1000 sequencers. We assembled the sequences using Phred-phrap-consed [21, 22] with default parameters after removal of vector and low-quality (< 100 bp) sequences. We annotated our unigene sequences (consensus sequences of assembled clusters, including contigs and singlets) based on protein sequences in NCBI non-redundant (nr) protein database and several other public databases including Drosophila melanogaster protein database (dmel-all-translation; release 4.3) from FlyBase  and silkworm Bombyx mori BGF protein database database from SilkDB , using BLAST  based tools, such as blastx and tblastx. We also collected available EST sequences of Insecta and Crustacean from NCBI's dbEST. To classify the Artemia unigenes into Gene Ontology (GO) categories, we compared them with proteins in UniProt (uniprot_sprot and uniprot_trembl) and assigned GO terms according to their best matches, by using a UniProt2GO data provided by the European Bioinformatics Institute (EBI).
Identification of differentially-expressed genes (DEGs) and quantitative PCR validation
We used a program IDEG6 [9, 26] to identify genes that are differentially expressed among libraries. A unigene (or gene) is said to be differentially expressed when it produces a P < 0.05 using Chi-Square algorithm. We further validate the expression profiles of nine DEGs, using quantitative real-time PCR (qRT-PCR). We designed primers for these genes, using a program OLIGO6 with the following parameters: Tm, 60 ± 2°C, difference < = 2 between a pair of primers; primer length, 17–22 bp; GC content, 40–60%; PCR product length, 150–200 bp. Whenever possible, we chose primers that locate at the 3' end of transcripts. A full list of designed primers is shown in Table 2 and their sequences in Additional file 3.
We used the same RNA samples as what for constructing the cDNA libraries. The RNA samples were treated with RNase-free DNase I (Promega) to remove possible DNA contaminations. The first-strand cDNA was synthesized by using 500 ng total RNAs, poly(T) primers and SSII reverse polymerase (Invitrogen). qRT-PCR was conducted by using a Quant SYBR Green PCR kit (Tiangen, China). We chose the GADPH gene as an external reference for data normalization. The PCR reaction parameters were as follows: 95°C for 2 min; 40 cycles of 3-temperature of 95°C for 15 s, 60°C for 20 s, 72°C for 30 s. Three replicates for each pair of primers per template were included. qRT-PCR data were analyzed by using Opticon Monitor ® software. Melting curves for each PCR were carefully analyzed to avoid non-specific amplifications. Gene expressions were quantified and transformed by using the ΔCt formula normalized with the expression of GAPDH.
We would like to thank Dr. Jianhai Xiang of Institute of Oceanology, Chinese Academy of Sciences, for kindly providing Artemia cysts. The work is supported by grants from Chinese Academy of Sciences (KSCX2-SW-223 to Jun Yu), Chinese Natural Science Foundation (30270748 to Jun Yu), and Ministry of Sciences and Technology (2005AA235110 to Jun Yu).
- Anderson E, Lochhead JH, Lochhead MS, Huebner E: The origin and structure of the tertiary envelope in thick-shelled eggs of the brine shrimp, Artemia. J Ultrastr Res. 1970, 32 (5–6): 497-525. 10.1016/S0022-5320(70)80025-9.View ArticleGoogle Scholar
- Morris JE, Afzelius BA: The structure of the shell and outer membranes in encysted Artemianext term salina embryos during cryptobiosis and development. J Ultrastr Res. 1967, 20 (3–4): 244-259. 10.1016/S0022-5320(67)90285-7.View ArticleGoogle Scholar
- Taylor F: Insect Life Histories: Seasonal Adaptations of Insects. Science. 1986, 232 (4754): 1152-10.1126/science.232.4754.1152-a.View ArticlePubMedGoogle Scholar
- Abatzopoulos TJB, J A, Clegg JS, Sorgeloos P: 2002, Artemia: Basic and Applied Biology: Springer, [http://www.springer.com/life+sci/zoology/book/978-1-4020-0746-0]
- Warner AH, Clegg JS: Diguanosine nucleotide metabolism and the survival of artemia embryos during years of continuous anoxia. Eur J Biochem. 2001, 268 (6): 1568-1576. 10.1046/j.1432-1327.2001.01993.x.View ArticlePubMedGoogle Scholar
- Clegg J: Embryos of Artemia franciscana survive four years of continuous anoxia: the case for complete metabolic rate depression. J Exp Biol. 1997, 200 (Pt 3): 467-475.PubMedGoogle Scholar
- Hontoria F, Crowe JH, Crowe LM, Amat F: METABOLIC HEAT PRODUCTION BY ARTEMIA EMBRYOS UNDER ANOXIC CONDITIONS. J Exp Biol. 1993, 178 (1): 149-159. [http://jeb.biologists.org/cgi/content/abstract/178/1/149]Google Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16 (6): 276-277. 10.1016/S0168-9525(00)02024-2.View ArticlePubMedGoogle Scholar
- Romualdi C, Bortoluzzi S, Danieli GA: Detecting differentially expressed genes in multiple tag sampling experiments: comparative evaluation of statistical tests. Hum Mol Genet. 2001, 10 (19): 2133-2141. 10.1093/hmg/10.19.2133.View ArticlePubMedGoogle Scholar
- Pierandrei-Amaldi P, Campioni N: Messenger RNA for ribosomal proteins in dormant and developing Artemia salina embryos. Biochim Biophys Acta. 1981, 655 (3): 359-365.View ArticlePubMedGoogle Scholar
- Warner AH, Perz MJ, Osahan JK, Zielinski BS: Potential role in development of the major cysteine protease in larvae of the brine shrimp Artemia franciscana. Cell Tissue Res. 1995, 282 (1): 21-31.PubMedGoogle Scholar
- Bray EA: Molecular Responses to Water Deficit. Plant Physiol. 1993, 103 (4): 1035-1040.PubMed CentralPubMedGoogle Scholar
- Abrahamson M, Alvarez-Fernandez M, Nathanson CM: Cystatins. Biochem Soc Symp. 2003, 179-199. 70Google Scholar
- Day RM, Gupta JS, MacRae TH: A small heat shock/alpha-crystallin protein from encysted Artemia embryos suppresses tubulin denaturation. Cell Stress Chaperones. 2003, 8 (2): 183-193. 10.1379/1466-1268(2003)008<0183:ASHCPF>2.0.CO;2.PubMed CentralView ArticlePubMedGoogle Scholar
- Collins CH, Clegg JS: A small heat-shock protein, p26, from the crustacean Artemia protects mammalian cells (Cos-1) against oxidative damage. Cell Biol Int. 2004, 28 (6): 449-455. 10.1016/j.cellbi.2004.03.014.View ArticlePubMedGoogle Scholar
- Willsie JK, Clegg JS: Small heat shock protein p26 associates with nuclear lamins and HSP70 in nuclei and nuclear matrix fractions from stressed cells. J Cell Biochem. 2002, 84 (3): 601-614. 10.1002/jcb.10040.View ArticlePubMedGoogle Scholar
- Ingram J, Bartels D: The Molecular Basis Of Dehydration Tolerance In Plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 377-403. 10.1146/annurev.arplant.47.1.377.View ArticlePubMedGoogle Scholar
- Hand Steven, J D, M M, W T: Life without water: expression of plant LEA genes by an anhydrobiotic arthropod. Journal of Experimental Zoology Part A: Ecological Genetics and Physiology. 2007, 307A (1): 62-66. 10.1002/jez.a.343.View ArticleGoogle Scholar
- Qiu Z, Tsoi SC, MacRae TH: Gene expression in diapause-destined embryos of the crustacean, Artemia franciscana. Mech Dev. 2007, 124 (11–12): 856-867. 10.1016/j.mod.2007.09.001.View ArticlePubMedGoogle Scholar
- Wang W, Meng B, Chen W, Ge X, Liu S, Yu J: A proteomic study on postdiapaused embryonic development of brine shrimp (Artemia franciscana). Proteomics. 2007, 7 (19): 3580-3591. 10.1002/pmic.200700259.View ArticlePubMedGoogle Scholar
- Ewing B, Green P: Base-Calling of Automated Sequencer Traces Using Phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.View ArticlePubMedGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-Calling of Automated Sequencer Traces Using Phred. I. Accuracy assessment. Genome Res. 1998, 8 (3): 175-185.View ArticlePubMedGoogle Scholar
- Drysdale RA, Crosby MA: FlyBase: genes and gene models. Nucleic Acids Res. 2005, D390-395. 33 DatabaseGoogle Scholar
- Wang J, Xia Q, He X, Dai M, Ruan J, Chen J, Yu G, Yuan H, Hu Y, Li R: SilkDB: a knowledgebase for silkworm biology and genomics. Nucleic Acids Res. 2005, D399-402. 33 DatabaseGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Romualdi C, Bortoluzzi S, d'Alessi F, Danieli GA: IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol Genomics. 2003, 12 (2): 159-162.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.