SEARCH · Search PubMed

Results for “transcriptome analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76Linked to original sources

PubMed2026

Increased expression of Ribonucleic acid export 1 (RAE1) gene promotes gastric carcinogenesis and is associated with Hippo signaling pathway.

BACKGROUND: Ribonucleic acid export 1 (RAE1) autoantibody may have good potential for early detection of gastric cancer (GC). However, the carcinogenicity of RAE1 in GC remains unknown. We aimed to explore the role and the potential mechanism for RAE1 in the carcinogenesis of GC. METHODS: Immunohistochemical assay was applied to analyze the expression of RAE1 in GC and precancerous lesion (PL) tissues and its relationship with clinical characteristics. The effects of RAE1 on proliferation, migration, apoptosis, and cell cycle were explored by constructing RAE1 knockdown and overexpression in GC cells. The effects of RAE1 knockdown on tumor growth were observed in a murine xenograft model. The signaling pathways involved in GC development that may be affected by RAE1 were investigated by transcriptome sequencing and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. RESULTS: Immunohistochemical expression of RAE1 was significantly higher in early GC patients than in PL and normal tissues, and the RAE1 expression was correlated with clinical stage (P<0.001). In vitro, knockdown of RAE1 inhibited the proliferation and migration of GC cells with promoted apoptosis and arrested cell cycle in the S phase, while overexpression of RAE1 showed the opposite trend. In vivo, knockdown of RAE1 inhibited the growth of subcutaneous graft tumors in mice. Transcriptome sequencing and subsequent analysis of RAE1 knockdown cells revealed that the Hippo signaling pathway was activated by RAE1 knockdown. CONCLUSIONS: RAE1 promotes GC cells proliferation and migration and is associated with the inhibition of the Hippo signaling pathway, and may be a potential biomarker for early diagnosis and treatment of GC.

Gastric cancer (GC)↗

PubMed2006

p53 suppresses the self-renewal of adult neural stem cells.

There is increasing evidence that tumors are heterogeneous and that a subset of cells act as cancer stem cells. Several proto-oncogenes and tumor suppressors control key aspects of stem cell function, suggesting that similar mechanisms control normal and cancer stem cell properties. We show here that the prototypical tumor suppressor p53, which plays an important role in brain tumor initiation and growth, is expressed in the neural stem cell lineage in the adult brain. p53 negatively regulates proliferation and survival, and thereby self-renewal, of neural stem cells. Analysis of the neural stem cell transcriptome identified the dysregulation of several cell cycle regulators in the absence of p53, most notably a pronounced downregulation of p21 expression. These data implicate p53 as a suppressor of tissue and cancer stem cell self-renewal.

Animals↗

PubMed2025

Spatial transcriptomics-aided localization for single-cell transcriptomics with STALocator.

Single-cell RNA-sequencing (scRNA-seq) techniques can measure gene expression at single-cell resolution but lack spatial information. Spatial transcriptomics (ST) techniques simultaneously provide gene expression data and spatial information. However, the data quality of the spatial resolution or gene coverage is still much lower than the quality of the single-cell transcriptomics data. To this end, we develop a ST-Aided Locator for single-cell transcriptomics (STALocator) to localize single cells to corresponding ST data. Applications on simulated data showed that STALocator performed better than other localization methods. When applied to the human brain and squamous cell carcinoma data, STALocator could robustly reconstruct the relative spatial organization of critical cell populations. Moreover, STALocator could enhance gene expression patterns for Slide-seqV2 data and predict genome-wide gene expression data for fluorescence in situ hybridization (FISH) and Xenium data, leading to the identification of more spatially variable genes and more biologically relevant Gene Ontology (GO) terms compared with the raw data. A record of this paper's transparent peer review process is included in the supplemental information.

Single-Cell Analysis↗

PubMed2002

Oligonucleotide microarray analysis of gene expression in follicle-stimulating hormone-treated rat Sertoli cells.

Spermatogenesis requires the presence of functional somatic Sertoli cells in the seminiferous tubules of the testis. Sertoli cells provide support and factors necessary for the successful progression of germ cells into spermatozoa. Sertoli cells are regulated to a large degree by the glycoprotein hormone FSH, which is required for the testis to acquire full size and spermatogenic capacity. Signaling events initiated by the binding of FSH to its receptor lead to an alteration of Sertoli cell gene expression. To characterize the changes in gene expression in FSH-treated Sertoli cells, we used the mRNA from these cells to screen Affymetrix U34A rat GeneChip oligonucleotide microarrays. Sertoli cells from 20-d-old rats were cultured in the presence of 25 ng/ml ovine FSH. At 0, 2, 4, 8, and 24 h after the addition of FSH, total RNA was purified and used to prepare biotinylated target, which was hybridized to the U34A rat microarray containing approximately 9000 rat genes. Analysis identified 100-300 transcripts at each time point that were up-regulated or down-regulated by 2-fold or greater. Genes previously reported to be FSH or cAMP regulated in rat Sertoli cells were identified, in addition to numerous genes not reported to be expressed or FSH regulated in Sertoli cells. The expression patterns of five of these genes, encoding nerve growth factor inducible gene B, PRL-1, PC3 nerve growth factor-inducible antiproliferative putative secreted protein, diacylglycerol acyltransferase, and an expressed sequence tag, in FSH- and N,O'-dibutyryl cAMP-treated rat Sertoli cells were confirmed and characterized by Northern blot analysis. Thus, we have begun to define the transcriptome induced and repressed by FSH in rat Sertoli cells, and we have generated datasets of genes available for further analysis in regard to spermatogenesis and Sertoli cell signaling.

Acyltransferases↗

PubMed2026

Single-cell analysis of the human retina reveals stage-linked microglial states and neural-immune circuit rewiring in diabetic retinopathy.

Diabetic retinopathy (DR) is a major cause of vision loss worldwide. Here, we conduct single-cell RNA sequencing of twenty human retina samples (from living and post-mortem donors) across non-diabetic, diabetic, and DR states to create a comprehensive transcriptomic atlas. We identify two stable microglial populations-homeostatic and inflammatory-that exist along a functional continuum, plus a neutrophil cluster within C1QA+ myeloid cells with dynamic transitions occurring throughout disease progression. Module-level analysis reveals divergent transcriptional trajectories: homeostatic microglia maintain energetic programs while selectively upregulating stress elements, whereas inflammatory microglia layer additional pro-inflammatory programs onto preserved biosynthetic foundations. Eleven co-expression modules organize into two major axes: an inflammatory-stress axis, and a regulatory/metabolic-motility axis, with a stable translation module persisting across disease stages. Cell communication analysis further highlights sophisticated neural-immune interactions, particularly between photoreceptors and microglia. Our findings provide insights into the complex cellular dynamics of DR progression and suggest potential therapeutic targets for early intervention.

Humans↗

PubMed2026

ORFannotate: reproducible coding sequence annotation of transcriptome assemblies.

SUMMARY: Accurate annotation of coding sequences and translational features within transcript models is essential for interpreting assembled transcriptomes and their functional potential. Existing open reading frame (ORF) prediction tools typically operate on transcript FASTA files and do not reintegrate coding sequence (CDS) information back into transcript models, limiting their utility in long-read sequencing workflows where GTF/GFF annotations are the primary output. We present ORFannotate, a lightweight, GTF-native Python command-line tool that predicts ORFs from transcript annotations and reinserts precise, exon-aware CDS and UTR features into the original GTF/GFF file. In addition, ORFannotate provides biologically informative translational context by annotating Kozak sequence strength, detecting non-overlapping upstream ORFs (uORFs) with coding probabilities, characterising 5' and 3' untranslated regions (UTRs), and predicting nonsense-mediated decay (NMD) susceptibility. All annotations are consolidated in a transcript-level summary to support downstream analysis. By generating GTF files with accurate CDS annotations, ORFannotate facilitates reproducible analysis of both long- and short-read transcriptomes and integrates seamlessly with visualization tools, genome browsers, and comparative transcript analysis workflows. ORFannotate is fast, scalable and provides a practical solution for transcriptome annotation beyond coding potential prediction alone. AVAILABILITY AND IMPLEMENTATION: ORFannotate is implemented in Python and freely available under the GNU General Public License v3 (GPL-3.0) at: https://github.com/egustavsson/ORFannotate (DOI: https://doi.org/10.5281/zenodo.16812866).

Open Reading Frames↗

PubMed2003

A panoramic view of gene expression in the human kidney.

To gain a molecular understanding of kidney functions, we established a high-resolution map of gene expression patterns in the human kidney. The glomerulus and seven different nephron segments were isolated by microdissection from fresh tissue specimens, and their transcriptome was characterized by using the serial analysis of gene expression (SAGE) method. More than 400,000 mRNA SAGE tags were sequenced, making it possible to detect in each structure transcripts present at 18 copies per cell with a 95% confidence level. Expression of genes responsible for nephron transport and permeability properties was evidenced through transcripts for 119 solute carriers, 84 channels, 43 ion-transport ATPases, and 12 claudins. Searching for differences between the transcriptomes, we found 998 transcripts greatly varying in abundance from one nephron portion to another. Clustering analysis of these transcripts evidenced different extents of similarity between the nephron portions. Approximately 75% of the differentially distributed transcripts corresponded to cDNAs of known or unknown function that are accurately mapped in the human genome. This systematic large-scale analysis of individual structures of a complex human tissue reveals sets of genes underlying the function of well-defined nephron portions. It also provides quantitative expression data for a variety of genes mutated in hereditary diseases and helps in sorting candidate genes for renal diseases that affect specific portions of the human nephron.

Cluster Analysis↗

PubMed2003

A snapshot of the low temperature stress transcriptome of developing rice seedlings (Oryza sativa L.) via ESTs from subtracted cDNA library.

Rice (Oryza sativa L.) is sensitive to chilling particularly during early seedling development. Given the biochemical complexity of tolerance mechanisms, genetic potential for this trait depends on highly coordinated expression of many genes. We used a simple cDNA subtraction strategy to develop Expressed Sequence Tags (ESTs) that represent an important subset of cold stress-upregulated genes. The 3,084 subtracted cDNA clones represent a total of 1,967 unigenes from 1,354 singletons and 613 contigs. As expected in the developing seedlings, genes involved in basic cellular processes, i.e., metabolism, growth and development, protein synthesis, folding and destination, cellular transport, cell division and DNA replication were widely represented. Genes with stress-related and regulatory functions comprised 23.17% of the total ESTs. These categories included proteins with known function in cellular defenses against abiotic (drought, cold and salinity) and biotic (pathogen) stresses, and proteins involved in developmental and stress response signalling and transcription. Based on the types of genes represented, tolerance mechanisms rely on precise integration of developmental processes with stress-related responses. A large fraction of the ESTs (38.7%) represents unknown proteins. This EST library is a rich source of cold stress-related genes, and supplements for other publicly available libraries for comprehensive analysis of the stress-response transcriptome.

Cold Temperature↗

PubMed2006

Genome-wide analysis of group a streptococci reveals a mutation that modulates global phenotype and disease specificity.

Many human pathogens produce phenotypic variants as a means to circumvent the host immune system and enhance survival and, as a potential consequence, exhibit increased virulence. For example, it has been known for almost 90 y that clinical isolates of the human bacterial pathogen group A streptococci (GAS) have extensive phenotypic heterogeneity linked to variation in virulence. However, the complete underlying molecular mechanism(s) have not been defined. Expression microarray analysis of nine clinical isolates identified two fundamentally different transcriptomes, designated pharyngeal transcriptome profile (PTP) and invasive transcriptome profile (ITP). PTP and ITP GAS differed in approximately 10% of the transcriptome, including at least 23 proven or putative virulence factor genes. ITP organisms were recovered from skin lesions of mice infected subcutaneously with PTP GAS and were significantly more able to survive phagocytosis and killing by human polymorphonuclear leukocytes. Complete genome resequencing of a mouse-derived ITP GAS revealed that the organism differed from its precursor by only a 7-bp frameshift mutation in the gene (covS) encoding the sensor kinase component of a two-component signal transduction system implicated in virulence. Genetic complementation, and sequence analysis of covR/S in 42 GAS isolates confirmed the central role of covR/S in transcriptome, exoproteome, and virulence modulation. Genome-wide analysis provides a heretofore unattained understanding of phenotypic variation and disease specificity in microbial pathogens, resulting in new avenues for vaccine and therapeutics research.

Journal Article↗

PubMed2026

A latent activated olfactory stem cell state revealed by single-cell transcriptomic and epigenomic profiling.

The olfactory epithelium is one of the few regions of the nervous system that sustains neurogenesis throughout life. Its experimental accessibility makes it especially tractable for studying molecular mechanisms that drive neural regeneration in response to injury. In this study, we used single-cell sequencing to identify transcriptional and epigenetic processes involved in determining olfactory epithelial stem cell fate during injury-induced regeneration. By combining gene expression and accessible chromatin profiles of individual lineage-traced olfactory stem cells, we identified transcriptional heterogeneity among activated stem cells at a stage when cell fates are being specified. We further identified a subset of resting cells that appears poised for activation, characterized by accessible chromatin around silent genes prior to their expression in response to injury. These results provide evidence for a latent activated stem cell state in which a subset of quiescent olfactory epithelial stem cells are epigenetically primed to support injury-induced regeneration.

Animals↗

PubMed2025

PUS7-dependent Ψ reshapes specific synaptic gene exons to facilitate fear extinction memory formation.

RNA modifications serve as dynamic regulators of neural plasticity through their ability to fine-tune transcript stability and splicing. Pseudouridine (Ψ), an evolutionarily conserved RNA modification catalyzed by pseudouridine synthases, plays established roles in neurodevelopment, yet its functional significance in activity-dependent behavioral adaptation remains poorly defined. Here, we investigate Ψ-mediated epitranscriptomic regulation within the infralimbic prefrontal cortex (ILPFC), a brain region requiring precise synaptic remodeling for the clinically relevant form of fear extinction memory. Combining transcriptome-wide pseudouridylation profiling with behavioral analysis in mice, we identified selective Ψ enrichment at exons of synaptic regulatory genes within ILPFC during fear extinction learning. Fear extinction in the ILPFC drives concomitant exonic Ψ deposition and upregulation of synaptogenic transcripts, processes that involve pseudouridine synthase PUS7. Crucially, PUS7 knockdown in the ILPFC selectively impaired fear extinction memory formation without altering baseline fear expression, establishing a causal link between Ψ-dependent RNA processing and activity-dependent synaptic structural remodeling in this microcircuit. Our findings demonstrate that PUS7-mediated Ψ modification spatiotemporally regulates activity-dependent RNA dynamics in the ILPFC, providing the evidence that epitranscriptomic mechanisms precisely coordinate synaptic gene expression within behaviorally defined brain sub-region. This work bridges molecular RNA biology with systems neuroscience, revealing a novel mechanism for activity-dependent regulation of fear extinction in ILPFC.

Animals↗

PubMed2026

Molecular Evolution and Expression Analysis of the ADH Gene Family in Apple Bud Mutants.

Alcohol dehydrogenase (ADH) catalyzes the reduction of aldehydes to alcohols, key precursor substrates for volatile ester biosynthesis, which determines the characteristic aroma of apple fruit. However, a comprehensive genome-wide investigation of the ADH gene family in apple has been lacking. In this study, we systematically identified ADH genes in the apple genome using integrated bioinformatics approaches, including phylogenetic analysis, synteny evaluation, promoter cis-element prediction, codon usage bias assessment, and protein interaction network modeling. Expression patterns were examined through transcriptomic data and validated by RT-qPCR analysis across different organs and among 'Red Delicious' and its four bud mutant lines. We identified 44 ADH genes, with 12 forming a prominent cluster on chromosome 1. RT-qPCR analysis revealed that MdADH20 was dramatically upregulated in the 'Red Chief' mutant (relative expression of 59.38), suggesting its pivotal role. Phylogenetic analysis revealed a close evolutionary relationship with wild strawberry. The encoded proteins were generally stable and predominantly localized to the cytoplasm. Promoter analysis showed enrichment of growth/development-related and ARE elements, while codon usage analysis identified AGA, GCU, GUU, and CUU as preferred codons. Protein interaction prediction suggested MdADH19 and MdADH20 as hub proteins. Expression profiling and RT-qPCR further identified MdADH20 as a core candidate gene, characterized by its stable and high expression, particularly in the 'Red Delicious' mutant. Its central position in the predicted protein-protein interaction network suggests a potential regulatory role in the aroma biosynthesis pathway of apple fruit. This study provides the first systematic genome-wide characterization of the apple ADH gene family, establishing a theoretical groundwork for deciphering aroma biosynthesis mechanisms and offering potential target genes for flavor improvement through bud mutation breeding strategies.

ADH gene family↗

PubMed2025

GZMK+CD8+ T cells target a specific acinar cell type in Sjögren's disease.

OBJECTIVES: Sjögren's disease (SjD) is a systemic autoimmune disorder characterized by dysfunction of exocrine glands, particularly the salivary and lacrimal glands, with no clear etiology or effective therapy. This study explores the complex interplay of varied cell types in the salivary glands and their role in the pathology of Sjögren's disease. METHODS: Utilizing single-cell and spatial transcriptomics alongside spatial immunophenotyping to analyze human minor salivary glands, we developed a comprehensive understanding of the cellular landscape of non-SjD salivary glands and how that landscape changes in SjD patients. In vitro cellular assays and novel patient-derived primary epithelial cells were co-cultured with autologous T cells to confirm effector states and the delivery and effect of disease-associated granzymes. RESULTS: We identified previously unrecognized heterogeneity among acinar cells, including a PRR4⁺CST3⁺WFDC2⁻ seromucous acinar population that is selectively lost in Sjögren's disease. Expression and organizational changes were linked to clinical features: (i) T cells in the glands of SSA⁺, high-focus score patients showed increased transcriptional signatures of activation, antigen presentation, and apoptosis resistance compared with patients with mild or moderate disease, and (ii) patients with low immune infiltration exhibited distinct epithelial organization. Notably, GZMK⁺CD8⁺ T cells, which accumulate with disease severity, displayed a cytotoxic transcriptional program, degranulated upon stimulation ex vivo, and localized spatially with immune-engaged epithelial cells. Functional assays demonstrated that GZMK activates interferon signaling in vitro, and autologous co-cultures of patient-derived T cells and epithelial cells validated these findings. CONCLUSIONS: Using single-cell and spatial transcriptomics and proteomics, this study identifies a selective loss of PRR4⁺CST3⁺WFDC2⁻ seromucous acinar cells and a rise in GZMK⁺CD8⁺ T cells in Sjögren's disease, revealing distinct immune-mediated epithelial remodeling and interferon-driven dysfunction across diverse clinical presentations. These findings uncover a novel sub-cytolytic effector mechanism by which GZMK⁺CD8⁺ T cells impair mitochondrial integrity and activate innate immune signaling, linking epithelial injury to type I interferon responses and offering new therapeutic targets.

Humans↗

PubMed2026

Optimized Hot Phenol-Based RNA Extraction from Mycobacteria: A Robust Approach for Reliable Gene Expression Analysis.

Mycobacterium tuberculosis (Mtb) remains a major global health threat, underscoring the need for reliable transcriptomic studies to understand its biology and drug resistance mechanisms. Such analyses depend on obtaining high-quality, high-yield RNA. Although several RNA extraction methods are available, many require expensive reagents, large culture volumes, or specialized equipment, limiting their suitability for large-scale studies, particularly in resource-constrained settings. Here, an optimized Hot Phenol based RNA extraction method specifically tailored for mycobacteria is presented. The method uses minimal culture volume and commonly available reagents to consistently yield high-quality RNA suitable for high-throughput transcriptomic applications. RNA quantity and integrity were assessed by gel electrophoresis and RNA integrity analysis (RIN), and its suitability for downstream applications was confirmed by qPCR and Qubit 4. To benchmark the performance of the optimized method, a parallel RNA extraction using TRIzol and RNeasy under identical experimental conditions was carried out, including the same Mycobacterium species, culture volume, growth phase (logarithmic and stationary), and lysis conditions. This allowed a direct comparison of yield, quality, feasibility, and cost. The optimized Hot Phenol method demonstrated comparable or improved RNA yield and quality while significantly reducing reagent cost and dependence on specialized equipment. Owing to its efficiency, reproducibility, and affordability, this protocol provides a practical alternative for large-scale gene expression and transcriptomic studies in Mtb and other mycobacterial species.

RNA, Bacterial↗

PubMed2026

A reproducible computational transcriptomic framework for cell-type-resolved fibroinflammatory-AKT remodeling in human heart failure.

BACKGROUND: Human heart failure involves multicellular transcriptional remodeling, but public transcriptomic studies often remain disconnected from cell-type localization and perturbational interpretation. METHODS: We developed a reproducible computational workflow integrating human left-ventricular bulk transcriptomes, donor-level cell-type pseudobulk results from a human heart-failure single-cell/single-nucleus atlas, external snRNA-seq support, curated module scoring, focused ligand-receptor prioritization and LINCS/L1000 perturbational matching. RESULTS: Cross-cohort analysis identified 14,358 same-direction HF-associated genes, including 1633 replicated HF-up and 785 replicated HF-down genes. Donor-level pseudobulk analysis localized disease remodeling to cardiomyocyte, fibroblast and myeloid compartments. Activated fibroblast and inflammatory myeloid programs defined a fibroinflammatory remodeling axis connected to context-dependent AKT-associated transcriptional shifts. External snRNA-seq support was strongest for fibroblast activation and AKT-associated remodeling, with etiology-dependent heterogeneity across validation resources. L1000FWD screening prioritized safety-aware perturbational hypotheses, including glimepiride and simvastatin as interpretable candidates requiring experimental validation. CONCLUSIONS: This study provides a computational transcriptomic framework linking reproducible human HF signatures, cell-type-resolved fibroinflammatory remodeling and perturbational genomic prioritization without claiming drug efficacy or AKT causality.

Humans↗

PubMed2005

Serial analysis of gene expression study of a hybrid rice strain (LYP9) and its parental cultivars.

Using the serial analysis of gene expression technique, we surveyed transcriptomes of three major tissues (panicles, leaves, and roots) of a super-hybrid rice (Oryza sativa) strain, LYP9, in comparison to its parental cultivars, 93-11 (indica) and PA64s (japonica). We acquired 465,679 tags from the serial analysis of gene expression libraries, which were consolidated into 68,483 unique tags. Focusing our initial functional analyses on a subset of the data that are supported by full-length cDNAs and the tags (genes) differentially expressed in the hybrid at a significant level (P<0.01), we identified 595 up-regulated (22 tags in panicles, 228 in leaves, and 345 in roots) and 25 down-regulated (seven tags in panicles, 15 in leaves, and three in roots) in LYP9. Most of the tag-identified and up-regulated genes were found related to enhancing carbon- and nitrogen-assimilation, including photosynthesis in leaves, nitrogen uptake in roots, and rapid growth in both roots and panicles. Among the down-regulated genes in LYP9, there is an essential enzyme in photorespiration, alanine:glyoxylate aminotransferase 1. Our study adds a new set of data crucial for the understanding of molecular mechanisms of heterosis and gene regulation networks of the cultivated rice.

Chromosome Mapping↗

PubMed2005

Overview and perspectives the transcriptome of Paracoccidioides brasiliensis.

Paracoccidioides brasiliensis is a dimorphic and thermo-regulated fungus which is the causative agent of paracoccidioidomycosis, an endemic disease widespread in Latin America that affects 10 million individuals. Pathogenicity is assumed to be a consequence of the dimorphic transition from mycelium to yeast cells during human infection. This review shows the results of the P. brasiliensis transcriptome project which generated 6,022 assembled groups from mycelium and yeast phases. Computer analysis using the tools of bioinformatics revealed several aspects from the transcriptome of this pathogen such as: general and differential metabolism in mycelium and yeast cells; cell cycle, DNA replication, repair and recombination; RNA biogenesis apparatus; translation and protein fate machineries; cell wall; hydrolytic enzymes; proteases; GPI-anchored proteins; molecular chaperones; insights into drug resistance and transporters; oxidative stress response and virulence. The present analysis has provided a more comprehensive view of some specific features considered relevant for the understanding of basic and applied knowledge of P. brasiliensis.

Cell Wall↗

PubMed2006

Gene expression patterns in Euglena gracilis: insights into the cellular response to environmental stress.

To better understand Euglena gracilis gene expression under different stress conditions (Chromium, Streptomycin or darkness), we undertook a survey of the E. gracilis transcriptome by cDNA sequencing and microarray analysis. First, we constructed a non-normalized cDNA library from the E. gracilis UTEX strain and sequenced a total of 1000 cDNAs. Six hundred and ten of these ESTs were similar to either Plantae or Protistae genes (e-value<e(-10)). Second, microarrays were built by spotting all the ESTs onto mirror slides. Microarray expression analysis indicated that 90 out of those 610 ESTs changed their expression level in response to different stress treatments (p<0.05). In addition, we detected 10 ESTs that changed expression levels irrespective of the tested stress. These may be considered as part of a larger set of stress-related genes in E. gracilis. Finally, we identified 23 unknown ESTs (U-ESTs) following the expression profiles of these putative stress-related genes suggesting that they could be related to the cellular mechanism of stress response.

Animals↗