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.

Reduced transcription of TCOF1 in adult cells of Treacher Collins syndrome patients

Background: Treacher Collins syndrome (TCS) is an autosomal dominant craniofacial disorder caused by frameshift deletions or duplications in the TCOF1 gene. These mutations cause premature termination codons, which are predicted to lead to mRNA degradation by nonsense mediated mRNA decay (NMD). Haploinsufficiency of the gene product (treacle) during embryonic development is the proposed molecular mechanism underlying TCS. However, it is still unknown if TCOF1 expression levels are decreased in postembryonic human cells. Methods: We have estimated TCOF1 transcript levels through real time PCR in mRNA obtained from leucocytes and mesenchymal cells of TCS patients (n = 23) and controls (n = 18). Mutational screening and analysis of NMD were performed by direct sequencing of gDNA and cDNA, respectively. Results: All the 23 patients had typical clinical features of the syndrome and pathogenic mutations were detected in 19 of them. We demonstrated that the expression level of TCOF1 is 18-31% lower in patients than in controls (p < 0.05), even if we exclude the patients in whom we did not detect the pathogenic mutation. We also observed that the mutant allele is usually less abundant than the wild type one in mesenchymal cells. Conclusions: This is the first study to report decreased expression levels of TCOF1 in TCS adult human cells, but it is still unknown if this finding is associated to any phenotype in adulthood. In addition, as we demonstrated that alleles harboring the pathogenic mutations have lower expression, we herein corroborate the current hypothesis of NMD of the mutant transcript as the explanation for diminished levels of TCOF1 expression. Further, considering that TCOF1 deficiency in adult cells could be associated to pathologic clinical findings, it will be important to verify if TCS patients have an impairment in adult stem cell properties, as this can reduce the efficiency of plastic surgery results during rehabilitation of these patients.

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.

Increased Levels of NOTCH1, NF-kappa B, and Other Interconnected Transcription Factors Characterize Primitive Sets of Hematopoietic Stem Cells

As previously shown, higher levels of NOTCH1 and increased NF-kappa B signaling is a distinctive feature of the more primitive umbilical cord blood (UCB) CD34+ hematopoietic stem cells (HSCs), as compared to bone marrow ( BM). Differences between BM and UCB cell composition also account for this finding. The CD133 marker defines a more primitive cell subset among CD34+ HSC with a proposed hemangioblast potential. To further evaluate the molecular basis related to the more primitive characteristics of UCB and CD133+ HSC, immunomagnetically purified human CD34+ and CD133+ cells from BM and UCB were used on gene expression microarrays studies. UCB CD34+ cells contained a significantly higher proportion of CD133+ cells than BM (70% and 40%, respectively). Cluster analysis showed that BM CD133+ cells grouped with the UCB cells ( CD133+ and CD34+) rather than to BM CD34+ cells. Compared with CD34+ cells, CD133+ had a higher expression of many transcription factors (TFs). Promoter analysis on all these TF genes revealed a significantly higher frequency ( than expected by chance) of NF-kappa B-binding sites (BS), including potentially novel NF-kappa B targets such as RUNX1, GATA3, and USF1. Selected transcripts of TF related to primitive hematopoiesis and self-renewal, such as RUNX1, GATA3, USF1, TAL1, HOXA9, HOXB4, NOTCH1, RELB, and NFKB2 were evaluated by real-time PCR and were all significantly positively correlated. Taken together, our data indicate the existence of an interconnected transcriptional network characterized by higher levels of NOTCH1, NF-kappa B, and other important TFs on more primitive HSC sets.

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.

The CEBPA gene is down-regulated in acute promyelocytic leukemia and its upstream promoter, but not the core promoter, is highly methylated

impairment of CCAAT Enhancer Binding Protein alpha (CEBPA) function is a common finding in acute myeloid leukemia; nevertheless, its relevance for acute promyelocytic leukemia pathogenesis is unclear. We analyzed the expression and assessed the methylation status of the core and upstream promoters of CEBPA in acute promyelocytic leukemia at diagnosis. Patients with acute promyelocytic leukemia (n=18) presented lower levels of CEBPA expression compared to healthy controls (n=5), but higher levels than those in acute myeloid leukemia with t(8;21) (n=9) and with inv(16) (n=5). Regarding the core promoter, we detected no methylation in 39 acute promyelocytic leukemia samples or in 8 samples from controls. In contrast, analysis of the upstream promoter showed methylation in 37 of 39 samples, with 17 patients showing methylation levels over 30%. Our results corroborate data obtained in animal models showing that CEBPA is down-regulated in acute promyelocytic leukemia stem cells and suggest that epigenetic mechanisms may be involved.

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.

Group I introns and associated homing endonuclease genes reveals a clinal structure for Porphyra spiralis var. amplifolia (Bangiales, Rhodophyta) along the Eastern coast of South America

Background: Group I introns are found in the nuclear small subunit ribosomal RNA gene (SSU rDNA) of some species of the genus Porphyra (Bangiales, Rhodophyta). Size polymorphisms in group I introns has been interpreted as the result of the degeneration of homing endonuclease genes (HEG) inserted in peripheral loops of intron paired elements. In this study, intron size polymorphisms were characterized for different Porphyra spiralis var. amplifolia (PSA) populations on the Southern Brazilian coast, and were used to infer genetic relationships and genetic structure of these PSA populations, in addition to cox2-3 and rbcL-S regions. Introns of different sizes were tested qualitatively for in vitro self-splicing. Results: Five intron size polymorphisms within 17 haplotypes were obtained from 80 individuals representing eight localities along the distribution of PSA in the Eastern coast of South America. In order to infer genetic structure and genetic relationships of PSA, these polymorphisms and haplotypes were used as markers for pairwise Fst analyses, Mantel's test and median joining network. The five cox2-3 haplotypes and the unique rbcL-S haplotype were used as markers for summary statistics, neutrality tests Tajima's D and Fu's Fs and for median joining network analyses. An event of demographic expansion from a population with low effective number, followed by a pattern of isolation by distance was obtained for PSA populations with the three analyses. In vitro experiments have shown that introns of different lengths were able to self-splice from pre-RNA transcripts. Conclusion: The findings indicated that degenerated HEGs are reminiscent of the presence of a full-length and functional HEG, once fixed for PSA populations. The cline of HEG degeneration determined the pattern of isolation by distance. Analyses with the other markers indicated an event of demographic expansion from a population with low effective number. The different degrees of degeneration of the HEG do not refrain intron self-splicing. To our knowledge, this was the first study to address intraspecific evolutionary history of a nuclear group I intron; to use nuclear, mitochondrial and chloroplast DNA for population level analyses of Porphyra; and intron size polymorphism as a marker for population genetics.

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.

Characterization of global transcription profile of normal and HPV-immortalized keratinocytes and their response to TNF treatment

Background: Persistent infection by high risk HPV types (e.g. HPV-16, -18, -31, and -45) is the main risk factor for development of cervical intraepithelial neoplasia and cervical cancer. Tumor necrosis factor (TNF) is a key mediator of epithelial cell inflammatory response and exerts a potent cytostatic effect on normal or HPV16, but not on HPV18 immortalized keratinocytes. Moreover, several cervical carcinoma-derived cell lines are resistant to TNF anti-proliferative effect suggesting that the acquisition of TNF-resistance may constitute an important step in HPV-mediated carcinogenesis. In the present study, we compared the gene expression profiles of normal and HPV16 or 18 immortalized human keratinocytes before and after treatment with TNF for 3 or 60 hours. Methods: In this study, we determined the transcriptional changes 3 and 60 hours after TNF treatment of normal, HPV16 and HPV18 immortalized keratinocytes by microarray analysis. The expression pattern of two genes observed by microarray was confirmed by Northern Blot. NF-kappa B activation was also determined by electrophoretic mobility shift assay (EMSA) using specific oligonucleotides and nuclear protein extracts. Results: We observed the differential expression of a common set of genes in two TNF-sensitive cell lines that differs from those modulated in TNF-resistant ones. This information was used to define genes whose differential expression could be associated with the differential response to TNF, such as: KLK7 (kallikrein 7), SOD2 (superoxide dismutase 2), 100P (S100 calcium binding protein P), PI3 (protease inhibitor 3, skin-derived), CSTA (cystatin A), RARRES1 (retinoic acid receptor responder 1), and LXN (latexin). The differential expression of the KLK7 and SOD2 transcripts was confirmed by Northern blot. Moreover, we observed that SOD2 expression correlates with the differential NF-kappa B activation exhibited by TNF-sensitive and TNF-resistant cells. Conclusion: This is the first in depth analysis of the differential effect of TNF on normal and HPV16 or HPV18 immortalized keratinocytes. Our findings may be useful for the identification of genes involved in TNF resistance acquisition and candidate genes which deregulated expression may be associated with cervical disease establishment and/or progression.

Intercellular transport of epidermis-expressed MADS domain transcription factors and their effect on plant morphology and floral transition

P>During the lifetime of an angiosperm plant various important processes such as floral transition, specification of floral organ identity and floral determinacy, are controlled by members of the MADS domain transcription factor family. To investigate the possible non-cell-autonomous function of MADS domain proteins, we expressed GFP-tagged clones of AGAMOUS (AG), APETALA3 (AP3), PISTILLATA (PI) and SEPALLATA3 (SEP3) under the control of the MERISTEMLAYER1 promoter in Arabidopsis thaliana plants. Morphological analyses revealed that epidermal overexpression was sufficient for homeotic changes in floral organs, but that it did not result in early flowering or terminal flower phenotypes that are associated with constitutive overexpression of these proteins. Localisations of the tagged proteins in these plants were analysed with confocal laser scanning microscopy in leaf tissue, inflorescence meristems and floral meristems. We demonstrated that only AG is able to move via secondary plasmodesmata from the epidermal cell layer to the subepidermal cell layer in the floral meristem and to a lesser extent in the inflorescence meristem. To study the homeotic effects in more detail, the capacity of trafficking AG to complement the ag mutant phenotype was compared with the capacity of the non-inwards-moving AP3 protein to complement the ap3 mutant phenotype. While epidermal expression of AG gave full complementation, AP3 appeared not to be able to drive all homeotic functions from the epidermis, perhaps reflecting the difference in mobility of these proteins.

Transgenic Sweet Orange (Citrus sinensis L. Osbeck) Expressing the attacin A Gene for Resistance to Xanthomonas citri subsp citri

Genetic transformation with genes that code for antimicrobial peptides has been an important strategy used to control bacterial diseases in fruit crops, including apples, pears, and citrus. Asian citrus canker (ACC) caused by Xanthomonas citri subsp. citri Schaad et al. (Xcc) is a very destructive disease, which affects the citrus industry in most citrus-producing areas of the world. Here, we report the production of genetically transformed Natal, Pera, and Valencia sweet orange cultivars (Citrus sinensis L. Osbeck) with the insect-derived attacin A (attA) gene and the evaluation of the transgenic plants for resistance to Xcc. Agrobacterium tumefaciens Smith and Towns-mediated genetic transformation experiments involving these cultivars led to the regeneration of 23 different lines. Genetically transformed plants were identified by polymerase chain reaction, and transgene integration was confirmed by Southern blot analyses. Transcription of attA gene was detected by Northern blot analysis in all plants, except for one Natal sweet orange transformation event. Transgenic lines were multiplied by grafting onto Rangpur lime rootstock plants (Citrus limonia Osbeck) and spray-inoculated with an Xcc suspension (10(6) cfu mL(-1)). Experiments were repeated three times in a completely randomized design with seven to ten replicates. Disease severity was determined in all transgenic lines and in the control (non-transgenic) plants 30 days after inoculation. Four transgenic lines of Valencia sweet orange showed a significant reduction in disease severity caused by Xcc. These reductions ranged from 58.3% to 77.8%, corresponding to only 0.16-0.30% of leaf diseased area as opposed to 0.72% on control plants. One transgenic line of Natal sweet orange was significantly more resistant to Xcc, with a reduction of 45.2% comparing to the control plants, with only 0.14% of leaf diseased area. Genetically transformed Pera sweet orange plants expressing attA gene did not show a significant enhanced resistance to Xcc, probably due to its genetic background, which is naturally more resistant to this pathogen. The potential effect of attacin A antimicrobial peptide to control ACC may be related to the genetic background of each sweet orange cultivar regarding their natural resistance to the pathogen.

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.

Leishmania (Viannia) braziliensis transfectants overexpressing the miniexon gene lose virulence in vivo

The miniexon gene has a central role in the processing of polycistronic pre-mRNA of kinetoplastids. It is added to the 5` extremity of each mRNA, supplying the 5`-capped structure to the molecule. Previous studies in Leishmania (Leishmania) major showed that the overexpression of the miniexon array attenuates the Virulence of the parasite in in vivo assays. The results presented here extend those findings to Vionnia subgenus. Leishmania (Vionnia) braziliensis was transfected with a cosmid harboring a tandem array of one hundred miniexon gene copies and then characterized by Northern blot analysis. The overexpression of the exogenous gene was confirmed and its effect on the virulence of L (V.) braziliensis was investigated in hamsters. In BALB/c mice we could not detect parasites during the course of 15 weeks of infection. In addition, hamsters infected with transfectants overexpressing the miniexon gene exhibited only a minor footpad swelling of late onset and failed to develop progressive lesion, these attenuated parasites could be recovered from the inoculation site 1 year after infection. The persistence of parasites in the host indicates that a stable line overexpressing the miniexon may be tested as live vaccine against leishmaniasis. (C) 2008 Elsevier Ireland Ltd. All rights reserved.