Evolutionary dynamics of an Arabidopsis insect resistance quantitative trait locus

Evolutionary Systems Biology Identifies Genetic Trade-offs in Rice Defense against Aboveground and Belowground Attackers

Like other plants, wild and domesticated rice species (Oryza nivara, O. rufipogon, and O. sativa) evolve in environments with various biotic and abiotic stresses that fluctuate in intensity through space and time. Microbial pathogens and invertebrate herbivores such as plant-parasitic nematodes and caterpillars show geographical and temporal variation in activity patterns and may respond differently to certain plant-defensive mechanisms. As such, plant interactions with multiple community members may result in conflicting selection pressures on genetic polymorphisms. Here, through assays with different aboveground and belowground herbivores, the fall armyworm (Spodoptera frugiperda) and the southern root-knot nematode (Meloidogyne incognita), and comparison with rice responses to microbial pathogens, we identify potential genetic trade-offs at the KSL8 and MG1 loci on chromosome 11. KSL8 encodes the first committed step toward the biosynthesis of either stemarane-type or stemodane-type diterpenoids through the japonica (KSL8-jap) or indica (KSL8-ind) allele. Knocking out KSL8-jap and CPS4, encoding an enzyme that acts upstream in diterpenoid synthesis, in japonica rice cultivars increased resistance to S. frugiperda and decreased resistance to M. incognita. Furthermore, MG1 resides in a haplotype that provided resistance to M. incognita, while alternative haplotypes are involved in mediating resistance to the rice blast fungus Magnaporthe oryzae and other pests and pathogens. Finally, KSL8 and MG1 alleles are located within trans-species polymorphic haplotypes and may be evolving under long-term balancing selection. Our data are consistent with a hypothesis that polymorphisms at KSL8 and MG1 may be maintained through complex and diffuse community interactions.

Ectopic Gene Conversion Causing Quantitative Trait Variation

Why is there so much non-neutral genetic variation segregating in natural populations? We dissect function and evolution of a near-cryptic quantitative trait locus (QTL) for defense metabolites in Arabidopsis using the CRISPR/Cas9 system and nucleotide polymorphism patterns. The QTL is explained by genetic variation in a family of 4 tightly linked indole-glucosinolate O-methyltransferase genes. Some of this variation appears to be maintained by balancing selection, some appears to be generated by non-reciprocal transfer of sequence, also known as ectopic gene conversion (EGC), between functionally diverged gene copies. Here, we elucidate how EGC, as an inevitable consequence of gene duplication, could be a general mechanism for generating genetic variation for fitness traits.

Glucosinolate profile and specifier protein activity determine the glucosinolate hydrolysis product formation in kohlrabi (Brassica oleracea var. gongylodes) in a tissue-specific way

Glucosinolates, commonly found in Brassica vegetables, are hydrolyzed by myrosinase to form bioactive isothiocyanates, unless specifier proteins redirect the degradation to less bioactive nitriles and epithionitriles. Here, the tissue-specific impact of specifier proteins on the outcome of glucosinolate hydrolysis in nine kohlrabi tissues was investigated. Glucosinolates and their hydrolysis product profiles, epithiospecifier protein and myrosinase activity, and protein abundance patterns of key glucosinolate biosynthesis, transport and hydrolysis enzymes were determined and correlated to the metabolites in the kohlrabi tissues. Leaf tissues were rich in glucoraphanin, while bulb tissues contained more glucoerucin. Moreover, a higher proportion of isothiocyanates was formed in leaf stalk, bulb peel, stem and root, tissues with relatively higher epithiospecifier modifier 1 abundance and a higher ratio of myrosinase activity to ESP activity. This study shows that tissue-specific glucosinolate hydrolysis is mediated by glucosinolate profiles, myrosinase activity and specifier protein as well as modifier protein abundance.

Allopolyploidy expanded gene content but not pangenomic variation in the hexaploid oilseed Camelina sativa

Ancient whole-genome duplications are believed to facilitate novelty and adaptation by providing the raw fuel for new genes. However, it is unclear how recent whole-genome duplications may contribute to evolvability within recent polyploids. Hybridization accompanying some whole-genome duplications may combine divergent gene content among diploid species. Some theory and evidence suggest that polyploids have a greater accumulation and tolerance of gene presence-absence and genomic structural variation, but it is unclear to what extent either is true. To test how recent polyploidy may influence pangenomic variation, we sequenced, assembled, and annotated 12 complete, chromosome-scale genomes of Camelina sativa, an allohexaploid biofuel crop with 3 distinct subgenomes. Using pangenomic comparative analyses, we characterized gene presence-absence and genomic structural variation both within and between the subgenomes. We found over 75% of ortholog gene clusters are core in C. sativa and <10% of sequence space was affected by genomic structural rearrangements. In contrast, 19% of gene clusters were unique to one subgenome, and the majority of these were Camelina specific (no ortholog in Arabidopsis). We identified an inversion that may contribute to vernalization requirements in winter-type Camelina and an enrichment of Camelina-specific genes with enzymatic processes related to seed oil quality and Camelina's unique glucosinolate profile. Genes related to these traits exhibited little presence-absence variation. Our results reveal minimal pangenomic variation in this species and instead show how hybridization accompanied by whole-genome duplication may benefit polyploids by merging diverged gene content of different species.

Unraveling the specialized metabolic pathways in medicinal plant genomes: a review

Medicinal plants are important sources of bioactive specialized metabolites with significant therapeutic potential. Advances in multi-omics have accelerated the understanding of specialized metabolite biosynthesis and regulation. Genomics, transcriptomics, proteomics, and metabolomics have each contributed new insights into biosynthetic gene clusters (BGCs), metabolic pathways, and stress responses. However, single-omics approaches often fail to fully address these complex processes. Integrated multi-omics provides a holistic perspective on key regulatory networks. High-throughput sequencing and emerging technologies like single-cell and spatial omics have deepened our understanding of cell-specific and spatially resolved biosynthetic dynamics. Despite these advancements, challenges remain in managing large datasets, standardizing protocols, accounting for the dynamic nature of specialized metabolism, and effectively applying synthetic biology for sustainable specialized metabolite production. This review highlights recent progress in omics-based research on medicinal plants, discusses available bioinformatics tools, and explores future research trends aimed at leveraging integrated multi-omics to improve the medicinal quality and sustainable utilization of plant resources.

Integrated Analysis of Metabolome and Transcriptome Reveals the Effect of Burdock Fructooligosaccharide on the Quality of Chinese Cabbage (Brassica rapa L. ssp. Pekinensis)

Burdock fructooligosaccharide (BFO) is fructose with a low polymerization degree, which could improve the immunity to pathogens, quality, and stress resistance of vegetables. Still, there are no studies on applying BFO in Chinese cabbage. In this study, the effects of exogenous BFO sprayed with different concentrations (0, 5, 10, 20, 30 g·L-1) on the growth and soluble sugar content of Chinese cabbage seedlings were determined. The result showed that 10 g·L-1 was the appropriate spraying concentration. Based on metabolome analysis, a total of 220 differentially accumulated metabolites (DAMs) were found, among which flavonoid metabolites, glucosinolate metabolites, and soluble sugar-related metabolites were the key metabolites involved in improving the quality of Chinese cabbage caused by BFO. Further combination analysis with transcriptome, trans-cinnamate 4-monooxygenase (CYP73A5), and chalcone synthase 1 (CHS1) were more closely associated with the DAMs of flavonoid biosynthesis. Sulfotransferases 18 (SOT18), Branched-chain amino acid amino transferases 6 (BCAT6), and cytochrome P450 monooxygenase (CYP83A1) were the key genes in glucosinolate biosynthesis. Hexokinase (HxK1), beta-glucosidase 8 (BGL08), invertase 3 (INV3), beta-glucosidase 3B (BGL3B), and sucrose phosphate synthase 1 (SPS1) were significantly upregulated, potentially playing crucial roles in the soluble sugar metabolism. In conclusion, these results provided an understanding of the effects of BFO on the expression of genes and the accumulation of metabolites related to quality formation in Chinese cabbage.

Bivariate GWA mapping reveals associations between aliphatic glucosinolates and plant responses to thrips and heat stress

Although plants harbor a huge phytochemical diversity, only a fraction of plant metabolites is functionally characterized. In this work, we aimed to identify the genetic basis of metabolite functions during harsh environmental conditions in Arabidopsis thaliana. With machine learning algorithms we predicted stress-specific metabolomes for 23 (a)biotic stress phenotypes of 300 natural Arabidopsis accessions. The prediction models identified several aliphatic glucosinolates (GLSs) and their breakdown products to be implicated in responses to heat stress in siliques and herbivory by Western flower thrips, Frankliniella occidentalis. Bivariate GWA mapping of the metabolome predictions and their respective (a)biotic stress phenotype revealed genetic associations with MAM, AOP, and GS-OH, all three involved in aliphatic GSL biosynthesis. We, therefore, investigated thrips herbivory on AOP, MAM, and GS-OH loss-of-function and/or overexpression lines. Arabidopsis accessions with a combination of MAM2 and AOP3, leading to 3-hydroxypropyl dominance, suffered less from thrips feeding damage. The requirement of MAM2 for this effect could, however, not be confirmed with an introgression line of ecotypes Cvi and Ler, most likely due to other, unknown susceptibility factors in the Ler background. However, AOP2 and GS-OH, adding alkenyl or hydroxy-butenyl groups, respectively, did not have major effects on thrips feeding. Overall, this study illustrates the complex implications of aliphatic GSL diversity in plant responses to heat stress and a cell-content-feeding herbivore.

Facilitation of Sclerotinia sclerotiorum infestation by aphid feeding behaviour is not affected by aphid resistance in oilseed rape

The relation between aphids and Sclerotinia stem rot (SSR) in oilseed rape is rarely examined because they are often studied alone. We have observed a significant correlation between the number of aphids and the occurrence of SSR in our field studies. Electropenetrography (EPG) was used to evaluate the effects of Brevicoryne brassicae (Linnaeus) on two oilseed rape cultivars while acquiring, vectoring and inoculating of Sclerotinia sclerotiorum Lib. (de Bary) ascospores. The results demonstrated that aphid feeding followed by the application of an ascospore suspension significantly increased S. sclerotiorum incidence. Aphids were capable of adhering to ascospores and carrying them to healthy plants, thereby causing diseases. The results of the EPG analysis indicated that aphid feeding behaviour was significantly altered in all leaf tissue levels following infection with S. sclerotiorum. Aphids initiated their first puncture significantly sooner than the control group, began probing mesophyll cells earlier, significantly increased the frequency of both short probes and intracellular punctures and had a significantly shorter pathway duration. On infected aphid-susceptible cultivars, aphids secreted more saliva but had reduced ingestion compared with aphids feeding on non-infected oilseed rape. In addition, ascospores can affect aphid feeding behaviour by adhering to aphids. Aphids carrying ascospores punctured cells earlier, with a significant increase in the frequency and duration of short probes and cell punctures, shortened pathway durations, increased salivation and reduced ingestion compared with aphids not carrying ascospores. On aphid-susceptible cultivars, aphids carrying ascospores delayed puncture onset, but on resistant cultivars, puncture onset was shortened. There is a correlation between aphids and S. sclerotiorum. The impact of S. sclerotiorum on aphid feeding behaviour is directional, favouring the spread of the fungus. This promotion does not appear to be altered by the aphid resistance of the cultivar.

Chromosome-scale reference genome of broccoli (Brassica oleracea var. italica Plenck) provides insights into glucosinolate biosynthesis

Broccoli (Brassica oleracea var. italica Plenck) is an important vegetable crop, as it is rich in health-beneficial glucosinolates (GSLs). However, the genetic basis of the GSL diversity in Brassicaceae remains unclear. Here we report a chromosome-level genome assembly of broccoli generated using PacBio HiFi reads and Hi-C technology. The final genome assembly is 613.79 Mb in size, with a contig N50 of 14.70 Mb. The GSL profile and content analysis of different B. oleracea varieties, combined with a phylogenetic tree analysis, sequence alignment, and the construction of a 3D model of the methylthioalkylmalate synthase 1 (MAM1) protein, revealed that the gene copy number and amino acid sequence variation both contributed to the diversity of GSL biosynthesis in B. oleracea. The overexpression of BoMAM1 (BolI0108790) in broccoli resulted in high accumulation and a high ratio of C4-GSLs, demonstrating that BoMAM1 is the key enzyme in C4-GSL biosynthesis. These results provide valuable insights for future genetic studies and nutritive component applications of Brassica crops.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

Identification and evolution of a diterpenoid phytoalexin oryzalactone biosynthetic gene in the genus Oryza

The natural variation of plant-specialized metabolites represents the evolutionary adaptation of plants to their environments. However, the molecular mechanisms that account for the diversification of the metabolic pathways have not been fully clarified. Rice plants resist attacks from pathogens by accumulating diterpenoid phytoalexins. It has been confirmed that the composition of rice phytoalexins exhibits numerous natural variations. Major rice phytoalexins (momilactones and phytocassanes) are accumulated in most cultivars, although oryzalactone is a cultivar-specific compound. Here, we attempted to reveal the evolutionary trajectory of the diversification of phytoalexins by analyzing the oryzalactone biosynthetic gene in Oryza species. The candidate gene, KSLX-OL, which accounts for oryzalactone biosynthesis, was found around the single-nucleotide polymorphisms specific to the oryzalactone-accumulating cultivars in the long arm of chromosome 11. The metabolite analyses in Nicotiana benthamiana and rice plants overexpressing KSLX-OL indicated that KSLX-OL is responsible for the oryzalactone biosynthesis. KSLX-OL is an allele of KSL8 that is involved in the biosynthesis of another diterpenoid phytoalexin, oryzalexin S and is specifically distributed in the AA genome species. KSLX-NOL and KSLX-bar, which encode similar enzymes but are not involved in oryzalactone biosynthesis, were also found in AA genome species. The phylogenetic analyses of KSLXs, KSL8s, and related pseudogenes (KSL9s) indicated that KSLX-OL was generated from a common ancestor with KSL8 and KSL9 via gene duplication, functional differentiation, and gene fusion. The wide distributions of KSLX-OL and KSL8 in AA genome species demonstrate their long-term coexistence beyond species differentiation, suggesting a balancing selection between the genes.

Natural genetic variation and negative density effects in plant-nematode interactions

Arabidopsis thaliana is a suitable host for phytoparasitic nematodes of the genus Meloidogyne. Successful nematode infection leads to the formation of root galls. We tested for natural genetic variation and inoculation density effects on nematode reproductive success in the interaction between A. thaliana and Meloidogyne javanica. We inoculated different Arabidopsis genotypes with two sources of nematodes at two different doses, using a mild protocol for inoculum preparation. We counted root galls and egg masses 2 months after inoculation. We obtained a high number of successful nematode infections. Infection success differed among Arabidopsis genotypes in interaction with the nematode source. Overall, infection success and reproductive success of nematodes were lower at a higher inoculum dose of nematodes. Our results indicate that natural genetic variation in both host plants and nematodes, as well as short- and long-term negative density effects, shape nematode reproductive success.

Large-scale identification of novel transcriptional regulators of the aliphatic glucosinolate pathway in Arabidopsis

Aliphatic glucosinolates are a large group of plant secondary metabolites characteristic of Brassicaceae, including the model plant Arabidopsis. The diverse and complex degradation products of aliphatic glucosinolates contribute to plant responses to herbivory, pathogen attack, and environmental stresses. Most of the biosynthesis genes in the aliphatic glucosinolate pathway have been cloned in Arabidopsis, and the research focus has recently shifted to the regulatory mechanisms controlling aliphatic glucosinolate accumulation. Up till now, more than 40 transcriptional regulators have been identified as regulating the aliphatic glucosinolate pathway, but many more novel regulators likely remain to be discovered based on research evidence over the past decade. In the current study, we took a systemic approach to functionally test 155 candidate transcription factors in Arabidopsis identified by yeast one-hybrid assay, and successfully validated at least 30 novel regulators that could significantly influence the accumulation of aliphatic glucosinolates in our experimental set-up. We also showed that the regulators of the aliphatic glucosinolate pathway have balanced positive and negative effects, and glucosinolate metabolism and plant development can be coordinated. Our work is the largest scale effort so far to validate transcriptional regulators of a plant secondary metabolism pathway, and provides new insights into how the highly diverse plant secondary metabolism is regulated at the transcriptional level.

Identification and in vitro enzymatic activity analysis of the AOP2 gene family associated with glucosinolate biosynthesis in Tumorous stem mustard (Brassica juncea var. tumida)

The major enzyme encoded by the glucosinolate biosynthetic gene AOP2 is involved in catalyzing the conversion of glucoiberin (GIB) into sinigrin (SIN) in Brassicaceae crops. The AOP2 proteins have previously been identified in several Brassicaceae species, but not in Tumorous stem mustard. As per this research, the five identified members of the AOP2 family from the whole genome of Brassica juncea named BjuAOP2.1-BjuAOP2.5 were found to be evenly distributed on five chromosomes. The subcellular localization results implied that BjuAOP2 proteins were mainly concentrated in the cytoplasm. Phylogenetic analysis of the AOP2 proteins from the sequenced Brassicaceae species in BRAD showed that BjuAOP2 genes were more closely linked to Brassica carinata and Brassica rapa than Arabidopsis. In comparison with other Brassicaceae plants, the BjuAOP2 members were conserved in terms of gene structures, protein sequences, and motifs. The light response and hormone response elements were included in the BjuAOP2 genes' cis-regulatory elements. The expression pattern of BjuAOP2 genes was influenced by the different stages of development and the type of tissue being examined. The BjuAOP2 proteins were used to perform the heterologous expression experiment. The results showed that all the five BjuAOP2 proteins can catalyze the conversion of GIB to SIN with different catalytic activity. These results provide the basis for further investigation of the functional study of BjuAOP2 in Tumorous stem mustard glucosinolate biosynthesis.

Natural variation and artificial selection at the BnaC2.MYB28 locus modulate Brassica napus seed glucosinolate

The degradation products of glucosinolates (GSLs) greatly lower the nutritional value of rapeseed (Brassica napus) meal; thus, reduction of seed GSL content (SGC) has become an important objective of rapeseed breeding. In our previous study, we finely mapped a major QTL (qGSL-C2) for SGC to a 49-kb collinear region on B. rapa chromosome A2. Here, we experimentally validated that BnaC2.MYB28, encoding an R2R3-MYB transcription factor, is the causal gene of qGSL-C2. BnaC2.MYB28 is a nucleus-localized protein mainly expressed in vegetative tissues. Knockout of BnaC2.MYB28 in the high-SGC parent G120 reduced SGC to a value lower than that in the low-SGC parent ZY50, while overexpression of BnaC2.MYB28 in both parental lines (G120 and ZY50) led to extremely high SGC, indicating that BnaC2.MYB28 acts as a positive regulator of SGC in both parents. Molecular characterization revealed that BnaC2.MYB28 forms a homodimer and specifically interacts with BnaMYC3. Moreover, BnaC2.MYB28 can directly activate the expression of GSL biosynthesis genes. Differential expression abundance resulting from the polymorphic promoter sequences, in combination with the different capability in activating downstream genes involved in aliphatic GSL biosynthesis, caused the functional divergence of BnaC2.MYB28 in SGC regulation between the parents. Natural variation of BnaC2.MYB28 was highly associated with SGC in natural germplasm and has undergone artificial selection in modern low-GSL breeding. This study provides important insights into the core function of BnaC2.MYB28 in regulating SGC and a promising strategy for manipulating SGC in rapeseed.

Transcriptome Analysis of Persian Oak (Quercus brantii L.) Decline Using RNA-seq Technology

Since the late 1980s, the oak decline has affected the Zagros oak forests in western Iran. Persian oak (Quercus brantii L.) the most important tree species of these forests has been damaged more than any other plant species. In the present study, the RNA sequencing technique was used for the first time to identify key genes and molecular mechanisms involved in Persian oak decline. The RNA was extracted from the leaves of healthy and declined oak trees, and sequenced using the Illumina HiSeq 2500 platform (2 × 150 bp paired-end reads). De novo transcriptome assembly of Persian oak revealed 56,743 unigenes and 6049 differentially expressed genes (DEGs) between declined and control samples. The results of gene ontology analysis showed that most of the DEGs involved in oak decline belong to the group of stress-responsive genes. In general, oak decline samples showed significant reductions in gene expression associated with "photosynthesis and storage of sugar" and "protein synthesis and related processes." Additionally, DEGs related to the starch degradation pathway were up-regulated, whereas DEGs associated with acetate-mevalonate (MVA), biosynthesis of lignin, and lignases pathways were down-regulated. The present study's findings can be an effective step in identifying the genes involved in oak decline and deciphering the relationship between this phenomenon and biotic and abiotic stresses.

Experimentally heat-induced transposition increases drought tolerance in Arabidopsis thaliana

Eukaryotic genomes contain a vast diversity of transposable elements (TEs). Formerly often described as selfish and parasitic DNA sequences, TEs are now recognised as a source of genetic diversity and powerful drivers of evolution. However, because their mobility is tightly controlled by the host, studies experimentally assessing how fast TEs may mediate the emergence of adaptive traits are scarce. We exposed Arabidopsis thaliana high-copy TE lines (hcLines) with up to c. eight-fold increased copy numbers of the heat-responsive ONSEN TE to drought as a straightforward and ecologically highly relevant selection pressure. We provide evidence for increased drought tolerance in five out of the 23 tested hcLines and further pinpoint one of the causative mutations to an exonic insertion of ONSEN in the ribose-5-phosphate-isomerase 2 gene. The resulting loss-of-function mutation caused a decreased rate of photosynthesis, plant size and water consumption. Overall, we show that the heat-induced transposition of a low-copy TE increases phenotypic diversity and leads to the emergence of drought-tolerant individuals in A. thaliana. This is one of the rare empirical examples substantiating the adaptive potential of mobilised stress-responsive TEs in eukaryotes. Our work demonstrates the potential of TE-mediated loss-of-function mutations in stress adaptation.

Identification and Characterization of PSEUDO-RESPONSE REGULATOR (PRR) 1a and 1b Genes by CRISPR/Cas9-Targeted Mutagenesis in Chinese Cabbage (Brassica rapa L.)

The CRISPR/Cas9 site-directed gene-editing system offers great advantages for identifying gene function and crop improvement. The circadian clock measures and conveys day length information to control rhythmic hypocotyl growth in photoperiodic conditions, to achieve optimal fitness, but operates through largely unknown mechanisms. Here, we generated core circadian clock evening components, Brassica rapa PSEUDO-RESPONSE REGULATOR (BrPRR) 1a, 1b, and 1ab (both 1a and 1b double knockout) mutants, using CRISPR/Cas9 genome editing in Chinese cabbage, where 9–16 genetic edited lines of each mutant were obtained. The targeted deep sequencing showed that each mutant had 2–4 different mutation types at the target sites in the BrPRR1a and BrPRR1b genes. To identify the functions of BrPRR1a and 1b genes, hypocotyl length, and mRNA and protein levels of core circadian clock morning components, BrCCA1 (CIRCADIAN CLOCK-ASSOCIATED 1) and BrLHY (LATE ELONGATED HYPOCOTYL) a and b were examined under light/dark cycles and continuous light conditions. The BrPRR1a and 1ab double mutants showed longer hypocotyls, lower core circadian clock morning component mRNA and protein levels, and a shorter circadian rhythm than wildtype (WT). On the other hand, the BrPRR1b mutant was not significantly different from WT. These results suggested that two paralogous genes may not be associated with the same regulatory function in Chinese cabbage. Taken together, our results demonstrated that CRISPR/Cas9 is an efficient tool for achieving targeted genome modifications and elucidating the biological functions of circadian clock genes in B. rapa, for both breeding and improvement.

Identification and Quantification of Glucosinolates and Phenolics in a Large Panel of Brassica napus Highlight Valuable Genetic Resources for Chemical Ecology and Breeding

Glucosinolate (GLS) and phenolic contents in Brassicaceae contribute to biotic and abiotic stress responses. Breeding crop accessions harboring agroecologically relevant metabolic profiles require a characterization of the chemical diversity in Brassica germplasm. This work investigates the diversity of specialized metabolites in 281 accessions of B. napus. First, an LC-HRMS2-based approach allowed the annotation of 32 phenolics and 36 GLSs, revealing 13 branched and linear alkyl-GLSs and 4 isomers of hydroxyphenylalkyl-GLSs, many of which have been rarely reported in Brassica. Then, quantitative UPLC-UV-MS-based profiling was performed in leaves and roots for the whole panel. This revealed striking variations in the content of 1-methylpropyl-GLS (glucocochlearin) and a large variation of tetra- and penta-glucosyl kaempferol derivatives among accessions. It also highlighted two main chemotypes related to sinapoyl-O-hexoside and kaempferol-O-trihexoside contents. By offering an unprecedented overview of the phytochemical diversity in B. napus, this work provides a useful resource for chemical ecology and breeding.

A keystone gene underlies the persistence of an experimental food web

Genes encode information that determines an organism's fitness. Yet we know little about whether genes of one species influence the persistence of interacting species in an ecological community. Here, we experimentally tested the effect of three plant defense genes on the persistence of an insect food web and found that a single allele at a single gene promoted coexistence by increasing plant growth rate, which in turn increased the intrinsic growth rates of species across multiple trophic levels. Our discovery of a "keystone gene" illustrates the need to bridge between biological scales, from genes to ecosystems, to understand community persistence.

Microevolution of Pieris butterfly genes involved in host-plant adaptation along a host-plant community cline

Herbivorous insects have evolved counteradaptations to overcome the chemical defenses of their host plants. Several of these counteradaptations have been elucidated at the molecular level, in particular for insects specialized on cruciferous host plants. While the importance of these counteradaptations for host plant colonization is well established, little is known about their microevolutionary dynamics in the field. In particular, it is not known whether and how host plant diversity shapes diversity in insect counteradaptations. In this study, we examine patterns of host plant use and insect counteradaptation in three Pieris butterfly species across Japan. The larvae of these butterflies express nitrile-specifier protein (NSP) and its paralog major allergen (MA) in their gut to overcome the highly diversified glucosinolate-myrosinase defense system of their cruciferous host plants. Pieris napi and Pieris melete colonize wild Brassicaceae whereas Pieris rapae typically uses cultivated Brassica as a host, regardless of the local composition of wild crucifers. As expected, NSP and MA diversity was independent of the local composition of wild Brassicaceae in P. rapae. In contrast, NSP diversity correlated with local host plant diversity in both species that preferred wild Brassicaceae. P. melete and P. napi both revealed two distinct major NSP alleles, which shaped diversity among local populations, albeit with different evolutionary trajectories. In comparison, MA showed no indication for local adaptation. Altogether, MA appeared to be evolutionary more conserved than NSP, suggesting that both genes play different roles in diverting host plant chemical defense.

Multigenerational Exposure to Heat Stress Induces Phenotypic Resilience, and Genetic and Epigenetic Variations in Arabidopsis thaliana Offspring

Plants are sedentary organisms that constantly sense changes in their environment and react to various environmental cues. On a short-time scale, plants respond through alterations in their physiology, and on a long-time scale, plants alter their development and pass on the memory of stress to the progeny. The latter is controlled genetically and epigenetically and allows the progeny to be primed for future stress encounters, thus increasing the likelihood of survival. The current study intended to explore the effects of multigenerational heat stress in Arabidopsis thaliana. Twenty-five generations of Arabidopsis thaliana were propagated in the presence of heat stress. The multigenerational stressed lineage F25H exhibited a higher tolerance to heat stress and elevated frequency of homologous recombination, as compared to the parallel control progeny F25C. A comparison of genomic sequences revealed that the F25H lineage had a three-fold higher number of mutations [single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs)] as compared control lineages, suggesting that heat stress induced genetic variations in the heat-stressed progeny. The F25H stressed progeny showed a 7-fold higher number of non-synonymous mutations than the F25C line. Methylome analysis revealed that the F25H stressed progeny showed a lower global methylation level in the CHH context than the control progeny. The F25H and F25C lineages were different from the parental control lineage F2C by 66,491 and 80,464 differentially methylated positions (DMPs), respectively. F25H stressed progeny displayed higher frequency of methylation changes in the gene body and lower in the body of transposable elements (TEs). Gene Ontology analysis revealed that CG-DMRs were enriched in processes such as response to abiotic and biotic stimulus, cell organizations and biogenesis, and DNA or RNA metabolism. Hierarchical clustering of these epimutations separated the heat stressed and control parental progenies into distinct groups which revealed the non-random nature of epimutations. We observed an overall higher number of epigenetic variations than genetic variations in all comparison groups, indicating that epigenetic variations are more prevalent than genetic variations. The largest difference in epigenetic and genetic variations was observed between control plants comparison (F25C vs. F2C), which clearly indicated that the spontaneous nature of epigenetic variations and heat-inducible nature of genetic variations. Overall, our study showed that progenies derived from multigenerational heat stress displayed a notable adaption in context of phenotypic, genotypic and epigenotypic resilience.

Structural Studies of Aliphatic Glucosinolate Chain-Elongation Enzymes

Plants evolved specialized metabolic pathways through gene duplication and functional divergence of enzymes involved in primary metabolism. The results of this process are varied pathways that produce an array of natural products useful to both plants and humans. In plants, glucosinolates are a diverse class of natural products. Glucosinolate function stems from their hydrolysis products, which are responsible for the strong flavors of Brassicales plants, such as mustard, and serve as plant defense molecules by repelling insects, fighting fungal infections, and discouraging herbivory. Additionally, certain hydrolysis products such as isothiocyanates can potentially serve as cancer prevention agents in humans. The breadth of glucosinolate function is a result of its great structural diversity, which comes from the use of aliphatic, aromatic and indole amino acids as precursors and elongation of some side chains by up to nine carbons, which, after the formation of the core glucosinolate structure, can undergo further chemical modifications. Aliphatic methionine-derived glucosinolates are the most abundant form of these compounds. Although both elongation and chemical modification of amino acid side chains are important for aliphatic glucosinolate diversity, its elongation process has not been well described at the molecular level. Here, we summarize new insights on the iterative chain-elongation enzymes methylthioalkylmalate synthase (MAMS) and isopropylmalate dehydrogenase (IPMDH).

The Roles of Cruciferae Glucosinolates in Disease and Pest Resistance

With the expansion of the area under Cruciferae vegetable cultivation, and an increase in the incidence of natural threats such as pests and diseases globally, Cruciferae vegetable losses caused by pathogens, insects, and pests are on the rise. As one of the key metabolites produced by Cruciferae vegetables, glucosinolate (GLS) is not only an indicator of their quality but also controls infestation by numerous fungi, bacteria, aphids, and worms. Today, the safe and pollution-free production of vegetables is advocated globally, and environmentally friendly pest and disease control strategies, such as biological control, to minimize the adverse impacts of pathogen and insect pest stress on Cruciferae vegetables, have attracted the attention of researchers. This review explores the mechanisms via which GLS acts as a defensive substance, participates in responses to biotic stress, and enhances plant tolerance to the various stress factors. According to the current research status, future research directions are also proposed.

Genetic variation, environment and demography intersect to shape Arabidopsis defense metabolite variation across Europe

Plants produce diverse metabolites to cope with the challenges presented by complex and ever-changing environments. These challenges drive the diversification of specialized metabolites within and between plant species. However, we are just beginning to understand how frequently new alleles arise controlling specialized metabolite diversity and how the geographic distribution of these alleles may be structured by ecological and demographic pressures. Here, we measure the variation in specialized metabolites across a population of 797 natural Arabidopsis thaliana accessions. We show that a combination of geography, environmental parameters, demography and different genetic processes all combine to influence the specific chemotypes and their distribution. This showed that causal loci in specialized metabolism contain frequent independently generated alleles with patterns suggesting potential within-species convergence. This provides a new perspective about the complexity of the selective forces and mechanisms that shape the generation and distribution of allelic variation that may influence local adaptation.

Since plants cannot move, they have evolved chemical defenses to help them respond to changes in their surroundings. For example, where animals run from predators, plants may produce toxins to put predators off. This approach is why plants are such a rich source of drugs, poisons, dyes and other useful substances. The chemicals plants produce are known as specialized metabolites, and they can change a lot between, and even within, plant species. The variety of specialized metabolites is a result of genetic changes and evolution over millions of years.

Evolution is a slow process, yet plants are able to rapidly develop new specialized metabolites to protect them from new threats. Even different populations of the same species produce many distinct metabolites that help them survive in their surroundings. However, the factors that lead plants to produce new metabolites are not well understood, and it is not known how this affects genetic variation.

To gain a better understanding of this process, Katz et al. studied 797 European variants of a common weed species called Arabidopsis thaliana, which is widely studied. The investigation found that many factors affect the range of specialized metabolites in each variant. These included local geography and environment, as well as genetics and population history (demography). Katz et al. revealed a pattern of relationships between the variants that could mirror their evolutionary history as the species spread and adapted to new locations.

These results highlight the complex network of factors that affect plant evolution. Rapid diversification is key to plant survival in new and changing environments and has resulted in a wide range of specialized metabolites. As such they are of interest both for studying plant evolution and for understanding their ecology. Expanding similar work to more populations and other species will broaden the scope of our ability to understand how plants adapt to their surroundings.

Engineering and optimization of the 2-phenylethylglucosinolate production in Nicotiana benthamiana by combining biosynthetic genes from Barbarea vulgaris and Arabidopsis thaliana

2-Phenylethylglucosinolate (2PE) derived from homophenylalanine is present in plants of the Brassicales order as a defense compound. It is associated with multiple biological properties, including deterrent effects on pests and antimicrobial and health-promoting functions, due to its hydrolysis product 2-phenylethyl isothiocyanate, which confers 2PE as a potential application in agriculture and industry. In this study, we characterized the putative key genes for 2PE biosynthesis from Barbarea vulgaris W.T. Aiton and demonstrated the feasibility of engineering 2PE production in Nicotiana benthamiana Domin. We used different combinations of genes from B. vulgaris and Arabidopsis thaliana (L.) Heynh. to demonstrate that: (i) BvBCAT4 performed more efficiently than AtBCAT4 in biosynthesis of both homophenylalanine and dihomomethionine; (ii) MAM1 enzymes were critical for the chain-elongated profile, while CYP79F enzymes accepted both chain-elongated methionine and homophenylalanine; (iii) aliphatic but not aromatic core structure pathway catalyzed the 2PE biosynthesis; (iv) a chimeric pathway containing BvBCAT4, BvMAM1, AtIPMI and AtIPMDH1 resulted in a two-fold increase in 2PE production compared with the B. vulgaris-specific chain elongation pathway; and (v) profiles of chain-elongated products and glucosinolates partially mirrored the profiles in the gene donor plant, but were wider in N. benthamiana than in the native plants. Our study provides a strategy to produce the important homophenylalanine and 2PE in a heterologous host. Furthermore, chimeric engineering of the complex 2PE biosynthetic pathway enabled detailed understanding of catalytic properties of individual enzymes - a prerequisite for understanding biochemical evolution. The new-to-nature gene combinations have the potential for application in biotechnological and plant breeding.

Evolutionary changes in the glucosinolate biosynthetic capacity in species representing Capsella, Camelina and Neslia genera

Glucosinolates are unique thioglucosides that evolved in the order Brassicales. These compounds function in plant adaptation to the environment, including combating plant pathogens, herbivore deterrence and abiotic stress tolerance. In line with their defensive functions glucosinolates usually accumulate constitutively in relatively high amounts in all tissues of Brassicaceae plants. Here we performed glucosinolate analysis in different organs of selected species representing Capsella, Camelina and Neslia genera, which similarly as the model plant Arabidopsis thaliana belong to the Camelineae tribe. We also identified orthologs of A. thaliana glucosinolate biosynthetic genes in the published genomes of some of the investigated species. Subsequent gene expression and phylogenetic analyses enabled us an insight into the evolutionary changes in the transcription of these genes and in the sequences of respective proteins that occurred within the Camelineae tribe. Our results indicated that glucosinolates are highly abundant in siliques and roots of the investigated species but hardly, if at all, produced in leaves. In addition to this unusual tissular distribution we revealed reduced structural diversity of methionine-derived aliphatic glucosinolates (AGs) with elevated accumulation of rare long chain AGs. This preference seems to correlate with evolutionary changes in genes encoding methylthioalkylmalate synthases that are responsible for the elongation of AG side chains. Finally, our results indicate that the biosynthetic pathway for tryptophan-derived indolic glucosinolates likely lost its main functions in immunity and resistance towards sucking insects and is on its evolutionary route to be shut off in the investigated species.

A Comprehensive Gene Inventory for Glucosinolate Biosynthetic Pathway in Arabidopsis thaliana

Glucosinolates (GSLs) are plant secondary metabolites comprising sulfur and nitrogen mainly found in plants from the order of Brassicales, such as broccoli, cabbage, and Arabidopsis thaliana. The activated forms of GSL play important roles in fighting against pathogens and have health benefits to humans. The increasing amount of data on A. thaliana generated from various omics technologies can be investigated more deeply in search of new genes or compounds involved in GSL biosynthesis and metabolism. This review describes a comprehensive inventory of A. thaliana GSLs identified from published literature and databases such as KNApSAcK, KEGG, and AraCyc. A total of 113 GSL genes encoding for 23 transcription components, 85 enzymes, and five protein transporters were experimentally characterized in the past two decades. Continuous efforts are still on going to identify all molecules related to the production of GSLs. A manually curated database known as SuCCombase (http://plant-scc.org) was developed to serve as a comprehensive GSL inventory. Realizing lack of information on the regulation of GSL biosynthesis and degradation mechanisms, this review also includes relevant information and their connections with crosstalk among various factors, such as light, sulfur metabolism, and nitrogen metabolism, not only in A. thaliana but also in other crucifers.

In Arabidopsis thaliana Substrate Recognition and Tissue- as Well as Plastid Type-Specific Expression Define the Roles of Distinct Small Subunits of Isopropylmalate Isomerase

In Arabidopsis thaliana, the heterodimeric isopropylmalate isomerase (IPMI) is composed of a single large (IPMI LSU1) and one of three different small subunits (IPMI SSU1 to 3). The function of IPMI is defined by the small subunits. IPMI SSU1 is required for Leu biosynthesis and has previously also been proposed to be involved in the first cycle of Met chain elongation, the first phase of the synthesis of Met-derived glucosinolates. IPMI SSU2 and IPMI SSU3 participate in the Met chain elongation pathway. Here, we investigate the role of the three IPMI SSUs through the analysis of the role of the substrate recognition region spanning five amino acids on the substrate specificity of IPMI SSU1. Furthermore, we analyze in detail the expression pattern of fluorophore-tagged IPMI SSUs throughout plant development. Our study shows that the substrate recognition region that differs between IPMI SSU1 and the other two IMPI SSUs determines the substrate preference of IPMI. Expression of IPMI SSU1 is spatially separated from the expression of IPMI SSU2 and IPMI SSU3, and IPMI SSU1 is found in small plastids, whereas IMPI SSU2 and SSU3 are found in chloroplasts. Our data show a distinct role for IMPI SSU1 in Leu biosynthesis and for IMPI SSU2 and SSU3 in the Met chain elongation pathway.

AthCNV: A Map of DNA Copy Number Variations in the Arabidopsis Genome

Copy number variations (CNVs) greatly contribute to intraspecies genetic polymorphism and phenotypic diversity. Recent analyses of sequencing data for >1000 Arabidopsis (Arabidopsis thaliana) accessions focused on small variations and did not include CNVs. Here, we performed genome-wide analysis and identified large indels (50 to 499 bp) and CNVs (500 bp and larger) in these accessions. The CNVs fully overlap with 18.3% of protein-coding genes, with enrichment for evolutionarily young genes and genes involved in stress and defense. By combining analysis of both genes and transposable elements (TEs) affected by CNVs, we revealed that the variation statuses of genes and TEs are tightly linked and jointly contribute to the unequal distribution of these elements in the genome. We also determined the gene copy numbers in a set of 1060 accessions and experimentally validated the accuracy of our predictions by multiplex ligation-dependent probe amplification assays. We then successfully used the CNVs as markers to analyze population structure and migration patterns. Finally, we examined the impact of gene dosage variation triggered by a CNV spanning the SEC10 gene on SEC10 expression at both the transcript and protein levels. The catalog of CNVs, CNV-overlapping genes, and their genotypes in a top model dicot will stimulate the exploration of the genetic basis of phenotypic variation.

mGWAS Uncovers Gln-Glucosinolate Seed-Specific Interaction and its Role in Metabolic Homeostasis

Gln is a key player in plant metabolism. It is one of the major free amino acids that is transported into the developing seed and is central for nitrogen metabolism. However, Gln natural variation and its regulation and interaction with other metabolic processes in seeds remain poorly understood. To investigate the latter, we performed a metabolic genome-wide association study (mGWAS) of Gln-related traits measured from the dry seeds of the Arabidopsis (Arabidopsis thaliana) diversity panel using all potential ratios between Gln and the other members of the Glu family as traits. This semicombinatorial approach yielded multiple candidate genes that, upon further analysis, revealed an unexpected association between the aliphatic glucosinolates (GLS) and the Gln-related traits. This finding was confirmed by an independent quantitative trait loci mapping and statistical analysis of the relationships between the Gln-related traits and the presence of specific GLS in seeds. Moreover, an analysis of Arabidopsis mutants lacking GLS showed an extensive seed-specific impact on Gln levels and composition that manifested early in seed development. The elimination of GLS in seeds was associated with a large effect on seed nitrogen and sulfur homeostasis, which conceivably led to the Gln response. This finding indicates that both Gln and GLS play key roles in shaping the seed metabolic homeostasis. It also implies that select secondary metabolites might have key functions in primary seed metabolism. Finally, our study shows that an mGWAS performed on dry seeds can uncover key metabolic interactions that occur early in seed development.

Dissection of genetic architecture for glucosinolate accumulations in leaves and seeds of Brassica napus by genome-wide association study

Glucosinolates (GSLs), whose degradation products have been shown to be increasingly important for human health and plant defence, compose important secondary metabolites found in the order Brassicales. It is highly desired to enhance pest and disease resistance by increasing the leaf GSL content while keeping the content low in seeds of Brassica napus, one of the most important oil crops worldwide. Little is known about the regulation of GSL accumulation in the leaves. We quantified the levels of 9 different GSLs and 15 related traits in the leaves of 366 accessions and found that the seed and leaf GSL content were highly correlated (r = 0.79). A total of 78 loci were associated with GSL traits, and five common and eleven tissue-specific associated loci were related to total leaf and seed GSL content. Thirty-six candidate genes were inferred to be involved in GSL biosynthesis. The candidate gene BnaA03g40190D (BnaA3.MYB28) was validated by DNA polymorphisms and gene expression analysis. This gene was responsible for high leaf/low seed GSL content and could explain 30.62% of the total leaf GSL variation in the low seed GSL panel and was not fixed during double-low rapeseed breeding. Our results provide new insights into the genetic basis of GSL variation in leaves and seeds and may facilitate the metabolic engineering of GSLs and the breeding of high leaf/low seed GSL content in B. napus.

Independent evolution of ancestral and novel defenses in a genus of toxic plants (Erysimum, Brassicaceae)

Phytochemical diversity is thought to result from coevolutionary cycles as specialization in herbivores imposes diversifying selection on plant chemical defenses. Plants in the speciose genus Erysimum (Brassicaceae) produce both ancestral glucosinolates and evolutionarily novel cardenolides as defenses. Here we test macroevolutionary hypotheses on co-expression, co-regulation, and diversification of these potentially redundant defenses across this genus. We sequenced and assembled the genome of E. cheiranthoides and foliar transcriptomes of 47 additional Erysimum species to construct a phylogeny from 9868 orthologous genes, revealing several geographic clades but also high levels of gene discordance. Concentrations, inducibility, and diversity of the two defenses varied independently among species, with no evidence for trade-offs. Closely related, geographically co-occurring species shared similar cardenolide traits, but not glucosinolate traits, likely as a result of specific selective pressures acting on each defense. Ancestral and novel chemical defenses in Erysimum thus appear to provide complementary rather than redundant functions.

Using networks to identify and interpret natural variation

Studies on natural variation and network biology inherently work to summarize vast amounts of information and data. The combination of these two areas of study while creating datasets of immense complexity is critical to their mutual progress. Networks are necessary as a way to work to reduce the dimensionality inherent in natural variation with 100 s to 1000 s of genotypes. Correspondingly natural variation is essential for testing how networks may or may not be shared across individuals or species. Advances in this area of cross-fertilization including using networks directly as phenotypes and the use of networks to help in prioritizing candidate gene validation efforts. Interesting new observations on frequent presence-absence variation in gene content and adaptation is beginning to highlight the potential for natural variation in network presence-absence. This review attempts to delve into these new insights.

Overexpression of HMG-CoA synthase promotes Arabidopsis root growth and adversely affects glucosinolate biosynthesis

3-Hydroxy-3-methylglutaryl-CoA synthase (HMGS) catalyses the second step of the mevalonate (MVA) pathway. An HMGS inhibitor (F-244) has been reported to retard growth in wheat, tobacco, and Brassica juncea, but the mechanism remains unknown. Although the effects of HMGS on downstream isoprenoid metabolites have been extensively reported, not much is known on how it might affect non-isoprenoid metabolic pathways. Here, the mechanism of F-244-mediated inhibition of primary root growth in Arabidopsis and the relationship between HMGS and non-isoprenoid metabolic pathways were investigated by untargeted SWATH-MS quantitative proteomics, quantitative real-time PCR, and target metabolite analysis. Our results revealed that the inhibition of primary root growth caused by F-244 was a consequence of reduced stigmasterol, auxin, and cytokinin levels. Interestingly, proteomic analyses identified a relationship between HMGS and glucosinolate biosynthesis. Inhibition of HMGS activated glucosinolate biosynthesis, resulting from the induction of glucosinolate biosynthesis-related genes, suppression of sterol biosynthesis-related genes, and reduction in sterol levels. In contrast, HMGS overexpression inhibited glucosinolate biosynthesis, due to down-regulation of glucosinolate biosynthesis-related genes, up-regulation of sterol biosynthesis-related genes, and increase in sterol content. Thus, HMGS might represent a target for the manipulation of glucosinolate biosynthesis, given the regulatory relationship between HMGS in the MVA pathway and glucosinolate biosynthesis.

Transcriptional Variation in Glucosinolate Biosynthetic Genes and Inducible Responses to Aphid Herbivory on Field-Grown Arabidopsis thaliana

Recently, increasing attempts have been made to understand how plant genes function in natura. In this context, transcriptional profiles represent plant physiological status in response to environmental stimuli. Herein, we combined high-throughput RNA-Seq with insect survey data on 19 accessions of Arabidopsis thaliana grown at a field site in Switzerland. We found that genes with the gene ontology (GO) annotations of "glucosinolate biosynthetic process" and "response to insects" were most significantly enriched, and the expression of these genes was highly variable among plant accessions. Nearly half of the total expression variation in the glucosinolate biosynthetic genes (AOPs, ESM1, ESP, and TGG1) was explained by among-accession variation. Of these genes, the expression level of AOP3 differed among Col-0 accession individuals depending on the abundance of the mustard aphid (Lipaphis erysimi). We also found that the expression of the major cis-jasmone activated gene CYP81D11 was positively correlated with the number of flea beetles (Phyllotreta striolata and Phyllotreta atra). Combined with the field RNA-Seq data, bioassays confirmed that AOP3 was up-regulated in response to attack by mustard aphids. The combined results from RNA-Seq and our ecological survey illustrate the feasibility of using field transcriptomics to detect an inducible defense, providing a first step towards an in natura understanding of biotic interactions involving phenotypic plasticity.

Molecular Basis of the Evolution of Methylthioalkylmalate Synthase and the Diversity of Methionine-Derived Glucosinolates

The globally cultivated Brassica species possess diverse aliphatic glucosinolates, which are important for plant defense and animal nutrition. The committed step in the side chain elongation of methionine-derived aliphatic glucosinolates is catalyzed by methylthioalkylmalate synthase, which likely evolved from the isopropylmalate synthases of leucine biosynthesis. However, the molecular basis for the evolution of methylthioalkylmalate synthase and its generation of natural product diversity in Brassica is poorly understood. Here, we show that Brassica genomes encode multiple methylthioalkylmalate synthases that have differences in expression profiles and 2-oxo substrate preferences, which account for the diversity of aliphatic glucosinolates across Brassica accessions. Analysis of the 2.1 Å resolution x-ray crystal structure of Brassica juncea methylthioalkylmalate synthase identified key active site residues responsible for controlling the specificity for different 2-oxo substrates and the determinants of side chain length in aliphatic glucosinolates. Overall, these results provide the evolutionary and biochemical foundation for the diversification of glucosinolate profiles across globally cultivated Brassica species, which could be used with ongoing breeding strategies toward the manipulation of beneficial glucosinolate compounds for animal health and plant protection.

Glucosinolate variability between turnip organs during development

Turnip (Brassica rapa spp. rapa) is an important vegetable species, with a unique physiology. Several plant parts, including both the turnip tubers and leaves, are important for human consumption. During the development of turnip plants, the leaves function as metabolic source tissues, while the tuber first functions as a sink, while later the tuber turns into a source for development of flowers and seeds. In the present study, chemical changes were determined for two genotypes with different genetic background, and included seedling, young leaves, mature leaves, tuber surface, tuber core, stalk, flower and seed tissues, at seven different time points during plant development. As a basis for understanding changes in glucosinolates during plant development, the profile of glucosinolates was analysed using liquid chromatography (LC) coupled to mass spectrometry (MS). This analysis was complemented by a gene expression analysis, focussed on GLS biosynthesis, which could explain part of the observed variation, pointing to important roles of specific gene orthologues for defining the chemical differences. Substantial differences in glucosinolate profiles were observed between above-ground tissues and turnip tuber, reflecting the differences in physiological role. In addition, differences between the two genotypes and between tissues that were harvested early or late during the plant lifecycle. The importance of the observed differences in glucosinolate profile for the ecophysiology of the turnip and for breeding turnips with optimal chemical profiles is discussed.

Genetic Interactions and Molecular Evolution of the Duplicated Genes ICARUS2 and ICARUS1 Help Arabidopsis Plants Adapt to Different Ambient Temperatures

Understanding how plants adapt to ambient temperatures has become a major challenge prompted by global climate change. This has led to the identification of several genes regulating the thermal plasticity of plant growth and flowering time. However, the mechanisms accounting for the natural variation and evolution of such developmental plasticity remain mostly unknown. In this study, we determined that natural variation at ICARUS2 (ICA2), which interacts genetically with its homolog ICA1, alters growth and flowering time plasticity in relation to temperature in Arabidopsis (Arabidopsis thaliana). Transgenic analyses demonstrated multiple functional effects for ICA2 and supported the notion that structural polymorphisms in ICA2 likely underlie its natural variation. Two major ICA2 haplogroups carrying distinct functionally active alleles showed high frequency, strong geographic structure, and significant associations with climatic variables related to annual and daily fluctuations in temperature. Genome analyses across the plant phylogeny indicated that the prevalent plant ICA genes encoding two tRNAHis guanylyl transferase 1 units evolved ∼120 million years ago during the early divergence of mono- and dicotyledonous clades. In addition, ICA1/ICA2 duplication occurred specifically in the Camelineae tribe (Brassicaceae). Thus, ICA2 appears to be ubiquitous across plant evolution and likely contributes to climate adaptation through modifications of thermal developmental plasticity in Arabidopsis.

Evolution of Glucosinolate Diversity via Whole-Genome Duplications, Gene Rearrangements, and Substrate Promiscuity

Over several decades, glucosinolates have become a model system for the study of specialized metabolic diversity in plants. The near-complete identification of biosynthetic enzymes, regulators, and transporters has provided support for the role of gene duplication and subsequent changes in gene expression, protein function, and substrate specificity as the evolutionary bases of glucosinolate diversity. Here, we provide examples of how whole-genome duplications, gene rearrangements, and substrate promiscuity potentiated the evolution of glucosinolate biosynthetic enzymes, regulators, and transporters by natural selection. This in turn may have led to the repeated evolution of glucosinolate metabolism and diversity in higher plants.

Glucosinolates: Molecular structure, breakdown, genetic, bioavailability, properties and healthy and adverse effects

Glucosinolates are a large group of plant secondary metabolites with nutritional effects and biologically active compounds. Glucosinolates are mainly found in cruciferous plants such as Brassicaceae family, including common edible plants such as broccoli (Brassica oleracea var. italica), cabbage (B. oleracea var. capitata f. alba), cauliflower (B. oleracea var. botrytis), rapeseed (Brassica napus), mustard (Brassica nigra), and horseradish (Armoracia rusticana). If cruciferous plants are consumed without processing, myrosinase enzyme will hydrolyze the glucosinolates to various metabolites, such as isothiocyanates, nitriles, oxazolidine-2-thiones, and indole-3-carbinols. On the other hand, when cruciferous are cooked before consumption, myrosinase is inactivated and glucosinolates could be partially absorbed in their intact form through the gastrointestinal mucosa. This review paper summarizes the glucosinolate molecular breakdown, their genetic aspects from biosynthesis to precursors, their bioavailability (assimilation, absorption, and elimination of these molecules), their sensory properties, identified healthy and adverse effects, as well as the impact of processing on their bioavailability.

Drought adaptation in Arabidopsis thaliana by extensive genetic loss-of-function

Interdisciplinary syntheses are needed to scale up discovery of the environmental drivers and molecular basis of adaptation in nature. Here we integrated novel approaches using whole genome sequences, satellite remote sensing, and transgenic experiments to study natural loss-of-function alleles associated with drought histories in wild Arabidopsis thaliana. The genes we identified exhibit population genetic signatures of parallel molecular evolution, selection for loss-of-function, and shared associations with flowering time phenotypes in directions consistent with longstanding adaptive hypotheses seven times more often than expected by chance. We then confirmed predicted phenotypes experimentally in transgenic knockout lines. These findings reveal the importance of drought timing to explain the evolution of alternative drought tolerance strategies and further challenge popular assumptions about the adaptive value of genetic loss-of-function in nature. These results also motivate improved species-wide sequencing efforts to better identify loss-of-function variants and inspire new opportunities for engineering climate resilience in crops.

Water shortages caused by droughts lead to crop losses that affect billions of people around the world each year. By discovering how wild plants adapt to drought, it may be possible to identify traits and genes that help to improve the growth of crop plants when water is scarce. It has been suggested that plants have adapted to droughts by flowering at times of the year when droughts are less likely to occur. For example, if droughts are more likely to happen in spring, the plants may delay flowering until the summer.

Arabidopsis thaliana is a small plant that is found across Eurasia, Africa and North America, including in areas that are prone to drought at different times of the year. Individual plants of the same species may carry different versions of the same gene (known as alleles). Some of these alleles may not work properly and are referred to as loss-of-function alleles. Monroe et al. investigated whether A. thaliana plants carry any loss-of-function alleles that are associated with droughts happening in the spring or summer, and whether they are linked to when those plants will flower.

Monroe et al. analyzed satellite images collected over the last 30 years to measure when droughts have occurred. Next, they searched genome sequences of Arabidopsis thaliana for alleles that might help the plants to adapt to droughts in the spring or summer. Combining the two approaches revealed that loss-of-function alleles associated with spring droughts were strongly predicted to be associated with the plants flowering later in the year. Similarly, loss-of-function alleles associated with summer droughts were predicted to be associated with the plants flowering earlier in the year.

These findings support the idea that plants can adapt to drought by changing when they produce flowers, and suggest that loss-of-function alleles play a major role in this process. New techniques for editing genes mean it is easier than ever to generate new loss-of-function alleles in specific genes. Therefore, the results presented by Monroe et al. may help researchers to develop new varieties of crop plants that are better adapted to droughts.

Plant Secondary Metabolite Diversity and Species Interactions

Ever since the first plant secondary metabolites (PSMs) were isolated and identified, questions about their ecological functions and diversity have been raised. Recent advances in analytical chemistry and complex data computation, as well as progress in chemical ecology from mechanistic to functional and evolutionary questions, open a new box of hypotheses. Addressing these hypotheses includes the measurement of complex traits, such as chemodiversity, in a context-dependent manner and allows for a deeper understanding of the multifunctionality and functional redundancy of PSMs. Here we review a hypothesis framework that addresses PSM diversity on multiple ecological levels (α, β, and γ chemodiversity), its variation in space and time, and the potential agents of natural selection. We use the concept of chemical information transfer as mediator of antagonistic and mutualistic interaction to interpret functional and microevolutionary studies and create a hypothesis framework for understanding chemodiversity as a factor driving ecological processes.

Balancing selection for aflatoxin in Aspergillus flavus is maintained through interference competition with, and fungivory by insects

The role of microbial secondary metabolites in the ecology of the organisms that produce them remains poorly understood. Variation in aflatoxin production by Aspergillus flavus is maintained by balancing selection, but the ecological function and impact on fungal fitness of this compound are unknown. We hypothesize that balancing selection for aflatoxin production in A. flavus is driven by interaction with insects. To test this, we competed naturally occurring aflatoxigenic and non-aflatoxigenic fungal isolates against Drosophila larvae on medium containing 0-1750 ppb aflatoxin, using quantitative PCR to quantify A. flavus DNA as a proxy for fungal fitness. The addition of aflatoxin across this range resulted in a 26-fold increase in fungal fitness. With no added toxin, aflatoxigenic isolates caused higher mortality of Drosophila larvae and had slightly higher fitness than non-aflatoxigenic isolates. Additionally, aflatoxin production increased an average of 1.5-fold in the presence of a single larva and nearly threefold when the fungus was mechanically damaged. We argue that the role of aflatoxin in protection from fungivory is inextricably linked to its role in interference competition. Our results, to our knowledge, provide the first clear evidence of a fitness advantage conferred to A. flavus by aflatoxin when interacting with insects.

Comparison of the Relative Potential for Epigenetic and Genetic Variation To Contribute to Trait Stability

The theoretical ability of epigenetic variation to influence the heritable variation of complex traits is gaining traction in the study of adaptation. This theory posits that epigenetic marks can control adaptive phenotypes but the relative potential of epigenetic variation in comparison to genetic variation in these traits is not presently understood. To compare the potential of epigenetic and genetic variation in adaptive traits, we analyzed the influence of DNA methylation variation on the accumulation of chemical defense compounds glucosinolates from the order Brassicales. Several decades of work on glucosinolates has generated extensive knowledge about their synthesis, regulation, genetic variation and contribution to fitness establishing this pathway as a model pathway for complex adaptive traits. Using high-throughput phenotyping with a randomized block design of ddm1 derived Arabidopsis thaliana epigenetic Recombinant Inbred Lines, we measured the correlation between DNA methylation variation and mean glucosinolate variation and within line stochastic variation. Using this information, we identified epigenetic Quantitative Trait Loci that contained specific Differentially Methylated Regions associated with glucosinolate traits. This showed that variation in DNA methylation correlates both with levels and variance of glucosinolates and flowering time with trait-specific loci. By conducting a meta-analysis comparing the results to different genetically variable populations, we conclude that the influence of DNA methylation variation on these adaptive traits is much lower than the corresponding impact of standing genetic variation. As such, selective pressure on these traits should mainly affect standing genetic variation to lead to adaptation.

Reconstructing Biosynthetic Pathway of the Plant-Derived Cancer Chemopreventive-Precursor Glucoraphanin in Escherichia coli

Epidemiological data confirmed a strong correlation between regular consumption of cruciferous vegetables and lower cancer risk. This cancer preventive property is mainly attributed to the glucosinolate products, such as glucoraphanin found in broccoli that is derived from methionine. Here we report the first successful reconstruction of the complete biosynthetic pathway of glucoraphanin from methionine in Escherichia coli via gene selection, pathway design, and protein engineering. We used branched-chain amino transferase 3 to catalyze two transamination steps to ensure the purity of precursor molecules and used cysteine as a sulfur donor to simplify the synthesis pathway. Two chimeric cytochrome P450 enzymes were engineered and expressed in E. coli functionally. The original plant C-S lyase was replaced by the Neurospora crassa hercynylcysteine sulfoxide lyase. Other pathway enzymes were successfully mined from Arabidopsis thaliana, Brassica rapa, and Brassica oleracea. Biosynthesis of glucoraphanin upon coexpression of the optimized enzymes in vivo was confirmed by liquid chromatography-tandem mass spectrometry analysis. No other glucosinolate analogues (except for glucoiberin) were identified that could facilitate the downstream purification processes. Production of glucoraphanin in this study laid the foundation for microbial production of such health-beneficial glucosinolates in a large-scale.

A naturally occurring variation in the BrMAM-3 gene is associated with aliphatic glucosinolate accumulation in Brassica rapa leaves

Glucosinolate profiles significantly vary among Brassica rapa genotypes. However, the molecular basis of these variations is largely unknown. In this study, we investigated a major quantitative trait locus (QTL) controlling aliphatic glucosinolate accumulation in B. rapa leaves. The QTL, which encompasses three tandem MAM genes and two MYB genes, was detected in two BC₂DH populations. Among the five-candidate genes, only the expression level of BrMAM-3 (Bra013007) was significantly correlated with the accumulation of aliphatic glucosinolates in B. rapa leaves. We identified a naturally occurring insertion within exon 1 of BrMAM-3, which is predicted to be a loss-of-function mutation, as confirmed by qRT-PCR. We determined that the loss of function was associated with the low glucosinolate content in B. rapa accessions. Furthermore, overexpressing the BrMAM-3 gene resulted in an increase in total aliphatic glucosinolates in Arabidopsis transgenic lines. Our study provides insights into the molecular mechanism underlying the accumulation of aliphatic glucosinolates in B. rapa leaves, thereby facilitating in the manipulation of total aliphatic glucosinolate content in Brassica crops.

Network-Guided Discovery of Extensive Epistasis between Transcription Factors Involved in Aliphatic Glucosinolate Biosynthesis

Plants use diverse mechanisms influenced by vast regulatory networks of indefinite scale to adapt to their environment. These regulatory networks have an unknown potential for epistasis between genes within and across networks. To test for epistasis within an adaptive trait genetic network, we generated and tested 47 Arabidopsis thaliana double mutant combinations for 20 transcription factors, which all influence the accumulation of aliphatic glucosinolates, the defense metabolites that control fitness. The epistatic combinations were used to test if there is more or less epistasis depending on gene membership within the same or different phenotypic subnetworks. Extensive epistasis was observed between the transcription factors, regardless of subnetwork membership. Metabolite accumulation displayed antagonistic epistasis, suggesting the presence of a buffering mechanism. Epistasis affecting enzymatic estimated activity was highly conditional on the tissue and environment and shifted between both antagonistic and synergistic forms. Transcriptional analysis showed that epistasis shifts depend on how the trait is measured. Because the 47 combinations described here represent a small sampling of the potential epistatic combinations in this genetic network, there is potential for significantly more epistasis. Additionally, the main effect of the individual gene was not predictive of the epistatic effects, suggesting that there is a need for further studies.

The Effect of Single and Multiple SERAT Mutants on Serine and Sulfur Metabolism

The gene family of serine acetyltransferases (SERATs) constitutes an interface between the plant pathways of serine and sulfur metabolism. SERATs provide the activated precursor, O-acetylserine for the fixation of reduced sulfur into cysteine by exchanging the serine hydroxyl moiety by a sulfhydryl moiety, and subsequently all organic compounds containing reduced sulfur moieties. We investigate here, how manipulation of the SERAT interface results in metabolic alterations upstream or downstream of this boundary and the extent to which the five SERAT isoforms exert an effect on the coupled system, respectively. Serine is synthesized through three distinct pathways while cysteine biosynthesis is distributed over the three compartments cytosol, mitochondria, and plastids. As the respective mutants are viable, all necessary metabolites can obviously cross various membrane systems to compensate what would otherwise constitute a lethal failure in cysteine biosynthesis. Furthermore, given that cysteine serves as precursor for multiple pathways, cysteine biosynthesis is highly regulated at both, the enzyme and the expression level. In this study, metabolite profiles of a mutant series of the SERAT gene family displayed that levels of the downstream metabolites in sulfur metabolism were affected in correlation with the reduction levels of SERAT activities and the growth phenotypes, while levels of the upstream metabolites in serine metabolism were unchanged in the serat mutants compared to wild-type plants. These results suggest that despite of the fact that the two metabolic pathways are directly connected, there seems to be no causal link in metabolic alterations. This might be caused by the difference of their pool sizes or the tight regulation by homeostatic mechanisms that control the metabolite concentration in plant cells. Additionally, growth conditions exerted an influence on metabolic compositions.