A reliable measure of similarity based on dependency for short time series: an application to gene expression networks

Background: Microarray techniques have become an important tool to the investigation of genetic relationships and the assignment of different phenotypes. Since microarrays are still very expensive, most of the experiments are performed with small samples. This paper introduces a method to quantify dependency between data series composed of few sample points. The method is used to construct gene co-expression subnetworks of highly significant edges. Results: The results shown here are for an adapted subset of a Saccharomyces cerevisiae gene expression data set with low temporal resolution and poor statistics. The method reveals common transcription factors with a high confidence level and allows the construction of subnetworks with high biological relevance that reveals characteristic features of the processes driving the organism adaptations to specific environmental conditions. Conclusion: Our method allows a reliable and sophisticated analysis of microarray data even under severe constraints. The utilization of systems biology improves the biologists ability to elucidate the mechanisms underlying celular processes and to formulate new hypotheses.

Gene Expression Profiles in Radiation Workers Occupationally Exposed to Ionizing Radiation

Ionizing radiation OR) imposes risks to human health and the environment. IR at low doses and low (lose rates has the potency to initiate carcinogenesis. Genotoxic environmental agents such as IR trigger a cascade of signal transduction pathways for cellular protection. In this study, using cDNA microarray technique, we monitored the gene expression profiles in lymphocytes derived from radiation-ex posed individuals (radiation workers). Physical dosimetry records on these patients indicated that the absorbed dose ranged from 0.696 to 39.088 mSv. Gene expression analysis revealed statistically significant transcriptional changes in a total of 78 genes (21 up-regulated and 57 clown-regulated) involved in several biological processes such as ubiquitin cycle (UHRF2 and PIAS1), DNA repair (LIG3, XPA, ERCC5, RAD52, DCLRE1C), cell cycle regulation/proliferation (RHOA, CABLES2, TGFB2, IL16), and stress response (GSTP1, PPP2R5A, DUSP22). Some of the genes that showed altered expression profiles in this study call be used as biomarkers for monitoring the chronic low level exposure in humans. Additionally, alterations in gene expression patterns observed in chronically exposed radiation workers reinforces the need for defining the effective radiation dose that causes immediate genetic damage as well as the long-term effects on genomic instability, including cancer.

Global gene expression under nitrogen starvation in Xylella fastidiosa: contribution of the sigma(54) regulon

Background: Xylella fastidiosa, a Gram-negative fastidious bacterium, grows in the xylem of several plants causing diseases such as citrus variegated chlorosis. As the xylem sap contains low concentrations of amino acids and other compounds, X. fastidiosa needs to cope with nitrogen limitation in its natural habitat. Results: In this work, we performed a whole-genome microarray analysis of the X. fastidiosa nitrogen starvation response. A time course experiment (2, 8 and 12 hours) of cultures grown in defined medium under nitrogen starvation revealed many differentially expressed genes, such as those related to transport, nitrogen assimilation, amino acid biosynthesis, transcriptional regulation, and many genes encoding hypothetical proteins. In addition, a decrease in the expression levels of many genes involved in carbon metabolism and energy generation pathways was also observed. Comparison of gene expression profiles between the wild type strain and the rpoN null mutant allowed the identification of genes directly or indirectly induced by nitrogen starvation in a sigma(54)-dependent manner. A more complete picture of the sigma(54) regulon was achieved by combining the transcriptome data with an in silico search for potential sigma(54)-dependent promoters, using a position weight matrix approach. One of these sigma(54)-predicted binding sites, located upstream of the glnA gene (encoding glutamine synthetase), was validated by primer extension assays, confirming that this gene has a sigma(54)-dependent promoter. Conclusions: Together, these results show that nitrogen starvation causes intense changes in the X. fastidiosa transcriptome and some of these differentially expressed genes belong to the sigma(54) regulon.

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.

Global gene expression profile in myelodysplastic syndromes using SAGE

The molecular pathogenesis of myelodysplastic syndromes (MDS) is poorly understood. In order to expand our knowledge of genetic defects in MDS, we determined the overall profile of genes expressed in bone marrow from patients with refractory anemia with excess blasts ( RAEB) by serial analysis of gene expression ( SAGE). The present report describes a partial transcriptome of RAEB bone marrow derived from 56,694 sequenced tags that provides information about expressed gene products. This is the first attempt to determine an overall profile of gene expression specifically in RAEB at diagnosis using SAGE, which should be useful in the understanding of the physiopathology of MDS and in identifying the genes involved.

Hormone-regulated expression and distribution of versican in mouse uterine tissues

Background: Remodeling of the extracellular matrix is one of the most striking features observed in the uterus during the estrous cycle and after hormone replacement. Versican (VER) is a hyaluronan-binding proteoglycan that undergoes RNA alternative splicing, generating four distinct isoforms. This study analyzed the synthesis and distribution of VER in mouse uterine tissues during the estrous cycle, in ovariectomized (OVX) animals and after 17beta-estradiol (E2) and medroxyprogesterone (MPA) treatments, either alone or in combination. Methods: Uteri from mice in all phases of the estrous cycle, and animals subjected to ovariectomy and hormone replacement were collected for immunoperoxidase staining for versican, as well as PCR and quantitative Real Time PCR. Results: In diestrus and proestrus, VER was exclusively expressed in the endometrial stroma. In estrus and metaestrus, VER was present in both endometrial stroma and myometrium. In OVX mice, VER immunoreaction was abolished in all uterine tissues. VER expression was restored by E2, MPA and E2+MPA treatments. Real Time PCR analysis showed that VER expression increases considerably in the MPA-treated group. Analysis of mRNA identified isoforms V0, V1 and V3 in the mouse uterus. Conclusion: These results show that the expression of versican in uterine tissues is modulated by ovarian steroid hormones, in a tissue-specific manner. VER is induced in the myometrium exclusively by E2, whereas MPA induces VER deposition only in the endometrial stroma.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

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.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Immunological responses and cytokine gene expression analysis to Cooperia punctata infections in resistant and susceptible Nelore cattle

Cellular and humoral immune response, as well as cytokine gene expression, was assessed in Nelore cattle with different degrees of resistance to Cooperia punctata natural infection. One hundred cattle (male, weaned, 11-12 months old), kept together on pasture, were evaluated. Faecal and blood samples were collected for parasitological and immunological assays. Based on nematode faecal egg counts (FEC) and worm burden, the seven most resistant and the eight most susceptible animals were selected. Tissue samples of the small intestine were collected for histological quantification of inflammatory cells and analysis of cytokine gene expression (IL-2, IL-4, IL-8, IL-1 2p35, IL-13, TNF-alpha, IFN-gamma, MCP-1, MCP-2, and MUC- 1) using real-time RT-PCR. Mucus samples were also collected for IgA levels determination. Serum IgG1 mean levels against C. punctata antigens were higher in the resistant group, but significant differences between groups were only observed 14 days after the beginning of the experiment against infective larvae (1-3) and 14 and 84 days against adult antigens. The resistant group also presented higher IgA levels against C. punctata (L3 and adult) antigens with significant difference 14 days after the beginning of the trial (P < 0.05). In the small-intestine mucosa, levels of IgA anti-L3 and anti-adult C. punctata were higher in the resistant group, compared with the susceptible group (P < 0.05). Gene expression of both T(H)2 cytokines (IL-4 and IL-13) in the resistant group and T(H)1 cytokines (IL-2, IL-1 2p35, IFN-gamma and MCP-1) in the susceptible group was up-regulated. Such results suggested that immune response to C. punctata was probably mediated by TH2 cytokines in the resistant group and by T(H)1 cytokines in the susceptible group. (C) 2008 Elsevier B.V. All rights reserved.

Enteroinvasive Escherichia coli vs. Shigella flexneri : how different patterns of gene expression affect virulence

Important features of the enteroinvasive Escherichia coli (EIEC) phenotype and gene expression likely to confer EIEC with a lower ability to cause disease than Shigella flexneri were described here for the first time. To confirm the lower pathogenicity of EIEC, we have analyzed the keratoconjunctivitis developed in guinea-pigs with EIEC or S. flexneri. Shigella flexneri induced a more pronounced proinflammatory response, whereas EIEC induced a mild form of the disease. EIEC showed a significantly less efficient cell-to-cell Caco-2 dissemination when compared with S. flexneri. Plaques formed by EIEC during intercellular spreading were four times smaller than those formed by S. flexneri. At the molecular level, the lower expression of virulence genes by EIEC during infection of Caco-2 cells highlighted the importance of effective gene transcription for bacterial pathogenicity.

Atorvastatin effects on SREBF1a and SCAP gene expression in mononuclear cells and its relation with lowering-lipids response

Background: The transcription factors SREBP1 and SCAP are involved in intracellular cholesterol homeostasis. Polymorphisms of these genes have been associated with variations on serum lipid levels and response to statins that are potent cholesterol-lowering drugs. We evaluated the effects of atorvastatin on SREBF1a and SCAP mRNA expression in peripheral blood mononuclear cells (PBMC) and a possible association with gene polymorphisms and lowering-cholesterol response. Methods: Fifty-nine hypercholesterolemic patients were treated with atorvastatin (10 mg/day for 4 weeks). Serum lipid profile and mRNA expression in PBMC were assessed before and after the treatment. Gene expression was quantified by real-time PCR using GAPD as endogenous reference and mRNA expression in HepG2 cells as calibrator. SREBF1 -36delG and SCAP A2386G polymorphisms were detected by PCR-RFLP. Results: Our results showed that transcription of SREBF1a and SCAP was coordinately regulated by atorvastatin (r=0.595, p<0.001), and that reduction in SCAP transcription was associated with the 2386AA genotype (p=0.019). Individuals who responded to atorvastatin with a downregulation of SCAP had also a lower triglyceride compared to those who responded to atorvastatin with an upregulation of SCAP. Conclusion: Atorvastatin has differential effects on SREBF1a and SCAP mRNA expression in PBMC that are associated with baseline transcription levels, triglycerides response to atorvastatin and SCAP A2386G polymorphism. (c) 2008 Elsevier B.V. All rights reserved.

EIN3-like gene expression during fruit ripening of Cavendish banana (Musa acuminata cv. Grande naine)

Ethylene signal transduction initiates with ethylene binding at receptor proteins and terminates in a transcription cascade involving the EIN3/EIL transcription factors. Here, we have isolated four cDNAs homologs of the Arabidopsis EIN3/EIN3-like gene, MA-EILs (Musa acuminata ethylene insensitive 3-like) from banana fruit. Sequence comparison with other banana EIL gene already registered in the database led us to conclude that, at this day, at least five different genes namely MA-EIL1, MA-EIL2/AB266318, MA-EIL3/AB266319, MA-EIL4/AB266320 and AB266321 exist in banana. Phylogenetic analyses included all banana EIL genes within a same cluster consisting of rice OsEILs, a monocotyledonous plant as banana. However, MA-EIL1, MA-EIL2/AB266318, MA-EIL4/AB266320 and AB266321 on one side, and MA-EIL3/AB266319 on the other side, belong to two distant subclusters. MA-EIL mRNAs were detected in all examined banana tissues but at lower level in peel than in pulp. According to tissues, MA-EIL genes were differentially regulated by ripening and ethylene in mature green fruit and wounding in old and young leaves. MA-EIL2/AB266318 was the unique ripening- and ethylene-induced gene; MA-EIL1, MA-EIL4/Ab266320 and AB266321 genes were downregulated, while MA-EIL3/AB266319 presented an unusual pattern of expression. Interestingly, a marked change was observed mainly in MA-EIL1 and MA-EIL3/Ab266319 mRNA accumulation concomitantly with changes in ethylene responsiveness of fruit. Upon wounding, the main effect was observed in MA-EIL4/AB266320 and AB266321 mRNA levels, which presented a markedly increase in both young and old leaves, respectively. Data presented in this study suggest the importance of a transcriptionally step control in the regulation of EIL genes during banana fruit ripening.

Gene expression analysis of Paracoccidioides brasiliensis transition from conidium to yeast cell

Paracoccidioides brasiliensis infectious process relies on the initial expression of virulence faactors that are assumed to be controlled by molecular mechanisms through which the conidia and/or mycelial fragments convert to yeast cells. In order to analyze the profile of the thermally-induced dimorphic gene expression, 48 h C-L transition cultures which had been incubated at 36 degrees C were studied. By this time approximately 50% of the conidial population had already reverted to yeast form cells. At this transition time, an EST-Orestes library was constructed and characterized. As a result, 79 sequences were obtained, of which 39 (49.4%) had not been described previously in other libraries of this fungus and which could represent novel exclusive C-Y transition genes. Two of these sequences are, among others, cholestanol delta-isomerase, and electron transfer flavoprotein-ubiquinoneoxidoreductase (ETF-QO). The other 40 (50.6%) sequences were shared with Mycelia (M), Yeast (Y) or Mycelia to yest transition (M-Y) libraries. An important component of this group of sequences is a putative response regulator receiver SKN7, a protein of high importance in stress adaptation and a regulator of virulence in some bacteria and fungi. This is the first report identifying genes expressed during the C-Y transition process, the initial step required to understand the natural history of P brasiliensis conidia induced infection.