Search

Chloroplasts are the powerhouse of the plant cell, yet they are resource-intensive and will cause photooxidative damage if their activity overshoots the demands of growth. The adjustment of chloroplast activity to match growth is therefore vital for stress acclimation. Here we identify a novel post-translational mechanism linking the conserved eukaryotic TOR kinase that promotes growth and the guanosine tetraphosphate (ppGpp) signaling pathway of prokaryotic origin that regulates chloroplast activity, and photosynthesis in particular. We show that RelA SpoT Homologue 3 (RSH3), a nuclear-encoded chloroplastic enzyme responsible for ppGpp biosynthesis, interacts directly with the TOR complex via a plant-specific N-terminal region (NTR) which is hyper-phosphorylated in a TOR-dependent manner. Downregulation of TOR activity reduces NTR phosphorylation, enhances ppGpp synthesis by RSH3, and causes a ppGpp-dependent decrease in photosynthetic capacity. Altogether we demonstrate that the TOR-RSH3 signaling axis is a novel and direct post-translational mechanism that allows chloroplast activity to be matched with plant growth, setting a new precedent for the regulation of organellar function by TOR.One sentence summaryThe TOR kinase post-translationally controls guanosine tetraphosphate signaling to regulate plant photosynthetic activity.

Protein is biology workhorse. Since the recent break-through of novel folding methods, the amount of available structural data is increasing, closing the gap between data-driven sequence-based and structure-based methods. In this work, we focus on the inverse folding problem that consists in predicting an amino-acid primary sequence from protein 3D structure. For this purpose, we introduce a simple Transformer model from Natural Language Processing augmented 3D-structural data. We call the resulting model PeTriBERT: Proteins embedded in tridimensional representation in a BERT model. We train this small 40-million parameters model on more than 350 000 proteins sequences retrieved from the newly available AlphaFoldDB database. Using PetriBert, we are able to in silico generate totally new proteins with a GFP-like structure. These 9 of 10 of these GFP structural homologues have no ressemblance when blasted on the whole entry proteome database. This shows that PetriBert indeed capture protein folding rules and become a valuable tool for de novo protein design.

BackgroundHigh-throughput phenotyping is crucial for the genetic and molecular understanding of adaptive root system development. In recent years, imaging automata have been developed to acquire the root system architecture of many genotypes grown in Petri dishes to explore the Genetic x Environment (GxE) interaction. There is now an increasing interest in understanding the dynamics of the adaptive responses, such as the organ apparition or the growth rate. However, due to the increasing complexity of root architectures in development, the accurate description of the topology, geometry, and dynamics of a growing root system remains a challenge.ResultsWe designed a high-throughput phenotyping method, combining an imaging device and an automatic analysis pipeline based on registration and topological tracking, capable of accurately describing the topology and geometry of observed root systems in 2D+t. The method was tested on a challenging Arabidopsis seedling dataset, including numerous root occlusions and crossovers. Static phenes are estimated with high accuracy (R2 = 0.996 and 0, 923 for primary and second-order roots length, respectively). These performances are similar to state-of-the-art results obtained on root systems of equal or lower complexity. In addition, our pipeline estimates dynamic phenes accurately between two successive observations (R2 = 0. 938 for lateral root growth).ConclusionsWe designed a novel method of root tracking that accurately and automatically measures both static and dynamic RSA parameters from a novel high-throughput root phenotyping platform. It has been used to characterize developing patterns of root systems grown under various environmental conditions. It provides a solid basis to explore the GxE interaction controlling the dynamics of root system architecture adaptive responses. In future work, our approach will be adapted to a wider range of imaging configurations and species.

Guanosine pentaphosphate and tetraphosphate (together referred to as ppGpp) are hyperphosphorylated nucleotides found in bacteria and the chloroplasts of plants and algae. In plants and algae artificial ppGpp accumulation can inhibit chloroplast gene expression, and influence photosynthesis, nutrient remobilisation, growth, and immunity. However, it is so far unknown whether ppGpp is required for abiotic stress acclimation in plants. Here, we demonstrate that ppGpp biosynthesis is necessary for acclimation to nitrogen starvation in Arabidopsis. We show that ppGpp is required for remodeling the photosynthetic electron transport chain to downregulate photosynthetic activity and for protection against oxidative stress. Furthermore, we demonstrate that ppGpp is required for coupling chloroplastic and nuclear gene expression during nitrogen starvation. Altogether, our work indicates that ppGpp is a pivotal regulator of chloroplast activity for stress acclimation in plants.

Background The sequencing of the wheat (Triticum aestivum) genome has been a methodological challenge for many years owing to its large size (15.5 Gb), repeat content, and hexaploidy. Many initiatives aiming at obtaining a reference genome of cultivar Chinese Spring have been launched in the past years and it was achieved in 2018 as the result of a huge effort to combine short-read sequencing with many other resources. Reference-quality genome assemblies were then produced for other accessions, but the rapid evolution of sequencing technologies offers opportunities to reach high-quality standards at lower cost. Results Here, we report on an optimized procedure based on long reads produced on the Oxford Nanopore Technology PromethION device to assemble the genome of the French bread wheat cultivar Renan. Conclusions We provide the most contiguous chromosome-scale assembly of a bread wheat genome to date. Coupled with an annotation based on RNA-sequencing data, this resource will be valuable for the crop community and will facilitate the rapid selection of agronomically important traits. We also provide a framework to generate high-quality assemblies of complex genomes using ONT.

Viral infections impose extraordinary RNA stress on a cell, triggering cellular RNA surveillance pathways like RNA decapping, nonsense-mediated decay and RNA silencing. Viruses need to maneuver between these pathways to establish infection and succeed in producing high amounts of viral proteins. Processing bodies (PBs) are integral to RNA triage in eukaryotic cells with several distinct RNA quality control pathways converging for selective RNA regulation. In this study, we investigate the role of Arabidopsis thaliana PBs during Cauliflower Mosaic Virus (CaMV) infection. We find that several PB components are co-opted into viral replication factories and support virus multiplication. This pro-viral role was not associated with RNA decay pathways but instead, we could establish PB components as essential helpers in viral RNA translation. While CaMV is normally resilient to RNA silencing, PB dysfunctions expose the virus to this pathway, similar to previous observations on transgenes. Transgenes, however, undergo RNA Quality Control dependent RNA degradation, whereas CaMV RNA remains stable but becomes translationally repressed through decreased ribosome association, revealing a unique dependence between PBs, RNA silencing and translational repression. Together, our study shows that PB components are co-opted by the virus to maintain efficient translation, a mechanism not associated with canonical PB functions.

Since modern wheat varieties are grown with chemical inputs, we ignore if changes observed in rhizosphere microorganisms between ancient and modern varieties are due to i) breeding-induced changes in plant genotype, ii) modifications of the environment via synthetic chemical inputs, or (iii) phenotypic plasticity, defined as the interaction between the genotype and the environment. In the field, we evaluated the effects of various wheat varieties (modern and ancient) grown with or without chemical inputs (N-fertilizer, fungicide and herbicide together) in a crossed factorial design. We analysed rhizosphere bacteria and fungi by amplicons sequencing and mycorrhizal association by microscopic observations. When considered independently of plant genotype, chemical inputs were responsible for an increase in dominance for bacteria and decrease in evenness for bacteria and fungi. Independently of inputs, modern varieties had richer and more even bacterial communities compared to ancient varieties. Phenotypic plasticity had a significant effect: bacterial and fungal diversity decreased when inputs were applied in ancient varieties but not in modern ones. Mycorrhiza were more abundant in modern than ancient varieties, and less abundant when using chemical inputs. Although neglected, phenotypic plasticity is important to understand the evolution of plant-microbiota associations and a relevant target in breeding programs.

SummaryThe bimolecular fluorescence complementation (BiFC) assay has emerged as one of the most popular methods for analysing protein-protein interactions (PPIs) in plant biology. This includes its increasing use as a tool for dissecting the molecular mechanisms of chloroplast function. However, the construction of chloroplast fusion proteins for BiFC can be difficult, and the availability and selection of appropriate controls is not trivial. Furthermore, the challenges of performing BiFC in restricted cellular compartments has not been specifically addressed. Here we describe the development of a flexible modular cloning-based toolkit (MoBiFC) for chloroplast BiFC and proximity labelling using synthetic biology principles. The approach facilitates the cloning process for chloroplast-targeted proteins, allows robust ratiometric quantification, and the toolkit comes with model positive and negative controls. Our study also highlights many potential pitfalls including the choice of fluorescent protein (FP) split, negative controls, cell type, and reference FP. Finally, we provide an example of how users can enrich the toolset by providing functional proximity labelling modules, and we discuss how MoBiFC could be further improved and extended to other compartments of the plant cell.

Iron deficiency hampers photosynthesis and is associated with chlorosis. We recently showed that iron deficiency-induced chlorosis depends on phosphorus availability. How plants integrate these cues to control chlorophyll accumulation is unknown. Here, we show that iron limitation downregulates photosynthesis genes in a phosphorus-dependent manner. Using transcriptomics and genome-wide association analysis, we identify two genes, a chloroplastic ascorbate transporter (PHT4;4) and a nuclear transcription factor (bZIP58), which prevent the downregulation of photosynthesis genes leading to the stay-green phenotype under iron-phosphorus deficiency. Joint limitation of these nutrients induces ascorbate accumulation by activating expression of an ascorbate biosynthesis gene, VTC4, which requires bZIP58. Exogenous ascorbate prevents iron deficiency-induced chlorosis in vtc4 mutants, but not in bzip58 or pht4;4. Our study demonstrates chloroplastic ascorbate transport is essential for preventing the downregulation of photosynthesis genes under iron-phosphorus combined deficiency. These findings uncover a molecular pathway coordinating chloroplast-nucleus communication to adapt photosynthesis to nutrient availability.

International audience; Oryza sativa (rice) plays an essential food security role for more than half of the world's population. Obtaining crops with high levels of disease resistance is a major challenge for breeders, especially today, given the urgent need for agriculture to be more sustainable. Plant resistance genes are mainly encoded by three large leucine-rich repeat (LRR)-containing receptor (LRR-CR) families: the LRR-receptor-like kinase (LRR-RLK), LRR-receptor-like protein (LRR-RLP) and nucleotide-binding LRR receptor (NLR). Using LRRPROFILER, a pipeline that we developed to annotate and classify these proteins, we compared three publicly available annotations of the rice Nipponbare reference genome. The extended discrepancies that we observed for LRR-CR gene models led us to perform an in-depth manual curation of their annotations while paying special attention to nonsense mutations. We then transferred this manually curated annotation to Kitaake, a cultivar that is closely related to Nipponbare, using an optimized strategy. Here, we discuss the breakthrough achieved by manual curation when comparing genomes and, in addition to 'functional' and 'structural' annotations, we propose that the community adopts this approach, which we call 'comprehensive' annotation. The resulting data are crucial for further studies on the natural variability and evolution of LRR-CR genes in order to promote their use in breeding future resilient varieties.

Transcriptomic changes at the cessation of sugar accumulation in the pericarp of Vitis vinifera were addressed on single berries re-synchronized according to their individual growth patterns. The net rates of water, sugars and K+ accumulation inferred from individual growth and solute concentration confirmed that these inflows stopped simultaneously in the ripe berry, while the small amount of malic acid remaining at this stage was still being oxidized at low rate. Re-synchronized individual berries displayed negligible variations in gene expression among triplicates. RNA-Seq studies revealed sharp reprogramming of cell wall enzymes and structural proteins at the stop of phloem unloading, associated with an 80% repression of multiple sugar transporters and aquaporins on the plasma or tonoplast membranes, with the noticeable exception of H+/sugar symporters, that were rather weakly and constitutively expressed. This was verified in three genotypes placed in contrasted thermo-hydric conditions. The prevalence of SWEET suggests that electrogenic transporters would play a minor role on the plasma membranes of SE/CC complex and the one of the flesh, while sucrose/H+ exchangers dominate on its tonoplast. Cis-regulatory elements present in their promoters allowed to sort these transporters in different groups, also including specific TIPs and PIPs paralogs, and cohorts of cell wall related genes. Together with simple thermodynamic considerations, these results lead to propose that H+/sugar exchangers at the tonoplast, associated with a considerably acidic vacuolar pH, may exhaust cytosolic sugars in the flesh and alleviate the need for supplementary energization of sugar transport at the plasma membrane.

SummaryThe nucleotides guanosine tetraphosphate and pentaphosphate (or ppGpp) are implicated in the regulation of chloroplast function in plants. ppGpp signalling is best understood in the model vascular plant Arabidopsis thaliana where it acts to regulate plastid gene expression to influence photosynthesis, plant development and immunity. However, little is known about the conservation or diversity of ppGpp signaling in other land plants. Here, we studied the function of ppGpp in the moss Physcomitrium (previously Physcomitrella) patens using an inducible system for triggering ppGpp accumulation. We used this approach to investigate the effects of ppGpp on chloroplast function, photosynthesis and growth. We demonstrate that ppGpp accumulation causes a dramatic drop in photosynthetic capacity by inhibiting chloroplast gene expression. This was accompanied by the unexpected reorganisation of the thylakoid system into super grana. Surprisingly, these changes did not affect gametophore growth, suggesting that bryophytes and vascular plants may have different tolerances to defects in photosynthesis. Our findings point to the existence of both highly conserved and more specific targets of ppGpp signalling in the land plants that may reflect different growth strategies.

Among legumes (Fabaceae) capable of nitrogen-fixing nodulation, several Aeschynomene spp. use a unique symbiotic process that is independent of Nod factors and infection threads. They are also distinctive in developing root and stem nodules with photosynthetic bradyrhizobia. Despite the significance of these symbiotic features, their understanding remains limited. To overcome such limitations, we conducted genetic studies of nodulation in Aeschynomene evenia, supported by the development of a genome sequence for A. evenia and transcriptomic resources for 10 additional Aeschynomene spp. Comparative analysis of symbiotic genes substantiated singular mechanisms in the early and late nodulation steps. A forward genetic screen also showed that AeCRK, coding a novel receptor-like kinase, and the symbiotic signaling genes AePOLLUX, AeCCamK, AeCYCLOPS, AeNSP2 and AeNIN, are required to trigger both root and stem nodulation. This work demonstrates the utility of the A. evenia model and provides a cornerstone to unravel new mechanisms underlying the rhizobium-legume symbiosis.

Shoot branching, a major component of shoot architecture, is regulated by multiple signals. Previous studies have indicated that sucrose may promote branching through suppressing the inhibitory effect of the hormone strigolactone (SL). However, the molecular mechanisms underlying this effect are unknown.- Here we used molecular and genetic tools to identify the molecular targets underlying the antagonistic interaction between sucrose and SL.- We showed that sucrose antagonises the suppressive action of SL on tillering in rice and on the degradation of D53, a major target of SL signalling. Sucrose inhibits the expression of D3, the orthologue of the arabidopsis F-box protein MAX2 required for SL signalling. Over-expression of D3 prevents sucrose from inhibiting D53 degradation and enabled the SL inhibition of tillering under high sucrose. Sucrose also prevents SL-induced degradation of D14, the SL receptor involved in D53 degradation. Interestingly, D14 over-expression enhances D53 protein levels and sucrose-induced tillering.- Our results show that sucrose inhibits SL perception by targeting key components of SL signalling and, together with previous studies reporting the inhibition of SL synthesis by nitrate and phosphate, demonstrate the central role played by strigolactones in the regulation of plant architecture by nutrients.

Fruit colour, which is central to the sensorial and nutritional quality of the fruit, is a major breeding target in cultivated strawberry (Fragaria × ananassa). Red fruit colour is caused by the accumulation of anthocyanins, which are water-soluble flavonoids. Here, using pseudo F1progeny derived from the cv. ‘Capitola’ and the advanced line ‘CF1116’ and taking advantage of the available high density SNP array, we delineated fruit flavonoids mQTLs (13 compounds: anthocyanin, flavonols and flavan-3-ols) and colour-related QTLs (total anthocyanins and fruit colour) to narrow genomic regions corresponding to specific subgenomes of the cultivated strawberry. Findings showed that the overwhelming majority of the anthocyanin mQTLs and other colour-related QTLs are specific toF. vescasubgenome but that other subgenomes also contribute to colour variations. We then focused on two major homoeo-mQTLs for pelargonidin-3-glucoside (PgGs) localized on both male and female maps on linkage group LG3a (F. vescasubgenome). Combined high-resolution mapping of PgGs mQTLs and colour QTLs, transcriptome analysis of selected progeny individuals and whole genome sequencing of the parents led to the identification of several INDELS in thecis-regulatory region of aMYB102-like ODORANTgene and of the deletion of a putative MADS box binding motif in the 5’UTR upstream region of an anthocyanidin reductase (ANR)gene, which likely underlie significant colour variations in strawberry fruit. The implications of these findings are important for the functional analysis and genetic engineering of colour-related genes and for the breeding by Marker-Assisted-Selection of new strawberry varieties with improved colour and health-benefits.

Enhancers are important regulators of gene expression during numerous crucial processes including tissue differentiation across development. In plants, their recent molecular characterization revealed their capacity to activate the expression of several target genes through the binding of transcription factors. Nevertheless, identifying these target genes at a genome-wide level remains a challenge, in particular in species with large genomes, where enhancers and target genes can be hundreds of kilobases away. Therefore, the contribution of enhancers to regulatory network is still poorly understood in plants. In this study, we investigate the enhancer-driven regulatory network of two maize tissues at different stages: leaves at seedling stage and husks (bracts) at flowering. Using a systems biology approach, we integrate genomic, epigenomic and transcriptomic data to model the regulatory relationship between transcription factors and their potential target genes. We identify regulatory modules specific to husk and V2-IST, and show that they are involved in distinct functions related to the biology of each tissue. We evidence enhancers exhibiting binding sites for two distinct transcription factor families (DOF and AP2/ERF) that drive the tissue-specificity of gene expression in seedling immature leaf and husk. Analysis of the corresponding enhancer sequences reveals that two different transposable element families (TIR transposon Mutator and MITE Pif/Harbinger) have shaped the regulatory network in each tissue, and that MITEs have provided new transcription factor binding sites that are involved in husk tissue-specificity.SignificanceEnhancers play a major role in regulating tissue-specific gene expression in higher eukaryotes, including angiosperms. While molecular characterization of enhancers has improved over the past years, identifying their target genes at the genome-wide scale remains challenging. Here, we integrate genomic, epigenomic and transcriptomic data to decipher the tissue-specific gene regulatory network controlled by enhancers at two different stages of maize leaf development. Using a systems biology approach, we identify transcription factor families regulating gene tissue-specific expression in husk and seedling leaves, and characterize the enhancers likely to be involved. We show that a large part of maize enhancers is derived from transposable elements, which can provide novel transcription factor binding sites crucial to the regulation of tissue-specific biological functions.

In the course of their growth and development plants have to constantly perceive and react to their environment. This is achieved in cells, by the coordination of complex combinatorial signaling networks. However, how signal integration and specificity are achieved in this context is unknown. With a focus on the hyperosmotic stimulus, we use live super-resolution light imaging methods to demonstrate that a Rho GTPase, Rho-of-Plant 6 (ROP6), forms stimuli-dependant nanodomains within the PM. These nanodomains are necessary and sufficient to transduce production of reactive oxygen species (ROS),that act as secondary messengers and trigger several plant adaptive responses to osmotic constraints. Furthermore, ROP6 activation triggers the nanoclustering of two NADPH oxidases that subsequently generates ROS. ROP6 nanoclustering is also needed for cell surface auxin signaling, but short-time auxin treatment does not induce ROS accumulation. We show that auxin-induced ROP6 nanodomains, unlike osmotically-driven ROP6 clusters, do not recruit the NADPH oxidase, RBOHD. Together, our results suggest that Rho GTPase nano-partitioning at the PM ensure signal specificity downstream of independent stimuli.

Horizontal gene transfer is of major evolutionary importance as it allows for the redistribution of phenotypically important genes among lineages. Such genes with essential functions include those involved in resistance to antimicrobial compounds and virulence factors in pathogenic bacteria. Understanding gene turnover at microevolutionary scales is critical to assess the pace of this evolutionary process. Here, we characterized and quantified gene turnover for the epidemic lineage of a bacterial plant pathogen of major agricultural importance worldwide. Relying on a dense geographic sampling spanning 39 years of evolution, we estimated both the dynamics of single nucleotide polymorphism accumulation and gene content turnover. We identified extensive gene content variation among lineages even at the smallest phylogenetic and geographic scales. Gene turnover rate exceeded nucleotide substitution rate by three orders of magnitude. Accessory genes were found preferentially located on plasmids, but we identified a highly plastic chromosomal region hosting ecologically important genes such as transcription activator-like effectors. Whereas most changes in the gene content are probably transient, the rapid spread of a mobile element conferring resistance to copper compounds widely used for the management of plant bacterial pathogens illustrates how some accessory genes can become ubiquitous within a population over short timeframes.

17 Pág., A decade of genetic association studies in multiple organisms suggests that most complex traits are polygenic; that is, they have a genetic architecture determined by numerous loci, each with small effect-size. Thus, determining the degree of polygenicity and its variation across traits, environments and time is crucial to understand the genetic basis of phenotypic variation. We applied multilocus approaches to estimate the degree of polygenicity of fitness-related traits in a long-lived plant (Pinus pinaster Ait., maritime pine) and to analyse this variation across environments and years. We evaluated five categories of fitness-related traits (survival, height, phenology, functional, and biotic-stress response) in a clonal common-garden network planted in contrasted environments (over 20,500 trees). Most of the analysed traits showed evidence of local adaptation based on Qst -Fst comparisons. We further observed a remarkably stable degree of polygenicity, averaging 6% (range of 0%-27%), across traits, environments and years. We detected evidence of negative selection, which could explain, at least partially, the high degree of polygenicity. Because polygenic adaptation can occur rapidly, our results suggest that current predictions on the capacity of natural forest tree populations to adapt to new environments should be revised, especially in the current context of climate change., This study was funded by the Spanish Ministry of Economy and Competitiveness through projects RTA2010-00120-C02- 02 (CLONAPIN), CGL2011-30182- C02- 01 (AdapCon) and AGL2012-40151- C03- 02 (FENOPIN). The study was also supported by the “Initiative d’Excellence (IdEx) de l’Université de Bordeaux - Chaires d'installation 2015” (EcoGenPin) and the European Union's Horizon 2020 research and innovation programme under grant agreement No 773383 (B4EST)

NRT2.1 can be phosphorylated at five different sites within N- and C-terminus. Here, we provide a systematic functional characterization of phosphorylation at S21 and S28 within the N-terminus of NRT2.1. We used existing phosphoproteomic data sets of nitrate starvation and nitrate resupply to construct a site-specific correlation network identifying kinase candidates to phosphorylate NRT2.1. By this approach, we identified NITRATE UPTAKE REGULATORY KINASE 1 (AT5G49770) which itself was regulated by phosphorylation at S839 and S870 within its kinase domain. In the active state, when S839 was dephosphorylated and S870 was phosphorylated, NURK1 was found to interact with NRT2.1 at dephosphorylated S28. Upon that interaction, NURK1 can phosphorylate NRT2.1 at S21. Phosphorylation of NRT2.1 at S21 resulted in low interaction of NRT2.1 with its activator protein NAR2.1. By contrast, phosphorylation of NRT2.1 at S28 by a yet unknown kinase enhanced the interaction with NAR2.1, but inhibited the interaction with NURK1. We propose that serines S21 and S28 are involved in a phospho-switch mechanism and by which the interaction of NRT2.1 with its activator NAR2.1, and thus NRT2.1 activity, is modulated. NURK1 here was identified as the kinase affecting this phospho-switch through phosphorylation of NRT2.1 at S21 leading to inactivation of NRT2.1.

The preservation of the maximum of diversity within the smallest number of accessions is one of the challenges of germplasm management. To construct core-collections, the assessment of the population structure and the relationships between the accessions represents a key step and the choice of suitable molecular markers is the starting point. Since the expansion of available SNP-based genomics tools, a debate has emerged regarding the usefulness of the widely used microsatellites (SSRs) markers. In this study, we analysed a part of the INRAE walnut germplasm collection of 150 accessions, unique in Europe for walnut biodiversity conservation, by comparing the power of both types of marker. We found that the first level of structure is equally detected using 13 SSRs or the Axiom™ J. regia 700K SNP array, and is in relation with the geographical origin of the accessions. For K=2, there was no exchange of accession between the two groups when both markers were compared. We also highlighted empirically that approximately 100 SNPs are needed to obtain similar clustering to SSRs in Principal Coordinate Analysis (PCoA). The neighbor-joining trees constructed were also consistent between both types of marker. The main differences lied in the upper levels of structure from K=3 to K=6, more powerful using the SNPs, and in the percentage of the explained variation in PCoA for K=2, higher using SSRs. We then constructed core-collections of 50 accessions, a crucial step in genetic resources management to reduce the costs and preserve the allelic diversity. Using two different construction methods, both SSR and SNP markers were suitable and able to keep at least 88.57% of the alleles. 32/50 accessions were in common between the two markers, for both methods. We concluded that the use of either marker is dependent on the researcher’s goal.

BackgroundParasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biocontrol, and infecting aphids requires overcoming both aphid defenses and their defensive endosymbionts.ResultsWe present thede novogenome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids:Aphidius erviandLysiphlebus fabarum(Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp), highly syntenic, and the most AT-rich reported thus far for any arthropod (GC content: 25.8% and 23.8%). This nucleotide bias is accompanied by skewed codon usage, and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and age-specific energy demands. We identify expansions of F-box/Leucine-rich-repeat proteins, suggesting that diversification in this gene family may be associated with their broad host range or with countering defenses from aphids’ endosymbionts. The absence of some immune genes (Toll and Imd pathways) resembles similar losses in their aphid hosts, highlighting the potential impact of symbiosis on both aphids and their parasitoids.ConclusionsThese findings are of fundamental interest for insect evolution and beyond. This will provide a strong foundation for further functional studies including coevolution with respect to their hosts, the basis of successful infection, and biocontrol. Both genomes are available athttps://bipaa.genouest.org.

The vast majority of eukaryotic cells contain mitochondria, essential powerhouses and metabolic hubs1. These organelles have a bacterial origin and were acquired during an early endosymbiosis event2. Mitochondria possess specialized gene expression systems composed of various molecular machines including the mitochondrial ribosomes (mitoribosomes). Mitoribosomes are in charge of translating the few essential mRNAs still encoded by mitochondrial genomes3. While chloroplast ribosomes strongly resemble those of bacteria4,5, mitoribosomes have diverged significantly during evolution and present strikingly different structures across eukaryotic species6–10. In contrast to animals and trypanosomatides, plants mitoribosomes have unusually expanded ribosomal RNAs and conserved the short 5S rRNA, which is usually missing in mitoribosomes11. We have previously characterized the composition of the plant mitoribosome6revealing a dozen plant-specific proteins, in addition to the common conserved mitoribosomal proteins. In spite of the tremendous recent advances in the field, plant mitoribosomes remained elusive to high-resolution structural investigations, and the plant-specific ribosomal features of unknown structures. Here, we present a cryo-electron microscopy study of the plant 78S mitoribosome from cauliflower at near-atomic resolution. We show that most of the plant-specific ribosomal proteins are pentatricopeptide repeat proteins (PPR) that deeply interact with the plant-specific rRNA expansion segments. These additional rRNA segments and proteins reshape the overall structure of the plant mitochondrial ribosome, and we discuss their involvement in the membrane association and mRNA recruitment prior to translation initiation. Finally, our structure unveils an rRNA-constructive phase of mitoribosome evolution across eukaryotes.

Plants use nitrate and ammonium in the soil as their main nitrogen sources. Recently, ammonium has attracted attention due to evidence suggesting that, in C3 species, an elevated CO2 environment inhibits nitrate assimilation. However, high concentrations of ammonium as the sole nitrogen source for plants causes impaired growth, i.e. ammonium toxicity. Although ammonium toxicity has been studied for a long time, the primary cause remains to be elucidated. Here, we show that ammonium assimilation in plastids rather than ammonium accumulation is a primary cause for toxicity. Our genetic screen of ammonium-tolerant Arabidopsis lines with enhanced shoot growth identified plastidic GLUTAMINE SYNTHETASE 2 (GLN2) as the causal gene. Our reciprocal grafting of wild-type and GLN2 or GLN1;2-deficient lines suggested that shoot GLN2 activity results in ammonium toxicity, whilst root GLN1;2 activity prevents it. With exposure to toxic levels of ammonium, the shoot GLN2 reaction produced an abundance of protons within cells, thereby elevating shoot acidity and stimulating expression of acidic stress-responsive genes. Application of an alkaline ammonia solution to the toxic ammonium medium efficiently alleviated the ammonium toxicity with a concomitant reduction in shoot acidity. Consequently, we conclude that a primary cause of ammonium toxicity is acidic stress in the shoot. This fundamental insight provides a framework for enhanced understanding of ammonium toxicity in plants.

High-throughput sequencing (HTS), combined with bioinformatics forde novodiscovery and assembly of plant virus or viroid genome reads, has promoted the discovery of abundant novel DNA and RNA viruses and viroids. However, the elucidation of a viral population in a single plant is rarely reported. In five birch trees of German and Finnish origin exhibiting symptoms of birch leaf-roll disease (BRLD), we identified in total five viruses, among which three are novel. The number of identified virus variants in each transcriptome ranged from one to five. The novel species are genetically - fully or partially - characterized, they belong to the generaCarlavirus, IdaeovirusandCapillovirusand they are tentatively namedbirch carlavirus,birch idaeovirus, andbirch capillovirus, respectively. The only virus systematically detected by HTS in symptomatic trees affected by the BRLD was the recently discovered birch leafroll-associated virus. The role of the new carlavirus in BLRD etiology seems at best weak, as it was detected only in one of three symptomatic trees. Continuing studies have to clarify the impact of the carlavirus to the BLRD. The role of theCapillovirusand theIdaeoviruswithin the BLRD complex and whether they influence plant vitality need to be investigated. Our study reveals the viral population in single birch trees and provides a comprehensive overview for the diversities of the viral communities they harbor.

Grapevine fanleaf virus (GFLV) is a picorna-like plant virus transmitted by nematodes that affects vineyards worldwide. Nanobody (Nb)-mediated resistance against GFLV has been created recently and shown to be highly effective in plants including grapevine, but the underlying mechanism is unknown. Here we present the high-resolution cryo-EM structure of the GFLV-Nb23 complex which provides the basis for the molecular recognition by the nanobody. The structure reveals a composite binding site bridging over 3 domains of the capsid protein (CP) monomer. The structure provides a precise mapping of the Nb23 epitope on the GFLV capsid in which the antigen loop is accommodated through an induced fit mechanism. Moreover, we uncover and characterize several resistance-breaking GFLV isolates with amino acids mapping within this epitope, including C-terminal extensions of the CP, which would sterically interfere with Nb binding. Escape variants with such extended CP fail to be transmitted by nematodes linking Nb-mediated resistance to vector transmission. Together, these data provide insights into the molecular mechanism of Nb23-mediated recognition of GFLV and of virus resistance loss.SignificanceGrapevine fanleaf virus (GFLV) is a picorna-like plant virus that severely impacts vineyards worldwide. While Nanobodies (Nb) confer resistance to GFLV in plants the underlying molecular mechanism of action is unknown. Here we present the high-resolution cryo-EM structure of the GFLV-Nb complex. It uncovers the conformational epitope on the capsid surface which is a composite binding site into which the antigen loop is accommodated through an induced fit mechanism. Furthermore, we describe several resistance-breaking isolates of GFLV with reduced Nb binding capacity. Those that carry a C-terminal extension also fail to be transmitted by nematodes. Together, these data provide structure-function insights into the Nb-GFLV recognition and the molecular mechanism leading to loss of resistance.

Glutathione-dependent root growth in Arabidopsis is linked to polar auxin transport (PAT). Arabidopsis mutants with reduced glutathione (GSH) levels also show reduced PAT. To gain an insight into the relationship between PAT and GSH level, we analysed tomato polycotyledon mutant, pct1-2, which has enhanced PAT. Microarray analysis of gene expression in pct1-2 mutant revealed underexpression of several genes related to glutamate and glutathione metabolism. In consonance with microarray analysis, enzymatic as well as in-vivo assay revealed higher glutathione level in the early phase of pct1-2 seedling growth than WT. The inhibition of auxin transport by 2,3,5-triiodobenzoic acid (TIBA) reduced both GSH level and PIN1 expression in pct1-2 root tips. The reduction of in vivo GSH accumulation in pct1-2 root tips by buthionine sulfoximine (BSO) stimulated elongation of the short root of pct1-2 mutant akin to TIBA. The rescue of short root phenotype of pct1-2 mutant was restricted to TIBA and BSO. The other auxin transport inhibitors 1-N-naphthylphthalamic acid (NPA), 2-[4-(diethylamino)-2-hydroxybenzoyl] benzoic acid (BUM), 3-chloro-4-hydroxyphenylacetic acid (CHPAA), brefeldin and gravacin inhibited root elongation in both WT and pct1-2 mutant. Our results indicate a relationship between PAT and GSH level in tomato akin to Arabidopsis. Our work also highlights that TIBA rescues short root phenotype of the pct1-2 mutant by acting on a PAT component distinct from the site of action of other PAT inhibitors.

Soybean (Glycine max (L.) Merr.) has potential to improve sustainability of agricultural production systems. A higher focus on this crop is needed to re-launch its production in the EU. A better understanding of key determinants affecting soybean establishment represents a first step to facilitate its adoption in cropping systems. To this objective, we conducted laboratory and field experiments in order to better characterize seed germination and seedling growth in relation to temperatures, water content, and soil structure. We then used these data to parametrize the SIMPLE crop emergence model and to evaluate its prediction quality, by comparing observed field germination and emergence data with the predicted ones. Finally, we performed a simulation study over the 2020-2100 period, for three sowing dates, from mid-March to mid-April, in the northern climate of France to evaluate whether future climate change will help expand soybean from Southern to Northern part of the country. Soybean germination was very fast, taking only 15 °C days to reach 50% germination at optimal conditions. The base, optimum and maximum temperatures were determined as 4, 30 and 40°C, respectively while the base water potential was −0.7 MPa, indicating a high sensitivity to water stress. The SIMPLE model well-predicted germination and emergence courses and their final rates, compared with the observed field data. The simulation study showed average emergence rate ranging from 61 to 78% with little variability among sowing dates and periods, but a high variability between years. Main causes of non-emergence were seedling mortality due to clods or soil surface crust followed by non-germination and seedling mortality due to drought, especially for mid-April sowing. These results provide a better knowledge of soybean establishment that are encouraging to introduce soybean with early sowings to diversify current cropping systems.

Root-knot nematodes (genusMeloidogyne) are plant parasitic species that cause huge economic loss in the agricultural industry and affect the prosperity of communities in developing countries. Control methods against these plant pests are sparse and the current preferred method is deployment of plant cultivars bearing resistance genes againstMeloidogynespecies. However, some species such asM. enterolobiiare not controlled by the resistance genes deployed in the most important crop plants cultivated in Europe. The recent identification of this species in Europe is thus a major concern. Like the other most damaging Meloidogyne species (e.g.M. incognita,M. arenariaandM. javanica),M. enterolobiireproduces by obligatory mitotic parthenogenesis. Genomic singularities such as a duplicated genome structure and a relatively high proportion of transposable elements have previously been described in the above mentioned mitotic parthenogenetic Meloidogyne.To gain a better understanding of the genomic and evolutionary background we sequenced the genome ofM. enterolobiiusing high coverage short and long read technologies. The information contained in the long reads helped produce a highly contiguous genome assembly ofM. enterolobii, thus enabling us to perform high quality annotations of coding and non-coding genes, and transposable elements.The genome assembly and annotation reveals a genome structure similar to the ones described in the other mitotic parthenogenetic Meloidogyne, described as recent hybrids. Most of the genome is present in 3 different copies that show high divergence. Because most of the genes belong to these duplicated regions only few gene losses took place, which suggest a recent polyploidization. The most likely hypothesis to reconcile high divergence between genome copies despite few gene losses and translocations is also a recent hybrid origin. Consistent with this hypothesis, we found an abundance of transposable elements at least as high as the one observed in the mitotic parthenogenetic nematodesM. incognitaandM. javanica.

Pearl millet, unlike other cereals, is able to withstand dry and hot conditions and plays an important role for food security in arid and semi-arid areas of Africa and India. However, low soil fertility and drought constrain pearl millet yield. One of the main targets to address these constraints through agricultural practices or breeding is root system architecture. In this study, in order to easily phenotype the root system in field conditions, we developed a model to predict root length density (RLD) of pearl millet plants from root intersection densities (RID) counted on a trench profile in field conditions. We identified root orientation as an important parameter to improve the relationship between RID and RLD. Root orientation was notably found to differ between thick roots (more anisotropic with depth) and fine roots (isotropic at all depths). We used our model to study pearl millet root system response to drought and showed that pearl millet reorients its root growth toward deeper soil layers that retain more water in these conditions. Overall, this model opens ways for the characterization of the impact of environmental factors and management practices on pearl millet root system development.

The basidiomycete red yeast Rhodotorula toruloides is a promising platform organism for production of biooils. We present rhto-GEM, the first genome-scale model of R. toruloides metabolism, that was largely reconstructed using RAVEN toolbox. The model includes 4869 genes, 897 reactions, and 3334 metabolites. Lipid metabolism was described using the SLIMEr formalism, which allows direct integration of lipid class and acyl chain experimental distribution data. The simulation results confirmed that the R. toruloides model provides valid growth predictions on glucose, xylose and glycerol, while prediction of genetic engineering targets to increase production of linolenic acid and triacylglycerols highlighted genes that have previously been successfully used to increase lipid production. This renders rtho-GEM useful for future studies to improve the production of other industrially important oleochemicals including both value-added fatty acids and carotenoids, while it will also be valuable tool for system-wide omics-data analysis in R. toruloides. Expanding the portfolio of GEMs for lipid accumulating fungi contributes to both understanding of metabolic mechanisms of the oleaginous phenotype but also uncover particularities of lipid production machinery in R. toruloides.

Whole genome profiling (WGP) is a sequence-based physical mapping technology and uses sequence tags generated by next generation sequencing for construction of bacterial artificial chromosome (BAC) contigs of complex genomes. The physical map provides a framework for assembly of genome sequence and information for localization of genes that are difficult to find through positional cloning. To address the challenges of accurate assembly of the pea genome (~4.2 GB of which approximately 85% is repetitive sequences), we have adopted the WGP technology for assembly of a pea BAC library. Multi-dimensional pooling of 295,680 BAC clones and sequencing the ends of restriction fragments of pooled DNA generated 1,814 million high quality reads, of which 825 million were deconvolutable to 1.11 million unique WGP sequence tags. These WGP tags were used to assemble 220,013 BACs into contigs. Assembly of the BAC clones using the modified Fingerprinted Contigs (FPC) program has resulted in 13,040 contigs, consisting of 213,719 BACs, and 6,294 singleton BACs. The average contig size is 0.33 Mbp and the N50 contig size is 0.62 Mbp. WGPTM technology has proved to provide a robust physical map of the pea genome, which would have been difficult to assemble using traditional restriction digestion based methods. This sequence-based physical map will be useful to assemble the genome sequence of pea. Additionally, the 1.1 million WGP tags will support efficient assignment of sequence scaffolds to the BAC clones, and thus an efficient sequencing of BAC pools with targeted genome regions of interest.

In Arabidopsis thaliana, canonical auxin-dependent gene regulation is mediated by 23 transcription factors from the AUXIN RESPONSE FACTOR (ARF) family, most of which interact with 29 auxin/indole acetic acid (Aux/IAA) repressors, themselves forming, in the presence of auxin, coreceptor complexes with one of six TRANSPORT INHIBITOR1/AUXIN-SIGNALLING F-BOX PROTEINS (TIR1/AFB). Different combinations of co-receptors drive specific sensing outputs, allowing auxin to control a myriad of processes. Considerable efforts have been made to discern the temporal and spatial specificity of auxin action. However, owing to a lack of obvious phenotype in single loss-of-function mutants in Aux/IAA genes, most genetic studies have relied on gain-of-function mutants, which are highly pleiotropic. In this article, we describe a molecular framework for the role of several members of the auxin sensing machinery. Using loss-of-function mutants, we demonstrate that TIR1 and AFB2 are positive regulators, whereas IAA6, IAA9 and IAA17 are negative regulators of adventitious root (AR) formation. The three Aux/IAA proteins interact with ARF6 and/or ARF8, which we have previously shown to be positive regulators of AR formation upstream of jasmonate, and likely repress their activity. Our data also suggest a dual role for TIR1 in the control of JA biosynthesis and conjugation, as revealed by upregulation of several JA biosynthesis genes in the tir1-1 mutant. In conclusion, we propose that in the presence of auxin, TIR1 and AFB2 form specific sensing complexes with IAA6, IAA9 and/or IAA17 that modulate JA homeostasis to control AR initiation.

Alternative splicing (AS) of pre-mRNAs in plants is an important mechanism of gene regulation in environmental stress tolerance but plant signals involved are essentially unknown. Pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) is mediated by mitogen-activated protein kinases and the majority of PTI defense genes are regulated by MPK3, MPK4 and MPK6. These responses have been mainly analyzed at the transcriptional level, however many splicing factors are direct targets of MAPKs. Here, we studied alternative splicing induced by the PAMP flagellin in Arabidopsis. We identified 506 PAMP-induced differentially alternatively spliced (DAS) genes. Although many DAS genes are targets of nonsense-mediated degradation (NMD), only 19% are potential NMD targets. Importantly, of the 506 PAMP-induced DAS genes, only 89 overlap with the set of 1849 PAMP-induced differentially expressed genes (DEG), indicating that transcriptome analysis does not identify most DASevents. Global DAS analysis of mpk3, mpk4, and mpk6 mutants revealed that MPK4 is a key regulator of PAMP-induced differential splicing, regulating AS of a number of splicing factors and immunity-related protein kinases, such as the calcium-dependent protein kinase CPK28, the cysteine-rich receptor like kinases CRK13 and CRK29 or the FLS2 co-receptor SERK4/BKK1. These data suggest that MAP kinase regulation of splicing factors is a key mechanism in PAMP-induced AS regulation of PTI.Significance statementAlternative pre-mRNA splicing (AS) affects plant responses to environmental stresses. So far, however, the regulation of AS is little understood. Here, we studied AS induced by the pathogen-associated molecular pattern (PAMP) flagellin in Arabidopsis. We identified 506 differentially alternatively spliced (DAS) genes, 89 of which overlap with the 1849 DEG, indicating that the majority of DAS events go undetected by common transcriptome analysis. PAMP-triggered immunity is mediated by mitogen-activated protein kinases. Global DAS analysis of MAPK mutants revealed that MPK4 is a key regulator of AS by affecting splicing factors and a number of important protein kinases involved in immunity. Since PAMP-triggered phosphorylation of several splicing factors is directly mediated by MAPKs, we discovered a key mechanism of AS regulation.

SummaryThe slow-fast continuum describes how resource allocation constrains life-history strategies in many organisms. In plants, it is reflected by a trade-off at the leaf level between the rate of carbon assimilation and lifespan, the so-called Leaf Economics Spectrum (LES). However, it is still unclear how the LES is connected to the slow-fast syndrome, and reflects adaptation to climate. Here, we measured growth, morpho-physiological and life-history traits at both leaf and whole-plant levels in 384 natural accessions of Arabidopsis thaliana. We examined the extent to which the LES continuum parallels the slow-fast continuum, and compared trait variation to neutral genetic differentiation between lineages. We found that the LES is tightly linked to variation in whole-plant functioning, relative growth rate and life history. A genetic analysis further suggested that phenotypic differentiation is linked to the evolution of different slow-fast strategies in contrasted climates. Together, our findings shed light on the physiological bases of the slow-fast continuum, and its role for plant adaptation to climate.

Predicting the usefulness of crosses in terms of expected genetic gain and genetic diversity is of interest to secure performance in the progeny and to maintain long-term genetic gain in plant breeding. A wide range of crossing schemes are possible including large biparental crosses, backcrosses, four-way crosses, and synthetic populations.In silicoprogeny simulations together with genome-based prediction of quantitative traits can be used to guide mating decisions. However, the large number of multi-parental combinations can hinder the use of simulations in practice. Analytical solutions have been proposed recently to predict the distribution of a quantitative trait in the progeny of biparental crosses using information of recombination frequency and linkage disequilibrium between loci. Here, we extend this approach to obtain the progeny distribution of more complex crosses including two to four parents. Considering agronomic traits and parental genome contribution as jointly multivariate normally distributed traits, the usefulness criterion parental contribution (UCPC) enables to (i) evaluate the expected genetic gain for agronomic traits, and at the same time (ii) evaluate parental genome contributions to the selected fraction of progeny. We validate and illustrate UCPC in the context of multiple allele introgression from a donor into one or several elite recipients in maize (Zea maysL.). Recommendations regarding the interest of two-way, three-way, and backcrosses were derived depending on the donor performance. We believe that the computationally efficient UCPC approach can be useful for mate selection and allocation in many plant and animal breeding contexts.

One sentence summaryExcessive N nutrition leads to ROS accumulation, and requires the function of major transcriptional regulators to maintain plants under physiological conditions.Author contributionsAn.M. and A.G. conceived research plans and supervised the experiments; F.B, Am.M., J.B., L.L., L.B. and An.M performed most of the experiments; F.B, Am.M., J.B., G.K., L.L., L.B. and An.M analyzed the data; An.M. wrote the article with contributions of all the authors.Competing interestsThe authors declare no competing financial interests.SummaryReactive Oxygen Species (ROS) can accumulate in cells at excessive levels, leading to unbalanced redox status and to a potential oxidative stress, which can have damaging effects to the molecular components of plant cells. Several environmental conditions have been described as causing an elevation of ROS production in plants. Consequently, this requires the expression of detoxification responses in order to maintain ROS homeostasis at physiological levels. In case of mis-regulation of the detoxification systems, oxidative stress can lead ultimately to growth retardation and developmental defects. Here, we demonstrate that Arabidopsis plants growing under high nitrogen environment have to express a set of genes involved in detoxification of ROS in order to maintain ROS at physiological levels. We show that the chromatin factor HNI9 is an important actor of this response, required for the expression of these detoxification genes. Mutation in HNI9 leads to elevated ROS levels, and to ROS-dependent phenotypic defects under high but not low N provision. In addition, we identify HY5 as one of the major transcription factors also required for the expression of this detoxification program under high N condition. Our results demonstrate the requirement of a balance between N nutrition and ROS production, and identified the first major regulators required to control ROS homeostasis under excessive N nutrition.

Rho GTPases are master regulators of cell signaling, but how they are regulated depending on the cellular context is unclear. Here, we show that the phospholipid phosphatidylserine acts as a developmentally-controlled lipid rheostat that tunes Rho GTPase signaling in Arabidopsis. Live super-resolution single molecule imaging revealed that RHO-OF-PLANT6 (ROP6) is stabilized by phosphatidylserine into plasma membrane (PM) nanodomains, which is required for auxin signaling. Furthermore, we uncovered that the PM phosphatidylserine content varies during plant root development and that the level of phosphatidylserine modulates the quantity of ROP6 nanoclusters induced by auxin and hence downstream signaling, including regulation of endocytosis and gravitropism. Our work reveals that variations in phosphatidylserine levels are a physiological process that may be leveraged to regulate small GTPase signaling during development.One Sentence SummaryPhosphatidylserine acts as a developmentally-controlled lipid rheostat that regulates cellular auxin sensitivity and plant development.

Allopolyploidy, combining interspecific hybridization with whole genome duplication, has had significant impact on plant evolution. Its evolutionary success is related to the rapid and profound genome reorganizations that allow neoallopolyploids to form and adapt. Nevertheless, how neoallopolyploid genomes adapt to regulate their expression remains poorly understood. The hypothesis of a major role for small noncoding RNAs (sRNAs) in mediating the transcriptional response of neoallopolyploid genomes has progressively emerged. Generally, 21-nt sRNAs mediate posttranscriptional gene silencing by mRNA cleavage, whereas 24-nt sRNAs repress transcription (transcriptional gene silencing) through epigenetic modifications. Here, we characterize the global response of sRNAs to allopolyploidy in Brassica, using three independently resynthesized Brassica napus allotetraploids originating from crosses between diploid Brassica oleracea and Brassica rapa accessions, surveyed at two different generations in comparison with their diploid progenitors. Our results suggest an immediate but transient response of specific sRNA populations to allopolyploidy. These sRNA populations mainly target noncoding components of the genome but also target the transcriptional regulation of genes involved in response to stresses and in metabolism; this suggests a broad role in adapting to allopolyploidy. We finally identify the early accumulation of both 21- and 24-nt sRNAs involved in regulating the same targets, supporting a posttranscriptional gene silencing to transcriptional gene silencing shift at the first stages of the neoallopolyploid formation. We propose that reorganization of sRNA production is an early response to allopolyploidy in order to control the transcriptional reactivation of various noncoding elements and stress-related genes, thus ensuring genome stability during the first steps of neoallopolyploid formation.

Theory predicts that environmental challenges can shape the composition of populations, which is manifest at the molecular level. Previously, we demonstrated that oil pollution affected gene expression patterns and altered genetic variation in natural populations of the foundation salt marsh grass, Spartina alterniflora. Here, we used a reduced representation bisulfite sequencing approach, epigenotyping by sequencing (epiGBS), to examine relationships among DNA sequence, DNA methylation, gene expression, and exposure to oil pollution. We documented genetic and methylation differentiation between oil-exposed and unexposed populations, suggesting that the Deepwater Horizon oil spill may have selected on genetic variation, and either selected on epigenetic variation or induced particular epigenotypes and expression patterns in exposed compared to unexposed populations. In support of the potential for differential response to the Deepwater Horizon oil spill, we demonstrate genotypic differences in response to oil under controlled conditions. Overall, these findings demonstrate genetic variation, epigenetic variation and gene expression are correlated to exposure to oil pollution, which may all contribute to the response to environmental stress.

Morphogenesis often yields organs with robust size and shapes, whereas cell growth and deformation feature significant spatio-temporal variability. Here, we investigate whether tissue responses to mechanical signals contribute to resolve this apparent paradox. We built a model of growing tissues made of fiber-like material, which may account for the cytoskeleton, polar cell-cell adhesion, or the extracellular matrix in animals, and for the cell wall in plants. We considered the synthesis and remodeling of this material, as well as the modulation of synthesis by isotropic and anisotropic response to mechanical stress. Formally, our model describes an expanding, mechanoresponsive, nematic, and active fluid. We show that mechanical responses buffer localized perturbations, with two possible regimes - hypo-responsive and hyper-responsive, and the transition between the two corresponds to a minimum value of the relaxation time. Whereas robustness of shapes suggests that growth fluctuations are confined to small scales, our model yields growth fluctuations that have long-range correlations. This indicates that growth fluctuations are a significant source of heterogeneity in development. Nevertheless, we find that mechanical responses may dampen such fluctuations, with a specific magnitude of anisotropic response that minimizes heterogeneity of tissue contours. We finally discuss how our predictions might apply to the development of plants and animals. Altogether, our results call for the systematic quantification of fluctuations in growing tissues.

In temperate trees, optimal timing and quality of flowering directly depend on adequate winter dormancy progression, regulated by a combination of chilling and warm temperatures. Physiological, genetic and functional genomic studies have shown that hormones play a key role in bud dormancy establishment, maintenance and release. We combined physiological, transcriptional analyses, quantification of abscisic acid (ABA) and gibberellins (GAs), and modelling to further investigate how these signaling pathways are associated with dormancy progression in the flower buds of two sweet cherry cultivars.Our results demonstrated that GA-associated pathways have distinct functions and may be differentially related with dormancy. In addition, ABA levels rise at the onset of dormancy, associated with enhanced expression of ABA biosynthesisPavNCEDgenes, and decreased prior to dormancy release. Following the observations that ABA levels are correlated with dormancy depth, we identifiedPavUG71B6, a sweet cherryUDP-GLYCOSYLTRANSFERASEgene that up-regulates active catabolism of ABA to ABA-GE and may be associated with low ABA content in the early cultivar. Subsequently, we modelled ABA content and dormancy behavior in three cultivars based on the expression of a small set of genes regulating ABA levels. These results strongly suggest the central role of ABA pathway in the control of dormancy progression and open up new perspectives for the development of molecular-based phenological modelling.

The Anthropocene epoch is associated with the spreading of metals in the environment increasing oxidative and genotoxic stress on living organisms1,2. Once regarded as a curiosity, plants hyperaccumulating metals are now envisioned as an opportunity to remediate metal contaminated soils. About 500 plant species adapted to metalliferous soils acquired the capacity to hyperaccumulate (>0.1% of dry weight) nickel in their shoot3. The phylogenetic distribution of these hyperaccumulators in 50 families suggest that this complex trait evolved multiple times independently from basic mechanisms involved in metal homeostasis. However, the exact nature of these mechanisms and whether they are shared between various lineages is not known. Here, using cross-species transcriptomic analyses in different plant families, we have identified convergent functions that may represent major nodes in the evolution of nickel hyperaccumulation. In particular, our data point out that constitutive high expression of IREG/Ferroportin transporters recurrently emerged as a mechanism involved in nickel hyperaccumulation.

Downy mildews are obligate biotrophic oomycete pathogens that cause devastating plant diseases on economically important crops. Plasmopara viticola is the causal agent of grapevine downy mildew, a major disease in vineyards worldwide. We sequenced the genome of Pl. viticola with PacBio long reads and obtained a new 92.94 Mb assembly with high continuity (359 scaffolds for a N50 of 706.5 kb) due to a better resolution of repeat regions. This assembly presented a high level of gene completeness, recovering 1,592 genes encoding secreted proteins involved in plant-pathogen interactions. Pl. viticola had a two-speed genome architecture, with secreted protein-encoding genes preferentially located in gene-sparse, repeat-rich regions and evolving rapidly, as indicated by pairwise dN/dS values. We also used short reads to assemble the genome of Plasmopara muralis, a closely related species infecting grape ivy (Parthenocissus tricuspidata). The lineage-specific proteins identified by comparative genomics analysis included a large proportion of RxLR cytoplasmic effectors and, more generally, genes with high dN/dS values. We identified 270 candidate genes under positive selection, including several genes encoding transporters and components of the RNA machinery potentially involved in host specialization. Finally, the Pl. viticola genome assembly generated here will allow the development of robust population genomics approaches for investigating the mechanisms involved in adaptation to biotic and abiotic selective pressures in this species.DATA AVAILABILITYRaw reads and genome assemblies have been deposited in GenBank (BioProjects PRJNA329579 for Pl. viticola and PRJNA448661 for Pl. muralis). Genome assemblies, gene annotations and analysis files (e.g. orthology relationships, full tables for GO enrichment analyses, pairwise dN/dS values and branch-site tests) have been deposited in Dataverse (Pl. viticola assembly and annotation: doi.org/10.15454/4NYHD6, Pl. muralis assembly and annotation: doi.org/10.15454/Q1QJYK, analysis files: doi.org/10.15454/8NZ8X9). Links to the data and information about the grapevine downy mildew genome project can be found at http://grapevine-downy-mildew-genome.com/.

Many animal species are comprised of discrete phenotypic forms. Understanding the genetic mechanisms generating and maintaining such phenotypic variation within species is essential to comprehending morphological diversity. A common and conspicuous example of discrete phenotypic variation in natural populations of insects is the occurrence of different colour patterns, which has motivated a rich body of ecological and genetic research1–6. The occurrence of dark, i.e. melanic, forms, displaying discrete colour patterns, is found across multiple taxa, but the underlying genomic basis remains poorly characterized. In numerous ladybird species (Coccinellidae), the spatial arrangement of black and orange patches on adult elytra varies wildly within species, forming strikingly different complex colour patterns7,8. In the harlequin ladybirdHarmonia axyridis, more than 200 distinct colour forms have been described, which classic genetic studies suggest result from allelic variation at a single, unknown, locus9,10. Here, we combined whole-genome sequencing, population genomics, gene expression and functional analyses, to establish that the genepanniercontrols melanic pattern polymorphism inH. axyridis. We show thatpannier, which encodes an evolutionary conserved transcription factor, is necessary for the formation of melanic elements on the elytra. Allelic variation inpannierleads to protein expression in distinct domains on the elytra, and thus determines the distinct colour patterns inH. axyridis. Recombination betweenpannieralleles may be reduced by a highly divergent sequence of ca. 170 kb in thecis-regulatory regions ofpannierwith a 50 kb inversion between colour forms. This likely helps maintaining the distinct alleles found in natural populations. Thus we propose that highly variable discrete colour forms can arise in natural populations throughcis-regulatory allelic variation of a single gene.

Identifying genetic variation that increases crop yields is a primary objective in plant breeding. We have used association analyses of Brassica napus (oilseed rape/canola) accessions to identify variation in the expression of a HECT E3 ligase gene, BnaUPL3.C03, that influences seed size and final yield. We establish a mechanism in which UPL3 mediates the proteasomal degradation of LEC2, a master regulator of seed maturation. Reduced UPL3 expression increases LEC2 protein levels and prolongs expression of lipid biosynthetic genes and seed maturation. Natural variation in BnaUPL3.C03 expression has not yet been exploited in current Brassica napus breeding lines and can therefore be used as a new approach to maximize future yields in this important oil crop.

The bacterium Xylella fastidiosa (Xf) is a plant endophyte native to the Americas that causes worldwide concern. Xf has been recently detected in several regions outside its natural range including Europe. In that context, accurate estimates of its response to climate change are required to design cost-efficient and environment-friendly control strategies. In the present study, we collected data documenting the native and invaded ranges of the three main subspecies of Xf: fastidiosa, pauca and multiplex, as well as two strains of Xf subsp. multiplex recently detected in southern France (ST6 and ST7). We fitted bioclimatic species distribution models (SDMs) to forecast their potential geographic range and impact in Europe under current and future climate conditions. According to model predictions, the geographical range of Xf as presently reported in Europe is small compared to the large extent of suitable areas. The European regions most threatened by Xf encompass the Mediterranean coastal areas of Spain, Greece, Italy and France, the Atlantic coastal areas of France, Portugal and Spain as well as the south-western regions of Spain and lowlands in southern Italy. Potential distribution of the different subspecies / strains are contrasted but all are predicted to increase by 2050, which could threaten several of the most economically important wine-, olive- and fruit-producing regions of Europe, warranting the design of control strategies. Bioclimatic models also predict that the subspecies multiplex might represent a threat to most of Europe under current and future climate conditions. These results may serve as a basis for future design of a spatially informed European-scale integrated management strategy, including early detection surveys in plants and insect vectors, quarantine measures as well as agricultural practices.

The range of hosts that a parasite can infect in nature is a trait determined by its own evolutionary history and that of its potential hosts. However, knowledge on host range diversity and evolution at the family level is often lacking. Here, we investigate host range variation and diversification trends within theSclerotiniaceae, a family of Ascomycete fungi. Using a phylogenetic framework, we associate diversification rates, the frequency of host jump events, and host range variation during the evolution of this family. Variations in diversification rate during the evolution of the Sclerotiniaceae define three major macro-evolutionary regimes with contrasted proportions of species infecting a broad range of hosts. Host-parasite co-phylogenetic analyses pointed towards parasite radiation on distant hosts long after host speciation (host jump or duplication events) as the dominant mode of association with plants in theSclerotiniaceae. The intermediate macro-evolutionary regime showed a low diversification rate, high frequency of duplication events, and the highest proportion of broad host range species. Consistent with previous reports on oomycete parasites, our findings suggest that host jump and radiation, possibly combined with low speciation rates, could associate with the emergence of generalist pathogens. These results have important implications for our understanding of fungal parasites evolution and are of particular relevance for the durable management of disease epidemics.

This preprint has been reviewed and recommended by Peer Community In Evolutionary Biology (http://dx.doi.org/10.24072/pci.evolbiol.100041).The mating system of a species is expected to have important effects on its genetic diversity. In this paper, we explore the effects of partial selfing on the equilibrium genetic varianceVg, mutation loadLand inbreeding depressionδunder stabilizing selection acting on a arbitrary numbernof quantitative traits coded by biallelic loci with additive effects. Overall, our model predicts a decrease in the equilibrium genetic variance with increasing selfing rates; however, the relationship between self-fertilization and the variables of interest depends on the strength of associations between loci, and three different regimes are observed. When theU/nratio is low (whereUis the total haploid mutation rate on selected traits) and effective recombination rates are sufficiently high, genetic associations between loci are negligible and the genetic variance, mutation load and inbreeding depression are well predicted by approximations based on single-locus models. For higher values ofU/nand/or lower effective recombination, moderate genetic associations generated by epistasis tend to increaseVg,Landδ, this regime being well predicted by approximations including the effects of pairwise associations between loci. For yet higher values ofU/nand/or lower effective recombination, a different regime is reached under which the maintenance of coadapted gene complexes reducesVg,Landδ. Simulations indicate that the values ofVg,Landδare little affected by assumptions regarding the number of possible alleles per locus.

Heterosis, the superiority of hybrids over their parents for quantitative traits, represents a crucial issue in plant and animal breeding. Heterosis has given rise to countless genetic, genomic and molecular studies, but has rarely been investigated from the point of view of systems biology. We hypothesized that heterosis is an emergent property of living systems resulting from frequent concave relationships between genotypic variables and phenotypes, or between different phenotypic levels. We chose the enzyme-flux relationship as a model of the concave genotype-phenotype (GP) relationship, and showed that heterosis can be easily created in the laboratory. First, we reconstitutedin vitrothe upper part of glycolysis. We simulated genetic variability of enzyme activity by varying enzyme concentrations in test tubes. Mixing the content of “parental” tubes resulted in “hybrids”, whose fluxes were compared to the parental fluxes. Frequent heterotic fluxes were observed, under conditions that were determined analytically and confirmed by computer simulation. Second, to test this model in a more realistic situation, we modeled the glycolysis/fermentation network in yeast by considering one input flux, glucose, and two output fluxes, glycerol and acetaldehyde. We simulated genetic variability by randomly drawing parental enzyme concentrations under various conditions, and computed the parental and hybrid fluxes using a system of differential equations. Again we found that a majority of hybrids exhibited positive heterosis for metabolic fluxes. Cases of negative heterosis were due to local convexity between certain enzyme concentrations and fluxes. In both approaches, heterosis was maximized when the parents were phenotypically closeandwhen the distributions of parental enzyme concentrations were contrasted and constrained. These conclusions are not restricted to metabolic systems: they only depend on the concavity of the GP relationship, which is commonly observed at various levels of the phenotypic hierarchy, and could account for the pervasiveness of heterosis.