Discovery of molecular associations among aging, stem cells, and cancer based on gene expression profiling

The emergence of a huge volume of "omics" data enables a computational approach to the investigation of the biology of cancer. The cancer informatics approach is a useful supplement to the traditional experimental approach. I reviewed several reports that used a bioinformatics approach to analyze the associations among aging, stem cells, and cancer by microarray gene expression profiling. The high expression of aging- or human embryonic stem cell-related molecules in cancer suggests that certain important mechanisms are commonly underlying aging, stem cells, and cancer. These mechanisms are involved in cell cycle regulation, metabolic process, DNA damage response, apoptosis, p53 signaling pathway, immune/inflammatory response, and other processes, suggesting that cancer is a developmental and evolutional disease that is strongly related to aging. Moreover, these mechanisms demonstrate that the initiation, proliferation, and metastasis of cancer are associated with the deregulation of stem cells. These findings provide insights into the biology of cancer. Certainly, the findings that are obtained by the informatics approach should be justified by experimental validation. This review also noted that next-generation sequencing data provide enriched sources for cancer informatics study.

Gene expression profiling defined pathways correlated with fibroblast cell proliferation induced by Opisthorchis viverrini excretory /secretory product

AIM： To investigate the mechanism of fibroblast cell proliferation stimulated by the Opisthorchis viverrini excretory/secretory （ES） product. METHODS： NIH-3T3, mouse fibroblast cells were treated with O. viverrini ES product by non-contact co-cultured with the adult parasites. Total RNA from NIH-3T3 treated and untreated with O. viverrini was extracted, reverse transcribed and hybridized with the mouse 15K complementary DNA （cDNA） array. The result was analyzed by ArrayVision version 5 and GeneSpring version 5 softwares. After normalization, the ratios of gene expression of parasite treated to untreated NIH-3T3 cells of 2-and more-fold upregulated was defined as the differentially expressed genes. The expression levels of the signal transduction genes were validated by semiquantitative SYBR-based real-time RT-PCR. RESULTS： Among a total of 15 000 genes/ESTs, 239 genes with established cell proliferation-related function were 2 fold-and more-up-regulated by O. viverrini ES product compared to those in cells without exposure to the parasitic product. These genes were classified into groups including energy and metabolism, signal transduction, protein synthesis and translation, matrix and structural protein, transcription control, cell cycle and DNA replication. Moreover, the expressions of serinethreonine kinase receptor, receptor tyrosine kinase and collagen production-related genes were up-regulated by O.viverrini ES product, The expression level of signal transduction genes, pkC, pdgfra, jak 1, eps 8, tgfβ 1/4,strap and h ras measured by real-time RT-PCR confirmed their expression levels to those obtained from cDNA array. However, only the up-regulated expression of pkC,eps 8 and tgfβ3 1/4 which are the downstream signaling molecules of either epidermal growth factor （EGF） or transforming growth factor-β （TGF-β） showed statistical significance （P 〈 0.05）.CONCLUSION： O. viverrini ES product stimulates the significant changes of gene expression in several functional categories and these mainly include transcripts related to cell proliferation. The TGF-β and EGF signal transduction pathways are indicated as the possible pathways of O. viverrini-driven cell proliferation.

GENE EXPRESSION PROFILING OF HUMAN PROMYELOCYTIC LEUKEMIA HL-60 CELL TREATED BY AJOENE

Objective: Ajoene, a major compound extracted from crashed garlic, has been shown to have antitumor, antimycotic, antimicrobial, antimutagenic functions in vivo or in vitro and treated as a potential antitumor drug. However, the molecular mechanisms underlying the tumor cytotoxicity of ajoene and even garlic substances are poorly defined. In the present study, we aimed to generate gene expression profiles of HL-60 cell treated by ajoene. Methods: A cDNA microarray presenting 2400 of genes amplified from human leukocyte cDNA library was constructed and the gene expression profiles of HL-60 cell induced by ajoene were generated. Results: After data analysis, 28 differentially expressed genes were identified and sequenced. These genes include 21 known genes and 7 ESTs. Most of the known genes are related to cell apoptosis, such as secretory granule (PRG1), beta-2 microglobulin (B2M), 16S ribosomal RNA gene and ribosomal protein S12. Several genes are related to cell differentiation, including the genes similar to H3 histone and ribosomal protein L31. Northern blot analysis was used to verify and quantify the expression of selected genes. Conclusion: Ajoene can induce HL-60 cell apoptosis significantly and may play a role in differentiation. cDNA microarray technology can be a valuable tool to gain insight into molecular events of pharmacological mechanism of herbal medicine.

GENE EXPRESSION PROFILING OF GANGLIOGLIOMA MALIGNANT PROGRESSION BY cDNA ARRAY

Objective: To establish gene expression profiles associated with malignant progression of ganglioglioma. Methods: The primary and two recurrent glioma specimens were collected intraoperatively from the same patient who experienced tumor transformation into anaplastic astrocytoma and glioblastoma multiform for the first and second recurrence respectively. Gene expression was assayed through cDNA array and bioinformatics analysis. Results: A total of 197 differentially expressed genes with differential ratio value more than 3 compared with normal brain tissue were obtained. Among 109 functionally denned genes, those associated with development ranked the first by frequency, followed by genes associated with metabolism, differentiation, signal transduction and so on. As a result of cluster analysis among 368 genes, eleven genes were up regulated with malignant progression, while six genes were down regulated. Conclusion: Gene expression profiles associated with malignant progression of glioma were successfully established, which provides a powerful tool for research on molecular mechanisms of malignant progression of gliomas.

Gene Expression Profiling of the Mouse Pancreas during the Secondary Transition in the Organogenesis of the Pancreatic Gland*

Diabetes mellitus is a chronic disease that impacts the homeostasis of blood sugar levels caused by loss or defect of insulin-producing β-cells in the Islets of Langerhans. Type 1 diabetes (T1D) is caused by auto-immune mediated destruction of β-cells, whereas in T2D, insulin is produced but used inefficiently. T2D accounts for 90% of people with diabetes worldwide (WHO 1999) and is the fastest increasing disease worldwide (https://diabetesatlas.org/en/). For an improved understanding of the pathomechanism of diabetes, profound knowledge of pancreas organogenesis and the associated gene regulatory networks is required. Therefore, we dissected and profiled the pancreatic endodermal and non-endodermal compartment between the embryonic stages (E) 12.5 and E 15.5 when progenitor cells commit to their different pancreatic lineages. Our associated study mined the global mRNA expression profile to increase the understanding of the secondary transition, endodermal-non-endodermal tissue interaction, and diabetic-related gene regulation. Furthermore, we validated 635 regulated pancreatic genes using the publicly available GenePaint.org, respective gp3.mpg.de to evaluate genes associated with genetic variants in Single-nucleotide polymorphism (SNP) related to T2D.

Effect of Chuanzhifang component (ZGC) on macrophage inflammatory injury based on whole gene expression profile

Objective: The effect of Chuanzhi Fang (ZGC) on the whole genome expression profile of RAW264.7 cells activated by lipopolysaccharide (LPS) was analyzed, and to explore the possible mechanism of action and core target of this formula on macrophage inflammatory injury at the overall level. Methods: A model of LPS-induced inflammation in RAW264.7 cells was constructed, and the effect of ZGC intervention on the genome-wide expression of inflammatory macrophages 3was examined by gene microarray technology, GO/KEGG enrichment analysis was performed for significantly differentially expressed genes among each group. Results: The results of genome-wide expression profiling microarray analysis showed that the ZGC intervention group upregulated the expression of 5 genes including C4bp and inhibited the expression of 22 genes including Mgat3, Psma6, and Siglecg relative to the LPS model group. KEGG signaling pathway analysis results showed that ZGC mainly acted through cytokine receptor interaction and the C-type lectin receptor signaling pathway. Conclusion: ZGC can interfere with the abnormal expression of 27 genes in inflammatory macrophages, and the related genes may exert corresponding anti-inflammatory effects by affecting cytokine receptor interactions, C-type lectin receptor signaling pathway, and TLR4/ NF-κB signaling pathway.

Identification of biomarkers of human pancreatic adenocarcinomas by expression profiling and validation with gene expression analysis in endoscopic ultrasound-guided fine needle aspiration samples

AIM： To compare gene expression profiles of pancreatic adenocarcinoma tissue specimens, human pancreatic and colon adenocarcinoma and leukemia cell lines and normal pancreas samples in order to distinguish differentially expressed genes and to validate the differential expression of a subset of genes by quantitative real-time RT-PCR （RT-QPCR） in endoscopic ultrasound-guided fine needle aspiration （EUS-guided FNA） specimens.METHODS： Commercially dedicated cancer cDNA macroarrays （Atlas Human Cancer 1.2） containing 1176 genes were used. Different statistical approaches （hierarchical clustering, principal component analysis （PCA） and SAM） were used to analyze the expression data. RT-QPCR and immunohistochemical studies were used for validation of results.RESULTS： RT-QPCR validated the increased expression of LCN2 （lipocalin 2） and for the first time PLAT （tissue-type plasminogen activator or tPA） in malignant pancreas as compared with normal pancreas. Immunohistochemical analysis confirmed the increased expression of LCN2 protein localized in epithelial cells of ducts invaded by carcinoma. The analysis of PLAT and LCN2 transcripts in 12 samples obtained through EUS-guided FNA from patients with pancreatic adenocarcinoma showed significantly increased expression levels in comparison with those found in normal tissues, indicating that a sufficient amount of high quality RNA can be obtained with this technique.CONCLUSION： Expression profiling is a useful method to identify biomarkers and potential target genes. Molecular analysis of EUS-guided FNA samples in pancreatic cancer appears as a valuable strategy for the diagnosis of pancreatic adenocarcinomas.

Molecular subtypes identified by gene expression profiling in early stage endometrioid endometrial adenocarcinoma

Background Early stage （FIGO stage Ⅰ-Ⅱ） endometrioid endometrial adenocarcinoma （EEA） is very common in clinical practice.However,patients with the early stage EEA show various clinical behaviors due to biological heterogeneity.Hence,we aimed to discover distinct classes of tumors based on gene expression profiling,and analyze whether the molecular classification correlated with the histopathological stages or other clinical parameters.Methods Hierarchical clustering was performed for class discovery in 28 eady stage EEA samples using a special cDNA microarray chip containing 492 genes designed for endometrial cancer.Correlations between clinicopathologic parameters and our classification were analyzed.And the significance analysis of microarrays （SAM） array was used to identify the signature genes according to the tumor grade and myometrial invasion.Results Three tumor subtypes （subtypes Ⅰ,Ⅱ and Ⅲ） were identified by hierarchical clustering,each subtype had different clinicopathological factors,such as tumor grade,myometrial invasion status,and FIGO stage.Moreover,SAM analysis showed 34 up-regulated genes in high grade tumors,and 38 up-regulated genes and 1 down-regulated in deep myometrial invasive tumors.The overlap genes between these two high-risk factors were markedly up-regulated in subtype Ⅰ,but down-regulated in subtype Ⅲ.Conclusion We have identified novel molecular subtypes in early stage EEA.Differential gene signatures characterize each tumor subtype,which could be used for recognizing the tumor risk and providing a basis for further treatment stratification.

Comprehensive genomic analysis and expression profiling of cysteine-rich polycomb-like transcription factor gene family in tea tree

Cysteine-rich polycomb-like(CPP)is a small gene family in plants,which plays key role in plant development and stress response.Although CPP transcription factors have been characterized in several other plant species,a genome-wide characterization of the CPP gene family has been absent in Camellia sinensis.In this study,we totally identified 7,8,and 8 non-redundant CsCPP genes in three published genomes,including Camellia sinensis var.assamica cv.Yunkang-10(CSA-YK10),Camellia sinensis var.sinensis cv.Biyun(CSS-BY)and Camellia sinensis var.sinensis cv.Shuchazao(CSS-SCZ).CPP proteins from tea tree and other plant species were classified into three groups,which were further divided into four subgroups based on phylogenetic relationships.Most CPP genes in the same subgroup had similar gene structures and conserved motifs.The cis-acting elements analysis indicated that CPP genes might be involved in plant growth,development and stress responses.Analysis of gene expression using qRT-PCR experiments validated that CPP genes exhibited different expression patterns across the examined tissues.All the genes were expressed differentially in a range of tissues,indicating that CPPs were involved in a range of developmental and physiological processes.This study has obtained new insights into the evolution and function of the CPP gene family in the growth and development of tea plants,and also provide candidate genes for further functional characterization in tea tree.

MP3RNA-seq:Massively parallel 3’end RNA sequencing for high-throughput gene expression profiling and genotyping

Transcriptome deep sequencing(RNA‐seq)hasbecome a routine method for global geneexpression profiling.However,its application tolarge‐scale experiments remains limited by costand labor constraints.Here we describe amassively parallel 3′end RNA‐seq(MP3RNA‐seq)method that introduces unique samplebarcodes during reverse transcription to permitsample pooling immediately following this initialstep.MP3RNA‐seq allows for handling of hun-dreds of samples in a single experiment,at acost of about$6 per sample for libraryconstruction and sequencing.MP3RNA‐seq iseffective for not only high‐throughput geneexpression profiling,but also genotyping.Todemonstrate its utility,we applied MP3RNA‐seqto 477 double haploid lines of maize.We iden-tified 19,429 genes expressed in at least 50%ofthe lines and 35,836 high‐quality singlenucleotide polymorphisms for genotypinganalysis.Armed with these data,we performedexpression and agronomic trait quantitativetrait locus(QTL)mapping and identified 25,797expression QTLs for 15,335 genes and 21 QTLsfor plant height,ear height,and relative earheight.We conclude that MP3RNA‐seq is highlyreproducible,accurate,and sensitive forhigh‐throughput gene expression profiling andgenotyping,and should be generally applicableto most eukaryotic species.

Gene expression profiling as a prognostic tool in multiple myeloma

Multiple myeloma(MM)is an aggressive plasma cell malignancy with high degrees of variability in outcome,some patients experience long remissions,whilst others survive less than two years from diagnosis.Therapy refractoriness and relapse remain challenges in MM management,and there is a need for improved prognostication and targeted therapies to improve overall survival(OS).The past decade has seen a surge in gene expression profiling(GEP)studies which have elucidated the molecular landscape of MM and led to the identification of novel gene signatures that predict OS and outperform current clinical predictors.In this review,we discuss the limitations of current prognostic tools and the emerging role of GEP in diagnostics and in the development of personalised medicine approaches to combat drug resistance.

Gene signatures to therapeutics:Assessing the potential of ivermectin against t(4;14)multiple myeloma

BACKGROUND Multiple myeloma(MM)is a terminal differentiated B-cell tumor disease characterized by clonal proliferation of malignant plasma cells and excessive levels of monoclonal immunoglobulins in the bone marrow.The translocation,(t)(4;14),results in high-risk MM with limited treatment alternatives.Thus,there is an urgent need for identification and validation of potential treatments for this MM subtype.Microarray data and sequencing information from public databases could offer opportunities for the discovery of new diagnostic or therapeutic targets.AIM To elucidate the molecular basis and search for potential effective drugs of t(4;14)MM subtype by employing a comprehensive approach.METHODS The transcriptional signature of t(4;14)MM was sourced from the Gene Expression Omnibus.Two datasets,GSE16558 and GSE116294,which included 17 and 15 t(4;14)MM bone marrow samples,and five and four normal bone marrow samples,respectively.After the differentially expressed genes were identified,the Cytohubba tool was used to screen for hub genes.Then,the hub genes were analyzed using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis.Using the STRING database and Cytoscape,protein–protein interaction networks and core targets were identified.Potential small-molecule drugs were identified and validated using the Connectivity Map database and molecular docking analysis,respectively.RESULTS In this study,a total of 258 differentially expressed genes with enriched functions in cancer pathways,namely cytokine receptor interactions,nuclear factor(NF)-κB signaling pathway,lipid metabolism,atherosclerosis,and Hippo signaling pathway,were identified.Ten hub genes(cd45,vcam1,ccl3,cd56,app,cd48,btk,ccr2,cybb,and cxcl12)were identified.Nine drugs,including ivermectin,deforolimus,and isoliquiritigenin,were predicted by the Connectivity Map database to have potential therapeutic effects on t(4;14)MM.In molecular docking,ivermectin showed strong binding affinity to all 10 identified targets,especially cd45 and cybb.Ivermectin inhibited t(4;14)MM cell growth via the NF-κB pathway and induced MM cell apoptosis in vitro.Furthermore,ivermectin increased reactive oxygen species accumulation and altered the mitochondrial membrane potential in t(4;14)MM cells.CONCLUSION Collectively,the findings offer valuable molecular insights for biomarker validation and potential drug development in t(4;14)MM diagnosis and treatment,with ivermectin emerging as a potential therapeutic alternative.

Analysis of Pathway Activity in Primary Tumors and NCI60 Cell Lines Using Gene Expression Profiling Data

To determine cancer pathway activities in nine types of primary tumors and NCI60 cell lines, we applied an in silico approach by examining gene signatures reflective of consequent pathway activation using gene expression data. Supervised learning approaches predicted that the Ras pathway is active in -70% of lung adenocarcinomas but inactive in most squamous cell carcinomas, pulmonary carcinoids, and small cell lung carcinomas. In contrast, the TGF-β, TNF-α, Src, Myc, E2F3, and β-catenin pathways are inactive in lung adenocarcinomas. We predicted an active Ras, Myc, Src, and/or E2F3 pathway in significant percentages of breast cancer, colorectal carcinoma, and gliomas. Our results also suggest that Ras may be the most prevailing oncogenic pathway. Additionally, many NCI60 cell lines exhibited a gene signature indicative of an active Ras, Myc, and/or Src, but not E2F3, β-catenin, TNF-α, or TGF-β pathway. To our knowledge, this is the first comprehensive survey of cancer pathway activities in nine major tumor types and the most widely used NCI60 cell lines. The ＂gene expression pathway signatures＂ we have defined could facilitate the understanding of molecular mechanisms in cancer development and provide guidance to the selection of appropriate cell lines for cancer research and pharmaceutical compound screening.

A Gene Expression Profiling of Early Rice Stamen Development that Reveals Inhibition of Photosynthetic Genes by OsMADS58

Stamen is a unique plant organ wherein germ cells or microsporocytes that commit to meiosis are initiated from somatic cells during its early developmental process. While genes determining stamen identity are known according to the ABC model of floral development, little information is available on how these genes affect germ cell initiation. By using the Affymetrix GeneChip Rice Genome Array to assess 51 279 tran- scripts, we established a dynamic gene expression profile （GEP） of the early developmental process of rice （Oryza sativa） stamen. Systematic analysis of the GEP data revealed novel expression patterns of some developmentally important genes including meiosis-, tapetum-, and phytohormone-related genes. Following the finding that a substantial amount of nuclear genes encoding photosynthetic proteins are ex- pressed at the low levels in early rice stamen, through the ChlP-seq analysis we found that a C-class MADS box protein, OsMADS58, binds many nuclear-encoded genes participated in photosystem and light reac- tions and the expression levels of most of them are increased when expression of OsMADS58 is downre- gulated in the osmads58 mutant. Furthermore, more pro-chloroplasts are observed and increased signals of reactive oxygen species are detected in the osmads58 mutant anthers. These findings implicate a novel link between stamen identity determination and hypoxia status establishment.

Characterization of the 19 Novel Cotton FLA Genes and Their Expression Profiling in Fiber Development and in Response to Phytohormones and Salt Stress

Fasciclin-like arabinogalactan proteins(FLAs),a subclass of arabinogalactan proteins(AGPs),are usually involved in cell development in plants.To investigate the expression profiling as well

Expression Profiling Identifies Candidate Genes for Fiber Yield and Quality

Gene expression profiling at early stages(0～2 DPA) of fiber development in Gossypium hirsutum identified a number of transcription factors which were down regulated in fiberless mutants relative to wild type controls and which could play a role in controlling early fiber development.Chief among these was GhMYB25,a Mixta-like MYB gene.Transgenic GhMYB25-silenced cotton

Identifying genetic susceptibility to Aspergillus fumigatus infection using collaborative cross mice and RNA-Seq approach

Background:Aspergillus fumigatus(Af)is one of the most ubiquitous fungi and its infection potency is suggested to be strongly controlled by the host genetic back-ground.The aim of this study was to search for candidate genes associated with host susceptibility to Aspergillus fumigatus(Af)using an RNAseq approach in CC lines and hepatic gene expression.Methods:We studied 31 male mice from 25 CC lines at 8 weeks old;the mice were infected with Af.Liver tissues were extracted from these mice 5 days post-infection,and next-generation RNA-sequencing(RNAseq)was performed.The GENE-E analysis platform was used to generate a clustered heat map matrix.Results:Significant variation in body weight changes between CC lines was ob-served.Hepatic gene expression revealed 12 top prioritized candidate genes differ-entially expressed in resistant versus susceptible mice based on body weight changes.Interestingly,three candidate genes are located within genomic intervals of the previ-ously mapped quantitative trait loci(QTL),including Gm16270 and Stox1 on chromo-some 10 and Gm11033 on chromosome 8.Conclusions:Our findings emphasize the CC mouse model's power in fine mapping the genetic components underlying susceptibility towards Af.As a next step,eQTL analysis will be performed for our RNA-Seq data.Suggested candidate genes from our study will be further assessed with a human cohort with aspergillosis.

Expression Profile Changes of Genes Involved in Lipid Metabolism Pathway During Liver Regeneration in Mice

[ Objective ] The aim of the research was to study the expression profile changes of genes involved in lipid metabolism pathway during liver regeneration in mice. [ Method] The CCI4 induced mouse model of liver regeneration was established and the total RNA was isolated from liver tissue of mouse. Then the changes of genes involved in lipid metabolism pathway during different stages of liver regeneration were detected through micro-array chip gene technique and their specific functions were also analyzed. [ Result] Dudng the process of liver regeneration, the expression level of 98 genes involved in lipid metabolism pathway changed, which were divided into eight groups according to change trend. In the mass, the expression of genes was inhibited in the early stage and up-regulated in the late phase. And the gene expression associated with fatty acid synthesis pathway was mainly up-regulated while the catabolic pathway did not change significantly. Most of genes involved in bile acid synthesis pathway were suppressed before 4.5 d and up-regulated after 4.5 d or 7 d. [ Conclusion] During the process of liver regeneration, the genes associated with lipid metabolism are expressed in different trends, and this data should provide a specific range of genes for further studying the regulation effect of lipid metabolism related pathway on liver regeneration.

Discovery and validation of prognostic markers in gastric cancer by genome-wide expression profiling

AIM: To develop a prognostic gene set that can predict patient overall survival status based on the whole genome expression analysis. METHODS: Using Illumina HumanWG-6 BeadChip followed by semi-supervised analysis, we analyzed the expression of 47 296 transcripts in two batches of gastric cancer patients who underwent surgical resection. Thirty-nine samples in the first batch were used as the training set to discover candidate markers correlated to overall survival, and thirty-three samples in the second batch were used for validation. RESULTS: A panel of ten genes were identified as prognostic marker in the first batch samples and classified patients into a lowand a high-risk group with significantly different survival times (P = 0.000047). This prognostic marker was then verified in an independent validation sample batch (P = 0.0009). By comparing with the traditional Tumor-node-metastasis (TNM) staging system, this ten-gene prognostic marker showed consistent prognosis results. It was the only independent prognostic value by multivariate Cox regression analysis (P = 0.007). Interestingly, six of these ten genes are ribosomal proteins, suggesting a possible association between the deregulation of ribosome related gene expression and the poor prognosis. CONCLUSION: A ten-gene marker correlated with overall prognosis, including 6 ribosomal proteins, was identified and verified, which may complement the predictive value of TNM staging system.

Genes Expression Profile Difference in Peripheral Blood Between Esophageal Carcinoma Patients and Normal Subjects

Objective： To study the genes expression profile differences in the peripheral blood between esophageal carcinoma patients and normal subjects using the gene chip technique and screen out the esophageal early conceration associated genes. Methods： The total RNA was extracted and purified in the peripheral blood obtained from the patients with esophageal carcinoma and normal subjects. The first strand of cDNA was synthesized through retro-transcription and labeled with Cy5 and Cy3 fluorescence as probes. The mixed probes were hybridized with a piece of 4096 double dot human whole gene chip. The acquired image was analyzed by microarrav suite software using a digital computer, and the intensity of ttuorescence signal and its ratio were calculated. Results： A total of 92 genes were screened out and its expression difference was more than 2 times in the peripheral blood between the patients with esophageal carcinoma and normal subjects. Among these, the expression difference of 36 genes was more than 3 times. Two human urokinase plasminogen activator surface receptor （UPAR） genes, 80K-L protein gene, human protein tyrosine-phosphatase gent arid proto-oncogene protein mRNA were significantly up-regulated, while the collagen V type （α-2 gene was markedly down-regulated. Conclusion： 80K-L protein gene, tyrosinephophatase gene, proto-oncogene protein arid the collagen V type α-2 gene might be associated with the ontogenesis, development and its metastasis in the esophageal carcinoma. The UPAR gene may play important roles in the diagnosing the micrometastasis in the peripheral blood of esophageal carcinoma.