Abstract Nitrogen fertilizers are essential for crop productivity but cause environmental harm, necessitating the development of cultivars that thrive under limited nitrogen. This study investigates the transcriptomic response to nitrate in Arabidopsis thaliana (a model dicot), Brachypodium distachyon (a model Pooideae), and Hordeum vulgare (barley, a domesticated Pooideae) to identify conserved and species-specific molecular mechanisms. Using RNA-seq after 1.5 and 3 hours of nitrate treatment, we found that core nitrate-responsive biological processes – such as nitrate transport, assimilation, carbon metabolism, and hormone signaling – are largely conserved across species. However, comparative analysis at gene level based on orthology revealed specificities between the species. For instance, rRNA processing was uniquely stimulated in Arabidopsis, while cysteine biosynthesis from serine and gibberellin biosynthesis were specifically regulated in Brachypodium and barley. Orthologs of key nitrate-responsive genes (e.g., NRT , NLP , TCP20 ) exhibited variable regulation, reflecting potential adaptations linked to domestication or nutrient acquisition strategies. These findings highlight the importance of integrating model and crop species to uncover targets for improving nitrogen use efficiency in cereals. This study provides a general pipeline combining Gene Ontology enrichment and orthology analysis to compare transcriptomic responses across species. This workflow is highly useful for any comparative transcriptomic study in which differences in annotation quality and gene divergence between genomes may bias interpretation.

Abstract Background Plant genomes harbour a substantial proportion of dispensable genes – present only in a subset of individuals – that differ from ubiquitously shared core genes in multiple genomic and expression features. While these differences have been repeatedly documented, the factors shaping gene dispensability remain poorly understood. Results We assembled a pan-gene set from eight maize inbred lines from American and European germplasms, together with their transcriptomic profile across 22 tissues/conditions, revealing the genomic and transcriptomic determinants of maize gene dispensability. Multivariate analysis demonstrates that gene expression level and purifying selection – rather than gene size – are the primary factors distinguishing core from dispensable genes. Dispensable genes overlap Helitrons at 4.6 times the rate of core genes, implicating Helitron -mediated gene capture as a major mechanism of dispensable gene formation. Classifying genes into stably expressed, variably expressed, and on-off categories shows that all three classes contain dispensable genes, though in different proportions than for core genes. Contrary to previous assumptions, we show that dispensable genes can participate in basal biological functions just as core genes, and that gene duplication likely provides only a partial mechanism for functional complementation of accessory genes absence. Conclusions Our results provide novel insights into the molecular and evolutionary factors distinguishing core from dispensable genes and into the biological mechanisms shaping gene dispensability in maize, and demonstrate that classifying genes by transcriptional patterns provides a powerful framework for understanding the biological functions and evolutionary dynamics of both core and dispensable genes.

Monoterpene indole alkaloids (MIAs) are valuable metabolites produced in numerous medicinal plants from the Apocynaceae family such as Alstonia scholaris, which synthesizes strictamine, a MIA displaying neuropharmacological properties of a potential importance. To get insights into the MIA metabolism in A. scholaris, we studied here both the spatial and transcriptional regulations of MIA genes by performing a robust transcriptomics analysis of the main plant organs, leaf epidermis but also by sequencing RNA from leaves transiently overexpressing the master transcriptional regulator MYC2. These transcriptomic studies notably demonstrated that the first steps of the MIA pathway are successively distributed in the internal phloem associated parenchyma and epidermis, and that MYC2 exerts a remarkable transcriptional effect by modulating the expression of around 1000 genes. By combining these distinct datasets, we initiated the search for MIA-related genes encoding CYP71, based on the similarity of expression compared to already known MIA genes. Transient expression of these candidates in Nicotiana benthamiana leaves and yeast notably led to the identification of a related isoform of rhazimal synthase (RHS) capable of converting the MIA precursor geissoschizine into akuammicine, strictamine and 16-epi-pleiocarpamine. Investigating its catalytic mechanism revealed that strictamine results from rhazimal deformylation and that a similar mechanism may also explain 16-epi-pleiocarpamine synthesis. This prompted us to rename these enzymes geissoschizine cyclase due to their capacity of cyclizing geissoschizine into three different MIA scaffolds and to form both C-C and C-N bonds. This identification thus illustrates the potential of integrating spatial and transcriptional regulation analysis for MIA gene identification.

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 .

Plant mitochondria represent the largest group of respiring organelles on the planet. Plant mitochondrial messenger RNAs (mRNAs) lack Shine-Dalgarno-like ribosome-binding sites, so it is unknown how plant mitoribosomes recognize mRNA. We show that “mitochondrial translation factors” mTRAN1 and mTRAN2 are land plant–specific proteins, required for normal mitochondrial respiration chain biogenesis. Our studies suggest that mTRANs are noncanonical pentatricopeptide repeat (PPR)–like RNA binding proteins of the mitoribosomal “small” subunit. We identified conserved Adenosine (A)/Uridine (U)-rich motifs in the 5′ regions of plant mitochondrial mRNAs. mTRAN1 binds this motif, suggesting that it is a mitoribosome homing factor to identify mRNAs. We demonstrate that mTRANs are likely required for translation of all plant mitochondrial mRNAs. Plant mitochondrial translation initiation thus appears to use a protein-mRNA interaction that is divergent from bacteria or mammalian mitochondria.

Background Establishing relationships between transcription factors and target genes is essential for understanding the mechanisms regulating gene expression, which play a fundamental role in plant adaptation to the local environment. Despite the importance of this research area and the tremendous progress in sequencing methods such as ChIP-seq and DAP-seq, we are still far from a complete reconstruction of the cis -regulatory landscape. Only a small number of transcription factors can be assessed from experimental data to identify their cis -regulatory binding sites. This highlights the role that in silico approaches can play to complement experimental data. Results We have developed Plant-PLMview, a web-accessible database for detecting preferentially cis -regulatory sequences in the gene-proximal regions of 20 plant species. Users of Plant-PLMview can (i) access the proximal regions of controlled genes from 20 plant species, (ii) query their own DNA motifs or access 840 cis -regulatory sequences from various plant resources, and (iii) use a tool called PLMdetect to search for preferentially located motifs in the gene-proximal regions of a list of genes. Results are displayed via a web interface with a list of DNA motifs preferentially located in a region near the start or end of genes, the distribution of these motifs and associated annotations. In addition, a graphical map of the preferential locations of the motifs in the 5’ and 3’-proximal regions of genes provides an overview of all motifs and genes examined. Conclusion Plant-PLMview provides the opportunity to study the proximal landscape of gene regulation in 20 plant genomes. The originality of the database lies in its ease of use thanks to the curated data ( cis -regulatory sequences and proximal regions of genes) and the possibility to search through PLMdetect in the 5’-gene-proximal region, but also in the 3’-gene-proximal region, which is rarely explored. The web interface provides numerous graphical views that allow the users to get an overview and interpret the results more easily.

Partial resistance in plants generally exerts a low selective pressure on pathogens, and thus ensuring their durability in agrosystems. However, little is known about the effect of partial resistance on the molecular mechanisms of pathogenicity, a knowledge that could advance plant breeding for sustainable plant health. Here we investigate the gene expression of Phytophthora capsici during infection of pepper (Capsicum annuum L.), where only partial genetic resistance is reported, using Illumina RNA-seq. Comparison of transcriptomes of P. capsici infecting susceptible and partially resistant peppers identified a small number of genes that redirected its own resources into lipid biosynthesis to subsist on partially resistant plants. The adapted and non-adapted isolates of P. capsici differed in expression of genes involved in nucleic acid synthesis and transporters. Transient ectopic expression of the RxLR effector genes CUST_2407 and CUST_16519 in pepper lines differing in resistance levels revealed specific host-isolate interactions that either triggered local necrotic lesions (hypersensitive response or HR) or elicited leave abscission (extreme resistance or ER), preventing the spread of the pathogen to healthy tissue. Although these effectors did not unequivocally explain the quantitative host resistance, our findings highlight the importance of plant genes limiting nutrient resources to select pepper cultivars with sustainable resistance to P. capsici.

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 biosyn-thesis 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 ex-pressed 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.

Identification of cis -regulatory sequences controlling gene expression is an arduous challenge that is being actively explored to discover key genetic factors responsible for traits of agronomic interest. Here, we used a genome-wide de novo approach to investigate preferentially located motifs (PLMs) in the proximal cis -regulatory landscape of Arabidopsis thaliana and Zea mays . We report three groups of PLMs in both the 5’- and 3’-gene-proximal regions and emphasize conserved PLMs in both species, particularly in the 3’-gene-proximal region. Comparison with resources from transcription factor and microRNA binding sites shows that 79% of the identified PLMs are unassigned, although some are supported by MNase-defined cistrome occupancy analysis. Enrichment analyses further reveal that unassigned PLMs provide functional predictions that differ from those derived from transcription factor and microRNA binding sites. Our study provides a comprehensive map of PLMs and demonstrates their potential utility for future characterization of orphan genes in plants.

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.

Shoot branching mechanisms where branches arise in association with leaves – referred to as lateral or axillary branching – evolved by convergence in the sporophyte of vascular plants and the gametophyte of bryophytes, and accompanied independent events of plant architectural diversification 1 . Previously, we showed that three hormonal cues, including auxin, have been recruited independently to co-ordinate branch patterning in flowering plant leafy shoots and moss gametophores (Coudert, Palubicki et al ., 2015) 2–4 . Moreover, auxin-mediated apical dominance, which relies on local auxin production, has been proposed as a unifying molecular regulatory mechanism of branch development across land plants 5 . Whilst our previous work in the moss Physcomitrium patens has gathered indirect evidence supporting the notion that auxin synthesized in gametophore apices regulates branch formation at a distance 2 , direct genetic evidence for a role of auxin biosynthesis in gametophore branching control is still lacking. Here, we show that gametophore apex decapitation promotes branch emergence through massive and rapid transcriptional reprogramming of auxin-responsive genes and altering auxin biosynthesis gene activity. Specifically, we identify a subset of P. patens TRYPTOPHAN AMINO-TRANSFERASE ( TAR ) and YUCCA FLAVIN MONOOXYGENASE-LIKE ( YUC ) auxin biosynthesis genes expressed in apical and basal regions of the gametophore, and show that they are essential for branch initiation and outgrowth control. Our results demonstrate that local auxin biosynthesis coordinates branch patterning in moss and thus constitutes a shared and ancient feature of shoot architecture control in land plants.

Abstract The high mutational load of mitochondrial genomes combined with their uniparental inheritance and high polyploidy favors the maintenance of deleterious mutations within populations. How cells compose and adapt to the accumulation of disadvantageous mitochondrial alleles remains unclear. Most harmful changes are likely corrected by purifying selection, however, the intimate collaboration between mitochondria- and nuclear-encoded gene products offers theoretical potential for compensatory adaptive changes. In plants, cytoplasmic male sterilities are known examples of nucleo-mitochondrial coadaptation situations in which nuclear-encoded restorer of fertility (Rf) genes evolve to counteract the effect of mitochondria-encoded cytoplasmic male sterility (CMS) genes and restore fertility. Most cloned Rfs belong to a small monophyletic group, comprising 26 pentatricopeptide repeat genes in Arabidopsis, called Rf-like (RFL). In this analysis, we explored the functional diversity of RFL genes in Arabidopsis and found that the RFL8 gene is not related to CMS suppression but essential for plant embryo development. In vitro-rescued rfl8 plantlets are deficient in the production of the mitochondrial heme–lyase complex. A complete ensemble of molecular and genetic analyses allowed us to demonstrate that the RFL8 gene has been selected to permit the translation of the mitochondrial ccmFN2 gene encoding a heme–lyase complex subunit which derives from the split of the ccmFN gene, specifically in Brassicaceae plants. This study represents thus a clear case of nuclear compensation to a lineage-specific mitochondrial genomic rearrangement in plants and demonstrates that RFL genes can be selected in response to other mitochondrial deviancies than CMS suppression.

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.

Summary 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.

AbstractPlant 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).

As fixed and long living organisms subjected to repeated environmental stresses, trees have developed mechanisms such as phenotypic plasticity that help them to cope with fluctuating environmental conditions. Here, we tested the role DNA methylation as a hub of integration, linking plasticity and physiological response to water deficit in the shoot apical meristem of the model tree poplar (Populus). Using a reverse genetic approach, we compared hypomethylated RNAi-ddm1 lines to wild-type trees for drought tolerance. An integrative analysis was realized with phytohormone balance, methylomes, transcriptomes and mobilomes. Hypomethylated lines were more tolerant when subjected to moderate water deficit and were intrinsically more tolerant to drought-induced cavitation. The alteration of the DDM1 machinery induced variation in DNA methylation in a cytosine context dependent manner, both in genes and transposable elements. Hypomethylated lines subjected to water deficit showed altered expression of genes involved in phytohormone pathways, such as salicylic acid and modified hormonal balance. Several transposable elements showed stress- and/or line-specific patterns of reactivation, and we could detect copy number variations for two of them in stressed ddm1 lines. Overall, our data highlight two major roles for DNA methylation in the shoot apical meristem: control of stress response and plasticity through transduction of hormone signaling and maintenance of genome integrity through the control of transposable elements.

Seeds are complex biological systems comprising three genetically distinct tissues nested one inside another (embryo, endosperm and maternal tissues). However, the complexity of the kernel makes it difficult to understand inter compartment interactions without access to spatially accurate information. Here we took advantage of the large size of the maize kernel to characterize genome-wide expression profiles of tissues at embryo/endosperm interfaces. Our analysis identifies specific transcriptomic signatures in two interface tissues compared to whole seed compartments: The scutellar aleurone layer (SAL), and the newly named endosperm adjacent to scutellum (EAS). The EAS, which appears around 9 days after pollination and persists for around 11 days, is confined to one to three endosperm cell layers adjacent to the embryonic scutellum. Its transcriptome is enriched in genes encoding transporters. The absence of the embryo in an embryo specific ( emb ) mutant can alter the expression pattern of EAS marker genes. The detection of cell death in some EAS cells together with an accumulation of crushed cell walls suggests that the EAS is a dynamic zone from which cell layers in contact with the embryo are regularly eliminated, and to which additional endosperm cells are recruited as the embryo grows.

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.

Characterizing the molecular processes developed by plants to respond to environmental cues is a major task to better understand local adaptation. DNA methylation is a chromatin mark involved in the transcriptional silencing of transposable elements (TEs) and gene expression regulation. While the molecular bases of DNA methylation regulation are now well described, involvement of DNA methylation in plant response to environmental cues remains poorly characterized. Here, using the TE-rich maize genome and analyzing methylome response to prolonged cold at the chromosome and feature scales, we investigate how genomic architecture affects methylome response to stress in a cold-sensitive genotype. Interestingly, we show that cold stress induces a genome-wide methylation increase through the hypermethylation of TE sequences and centromeres. Our work highlights a cytosine context-specific response of TE methylation that depends on TE types, chromosomal location and proximity to genes. The patterns observed can be explained by the parallel transcriptional activation of multiple DNA methylation pathways that methylate TEs in the various chromatin locations where they reside. Our results open new insights into the possible role of genome-wide DNA methylation in phenotypic response to stress.