Genome wide investigation of Hsf gene family in Phoebe bournei:identification,evolution,and expression after abiotic stresses

Heat shock transcription factors(Hsfs)have important roles during plant growth and development and responses to abiotic stresses.The identification and func-tion of Hsf genes have been thoroughly studied in various herbaceous plant species,but not woody species,especially Phoebe bournei,an endangered,unique species in China.In this study,17 members of the Hsf gene family were identi-fied from P.bournei using bioinformatic methods.Phyloge-netic analysis indicated that PbHsf genes were grouped into three subfamilies:A,B,and C.Conserved motifs,three-dimensional structure,and physicochemical properties of the PbHsf proteins were also analyzed.The structure of the PbHsf genes varied in the number of exons and introns.Pre-diction of cis-acting elements in the promoter region indi-cated that PbHsf genes are likely involved in responses to plant hormones and stresses.A collinearity analysis dem-onstrated that expansions of the PbHsf gene family mainly take place via segmental duplication.The expression levels of PbHsf genes varied across different plant tissues.On the basis of the expression profiles of five representative PbHsf genes during heat,cold,salt,and drought stress,PbHsf pro-teins seem to have multiple functions depending on the type of abiotic stress.This systematic,genome-wide investigation of PbHsf genes in P.bournei and their expression patterns provides valuable insights and information for further func-tional dissection of Hsf proteins in this endangered,unique species.

Genome-wide identification,characterization and functional prediction of the SRS gene family in sesame(Sesamum indicum L.)

Sesame(Sesamum indicum L.)is an ancient oilseed crop of the Pedaliaceae family with high oil content and potential health benefits.SHI RELATED SEQUENCE(SRS)proteins are the transcription factors(TFs)specific to plants that contain RING-like zinc finger domain and are associated with the regulation of several physiological and biochemical processes.They also play vital roles in plant growth and development such as root formation,leaf development,floral development,hormone biosynthesis,signal transduction,and biotic and abiotic stress responses.Nevertheless,the SRS gene family was not reported in sesame yet.In this study,identification,molecular characterization,phylogenetic relationship,cis-acting regulatory elements,protein-protein interaction,syntenic relationship,duplication events and expression pattern of SRS genes were analyzed in S.indicum.We identified total six SiSRS genes on seven different linkage groups in the S.indicum genome by comparing with the other species,including the model plant Arabidopsis thaliana.The SiSRS genes showed variation in their structure like2–5 exons and 1–4 introns.Like other species,SiSRS proteins also contained‘RING-like zinc finger'and‘LRP1'domains.Then,the SiSRS genes were clustered into subclasses via phylogenetic analysis with proteins of S.indicum,A.thaliana,and some other plant species.The cis-acting regulatory elements analysis revealed that the promoter region of SiSRS4(SIN_1011561)showed the highest 13 and 16 elements for light-and phytohormone-responses whereas,SiSRS1(SIN_1015187)showed the highest 15 elements for stress-response.The ABREs,or ABA-responsive elements,were found in a maximum of 8 copies in the SiSRS3(SIN 1009100).Moreover,the available RNA-seq based expression of SiSRS genes revealed variation in expression patterns between stress-treated and non-treated samples,especially in drought and salinity conditions in.S.indicum.Two SiSRS genes like SiSRS1(SIN_1015187)and SiSRS5(SIN_1021065),also exhibited variable expression patterns between control vs PEG-treated sesame root samples and three SiSRS genes,including SiSRS1(SIN_1015187),SiSRS2(SIN_1003328)and SiSRS5(SIN_1021065)were responsive to salinity treatments.The present outcomes will encourage more research into the gene expression and functionality analysis of SiSRS genes in S.indicum and other related species.

Transcriptome profiling and RXR gene family identification reveals the molecular mechanism of rapid aging after spawning of cuttlefish Sepiella japonica

Sepiella japonica is a worldwide marine cuttlefish species of high economic value.S.japonica routinely modifying behaviors in reproductive life,such as rapid aging until death after spawning,has been recognized in artificial breeding.However,reproductive behavior at the level of genes is rarely reported,thus,the research on the genetic basis of behavior,reproduction,and artificial breeding was limited.We applied RNA-seq in different stages of reproduction to investigate the reason of rapid aging after spawning,pre-maturity,pre-spawning after maturity,and post-spawning.The retinoid X receptor(RXR)gene family in S.japonica was identified,and 1343–1452 differentially expressed genes(DEGs)in all 3 stages of reproductive life were identified from pairwise m RNA comparisons.Furthermore,through the GO term and KEGG analysis,S.japonica could handle neuronal development and network formation before maturity and have a functional degradation of neural communication,signal transduction,vision,and gene expression after spawning.Eight Sj RXRαs have been identified and they played different roles in growth development or reproduction.Therefore,the regulation of several channels and receptors is the intrinsic molecular mechanism of rapid aging after spawning in S.japonica.This study revealed the survival strategy and provided fundamental data on the level of genes for understanding the reproductive behavior and the reproduction of S.japonica.

Genome-Wide Identification of Tomato (Solanum lycopersicum L.) CKX Gene Family and Expression Analysis in the Callus Tissue under Zeatin Treatment

The cytokinin oxidase/dehydrogenase(CKX)enzyme is essential for controlling thefluctuating levels of endogen-ous cytokinin(CK)and has a significant impact on different aspects of plant growth and development.Nonethe-less,there is limited knowledge about CKX genes in tomato(Solanum lycopersicum L.).Here we performed genome-wide identification and analysis of nine SlCKX family members in tomatoes using bioinformatics tools.The results revealed that nine SlCKX genes were unevenly distributed onfive chromosomes(Chr.1,Chr.4,Chr.8,Chr.10,and Chr.12).The amino acid length,isoelectric points,and molecular weight of the nine SlCKX proteins ranged from 453 to 553,5.77 to 8.59,and 51.661 to 62.494 kD,respectively.Subcellular localization analysis indi-cated that SlCKX2 proteins were located in both the vacuole and cytoplasmic matrix;SlCKX3 and SlCKX5 pro-teins were located in the vacuole;and SlCKX1,4,6,7,8,and 9 proteins were located in the cytoplasmic matrix.Furthermore,we observed differences in the gene structures and phylogenetic relationships of SlCKX proteins among different members.SlCKX1-9 were positioned on two out of the three branches of the CKX phylogenetic tree in the multispecies phylogenetic tree construction,revealing their strong conservation within phylogenetic subgroups.Unique patterns of expression of CKX genes were noticed in callus cultures exposed to varying con-centrations of exogenous ZT,suggesting their roles in specific developmental and physiological functions in the regeneration system.These results may facilitate subsequent functional analysis of SlCKX genes and provide valu-able insights for establishing an efficient regeneration system for tomatoes.

Identification and Molecular Characterization of the Alkaloid Biosynthesis Gene Family in Dendrobium catenatum

As one of the main active components of Dendrobium catenatum, alkaloids have high medicinal value. The physicochemicalproperties, conserved domains and motifs, phylogenetic analysis, and cis-acting elements of the genefamily members in the alkaloid biosynthesis pathway of D. catenatum were analyzed by bioinformatics, and theexpression of the genes in different years and tissues was analyzed by qRT-PCR. There are 16 gene families,including 25 genes, in the D. catenatum alkaloid biosynthesis pathway. The analysis of conserved domains andmotifs showed that the types, quantities, and orders of domains and motifs were similar among members ofthe same family, but there were significant differences among families. Phylogenetic analysis indicated that thegene family members showed some evolutionary conservation. Cis-acting element analysis revealed that therewere a large number of light-responsive elements and MYB (v-myb avian myeloblastosis viral oncogene homolog)-related elements in these genes. qRT-PCR showed that expressions of gene family members involved in alkaloidsynthesis were different in different years and tissues of D. catenatum. This study provides a theoretical basisfor further exploration of the regulatory mechanisms of these genes in the alkaloid biosynthesis of D. catenatum.

Identification and Analysis of the WRKY Transcription Factor Gene Family in Verbena bonariensis

The WRKY transcription factor gene family is one of the unique gene families in plants.It plays an important role in response to abiotic stresses such as cold and drought,hormone signal transduction,regulation of biosynthesis,leaf senescence seed germination,etc.However,little information is available about WRKY transcription factors in Verbena bonariensis.In this study,70 VbWRKY genes were identified from the whole genome.The phylogenetic analysis of the WRKY gene family in V.bonariensis and Arabidopsis shows that the WRKY genes in V.bonariensis can be divided into three groups:I,II,and III,which contain 13,47,and 10 members,respectively.Group II can be further divided into five subclasses:IIa(5),IIb(10),IIc(18),IId(6),and IIe(8).Conservative motif analysis showed that 64 proteins encoded by the VbWRKY gene had conserved motifs 1,2 and 3,and the same subclass motif elements were approximately the same.The collinearity analysis showed that there were 44 homologous gene pairs among the VbWRKYs,and these homologous gene pairs may have the same function.Promoter sequence analysis showed that the VbWRKY gene has multiple cis-acting elements,including not only cis-acting elements related to low-temperature and light responses,but also cis-acting elements related to hormone regulation,Among them,most VbWRKY genes contain response elements about low-temperature,and 30 VbWRKY genes contain low-temperature response elements(LTR),and 61 VbWRKY genes contain abscisic acid response elements(ABRE),indicating that VbWRKY plays a crucial role in plant growth and abiotic stress.According to the expression of VbWRKY in the cold stress and different tissues transcriptome,70 VbWRKY genes played their respective roles in various tissues and stages to regulate plant growth,Also,some of them participated in the process of cold stress tolerance,52 VbWRKYs showed significant differences in expression under cold stress,and 37 VbWRKY genes were up-regulated under cold stress.9 VbWRKY genes were selected for quantitative real-time PCR(qRT-PCR)analysis under low-temperature stress,and the results showed that all 9 genes were upregulated under low-temperature stress.Ultimately,the present study provides a comprehensive analysis of the predicted V.bonariensis WRKY genes family,which provided a theoretical basis for the study of low-temperature resistance and growth and development of V.bonariensis.

Molecular Composition and Evolution of the Chalcone Synthase (CHS) Gene Family in Five Species of Camellia (Theaceae)

The molecular composition and evolution of the chalcone synthase (CHS) gene family from five species in Camellia (Theaceae) are explored in this study. Sixteen CHS exon 2 from four Camellia species were amplified from total DNA by PCR method. Three sequences of the fifth species in Camellia and two sequences of Glycine max as the designated outgroups were obtained from GenBank. Our results indicated that CHS gene family in Camellia was differentiated to three subfamilies (A, B, C) during the evolutionary history with six groups (A1, A2, A3, BI, B2, C). Among them, only group A2 was possessed by all five species in this study. However, the other five groups were detected only in some species of the plants studied. All members of CHS gene family in this study had high sequence similarity, more than 90% among the members in the same subfamily and more than 78% among different subfamilies at nucleotide level., According to the estimated components of amino acids, the function of CHS genes in Camellia had been diverged. The nucleotide substitutions of the different groups were not identical. Based on phylogenetic analyse inferred from sequences of CHS genes and their deduced amino acid sequences, we concluded that the CHS genes with new function in this genus were evolved either by mutations on several important sites or by accumulation of the mutations after the gene duplication. A further analysis showed that the diversification of CHS genes in Camellia still occurred recently, and the evolutionary models were different to some extant among different species. So we assumed that the different evolutionary models resulted from the impacts of variable environmental elements after the events of speciation.

Expansion of the Actin Gene Family in Amphioxus

Actins are a small family of ubiquitous proteins that are essential cytoskeletal components and are highly conserved during evolution. Actins are usually divided into two classes, the cytoplasmic and muscle actins, which have different functional roles. Here we systematically analyzed the actin genes in the genome of the primitive chordate amphioxus （Branchiostoma floridae）. We found that amphioxus contains more than 30 actin genes, many of which are linked. Phylogenetic analysis suggests the amphioxus actin genes have clearly undergone extensive expansion through tandem duplications. The actin genes＇ structure also varies a lot, containing 2 to 7 exons. We also cloned two muscle type of actin genes from the amphioxus （B. belcheri） and compared their expression patterns during early development. The slight difference in their expression suggests functional diversification of these actin genes. Our results shed light on the evolution both of actin genes themselves and their functional roles in chordate development.

In silico genome-wide identification,phylogeny and expression analysis of the R2R3-MYB gene family in Medicago truncatula

The R2R3-MYB genes make up one of the largest transcription factor families in plants, and play regulatory roles in various biological processes such as development, metabolism and defense response. Although genome-wide analyses of this gene family have been conducted in several species, R2R3-MYB genes have not been systematically analyzed in Medicago truncatula, a sequenced model legume plant. Here, we performed a comprehensive, genome-wide computational analysis of the structural characteristics, phylogeny, functions and expression patterns of M. truncatula R2R3-MYB genes. DNA binding domains are highly conserved among the 155 putative MtR2R3-MYB proteins that we identified. Chromosomal location analysis revealed that these genes were distributed across all eight chromosomes. Results showed that the expansion of the MtR2R3-MYB family was mainly attributable to segmental duplication and tandem duplication. A comprehensive classification was performed based on phylogenetic analysis of the R2R3-MYB gene families in M. truncatula, Arabidopsis thaliana and other plant species. Evolutionary relationships within clades were supported by clade-specific conserved motifs outside the MYB domain. Species-specific clades have been gained or lost during evolution, resulting in functional divergence. Also, tissue-specific expression patterns were investigated. The functions of stress response-related clades were further verified by the changes in transcript levels of representative R2R3-MYB genes upon treatment with abiotic and biotic stresses. This study is the first report on identification and characterization of R2R3-MYB gene family based on the genome of M. truncatula, and will facilitate functional analysis of this gene family in the future.

Genome-wide identification and characterization of the bHLH gene family in an ornamental woody plant Prunus mume

The basic helix-loop-helix(bHLH)transcription factor family is the second-largest family in plants,where it plays essential roles in development,and the responses to multiple abiotic and biotic stressors.However,little information is available about this gene family in Prunus mume,which is widely cultivated in East Asia as an ornamental fruit tree.Here,100 PmbHLH genes were identified,and their evolution and functions were explored in P.mume for the first time.The PmbHLH genes were classified into 21 subfamilies.The chromosomal distribution,physicochemical properties,bHLH domain,conserved motif,and intron/exon compositions were also analyzed.Furthermore,the evolutionary pattern,divergence time of the PmbHLH family,and genetic relationships among P.mume,Arabidopsis thaliana,and Prunus persica and Fragaria vesca of Rosaceae were explored.The functional prediction analysis of these PmbHLHs indicated that their functions varied,and included participating in the formation of organs and tissues,responding to stress,and the biosynthesis and metabolism of hormones and other secondary metabolites.Interestingly,expression analyses of PmbHLHs also revealed diverse expression patterns.Most of the PmbHLH genes were highly expressed in roots and stems,and a few were highly expressed in leaves,buds,and fruits,indicating tissue expression specificity.Eight PmbHLH genes,which were upregulated during low-temperature stress,may have critical roles in the response to cold stress.Ten PmbHLHs were differentially expressed between weeping and upright branches in a P.mume F_(1) population.These results shed light on the structure and evolution of the PmbHLH gene family,and lay a foundation for further functional studies of the bHLH genes.

Genome-wide analysis of the calcium-dependent protein kinase gene family in Gossypium raimondii

Plant calcium-dependent protein kinases （CDPKs） play important roles in diverse physiological processes by regulating the downstream components of calcium signaling. To date, only a few species of the plant CDPK gene family have been functionally identified. In addition, there has been no systematic analysis of the CDPK family in cotton. Here, 41 putative cotton CDPK （GrCDPK） genes were identified via bioinformatics analysis of the entire genome of Gossypium raimondii and were classified into four groups based on evolutionary relatedness. Gene structure analysis indicated that most of these GrCDPK genes share a similar intron-exon structure （7 or 8 exons）, strongly supporting their close evolutionary relationships. Chromosomal distributions and phylogenetics analysis showed that 13 pairs of GrCDPK genes arose via segmental duplication events. Furthermore, using microarray data of upland cotton （G. hirsutum L.）, comparative profiles analysis of these GhCDPKs indicated that some of the encoding genes might be involved in the responses to multiple abiotic stresses and play important regulatory roles during cotton fiber development. This study is the first genome-wide analysis of the CDPK family in cotton, and it will provide valuable information for the further functional characterization of cotton CDPK genes.

Comprehensive identification and analyses of the Hsf gene family in the whole-genome of three Apiaceae species

Apiaceae is a major family from Apiales and includes many important vegetable and medicinal crops.Heat shock transcription factors(Hsf)play important roles in heat tolerance during plant development.Here,we conducted systematic analyses of the Hsf gene family in three Apiaceae species,including 17 Apium graveolens(celery),32 Coriandrum sativum(coriander),and 14 Daucus carota(carrot).A total of 73 Hsf genes were identified in three representative species,including Arabidopsis thaliana,Vitis vinifera,and Lactuca sativa.Whole-genome duplication played important roles in the Hsf gene family’s expansion within Apiaceae.Interestingly,we found that coriander had more Hsf genes than celery and carrot due to greater expansion and fewer losses.Twenty-seven branches of the phylogenetic tree underwent considerable positive selection in these Apiaceae species.We also explored the expression patterns of Hsf genes in three plant organs.Collectively,this study will serve as a rich gene resource for exploring the molecular mechanisms of heat tolerance.Additionally,this is the first study to report on the Hsf gene family in Apiaceae;thus,our research will provide guidance for future comparative and functional genomic studies on the Hsf gene family and others in Apiaceae.

The tyrosinase gene family and albinism in fish

Tyrosinase exists universally in organisms and is a characteristic enzyme of melanocytes. Tyrosinase family genes in vertebrates consist of 3 related members; tyrosinase （TYR, Tyr）, tyro sinase-related protein-1 （TRP-1, Tyrp 1）, and tyro sinase-related protein-2 （TRP-2, Tyrp2, Dct）. These proteins catalyze melanin biosynthesis in pigment cells and play important roles in determining vertebrate coloration. Transcription of the TYR and TRP genes is useful for studying neural crest and optic vesicle cell migration and differentiation during emblyogenesis and important in pigment rescue in fish. In this paper, the structure of gene and protein molecular evolution, function and roles of the TYR family in fish were reviewed.

Genome-wide identification of the MATE gene family and functional characterization of PbrMATE9 related to anthocyanin in pear

Plant multidrug and toxic compound extrusion(MATE) genes play an important role in the process of detoxification, plant morphogenesis, and anthocyanin accumulation. However, whether the MATE gene family functions in pear peel coloration is still unknown. To evaluate and identify the MATE gene family members which are involving in anthocyanin accumulation and coloration in pear. In this study, 85 MATE genes were identified in the reference pear genome of ‘Dangshansuli’ through genome-wide identification. Based on gene structure and phylogenetic tree analysis, the MATE family was divided into five subfamilies. RNA sequencing and quantitative real-time polymerase chain reaction(qRTPCR) indicated that the expression patterns of PbrMATEs were tissue-specific. 28.24%(24) of PbrMATE genes were expressed in the fruits, and44.71%(38) of PbrMATE genes were expressed in the leaves. Additionally, we found that the expression levels of PbrMATE9, PbrMATE26,PbrMATE50, and PbrMATE69 in debagged fruits with red peel were significantly higher than those in bagged fruits without red peel, according to our bagging/debagging treatment of ‘Mantianhong’. The expression pattern of PbrMATE9 was consistent with the variation trend in anthocyanin content, suggesting that it might play an important role in anthocyanin accumulation in response to light exposure. Subcellular localization showed that PbrMATE9 was a membrane protein. More strikingly, the transient overexpression of PbrMATE9 promoted anthocyanin accumulation in the peel of pear, and the expression of structural genes(PbrCHI, PbrANS, PbrDFR, and PbrUFGT) in the anthocyanin biosynthesis pathway also increased significantly. Through co-expression network analysis, the transcription factors were identified, such as WRKY, COL,GATA, and BBX, which might be involved in the regulation of PbrMATE9. The study has enriched the genetic resources and improved the understanding of the regulation network of anthocyanin accumulation in pear.

Genome-wide identification and characterization of the JAZ gene family and its expression patterns under various abiotic stresses in Sorghum bicolor

The jasmonate ZIM domain(JAZ)protein belongs to the TIFY((TIF[F/Y]XG)domain protein)family,which is composed of several plant-specific proteins that play important roles in plant growth,development,and defense responses.However,the mechanism of the sorghum JAZ family in response to abiotic stress remains unclear.In the present study,a total of 17 JAZ genes were identified in sorghum using a Hidden Markov Model search.In addition,real-time quantification polymerase chain reaction(RT-qPCR)was used to analyze the gene expression patterns under abiotic stress.Based on phylogenetic tree analysis,the sorghum JAZ proteins were mainly divided into nine subfamilies.A promoter analysis revealed that the SbJAZ family contains diverse types of promoter cis-acting elements,indicating that JAZ proteins function in multiple pathways upon stress stimulation in plants.According to RT-qPCR,SbJAZ gene expression is tissuespecific.Additionally,under cold,hot,polyethylene glycol,jasmonic acid,abscisic acid,and gibberellin treatments,the expression patterns of SbJAZ genes were distinctly different,indicating that the expression of SbJAZ genes may be coordinated with different stresses.Furthermore,the overexpression of SbJAZ1 in Escherichia coli was found to promote the growth of recombinant cells under abiotic stresses,such as PEG 6000,NaCl,and 40℃ treatments.Altogether,our findings help us to better understand the potential molecular mechanisms of the SbJAZ family in sorghum in response to abiotic stresses.

Bioinformatics analysis of structure and function in the MRP gene family and its expression in response to various drugs in Bursaphelenchus xylophilus

Genes homologous to members of the MRP gene family in Caenorhabditis elegans are important in drug resistance.To further explore the molecular mechanism of drug resistance in pine wood nematode(Bursaphelenchus xylophilus),we used bioinformatics approaches to analyze genomic data for B.xylophilus and identified Bx-MRP genes.We predicted the structure and function of the genes and encoded proteins.Using bioinformatics programs to predict and analyze various properties of the predicted proteins,including hydrophobicity,transmembrane regions,phosphorylation sites,and topologically isomeric structures,of these Bx-MRP genes,we determined that they function in transmembrane transport.From the results of RT-qPCR,the Bx-MRP family members confer significant differential resistance to different drug treatments.After treatment with different concentrations of emamectin benzoate,avermectin and matrine,the expression of each gene increased with increasing drug concentrations,indicating that the family members play a positive role in the regulation of multidrug resistance.

Genome-wide identification and expression analysis of NtbHLH gene family in tobacco(Nicotiana tabacum)and the role of NtbHLH86 in drought adaptation

The bHLH transcription factors play pivotal roles in plant growth and development,production of secondary metabolites and responses to various environmental stresses.Although the bHLH genes have been well studied in model plant species,a comprehensive investigation of the bHLH genes is required for tobacco with newly obtained high-quality genome.In the present study,a total of 309 NtbHLH genes were identified and can be divided into 23 subfamilies.The conserved amino acids which are essential for their function were predicted for the NtbHLH proteins.Moreover,the NtbHLH genes were conserved during evolution through analyzing the gene structures and conserved motifs.A total of 265 NtbHLH genes were localized in the 24 tobacco chromosomes while the remained 44 NtbHLH genes were mapped to the scaffolds due to the complexity of tobacco genome.Moreover,transcripts of NtbHLH genes were obviously tissue-specific expressed from the gene-chip data from 23 tobacco tissues,and expressions of 20 random selected NtbHLH genes were further confirmed by quantitative real-time PCR,indicating their potential functions in the plant growth and development.Importantly,overexpressed NtbHLH86 gene confers improve drought tolerance in tobacco indicating that it might be involved in the regulation of drought stress.Therefore,our findings here provide a valuable information on the characterization of NtbHLH genes and further investigation of their functions in tobacco.

Prognostic and biological role of the N-Myc downstream-regulated gene family in hepatocellular carcinoma

BACKGROUND The N-Myc downstream-regulated gene(NDRG)family is comprised of four members(NDRG1-4)involved in various important biological processes.However,there is no systematic evaluation of the prognostic of the NDRG family in hepatocellular carcinoma(HCC).AIM To analyze comprehensively the biological role of the NDRG family in HCC.METHODS The NDRG family expression was explored using The Cancer Genome Atlas.DNA methylation interactive visualization database was used for methylation analysis of the NDRG family.The NDRG family genomic alteration was assessed using the cBioPortal.Single-sample Gene Set Enrichment Analysis was used to determine the degree of immune cell infiltration in tumors.RESULTS NDRG1 and NDRG3 were up-regulated in HCC,while NDRG2 was down-regulated.Consistent with expression patterns,high expression of NDRG1 and NDRG3 was associated with poor survival outcomes(P<0.05).High expression of NDRG2 was associated with favorable survival(P<0.005).An NDRG-based signature that statistically stratified the prognosis of the patients was constructed.The percentage of genetic alterations in the NDRG family varied from 0.3%to 11.0%,and the NDRG1 mutation rate was the highest.NDRG 1-3 expression was associated with various types of infiltrated immune cells.Gene ontology analysis revealed that organic acid catabolism was the most important biological process related to the NDRG family.Gene Set Enrichment Analysis showed that metabolic,proliferation,and immune-related gene sets were enriched during NDRG1 and NDRG3 high expression and NDRG2 low expression.CONCLUSION Overexpression of NDRG1 and NDRG3 and down-expression of NDRG2 are correlated with poor overall HCC prognosis.Our results may provide new insights into the indispensable role of NDRG1,2,and 3 in the development of HCC and guide a promising new strategy for treating HCC.

Identification and expression of the CEP gene family in apple (Malus×domestica)

Plant peptide hormones play important roles in plant growth and development. Among these hormones, the C-TERMINALLYENCODED PEPTIDE（CEP） belongs to a newly found peptide family that regulates root development in Arabidopsis as well as in other species. However, nothing is known about the CEP genes in apple（Malus×domestica, MdCEP）. In this study, a total of 27 apple CEP genes were identified through a genome-wide analysis and were phylogenetically divided into three classes（Ⅰ, Ⅱ and Ⅲ）. The predicted MdCEP genes were distributed across 10 of 17 chromosomes with different densities. Next, the gene structures and motif compositions of the MdCEP genes were analyzed. Subsequently, the expression analysis suggested that the MdCEP genes were highly activated in roots and that MdCEP23 may play an important role in regulating the growth and development of roots. Moreover, all of the MdCEP genes were responsive to multiple abiotic stresses, indicating that MdCEP genes may be involved with various aspects of physiological processes in apple. Nearly one-third of MdCEP genes had a significant response to low nitrogen treatment. Most of the MdCEP genes were up-regulated under stress, including mannitol, polyethylene glycol（PEG） and abscisic acid（ABA）, suggesting that MdCEP genes may be involved in the drought stress response. This study provides insight into the putative functions of the MdCEP genes using a genome-wide analysis of the CEP gene family.

Genome-wide Identification and Analysis of MVD Gene Family in Euphorbiaceae Plants

The mevalonate diphosphate deearboxylase （MVD） is an essential enzyme in mevalonate （MVA） pathway that catalyzes the irreversible Mg2＋ -ATP de- pendent decarboxylation of 6-carben compound mevalonate-5-diphosphate （MVAPP） into 5-carbon isopentenyl diphosphate （ IPP）, the building block of sterol and isoprenoid biosynthesis. In this study, based on the published geanme sequences and ESTs, a genome-wide search was carried out for the first time to identify MVD gene family in four genome-sequenced Euphorbiaceae plants, i.e. castor bean （ Ricinus communis）, physic nut （ Jatropha curcas）, cassava （Manihot esculenta） and rubber tree （Hevea brasiliensis）, and to analyze the gene structure and phylogenetic characteristics. According to the experimental results, 1, 1,2 and 2 MVD genes, which all contain 9 introns, were identh＇ied from castor bean, physic nut, cassava and rubber tree, respectively. Homology analysis indicates that MVD genes are widely distributed in eukaryotes, some archaea and eubacteria, which suggests an early origin of this gerte family. Although MVD genes were identified in most green plants, no homologous genes were found in unicellular green algae. In most genome-sequenced plants including castor bean and physic nut, a single copy of MVD gene was found; however, in cassava and rubber tree, two copies were identified just like that in moss, maize, Arabidopsis and poplar. ＂In castor bean, digital expression profiling suggests that in five examined tissues, i.e. leaf, flower, II/III stage endosperm, V/VI stage endosperm and seed, RcPMK was expressed strongly in flower and II/III stage endosperm, moderately in V/VI stage endosperm and leaf, and weakly in seed.