Identification of key genes underlying clinical features of hepatocellular carcinoma based on weighted gene co‑expression network analysis and bioinformatics analysis

Objective: To identify module genes that are closely related to clinical features of hepatocellular carcinoma (HCC) by weighted gene co‑expression network analysis, and to provide a reference for early clinical diagnosis and treatment. Methods: GSE84598 chip data were downloaded from the GEO database, and module genes closely related to the clinical features of HCC were extracted by comprehensive weighted gene co‑expression network analysis. Hub genes were identified through protein interaction network analysis by the maximum clique centrality (MCC) algorithm;Finally, the expression of hub genes was validated by TCGA database and the Kaplan Meier plotter online database was used to evaluate the prognostic relationship between hub genes and HCC patients. Results: By comparing the gene expression data between HCC tissue samples and normal liver tissue samples, a total of 6 262 differentially expressed genes were obtained, of which 2 207 were upregulated and 4 055 were downregulated. Weighted gene co‑expression network analysis was applied to identify 120 genes of key modules. By intersecting with the differentially expressed genes, 115 candidate hub genes were obtained. The results of enrichment analysis showed that the candidate hub genes were closely related to cell mitosis, p53 signaling pathway and so on. Further application of the MCC algorithm to the protein interaction network of 115 candidate hub genes identified five hub genes, namely NUF2, RRM2, UBE2C, CDC20 and MAD2L1. Validation of hub genes by TCGA database revealed that all five hub genes were significantly upregulated in HCC tissues compared to normal liver tissues;Moreover, survival analysis revealed that high expression of hub genes was closely associated with poor prognosis in HCC patients. Conclusions: This study identifies five hub genes by combining multiple databases, which may provide directions for the clinical diagnosis and treatment of HCC.

Identification and expression analysis of group Ⅲ WRKY transcription factors in cotton

The WRKY proteins constitute a large family of transcription factors in plants containing highly conserved WRKYGQK sequences and zinc-finger-like motifs. To comprehensively study WRKY III genes in cotton, we analyzed the genome sequences of Gossypium hirsutum, G. raimondii and G. arboreum. According to the three genome sequences, 18 group III Gh WRKY genes were identified in G. hirsutum, 12 both in G. raimondii and G. arboreum. Phylogenetic and motif analysis showed that proteins with high similarities could be clustered together and had the same motif components. The ratios of non-synonymous（Ka） to synonymous（Ks） of the Gh WRKY to Gr WRKY or Ga WRKY were lower than 1, which indicated that group III WRKY genes in Gossypium species are under purifying selection. Expression analysis revealed that group III Gh WRKY genes expressed during fiber development and leaf senescence, and most of them could be induced by salicylic acid（SA）, jasmonic acid（JA）, ethylene, abscisic acid（ABA）, mannitol, and Na Cl both in roots and cotyledons. Our study gives a briefly introduction on cotton group III WRKY genes and implicates their potential function in cotton fiber development, leaf senescence and abiotic stresses.

Transcriptome-based identification and expression analysis of the glutathione S-transferase(GST)family in tree peony reveals a likely role in anthocyanin transport

For red-flowered cultivars of tree peony(Paeonia suffruticosa),anthocyanin content is a critical factor determining the different petal pigmentations.Glutathione S-transferases(GSTs)are ubiquitous and multifunctional conjugating proteins that may be responsible for the transport of anthocyanin pigments from the cytoplasm to vacuole.The underlying function of the GST family in tree peony,however,remains unclear.In this study,we systematically isolated and identified a total of 54 putative full-length Ps GST genes through a combination of bioinformatics approaches from transcriptome databases.Intraspecific phylogenetic analyses revealed extensive differentiation in their coding sequences and divided them into 10 of the 14 known classes of plant GSTs.The phylogenetic relationships,evolutionary characteristics,protein domain,and motif organization were clearly conserved among the different phylogenetic subclasses.The results of the RNA-seq and quantitative real-time polymerase chain reaction experiments exhibited extensive variation in gene expression profiles among different developmental stages and varieties.Furthermore,the phylogenetic relationships,expression profiles,protein interactions,weighted gene co-expression network analysis,and correlation analysis results suggested that PsGSTF3(Unigene 0064200)is a candidate participant in anthocyanin transport and the promotion of pigment accumulation,exhibiting a strong positive correlation with anthocyanin content among different tissues(r=0.908**)and an increasing rate of anthocyanin content during the flower developmental process(r=0.961*).These results furthered our understanding of the transport and accumulation functions of the GST family as well as the enhancement of tree peony breeding through molecular biology techniques.

Genome-wide Identification,Classification and Expression Analysis of Lhc Supergene Family in Castor Bean(Ricinus communis L.)

Light-harvesting chlorophyll a/b-binding （LHC） proteins are a group of nuclear-encoded thylakoid proteins that play a key role in plant photosynthesis and are widely involved in light harvesting, energy transfer to the reaction center, maintenance of thylakoid membrane structure, photoprotection and response to en- vironmental conditions, etc. Although/dw supergene family is well characterized in model plants such as Arabidopsis, rice and poplar, little information is available in castor bean （Ricinus communis L. ）. In this study, a genome-wide search was carried out for the first time to identify castor bean L/w genes and analyze the gene structures, biochemical properties, evolutionary relationships and expression characteristics based on the published data of castor bean genome and ESTs. According to the results, a total of 28 Rclhcs genes representing 13 gene families （ l_hca , l_hcb , Elip , Ohpl , Ohp2 , SEP1, SEP2 , SEP3 , SEP4 , SEP5 , PsbS , Rieske and FCII） and 25 subgene families were identified in castor bean genome; to be specific, 25 and 5 genes were found to have corresponding ESTs in NCBI and have al- ternative splicing isoforlns, respectively. These RcLhcs contain 0 to 9 introns and distribute on 26 of the 25 878 released scaffolds. All RcLhcs genes were found to be expressed in all examined tissues, i.e. leaf, flower, II/III stage endosperm, V/VI stage endosperm and seed, with the highest expression level in leaf tissue.

Molecular cloning and expression analysis of Crustin-like gene from Chinese shrimp Fenneropenaeus chinensis

A new member of antimicrobial protein genes of the Crustin family was cloned from haemocytes of the Chinese shrimp Fennero- penaeus chinensis by 3 ′and 5′ RACE. The full-length cDNA of Crustin-like gene contains a 390 bp open reading frame, encoding 130 amino acids. The deduced peptide contains a putative signal peptide of 17 amino acids and mature peptide of 113 amino acids. The molecular mass of the deduced mature peptide is 12. 3 ku. It is highly cationic with a theoretical isoelectric point of 8.5. The deduced amino acids sequence of this Crustin showed high homology with those of Penaeus （ Litopenaeus ） setferus. Northern blotting showed that the cloned Crustin gene was mainly expressed in haemocytes, gill, intestine, and RNA in situ hybridization indicated that the Crustin gene was constitutively expressed exclusively in haemocytes of these tissues. Capillary elec- trophoresis RT-PCR analysis showed that Crustin was up-regulated dramatically from 12 to 48 h after a brief decrease of mRNA during first 6 h in response to microbe infection. The level of Crustin mRNA began to restore at 72 h post-challenge. This indicated that Crustin gene might play an important role when shrimps are infected by bacterial pathogen.

Genome-Wide Survey and Expression Analysis of P_(1B)-ATPases in Rice, Maize and Sorghum

P1B-type ATPase ion pumps that transport heavy metal ions across cellular membranes are essential for plant growth and development. To date, a genomic comparison overview of the family in rice, maize and sorghum is not yet available. In this study, a total of 31 heavy metal P1B-type ATPase （HMA） genes were identified, including 9 in rice, 11 each from maize and sorghum. They were classified into two distinct subfamilies based on their sequence composition and phylogenetic relationship. Four pairs of HMA genes were expanded via gene duplication with tandemly duplicated. Comprehensive analyses were performed to investigate the expression profiles of HMA genes in various tissues by using quantitative real-time PCR. Some HMA members exhibited abundant and tissue-specific expression patterns. Moreover, most of the genes were found to be differentially expressed under the Cu/Cd treatment. This study will facilitate further studies on P1B-type ATPase family and provide valuable hints for the functional validation in rice, maize and sorghum.

Novel structural annotation and functional expression analysis of GTP_EFTU conserved genes in pepper based on the PacBio sequencing data

Genes containing GTP_EFTU domain mainly express elongation factors(EF),Small GTPases,and GTP-binding proteins,which are closely related to protein synthesis,extension and ATP synthesis.In this study,we identified 39 genes containing GTP_EFTU domains from peppers.The evolutionary trees constructed from capsicum,Arabidopsis,rice,and tomato are mainly divided into 7 subfamilies.Using PacBio(Pacific Biosciences)sequencing and assembly data,we extracted these 39 gene sequences,fromwhich 25 genes had alternative splicing.Particularly,the Capana08g000545 had 16 alternative splicing processes.Accordingly,we performed promoter sequence analysis,subcellular location prediction,the expression analysis of different tissues and periods,and also the GO(Gene ontology)analysis of co-expressed genes.Lastly we did the qRTPCR analysis in 5 stages of pepper fruit development.These analyses revealed important structural and functional information for the identified 39 genes that contain GTP_EFTU domains,providing important references for further follow-up experiments to verify the genes function on plants or their unique roles in peppers.

Cloning and Expression Analysis of a Mitogen-Activated Protein Kinase Gene OsMPK14 in Rice

Mitogen activated-protein kinases （MAPKs） are important components in signal transduction pathways responding to various biotic and abiotic stresses. An MAPK gene, OsMPK14 （GenBank Accession No. GQ265780） from rice （Oryza sativa L.）, was cloned by RT-PCR. The full-length cDNA of OsMPK14 consists of 1660 bp in size, containing an open reading frame of 1629 bp, which encodes a 542-amino-acid polypeptide and has a typical protein kinase domain and a phosphorylation activation motif TDY. Sequence alignment and analysis revealed that OsMPK14 was located on rice chromosome 5, and composed of nine exons and eight introns in the coding region. Semi-quantitative RT-PCR was performed to detect the expression patterns of OsMPK14 in rice shoots and roots under darkness, drought, high salinity, low temperature and abscisic acid treatments. The OsMPK14 mRNA was induced by abscisic acid, low temperature and high salinity, but weakly inhibited by drought. In addition, the expression of OsMPK14 was up-regulated in roots, but down-regulated in shoots by light. The results indicate that OsMPK14 could be implicated in diverse rice stimuli-responsive signaling cascades, and its expression might be regulated by multiple factors.

Genome-wide identification and expression analysis of Gossypium RING-H2 finger E3 ligase genes revealed their roles in fiber development,and phytohormone and abiotic stress responses

Background： RING H2 finger E3 ligase （RH2FE3） genes encode cysteine rich proteins that mediate E3 ubiquitin ligase activity and degrade target substrates. The roles of these genes in plant responses to phytohormones and abiotic stresses are well documented in various species, but their roles in cotton fiber development are poorly understood. To date, genome wide identification and expression analyses of Gossypium hirsutum RH2FE3 genes have not been reported. Methods： We performed computational identification, structural and phylogenetic analyses, chromosomal distribution analysis and estimated KJKs values of G hirsutum RH2FE3 genes. Orthologous and paralogous gene pairs were identified by all versus all BLASTP searches. We predicted cis regulatory elements and analyzed microarray data sets to generate heatmaps at different development stages. Tissue specific expression in cotton fiber, and hormonal and abiotic stress responses were determined by quantitative real time polymerase chain reaction （qRT PCR） analysis. Results： We investigated 140 G hirsutum, 80 G. orboreum, and evolutionary mechanisms and compared them with orthologs 89 G. roimondii putative RH2FB genes and their in Arobidopsis and rice. A domain based analysis of the G hirsutum RH2FE3 genes predicted conserved signature motifs and gene structures. Chromosomal localization showed the genes were distributed across all G hirsutum chromosomes, and 60 duplication events （4 tandem and 56 segmental duplications） and 98 orthologs were detected, cis elements were detected in the promoter regions of G hirsutum RH2FE3 genes. Microarray data and qRT PCR analyses showed that G hirsutum RH2FE3 genes were strongly correlated with cotton fiber development. Additionally, almost all the （brassinolide, gibberellic acid （GA）, indole 3-acetic acid drought, and salt）. dentified genes were up regulated in response to phytohormones （IAA）, and salicylic acid （SA）） and abiotic stresses （cold, heat, Conclusions： The genome wide identification, comprehensive analysis, and characterization of conserved domains and gene structures, as well as phylogenetic analysis, cis element prediction, and expression profile analysis of G hirsutum RH2FE3 genes and their roles in cotton fiber development and responses to plant hormones and abiotic stresses are reported here for the first time. Our findings will contribute to the genome wide analysis of putative RH2FE3 genes in other species and lay a foundation for future physiological and functional research on G hirsutum RH2FE3 genes.

Molecular cloning, characterization and expression analysis of a catalase gene in Paphia textile

Catalase is an important antioxidant protein that can protect organisms against various forms of oxidative damage by eliminating hydrogen peroxide. In this study, the catalase c DNA of Paphia textile（Pt CAT） was cloned using RTPCR and rapid amplification of c DNA ends（RACE）. Pt CAT is 1 921 bp long and consists of a 5′-UTR of 50 bp, a 3′-UTR of 349 bp, and an ORF of 1 542 bp that encodes 513 amino acids with a molecular weight of 58.4 k D and an estimated isoelectric point of 8.2. Sequence alignment indicated that Pt CAT contained a highly conserved catalytic signature motif（^（61）FNRERIPERVVHAKGAG^（77））, a proximal heme-ligand signature sequence（^（352）RLFSYSDP^（359））, and three catalytic amino acid residues（H^（72）, N^（145）, and Y^（356））. Pt CAT also contains two putative N-glycosylation sites（^（34）NKT^（36） and ^（437）NFT^（439）） and a peroxisome-targeting signal（^（511）AQL^（513））. Furthermore, Pt CAT shares 53%–88% identity and 29%–89% similarity with other catalase amino acid sequences. Pt CAT m RNA was present in all tested organs, including the heart, digestive gland, adductor muscle, gonad, gill, and mantle, but its expression was highest in the digestive gland. High-temperature-induced stress produced two expression patterns of Pt CAT m RNA： first, an initial up-regulation followed by a down-regulation in the heart, digestive gland, and gonad and, second, consistent down-regulation in all other organs. These results demonstrate that Pt CAT is a typical member of the catalase family and might be involved in the responses to harmful environmental factors.

Identification of potential immune-related prognostic biomarkers of lung cancer using gene co-expression network analysis

Objective The objective of this study was to identify new carcinogenetic hub genes and develop the integration of differentially expressed genes to predict the prognosis of lung cancer.Methods GSE139032 microarray data packages were downloaded from the Gene Expression Omnibus for planning,testing,and review of data.We identified KRT6C,LAMC2,LAMB3,KRT6A,and MYEOV from a key module for validation.Results We found that the five genes were related to a poor prognosis,and the expression levels of these genes were associated with tumor stage.Furthermore,Kaplan-Meier plotter showed that the five hub genes had better prognostic values.The mean levels of methylation in lung adenocarcinoma(LUAD)were significantly lower than those in healthy lung tissues for the hub genes.However,gene set enrichment analysis(GSEA)for single hub genes showed that all of them were immune-related.Conclusion Our findings demonstrated that KRT6C,LAMC2,LAMB3,KRT6A,and MYEOV are all candidate diagnostic and prognostic biomarkers for LUAD.They may have clinical implications in LUAD patients not only for the improvement of risk stratification but also for therapeutic decisions and prognosis prediction.

Genome-wide identification and expression analysis of the GhIQD gene family in upland cotton(Gossypium hirsutum L.)

Background:Calmodulin(CaM)is one of the most important Ca^(2+)signaling receptors because it regulates diverse physiological and biochemical reactions in plants.CaM functions by interacting with CaM-binding proteins(CaMBPs)to modulate Ca^(2+)signaling.IQ domain(IQD)proteins are plant-specific CaMBPs that bind to CaM by their specific CaM binding sites.Results:In this study,we identified 102 GhIQD genes in the Gossypium hirsutum L.genome.The GhIQD gene family was classified into four clusters(Ⅰ,Ⅱ,Ⅲ,andⅣ),and we then mapped the GhIQD genes to the G.hirsutum L.chromosomes.Moreover,we found that 100 of the 102 GhIQD genes resulted from segmental duplication events,indicating that segmental duplication is the main force driving GhIQD gene expansion.Gene expression pattern analysis showed that a total of 89 GhIQD genes expressed in the elongation stage and second cell wall biosynthesis stage of the fiber cells,suggesting that GhIQD genes may contribute to fiber cell development in cotton.In addition,we found that 20 selected GhIQD genes were highly expressed in various tissues.Exogenous application of MeJA significantly enhanced the expression levels of GhIQD genes.Conclusions:Our study shows that GhIQD genes are involved in fiber cell development in cotton and are also widely induced by MeJA.Thw results provide bases to systematically characterize the evolution and biological functions of GhIQD genes,as well as clues to breed better cotton varieties in the future.

Characterization and expression analysis of genes encoding Taxol biosynthetic enzymes in Taxus spp.

Taxol(Paclitaxel),an important anticancer drug,is derived at very low yields from Taxus(yew)species that grow very slowly.In the present study,thirteen genes that encode enzymes involved in Taxol biosynthesis in Taxus spp.were analyzed with bioinformatics methods,and their expression levels in different tissues and after cold and hormone treatments were also analyzed.The results indicated that many cis-elements related to abiotic stresses and hormones were found in the promoter sequences of the 8 genes involved in Taxol biosynthesis.Moreover,the 13 enzymes encoded by the target genes were located in different organelles and had many phosphorylation sites in the response proteins.The 13 genes were expressed highly either in roots or in stems,with lower transcripts in needles,and they were highly expressed after treatment with cold,gibberellin,methyl jasmonate or coronatine,consistent with predictions based on the bioinformatics analysis.These results suggest that the factors such as hormones and abiotic stresses stimulate taxane biosynthesis in yews,providing an important way to sustainably generate taxanes from yew trees or their cell cultures to improve Taxol yields.

Transcriptome analysis shows nicotinamide seed treatment alters expression of genes involved in defense and epigenetic processes in roots of seedlings of Picea abies

Nicotinamide treatment of plants and plant cell cultures has been shown to promote defense and decrease levels of DNA methylation.In the present study,we used RNA-seq technology to study overall changes in gene expression induced in roots of 3-month-old spruce(Picea abies)seedlings grown from nicotinamide-treated seeds to examine the molecular mechanisms underlying the defense promotion.Approximately 350 genes were identified as differentially expressed in roots after the seed treatment.Stress response genes,including transcription factors MYB77 and LHY and two chitinase enzymes,were generally upregulated,whereas genes thought to be involved in epigenetic regulation such as DDM1,known to promote DNA methylation,were present at high frequency among the downregulated genes.Across all samples,the expression of downregulated epigenetic-related genes was highly correlated with the nicotinamide treatment,indicating a common regulation.Our results support an earlier hypothesis regarding a potential role of nicotinamide as a defense-signal mediator.

Cloning and Expression Analysis of Sucrose Non-fermenting-1 Related Protein Kinase (SnRK) in Cucumis sativus L. Under Low Nitrogen Conditions

The sucrose non-fermenting-1 related protein kinase（SnRK）, whose expression is induced by kinds of hyperosmotic stresses, plays a key role in improving stress resistance of plants. In order to investigate the molecular mechanism of low nitrogen resistance in cucumber, the full-length cDNA of SnRK gene was cloned in this study. The result showed that SnRK gene was 1 548 bp in length, encoded 515 amino acids, and had more than 80% homology with other crops. The protein encoded by this gene was an unstable and hydrophilic protein with no transmembrane structure and no signal peptide. Under nitrogen-free conditions and low nitrogen conditions, the expression pattern analysis of SnRK gene showed that this gene was up-regulated and its expression increased and was significantly higher than the normal level as the nitrogen concentration decreased. In addition, the expression of SnRK gene was also inhibited in the high nitrogen level and was significantly lower than the normal level. The result of this study would help us understand the molecular mechanism of low nitrogen resistance in cucumber.

Alternative RNA splicing in stem cells and cancer stem cells:Importance of transcript-based expression analysis

Alternative ribonucleic acid(RNA)splicing can lead to the assembly of different protein isoforms with distinctive functions.The outcome of alternative splicing(AS)can result in a complete loss of function or the acquisition of new functions.There is a gap in knowledge of abnormal RNA splice variants promoting cancer stem cells(CSCs),and their prospective contribution in cancer progression.AS directly regulates the self-renewal features of stem cells(SCs)and stem-like cancer cells.Notably,octamer-binding transcription factor 4A spliced variant of octamerbinding transcription factor 4 contributes to maintaining stemness properties in both SCs and CSCs.The epithelial to mesenchymal transition pathway regulates the AS events in CSCs to maintain stemness.The alternative spliced variants of CSCs markers,including cluster of differentiation 44,aldehyde dehydrogenase,and doublecortin-like kinase,α6β1 integrin,have pivotal roles in increasing selfrenewal properties and maintaining the pluripotency of CSCs.Various splicing analysis tools are considered in this study.LeafCutter software can be considered as the best tool for differential splicing analysis and identification of the type of splicing events.Additionally,LeafCutter can be used for efficient mapping splicing quantitative trait loci.Altogether,the accumulating evidence re-enforces the fact that gene and protein expression need to be investigated in parallel with alternative splice variants.

Genome-Wide Identification of the F-box Gene Family and Expression Analysis under Drought and Salt Stress in Barley

The F-box protein-encoding gene family plays an essential role in plant stress resistance.In present study,126 non-redundant F-box genes were identified in barley(Hordeum vulgare L.,Hv).The corresponding proteins contained 165–887 amino acid residues and all were amphiphilic,except 5 proteins.Phylogenetic analysis of F-box protein sequences in barley and stress-related F-box protein sequences in wheat and Arabidopsis thaliana(At)was used to classify barley F-box genes are divided into 9 subfamilies(A–I).A structure-based sequence alignment demonstrated that F-box proteins were highly conserved with a total of 10 conserved motifs.In total,124 F-box genes were unevenly distributed on 7 chromosomes;another 2 genes have not been anchored yet.The gene structure analysis revealed high variability in the number of exons and introns in F-box genes.Comprehensive analysis of expression profiles and phylogenetic tree analysis,a total of 12 F-box genes that may be related to stress tolerance in barley were screened.Of the 12 detected F-box genes,8 and 10 were upregulated after drought and salt stress treatments,respectively,using quantitative real-time polymerase chain reaction(qRT-PCR).This study is the first systematic analysis conducted on the F-box gene family in barley,which is of great importance for clarifying this family’s bioinformatic characteristics and elucidating its function in barley stress resistance.These results will serve as a theoretical reference for subsequent research on molecular regulation mechanisms,genetic breeding,and improvement.

Characterization and Expression Analysis of Peroxiredoxin Genes in NNK-induced V79 Cells

4-（Methylnitrosamino）-1-（3-pyridyl）-1-butanone （NNK） is a potent and prevalent nitrosamine procarcinogen found in cigarette smoke. The aim of this work is to study alterations in peroxiredoxin （Prx） expression induced by NNK during carcinogenesis. Characterization of Prx genes from hamster was performed using bioinformatics.

Cloning and Expression Analysis of Mlo Gene from Pericallis hybrida B. Nord.

The full-length Mlo gene was obtained by reverse transcription polymerase chain reaction （RT-PCR） and RACE. The result of sequence analysis indicated that M/o gene from Pericallis hybrida B. Nord. contained about 1296bp open reading frame and encoded 431 amino acids. According to the comparison of the exogenous gene sequences by BLAST analysis and phylogenetic analysis, Mlo gene shared over 85% nucleotide homology and 98% amino acid homology. Finally, through semi-quantitative-PCR and fluorescence quantitative analysis, we found that Mlo gene showed the highest expression levels in leaves and the lowest in roots after inoculated with powdery mildew pathogen for different days.