Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns.

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

The giant diploid faba genome unlocks variation in a global protein crop.

Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

A Solanum lycopersicoides reference genome facilitates insights into tomato specialized metabolism and immunity.

Wild relatives of tomato are a valuable source of natural variation in tomato breeding, as many can be hybridized to the cultivated species (Solanum lycopersicum). Several, including Solanum lycopersicoides, have been crossed to S. lycopersicum for the development of ordered introgression lines (ILs), facilitating breeding for desirable traits. Despite the utility of these wild relatives and their associated ILs, few finished genome sequences have been produced to aid genetic and genomic studies. Here we report a chromosome-scale genome assembly for S. lycopersicoides LA2951, which contains 37 938 predicted protein-coding genes. With the aid of this genome assembly, we have precisely delimited the boundaries of the S. lycopersicoides introgressions in a set of S. lycopersicum cv. VF36 × LA2951 ILs. We demonstrate the usefulness of the LA2951 genome by identifying several quantitative trait loci for phenolics and carotenoids, including underlying candidate genes, and by investigating the genome organization and immunity-associated function of the clustered Pto gene family. In addition, syntenic analysis of R2R3MYB genes sheds light on the identity of the Aubergine locus underlying anthocyanin production. The genome sequence and IL map provide valuable resources for studying fruit nutrient/quality traits, pathogen resistance, and environmental stress tolerance. We present a new genome resource for the wild species S. lycopersicoides, which we use to shed light on the Aubergine locus responsible for anthocyanin production. We also provide IL boundary mappings, which facilitated identifying novel carotenoid quantitative trait loci of which one was likely driven by an uncharacterized lycopene ß-cyclase whose function we demonstrate.

GXP: Analyze and Plot Plant Omics Data in Web Browsers.

Next-generation sequencing and metabolomics have become very cost and work efficient and are integrated into an ever-growing number of life science research projects. Typically, established software pipelines analyze raw data and produce quantitative data informing about gene expression or concentrations of metabolites. These results need to be visualized and further analyzed in order to support scientific hypothesis building and identification of underlying biological patterns. Some of these tools already exist, but require installation or manual programming. We developed "Gene Expression Plotter" (GXP), an RNAseq and Metabolomics data visualization and analysis tool entirely running in the user's web browser, thus not needing any custom installation, manual programming or uploading of confidential data to third party servers. Consequently, upon receiving the bioinformatic raw data analysis of RNAseq or other omics results, GXP immediately enables the user to interact with the data according to biological questions by performing knowledge-driven, in-depth data analyses and candidate identification via visualization and data exploration. Thereby, GXP can support and accelerate complex interdisciplinary omics projects and downstream analyses. GXP offers an easy way to publish data, plots, and analysis results either as a simple exported file or as a custom website. GXP is freely available on GitHub (see introduction).

Crop genetic diversity uncovers metabolites, elements, and gene networks predicted to be associated with high plant biomass yields in maize.

Rapid population growth and increasing demand for food, feed, and bioenergy in these times of unprecedented climate change require breeding for increased biomass production on the world's croplands. To accelerate breeding programs, knowledge of the relationship between biomass features and underlying gene networks is needed to guide future breeding efforts. To this end, large-scale multiomics datasets were created with genetically diverse maize lines, all grown in long-term organic and conventional cropping systems. Analysis of the datasets, integrated using regression modeling and network analysis revealed key metabolites, elements, gene transcripts, and gene networks, whose contents during vegetative growth substantially influence the build-up of plant biomass in the reproductive phase. We found that S and P content in the source leaf and P content in the root during the vegetative stage contributed the most to predicting plant performance at the reproductive stage. In agreement with the Gene Ontology enrichment analysis, the cis-motifs and identified transcription factors associated with upregulated genes under phosphate deficiency showed great diversity in the molecular response to phosphate deficiency in selected lines. Furthermore, our data demonstrate that genotype-dependent uptake, assimilation, and allocation of essential nutrient elements (especially C and N) during vegetative growth under phosphate starvation plays an important role in determining plant biomass by controlling root traits related to nutrient uptake. These integrative multiomics results revealed key factors underlying maize productivity and open new opportunities for efficient, rapid, and cost-effective plant breeding to increase biomass yield of the cereal crop maize under adverse environmental factors.

Plant PhysioSpace: a robust tool to compare stress response across plant species.

Generalization of transcriptomics results can be achieved by comparison across experiments. This generalization is based on integration of interrelated transcriptomics studies into a compendium. Such a focus on the bigger picture enables both characterizations of the fate of an organism and distinction between generic and specific responses. Numerous methods for analyzing transcriptomics datasets exist. Yet, most of these methods focus on gene-wise dimension reduction to obtain marker genes and gene sets for, for example, pathway analysis. Relying only on isolated biological modules might result in missing important confounders and relevant contexts. We developed a method called Plant PhysioSpace, which enables researchers to compute experimental conditions across species and platforms without a priori reducing the reference information to specific gene sets. Plant PhysioSpace extracts physiologically relevant signatures from a reference dataset (i.e. a collection of public datasets) by integrating and transforming heterogeneous reference gene expression data into a set of physiology-specific patterns. New experimental data can be mapped to these patterns, resulting in similarity scores between the acquired data and the extracted compendium. Because of its robustness against platform bias and noise, Plant PhysioSpace can function as an inter-species or cross-platform similarity measure. We have demonstrated its success in translating stress responses between different species and platforms, including single-cell technologies. We have also implemented two R packages, one software and one data package, and a Shiny web application to facilitate access to our method and precomputed models.

Photoreceptor Activity Contributes to Contrasting Responses to Shade in Cardamine and Arabidopsis Seedlings.

Plants have evolved two major ways to deal with nearby vegetation or shade: avoidance and tolerance. Moreover, some plants respond to shade in different ways; for example, Arabidopsis (Arabidopsis thaliana) undergoes an avoidance response to shade produced by vegetation, but its close relative Cardamine hirsuta tolerates shade. How plants adopt opposite strategies to respond to the same environmental challenge is unknown. Here, using a genetic strategy, we identified the C. hirsuta slender in shade1 mutants, which produce strongly elongated hypocotyls in response to shade. These mutants lack the phytochrome A (phyA) photoreceptor. Our findings suggest that C. hirsuta has evolved a highly efficient phyA-dependent pathway that suppresses hypocotyl elongation when challenged by shade from nearby vegetation. This suppression relies, at least in part, on stronger phyA activity in C. hirsuta; this is achieved by increased ChPHYA expression and protein accumulation combined with a stronger specific intrinsic repressor activity. We suggest that modulation of photoreceptor activity is a powerful mechanism in nature to achieve physiological variation (shade tolerance versus avoidance) for species to colonize different habitats.

Cell type specific transcriptional reprogramming of maize leaves during Ustilago maydis induced tumor formation.

Ustilago maydis is a biotrophic pathogen and well-established genetic model to understand the molecular basis of biotrophic interactions. U. maydis suppresses plant defense and induces tumors on all aerial parts of its host plant maize. In a previous study we found that U. maydis induced leaf tumor formation builds on two major processes: the induction of hypertrophy in the mesophyll and the induction of cell division (hyperplasia) in the bundle sheath. In this study we analyzed the cell-type specific transcriptome of maize leaves 4 days post infection. This analysis allowed identification of key features underlying the hypertrophic and hyperplasic cell identities derived from mesophyll and bundle sheath cells, respectively. We examined the differentially expressed (DE) genes with particular focus on maize cell cycle genes and found that three A-type cyclins, one B-, D- and T-type are upregulated in the hyperplasic tumorous cells, in which the U. maydis effector protein See1 promotes cell division. Additionally, most of the proteins involved in the formation of the pre-replication complex (pre-RC, that assure that each daughter cell receives identic DNA copies), the transcription factors E2F and DPa as well as several D-type cyclins are deregulated in the hypertrophic cells.

Recently duplicated sesterterpene (C25) gene clusters in Arabidopsis thaliana modulate root microbiota.

Land plants co-speciate with a diversity of continually expanding plant specialized metabolites (PSMs) and root microbial communities (microbiota). Homeostatic interactions between plants and root microbiota are essential for plant survival in natural environments. A growing appreciation of microbiota for plant health is fuelling rapid advances in genetic mechanisms of controlling microbiota by host plants. PSMs have long been proposed to mediate plant and single microbe interactions. However, the effects of PSMs, especially those evolutionarily new PSMs, on root microbiota at community level remain to be elucidated. Here, we discovered sesterterpenes in Arabidopsis thaliana, produced by recently duplicated prenyltransferase-terpene synthase (PT-TPS) gene clusters, with neo-functionalization. A single-residue substitution played a critical role in the acquisition of sesterterpene synthase (sesterTPS) activity in Brassicaceae plants. Moreover, we found that the absence of two root-specific sesterterpenoids, with similar chemical structure, significantly affected root microbiota assembly in similar patterns. Our results not only demonstrate the sensitivity of plant microbiota to PSMs but also establish a complete framework of host plants to control root microbiota composition through evolutionarily dynamic PSMs.

MapMan4: A Refined Protein Classification and Annotation Framework Applicable to Multi-Omics Data Analysis.

Genome sequences from over 200 plant species have already been published, with this number expected to increase rapidly due to advances in sequencing technologies. Once a new genome has been assembled and the genes identified, the functional annotation of their putative translational products, proteins, using ontologies is of key importance as it places the sequencing data in a biological context. Furthermore, to keep pace with rapid production of genome sequences, this functional annotation process must be fully automated. Here we present a redesigned and significantly enhanced MapMan4 framework, together with a revised version of the associated online Mercator annotation tool. Compared with the original MapMan, the new ontology has been expanded almost threefold and enforces stricter assignment rules. This framework was then incorporated into Mercator4, which has been upgraded to reflect current knowledge across the land plant group, providing protein annotations for all embryophytes with a comparably high quality. The annotation process has been optimized to allow a plant genome to be annotated in a matter of minutes. The output results continue to be compatible with the established MapMan desktop application.

Heterologous Expression, Biosynthetic Studies, and Ecological Function of the Selective Gq-Signaling Inhibitor FR900359.

The cyclic depsipeptide FR900359 (FR), isolated from the tropical plant Ardisia crenata, is a strong and selective inhibitor of Gq proteins, making it an indispensable pharmacological tool to study Gq-related processes, as well as a promising drug candidate. Gq inhibition is a novel mode of action for defense chemicals and crucial for the ecological function of FR, as shown by in vivo experiments in mice, its affinity to insect Gq proteins, and insect toxicity studies. The uncultured endosymbiont of A. crenata was sequenced, revealing the FR nonribosomal peptide synthetase (frs) gene cluster. We here provide a detailed model of FR biosynthesis, supported by in vitro enzymatic and bioinformatic studies, and the novel analogue AC-1, which demonstrates the flexibility of the FR starter condensation domains. Finally, expression of the frs genes in E. coli led to heterologous FR production in a cultivable, bacterial host for the first time.

The Cardamine hirsuta genome offers insight into the evolution of morphological diversity.

Finding causal relationships between genotypic and phenotypic variation is a key focus of evolutionary biology, human genetics and plant breeding. To identify genome-wide patterns underlying trait diversity, we assembled a high-quality reference genome of Cardamine hirsuta, a close relative of the model plant Arabidopsis thaliana. We combined comparative genome and transcriptome analyses with the experimental tools available in C. hirsuta to investigate gene function and phenotypic diversification. Our findings highlight the prevalent role of transcription factors and tandem gene duplications in morphological evolution. We identified a specific role for the transcriptional regulators PLETHORA5/7 in shaping leaf diversity and link tandem gene duplication with differential gene expression in the explosive seed pod of C. hirsuta. Our work highlights the value of comparative approaches in genetically tractable species to understand the genetic basis for evolutionary change.

Structural basis of head to head polyketide fusion by CorB.

Corallopyronin A is a polyketide derived from the myxobacterium Corallococcus coralloides with potent antibiotic features. The gene cluster responsible for the biosynthesis of corallopyronin A has been described recently, and it was proposed that CorB acts as a ketosynthase to interconnect two polyketide chains in a rare head-to-head condensation reaction. We determined the structure of CorB, the interconnecting polyketide synthase, to high resolution and found that CorB displays a thiolase fold. Site-directed mutagenesis showed that the catalytic triad consisting of a cysteine, a histidine and an asparagine is crucial for catalysis, and that this triad shares similarities with the triad found in HMG-CoA synthases. We synthesized a substrate mimic to derivatize purified CorB and confirmed substrate attachment by ESI-MS. Structural analysis of the complex yielded an electron density-based model for the polyketide chain and showed that the unusually wide, T-shaped active site is able to accommodate two polyketides simultaneously. Our structural analysis provides a platform for understanding the unusual head-to-head polyketide-interconnecting reaction catalyzed by CorB.