A cell wall-associated gene network shapes leaf boundary domains.

Boundary domains delimit and organize organ growth throughout plant development almost relentlessly, building plant architecture and morphogenesis. Boundary domains display reduced growth and orchestrate development of adjacent tissues in a non-cell-autonomous manner. How these two functions are achieved remains elusive despite the identification of several boundary-specific genes. Here, we show using morphometrics at the organ and cellular levels that leaf boundary domain development requires SPINDLY (SPY), an O-fucosyltransferase, to act as cell growth repressor. Furthermore, we show that SPY acts redundantly with the CUP-SHAPED COTYLEDON transcription factors (CUC2 and CUC3), which are major determinants of boundaries development. Accordingly, at the molecular level CUC2 and SPY repress a common set of genes involved in cell wall loosening, providing a molecular framework for the growth repression associated with boundary domains. Atomic force microscopy confirmed that young leaf boundary domain cells have stiffer cell walls than marginal outgrowth. This differential cell wall stiffness was reduced in spy mutant plants. Taken together, our data reveal a concealed CUC2 cell wall-associated gene network linking tissue patterning with cell growth and mechanics.

Topoisomerase VI participates in an insulator-like function that prevents H3K9me2 spreading.

The organization of the genome into transcriptionally active and inactive chromatin domains requires well-delineated chromatin boundaries and insulator functions in order to maintain the identity of adjacent genomic loci with antagonistic chromatin marks and functionality. In plants that lack known chromatin insulators, the mechanisms that prevent heterochromatin spreading into euchromatin remain to be identified. Here, we show that DNA Topoisomerase VI participates in a chromatin boundary function that safeguards the expression of genes in euchromatin islands within silenced heterochromatin regions. While some transposable elements are reactivated in mutants of the Topoisomerase VI complex, genes insulated in euchromatin islands within heterochromatic regions of the Arabidopsis thaliana genome are specifically down-regulated. H3K9me2 levels consistently increase at euchromatin island loci and decrease at some transposable element loci. We further show that Topoisomerase VI physically interacts with S-adenosylmethionine synthase methionine adenosyl transferase 3 (MAT3), which is required for H3K9me2. A Topoisomerase VI defect affects MAT3 occupancy on heterochromatic elements and its exclusion from euchromatic islands, thereby providing a possible mechanistic explanation to the essential role of Topoisomerase VI in the delimitation of chromatin domains.

DIACYLGLYCEROL KINASE 5 participates in flagellin-induced signaling in Arabidopsis.

Flagellin perception is a keystone of pattern-triggered immunity in plants. The recognition of this protein by a plasma membrane (PM) receptor complex is the beginning of a signaling cascade that includes protein phosphorylation and the production of reactive oxygen species (ROS). In both Arabidopsis (Arabidopsis thaliana) seedlings and suspension cells, we found that treatment with flg22, a peptide corresponding to the most conserved domain of bacterial flagellin, caused a rapid and transient decrease in the level of phosphatidylinositol (PI) 4,5-bisphosphate along with a parallel increase in phosphatidic acid (PA). In suspension cells, inhibitors of either phosphoinositide-dependent phospholipases C (PLC) or diacylglycerol kinases (DGKs) inhibited flg22-triggered PA production and the oxidative burst. In response to flg22, receptor-like kinase-deficient fls2, bak1, and bik1 mutants (FLAGELLIN SENSITIVE 2, BRASSINOSTEROID INSENSITIVE 1-associated kinase 1, and BOTRYTIS-INDUCED KINASE 1, respectively) produced less PA than wild-type (WT) plants, whereas this response did not differ in NADPH oxidase-deficient rbohD (RESPIRATORY BURST OXIDASE HOMOLOG D) plants. Among the DGK-deficient lines tested, the dgk5.1 mutant produced less PA and less ROS after flg22 treatment compared with WT seedlings. In response to flg22, dgk5.1 plants showed lower callose accumulation and impaired resistance to Pseudomonas syringae pv. tomato DC3000 hrcC-. Transcriptomics revealed that the basal expression of defense-related genes was altered in dgk5.1 seedlings compared with the WT. A GFP-DGK5 fusion protein localized to the PM, where RBOHD and PLC2 (proteins involved in plant immunity) are also located. The role of DGK5 and its enzymatic activity in flagellin signaling and fine-tuning of early immune responses in plant-microbe interactions is discussed.

Wide cross-species RNA-Seq comparison reveals convergent molecular mechanisms involved in nickel hyperaccumulation across dicotyledons.

The Anthropocene epoch is associated with the spreading of metals in the environment increasing oxidative and genotoxic stress on organisms. Interestingly, c. 520 plant species growing on metalliferous soils acquired the capacity to accumulate and tolerate a tremendous amount of nickel in their shoots. The wide phylogenetic distribution of these species suggests that nickel hyperaccumulation evolved multiple times independently. However, the exact nature of these mechanisms and whether they have been recruited convergently in distant species is not known. To address these questions, we have developed a cross-species RNA-Seq approach combining differential gene expression analysis and cluster of orthologous group annotation to identify genes linked to nickel hyperaccumulation in distant plant families. Our analysis reveals candidate orthologous genes encoding convergent function involved in nickel hyperaccumulation, including the biosynthesis of specialized metabolites and cell wall organization. Our data also point out that the high expression of IREG/Ferroportin transporters recurrently emerged as a mechanism involved in nickel hyperaccumulation in plants. We further provide genetic evidence in the hyperaccumulator Noccaea caerulescens for the role of the NcIREG2 transporter in nickel sequestration in vacuoles. Our results provide molecular tools to better understand the mechanisms of nickel hyperaccumulation and study their evolution in plants.

RNAi suppression of DNA methylation affects the drought stress response and genome integrity in transgenic poplar.

Trees are long-lived organisms that continuously adapt to their environments, a process in which epigenetic mechanisms are likely to play a key role. Via downregulation of the chromatin remodeler DECREASED IN DNA METHYLATION 1 (DDM1) in poplar (Populus tremula × Populus alba) RNAi lines, we examined how DNA methylation coordinates genomic and physiological responses to moderate water deficit. We compared the growth and drought response of two RNAi-ddm1 lines to wild-type (WT) trees under well-watered and water deficit/rewatering conditions, and analyzed their methylomes, transcriptomes, mobilomes and phytohormone contents in the shoot apical meristem. The RNAi-ddm1 lines were more tolerant to drought-induced cavitation but did not differ in height or stem diameter growth. About 5000 differentially methylated regions were consistently detected in both RNAi-ddm1 lines, colocalizing with 910 genes and 89 active transposable elements. Under water deficit conditions, 136 differentially expressed genes were found, including many involved in phytohormone pathways; changes in phytohormone concentrations were also detected. Finally, the combination of hypomethylation and drought led to the mobility of two transposable elements. Our findings suggest major roles for DNA methylation in regulation of genes involved in hormone-related stress responses, and the maintenance of genome integrity through repression of transposable elements.

Gene expression predictions and networks in natural populations supports the omnigenic theory.

BACKGROUND: Recent literature on the differential role of genes within networks distinguishes core from peripheral genes. If previous works have shown contrasting features between them, whether such categorization matters for phenotype prediction remains to be studied. RESULTS: We measured 17 phenotypic traits for 241 cloned genotypes from a Populus nigra collection, covering growth, phenology, chemical and physical properties. We also sequenced RNA for each genotype and built co-expression networks to define core and peripheral genes. We found that cores were more differentiated between populations than peripherals while being less variable, suggesting that they have been constrained through potentially divergent selection. We also showed that while cores were overrepresented in a subset of genes statistically selected for their capacity to predict the phenotypes (by Boruta algorithm), they did not systematically predict better than peripherals or even random genes. CONCLUSION: Our work is the first attempt to assess the importance of co-expression network connectivity in phenotype prediction. While highly connected core genes appear to be important, they do not bear enough information to systematically predict better quantitative traits than other gene sets.

Synthetic data sets for the identification of key ingredients for RNA-seq differential analysis.

Numerous statistical pipelines are now available for the differential analysis of gene expression measured with RNA-sequencing technology. Most of them are based on similar statistical frameworks after normalization, differing primarily in the choice of data distribution, mean and variance estimation strategy and data filtering. We propose an evaluation of the impact of these choices when few biological replicates are available through the use of synthetic data sets. This framework is based on real data sets and allows the exploration of various scenarios differing in the proportion of non-differentially expressed genes. Hence, it provides an evaluation of the key ingredients of the differential analysis, free of the biases associated with the simulation of data using parametric models. Our results show the relevance of a proper modeling of the mean by using linear or generalized linear modeling. Once the mean is properly modeled, the impact of the other parameters on the performance of the test is much less important. Finally, we propose to use the simple visualization of the raw P-value histogram as a practical evaluation criterion of the performance of differential analysis methods on real data sets.

Multi-omics Analysis Reveals Sequential Roles for ABA during Seed Maturation.

Abscisic acid (ABA) is an important hormone for seed development and germination whose physiological action is modulated by its endogenous levels. Cleavage of carotenoid precursors by 9-cis epoxycarotenoid dioxygenase (NCED) and inactivation of ABA by ABA 8'-hydroxylase (CYP707A) are key regulatory metabolic steps. In Arabidopsis (Arabidopsis thaliana), both enzymes are encoded by multigene families, having distinctive expression patterns. To evaluate the genome-wide impact of ABA deficiency in developing seeds at the maturation stage when dormancy is induced, we used a nced2569 quadruple mutant in which ABA deficiency is mostly restricted to seeds, thus limiting the impact of maternal defects on seed physiology. ABA content was very low in nced2569 seeds, similar to the severe mutant aba2; unexpectedly, ABA Glc ester was detected in aba2 seeds, suggesting the existence of an alternative metabolic route. Hormone content in nced2569 seeds compared with nced259 and wild type strongly suggested that specific expression of NCED6 in the endosperm is mainly responsible for ABA production. In accordance, transcriptome analyses revealed broad similarities in gene expression between nced2569 and either wild-type or nced259 developing seeds. Gene ontology enrichments revealed a large spectrum of ABA activation targets involved in reserve storage and desiccation tolerance, and repression of photosynthesis and cell cycle. Proteome and metabolome profiles in dry nced2569 seeds, compared with wild-type and cyp707a1a2 seeds, also highlighted an inhibitory role of ABA on remobilization of reserves, reactive oxygen species production, and protein oxidation. Down-regulation of these oxidative processes by ABA may have an essential role in dormancy control.

Extensive nuclear reprogramming and endoreduplication in mature leaf during floral induction.

BACKGROUND: The floral transition is a complex developmental event, fine-tuned by various environmental and endogenous cues to ensure the success of offspring production. Leaves are key organs in sensing floral inductive signals, such as a change in light regime, and in the production of the mobile florigen. CONSTANS and FLOWERING LOCUS T are major players in leaves in response to photoperiod. Morphological and molecular events during the floral transition have been intensively studied in the shoot apical meristem. To better understand the concomitant processes in leaves, which are less described, we investigated the nuclear changes in fully developed leaves during the time course of the floral transition. RESULTS: We highlighted new putative regulatory candidates of flowering in leaves. We observed differential expression profiles of genes related to cellular, hormonal and metabolic actions, but also of genes encoding long non-coding RNAs and new natural antisense transcripts. In addition, we detected a significant increase in ploidy level during the floral transition, indicating endoreduplication. CONCLUSIONS: Our data indicate that differentiated mature leaves, possess physiological plasticity and undergo extensive nuclear reprogramming during the floral transition. The dynamic events point at functionally related networks of transcription factors and novel regulatory motifs, but also complex hormonal and metabolic changes.

On a Cold Night: Transcriptomics of Grapevine Flower Unveils Signal Transduction and Impacted Metabolism.

Low temperature is a critical environmental factor limiting plant productivity, especially in northern vineyards. To clarify the impact of this stress on grapevine flower, we used the Vitis array based on Roche-NimbleGen technology to investigate the gene expression of flowers submitted to a cold night. Our objectives were to identify modifications in the transcript levels after stress and during recovery. Consequently, our results confirmed some mechanisms known in grapes or other plants in response to cold stress, notably, (1) the pivotal role of calcium/calmodulin-mediated signaling; (2) the over-expression of sugar transporters and some genes involved in plant defense (especially in carbon metabolism), and (3) the down-regulation of genes encoding galactinol synthase (GOLS), pectate lyases, or polygalacturonases. We also identified some mechanisms not yet known to be involved in the response to cold stress, i.e., (1) the up-regulation of genes encoding G-type lectin S-receptor-like serine threonine-protein kinase, pathogen recognition receptor (PRR5), or heat-shock factors among others; (2) the down-regulation of Myeloblastosis (MYB)-related transcription factors and the Constans-like zinc finger family; and (3) the down-regulation of some genes encoding Pathogen-Related (PR)-proteins. Taken together, our results revealed interesting features and potentially valuable traits associated with stress responses in the grapevine flower. From a long-term perspective, our study provides useful starting points for future investigation.

Accuracy of RNAseq based SNP discovery and genotyping in Populusnigra.

BACKGROUD: Populus nigra is a major tree species of ecological and economic importance for which several initiatives have been set up to create genomic resources. In order to access the large number of Single Nucleotide Polymorphisms (SNPs) typically needed to carry out a genome scan, the present study aimed at evaluating RNA sequencing as a tool to discover and type SNPs in genes within natural populations of P. nigra. RESULTS: We have devised a bioinformatics pipeline to call and type SNPs from RNAseq reads and applied it to P. nigra transcriptomic data. The accuracy of the resulting RNAseq-based SNP calling and typing has been evaluated by (i) comparing their position and alleles to those previously reported in candidate genes, (ii) assessing their genotyping accuracy with respect to a previously available SNP chip and (iii) evaluating their inter-annual repeatability. We found that a combination of several callers yields a good compromise between the number of variants type and the accuracy of genotyping. We further used the resulting genotypic data to carry out basic genetic analyses whose results confirm the quality of the RNAseq-based SNP dataset. CONCLUSIONS: We demonstrated the potential and accuracy of RNAseq as an efficient way to genotype SNPs in P. nigra. These results open prospects towards the use of this technology for quantitative and population genomics studies.

Transcriptome reprogramming due to the introduction of a barley telosome into bread wheat affects more barley genes than wheat.

Despite a long history, the production of useful alien introgression lines in wheat remains difficult mainly due to linkage drag and incomplete genetic compensation. In addition, little is known about the molecular mechanisms underlying the impact of foreign chromatin on plant phenotype. Here, a comparison of the transcriptomes of barley, wheat and a wheat-barley 7HL addition line allowed the transcriptional impact both on 7HL genes of a non-native genetic background and on the wheat gene complement as a result of the presence of 7HL to be assessed. Some 42% (389/923) of the 7HL genes assayed were differentially transcribed, which was the case for only 3% (960/35 301) of the wheat gene complement. The absence of any transcript in the addition line of a suite of chromosome 7A genes implied the presence of a 36 Mbp deletion at the distal end of the 7AL arm; this deletion was found to be in common across the full set of Chinese Spring/Betzes barley addition lines. The remaining differentially transcribed wheat genes were distributed across the whole genome. The up-regulated barley genes were mostly located in the proximal part of the 7HL arm, while the down-regulated ones were concentrated in the distal part; as a result, genes encoding basal cellular functions tended to be transcribed, while those encoding specific functions were suppressed. An insight has been gained into gene transcription in an alien introgression line, thereby providing a basis for understanding the interactions between wheat and exotic genes in introgression materials.

MYB118 represses endosperm maturation in seeds of Arabidopsis.

In the exalbuminous species Arabidopsis thaliana, seed maturation is accompanied by the deposition of oil and storage proteins and the reduction of the endosperm to one cell layer. Here, we consider reserve partitioning between embryo and endosperm compartments. The pattern of deposition, final amount, and composition of these reserves differ between the two compartments, with the embryo representing the principal storage tissue in mature seeds. Complex regulatory mechanisms are known to prevent activation of maturation-related programs during embryo morphogenesis and, later, during vegetative growth. Here, we describe a regulator that represses the expression of maturation-related genes during maturation within the endosperm. MYB118 is transcriptionally induced in the maturing endosperm, and seeds of myb118 mutants exhibit an endosperm-specific derepression of maturation-related genes associated with a partial relocation of storage compounds from the embryo to the endosperm. Moreover, MYB118 activates endosperm-induced genes through the recognition of TAACGG elements. These results demonstrate that the differential partitioning of reserves between the embryo and endosperm in exalbuminous Arabidopsis seeds does not only result from developmental programs that establish the embryo as the preponderant tissue within seeds. This differential partitioning is also regulated by MYB118, which regulates the biosynthesis of reserves at the spatial level during maturation.

Germination Potential of Dormant and Nondormant Arabidopsis Seeds Is Driven by Distinct Recruitment of Messenger RNAs to Polysomes.

Dormancy is a complex evolutionary trait that temporally prevents seed germination, thus allowing seedling growth at a favorable season. High-throughput analyses of transcriptomes have led to significant progress in understanding the molecular regulation of this process, but the role of posttranscriptional mechanisms has received little attention. In this work, we have studied the dynamics of messenger RNA association with polysomes and compared the transcriptome with the translatome in dormant and nondormant seeds of Arabidopsis (Arabidopsis thaliana) during their imbibition at 25 °C in darkness, a temperature preventing germination of dormant seeds only. DNA microarray analysis revealed that 4,670 and 7,028 transcripts were differentially abundant in dormant and nondormant seeds in the transcriptome and the translatome, respectively. We show that there is no correlation between transcriptome and translatome and that germination regulation is also largely translational, implying a selective and dynamic recruitment of messenger RNAs to polysomes in both dormant and nondormant seeds. The study of 5' untranslated region features revealed that GC content and the number of upstream open reading frames could play a role in selective translation occurring during germination. Gene Ontology clustering showed that the functions of polysome-associated transcripts differed between dormant and nondormant seeds and revealed actors in seed dormancy and germination. In conclusion, our results demonstrate the essential role of selective polysome loading in this biological process.

The Polyadenylation Factor Subunit CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR30: A Key Factor of Programmed Cell Death and a Regulator of Immunity in Arabidopsis.

Programmed cell death (PCD) is essential for several aspects of plant life, including development and stress responses. Indeed, incompatible plant-pathogen interactions are well known to induce the hypersensitive response, a localized cell death. Mutational analyses have identified several key PCD components, and we recently identified the mips1 mutant of Arabidopsis (Arabidopsis thaliana), which is deficient for the key enzyme catalyzing the limiting step of myoinositol synthesis. One of the most striking features of mips1 is the light-dependent formation of lesions on leaves due to salicylic acid (SA)-dependent PCD, revealing roles for myoinositol or inositol derivatives in the regulation of PCD. Here, we identified a regulator of plant PCD by screening for mutants that display transcriptomic profiles opposing that of the mips1 mutant. Our screen identified the oxt6 mutant, which has been described previously as being tolerant to oxidative stress. In the oxt6 mutant, a transfer DNA is inserted in the CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR30 (CPSF30) gene, which encodes a polyadenylation factor subunit homolog. We show that CPSF30 is required for lesion formation in mips1 via SA-dependent signaling, that the prodeath function of CPSF30 is not mediated by changes in the glutathione status, and that CPSF30 activity is required for Pseudomonas syringae resistance. We also show that the oxt6 mutation suppresses cell death in other lesion-mimic mutants, including lesion-simulating disease1, mitogen-activated protein kinase4, constitutive expressor of pathogenesis-related genes5, and catalase2, suggesting that CPSF30 and, thus, the control of messenger RNA 3' end processing, through the regulation of SA production, is a key component of plant immune responses.

Scavenging iron: a novel mechanism of plant immunity activation by microbial siderophores.

Siderophores are specific ferric iron chelators synthesized by virtually all microorganisms in response to iron deficiency. We have previously shown that they promote infection by the phytopathogenic enterobacteria Dickeya dadantii and Erwinia amylovora. Siderophores also have the ability to activate plant immunity. We have used complete Arabidopsis transcriptome microarrays to investigate the global transcriptional modifications in roots and leaves of Arabidopsis (Arabidopsis thaliana) plants after leaf treatment with the siderophore deferrioxamine (DFO). Physiological relevance of these transcriptional modifications was validated experimentally. Immunity and heavy-metal homeostasis were the major processes affected by DFO. These two physiological responses could be activated by a synthetic iron chelator ethylenediamine-di(o-hydroxyphenylacetic) acid, indicating that siderophores eliciting activities rely on their strong iron-chelating capacity. DFO was able to protect Arabidopsis against the pathogenic bacterium Pseudomonas syringae pv tomato DC3000. Siderophore treatment caused local modifications of iron distribution in leaf cells visible by ferrocyanide and diaminobenzidine-H2O2 staining. Metal quantifications showed that DFO causes a transient iron and zinc uptake at the root level, which is presumably mediated by the metal transporter iron regulated transporter1 (IRT1). Defense gene expression and callose deposition in response to DFO were compromised in an irt1 mutant. Consistently, plant susceptibility to D. dadantii was increased in the irt1 mutant. Our work shows that iron scavenging is a unique mechanism of immunity activation in plants. It highlights the strong relationship between heavy-metal homeostasis and immunity.

Carotenoid oxidation products are stress signals that mediate gene responses to singlet oxygen in plants.

(1)O(2) (singlet oxygen) is a reactive O(2) species produced from triplet excited chlorophylls in the chloroplasts, especially when plants are exposed to excess light energy. Similarly to other active O(2) species, (1)O(2) has a dual effect: It is toxic, causing oxidation of biomolecules, and it can act as a signal molecule that leads to cell death or to acclimation. Carotenoids are considered to be the main (1)O(2) quenchers in chloroplasts, and we show here that light stress induces the oxidation of the carotenoid ß-carotene in Arabidopsis plants, leading to the accumulation of different volatile derivatives. One such compound, ß-cyclocitral, was found to induce changes in the expression of a large set of genes that have been identified as (1)O(2) responsive genes. In contrast, ß-cyclocitral had little effect on the expression of H(2)O(2) gene markers. ß-Cyclocitral-induced reprogramming of gene expression was associated with an increased tolerance to photooxidative stress. The results indicate that ß-cyclocitral is a stress signal produced in high light that is able to induce defense mechanisms and represents a likely messenger involved in the (1)O(2) signaling pathway in plants.

γ-Aminobutyric acid transaminase deficiency impairs central carbon metabolism and leads to cell wall defects during salt stress in Arabidopsis roots.

Environmental constraints challenge cell homeostasis and thus require a tight regulation of metabolic activity. We have previously reported that the γ-aminobutyric acid (GABA) metabolism is crucial for Arabidopsis salt tolerance as revealed by the NaCl hypersensitivity of the GABA transaminase (GABA-T, At3g22200) gaba-t/pop2-1 mutant. In this study, we demonstrate that GABA-T deficiency during salt stress causes root and hypocotyl developmental defects and alterations of cell wall composition. A comparative genome-wide transcriptional analysis revealed that expression levels of genes involved in carbon metabolism, particularly sucrose and starch catabolism, were found to increase upon the loss of GABA-T function under salt stress conditions. Consistent with the altered mutant cell wall composition, a number of cell wall-related genes were also found differentially expressed. A targeted quantitative analysis of primary metabolites revealed that glutamate (GABA precursor) accumulated while succinate (the final product of GABA metabolism) significantly decreased in mutant roots after 1 d of NaCl treatment. Furthermore, sugar concentration was twofold reduced in gaba-t/pop2-1 mutant roots compared with wild type. Together, our results provide strong evidence that GABA metabolism is a major route for succinate production in roots and identify GABA as a major player of central carbon adjustment during salt stress.

Cell specialization and coordination in Arabidopsis leaves upon pathogenic attack revealed by scRNA-seq.

Plant defense responses involve several biological processes that allow plants to fight against pathogenic attacks. How these different processes are orchestrated within organs and depend on specific cell types is poorly known. Here, using single-cell RNA sequencing (scRNA-seq) technology on three independent biological replicates, we identified several cell populations representing the core transcriptional responses of wild-type Arabidopsis leaves inoculated with the bacterial pathogen Pseudomonas syringae DC3000. Among these populations, we retrieved major cell types of the leaves (mesophyll, guard, epidermal, companion, and vascular S cells) with which we could associate characteristic transcriptional reprogramming and regulators, thereby specifying different cell-type responses to the pathogen. Further analyses of transcriptional dynamics, on the basis of inference of cell trajectories, indicated that the different cell types, in addition to their characteristic defense responses, can also share similar modules of gene reprogramming, uncovering a ubiquitous antagonism between immune and susceptible processes. Moreover, it appears that the defense responses of vascular S cells, epidermal cells, and mesophyll cells can evolve along two separate paths, one converging toward an identical cell fate, characterized mostly by lignification and detoxification functions. As this divergence does not correspond to the differentiation between immune and susceptible cells, we speculate that this might reflect the discrimination between cell-autonomous and non-cell-autonomous responses. Altogether our data provide an upgraded framework to describe, explore, and explain the specialization and the coordination of plant cell responses upon pathogenic challenge.

RNAseq based variant dataset in a black poplar association panel.

OBJECTIVE: Black poplar (Populus nigra L.) is a species native to Eurasia with a wide distribution area. It is an ecologically important species from riparian ecosystems, that is used as a parent of interspecific (P. deltoides x P. nigra) cultivated poplar hybrids. Variant detection from transcriptomics sequences of 241 P. nigra individuals, sampled in natural populations from 11 river catchments (in four European countries) is described here. These data provide new valuable resources for population structure analysis, population genomics and genome-wide association studies. DATA DESCRIPTION: We generated transcriptomics data from a mixture of young differentiating xylem and cambium tissues of 480 Populus nigra trees sampled in a common garden experiment located at Orléans (France), corresponding to 241 genotypes (2 clonal replicates per genotype, at maximum) by using RNAseq technology. We launched on the resulting sequences an in-silico pipeline that allowed us to obtain 878,957 biallelic polymorphisms without missing data. More than 99% of these positions are annotated and 98.8% are located on the 19 chromosomes of the P. trichocarpa reference genome. The raw RNAseq sequences are available at the NCBI Sequence Read Archive SPR188754 and the variant dataset at the Recherche Data Gouv repository under https://doi.org/10.15454/8DQXK5 .