Construction of a cDNA library and preliminary analysis of the expressed sequence tags of the earthworm Eisenia fetida (Savigny, 1826)
- Authors:
- Published online on: February 5, 2019 https://doi.org/10.3892/mmr.2019.9938
- Pages: 2537-2550
-
Copyright: © Liu et al. This is an open access article distributed under the terms of Creative Commons Attribution License.
Abstract
Introduction
Earthworms are terrestrial annelids in the oligochaeta subclass, with a generally preferred habitat of damp and loose soil. They include ~3,000 species worldwide, with 229 species in China (1,2). In a wide variety of soil types, earthworms serve vital roles in converting large pieces of organic matter into rich humus to enrich soil fertility. The earthworms are the highest evolutionary species capable of regenerating an anterior portion containing the central nerve system, heart and clitellum (3). The anterior regeneration is a unique developmental process that requires cell proliferation, re-differentiation and sophisticated cell-cell communication. This process can serve as a useful model for investigating normal development and differentiation (4).
Over the past several years, cDNA library construction and analysis have become established as indispensable methods for functional genome analysis since they provide detailed information about the genomic mechanisms underlying the diverse processes of an organism (5). However, conventionally generated cDNA libraries contain a high percentage of 5′-truncated clones, limiting the utility of such libraries. The Switching Mechanism at 5′End of RNA Template (SMART) technique (6) amplifies and enriches the full-length mRNA, and thus generates cDNA libraries with a significantly improved ratio of full-length to partial cDNA sequences. In the present study, the SMART technique was adopted to construct a high quality library of full-length cDNAs representative of adult earthworms, namely of the earthworm Eisenia fetida (Savigny, 1826).
Unlike other model organisms, none of the oligochaete genomes have been sequenced to the best of our knowledge, and genomic research on earthworms lags behind that of other model species such as Mus musculus. In the absence of the full genome sequences, expressed sequence tags (ESTs) aid the rapid detection of expressed genes via sequence analysis, and are a significant resource for comparative and functional genomic studies (2). In addition, among the biological techniques for transcriptome analysis, the determination of ESTs is considered the simplest method for profiling the transcriptome, which is also particularly useful in the development of cDNA microarrays for systematic identification of differentially expressed genes (7). Analysis of ESTs is an effective method for rapidly analyzing gene expression, characterizing gene functions and discovering new genes that are important for specific developmental and physiological processes (8). The present study established 593 ESTs, representing 168 genes and 425 unknown tags, providing a gene expression profile of earthworm development. This collection of ESTs may provide a valuable basis for future research on the physiology of earthworms.
Materials and methods
Isolation of total RNA and mRNA
Eisenia fetida earthworms were purchased from Beijing Shuangqiao Farm (Beijing, China). Fully developed adult Eisenia fetida earthworms weighing 0.3–0.6 g (live weight) were selected for all experiments. All earthworm tissues were harvested and total RNA was isolated using TRIzol reagent (Invitrogen; Thermo Fisher Scientific, Inc., Waltham, MA, USA). The integrity of total RNA was analyzed by electrophoresis using 1% agarose gels. Isolation of poly(A) mRNA from total RNA was carried out using an MN-NucleoTrap®mRNA kit (Machery-Nagel GmbH & Co. KG, Düren, Germany) according to the manufacturer's protocol. Oligo(dT) beads suspension was applied to total RNA and incubated at 68°C for 5 min before eluting mRNA. Then ethidium bromide (EB) staining was applied and 1% agarose gel was used to visualize the result. The isolated mRNA was further vacuum concentrated using Concentrator plus™ (Eppendorf, Hamburg, Germany). The quantity and integrity of isolated mRNA were determined using a nanodrop spectrophotometer and agarose gel electrophoresis, respectively.
cDNA library construction
A total of ~8,048.4 ng mRNA was used for single-stranded cDNA synthesis. The purified mRNA was used as a template, Oligod(T)18 with XhoI cleavage site was used as the primer, and first strand cDNA was transcribed at 42°C using SuperScript™ II RnaseH-Reverse (Thermo Fisher Scientific, Inc.). Then the mRNA was digested using RNaseH, and the resultant mRNA fragments were used as further primers. The first cDNA chain was used as a template for double-stranded cDNA synthesis, using DNA Polymerase I (Takara Biotechnology Co., Ltd., Beijing, China). The ends of the double-stranded cDNA were ligated by T4 DNA polymerase and the ligation products were purified by phenol/chloroform/isoamyl alcohol to remove excess impurities such as protein. Subsequently, the double-stranded DNA fragments were ligated into EcoRI Adaptor using T4 DNA ligase at 4°C overnight. Then the double-stranded DNA fragments were phosphorylated with T4 Polynucleotide Kinase and digested with XhoI. Following XhoI digestion of the double-stranded cDNA, producing XhoI sticky ends, a QIAEXII Gel Extraction kit (Beijing BioDev-Tech, Beijing, China) was used to recycle 0.5–4 Kb fragments. The recycled cDNA was preserved at −20°C. Then the cDNA was ligated into the pBluescript II SK(+) XR vector (Promega Corporation, Madison, WI, USA) in a 3:1 molar ratio with T4 DNA ligase at 4°C overnight. To reduce the redundancy and avoid the underrepresentation of different transcript species, cDNA fragments with different fractionated sizes were balanced and subjected to library construction (9,10). Prior to transformation, mixing of all ligated products with microporous membranes was performed to remove salt ions. Subsequently the products were transformed into 5×107/ml DH10B competent cells (Thermo Fisher Scientific, Inc.), plated on agar plates (10 cm diameter) by pipetting the cells onto the middle of the plate and spreading, and monoclonal colonies were selected for PCR amplification. The inserted sequences in the plasmids were amplified by PCR using T3 primers (5′-ATTAACCCTCACTAAAGGGA-3′) and T7 primers (5′-TAATACGACTCACTATAGGG-3′). The total volume of PCR reaction mixture was 20 µl, containing 1 µl template, 10 µl 2XTaq MasterMix (CWBIO, Beijing, China), 1 µl T3 primers (10 pmol), 1 µl T7 primers (10 pmol) and 7 µl ddH2O (CWBIO). Cycling conditions were: 94°C for 5 min, followed by 30 cycles of 94°C for 30 sec, 55°C for 40 sec and 72°C for 60 sec, followed by 72°C for 5 min.
Bioinformatic analysis
cDNA clones were selected randomly from the cDNA library and the vector sequences were trimmed from the raw sequence data using Vecscreen tool (www.ncbi.nlm.nih.gov/tools/vecscreen/) from the National Center for Biotechnology Information. The sequence of each EST was also edited, mainly to remove ambiguous bases and poor-quality sequences (nucleotide sequences <100 bp). All edited sequences were assembled into groups using SeqMan software version 8.0 (DNASTAR, Madison, WI, USA). The processed cDNA sequences were used to perform a BLAST search in the GenBank database to compare all available ESTs and genes to date (11). The Basic Local Alignment Search Tool (BLAST; blast.ncbi.nlm.nih.gov/Blast.cgi) results with P-values <1×10−8 were generally regarded as a significant match (12,13). A large-scale Unigene assembly of the ESTs was initiated to identify and functionally annotate as many unique transcripts as possible. BLAST analysis against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and protein and nucleic acid databases was conducted for examination of biological functions. The ESTs homologous to known proteins were further annotated for Gene Ontology (GO; www.geneontology.org) terms and the GO analysis was carried out using WEB-based GEne SeT AnaLysis Toolkit (WebGestalt; www.webgestalt.org/option.php) (14,15).
Results
Construction of cDNA library
Obtaining an adequate quantity of high quality mRNA initially is the key to yielding a sufficient quantity of first-strand full-length cDNA by reverse transcription. In the present study, total RNA was extracted from the tissues of earthworms. As shown in Fig. 1, 28s and 18s bands were clearly visible in the electrophoresis gel of total RNA, indicating that the total RNA was obtained. The optical density (OD)260/OD280 ratio for the total RNA was 2.04, well within the range of 1.8–2.1, indicating that the isolated total RNA was suitable for cDNA library construction.
Once the double-stranded cDNA was synthesized as described in the Materials and methods, the present study determined the size distribution of the products. Diffuse strips between 0.5–4.0 kb could be detected by 1% agarose gel electrophoresis, which indicated that double-stranded cDNAs were successfully synthesized. A cDNA library of 4.12×105 clones was obtained and half of the bacteria were cultured for amplification, which produced a total of 1.4×1011 clones. Several colonies were selected, and the inserted sequences in the plasmids were amplified by PCR using T3 and T7 primers. The PCR products were detected using 1% agarose gel electrophoresis as clear bands. No nonspecific bands were identified and the recombinant rate was 97% (Fig. 2).
To investigate the quality of the full-length cDNA library, the lengths of the cDNA inserts were assessed. Sequence outputs were manually edited to remove vector and ambiguous sequences. Then, the sequence data of the cDNA clones obtained by random partial sequencing were searched in the NCBI GenBank using BLAST to identify similarities with sequences in the nucleic acid databases. An evaluation of cDNA insert size and its distribution revealed a low level of insert size bias in the final cDNA library. The majority of the cDNA inserts were larger than 500 bp.
EST analysis
Instead of the amplified library, the primary cDNA library was used to generate ESTs to reduce the redundancy of cDNA clones. Following the removal of the redundant sequences and low-quality sequences (<100 bp), 1,148 effective sequences (>100 bp) from the total cDNA sequences were obtained. As shown in Fig. 3, 53 ESTs were 100–300 bp, 261 ESTs were 300–500 bp, 828 ESTs were 500–700 bp and 6 ESTs were larger than 700 bp. Taken together, 1,148 ESTs were larger than 100 bp. Among them, the shortest sequence was 100 bp, the longest was 718 bp and the average length was ~452 bp. Following sequencing, a homology BLAST search and assembling of the data, 368 singletons and 225 contigs were obtained out of the 1,140 high-quality ESTs, as shown in Table I. Additionally, a total of 593 individual ESTs were analyzed and 168 ESTs annotated in GenBank with nematode homology (Table I).
GO annotation and bioinformatic analysis
The cDNA functions were classified using the GO database into the three main categories of molecular functions, cell components and biological processes.
GO annotation of genes associated with molecular functions indicated that among the 168 ESTs, 46% (77/168 ESTs) were associated with growth and metabolic pathways, with the distribution of the 77 ESTs shown in Fig. 4A. Out of the 77 ESTs, 21% (16 ESTs) were associated with ‘proteolytic enzymes’, 16% (12 ESTs) with ‘protein ligases’, 14% (11 ESTs) with ‘oxido-reductases’, 13% (10 ESTs) with ‘energy release’, 10% (8 ESTs) with ‘signal transduction and cell communication’, 5% (4 ESTs) with ‘transport’ and only 3% (2 ESTs) with ‘post-translational modification’, ‘protein turnover’ and ‘chaperones’ (Fig. 4B).
Cellular components associated with the cDNAs included ‘myosin’, the ‘citrate lyase compound’, the ‘mitochondrial inner membrane translocase compound’, ‘microtubules’, the ‘mitochondrial inner membrane’ and ‘nucleosomes’. The proportions of cellular components are presented in Fig. 5. It can be observed that the proportion of ‘myosin’ among cellular components was the largest. Myosin, an actin-dependent molecular motor, is involved in a number of important functions in earthworms. In particular, the myosin network can drive movement and support different moving speeds of earthworms. This specific feature is closely related to the free moving ability of earthworms (16).
Regarding biological processes, known genes were determined as those presenting significant matches to protein sequences with known functions in non-redundant nucleotide databases. According to these biological functions, the biological processes component was divided into different functions including ‘larval development’ (46%), ‘changes of cell morphogenesis’ (6%), ‘the process of cytokinesis’ (12%), ‘post-translational protein modification’ (15%), ‘stress response’ (5%), ‘cell redox homeostasis’ (2%), ‘protein polymerization’ (2%) and ‘protein catabolism’ (12%), as shown in Fig. 6. It was concluded from the above data that promoting growth was considered to be an important biological function of genes associated with biological processes.
KEGG pathway annotation
KEGG is a collection of online databases describing pathways associated with biochemical, genomic and enzymatic processes. Furthermore, it provides annotations of biochemical pathways for the species in which the genome has been sequenced (17). In this analysis, proteins are not viewed as individual gene products but are organized into pathways and networks according to their biological function(s). In the present study, from the data in Table II, 15 of the 168 ESTs were revealed to be involved in metabolism. Notably, 4 of the 15 ESTs (27%) were involved in the ubiquitin-mediated proteolysis pathway, which was the most represented. The glutathione metabolism (3 ESTs) and arachidonic acid metabolism (2 ESTs) pathways were the second and third most represented pathways, respectively. Additionally, chondroitin sulfate biosynthesis, heparan sulfate biosynthesis, riboflavin metabolism, selenoamino acid metabolism, γ-hexachlorocyclohexane degradation, and the fructose and mannose metabolism pathways were also represented.
Mutual relationship between 287 genes
Relationships among the 168 ESTs were analyzed using the STRING database, and a functional association network was determined with 287 nodes, as shown in Table III. Each node corresponds to a gene and each (weighted) edge represents the evidence of a functional association between the gene pair. Predicted potential regulators are presented in Fig. 7. It can be seen from the STRING results that the most associated nodes (blue) included the 19 genes: Ribosomal protein L (rpl)-1, ribosomal protein S (rps)-0, rpl-4, rpl-5, rps-13, rps-2, acidic ribosomal protein (rla)-1, translocon-associated protein-4, transcription factor BTF3 homolog, rpl-7A, iff-1, rps-17, elongation factor 1α3, polyadenylate-binding protein-1, rps-28, translationally-controlled tumor protein homolog-1, rpl-18, rps-4 and adenosine triphosphate synthase subunit β mitochondrial precursor-2. Of note, the genes rpl-4, rps-13, rpl-5, rpl-1, rps-0, rla-1 and rpl-7A served critical roles in overall gene expression networks.
Discussion
The present study reported on the construction of a high-quality cDNA library from tissues of the earthworm Eisenia fetida (Savigny, 1826), following preliminary analysis of ESTs, putative functional analysis of the ESTs and the gene expression pattern associated with the physiological functions of this organism. cDNA libraries are widely used to identify genes and splice variants and are considered to be a physical resource for the construction of full-length clones (18,19). In the present study, a cDNA library was utilized to provide a molecular resource for the analysis of genes involved in the specific biology of earthworms in terms of their development, survival, pathogenicity and virulence. There are two main factors to consider when assessing the quality of a cDNA library: Representation and cDNA lengths. According to Clarke-Carbon's formula (20), a cDNA library should contain at least 1.7×105 independent clones to ensure that 99% of low-abundance mRNAs will be represented in the library (21). Furthermore, the average length of the inserted cDNAs should be no less than 1.0 kb to ensure the integrity of cDNAs, indicating that in the present study the fragment sizes were effective for ensuring full-length cDNAs in the cDNA library. Since selection bias could favor the smaller cDNAs, the present study used fewer PCR cycles to minimize such bias as previously suggested (10). In addition, up to 25 PCR amplification cycles were used to generate an adequate amount of cDNA for cloning.
The generation of ESTs is an effective and unique approach in molecular studies as it allows for the analysis and measurement of gene expression, as well as simultaneous discovery of new genes. As each EST represents a copy of the functional part of a genome, the study of ESTs is believed to be a more effective way to discover functional genes (22). Furthermore, analysis of the expression of a large number of genes combined with the knowledge of their functions enables insight into the overall situation in terms of biological processes (23), for the current purposes in earthworms. In the present study, ~91% of the ESTs generated were sequences with known or putative functions, while the remainder represented unknown proteins or sequences with no similarities to those in the databases. Although close to 600 ESTs were reported, this is actually far from what could be considered as a ‘complete’ transcriptome (which usually includes between 15,000–20,000 to >100,000 ESTs). Therefore, the present characterization of this seemingly partial transcriptome may far from reflect the full transcriptomic profile of tissues in earthworms.
A comparison of the classification of ESTs with a C. elegans cDNA library based on their putative functions was conducted. Based on identification of clusters via GO analysis, 168 ESTs were matched to C. elegans genes by BLASTx. It is well known that earthworms serve significant roles in organic matter decomposition and mineral cycling, and thus are considered to be important contributors to soil fertility and humification processes (24). In the present study, hydrolytic enzyme activity, conjugating protein ligase activity, oxidation reduction and energy release activity of metabolic enzymes accounted for a large proportion of the molecular functions component. These molecular functions are considered to be a key part of the physiological functions of earthworms, which allow earthworms to survive in different soil environments.
Cell components, as part of GO annotation, are mainly categorized based on subcellular location (including cell cytoplasm, mitochondria, lysosomes, nucleus, microtubules, plasma membrane and myosin), which is highly important for the study of protein functions (25). The results of subcellular localization analysis can provide significant clues to aid the understanding of protein functions. In the present study, several cell components were determined through GO analysis of the 168 ESTs (annotated genes). It was evident that these cell components had strong associations with the regulation of gene expression during the biological development of earthworms, enabling the regeneration of the anterior portion, alterations in movement ability and tissue differentiation.
Due to the temporal specificity of gene expression and interactions with other gene products, the specific pathway undertaken, sequence of gene expression and expression pattern may ultimately change the effects of multiple pathways in earthworms (2). Therefore, the cDNA library and transcription profile of genes representative of fully developed adult earthworms may differ markedly to those of juvenile earthworms. In the present study, gene expression profiles representative of adult earthworm development were generated; however, gene expression profiling of juvenile earthworms was not performed, nor were analyses of the expected differential gene expression between juveniles and adults. In the present study, 168 individual ESTs of earthworms were analyzed in terms of KEGG pathway annotation, which identified 9 corresponding categories. Among them, glutathione metabolism is involved in antioxidant defense systems in Eisenia Andrei, and the associated enzymes are mainly identified in cytosolic fractions (26). Chondroitin sulfate and heparan sulfate biosynthesis are involved in biosynthesis pathways, which have important effects on growth and regeneration. In addition, heparan sulfate is also considered to be a type of anticoagulant, which may also be a function of arachidonic acid (27). Therefore, understanding the molecular function of earthworms may provide some basis for the treatment of thrombotic disorders.
Information on functional annotation and relevant biological interactions associated with a particular gene is available from many online resources. The gene network comprises a collection of genes that cooperate with each other to control the main biological processes. The STRING database suggested a functional context for earthworm lumbrokinase with unknown specificity. A previous study revealed that lumbrokinase and dilong administration can efficiently reduce the incidence of cardiac disease among nonsmokers exposed to second-hand smoke (28). In this regard, the discovery of genes and protein interactions in earthworms has provided a basis for further investigation into human diseases.
Acknowledgements
Not applicable.
Funding
The present study was supported by the National Natural Science Foundation of China (grant nos. 30472251, 30901821 and 81370312), the Natural Science Foundation for Young Scientists of Shanxi Province, China (grant no. 2010021035-2), Natural Science Foundation of Shanxi Province, China (grant no. 2015011113), International Scientific and Technological Cooperative Foundation of Shanxi Province, China (grant no. 201703D421023), the Fund for Shanxi ‘1331 Project’ Key Subjects Construction and the Fund for Shanxi Key Subjects Construction.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Authors' contributions
BFY, BN and JX conceived and designed the experiments and bioinformatics analysis. CL and FXM performed bioinformatics analysis, and wrote the manuscript. XW, PYM, QZ and JBT performed the experiments. RG, ZZL, HLW and NLC contributed to designing the present study and revising the manuscript. JHW and GQS analyzed and interpreted the data.
Ethics approval and consent to participate
Not applicable.
Patient consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Glossary
Abbreviations
Abbreviations:
EST |
expressed sequence tags |
GO |
Gene Ontology |
KEGG |
Kyoto Encyclopedia of Genes and Genomes |
References
Torres-Leguizamon M, Mathieu J, Decaens T and Dupont L: Genetic structure of earthworm populations at a regional scale: Inferences from mitochondrial and microsatellite molecular markers in Aporrectodea icterica (Savigny 1826). PLoS One. 9:e1015972014. View Article : Google Scholar : PubMed/NCBI | |
Pirooznia M, Gong P, Guan X, Inouye LS, Yang K, Perkins EJ and Deng Y: Cloning, analysis and functional annotation of expressed sequence tags from the Earthworm Eisenia fetida. BMC Bioinformatics. 8 (Suppl 7):S72007. View Article : Google Scholar : PubMed/NCBI | |
Cho SJ, Koh KS, Lee E and Park SC: Differential expression of threelabial genes during earthworm head regeneration. Biosci Biotechnol Biochem. 73:2609–2614. 2014. View Article : Google Scholar | |
Zheng P, Shao Q, Diao X, Li Z and Han Q: Expression of stem cell pluripotency factors during regeneration in the earthworm Eisenia foetida. Gene. 575:58–65. 2016. View Article : Google Scholar : PubMed/NCBI | |
Kim TH, Kim NS, Lim D, Lee KT, Oh JH, Park HS, Jang GW, Kim HY, Jeon M, Choi BH, et al: Generation and analysis of large-scale expressed sequence tags (ESTs) from a full-length enriched cDNA library of porcine backfat tissue. BMC Genomics. 7:362006. View Article : Google Scholar : PubMed/NCBI | |
Wellenreuther R, Schupp I, Poustka A and Wiemann S; German cDNA Consortium, : SMART amplication combined with cDNA size fractionation in order to obtain large full-length clones. BMC Genomics. 5:362004. View Article : Google Scholar : PubMed/NCBI | |
Dai Y, Su W, Yang C, Song B, Li Y and Fu Y: Development of novel polymorphic EST-SSR markers in bailinggu (Pleurotus tuoliensis) for crossbreeding. Genes (Basel). 8(pii): E3252017. View Article : Google Scholar : PubMed/NCBI | |
Liu C, Liu D, Guo Y, Lu T, Li X, Zhang M, Ma J, Ma Y and Guan W: Construction of a full-length enriched cDNA library and preliminary analysis of expressed sequence tags from Bengal Tiger Panthera tigris tigris. Int J Mol Sci. 14:11072–11083. 2013. View Article : Google Scholar : PubMed/NCBI | |
Carninci P, Shibata Y, Hayatsu N, Itoh M, Shiraki T, Hirozane T, Watahiki A, Shibata K, Konno H, Muramatsu M and Hayashizaki Y: Balanced-size and long-size cloning of full-length, cap-trapped cDNAs into vectors of the novel lambda-FLC family allows enhanced gene discovery rate and functional analysis. Genomics. 77:79–90. 2001. View Article : Google Scholar : PubMed/NCBI | |
Ling P, Wang M, Chen X and Campbell KG: Construction and characterization of a full-length cDNA library for the wheat stripe rust pathogen (Puccinia striiformis f. sp. tritici). BMC Genomics. 8:1452007. View Article : Google Scholar : PubMed/NCBI | |
Ye J, McGinnis S and Madden TL: BLAST: Improvements for better sequence analysis. Nucleic Acids Res 34 (Web Server Issue). W6–W9. 2006. View Article : Google Scholar | |
NCBI Resource Coordinators: Database resources of the national center for biotechnology information. Nucleic Acids Res. 41:(Database Issue). D8–D20. 2013. View Article : Google Scholar : PubMed/NCBI | |
Thompson JD, Higgins DG and Gibson TJ: CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:4673–4680. 1994. View Article : Google Scholar : PubMed/NCBI | |
Wang J, Duncan D, Shi Z and Zhang B: WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): Update 2013. Nucleic Acids Res 41 (Web Server Issue). W77–W83. 2013. View Article : Google Scholar | |
Zhang B, Kirov S and Snoddy J: WebGestalt: An integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res 33 (Web Server Issue). W741–W748. 2005. View Article : Google Scholar | |
Ravaux J, Hassanin A, Deutsch J, Gaill F and Markmann- Mulisch U: Sequence analysis of the myosin regulatory light chain gene of the vestimentiferan Riftia pachyptila. Gene. 263:141–149. 2001. View Article : Google Scholar : PubMed/NCBI | |
Kanehisa M, Goto S, Kawashima S and Nakaya A: The KEGG databases at GenomeNet. Nucleic Acids Res. 30:42–46. 2002. View Article : Google Scholar : PubMed/NCBI | |
Wiemann S, Mehrle A, Bechtel S, Wellenreuther R, Pepperkok R and Poustka A; German cDNA Consortium, : cDNAs for functional genomics and proteomics: The German Consortium. Comptes Rendus Biologies. 326:1003–1009. 2003. View Article : Google Scholar : PubMed/NCBI | |
Blair MW, Fernandez AC, Ishitani M, Moreta D, Seki M, Ayling S and Shinozaki K: Construction and EST sequencing of full-length, drought stress cDNA libraries for common beans (Phaseolus vulgaris L.). BMC Plant Biol. 11:1712011. View Article : Google Scholar : PubMed/NCBI | |
Clarke L and Carbon J: A colony bank containing synthetic CoI EI hybrid plasmids representative of the entire E. coli genome. Cell. 9:91–99. 1976. View Article : Google Scholar : PubMed/NCBI | |
Li YP, Xia RX, Wang H, Li XS, Liu YQ, Wei ZJ, Lu C and Xiang ZH: Construction of a full-length cDNA Library from Chinese oak silkworm pupa and identification of a KK-42-binding protein gene in relation to pupa-diapause termination. Int J Biol Sci. 5:451–457. 2009. View Article : Google Scholar : PubMed/NCBI | |
Hatey F, Tosser-Klopp G, Clouscard-Martinato K, Mulsant P and Gasser F: Expressed sequence tags for genes: A review. Genet Selec Evol. 30:521–524. 1998. View Article : Google Scholar | |
Thanh T, Chi VT, Abdullah MP, Omar H, Noroozi M, Ky H and Napis S: Construction of cDNA library and preliminary analysis of expressed sequence tags from green microalga Ankistrodesmus convolutus Corda. Mol Biol Rep. 38:177–182. 2011. View Article : Google Scholar : PubMed/NCBI | |
Duo L, Yin L, Zhang C and Zhao S: Ecotoxicological responses of the earthworm Eisenia fetida to EDTA addition under turfgrass growing conditions. Chemosphere. 220:56–60. 2018. View Article : Google Scholar : PubMed/NCBI | |
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 25:25–29. 2000. View Article : Google Scholar : PubMed/NCBI | |
Homa J, Stalmach M, Wilczek G and Kolaczkowska E: Effective activation of antioxidant system by immune-relevant factors reversely correlates with apoptosis of Eisenia andrei coelomocytes. J Comp Physiol B. 186:417–430. 2016. View Article : Google Scholar : PubMed/NCBI | |
Li P, Sheng J, Liu Y, Li J, Liu J and Wang F: Heparosan-derived heparan sulfate/heparin-like compounds: One kind of potential therapeutic agents. Med Res Rev. 33:665–692. 2013. View Article : Google Scholar : PubMed/NCBI | |
Lai CH, Han CK, Shibu MA, Pai PY, Ho TJ, Day CH, Tsai FJ, Tsai CH, Yao CH and Huang CY: Lumbrokinase from earthworm extract ameliorates second-hand smoke-induced cardiac fibrosis. Environ Toxicol. 30:1216–1225. 2015. View Article : Google Scholar : PubMed/NCBI |