Identifying the connections between molecular and physiological processes underlying the diversity of drought stress responses in plants is key for basic and applied science. Drought stress response involves a large number of molecular pathways and subsequent physiological processes. Therefore, it constitutes an archetypical systems biology model. We first inferred a gene-phenotype network exploiting differences in drought responses of eight sunflower (Helianthus annuus) genotypes to two drought stress scenarios. Large transcriptomic data were obtained with the sunflower Affymetrix microarray, comprising 32423 probesets, and were associated to nine morpho-physiological traits (integrated transpired water, leaf transpiration rate, osmotic potential, relative water content, leaf mass per area, carbon isotope discrimination, plant height, number of leaves and collar diameter) using sPLS regression. Overall, we could associate the expression patterns of 1263 probesets to six phenotypic traits and identify if correlations were due to treatment, genotype and/or their interaction. We also identified genes whose expression is affected at moderate and/or intense drought stress together with genes whose expression variation could explain phenotypic and drought tolerance variability among our genetic material. We then used the network model to study phenotypic changes in less tractable agronomical conditions, i.e. sunflower hybrids subjected to different watering regimes in field trials. Mapping this new dataset in the gene-phenotype network allowed us to identify genes whose expression was robustly affected by water deprivation in both controlled and field conditions. The enrichment in genes correlated to relative water content and osmotic potential provides evidence of the importance of these traits in agronomical conditions.

Profiling of DNA and histone modifications has recently allowed the establishment of reference epigenomes from several model organisms. This identified a major chromatin state for active genes that contains monoubiquitinated H2B (H2Bub), a mark linked to transcription elongation. However, assessment of dynamic chromatin changes during the reprogramming of gene expression in response to extrinsic or developmental signals has been more difficult. Here we used the major developmental switch that Arabidopsis thaliana plants undergo upon their initial perception of light, known as photomorphogenesis, as a paradigm to assess spatial and temporal dynamics of monoubiquitinated H2B (H2Bub) and its impact on transcriptional responses. The process involves rapid and extensive transcriptional reprogramming and represents a developmental window well suited to studying cell division–independent chromatin changes. Genome-wide H2Bub distribution was determined together with transcriptome profiles at three time points during early photomorphogenesis. This revealed de novo marking of 177 genes upon the first hour of illumination, illustrating the dynamic nature of H2Bub enrichment in a genomic context. Gene upregulation was associated with H2Bub enrichment, while H2Bub levels generally remained stable during gene downregulation. We further report that H2Bub influences the modulation of gene expression, as both gene up- and downregulation were globally weaker in hub1 mutant plants that lack H2Bub. H2Bub-dependent regulation notably impacted genes with fast and transient light induction, and several circadian clock components whose mRNA levels are tightly regulated by sharp oscillations. Based on these findings, we propose that H2B monoubiquitination is part of a transcription-coupled, chromatin-based mechanism to rapidly modulate gene expression.

Author Summary

In eukaryotes, chromatin-based mechanisms overlay with DNA sequence information to determine the transcriptional output of the genome. Evaluating the role of chromatin state variations in the regulation of gene expression is therefore key to understanding their contribution to development. Several transcriptional coactivators contribute to the selective regulation of cellular pathways by coordinating histone H2B monoubiquitination (H2Bub) with other histone modifications. Although H2Bub is present on a large number of genes, its loss was shown to affect RNA levels for only a small subset of genes, and therefore its influence on gene expression is not well understood. Here we assessed the impact of H2Bub on expression changes during a rapid developmental transition that initiates upon exposure of plants to light. This revealed that H2Bub marking is highly dynamic in a genomic context. Furthermore, a large repertoire of light-responsive genes was impaired for rapid up- or downregulation, indicating that H2Bub is important for attaining appropriate expression levels. Regulatory factors and circadian clock components are well represented within the set of genes impacted by H2Bub dynamics for rapid changes in RNA levels, indicating that some genes whose mRNAs need tight and rapid control are particularly sensitive to chromatin-based mechanisms linked to H2Bub deposition.

Background

Expression traits can vary quantitatively between individuals and have a complex inheritance. Identification of the genetics underlying transcript variation can help in the understanding of phenotypic variation due to genetic factors regulating transcript abundance and shed light into divergence patterns. So far, only a limited number of studies have addressed this subject in Arabidopsis, with contrasting results due to dissimilar statistical power. Here, we present the transcriptome architecture in leaf tissue of two RIL sets obtained from a connected-cross design involving 3 commonly used accessions. We also present the transcriptome architecture observed in developing seeds of a third independent cross.

Results

The utilisation of the novel R/eqtl package (which goal is to automatize and extend functions from the R/qtl package) allowed us to map 4,290 and 6,534 eQTLs in the Cvi-0 × Col-0 and Bur-0 × Col-0 recombinant populations respectively. In agreement with previous studies, we observed a larger phenotypic variance explained by eQTLs in linkage with the controlled gene (potentially cis-acting), compared to distant loci (acting necessarily indirectly or in trans). Distant eQTLs hotspots were essentially not conserved between crosses, but instead, cross-specific. Accounting for confounding factors using a probabilistic approach (VBQTL) increased the mapping resolution and the number of significant associations. Moreover, using local eQTLs obtained from this approach, we detected evidence for a directional allelic effect in genes with related function, where significantly more eQTLs than expected by chance were up-regulated from one of the accessions. Primary experimental data, analysis parameters, eQTL results and visualisation of LOD score curves presented here are stored and accessible through the QTLstore service database http://qtlstore.versailles.inra.fr/.

Conclusions

Our results demonstrate the extensive diversity and moderately conserved eQTL landscape between crosses and validate the utilisation of expression traits to explore for candidates behind phenotypic variation among accessions. Furthermore, this stresses the need for a wider spectrum of diversity to fully understand expression trait variation within a species.

Cultivated for centuries, the varieties of rose have been selected based on a number of flower traits. Understanding the genetic and molecular basis that contributes to these traits will impact on future improvements for this economically important ornamental plant. In this study, we used scanning electron microscopy and sections of meristems and flowers to establish a precise morphological calendar from early rose flower development stages to senescing flowers. Global gene expression was investigated from floral meristem initiation up to flower senescence in three rose genotypes exhibiting contrasted floral traits including continuous versus once flowering and simple versus double flower architecture, using a newly developed Affymetrix microarray (Rosa1_Affyarray) tool containing sequences representing 4765 unigenes expressed during flower development. Data analyses permitted the identification of genes associated with floral transition, floral organs initiation up to flower senescence. Quantitative real time PCR analyses validated the mRNA accumulation changes observed in microarray hybridizations for a selection of 24 genes expressed at either high or low levels. Our data describe the early flower development stages in Rosa sp, the production of a rose microarray and demonstrate its usefulness and reliability to study gene expression during extensive development phases, from the vegetative meristem to the senescent flower.

Background

Whole genome duplication is a common evolutionary event in plants. Bread wheat (Triticum aestivum L.) is a good model to investigate the impact of paleo- and neoduplications on the organization and function of modern plant genomes.

Results

We performed an RNA sequencing-based inference of the grain filling gene network in bread wheat and identified a set of 37,695 non-redundant sequence clusters, which is an unprecedented resolution corresponding to an estimated half of the wheat genome unigene repertoire. Using the Brachypodium distachyon genome as a reference for the Triticeae, we classified gene clusters into orthologous, paralogous, and homoeologous relationships. Based on this wheat gene evolutionary classification, older duplicated copies (dating back 50 to 70 million years) exhibit more than 80% gene loss and expression divergence while recent duplicates (dating back 1.5 to 3 million years) show only 54% gene loss and 36 to 49% expression divergence.

Conclusions

We suggest that structural shuffling due to duplicated gene loss is a rapid process, whereas functional shuffling due to neo- and/or subfunctionalization of duplicates is a longer process, and that both shuffling mechanisms drive functional redundancy erosion. We conclude that, as a result of these mechanisms, half the gene duplicates in plants are structurally and functionally altered within 10 million years of evolution, and the diploidization process is completed after 45 to 50 million years following polyploidization.

Cinnamoyl-CoA reductase (CCR), which catalyses the first committed step of the lignin-specific branch of monolignol biosynthesis, has been extensively characterized in dicot species, but few data are available in monocots. By screening a Mu insertional mutant collection in maize, a mutant in the CCR1 gene was isolated named Zmccr1–. In this mutant, CCR1 gene expression is reduced to 31% of the residual wild-type level. Zmccr1– exhibited enhanced digestibility without compromising plant growth and development. Lignin analysis revealed a slight decrease in lignin content and significant changes in lignin structure. p-Hydroxyphenyl units were strongly decreased and the syringyl/guaiacyl ratio was slightly increased. At the cellular level, alterations in lignin deposition were mainly observed in the walls of the sclerenchymatic fibre cells surrounding the vascular bundles. These cell walls showed little to no staining with phloroglucinol. These histochemical changes were accompanied by an increase in sclerenchyma surface area and an alteration in cell shape. In keeping with this cell type-specific phenotype, transcriptomics performed at an early stage of plant development revealed the down-regulation of genes specifically associated with fibre wall formation. To the present authors’ knowledge, this is the first functional characterization of CCR1 in a grass species.

Background

Comparative genomics has emerged as a promising means of unravelling the molecular networks underlying complex traits such as drought tolerance. Here we assess the genotype-dependent component of the drought-induced transcriptome response in two poplar genotypes differing in drought tolerance. Drought-induced responses were analysed in leaves and root apices and were compared with available transcriptome data from other Populus species.

Results

Using a multi-species designed microarray, a genomic DNA-based selection of probesets provided an unambiguous between-genotype comparison. Analyses of functional group enrichment enabled the extraction of processes physiologically relevant to drought response. The drought-driven changes in gene expression occurring in root apices were consistent across treatments and genotypes. For mature leaves, the transcriptome response varied weakly but in accordance with the duration of water deficit. A differential clustering algorithm revealed similar and divergent gene co-expression patterns among the two genotypes. Since moderate stress levels induced similar physiological responses in both genotypes, the genotype-dependent transcriptional responses could be considered as intrinsic divergences in genome functioning. Our meta-analysis detected several candidate genes and processes that are differentially regulated in root and leaf, potentially under developmental control, and preferentially involved in early and long-term responses to drought.

Conclusions

In poplar, the well-known drought-induced activation of sensing and signalling cascades was specific to the early response in leaves but was found to be general in root apices. Comparing our results to what is known in arabidopsis, we found that transcriptional remodelling included signalling and a response to energy deficit in roots in parallel with transcriptional indices of hampered assimilation in leaves, particularly in the drought-sensitive poplar genotype.

Background

The physiological function of the prion protein remains largely elusive while its key role in prion infection has been expansively documented. To potentially assess this conundrum, we performed a comparative transcriptomic analysis of the brain of wild-type mice with that of transgenic mice invalidated at this locus either at the zygotic or at the adult stages.

Results

Only subtle transcriptomic differences resulting from the Prnp knockout could be evidenced, beside Prnp itself, in the analyzed adult brains following microarray analysis of 24 109 mouse genes and QPCR assessment of some of the putatively marginally modulated loci. When performed at the adult stage, neuronal Prnp disruption appeared to sequentially induce a response to an oxidative stress and a remodeling of the nervous system. However, these events involved only a limited number of genes, expression levels of which were only slightly modified and not always confirmed by RT-qPCR. If not, the qPCR obtained data suggested even less pronounced differences.

Conclusions

These results suggest that the physiological function of PrP is redundant at the adult stage or important for only a small subset of the brain cell population under classical breeding conditions. Following its early reported embryonic developmental regulation, this lack of response could also imply that PrP has a more detrimental role during mouse embryogenesis and that potential transient compensatory mechanisms have to be searched for at the time this locus becomes transcriptionally activated.

Background

The initial stages of development depend on mRNA and proteins accumulated in the oocyte, and during these stages, certain genes are essential for fertilization, first cleavage and embryonic genome activation. The aim of this study was first to search for avian oocyte-specific genes using an in silico and a microarray approaches, then to investigate the temporal and spatial dynamics of the expression of some of these genes during follicular maturation and early embryogenesis.

Results

The in silico approach allowed us to identify 18 chicken homologs of mouse potential oocyte genes found by digital differential display. Using the chicken Affymetrix microarray, we identified 461 genes overexpressed in granulosa cells (GCs) and 250 genes overexpressed in the germinal disc (GD) of the hen oocyte. Six genes were identified using both in silico and microarray approaches. Based on GO annotations, GC and GD genes were differentially involved in biological processes, reflecting different physiological destinations of these two cell layers. Finally we studied the spatial and temporal dynamics of the expression of 21 chicken genes. According to their expression patterns all these genes are involved in different stages of final follicular maturation and/or early embryogenesis in the chicken. Among them, 8 genes (btg4, chkmos, wee, zpA, dazL, cvh, zar1 and ktfn) were preferentially expressed in the maturing occyte and cvh, zar1 and ktfn were also highly expressed in the early embryo.

Conclusion

We showed that in silico and Affymetrix microarray approaches were relevant and complementary in order to find new avian genes potentially involved in oocyte maturation and/or early embryo development, and allowed the discovery of new potential chicken mature oocyte and chicken granulosa cell markers for future studies. Moreover, detailed study of the expression of some of these genes revealed promising candidates for maternal effect genes in the chicken. Finally, the finding concerning the different state of rRNA compared to that of mRNA during the postovulatory period shed light on some mechanisms through which oocyte to embryo transition occurs in the hen.

Background

Since the finishing of the sequencing of the Arabidopsis thaliana genome, the Arabidopsis community and the annotator centers have been working on the improvement of gene annotation at the structural and functional levels. In this context, we have used the large CATMA resource on the Arabidopsis transcriptome to search for genes missed by different annotation processes. Probes on the CATMA microarrays are specific gene sequence tags (GSTs) based on the CDS models predicted by the Eugene software. Among the 24 576 CATMA v2 GSTs, 677 are in regions considered as intergenic by the TAIR annotation. We analyzed the cognate transcriptome data in the CATMA resource and carried out data-mining to characterize novel genes and improve gene models.

Results

The statistical analysis of the results of more than 500 hybridized samples distributed among 12 organs provides an experimental validation for 465 novel genes. The hybridization evidence was confirmed by RT-PCR approaches for 88% of the 465 novel genes. Comparisons with the current annotation show that these novel genes often encode small proteins, with an average size of 137 aa. Our approach has also led to the improvement of pre-existing gene models through both the extension of 16 CDS and the identification of 13 gene models erroneously constituted of two merged CDS.

Conclusion

This work is a noticeable step forward in the improvement of the Arabidopsis genome annotation. We increased the number of Arabidopsis validated genes by 465 novel transcribed genes to which we associated several functional annotations such as expression profiles, sequence conservation in plants, cognate transcripts and protein motifs.

CATdb is a free resource available at http://urgv.evry.inra.fr/CATdb that provides public access to a large collection of transcriptome data for Arabidopsis thaliana produced by a single Complete Arabidopsis Transcriptome Micro Array (CATMA) platform. CATMA probes consist of gene-specific sequence tags (GSTs) of 150–500 bp. The v2 version of CATMA contains 24 576 GST probes representing most of the predicted A. thaliana genes, and 615 probes tiling the chloroplastic and mitochondrial genomes. Data in CATdb are entirely processed with the same standardized protocol, from microarray printing to data analyses. CATdb contains the results of 53 projects including 1724 hybridized samples distributed between 13 different organs, 49 different developmental conditions, 45 mutants and 63 environmental conditions. All the data contained in CATdb can be downloaded from the web site and subsets of data can be sorted out and displayed either by keywords, by experiments, genes or lists of genes up to 100. CATdb gives an easy access to the complete description of experiments with a picture of the experiment design.