Juan Carlos Oliveros

Genome-wide mapping of Arabidopsis thaliana origins of DNA replication and their associated epigenetic marks

Genome integrity requires faithful chromosome duplication. Origins of replication, the genomic sites at which DNA replication initiates, are scattered throughout the genome. Their mapping at a genomic scale in multicellular organisms has been challenging. In this study we profiled origins in Arabidopsis thaliana by high-throughput sequencing of newly synthesized DNA and identified ~1,500 putative origins genome-wide. This was supported by chromatin immunoprecipitation and microarray (ChIP-chip) experiments to identify ORC1- and CDC6-binding sites. We validated origin activity independently by measuring the abundance of nascent DNA strands. The midpoints of most A. thaliana origin regions are preferentially located within the 5′ half of genes, enriched in G+C, histone H2A.Z, H3K4me2, H3K4me3 and H4K5ac, and depleted in H3K4me1 and H3K9me2. Our data help clarify the epigenetic specification of DNA replication origins in A. thaliana and have implications for other eukaryotes.

Improved protein‐binding microarrays for the identification of DNA‐binding specificities of transcription factors

Transcriptional regulation depends on the specificity of transcription factors (TFs) recognizing cis regulatory sequences in the promoters of target genes. Current knowledge about DNA-binding specificities of TFs is based mostly on low- to medium-throughput methodologies, revealing DNA motifs bound by a TF with high affinity. These strategies are time-consuming and often fail to identify DNA motifs recognized by a TF with lower affinity but retaining biological relevance. Here we report on the development of a protein-binding microarray (PBM11) containing all possible double-stranded 11-mers for the determination of DNA-binding specificities of TFs. The large number of sequences in the PBM11 allows accurate and high-throughput quantification of TF-binding sites, outperforming previous methods. We applied this tool to determine binding site specificities of two Arabidopsis TFs, MYC2 and ERF1, rendering the G-box and the GCC-box, respectively, as their highest-affinity binding sites. In addition, we identified variants of the G-box recognized by MYC2 with high and medium affinity, whereas ERF1 only recognized GCC variants with low affinity, indicating that ERF1 binding to DNA has stricter base requirements than MYC2. Analysis of transcriptomic data revealed that high- and medium-affinity binding sites have biological significance, probably representing relevant cis-acting elements in vivo. Comparison of promoter sequences with putative orthologs from closely related species demonstrated a high degree of conservation of all the identified DNA elements. The combination of PBM11, transcriptomic data and phylogenomic footprinting provides a straightforward method for the prediction of biologically active cis-elements, and thus for identification of in vivo DNA targets of TFs.

Mining functional information associated with expression arrays

Abstract. Deciphering the networks of interactions between molecules in biological systems has gained momentum with the monitoring of gene expression patterns at the genomic scale. Expression array experiments provide vast amounts of experimental data about these networks, the analysis of which requires new computational methods. In particular, issues related to the extraction of biological information are key for the end users. We propose here a strategy, implemented in a system called GEISHA (gene expression information system for human analysis) and able to detect biological terms significantly associated to different gene expression clusters by mining collections of Medline abstracts. GEISHA is based on a comparison of the frequency of abstracts linked to different gene clusters and containing a given term. Interpretation by the end user of the biological meaning of the terms is facilitated by embedding them in the corresponding significant sentences and abstracts and by establishing relations with other, equally significant terms. The information provided by GEISHA for the available yeast expression data compares favorably with the functional annotations provided by human experts, demonstrating the potential value of GEISHA as an assistant for the analysis of expression array experiments.

Berry Flesh and Skin Ripening Features in Vitis vinifera as Assessed by Transcriptional Profiling

Background Ripening of fleshy fruit is a complex developmental process involving the differentiation of tissues with separate functions. During grapevine berry ripening important processes contributing to table and wine grape quality take place, some of them flesh- or skin-specific. In this study, transcriptional profiles throughout flesh and skin ripening were followed during two different seasons in a table grape cultivar ‘Muscat Hamburg’ to determine tissue-specific as well as common developmental programs. Methodology/Principal Findings Using an updated GrapeGen Affymetrix GeneChip® annotation based on grapevine 12×v1 gene predictions, 2188 differentially accumulated transcripts between flesh and skin and 2839 transcripts differentially accumulated throughout ripening in the same manner in both tissues were identified. Transcriptional profiles were dominated by changes at the beginning of veraison which affect both pericarp tissues, although frequently delayed or with lower intensity in the skin than in the flesh. Functional enrichment analysis identified the decay on biosynthetic processes, photosynthesis and transport as a major part of the program delayed in the skin. In addition, a higher number of functional categories, including several related to macromolecule transport and phenylpropanoid and lipid biosynthesis, were over-represented in transcripts accumulated to higher levels in the skin. Functional enrichment also indicated auxin, gibberellins and bHLH transcription factors to take part in the regulation of pre-veraison processes in the pericarp, whereas WRKY and C2H2 family transcription factors seems to more specifically participate in the regulation of skin and flesh ripening, respectively. Conclusions/Significance A transcriptomic analysis indicates that a large part of the ripening program is shared by both pericarp tissues despite some components are delayed in the skin. In addition, important tissue differences are present from early stages prior to the ripening onset including tissue-specific regulators. Altogether, these findings provide key elements to understand berry ripening and its differential regulation in flesh and skin.

Distinct and conserved transcriptomic changes during nematode‐induced giant cell development in tomato compared with Arabidopsis: a functional role for gene repression

Root-knot nematodes (RKNs) induce giant cells (GCs) from root vascular cells inside the galls. Accompanying molecular changes as a function of infection time and across different species, and their functional impact, are still poorly understood. Thus, the transcriptomes of tomato galls and laser capture microdissected (LCM) GCs over the course of parasitism were compared with those of Arabidopsis, and functional analysis of a repressed gene was performed. Microarray hybridization with RNA from galls and LCM GCs, infection-reproduction tests and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) transcriptional profiles in susceptible and resistant (Mi-1) lines were performed in tomato. Tomato GC-induced genes include some possibly contributing to the epigenetic control of GC identity. GC-repressed genes are conserved between tomato and Arabidopsis, notably those involved in lignin deposition. However, genes related to the regulation of gene expression diverge, suggesting that diverse transcriptional regulators mediate common responses leading to GC formation in different plant species. TPX1, a cell wall peroxidase specifically involved in lignification, was strongly repressed in GCs/galls, but induced in a nearly isogenic Mi-1 resistant line on nematode infection. TPX1 overexpression in susceptible plants hindered nematode reproduction and GC expansion. Time-course and cross-species comparisons of gall and GC transcriptomes provide novel insights pointing to the relevance of gene repression during RKN establishment.

Severe acute respiratory syndrome coronavirus envelope protein regulates cell stress response and apoptosis.

Severe acute respiratory syndrome virus (SARS-CoV) that lacks the envelope (E) gene (rSARS-CoV-ΔE) is attenuated in vivo. To identify factors that contribute to rSARS-CoV-ΔE attenuation, gene expression in cells infected by SARS-CoV with or without E gene was compared. Twenty-five stress response genes were preferentially upregulated during infection in the absence of the E gene. In addition, genes involved in signal transduction, transcription, cell metabolism, immunoregulation, inflammation, apoptosis and cell cycle and differentiation were differentially regulated in cells infected with rSARS-CoV with or without the E gene. Administration of E protein in trans reduced the stress response in cells infected with rSARS-CoV-ΔE or with respiratory syncytial virus, or treated with drugs, such as tunicamycin and thapsigargin that elicit cell stress by different mechanisms. In addition, SARS-CoV E protein down-regulated the signaling pathway inositol-requiring enzyme 1 (IRE-1) of the unfolded protein response, but not the PKR-like ER kinase (PERK) or activating transcription factor 6 (ATF-6) pathways, and reduced cell apoptosis. Overall, the activation of the IRE-1 pathway was not able to restore cell homeostasis, and apoptosis was induced probably as a measure to protect the host by limiting virus production and dissemination. The expression of proinflammatory cytokines was reduced in rSARS-CoV-ΔE-infected cells compared to rSARS-CoV-infected cells, suggesting that the increase in stress responses and the reduction of inflammation in the absence of the E gene contributed to the attenuation of rSARS-CoV-ΔE.

A freezing-sensitive mutant of Arabidopsis , frs1 , is a new aba3 allele

Abstract. To investigate the molecular mechanisms controlling the process of cold acclimation and to identify genes involved in plant freezing tolerance, mutations that impaired the cold acclimation capability of Arabidopsis thaliana (L.) Heynh. were screened for. A new mutation, frs1 (freezing sensitive 1), that reduced both the constitutive freezing tolerance as well as the freezing tolerance of Arabidopsis after cold acclimation was characterized. This mutation also produced a wilty phenotype and excessive water loss. Plants with the frs1 mutation recovered their wild-type phenotype, their capability to tolerate freezing temperatures and their capability to retain water after an exogenous abscisic acid (ABA) treatment. Measurements of ABA revealed that frs1 mutants were ABA deficient, and complementation tests indicated that frs1 mutation was a new allele of the ABA3 locus showing that a mutation in this locus leads to an impairment of freezing tolerance. These results constitute the first report showing that a mutation in ABA3 leads to an impairment of freezing tolerance, and not only strengthen the conclusion that ABA is required for full development of freezing tolerance in cold-acclimated plants, but also demonstrate that ABA mediates the constitutive freezing tolerance of Arabidopsis. Gene expression in frs1 mutants was altered in response to dehydration, suggesting that freezing tolerance in Arabidopsis depends on ABA-regulated proteins that allow plants to survive the challenges imposed by subzero temperatures, mainly freeze-induced cellular dehydration.

Analysis of Genome-Wide Changes in the Translatome of Arabidopsis Seedlings Subjected to Heat Stress

Heat stress is one of the most prominent and deleterious environmental threats affecting plant growth and development. Upon high temperatures, plants launch specialized gene expression programs that promote stress protection and survival. These programs involve global and specific changes at the transcriptional and translational levels. However, the coordination of these processes and their specific role in the establishment of the heat stress response is not fully elucidated. We have carried out a genome-wide analysis to monitor the changes in the translation efficiency of individual mRNAs of Arabidopsis thaliana seedlings after the exposure to a heat shock stress. Our results demonstrate that translation exerts a wide but dual regulation of gene expression. For the majority of mRNAs, translation is severely repressed, causing a decreased of 50% in the association of the bulk of mRNAs to polysomes. However, some relevant mRNAs involved in different aspects of homeostasis maintenance follow a differential pattern of translation. Sequence analyses of the differentially translated mRNAs unravels that some features, such as the 5′UTR G+C content and the cDNA length, may take part in the discrimination mechanisms for mRNA polysome loading. Among the differentially translated genes, master regulators of the stress response stand out, highlighting the main role of translation in the early establishment of the physiological response of plants to elevated temperatures.

Bioinformatics methods for the analysis of expression arrays: data clustering and information extraction

Expression arrays facilitate the monitoring of changes in the expression patterns of large collections of genes. The analysis of expression array data has become a computationally-intensive task that requires the development of bioinformatics technology for a number of key stages in the process, such as image analysis, database storage, gene clustering and information extraction. Here, we review the current trends in each of these areas, with particular emphasis on the development of the related technology being carried out within our groups.

Transcriptomic fingerprinting of Pseudomonas putida under alternative physiological regimes

Pseudomonas putida KT2440 is a metabolically versatile soil bacterium useful both as a model biodegradative organism and as a host of catalytic activities of biotechnological interest. In this report, we present the high-resolution transcriptome of P. putida cultured on different carbon sources as revealed by deep sequencing of the corresponding RNA pools. Examination of the data from growth on substrates that are processed through distinct pathways (glucose, fructose, succinate and glycerol) revealed that ≥ 20% of the P. putida genome is differentially expressed depending on the ensuing physiological regime. Changes affected not only metabolic genes but also a suite of global regulators, e.g. the rpoS sigma subunit of RNA polymerase, various cold-shock proteins and the three HU histone-like proteins. Specifically, the genes encoding HU subunit variants hupA, hupB and hupN drastically altered their expression levels (and thus their ability to form heterodimeric combinations) under the diverse growth conditions. Furthermore, we found that two small RNAs, crcZ and crcY, known to inhibit the Crc protein that mediates catabolite repression in P. putida, were both down-regulated by glucose. The raw transcriptomic data generated in this work is made available to the community through the Gene Expression Omnibus database.