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 .

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.

Evolutionary analyses and expression patterns of TCP genes in Ranunculales.

TCP transcription factors play a role in a large number of developmental processes and are at the crossroads of numerous hormonal biosynthetic and signaling pathways. The complete repertoire of TCP genes has already been characterized in several plant species, but not in any species of early diverging eudicots. We focused on the order Ranunculales because of its phylogenetic position as sister group to all other eudicots and its important morphological diversity. Results show that all the TCP genes expressed in the floral transcriptome of Nigella damascena (Ranunculaceae) are the orthologs of the TCP genes previously identified from the fully sequenced genome of Aquilegia coerulea. Phylogenetic analyses combined with the identification of conserved amino acid motifs suggest that six paralogous genes of class I TCP transcription factors were present in the common ancestor of angiosperms. We highlight independent duplications in core eudicots and Ranunculales within the class I and class II subfamilies, resulting in different numbers of paralogs within the main subclasses of TCP genes. This has most probably major consequences on the functional diversification of these genes in different plant clades. The expression patterns of TCP genes in Nigella damascena were consistent with the general suggestion that CIN and class I TCP genes may have redundant roles or take part in same pathways, while CYC/TB1 genes have more specific actions. Our findings open the way for future studies at the tissue level, and for investigating redundancy and subfunctionalisation in TCP genes and their role in the evolution of morphological novelties.

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.

Transcriptional reprogramming during floral fate acquisition.

Coordinating growth and patterning is essential for eukaryote morphogenesis. In plants, auxin is a key regulator of morphogenesis implicated throughout development. Despite this central role, our understanding of how auxin coordinates cell fate and growth changes is still limited. Here, we addressed this question using a combination of genomic screens to delve into the transcriptional network induced by auxin at the earliest stage of flower development, prior to morphological changes. We identify a shoot-specific network suggesting that auxin initiates growth through an antagonistic regulation of growth-promoting and growth-repressive hormones, quasi-synchronously to floral fate specification. We further identify two DNA-binding One Zinc Finger (DOF) transcription factors acting in an auxin-dependent network that could interface growth and cell fate from the early stages of flower development onward.

Apical dominance control by TAR-YUC-mediated auxin biosynthesis is a deep homology of land plants.

A key aim in biology is to identify which genetic changes contributed to the evolution of form through time. Apical dominance, the inhibitory effect exerted by shoot apices on the initiation or outgrowth of distant lateral buds, is a major regulatory mechanism of plant form.1 Nearly a century of studies in the sporophyte of flowering plants have established the phytohormone auxin as a front-runner in the search for key factors controlling apical dominance,2,3 identifying critical roles for long-range polar auxin transport and local auxin biosynthesis in modulating shoot branching.4-10 A capacity for lateral branching evolved by convergence in the gametophytic shoot of mosses and primed its diversification;11 however, polar auxin transport is relatively unimportant in this developmental process,12 the contribution of auxin biosynthesis genes has not been assessed, and more generally, the extent of conservation in apical dominance regulation within the land plants remains largely unknown. To fill this knowledge gap, we sought to identify genetic determinants of apical dominance in the moss Physcomitrium patens. Here, we show that leafy shoot apex decapitation releases apical dominance through massive and rapid transcriptional reprogramming of auxin-responsive genes and altering auxin biosynthesis gene activity. We pinpoint a subset of P. patens TRYPTOPHAN AMINO-TRANSFERASE (TAR) and YUCCA FLAVIN MONOOXYGENASE-LIKE (YUC) auxin biosynthesis genes expressed in the main and lateral shoot apices and show that they are essential for coordinating branch initiation and outgrowth. Our results demonstrate that local auxin biosynthesis acts as a pivotal regulator of apical dominance in moss and constitutes a shared mechanism underpinning shoot architecture control in land plants.

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.

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.

An Arabidopsis mutant deficient in phosphatidylinositol-4-phosphate kinases ß1 and ß2 displays altered auxin-related responses in roots.

Phosphatidylinositol 4-kinases (PI4Ks) are the first enzymes that commit phosphatidylinositol into the phosphoinositide pathway. Here, we show that Arabidopsis thaliana seedlings deficient in PI4Kß1 and ß2 have several developmental defects including shorter roots and unfinished cytokinesis. The pi4kß1ß2 double mutant was insensitive to exogenous auxin concerning inhibition of root length and cell elongation; it also responded more slowly to gravistimulation. The pi4kß1ß2 root transcriptome displayed some similarities to a wild type plant response to auxin. Yet, not all the genes displayed such a constitutive auxin-like response. Besides, most assessed genes did not respond to exogenous auxin. This is consistent with data with the transcriptional reporter DR5-GUS. The content of bioactive auxin in the pi4kß1ß2 roots was similar to that in wild-type ones. Yet, an enhanced auxin-conjugating activity was detected and the auxin level reporter DII-VENUS did not respond to exogenous auxin in pi4kß1ß2 mutant. The mutant exhibited altered subcellular trafficking behavior including the trapping of PIN-FORMED 2 protein in rapidly moving vesicles. Bigger and less fragmented vacuoles were observed in pi4kß1ß2 roots when compared to the wild type. Furthermore, the actin filament web of the pi4kß1ß2 double mutant was less dense than in wild-type seedling roots, and less prone to rebuilding after treatment with latrunculin B. A mechanistic model is proposed in which an altered PI4K activity leads to actin filament disorganization, changes in vesicle trafficking, and altered auxin homeostasis and response resulting in a pleiotropic root phenotypes.

Transcriptome Analysis Reveals Putative Target Genes of APETALA3-3 During Early Floral Development in Nigella damascena L.

Even though petals are homoplastic structures, their identity consistently involves genes of the APETALA3 (AP3) lineage. However, the extent to which the networks downstream of AP3 are conserved in species with petals of different evolutionary origins is unknown. In Ranunculaceae, the specificity of the AP3-III lineage offers a great opportunity to identify the petal gene regulatory network in a comparative framework. Using a transcriptomic approach, we investigated putative target genes of the AP3-III ortholog NdAP3-3 in Nigella damascena at early developmental stages when petal identity is determined, and we compared our data with that from selected eudicot species. We generated a de novo reference transcriptome to carry out a differential gene expression analysis between the wild-type and mutant NdAP3-3 genotypes differing by the presence vs. absence of petals at early stages of floral development. Among the 1,620 genes that were significantly differentially expressed between the two genotypes, functional annotation suggested a large involvement of nuclear activities, including regulation of transcription, and enrichment in processes linked to cell proliferation. Comparing with Arabidopsis data, we found that highly conserved genes between the two species are enriched in homologs of direct targets of the AtAP3 protein. Integrating AP3-3 binding site data from another Ranunculaceae species, Aquilegia coerulea, allowed us to identify a set of 18 putative target genes that were conserved between the three species. Our results suggest that, despite the independent evolutionary origin of petals in core eudicots and Ranunculaceae, a small conserved set of genes determines petal identity and early development in these taxa.

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.

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.

Plant nitrogen supply affects the Botrytis cinerea infection process and modulates known and novel virulence factors.

Plant nitrogen (N) fertilization is known to affect disease; however, the underlying mechanisms remain mostly unknown. We investigated the impact of N supply on the Arabidopsis thaliana-Botrytis cinerea interaction. A. thaliana plants grown in low nitrate were more tolerant to all wild-type B. cinerea strains tested. We determined leaf nitrate concentrations and showed that they had a limited impact on B. cinerea growth in vitro. For the first time, we performed a dual RNA-Seq of infected leaves of plants grown with different nitrate concentrations. Transcriptome analysis showed that plant and fungal transcriptomes were marginally affected by plant nitrate supply. Indeed, only a limited set of plant (182) and fungal (22) genes displayed expression profiles altered by nitrate supply. The expression of selected genes was confirmed by quantitative reverse transcription PCR at 6 hr postinfection (hpi) and analysed at a later time point (24 hpi). We selected three of the 22 B. cinerea genes identified for further analysis. B. cinerea mutants affected in these genes were less aggressive than the wild-type strain. We also showed that plants grown in ammonium were more tolerant to B. cinerea. Furthermore, expression of the selected B. cinerea genes in planta was altered when plants were grown with ammonium instead of nitrate, demonstrating an impact of the nature of N supplied to plants on the interaction. Identification of B. cinerea genes expressed differentially in planta according to plant N supply unveils two novel virulence functions required for full virulence in A. thaliana: a secondary metabolite (SM) and an acidic protease (AP).

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.

Mitochondrial Transcriptome Control and Intercompartment Cross-Talk During Plant Development.

We address here organellar genetic regulation and intercompartment genome coordination. We developed earlier a strategy relying on a tRNA-like shuttle to mediate import of nuclear transgene-encoded custom RNAs into mitochondria in plants. In the present work, we used this strategy to drive trans-cleaving hammerhead ribozymes into the organelles, to knock down specific mitochondrial RNAs and analyze the regulatory impact. In a similar approach, the tRNA mimic was used to import into mitochondria in Arabidopsis thaliana the orf77, an RNA associated with cytoplasmic male sterility in maize and possessing sequence identities with the atp9 mitochondrial RNA. In both cases, inducible expression of the transgenes allowed to characterise early regulation and signaling responses triggered by these respective manipulations of the organellar transcriptome. The results imply that the mitochondrial transcriptome is tightly controlled by a "buffering" mechanism at the early and intermediate stages of plant development, a control that is released at later stages. On the other hand, high throughput analyses showed that knocking down a specific mitochondrial mRNA triggered a retrograde signaling and an anterograde nuclear transcriptome response involving a series of transcription factor genes and small RNAs. Our results strongly support transcriptome coordination mechanisms within the organelles and between the organelles and the nucleus.

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.

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.

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.

A systems biology approach uncovers a gene co-expression network associated with cell wall degradability in maize.

Understanding the mechanisms triggering variation of cell wall degradability is a prerequisite to improving the energy value of lignocellulosic biomass for animal feed or biorefinery. Here, we implemented a multiscale systems approach to shed light on the genetic basis of cell wall degradability in maize. We demonstrated that allele replacement in two pairs of near-isogenic lines at a region encompassing a major quantitative trait locus (QTL) for cell wall degradability led to phenotypic variation of a similar magnitude and sign to that expected from a QTL analysis of cell wall degradability in the F271 × F288 recombinant inbred line progeny. Using DNA sequences within the QTL interval of both F271 and F288 inbred lines and Illumina RNA sequencing datasets from internodes of the selected near-isogenic lines, we annotated the genes present in the QTL interval and provided evidence that allelic variation at the introgressed QTL region gives rise to coordinated changes in gene expression. The identification of a gene co-expression network associated with cell wall-related trait variation revealed that the favorable F288 alleles exploit biological processes related to oxidation-reduction, regulation of hydrogen peroxide metabolism, protein folding and hormone responses. Nested in modules of co-expressed genes, potential new cell-wall regulators were identified, including two transcription factors of the group VII ethylene response factor family, that could be exploited to fine-tune cell wall degradability. Overall, these findings provide new insights into the regulatory mechanisms by which a major locus influences cell wall degradability, paving the way for its map-based cloning in maize.

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.