Gene Ontology Functional Classification Research Articles

PGD: Shared gene linking polycystic ovary syndrome and endometrial cancer, influencing proliferation and migration through glycometabolism.

The relationship among polycystic ovary syndrome (PCOS), endometrial cancer (EC), and glycometabolism remains unclear. We explored shared genes between PCOS and EC, using bioinformatics to unveil their pathogenic connection and influence on EC prognosis. Gene Expression Omnibus datasets GSE226146 (PCOS) and GSE196033 (EC) were used. A protein-protein interaction (PPI) network was constructed to identify the central genes. Candidate markers were screened using dataset GSE54250. Differences in marker expression were confirmed in mouse PCOS and human EC tissues using RT-PCR and immunohistochemistry. The effect of PGD on EC proliferation and migration was explored using Ki-67 and Transwell assays. PGD's impact on the glycometabolic pathway within carbon metabolism was assessed by quantifying glucose content and lactic acid production. R software identified 31 common genes in GSE226146 and GSE196033. Gene Ontology functional classification revealed enrichment in the "purine nucleoside triphosphate metabolism process," with key Kyoto Encyclopedia of Genes and Genomes pathways related to "carbon metabolism." The PPI network identified 15 hub genes. HK2, NDUFS8, PHGDH, PGD, and SMAD3 were confirmed as candidate markers. The RT-PCR analysis validated distinct HK2 and PGD expression patterns in mouse PCOS ovarian tissue and human EC tissue, as well as in normal and EC cells. Transfection experiments with Ishikawa cells further confirmed PGD's influence on cell proliferation and migration. Suppression of PGD expression impeded glycometabolism within the carbon metabolism of EC cells, suggesting PGD as a significant PCOS risk factor impacting EC proliferation and migration through modulation of single carbon metabolism. These findings highlight PGD's pivotal role in EC onset and prognosis.

Identification of hub genes and immune infiltration in ulcerative colitis using bioinformatics

Ulcerative colitis (UC) is a chronic inflammatory disease of the intestine, whose pathogenesis is not fully understood. Given that immune infiltration plays a key role in UC progression, our study aimed to assess the level of immune cells in UC intestinal mucosal tissues and identify potential immune-related genes. The GSE65114 UC dataset was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between healthy and UC tissues were identified using the “limma” package in R, while their Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were determined with the clusterProfiler package. Protein–protein interaction network analysis and visualization were performed with STRING and Cytoscape. Immune cell infiltration was calculated with CIBERSORT. The relationship between hub genes and immune-infiltrated cells in UC was determined by Pearson correlation. A total of 206 DEGs were identified, of which 174 were upregulated and 32 downregulated. GO and KEGG functional classification indicated DEG enrichment in immune response pathways, including Toll-like receptor signaling, IL-17 signaling, and immune system process and chemokine signaling. 13 hub genes were identified. Infiltration matrix analysis of immune cells showed abundant plasma cells, memory B cells, resting CD4 memory T cells, γδ T cells, M0 and M1 macrophages, and neutrophils in UC intestinal tissues. Correlation analysis revealed 13 hub genes associated with immune-infiltrated cells in UC. 13 hub genes associated with immune-infiltrated cells in UC were identified; they included CXCL13, CXCL10, CXCL9, CXCL8, CCL19, CTLA4, CCR1, CD69, CD163, IL7R, PECAM1, TLR8 and TLR2. These genes could potentially serve as markers for the diagnosis and treatment of UC.

Transcriptional analysis reveals formation of axillary solitary hook in vine plant Uncaria rhynchophylla

The hook of Uncaria rhynchophylla is an important external phenotype representing the quality of Gou-Teng. The hook formation pattern of U. rhynchophylla displays a unique feature that solitary hook and opposite hooks alternately occur in the leaf axils forming a liner series. The solitary hook only grows on the lower side of the leaf axil, resulting from the decline in the total yield of hooks. However, the reasons for solitary hook formation in U. rhynchophylla are not clear. Therefore, a survey on the molecular mechanisms of solitary hook formation needs to be performed urgently. We obtained 250.67 million high-quality clean reads from the lower side of axillary stem segments (HS) and the upper side of axillary stem segments (HSCK) libraries. We detected 42 DEGs (differentially expressed genes) between HS and HSCK libraries; 11 were down-regulated, and 31 were up-regulated. Gene ontology functional classification of the DEGs indicated that they included two genes that encoded LAZY1 related to the regulation of polar auxin transport, two genes that encoded AUX/IAA and AP2/ERF related to plant hormone signaling transduction, and two genes that encoded F3H related to the flavone synthase activity. These six genes were associated with lateral polar auxin transport (PAT) regulation, and their expression level had a high correlation with the pattern of IAA distribution upon gravity stimulation. And the inhibition of solitary hook development was observed after treatment with N-1-naphthylphthalamic acid (NPA) or quercetin, which was known as the polar auxin inhibitor. The above findings might provide new information about the molecular mechanisms of regulation of lateral PAT in the axillary bud development process under gravity.

Identification of hub ferroptosis-related genes and immune infiltration in lupus nephritis using bioinformatics

Lupus nephritis (LN) is one of the most severe and more common organ manifestations of the autoimmune disease, systemic lupus erythematosus. Ferroptosis, a novel type of programmed cell death, so far its role in LN remains uncertain. In the present study, we explored the role of ferroptosis in LN and its relationship with the immune response. The GSE112943 LN dataset was downloaded from the Gene Expression Omnibus database. Ferroptosis-Related Genes (FRGs) that drive, suppress or mark ferroptosis were retrieved from the public FerrDb database. The gene expression matrix of the GSE112943 dataset was analyzed with the “limma” package in R to obtain differentially expressed genes (DEGs) between LN and healthy samples. Subsequently, the crossover genes between DEGs and FRGs were identified as differentially expressed ferroptosis-related genes (DE-FRGs). Protein–protein interaction (PPI) network analysis, visualization, and identification of hub lupus nephritis ferroptosis-related genes (LN-FRGs) were performed with STRING and Cytoscape, while their Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were determined with the clusterProfiler package. Immune cell infiltration was calculated with CIBERSORT. The relationship between hub LN-FRGs and immune-infiltrated cells in LN was determined by Pearson correlation. A total of 96 DE-FRGs and 8 hub LN-FRGs (KRAS, PIK3CA, EGFR, MAPK14, SRC, MAPK3, VEGFA, and ATM) were identified. GO and KEGG functional classification indicated these genes enrichment in apoptotic process, programmed cell death, autophagy-animal, FoxO signaling pathway, relaxin signaling pathway, and VEGF signaling pathway. Infiltration matrix analysis of immune cells showed abundant Monocytes and M0/M1/M2 macrophages in LN kidney tissues. Correlation analysis revealed 8 hub LN-FRGs associated with immune-infiltrated cells in LN. In summary, overproduction of ROS and abnormal infiltration of immune cells would be implicated in the LN caused by ferroptosis. 8 hub lupus nephritis ferroptosis-related genes (LN-FRGs) which might be good biomarkers of ferroptosis in LN were identified in this study. These findings point to the immune response playing an important role in LN caused by ferroptosis via mutual regulation between hub LN-FRGs and immune-infiltrated cells.

Molecular mechanism of locus coeruleus involved in acupuncture anti-myocardial ischemia based on transcriptome sequencing

To explore the molecular mechanism of locus coeruleus（LC） involved in electroacupuncture (EA) anti myocardial ischemia. Twenty-four SD rats were randomly divided into sham-operation, model, EA and EA +lesion groups, with 6 rats in each group. The acute myocardial ischemia (AMI) model was established by ligation of the left anterior descending branch of coronary artery. EA (2 Hz/15 Hz, 1 mA) was applied to bilateral "Shenmen" (HT7) -"Tongli" (HT5) and the middle-point between HT7 and HT5 for 30 min, once daily for 3 days. For rats of the EA +lesion group, the virus (300 nL) was injected into bilateral LC before EA treatment. Serum aspartate aminotransferase (AST) was detected by ELISA. The gene expression profiles of rat heart were detected by transcriptome sequencing, the differentially expressed genes were screened, and Gene Ontology (GO) functional classification and Kyoto Encyclopedia of genes and genomes (KEGG) metabolic pathway enrichment analysis were performed. Compared with the sham-operation group, serum AST content was significantly increased in the model group (P<0.01). Following the intervention, serum AST was significantly reduced in the EA group (P<0.01), while the serum AST in the EA + lesion group was significantly higher compared with the EA group (P<0.05). Differential expression analysis showed that 1 138 differentially expressed genes were screened out between the model group and the sham-operation group, 1 330 differentially expressed genes between model and EA group, and 804 differentially expressed genes between EA and EA + lesion group. Among them, 218 differential genes were involved in the regulation of EA anti-myocardial ischemia in LC. GO functional classification analysis showed that these differentially expressed genes mainly involved in cell processes, metabolic processes and biological regulation in biological processes. KEGG pathway analysis showed that these differentially expressed genes were enriched in sulfur relay system, thiamine metabolism, glutathione metabolism, C5 branch dicarboxylic acid metabolism, cell adhesion molecules and Th1 and Th2 cell differentiation. EA intervention has a positive effect in anti-myocardial ischemia, which may be related to the sulfur relay system, thiamine metabolism, glutathione metabolism, C5 branch dicarboxylic acid metabolism, cell adhesion molecules and Th1 and Th2 cell differentiation involved in LC.

Isobaric Tags for Relative and Absolute Quantitation-Based Quantitative Serum Proteomics Analysis in Ischemic Stroke Patients With Hemorrhagic Transformation

Hemorrhagic transformation (HT), which occurs with or without reperfusion treatments (thrombolysis and/or thrombectomy), deteriorates the outcomes of ischemic stroke patients. It is essential to find clinically reliable biomarkers that can predict HT. In this study, we screened for potential serum biomarkers from an existing blood bank and database with 243 suspected acute ischemic stroke (AIS) patients. A total of 37 patients were enrolled, who were diagnosed as AIS without receiving reperfusion treatment. They were divided into two groups based on whether they were accompanied with HT or not (five HT and 32 non-HT). Serum samples were labeled by isobaric tags for relative and absolute quantitation (iTRAQ) and analyzed by liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) and compared under NCBInr database. A total of 647 proteins in sera samples were captured, and the levels of 17 proteins (12 upregulated and five downregulated) were significantly different. These differentially expressed proteins were further categorized with Gene Ontology functional classification annotation and Kyoto Encyclopedia of Genes and Genomes metabolic pathway analysis into biological processes. Further protein–protein interaction analysis using String database discovered that, among the differentially expressed proteins, 10 pairs of proteins were found to have crosstalk connections, which may have direct (physical) and indirect (functional) interactions for the development of HT. Our findings suggest that these differentially expressed proteins could serve as potential biomarkers for predicting HT after ischemic stroke.

Survival, Virulent Characteristics, and Transcriptomic Analyses of the Pathogenic Vibrio anguillarum Under Starvation Stress.

Many bacteria have developed strategies for metamorphosis into sophisticated survival forms to survive extended periods of environmental stress. As a global causative agent of vibriosis in marine fish farming, Vibrio anguillarum (V. anguillarum) can efficiently grow and proliferate under environmental stress, but the specific mechanism is not clear. In the present study, survival, virulent characteristics, and transcriptomic analysis of the V. anguillarum BH1 were performed under starvation stress. The results demonstrated that V. anguillarum was still culturable and showed rippled surface after 6 months of starvation. Starved cells maintained their infectivity in half-smooth tongue sole (Cynoglossus semilaevi). Detection of virulence factors and virulence-associated genes in starved cells showed that the starved strain still produced β-hemolysis on rabbit blood agar, caseinase, dnase, and gelatinase, and possessed empA, vah1, vah2, vah3, vah4, vah5, rtxA, flaA, flaD, flaE, virC, tonB, mreB, toxR, rpoS, and ftsZ virulence-related genes. In addition, we first reported the RNA-seq study for V. anguillarum with and without starvation treatment for a period of 6 months and emphasized the regulation of gene expression at the whole transcriptional level. It indicated that V. anguillarum expressed 3,089 and 3,072 genes in the control group and starvation stress group, respectively. The differently expressed genes (DEGs) of the starved strain were thereby identified, including 251 up-regulated genes and 272 down-regulated genes in comparison with the non-starved strain. Gene Ontology (GO) analysis and Kyto Encyclopedia Genes and Genomes (KEGG) enrichment analysis of DEGs were also analyzed. GO functional classification revealed that among the significantly regulated genes with known function categories, more genes affiliated with signal transducer activity, molecular transducer activity, and cell communication were significantly up-regulated, and more genes affiliated with cellular macromolecule, cellular component, and structural molecule activity were significantly down-regulated. In addition, the DEGs involved in the pathway of two-component system was significantly up-regulated, and the pathways of ribosome and flagellar assembly were significantly down-regulated. This study provides valuable insight into the survival strategies of V. anguillarum and suggests that a portion of the bacterial populations may remain pathogenic while persisting under starvation stress by up-regulating or down-regulating a series of genes.

The early transcriptome response of cassava (Manihot esculenta Crantz) to mealybug (Phenacoccus manihoti) feeding.

The mealybug, Phenacoccus manihoti, is a leading pest of cassava (Manihot esculenta Crantz), damaging this crop globally. Although the biological control of this mealybug using natural predators has been established, resistance breeding remains an important means of control. Understanding plant responses to insect herbivory, by determining and identifying differentially expressed genes (DEGs), is a vital step towards the understanding of molecular mechanisms of defence responses in plants and the development of resistant cultivars by gene editing. Morphological and molecular analysis confirmed the mealybug identity as Phenacoccus manihoti (Matile-Ferrero). The transcriptome response of the green mite resistant cassava genotype AR23.1 was compared to P40/1 with no known resistance at 24 and 72 hours of mealybug infestation compared to non-infested mock. A total of 301 and 206 genes were differentially expressed at 24 and 72 of mealybug infestation for AR23.1 and P40/1 genotypes respectively, using a log2 fold change and P-value ≤ 0.05. Gene ontology functional classification revealed an enrichment of genes in the secondary metabolic process category in AR23.1 in comparison with P40/1, while genes in the regulation of molecular function, cellular component biogenesis and electron carrier categories were more significantly enriched in P40/1 than in AR23.1. Biological pathway analysis, based on KEGG, revealed a significant enrichment of plant-pathogen interaction and plant hormonal signal transduction pathways for a cohort of up-regulated and down-regulated DEGs in both genotypes. Defence-related genes such as 2-oxogluterate, gibberellin oxidase and terpene synthase proteins were only induced in genotype AR23.1 and not in P40/1, and subsequently validated by RT-qPCR. The study revealed a difference in response to mealybug infestation in the two genotypes studied, with AR23.1 showing a higher number of differentially expressed transcripts post mealybug infestation at 24 and 72 hours. Candidate defence-related genes that were overexpressed in the AR23.1 genotype post mealybug infestation will be useful in future functional studies towards the control of mealybugs.

Transcriptome sequencing of active buds from Populus deltoides CL. and Populus × zhaiguanheibaiyang reveals phytohormones involved in branching

Branching in woody plants affects their ecological benefits and impacts wood formation. To obtain genome-wide insights into the transcriptome changes and regulatory mechanisms associated with branching, we performed high-throughput RNA sequencing to characterize cDNA libraries generated from active buds of Populus deltoides CL. ‘zhonglin2025’ (BC) and Populus × zhaiguanheibaiyang (NC). NC has more branches than BC and rapid growth. We obtained a total of 198.2 million high-quality clean reads from the NC and BC libraries. We detected 3543 differentially expressed genes (DEGs) between the NC and BC libraries; 1418 were down-regulated and 2125 were up-regulated. Gene ontology functional classification of the DEGs indicated that they included 89 genes that encoded proteins related to hormone biosynthesis, 364 genes related to hormone signaling transduction, and 104 related to the auxin efflux transmembrane transporter. We validated the expression profiles of 16° by real-time quantitative PCR and found that their expression patterns were similar to those obtained from the high-throughput RNA sequencing data. We also measured the hormone content in young buds of BC and NC by high-pressure liquid chromatography. In this study, we identified global hormone regulatory patterns and differences in gene expression between NC and BC, and constructed a hormone regulatory network to explain branching in Populus buds. In addition, candidate genes that may be useful for molecular breeding of particular plant types were identified. Our results will provide a starting point for future investigations into the molecular mechanisms of branching in Populus.

Detection of Genes Associated with Follicle Development Through Transcriptome Analysis of Bovine Ovarian Follicles GCs

Background: The ovarian follicle is an essential component of the reproductive process. It plays an important role in controlling the estrous cycle, determining estrous behaviour, ensuring oocyte competency and subsequent embryo survival rate, and determining both post ovulation corpus luteum function and progesterone synthesis. Traditionally, gene expression studies in the field of follicular development focus on the study of expression of candidate genes of interest. With the development of next-generation sequencing technologies, transcriptome profiling has become a powerful approach for identification of genes globally expressed in various tissues including ovarian follicles. Objective: To investigate potential differentially expressed genes associated with follicular development and gene expression profiles of bovine ovarian follicles at onset of deviation stages of a follicular wave. Method: Ovaries were removed, greatest (ODF1) and second greatest (ODF2) diameter follicles were isolated, and GCs were isolated from the two types of follicles (n=4). RNA samples of follicles were pooled within ODF1 and ODF2. The two cDNA libraries were sequenced, and mapped to the bovine RefSeq database. Gene ontology (GO) functional enrichment analysis was subjected and differentially expressed genes were analysed using visualization plug BiNGO of Cytoscape 3.0.0 software. KEGG pathways were performed and all the pathways and protein equivalent convert information of UniProt database of cattle. Twelve genes were selected randomly to verify their expression by QRT-PCR technology. Results: A total of 43,708,132 and 43,826,914 clean reads were obtained for ODF1 and ODF2, respectively. After being mapped to the bovine RefSeq database, a total of 15,519 genes were expressed (cut-off RPKM>0.5), of which 761 were highly expressed (cut-off RPKM>100) in both types of follicles. GO functional classification of these highly expressed genes revealed that many of them are involved in metabolic processes, multicellular organismal processes, and binding. Furthermore, 831 genes were identified as differentially expressed in ODF1 versus ODF2, in which 384 genes were up-regulated, and 447 genes were down-regulated. KEGG pathway analysis revealed eight differentially expressed genes in ODF1 vs. ODF2, QRT-PCR results showed that 9 of the 12 selected genes mRNA amounts were significantly greater in ODF1 than that ODF2 (P<0.05 or 0.01), which was consistent with the deep sequencing data. Conclusion: This study represents the first transcriptome investigation of bovine follicles at onset of deviation of the largest vs. second largest follicles. It provides a foundation for future investigation of the regulatory mechanisms involved in follicular development in cattle. Keywords: Follicular development, transcriptome, gene expression, GO analysis, onset of deviation.

Transcriptome profiling of Staphylococci-infected cow mammary gland parenchyma

BackgroundGenome-wide gene expression profiling allows for identification of genes involved in the defense response of the host against pathogens. As presented here, transcriptomic analysis and bioinformatics tools were applied in order to identify genes expressed in the mammary gland parenchyma of cows naturally infected with coagulase-positive and coagulase-negative Staphylococci.ResultsIn cows infected with coagulase-positive Staphylococci, being in 1st or 2nd lactation, 1700 differentially expressed genes (DEGs) were identified. However, examination of the 3rd or 4th lactations revealed 2200 DEGs. Gene ontology functional classification showed the molecular functions of the DEGs overrepresented the activity of cytokines, chemokines, and their receptors. In cows infected with coagulase-negative Staphylococci, in the 1st or 2nd lactations 418 DEGs, while in the 3rd or 4th lactations, 1200 DEGs were identified that involved in molecular functions such as protein, calcium ion and lipid binding, chemokine activity, and protein homodimerization. Gene network analysis showed DEGs associated with inflammation, cell migration, and immune response to infection, development of cells and tissues, and humoral responses to infections caused by both types of Staphylococci.ConclusionA coagulase-positive Staphylococci infection caused a markedly stronger host response than that of coagulase-negative, resulting in vastly increased DEGs. A significant increase in the expression of the FOS, TNF, and genes encoding the major histocompatibility complex proteins (MHC) was observed. It suggests these genes play a key role in the synchronization of the immune response of the cow’s parenchyma against mastitis-causing bacteria. Moreover, the following genes that belong to several physiological pathways (KEGG pathways) were selected for further studies as candidate genes of mammary gland immune response for use in Marker Assisted Selection (MAS): chemokine signaling pathway (CCL2, CXCL5, HCK, CCR1), cell adhesion molecules (CAMs) pathway (BOLA-DQA2, BOLA-DQA1, F11R, ITGAL, CD86), antigen processing and presentation pathway (CD8A, PDIA3, LGMN, IFI30, HSPA1A), and NOD-like receptor signaling pathway (TNF, IL8, IL18, NFKBIA).

Alteration of gene expression by zinc oxide nanoparticles or zinc sulfate in vivo and comparison with in vitro data: A harmonious case

Granulosa cells (GCs) are those somatic cells closest to the female germ cell. GCs play a vital role in oocyte growth and development, and the oocyte is necessary for multiplication of a species. Zinc oxide (ZnO) nanoparticles (NPs) readily cross biologic barriers to be absorbed into biologic systems that make them promising candidates as food additives. The objective of the present investigation was to explore the impact of intact NPs on gene expression and the functional classification of altered genes in hen GCs in vivo, to compare the data from in vivo and in vitro studies, and finally to point out the adverse effects of ZnO NPs on the reproductive system. After a 24-week treatment, hen GCs were isolated and gene expression was quantified. Intact NPs were found in the ovary and other organs. Zn levels were similar in ZnO-NP-100 mg/kg– and ZnSO4-100 mg/kg–treated hen ovaries. ZnO-NP-100 mg/kg and ZnSO4-100 mg/kg regulated the expression of the same sets of genes, and they also altered the expression of different sets of genes individually. The number of genes altered by the ZnO-NP-100 mg/kg and ZnSO4-100 mg/kg treatments was different. Gene Ontology (GO) functional analysis reported that different results for the two treatments and, in Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, 12 pathways (out of the top 20 pathways) in each treatment were different. These results suggested that intact NPs and Zn2+ had different effects on gene expression in GCs in vivo. In our recent publication, we noted that intact NPs and Zn2+ differentially altered gene expression in GCs in vitro. However, GO functional classification and KEGG pathway enrichment analyses revealed close similarities for the changed genes in vivo and in vitro after ZnO NP treatment. Furthermore, close similarities were observed for the changed genes after ZnSO4 treatments in vivo and in vitro by GO functional classification and KEGG pathway enrichment analyses. Therefore, the effects of ZnO NPs on gene expression in vitro might represent their effects on gene expression in vivo. The results from this study and our earlier studies support previous findings indicating ZnO NPs promote adverse effects on organisms. Therefore, precautions should be taken when ZnO NPs are used as diet additives for hens because they might cause reproductive issues.

Identification of Differentially Expressed Genes Associated with Apple Fruit Ripening and Softening by Suppression Subtractive Hybridization.

Apple is one of the most economically important horticultural fruit crops worldwide. It is critical to gain insights into fruit ripening and softening to improve apple fruit quality and extend shelf life. In this study, forward and reverse suppression subtractive hybridization libraries were generated from ‘Taishanzaoxia’ apple fruits sampled around the ethylene climacteric to isolate ripening- and softening-related genes. A set of 648 unigenes were derived from sequence alignment and cluster assembly of 918 expressed sequence tags. According to gene ontology functional classification, 390 out of 443 unigenes (88%) were assigned to the biological process category, 356 unigenes (80%) were classified in the molecular function category, and 381 unigenes (86%) were allocated to the cellular component category. A total of 26 unigenes differentially expressed during fruit development period were analyzed by quantitative RT-PCR. These genes were involved in cell wall modification, anthocyanin biosynthesis, aroma production, stress response, metabolism, transcription, or were non-annotated. Some genes associated with cell wall modification, anthocyanin biosynthesis and aroma production were up-regulated and significantly correlated with ethylene production, suggesting that fruit texture, coloration and aroma may be regulated by ethylene in ‘Taishanzaoxia’. Some of the identified unigenes associated with fruit ripening and softening have not been characterized in public databases. The results contribute to an improved characterization of changes in gene expression during apple fruit ripening and softening.

In silico characterization of hypothetical proteins from Paracoccidioides lutzii.

Nearly 60% of Paracoccidioides lutzii genes encode products annotated as hypothetical or predicted proteins (HPs). In this study, we describe the global detection and functional inference of HPs, using computational methods based on sequence similarity, identification of targeting signals, presence of known protein domains, and use of the Gene Ontology functional classification scheme. Our analysis enabled a high-throughput characterization of predicted cellular localization and presence of protein domains, clustering HPs into different functional categories including metabolism, localization, cell cycle, response to stimulus, and signaling. To investigate P. lutzii HP expression profiles, we used data obtained from the expressed sequence tag database (dbEST). These analyses revealed 2364 HPs expressed in different situations, namely in mycelial and yeast forms, during the transition from mycelium to yeast, and under conditions mimicking infection. Based on this transcriptomic data, we performed a functional enrichment analysis according to the domains present in the HPs expressed in each condition. The most overrepresented functional domains were those involved in the regulation of gene expression, suggesting important and as yet undescribed roles for these HPs in the adaptation of P. lutzii to environmental conditions. In addition, the expression profiles of six randomly selected HPs were analyzed by quantitative real-time polymerase chain reaction in order to verify their expression in the complementary DNA libraries analyzed in this investigation. The approach used in this study should improve functional characterization of P. lutzii HPs.

Identification of heat-related ESTs in moth bean through suppression subtraction hybridization.

Moth bean (Vigna aconitifolia (Jacq.) Marechal), an important grain-legume crop grown in hot desert regions of Thar, under scorching sun rays, was investigated for heat tolerance at molecular level. In the present study, we constructed a forward suppression subtractive hybridization (SSH) cDNA library of heat tolerant genotype RMO-40 to identify genes expressing under delayed response to elevated temperature. Heat induction was carried out by exposing 14-day-old seedlings to elevated temperature of 42 °C for 30 min. A total of 125 unigenes (33 contigs and 92 singletons) were derived by cluster assembly and sequence alignment of 200 ESTs; out of 125 unigenes, 21 (16 %) were found to be novel to moth bean. Gene ontology functional classification terms were retrieved for 98 (78.4 %) unigenes of which 73 (58.4 %) ESTs were functionally annotated (GO consensus) where 19 unigenes were annotated with 11 enzyme commission (EC) codes and were mapped to 25 different KEGG pathways. We have identified a majority of heat-shock proteins (constituting 35 % of the present library) aiding heat stress tolerance to moth bean. An expression level of 22 ESTs generated from the above SSH cDNA library was studied through semiquantitative RT-PCR assay simultaneously under 5 and 30 min of heat stress at 42 °C.

Elucidating cold acclimation pathway in blueberry by transcriptome profiling

A fundamental goal of cold acclimation research is to understand the mechanisms responsible for the increase in freezing tolerance in response to environmental cues. Changes in gene expression underlie some of the biochemical and physiological changes that occur during cold acclimation. Detailed and comprehensive transcriptome annotation can be considered a prerequisite for effective analysis and a fast and cost-effective way to rapidly obtain information in the context of a given physiological condition. By computational predictions and manual curation, we have annotated 454 sequence assemblies from two blueberry cDNA libraries that represent flower buds in the first and second stages of cold acclimation. Gene ontology functional classification terms were retrieved for 4343 (80.0%) sequences. GO annotation files compatible with a commonly used annotation tool have been generated and are publicly available. By mining the dataset further, it was possible to associate presence of certain transcripts related to carbohydrate metabolism and lipid metabolism with different stages of cold acclimation. This was concomitant with differential presence of Zn finger functional domains and C3H-family transcription factors. The expression of a few selected genes was validated by quantitative real-time PCR assay. Results demonstrate that our transcriptome database is a rich resource for mining cold acclimation-responsive genes.

Generation and analysis of expressed sequence tags for microsatellite marker development in Calamus simplicifolius C. F. Wei

Rattans serve as an important source of raw non-wood materials for furniture and handicraft industries worldwide. However, their genomic sequence information in public databases is very limited. In this study, a set of 2,528 good-quality expressed sequence tags (ESTs) were generated from a full-length cDNA library constructed previously with root, stem and male inflorescence tissues of Calamus simplicifolius C. F. Wei, a rattan species native to Hainan Island, China. The ESTs were assembled into 1,588 unigenes, including 1,221 singletons and 367 contigs. BlastX searches against the GenBank non-redundant protein database revealed that 1,248 (78.6 %) unigenes had at least one significant match (E ≤ 10−5). The gene ontology functional classification assigned 991, 669 and 977 of the unigenes to the cellular component, molecular function and biological process categories, respectively. A total of 71 simple sequence repeat (SSR) loci were developed among these ESTs, including 65 polymorphic across 19 rattan species representing three genera. High levels of cross-species/genus transferability were observed for the EST-SSRs. For the polymorphic EST-SSR markers, the number of alleles per locus and polymorphic information content ranged from 2 to 25 (mean 11.1) and from 0.135 to 0.949 (mean 0.695), respectively. The EST sequences and the EST-SSR primers have been deposited in GenBank databases of EST (IDs JK838364–40891) and Probe (IDs Pr16718978–9048, to be assigned).

Construction and analysis of an SSH cDNA library of early heat-induced genes ofVigna aconitifoliavariety RMO-40

Moth bean ( Vigna aconitifolia (Jacq.) Marechal) is an important grain legume crop grown in rain fed areas of hot desert regions of Thar, India, under scorching sun rays with very little supplementation of water. An SSH cDNA library was generated from leaf tissues of V. aconitifolia var. RMO-40 exposed to an elevated temperature of 42°C for 5min to identify early-induced genes. A total of 488 unigenes (114 contigs and 374 singletons) were derived by cluster assembly and sequence alignment of 738 ESTs; out of 206 ESTs (28%) of unknown proteins, 160 ESTs (14%) were found to be novel to moth bean. Only 578 ESTs (78%) showed significant BLASTX similarity (<1 × 10(-6)) in the NCBI non-redundant database. Gene ontology functional classification terms were retrieved for 479 (65%) sequences, and 339 sequences were annotated with 165 EC codes and mapped to 68 different KEGG pathways. Four hundred and fifty-two ESTs were further annotated with InterProScan (IPS), and no IPS was assigned to 153 ESTs. In addition, the expression level of 27 ESTs in response to heat stress was evaluated through semiquantitative RT-PCR assay. Approximately 20 different signaling genes and 16 different transcription factors have been shown to be associated with heat stress in moth bean for the first time.

Determination of the mechanism of action of repetitive halothane exposure on rat brain tissues using a combined method of microarray gene expression profiling and bioinformatics analysis

The present study aimed to investigate the gene expression profiles of rats brain tissues treated with halothane compared with untreated controls to improve current understanding of the mechanism of action of the inhaled anesthetic. The GSE357 gene expression profile was dowloaded from the Gene Expression Omnibus database, and included six gene chips of samples repeatedly exposed to halothane and 12 gene chips of untreated controls. The differentially expressed genes (DEGs) between these two groups were identified using the Limma package in R language. Subsequently, the Database for Annotation, Visualization and Integrated Discovery was used to annotate the function of these DEGs. In addition, the most significantly upregulated gene and downregulated gene were annotated, to reveal the functional interactions with other associated genes, in FuncBase database. A total of 44DEGs were obtained between The control and halothane exposure samples. Following Gene Ontology functional classification, these DEGs were found to be involved predominantly in the circulatory system, regulation of cell proliferation and response to endogenous stimulus and corticosteroid stimulus processes. KRT31 and HMGCS2, which were identified as the most significantly downregulated and upregulated DEGs, respectively, were associated with the lipid metabolic process and Tcell activation, respectively. These results provided a basis for the development of improved inhalational anesthetics with minimal side effects and are essential for optimization of inhaled anesthetic techniques for advanced surgical procedures.

Comparative analysis of expressed sequence tags (ESTs) between drought-tolerant and -susceptible genotypes of chickpea under terminal drought stress

BackgroundChickpea (Cicer arietinum L.) is an important grain-legume crop that is mainly grown in rainfed areas, where terminal drought is a major constraint to its productivity. We generated expressed sequence tags (ESTs) by suppression subtraction hybridization (SSH) to identify differentially expressed genes in drought-tolerant and -susceptible genotypes in chickpea.ResultsEST libraries were generated by SSH from root and shoot tissues of IC4958 (drought tolerant) and ICC 1882 (drought resistant) exposed to terminal drought conditions by the dry down method. SSH libraries were also constructed by using 2 sets of bulks prepared from the RNA of root tissues from selected recombinant inbred lines (RILs) (10 each) for the extreme high and low root biomass phenotype. A total of 3062 unigenes (638 contigs and 2424 singletons), 51.4% of which were novel in chickpea, were derived by cluster assembly and sequence alignment of 5949 ESTs. Only 2185 (71%) unigenes showed significant BLASTX similarity (<1E-06) in the NCBI non-redundant (nr) database. Gene ontology functional classification terms (BLASTX results and GO term), were retrieved for 2006 (92.0%) sequences, and 656 sequences were further annotated with 812 Enzyme Commission (EC) codes and were mapped to 108 different KEGG pathways. In addition, expression status of 830 unigenes in response to terminal drought stress was evaluated using macro-array (dot blots). The expression of few selected genes was validated by northern blotting and quantitative real-time PCR assay.ConclusionOur study compares not only genes that are up- and down-regulated in a drought-tolerant genotype under terminal drought stress and a drought susceptible genotype but also between the bulks of the selected RILs exhibiting extreme phenotypes. More than 50% of the genes identified have been shown to be associated with drought stress in chickpea for the first time. This study not only serves as resource for marker discovery, but can provide a better insight into the selection of candidate genes (both up- and downregulated) associated with drought tolerance. These results can be used to identify suitable targets for manipulating the drought-tolerance trait in chickpea.