Differential gene expression profiling in mucus glands of honey bee (Apis mellifera) drones during sexual maturation

The mating sign that each drone leaves when mating with a queen essentially consists of mucus gland proteins. We employed a Representational Difference Analysis (RDA) methodology to identify genes that are differentially expressed in mucus glands during sexual maturation of drones. The RDA library for mucus glands of newly emerged drones was more complex than that of 8 day-old drones, with matches to 20 predicted genes. Another 26 reads matched to the Apis genome but not to any predicted gene. Since these ESTs were located within ORFs they may represent novel honey bee genes, possibly fast evolving mucus gland proteins. In the RDA library for mucus glands of 8 day-old drones, most reads corresponded to a capsid protein of deformed wing virus, indicating high viral loads in these glands. The expression of two genes encoding venom allergens, acid phosphatase-1 and hyaluronidase, in drone mucus glands argues for their homology with the female venom glands, both associated with the reproductive system.

A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Genes up- and down-regulated by dermcidin in breast cancer: a microarray analysis

Dermcidin (DCD) is a human gene mapped to chromosome 12q13 region, which is co-amplified with multiple oncogenes with a well-established role in the growth, survival and progression of breast cancers. Here, we present a summary of a DNA microarray-based study that identified the genes that are up- and down-regulated in a human MDA-361 pLKO control clone and three clones expressing short hairpin RNA against three different regions of DCD mRNA. A list of 235 genes was differentially expressed among independent clones (> 3-fold change and P < 0.005). The gene expression of 208 was reduced and of 27 was increased in the three DCD-RNAi clones compared to pLKO control clone. The expression of 77 genes (37%) encoding for enzymes involved in amino acid metabolism, glucose metabolism and oxidoreductase activity and several genes required for cell survival and DNA repair were decreased. The expression of EGFR/ErbB-1 gene, an important predictor of outcome in breast cancer, was reduced together with the genes for betacellulin and amphiregulin, two known ligands of EGFR/ErbB receptors. Many of the 27 genes up-regulated by DCD-RNAi expression have not yet been fully characterized; among those with known function, we identified the calcium-calmodulin-dependent protein kinase-II delta and calcineurin A alpha. We compared 132 up-regulated and 12 down-regulated genes in our dataset with those genes up- and down-regulated by inhibitors targeting various signaling pathway components. The analysis showed that the genes in the DCD pathway are aligned with those functionally influenced by the drugs sirolimus, LY-294002 and wortmannin. Therefore, DCD may exert its function by activating the PI3K/AKT/mTOR signaling pathway. Together, these bioinformatic approaches suggest the involvement of DCD in the regulation of genes for breast cancer cell metabolism, proliferation and survival.

Var transcription profiling of Plasmodium falciparum 3D7: assignment of cytoadherent phenotypes to dominant transcripts

Background: Cytoadherence of Plasmodium falciparum-infected red blood cells is mediated by var gene-encoded P. falciparum erythrocyte membrane protein-1 and host receptor preference depends in most cases on which of the 50-60 var genes per genome is expressed. Enrichment of phenotypically homogenous parasites by panning on receptor expressing cells is fundamental for the identification of the corresponding var transcript. Methods: P. falciparum 3D7 parasites were panned on several transfected CHO-cell lines and their var transcripts analysed by i) reverse transcription/PCR/cloning/sequencing using a universal DBL alpha specific oligonucleotide pair and ii) by reverse transcription followed by quantitative PCR using 57 different oligonucleotide pairs. Results: Each cytoadherence selected parasite line also adhered to untransfected CHO-745 cells and upregulation of the var gene PFD995/PFD1000c was consistently associated with cytoadherence to all but one CHO cell line. In addition, parasites panned on different CHO cell lines revealed candidate var genes which reproducibly associated to the respective cytoadherent phenotype. The transcription profile obtained by RT-PCR/cloning/sequencing differed significantly from that of RT-quantitative PCR. Conclusion: Transfected CHO cell lines are of limited use for the creation of monophenotypic cytoadherent parasite lines. Nevertheless, 3D7 parasites can be reproducibly selected for the transcription of different determined var genes without genetic manipulation. Most importantly, var transcription analysis by RT-PCR/cloning/sequencing may lead to erroneous interpretation of var transcription profiles.

Functional microarray analysis suggests repressed cell-cell signaling and cell survival-related modules inhibit progression of head and neck squamous cell carcinoma

Background: Cancer shows a great diversity in its clinical behavior which cannot be easily predicted using the currently available clinical or pathological markers. The identification of pathways associated with lymph node metastasis (N+) and recurrent head and neck squamous cell carcinoma (HNSCC) may increase our understanding of the complex biology of this disease. Methods: Tumor samples were obtained from untreated HNSCC patients undergoing surgery. Patients were classified according to pathologic lymph node status (positive or negative) or tumor recurrence (recurrent or non-recurrent tumor) after treatment (surgery with neck dissection followed by radiotherapy). Using microarray gene expression, we screened tumor samples according to modules comprised by genes in the same pathway or functional category. Results: The most frequent alterations were the repression of modules in negative lymph node (N0) and in non-recurrent tumors rather than induction of modules in N+ or in recurrent tumors. N0 tumors showed repression of modules that contain cell survival genes and in non-recurrent tumors cell-cell signaling and extracellular region modules were repressed. Conclusions: The repression of modules that contain cell survival genes in N0 tumors reinforces the important role that apoptosis plays in the regulation of metastasis. In addition, because tumor samples used here were not microdissected, tumor gene expression data are represented together with the stroma, which may reveal signaling between the microenvironment and tumor cells. For instance, in non-recurrent tumors, extracellular region module was repressed, indicating that the stroma and tumor cells may have fewer interactions, which disable metastasis development. Finally, the genes highlighted in our analysis can be implicated in more than one pathway or characteristic, suggesting that therapeutic approaches to prevent tumor progression should target more than one gene or pathway, specially apoptosis and interactions between tumor cells and the stroma.

Methylation patterns revealed by MSAP profiling in genetically stable somatic embryogenic cultures of Ocotea catharinensis (Lauraceae)

Ocotea catharinensis is a rare tree species indigenous to the Atlantic rainforest of South America. In spite of its value as a hardwood species, it is in danger of extinction. The species erratically produces seeds showing irregular flowering and slow growth. Therefore, plants are not easily replaced. Tissue culture-based techniques are commonly used for obtaining living material for tree propagation and in vitro preservation. Therefore, a high-frequency somatic embryogenic system was developed for the species. In the present work, the genetic fidelity of cell aggregates and somatic embryos at various stages of in vitro development of O. catharinensis was investigated using RAPD and AFLP markers. Both analyses confirmed the absence of genetic variation in all developmental stages of O. catharinensis embryogenic cultures, verifying that the in vitro system is genetically stable. The cultures were also analyzed for their methylation profiles at 5`-CCGG-3` sites by identifying methylation-sensitive amplification polymorphisms. Some of these markers differentiated cell aggregates from embryo bodies. The sequencing of ten MSAP markers revealed that four sequences showed significant similarity to genes encoding plant proteins. Particularly, the predicted amino acid sequence of the fragment designated as OcEaggHMttc155 was similar to the enzyme 1-aminocyclopropane-1-carboxylate oxidase (ACO), which is involved in the biosynthesis of ethylene, and its expression was reported to occur from the beginning to the intermediate stages of plant embryo development. Here, we suggest that this enzyme is possibly involved in the control of the earliest stages of somatic embryogenesis of O. catharinensis, and an approach to study ACO expression during somatic embryogenesis is proposed.

SAGE transcript profiling of the juvenile cambial region of Eucalyptus grandis

Despite the importance of Eucalyptus spp. in the pulp and paper industry, functional genomic approaches have only recently been applied to understand wood formation in this genus. We attempted to establish a global view of gene expression in the juvenile cambial region of Eucalyptus grandis Hill ex Maiden. The expression profile was obtained from serial analysis of gene expression (SAGE) library data produced from 3- and 6-year-old trees. Fourteen-base expressed sequence tags (ESTs) were searched against public Eucalyptus ESTs and annotated with GenBank. Altogether 43,304 tags were generated producing 3066 unigenes with three or more copies each, 445 with a putative identity, 215 with unknown function and 2406 without an EST match. The expression profile of the juvenile cambial region revealed the presence of highly frequent transcripts related to general metabolism and energy metabolism, cellular processes, transport, structural components and information pathways. We made a quantitative analysis of a large number of genes involved in the biosynthesis of cellulose, pectin, hemicellulose and lignin. Our findings provide insight into the expression of functionally related genes involved in juvenile wood formation in young fast-growing E. grandis trees.

Transcriptional responses of host peripheral blood cells to tuberculosis infection

Host responses following exposure to Mycobacterium tuberculosis (TB) are complex and can significantly affect clinical outcome. These responses, which are largely mediated by complex immune mechanisms involving peripheral blood cells (PBCs) such as T-lymphocytes, NK cells and monocyte-derived macrophages, have not been fully characterized. We hypothesize that different clinical outcome following TB exposure will be uniquely reflected in host gene expression profiles, and expression profiling of PBCs can be used to discriminate between different TB infectious outcomes. In this study, microarray analysis was performed on PBCs from three TB groups (BCG-vaccinated, latent TB infection, and active TB infection) and a control healthy group. Supervised learning algorithms were used to identify signature genomic responses that differentiate among group samples. Gene Set Enrichment Analysis was used to determine sets of genes that were co-regulated. Multivariate permutation analysis (p < 0.01) gave 645 genes differentially expressed among the four groups, with both distinct and common patterns of gene expression observed for each group. A 127-probeset, representing 77 known genes, capable of accurately classifying samples into their respective groups was identified. In addition, 13 insulin-sensitive genes were found to be differentially regulated in all three TB infected groups, underscoring the functional association between insulin signaling pathway and TB infection. Published by Elsevier Ltd.

Transcript profiling of papaya fruit reveals differentially expressed genes associated with fruit ripening

Papaya (Carica papaya L) fruit has a short shelf life due to fast ripening induced by ethylene, but little is known about the genetic control of ripening and attributes of fruit quality. Therefore, we identified ripening-related genes affected by ethylene using cDNA-AFLP (Amplified Fragment Length Polymorphism of cDNA). Transcript profiling of non-induced and ethylene-induced fruit samples was performed, and 71 differentially expressed genes were identified. Among those genes some involved in ethylene biosynthesis, regulation of transcription, and stress responses or plant defence were found (heat shock proteins, polygalacturonase-inhibiting protein, and acyl-CoA oxidases). Several transcription factors were isolated, and except for a 14-3-3 protein, an AP2 domain-containing factor, a salt-tolerant zinc finger protein, and a suppressor of PhyA-105 1, most of them were negatively affected by ethylene, including fragments of transcripts similar to VRN1, and ethylene responsive factors (ERF). With respect to fruit quality, genes related to cell wall structure or metabolism, volatiles or pigment precursors, and vitamin biosynthesis were also found. (C) 2010 Elsevier Ireland Ltd. All rights reserved.

Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

Coordinated plant defense responses in Arabidopsis revealed by microarray analysis

Disease resistance is associated with a plant defense response that involves an integrated set of signal transduction pathways. Changes in the expression patterns of 2.375 selected genes were examined simultaneously by cDNA microarray analysis in Arabidopsis thaliana after inoculation with an incompatible fungal pathogen Alternaria brassicicola or treatment with the defense-related signaling molecules salicylic acid (SA), methyl jasmonate (MJ), or ethylene, Substantial changes (up- and down-regulation) in the steady-state abundance of 705 mRNAs were observed in response to one or more of the treatments, including known and putative defense-related genes and 106 genes with no previously described function or homology, In leaf tissue inoculated with A. brassicicola, the abundance of 168 mRNAs was increased more than 2.5-fold, whereas that of 39 mRNAs was reduced. Similarly, the abundance of 192, 221, and 55 mRNAs was highly (>2.5-fold) increased after treatment with SA, MJ, and ethylene, respectively. Data analysis revealed a surprising level of coordinated defense responses, including 169 mRNAs regulated by multiple treatments/defense pathways. The largest number of genes coinduced (one of four induced genes) and corepressed was found after treatments with SA and MJ. In addition, 50% of the genes induced by ethylene treatment were also induced by MJ treatment. These results indicated the existence of a substantial network of regulatory interactions and coordination occurring during plant defense among the different defense signaling pathways, notably between the salicylate and jasmonate pathways that were previously thought to act in an antagonistic fashion.

Gene expression profiling of papillary thyroid carcinoma identifies transcripts correlated with BRAF mutational status and lymph node metastasis

Purpose: To identify papillary thyroid carcinoma (PTC)-associated transcripts, we compared the gene expression profiles of three Serial Analysis of Gene Expression libraries generated from thyroid tumors and a normal thyroid tissue. Experimental Design: Selected transcripts were validated in a panel of 57 thyroid tumors using quantitative PCR (qPCR). An independent set of 71 paraffin-embedded sections was used for validation using immunohistochemical analysis. To determine if PTC-associated gene expression could predict lymph node involvement, a separate cohort of 130 primary PTC (54 metastatic and 76 nonmetastatic) was investigated. The BRAF(V600E) mutational status was compared with qPCR data to identify genes that might be regulated by abnormal BRAF/MEK/extracellular signal-regulated kinase signaling. Results: We identified and validated new PTC-associated transcripts. Three genes (CST6, CXCL14, and DHRS3) are strongly associated with PTC. Immunohistochemical analysis of CXCL14 confirmed the qPCR data and showed protein expression in PTC epithelial cells. We also observed that CST6, CXCL14, DHRS3, and SPP1 were associated with PTC lymph node metastasis, with CST6, CXCL14, and SPP1 being positively correlated with metastasis and DHRS3 being negatively correlated. Finally, we found a strong correlation between CST6 and CXCL14 expression and BRAF(V600E) mutational status, suggesting that these genes may be induced subsequently to BRAF activation and therefore may be downstream in the BRAF/MEK/extracellular signal-regulated kinase signaling pathway. Conclusion: CST6, CXCL14, DHRS3, and SPP1 may play a role in PTC pathogenesis and progression and are possible molecular targets for FTC therapy.

cDNA microarray analysis of differentially expressed genes in penile tissue after treatment with tadalafil

To evaluate differential gene expression in penile tissue after treatment with the phosphodiesterase 5 (PDE5) inhibitor tadalafil, as of the three clinically available PDE5 inhibitors (sildenafil, tadalafil, and vardenafil) used for the treatment of erectile dysfunction (ED), tadalafil has a long half-life and low incidence of side-effects. In all, 32 adult rats were divided into two groups. The control group received 0.5 mL of drinking water alone, while the tadalafil group was treated with tadalafil at a dose of 0.27 mg/kg. At 4 h after treatment with water or tadalafil the rats were killed and the penile tissue was removed. The total RNA was isolated from the penile tissue from both groups and differentially expressed genes were identified by cDNA microarray analysis. To validate the expression data from the microarray analysis, quantitative real-time polymerase chain reaction (PCR) and immunohistochemistry were used. In all, 153 genes were differentially expressed between the control group and the tadalafil group. We validated the microarray results by quantitative PCR for the insulin-like growth factor binding protein 6 (IGFBP-6) gene and the neuronal calcium sensor 1 (NCS-1) gene, both of which were up-regulated in the tadalafil group, and for the natriuretic peptide receptor 1 (NPR-1) gene that was down-regulated in this group. Immunohistochemistry showed localization of the NCS-1 protein in sinusoid trabeculae of the corpus cavernosum in control and tadalafil-treated rats. There was differential expression in 153 genes after tadalafil treatment. Some of these genes such as IGFBP-6, NPR-1 and NCS-1, might result in new targets in the treatment of ED.

Profiling the Proteome of the Venom from the Social Wasp Polybia paulista: A Clue to Understand the Envenoming Mechanism

The study reported here is a classical bottom-up proteomic approach where proteins from wasp venom were extracted and separated by 2-DE; the individual protein spots were proteolytically digested and subsequently identified by using tandem mass spectrometry and database query with the protein search engine MASCOT. Eighty-four venom proteins belonging to 12 different molecular functions were identified. These proteins were classified into three groups; the first is constituted of typical venom proteins: antigens-5, hyaluronidases, phospholipases, heat shock proteins, metalloproteinases, metalloproteinase-desintegrin like proteins, serine proteinases, proteinase inhibitors, vascular endothelial growth factor-related protein, arginine kinases, Sol i-II and -II like proteins, alpha-glucosidase, and superoxide dismutases. The second contained proteins structurally related to the muscles that involves the venom reservoir. The third group, associated with the housekeeping of cells from venom glands, was composed of enzymes, membrane proteins of different types, and transcriptional factors. The composition of P. paulista venom permits us to hypothesize about a general envenoming mechanism based on five actions: (i) diffusion of venom through the tissues and to the blood, (ii) tissue, (iii) hemolysis, (iv) inflammation, and (v) allergy-played by antigen-5, PLA1, hyaluronidase, HSP 60, HSP 90, and arginine kinases.