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.

Sequential analysis of global gene expression profiles in immature and in vitro matured bovine oocytes: potential molecular markers of oocyte maturation

Background: Without intensive selection, the majority of bovine oocytes submitted to in vitro embryo production (IVP) fail to develop to the blastocyst stage. This is attributed partly to their maturation status and competences. Using the Affymetrix GeneChip Bovine Genome Array, global mRNA expression analysis of immature (GV) and in vitro matured (IVM) bovine oocytes was carried out to characterize the transcriptome of bovine oocytes and then use a variety of approaches to determine whether the observed transcriptional changes during IVM was real or an artifact of the techniques used during analysis. Results: 8489 transcripts were detected across the two oocyte groups, of which similar to 25.0% (2117 transcripts) were differentially expressed (p < 0.001); corresponding to 589 over-expressed and 1528 under-expressed transcripts in the IVM oocytes compared to their immature counterparts. Over expression of transcripts by IVM oocytes is particularly interesting, therefore, a variety of approaches were employed to determine whether the observed transcriptional changes during IVM were real or an artifact of the techniques used during analysis, including the analysis of transcript abundance in oocytes in vitro matured in the presence of a-amanitin. Subsets of the differentially expressed genes were also validated by quantitative real-time PCR (qPCR) and the gene expression data was classified according to gene ontology and pathway enrichment. Numerous cell cycle linked (CDC2, CDK5, CDK8, HSPA2, MAPK14, TXNL4B), molecular transport (STX5, STX17, SEC22A, SEC22B), and differentiation (NACA) related genes were found to be among the several over-expressed transcripts in GV oocytes compared to the matured counterparts, while ANXA1, PLAU, STC1and LUM were among the over-expressed genes after oocyte maturation. Conclusion: Using sequential experiments, we have shown and confirmed transcriptional changes during oocyte maturation. This dataset provides a unique reference resource for studies concerned with the molecular mechanisms controlling oocyte meiotic maturation in cattle, addresses the existing conflicting issue of transcription during meiotic maturation and contributes to the global goal of improving assisted reproductive technology.

Identification of the Schistosoma mansoni TNF-Alpha Receptor Gene and the Effect of Human TNF-Alpha on the Parasite Gene Expression Profile

Background: Schistosoma mansoni is the major causative agent of schistosomiasis. The parasite takes advantage of host signals to complete its development in the human body. Tumor necrosis factor-alpha (TNF-alpha) is a human cytokine involved in skin inflammatory responses, and although its effect on the adult parasite's metabolism and egg-laying process has been previously described, a comprehensive assessment of the TNF-alpha pathway and its downstream molecular effects is lacking. Methodology/Principal Findings: In the present work we describe a possible TNF-alpha receptor (TNFR) homolog gene in S. mansoni (SmTNFR). SmTNFR encodes a complete receptor sequence composed of 599 amino acids, and contains four cysteine-rich domains as described for TNFR members. Real-time RT-PCR experiments revealed that SmTNFR highest expression level is in cercariae, 3.5 (+/- 0.7) times higher than in adult worms. Downstream members of the known human TNF-alpha pathway were identified by an in silico analysis, revealing a possible TNF-alpha signaling pathway in the parasite. In order to simulate parasite's exposure to human cytokine during penetration of the skin, schistosomula were exposed to human TNF-alpha just 3 h after cercariae-to-schistosomula in vitro transformation, and large-scale gene expression measurements were performed with microarrays. A total of 548 genes with significantly altered expression were detected, when compared to control parasites. In addition, treatment of adult worms with TNF-alpha caused a significantly altered expression of 1857 genes. Interestingly, the set of genes altered in adults is different from that of schistosomula, with 58 genes in common, representing 3% of altered genes in adults and 11% in 3 h-old early schistosomula. Conclusions/Significance: We describe the possible molecular elements and targets involved in human TNF-alpha effect on S. mansoni, highlighting the mechanism by which recently transformed schistosomula may sense and respond to this host mediator at the site of cercarial penetration into the skin.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

Multidimensional cluster stability analysis from a Brazilian Bradyrhizobium sp RFLP/PCR data set

The taxonomy of the N(2)-fixing bacteria belonging to the genus Bradyrhizobium is still poorly refined, mainly due to conflicting results obtained by the analysis of the phenotypic and genotypic properties. This paper presents an application of a method aiming at the identification of possible new clusters within a Brazilian collection of 119 Bradryrhizobium strains showing phenotypic characteristics of B. japonicum and B. elkanii. The stability was studied as a function of the number of restriction enzymes used in the RFLP-PCR analysis of three ribosomal regions with three restriction enzymes per region. The method proposed here uses Clustering algorithms with distances calculated by average-linkage clustering. Introducing perturbations using sub-sampling techniques makes the stability analysis. The method showed efficacy in the grouping of the species B. japonicum and B. elkanii. Furthermore, two new clusters were clearly defined, indicating possible new species, and sub-clusters within each detected cluster. (C) 2008 Elsevier B.V. All rights reserved.

A hybrid-Trefftz formulation for plane elasticity with selective enrichment of the approximations

This work is related to the so-called non-conventional finite element formulations. Essentially, a methodology for the enrichment of the initial approximation which is typical of the meshless methods and based on the clouds concept is introduced in the hybrid-Trefftz formulation for plane elasticity. The formulation presented allows for the approximation and direct enrichment of two independent fields: stresses in the domains and displacements on the boundaries of the elements. Defined by a set of elements and interior boundaries sharing a common node, the cloud notion is employed to select the enrichment support for the approximation fields. The numerical analysis performed reveals an excellent performance of the resulting formulation, characterized by the good approximation ability and a reduced computational effort. Copyright (C) 2009 John Wiley & Sons, Ltd.

Inhibition of Auxin Transport from the Ovary or from the Apical Shoot Induces Parthenocarpic Fruit-Set in Tomato Mediated by Gibberellins

Fruit-set in tomato (Solanum lycopersicum) depends on gibberellins and auxins (GAs). Here, we show, using the cv MicroTom, that application of N-1-naphthylphthalamic acid (NPA; an inhibitor of auxin transport) to unpollinated ovaries induced parthenocarpic fruit-set, associated with an increase of indole-3-acetic acid (IAA) content, and that this effect was negated by paclobutrazol (an inhibitor of GA biosynthesis). NPA-induced ovaries contained higher content of GA(1) (an active GA) and transcripts of GA biosynthetic genes (SlCPS, SlGA20ox1, and -2). Interestingly, application of NPA to pollinated ovaries prevented their growth, potentially due to supraoptimal IAA accumulation. Plant decapitation and inhibition of auxin transport by NPA from the apical shoot also induced parthenocarpic fruit growth of unpollinated ovaries. Application of IAA to the severed stump negated the plant decapitation effect, indicating that the apical shoot prevents unpollinated ovary growth through IAA transport. Parthenocarpic fruit growth induced by plant decapitation was associated with high levels of GA(1) and was counteracted by paclobutrazol treatment. Plant decapitation also produced changes in transcript levels of genes encoding enzymes of GA biosynthesis (SlCPS and SlGA20ox1) in the ovary, quite similar to those found in NPA-induced fruits. All these results suggest that auxin can have opposing effects on fruit-set, either inducing (when accumulated in the ovary) or repressing (when transported from the apical shoot) that process, and that GAs act as mediators in both cases. The effect of NPA application and decapitation on fruit-set induction was also observed in MicroTom lines bearing introgressed DWARF and SELF-PRUNING wild-type alleles.

A gene (Cargl) that regulates tissue resistance to Candida albicans maps to chromosome 14 of the mouse

Tissue susceptibility and resistance to infection with the yeast Candida albicans is genetically regulated. Analysis of the strain distribution pattern of the C. albicans resistance gene (Carg1) and additional gene and DNA segment markers in the AKXL recombinant inbred (RI) set showed that 13/15 RI strains were concordant for Carg1, Tcra and Rib1. Therefore, Carg1 is probably located within a 17 cM segment of chromosome 14, within approximately 4 cM of the other two genes. (C) 1998 Academic Press.

Cytokine-related genes identified from the RIKEN full-length mouse cDNA data set

To identify novel cytokine-related genes, we searched the set of 60,770 annotated RIKEN mouse cDNA clones (FANTOM2 clones), using keywords such as cytokine itself or cytokine names (such as interferon, interleukin, epidermal growth factor, fibroblast growth factor, and transforming growth factor). This search produced 108 known cytokines and cytokine-related products such as cytokine receptors, cytokine-associated genes, or their products (enhancers, accessory proteins, cytokine-induced genes). We found 15 clusters of FANTOM2 clones that are candidates for novel cytokine-related genes. These encoded products with strong sequence similarity to guanylate-binding protein (GBP-5), interleukin-1 receptor-associated kinase 2 (IRAK-2), interleukin 20 receptor alpha isoform 3, a member of the interferon-inducible proteins of the Ifi 200 cluster, four members of the membrane-associated family 1-8 of interferon-inducible proteins, one p27-like protein, and a hypothetical protein containing a Toll/Interleukin receptor domain. All four clones representing novel candidates of gene products from the family contain a novel highly conserved cross-species domain. Clones similar to growth factor-related products included transforming growth factor beta-inducible early growth response protein 2 (TIEG-2), TGFbeta-induced factor 2, integrin beta-like 1, latent TGF-binding protein 4S, and FGF receptor 4B. We performed a detailed sequence analysis of the candidate novel genes to elucidate their likely functional properties.

Identification of a juvenile hormone esterase-like gene in the honey bee, Apis mellifera L. - Expression analysis and functional assays

Tight control over circulating juvenile hormone (JH) levels is of prime importance in an insect`s life cycle. Consequently, enzymes involved in JH metabolism, especially juvenile hormone esterases (JHEs), play major roles during metamorphosis and reproduction. In the highly eusocial Hymenoptera, JH has been co-opted into additional functions, primarily in the development of the queen and worker castes and in age-related behavioral development of workers. Within a set of 21 carboxylesterases predicted in the honey bee genome we identified one gene (Amjhe-like) that contained the main functional motifs of insect JHEs. Its transcript levels during larval development showed a maximum at the switch from feeding to spinning behavior, coinciding with a JH titer minimum. In adult workers, the highest levels were observed in nurse bees, where a low JH titer is required to prevent the switch to foraging. Functional assays showed that Amjhe-like expression is induced by JH-III and suppressed by 20-hydroxyecdysone. RNAi-mediated silencing of Amjhe-like gene function resulted in a six-fold increase in the JH titer in adult worker bees. The temporal profile of Amjhe-like expression in larval and adult workers, the pattern of hormonal regulation and the knockdown phenotype are consistent with the function of this gene as an authentic JHE. (C) 2008 Elsevier Inc. All rights reserved.

Mutations of the KISS1 Gene in Disorders of Puberty

Context: Kisspeptin, encoded by the KISS1 gene, is a key stimulatory factor of GnRH secretion and puberty onset. Inactivating mutations of its receptor (KISS1R) cause isolated hypogonadotropic hypogonadism (IHH). A unique KISS1R-activating mutation was described in central precocious puberty (CPP). Objective: Our objective was to investigate KISS1 mutations in patients with idiopathic CPP and normosmic IHH. Patients: Eighty-three children with CPP (77 girls) and 61 patients with IHH (40 men) were studied. The control group consisted of 200 individuals with normal pubertal development. Methods: The promoter region and the three exons of KISS1 were amplified and sequenced. Cells expressing KISS1R were stimulated with synthetic human wild-type or mutant kisspeptin-54 (kp54), and inositol phosphate accumulation was measured. In a second set of experiments, kp54 was preincubated in human serum before stimulation of the cells. Results: Two novel KISS1 missense mutations, p.P74S and p.H90D, were identified in three unrelated children with idiopathic CPP. Both mutations were absent in 400 control alleles. The p.P74S mutation was identified in the heterozygous state in a boy who developed CPP at 1 yr of age. The p.H90D mutation was identified in the homozygous state in two unrelated girls with CPP. In vitro studies revealed that the capacity of the P74S and H90D mutants to stimulate IP production was similar to the wild type. After preincubation of wild-type and mutant kp54 in human serum, the capacity to stimulate signal transduction was significantly greater for P74S compared with the wild type, suggesting that the p.P74S variant is more stable. Only polymorphisms were found in the IHH group. Conclusion: Two KISS1 mutations were identified in unrelated patients with idiopathic CPP. The p.P74S variant was associated with higher kisspeptin resistance to degradation in comparison with the wild type, suggesting a role for this mutation in the precocious puberty phenotype. (J Clin Endocrinol Metab 95: 2276-2280, 2010)

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.

PVALB, a New Hurthle Adenoma Diagnostic Marker Identified through Gene Expression

Context: A better means to accurately identify malignant thyroid nodules and to distinguish them from benign tumors is needed. We previously identified markers for detecting thyroid malignancy, with sensitivity estimated at or close to 100%. One lingering problem with these markers was that false positives occurred with Hurthle cell adenomas (HCA) which lowered test specificity. Methods: To locate accurate diagnostic markers, we profiled in depth the transcripts of a HCA and a Hurthle cell carcinoma (HCC). From 1146 differentially expressed genes, 18 transcripts specifically expressed in HCA were tested by quantitative PCR in a wide range of thyroid tumors (n = 76). Sensibility and specificity were calculated using receiver operating characteristic (ROC). Selected markers were further validated in an independent set of thyroid tumors (n = 82) by immunohistochemistry. To define the panel that would yield best diagnostic accuracy, these markers were tested in combination with our previous identified markers. Results: Seventeen of the 18 genes showed statistical significance based on a mean relative level of expression (P < 0.05). KLK1 (sensitivity = 0.97) and PVALB (sensitivity = 0.94) were the best candidate markers. The combination of PVALB and C1orf24 increased specificity to > 97% and maintained sensitivity for detection of carcinoma. Conclusion: We identified tumor markers that can be used in combination for a more accurate preoperative diagnosis of thyroid nodules and for postoperative diagnosis of thyroid carcinoma in tumor sections. This improved test would help physicians rapidly focus treatment on true malignancies and avoid unnecessary treatment of benign tumors, simultaneously improving medical care and reducing costs. (J Clin Endocrinol Metab 96: E151-E160, 2011)

Affected and non-affected skin fibroblasts from systemic sclerosis patients share a gene expression profile deviated from the one observed in healthy individuals

Objectives To evaluate the gene expression profile of fibroblasts from affected and non-affected skin of systemic sclerosis (SSc) patients and from controls. Materials and methods Labeled cDNA from fibroblast cultures from forearm (affected) and axillary (non-affected) skin from six diffuse SSc patients, from three normal controls, and from MOLT-4/HEp-2/normal fibroblasts (reference pool) was probed in microarrays generated with 4193 human cDNAs from the IMAGE Consortium. Microarray images were converted into numerical data and gene expression was calculated as the ratio between fibroblast cDNA (Cy5) and reference pool cDNA (Cy3) data and analyzed by R environment/Aroma, Cluster, Tree View, and SAM softwares. Differential expression was confirmed by real time PCR for a set of selected genes. Results Eighty-eight genes were up- and 241 genes down-regulated in SSc fibroblasts. Gene expression correlation was strong between affected and non-affected fibroblast samples from the same patient (r>0.8), moderate among fibroblasts from all patients (r=0.72) and among fibroblasts from all controls (r=0.70), and modest among fibroblasts from patients and controls (r=0.55). The differential expression was confirmed by real time PCR for all selected genes. Conclusions Fibroblasts from affected and non-affected skin of SSc patients shared a similar abnormal gene expression profile, suggesting that the widespread molecular disturbance in SSc fibroblasts is more sensitive than histological and clinical alterations. Novel molecular elements potentially involved in SSc pathogenesis were identified.

Detection of Clonal Immunoglobulin and T-Cell Receptor Gene Rearrangements in Childhood Acute Lymphoblastic Leukemia Using a Low-Cost PCR Strategy

Background Imunoglobulin (Ig) and T cell receptor (TCR) gene rearrangements function as specific markers for minimal residual disease (MRD) which is one of the best predictors of outcome in childhood acute lymphoblastic leukemia (ALL) We recently reported on the prognostic value of MRD during the induction of remission through a simplified PCR method Here we report on gene rearrangement frequencies and offer guidelines for the application of the technique Procedure Two hundred thirty three children had DNA extracted from bone marrow Ig and TCR gene rearrangements were amplified using consensus primers and conventional PCR PCR products were submitted to homo/heteroduplex analysis A computer program was designed to define combinations of targets for clonal detection using a minimum set of primers and reactions Results At least one clonal marker could be detected in 98% of the patients and two markers in approximately 80% The most commonly rear ringed genes in precursor B cell ALL were IgH (75%) TCRD (59%) IgK (55%), and TCRG (54%) The most commonly rearranged genes for TALL were TCRG (100%) and TCRD (24%) The sensitivity of primers was limited to the detection of 1 leukemic cell among 100 normal cells Conclusions We propose that eight PCR reactions per ALL subtype would allow for the detection of two markers in most cases In addition these reactions ire suitable for MRD monitoring especially when aiming the selection of patients with high MRD levels (>= 10(-2)) at the end of induction therapy Such an approach would be very useful in centers with limited financial resources Pediatr Blood Cancer 2010 55 1278-1286 (C) 2010 Wiley Liss Inc