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.

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.

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 gene expression profiles in pancreatic carcinoma by using cDNA microarray

OBJECTIVES: To survey the gene expression profiles in pancreatic carcinoma by using cDNA microarray and detect target genes for further study. METHODS: Three mixed samples from 2 cases of normal pancreatic tissue and 4 cases of moderate-differentiated pancreatic carcinoma were studied by means of cDNA microarray consisting of 18 000 genes. RESULTS: 1484 and 1353 different expressed genes were observed in two cancer samples respectively. We identified 455 genes altered with the same tendency in both samples, including 102 up-regulated and 353 down-regulated genes. There were 274 known genes and 181 unknown genes; 27.8% and 52.0% genes respectively had an expression level in cancer that was 2-fold higher or lower than that in normal samples. Tumor suppressor genes, growth factors and receptor genes, signal conduction genes, transcription factor genes were identified. CONCLUSIONS: cDNA microarray is an efficient and high-throughout method to investigate gene expression profiles in pancreatic carcinoma. MBD1, EDG1 and gene hypermethylation mechanism would play an important role in the pathogenesis of pancreatic carcinoma.

Differential Gene Expression Profiles in Acute Hepatic Failure Model in Mice Infected with MHV-3 Virus Intervened by Anti-hepatic Failure Compound

Differential gene expression profiles in Balb/cJ mouse model of acute hepatic failure infected with MHV-3 virus intervened by anti-hepatic failure compound （AHFC） and the changes of cytokines regulated by genes were investigated. The Balb/cj mice were divided into AHFC-intervened group and control group randomly. Acute hepatic failure model of Balb/cJ mice infected with MHV-3 virus was established. The survival rate in the two groups was observed. It was found that the survival rate in the AHFC-intervened group and control group was 90% and 50% respectively 48 h after intraperitoneal injection of MHV-3 （P〈0.05）. Before and after the experiment, the cytokines in peripheral blood of the survival mice were determined, and RNA was extracted from survival mouse liver tissue for the analysis of the differential gene expression by a 36 kb mouse oligonuleotide DNA array. In all the genes of microarray there were 332 genes expressed differently in the two groups, in which 234 genes were up-regulated and 78 genes down-regulated. Through clustering analysis, the differential expression of immune related genes, including TNF receptor superfamily, Kctd9, Bcl-2, Fgl2, IL-8, IL-6, IFN-7, TNF-α etc. might be related with the curative effectiveness of AHFC. It was suggested that AHFC can balance the immune state of mouse model of acute hepatic failure infected with MHV-3 virus mainly through regulating the expression of immune related genes, decrease the immune damage and inhibit liver cell apoptosis of mouse acute hepatic failure model obviously so as to increase the survival rate of mouse models of acute hepatic failure.

Analysis of hippocampal gene expression profile of Alzheimer's disease model rats using genome chip bioinformatics

In this study, an Alzheimer＇s disease model was established in rats through stereotactic injection of condensed amyloid beta 1-40 into the bilateral hippocampus, and the changes of gene expression profile in the hippocampus of rat models and sham-operated rats were compared by genome expression profiling analysis. Results showed that the expression of 50 genes was significantly up-regulated （fold change 〉 2）, while 21 genes were significantly down-regulated in the hippocampus of Alzheimer＇s disease model rats （fold change 〈 0.5） compared with the sham-operation group. The differentially expressed genes are involved in many functions, such as brain nerve system development, neuronal differentiation and functional regulation, cellular growth, differentiation and apoptosis, synaptogenesis and plasticity, inflammatory and immune responses, ion channels/transporters, signal transduction, cell material/energy metabolism. Our findings indicate that several genes were abnormally expressed in the metabolic and signal transduction pathways in the hippocampus of amyloid beta 1 40-induced rat model of Alzheimer＇s disease, thereby affecting the hippocampal and brain functions.

Expression Profile of Metastasis-associated Genes in Esophageal Squamous Cell Carcinoma

The differentially expressed genes between esophageal squamous cell carcinoma （ESCC） with or without lymphatic metastasis were investigated by gene chip, and the lymphatic metastasisassociated genes were screened out. Expression array was used to detect the mRNA from both the primary carcinoma and the corresponding esophageal epithelium in 15 cases of human ESCC. The lymphatic metastasis-associated genes were screened by bioinformatics between ESCC with or without lymphatic metastasis. The results showed that 43 （4. 85 %） genes significantly differed between the ESCC with and without lymphatic metastasis （P〈0.05）, of which 18（2.03 %）were upregulated and 25 （2.82 %） down-regulated. The up-regulated genes were involved in cell adhesion molecules and cell membrane receptors and the down-regulated genes were mostly cell cycle regulators and intracellular signaling molecules. It was suggested that lymphatic metastasis-associated genes were screened by gene chip, which was helpful to understand the molecular mechanism of ESCC lymphatic metastasis and lymphatic metastasis-associated genes might be used as diagnostic markers and therapeutic targets for lymphatic metastasis.

Immune Responses to Trichloroethylene and Skin Gene Expression Profiles in Sprague Dawley Rats

Objective To characterize the immune reaction in SD rats exposed to trichloroethylene （TCE） and to identify the gene expression profiles involved in skin after TCE exposure. Methods Fifteen percent of TCE was injected intradermally into the rat back （100 μL/120 g） at intervals of 7 days. Whole blood was collected 24 h after the fifth or seventh intradermic administration of TCE. The percentages of CD4＋ and CD8＋ of T lymphocytes were measured by a flow cytometer. The concentrations of IFN-gamma and 1L-4 in the serum were semi-quantified by ELISA. Total RNAs of skin samples at 3 h or 24 h after the seventh dose of TCE in SD rats were extracted, and gene expression proftles of these tissues were analyszed by rat toxicology U34 array of Affymetrix. Results Obvious decline of CD4＋ in T lymphocytes was observed in the TCE-administer group. No significant concentration differences in IFN-gamma and IL-4 were found between TCE-treated and control rats. Gadd45a and Mel were significantly up regulated in skin tissue 24 h after TCE exposure. The expression regulation of immune response factors was as active as proteins associated with lipid metabolism and synthesis process in these skin samples of SD rats exposed to TCE. Conclusion T-helper type 1 cells mediate immune response can not be elicited in TCE-treated SD rats, but certain immune disorder can be induced.

THE GENE EXPRESSION PROFILE OF HIGHLY METASTATIC HUMAN OVARIAN CANCER CELL LINE BY GENE CHIP

Objective: To study the gene expression of high metastatic human ovarian carcinoma cell line (HO-8910PM) and to screen for novel metastasis- associated genes by cDNA microarray. Methods: The cDNA was retro-transcribed from equal quantity mRNA derived from tissues of highly metastatic ovarian carcinoma cell line and normal ovarian, and was labeled with Cy5 and Cy3 fluorescence as probes. The mixed probes were hybridized with BioDoor 4096 double dot human whole gene chip. The chip was scanned by scanArray 3000 laser scanner. The acquired image was analyzed by ImaGene 3.0 software. Results: By applying the cDNA microarray we found: A total of 323 genes whose expression level were 3 times higher or lower in HO-8910PM cell than normal ovarian epithelium cell were screened out, with 71 higher and 252 lower respectively. Among these 10 were new genes. 67 genes showed expression difference bigger than 6 times between HO-8910PM cell and normal ovarian epithelium cell, among these genes 12 were higher, 55 lower, and two new genes were found. Conclusion: cDNA microarray technique is effective in screening the differentially expressed genes between human ovarian cancer cell line (HO-8910PM) and normal ovarian epithelium cell. Using the cDNA microarray to analyze of human ovarian cancer cell line gene expression profile difference will help the gene diagnosis, treatment and protection.

Analysis of Gene Expression Profile in Lung Adenosquamous Carcinoma Using cDNA Microarray

Gene expression profile of the lung adenosquamous carcinoma was characterized by using cDNA microarray chip containing 4 096 human genes. Among target genes, 508 differentially expressed genes were identified in adenosquamous carcinoma of the lung, 232 genes were overexpressed and 276 genes were underexpressed. Among them, 92 genes are cell signals transduction genes, 34 genes are proto-oncogenes and tumor suppressor genes or cell cycle related genes or cell apoptosis related genes, 29 genes are cell skeleton genes, 28 genes are DNA synthesis, repair and recombination genes, 12 genes are DNA binding and transcription genes. These genes may be associated with the occurence and development of adenosquamous carinome of the lung. Key words lung carcinoma - adenosquamous carcinoma - microarray - gene expression profile CLC number R 734.2 Foundation item: Supported by the National Natural Science Foundation of China (39870305)Biography: YANG Fei(1972-), female, Ph. D candidate, research direction: etiology and pathogenesis of lung cancer.

Gene expression profile analysis reveals the effect of metformin treatment on HepG2 cells

Metformin is a first-line drug in the fight against type 2 diabetes.In recent years,studies have shown that metformin has some preventive and therapeutic effects on liver cancer,but the effects of metformin on the gene expression of liver cancer cells are not fully known.This study focused on the differences in the gene expression profiles in liver cancer cells treated with or without metformin.A total of 153 differentially expressed genes(DEGs)(FC>2 and q-values<0.001)were found,including 77 upregulated genes and 76 downregulated genes.These DEGs are involved in mitogen-activated protein kinase(MAPK),nuclear factor-kappa B(NF-κB),cell adhesion molecules(CAMs),and leukocyte transendothelial migration signaling pathways.These findings reveal the effects of metformin treatment on gene expression profiles in liver cancer cells and provide new clues for unveiling the mechanism of the antitumor effects of metformin.

Difference of Gene Expression Profiles between Barrett’s Esophagus and Cardia Intestinal Metaplasia by Gene Chip

The difference of gene expression profile changes in Barrett＇s esophagus （BE） and cardia intestinal metaplasia （CIM） epithelium was studied and the novel associated genes were screened in the early stage by cDNA microarray. The cDNA retro-transcribed from equal quantity mRNA from BE and CIM epithelial tissues were labeled with Cy3 and Cy5 fluorescence as probes. The mixed probe was hybridized with three pieces BiostarH-40s double dot human whole gene chip. The chips were scanned with a ScanArray 4000. The acquired images were analyzed using GenePix Pro 3.0 software. It was found a total of 141 genes were screened out that exhibited differentially expression more than 2 times in all three chips. It was identified that in gene expression profiles of BE, 74 genes were up-regulated and 67 down-regulated as compared with CIM. The comparison between the difference of gene expression profile changes in BE and CIM epithelia revealed that there existed the difference between BE and CIM at gene level. 141 genes with the expression more than two time were probably related to the occurrence and development of BE and the promotion or progress in adenocarcinoma.

Study on the effect of Baiban Ointment on the whole gene expression profile of vulvar sclerosing lichen

Objective:To study the effect of Baiban Ointment on gene expression profile of vulvar lichen sclerosus before and after treatment,and to provide theoretical basis for the early diagnosis and clinical treatment of vulvar lichen sclerosus.Methods:Nine patients with vulvar lichen sclerosus diagnosed pathologically were selected as the study object,and the Baiban ointment was applied locally for 3 months.Gene chip technology was used to detect the vulva skin tissue in the same area before and after treatment,and the peripheral normal skin was used as the normal control group to analyze the change of gene expression profile.FC and P values were used as research indicators,and the standard of differential expression was FC value≥1.5,P value<0.05.The differentially expressed genes were screened,and GO analysis and KEGG database gene pathway analysis were carried out.Results:Compared with the normal control group,there were 22 differentially expressed genes in Group A before treatment,of which 21 were up-regulated and 1 was down regulated;compared with Group B before treatment,there were 23 differentially expressed genes,all of which were up regulated.The differentially expressed genes related to Baiban ointment treatment included TBK1、STAT1、ITGAM、VCAN、PRKACB、PROM1、PLAT、SERPINA1.Go analysis showed that the up-regulated differential genes were mainly concentrated in the extracellular exosome,cytosol,extracellular space,Golgi apparatus and other cellular component,using the biological process such as positive regulation of cell metabolism,signal transduction regulation,cell adhesion,etc.,to play a role in protein binding,enzyme activity regulation and other molecular functions.KEGG signaling pathway showed that the differentially expressed genes were significantly enriched in Toll like receptor signaling pathway,Nod like receptor signaling pathway and Fc RI signaling pathway,all of which were up-regulated signaling pathways.Among them,Toll like receptor signaling pathway and Nod like receptor signaling pathway are most closely related to the disease.Conclusion:Baiban ointment may play a role in regulating metabolism,inflammation and immune response by regulating the expression of related genes,affecting the signal transduction such as Toll like receptor signal pathway.The pathway and the genes screened in this study will provide a direction for the future study of this disease.

Expression Profile Analysis of Genes Involved in Brassinosteroid Biosynthesis Pathway in Cotton Fiber Development

Cotton(Gossypium hirsutum L.) is the leading fiber crop and one of the mainstays of the economy in the world.Cotton fibers,as the main product of cotton plants,are unicellular,linear

Transcriptional regulatory network during axonal regeneration of dorsal root ganglion neurons:laser-capture microdissection and deep sequencing

The key regulators and regeneration-associated genes involved in axonal regeneration of neurons after injury have not been clarified.In high-throughput sequencing,various factors influence the final sequencing results,including the number and size of cells,the depth of sequencing,and the method of cell separation.There is still a lack of research on the detailed molecular expression profile during the regeneration of dorsal root ganglion neuron axon.In this study,we performed lase r-capture microdissection coupled with RNA sequencing on dorsal root ganglion neurons at 0,3,6,and 12 hours and 1,3,and 7 days after sciatic nerve crush in rats.We identified three stages after dorsal root ganglion injury:early(3-12 hours),pre-regeneration(1 day),and regeneration(3-7 days).Gene expression patterns and related function enrichment res ults showed that one module of genes was highly related to axonal regeneration.We verified the up-regulation of activating transcription factor 3(Atf3),Kruppel like factor 6(Klf6),AT-rich inte raction domain 5A(Arid5α),CAMP responsive element modulator(Crem),and FOS like 1,AP-1 transcription factor Subunit(Fosl1) in dorsal root ganglion neurons after injury.Suppressing these transcription factors(Crem,Arid5o,Fosl1 and Klf6) reduced axonal regrowth in vitro.As the hub transcription factor,Atf3 showed higher expression and activity at the preregeneration and regeneration stages.G protein-coupled estrogen receptor 1(Gper1),inte rleukin 12a(Il12α),estrogen receptor 1(ESR1),and interleukin 6(IL6) may be upstream factors that trigger the activation of Atf3 during the repair of axon injury in the early stage.Our study presents the detailed molecular expression profile during axonal regeneration of dorsal root ganglion neurons after peripheral nerve injury.These findings may provide reference for the clinical screening of molecular targets for the treatment of peripheral nerve injury.

Hippocampal gene expression in a rat model of depression after electroacupuncture at the Baihui and Yintang acupoints

Preliminary basic research and clinical findings have demonstrated that electroacupuncture ther- apy exhibits positive effects in ameliorating depression. However, most studies of the underlying mechanism are at the single gene level; there are few reports regarding the mechanism at the whole-genome level. Using a rat genomic gene-chip, we profiled hippocampal gene expression changes in rats after electroacupuncture therapy. Electroacupuncture therapy alleviated depres- sion-related manifestations in the model rats. Using gene-chip analysis, we demonstrated that electroacupuncture at Baihui （DU20） and Yintang （EX-HN3） regulates the expression of 21 genes. Real-time PCR showed that the genes Vgf, lgf2, Trnp32, Loc500373, Hifla, Folrl, Nrnb, and Rtn were upregulated or downregulated in depression and that their expression tended to nor- malize after electroacupuncture therapy. These results indicate that electroacupuncture at Baihui and Yintang modulates depression by regulating the expression of particular genes.

Genome-wide identification,phylogeny and expression analysis of the SBP-box gene family in maize(Zea mays)

The SQUAMOSA promoter binding protein （SBP）-box genes encode a kind of plant-specific transcription factors （TFs） and play important roles in the regulation of plant development. In this study, a genome-wide characterization of this family was conducted in maize （Zea mays）. Thirty-one SBP-box genes were identified to be distributed in nine chromosomes and 16 of them were complementary to the mature ZmmiR156 sequences. All the Z. mays SBP （ZmSBP） genes were classified into two clusters with eight subgroups according to the phylogenetic analysis of proteins, which were consistent with the pattern of exon-intron structures. The phylogenetic tree of the ZmSBP, Oryza sativa SBP-like （OsSPL） and Arabidopsis thaliana SBP-like （AtSPL） genes were constructed and all the SBP-box genes were divided into eight groups, which was the same as the classification of ZmSBP genes. The comparision of the expression profiles of all SBP-box genes in these three species indicated that most orthologous genes had similar expression patterns. The results from this study provided a basic understanding of the ZmSBP genes and might facilitate future researches for elucidating the SBP-box genes function in maize.

Methotrexate combined with methylprednisolone for the recovery of motor function and differential gene expression in rats with spinal cord injury

Methylprednisolone is a commonly used drug for the treatment of spinal cord injury, but high doses of methylprednisolone can increase the incidence of infectious diseases. Methotrexate has anti-inflammatory activity and immunosuppressive effects, and can reduce in- flammation after spinal cord injury. To analyze gene expression changes and the molecular mechanism of methotrexate combined with methylprednisolone in the treatment of spinal cord injury, a rat model of spinal cord contusion was prepared using the PinPointTM preci- sion cortical impactor technique. Rats were injected with methylprednisolone 30 mg/kg 30 minutes after injury, and then subcutaneously injected with 0.3 mg/kg methotrexate 1 day after injury, once a day, for 2 weeks. TreadScan gait analysis found that at 4 and 8 weeks after injury, methotrexate combined with methylprednisolone significantly improved hind limb swing time, stride time, minimum longitudinal deviation, instant speed, footprint area and regularity index. Solexa high-throughput sequencing was used to analyze differential gene ex- pression. Compared with methylprednisolone alone, differential expression of 316 genes was detected in injured spinal cord treated with methotrexate and methylprednisolone. The 275 up-regulated genes were mainly related to nerve recovery, anti-oxidative, anti-inflammatory and anti-apoptotic functions, while 41 down-regulated genes were mainly related to proinflammatory and pro-apoptotic functions. These results indicate that methotrexate combined with methylprednisolone exhibited better effects on inhibiting the activity of inflammatory cytokines and enhancing antioxidant and anti-apoptotic effects and thereby produced stronger neuroprotective effects than methotrexate alone. The 316 differentially expressed genes play an important role in the above processes.

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.

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.