Lineage-specific gene duplication and expansion of DUF1216 gene family in Brassicaceae.

Proteins containing domain of unknown function (DUF) are prevalent in eukaryotic genome. The DUF1216 proteins possess a conserved DUF1216 domain resembling to the mediator protein of Arabidopsis RNA polymerase II transcriptional subunit-like protein. The DUF1216 family are specifically existed in Brassicaceae, however, no comprehensive evolutionary analysis of DUF1216 genes have been performed. We performed a first comprehensive genome-wide analysis of DUF1216 proteins in Brassicaceae. Totally 284 DUF1216 genes were identified in 27 Brassicaceae species and classified into four subfamilies on the basis of phylogenetic analysis. The analysis of gene structure and conserved motifs revealed that DUF1216 genes within the same subfamily exhibited similar intron/exon patterns and motif composition. The majority members of DUF1216 genes contain a signal peptide in the N-terminal, and the ninth position of the signal peptide in most DUF1216 is cysteine. Synteny analysis revealed that segmental duplication is a major mechanism for expanding of DUF1216 genes in Brassica oleracea, Brassica juncea, Brassica napus, Lepidium meyneii, and Brassica carinata, while in Arabidopsis thaliana and Capsella rubella, tandem duplication plays a major role in the expansion of the DUF1216 gene family. The analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios for DUF1216 paralogous indicated that most of gene pairs underwent purifying selection. DUF1216 genes displayed a specifically high expression in reproductive tissues in most Brassicaceae species, while its expression in Brassica juncea was specifically high in root. Our studies offered new insights into the phylogenetic relationships, gene structures and expressional patterns of DUF1216 members in Brassicaceae, which provides a foundation for future functional analysis.

A radiation of Psylliodes flea beetles on Brassicaceae is associated with the evolution of specific detoxification enzymes.

Flea beetles of the genus Psylliodes have evolved specialized interactions with plant species belonging to several distantly related families, mainly Brassicaceae, Solanaceae, and Fagaceae. This diverse host use indicates that Psylliodes flea beetles are able to cope with different chemical defense metabolites, including glucosinolates, the characteristic defense metabolites of Brassicaceae. Here we investigated the evolution of host use and the emergence of a glucosinolate-specific detoxification mechanism in Psylliodes flea beetles. In phylogenetic analyses, Psylliodes species clustered into four major clades, three of which contained mainly species specialized on either Brassicaceae, Solanaceae, or Fagaceae. Most members of the fourth clade have broader host use, including Brassicaceae and Poaceae as major host plant families. Ancestral state reconstructions suggest that Psylliodes flea beetles were initially associated with Brassicaceae and then either shifted to Solanaceae or Fagaceae, or expanded their host repertoire to Poaceae. Despite a putative ancestral association with Brassicaceae, we found evidence that the evolution of glucosinolate-specific detoxification enzymes coincides with the radiation of Psylliodes on Brassicaceae, suggesting that these are not required for using Brassicaceae as hosts but could improve the efficiency of host use by specialized Psylliodes species.

Overexpression of gamma-glutamyl cyclotransferase 2;1 (CsGGCT2;1) reduces arsenic toxicity and accumulation in Camelina sativa (L.).

KEY MESSAGE: Overexpressing CsGGCT2;1 in Camelina enhances arsenic tolerance, reducing arsenic accumulation by 40-60%. Genetically modified Camelina can potentially thrive on contaminated lands and help safeguard food quality and sustainable food and biofuel production. Environmental arsenic contamination is a serious global issue that adversely affects human health and diminishes the quality of harvested produce. Glutathione (GSH) is known to bind and detoxify arsenic and other toxic metals. A steady level of GSH is maintained within cells via the γ-glutamyl cycle. The γ-glutamyl cyclotransferases (GGCTs) have previously been shown to be involved in GSH degradation and increased tolerance to toxic metals in plants. In this study, we characterized the GGCT2;1 homolog from Camelina sativa for its role in arsenic tolerance and accumulation. Overexpression of CsGGCT2;1 in Camelina under CaMV35S constitutive promoter resulted in strong tolerance to arsenite (AsIII). The overexpression (OE) lines had 2.6-3.5-fold higher shoots and sevenfold to tenfold enhanced root biomass on media supplemented with AsIII, relative to wild-type plants. The CsGGCT2;1 OE lines accumulated 40-60% less arsenic in root and shoot tissues compared to wild-type plants. Further, the OE lines had ~ twofold higher chlorophyll content and 35% lesser levels of malondialdehyde (MDA), an indicator of membrane damage via lipid peroxidation. There was a slight but non-significant increase in 5-oxoproline (5-OP), a product of GSH degradation, in OE lines. However, the transcript levels of Oxoprolinase 1 (OXP1) were upregulated, indicating the accelerated conversion of 5-OP to glutamate, which is further utilized for the resynthesis of GSH to maintain GSH homeostasis. Overall, this research suggests that genetically modified Camelina may have the potential for cultivation on contaminated marginal lands to reduce As accumulation; thereby could help in addressing food safety issues as well as future food and biofuel needs.

Antagonistic RALF peptides control an intergeneric hybridization barrier on Brassicaceae stigmas.

Pollen-pistil interactions establish interspecific/intergeneric pre-zygotic hybridization barriers in plants. The rejection of undesired pollen at the stigma is crucial to avoid outcrossing but can be overcome with the support of mentor pollen. The mechanisms underlying this hybridization barrier are largely unknown. Here, in Arabidopsis, we demonstrate that receptor-like kinases FERONIA/CURVY1/ANJEA/HERCULES RECEPTOR KINASE 1 and cell wall proteins LRX3/4/5 interact on papilla cell surfaces with autocrine stigmatic RALF1/22/23/33 peptide ligands (sRALFs) to establish a lock that blocks the penetration of undesired pollen tubes. Compatible pollen-derived RALF10/11/12/13/25/26/30 peptides (pRALFs) act as a key, outcompeting sRALFs and enabling pollen tube penetration. By treating Arabidopsis stigmas with synthetic pRALFs, we unlock the barrier, facilitating pollen tube penetration from distantly related Brassicaceae species and resulting in interspecific/intergeneric hybrid embryo formation. Therefore, we uncover a "lock-and-key" system governing the hybridization breadth of interspecific/intergeneric crosses in Brassicaceae. Manipulating this system holds promise for facilitating broad hybridization in crops.

Overexpression of BoLSU1 and BoLSU2 Confers Tolerance to Sulfur Deficiency in Arabidopsis by Manipulating Glucosinolate Metabolism.

Sulfur is an essential element for plant growth, development and resistance to environmental stresses. Glucosinolates (GSLs), a group of sulfur rich secondary metabolites found in Brassicaceae plants, are known for their defensive properties against pathogens and herbivores. Due to their integration of a large proportion of total sulfur, their biosynthesis and degradation are closely linked to sulfur metabolism. It has been demonstrated that GSLs can be broken down to release sulfur and facilitate the production of other thio-metabolites when the plant is under stress. However, the regulation of this process is still not fully understood. In this study, we constructed two broccoli LSU (low sulfur responsive) gene overexpressing lines, 35S::BoLSU1 and 35S::BoLSU2, to detect changes in GSL metabolism after sulfur deficiency treatment. The results showed that BoLSU1 and BoLSU2 inhibit the biosynthesis of aliphatic GSLs, while also promoting their degradation and increasing the content of glutathione (GSH), leading to the reallocation of sulfur from the GSL pool to other thio-metabolites such as GSH. Furthermore, this regulation of GSL metabolism mediated by BoLSU1 and BoLSU2 is found to be dependent on myrosinases BGLU28 and BGLU30. Our study provides insight into the physiological role of LSU proteins and their regulation of sulfur metabolism.

Cadmium and Zn hyperaccumulation provide efficient constitutive defense against Turnip yellow mosaic virus infection in Noccaea caerulescens.

To understand the role of Zn and Cd in anti-viral defence, Zn/Cd hyperaccumulator Noccaea caerulescens plants grown with deficient (0.3 µM), replete (10 µM) and excess (100 µM) Zn2+ and Cd (10 µM Zn2+ + 1 µM Cd2+) were infected with Turnip yellow mosaic virus (TYMV). Gas exchange and chlorophyll fluorescence kinetics analyses demonstrated direct TYMV effects on photosynthetic light reactions but N. caerulescens was more resistant against TYMV than the previously studied non-hyperaccumulator N. ochroleucum. Virus abundance and photosynthesis inhibition were the lowest in the high Zn and Cd treatments. RNAseq analysis of 10 µM Zn2+ plants revealed TYMV-induced upregulation of Ca transporters, chloroplastic ZTP29 and defence genes, but none of those that are known to be strongly involved in hyperaccumulation. Synchrotron µ-XRF tomography, however, showed that Zn hyperaccumulation remained strongest in vacuoles of epidermal storage cells regardless of infection. This was in contrast to N. ochroleucum, where apoplastic Zn drastically increased in response to TYMV. These results suggest that the antiviral response of N. caerulescens is less induced by the onset of this biotic stress, but it is rather a permanent resistant state of the plant. Real-time qPCR revealed upregulation of ferritin in Zn10 infected plants, suggesting Fe deprivation as a virus defence strategy under suboptimal Zn supply.

The meso-octoploid Heliophila variabilis genome sheds a new light on the impact of polyploidization and diploidization on the diversity of the Cape flora.

Although the South African Cape flora is one of the most remarkable biodiversity hotspots, its high diversity has not been associated with polyploidy. Here, we report the chromosome-scale genome assembly of an ephemeral cruciferous species Heliophila variabilis (~334 Mb, n = 11) adapted to South African semiarid biomes. Two pairs of differently fractionated subgenomes suggest an allo-octoploid origin of the genome at least 12 million years ago. The ancestral octoploid Heliophila genome (2n = 8x = ~60) has probably originated through hybridization between two allotetraploids (2n = 4x = ~30) formed by distant, intertribal, hybridization. Rediploidization of the ancestral genome was marked by extensive reorganization of parental subgenomes, genome downsizing, and speciation events in the genus Heliophila. We found evidence for loss-of-function changes in genes associated with leaf development and early flowering, and over-retention and sub/neofunctionalization of genes involved in pathogen response and chemical defense. The genomic resources of H. variabilis will help elucidate the role of polyploidization and genome diploidization in plant adaptation to hot arid environments and origin of the Cape flora. The sequenced H. variabilis represents the first chromosome-scale genome assembly of a meso-octoploid representative of the mustard family.

Metabolic and transcriptomic study of pennycress natural variation identifies targets for oil improvement.

Pennycress (Thlaspi arvense L.), a member of the Brassicaceae family, produces seed oil high in erucic acid, suitable for biodiesel and aviation fuel. Although pennycress, a winter annual, could be grown as a dedicated bioenergy crop, an increase in its seed oil content is required to improve its economic competitiveness. The success of crop improvement relies upon finding the right combination of biomarkers and targets, and the best genetic engineering and/or breeding strategies. In this work, we combined biomass composition with metabolomic and transcriptomic studies of developing embryos from 22 pennycress natural variants to identify targets for oil improvement. The selected accession collection presented diverse levels of fatty acids at maturity ranging from 29% to 41%. Pearson correlation analyses, weighted gene co-expression network analysis and biomarker identifications were used as complementary approaches to detect associations between metabolite level or gene expression and oil content at maturity. The results indicated that improving seed oil content can lead to a concomitant increase in the proportion of erucic acid without affecting the weight of embryos. Processes, such as carbon partitioning towards the chloroplast, lipid metabolism, photosynthesis, and a tight control of nitrogen availability, were found to be key for oil improvement in pennycress. Besides identifying specific targets, our results also provide guidance regarding the best timing for their modification, early or middle maturation. Thus, this work lays out promising strategies, specific for pennycress, to accelerate the successful development of lines with increased seed oil content for biofuel applications.

Functional Characterization of Lobularia maritima LmTrxh2 Gene Involved in Cold Tolerance in Tobacco through Alleviation of ROS Damage to the Plasma Membrane.

Cold stress is a key environmental factor affecting plant growth and development, crop productivity, and geographic distribution. Thioredoxins (Trxs) are small proteins that are ubiquitously expressed in all organisms and implicated in several cellular processes, including redox reactions. However, their role in the regulation of cold stress in the halophyte plant Lobularia maritima remains unknown. We recently showed that overexpression of LmTrxh2, which is the gene that encodes the h-type Trx protein previously isolated from L. maritima, led to an enhanced tolerance to salt and osmotic stress in transgenic tobacco. This study functionally characterized the LmTrxh2 gene via its overexpression in tobacco and explored its cold tolerance mechanisms. Results of the RT-qPCR and western blot analyses indicated differential temporal and spatial regulation of LmTrxh2 in L. maritima under cold stress at 4 °C. LmTrxh2 overexpression enhanced the cold tolerance of transgenic tobacco, as evidenced by increased germination rate, fresh weight and catalase (CAT), superoxide dismutase (SOD) and peroxidase (POD) activities; reduced malondialdehyde levels, membrane leakage, superoxide anion (O2-), and hydrogen peroxide (H2O2) levels; and higher retention of chlorophyll than in non-transgenic plants (NT). Furthermore, the transcript levels of reactive oxygen species (ROS)-related genes (NtSOD and NtCAT1), stress-responsive late embryogenis abundant protein 5 (NtLEA5), early response to dehydration 10C (NtERD10C), DRE-binding proteins 1A (NtDREB1A), and cold-responsive (COR) genes (NtCOR15A, NtCOR47, and NtKIN1) were upregulated in transgenic lines compared with those in NT plants under cold stress, indicating that LmTrxh2 conferred cold stress tolerance by enhancing the ROS scavenging ability of plants, thus enabling them to maintain membrane integrity. These results suggest that LmTrxh2 promotes cold tolerance in tobacco and provide new insight into the improvement of cold-stress resistance to cold stress in non-halophyte plants and crops.

Phylogenetic challenges in a recently diversified and polyploid-rich Alyssum (Brassicaceae) lineage: low divergence, reticulation, and parallel polyploid speciation.

Elucidating the evolution of recently diverged and polyploid-rich plant lineages may be challenging even with high-throughput sequencing, both for biological reasons and bioinformatic difficulties. Here, we apply target enrichment with genome skimming (Hyb-Seq) to unravel the evolutionary history of the Alyssum montanum-A. repens species complex. Reconstruction of phylogenetic relationships in diploids supported recent and rapid diversification accompanied by reticulation events. Of the 4 main clades identified among the diploids, 3 clades included species from the Alps, Apennine, and Balkan peninsulas, indicating close biogeographic links between these regions. We further focused on the clade distributed from the Western Alps to the Iberian Peninsula, which comprises numerous polyploids as opposed to a few diploids. Using a recently developed PhyloSD (phylogenomic subgenome detection) pipeline, we successfully tracked the ancestry of all polyploids. We inferred multiple polyploidization events that involved 2 closely related diploid progenitors, resulting into several sibling polyploids: 2 autopolyploids and 6 allopolyploids. The skewed proportions of major homeolog-types and the occurrence of some minor homeolog-types, both exhibiting geographic patterns, suggest introgression with the progenitors and other related diploids. Our study highlights a unique case of parallel polyploid speciation that was enhanced by ecological and geographic separation and provides an excellent resource for future studies of polyploid evolution.

Identification and functional analysis of calcium sensor calmodulins from heavy metal hyperaccumulator Noccaea caerulescens.

Noccaea caerulescens (J. Presl & C. Presl) F. K. Mey. is a heavy metal hyperaccumulator exhibiting extreme tolerance to various environmental stresses. To date, the functional role of Ca2+ -binding protein in this plant is largely unknown. To investigate the function of calmodulins (CaMs) in N. caerulescens , CaM2 , a Ca2+ sensor encoding gene, was identified and functionally characterised. Protein structure analysis showed that NcCaM2 contains four classic exchange factor (EF)-hand motifs with high sequence similarity to the CaM proteins from model plant Arabidopsis thaliana L. Tissue specific expression analysis showed that NcCaM2 is constitutively expressed in stems, leaves, and roots. Expression level of NcCaM2 was significantly upregulated under various environmental stimulus, indicating a potential involvement of NcCaM2 in the tolerance to abiotic stresses. The heterologous expression of NcCaM2 in a yeast mutant strain increased the heavy metal tolerance in yeast cells. Furthermore, the constitutive expression of NcCaM2 enhanced the heavy metal tolerance capability of transgenic tobacco (Nicotiana tabacum L.) plants. Our data suggested an important role of NcCaM2 in the responses to environmental stresses and provided a potential target gene to enhance of the ability to hyperaccumulate metals.

The Gynandropsis gynandra genome provides insights into whole-genome duplications and the evolution of C4 photosynthesis in Cleomaceae.

Gynandropsis gynandra (Cleomaceae) is a cosmopolitan leafy vegetable and medicinal plant, which has also been used as a model to study C4 photosynthesis due to its evolutionary proximity to C3 Arabidopsis (Arabidopsis thaliana). Here, we present the genome sequence of G. gynandra, anchored onto 17 main pseudomolecules with a total length of 740 Mb, an N50 of 42 Mb and 30,933 well-supported gene models. The G. gynandra genome and previously released genomes of C3 relatives in the Cleomaceae and Brassicaceae make an excellent model for studying the role of genome evolution in the transition from C3 to C4 photosynthesis. Our analyses revealed that G. gynandra and its C3 relative Tarenaya hassleriana shared a whole-genome duplication event (Gg-α), then an addition of a third genome (Th-α, +1×) took place in T. hassleriana but not in G. gynandra. Analysis of syntenic copy number of C4 photosynthesis-related gene families indicates that G. gynandra generally retained more duplicated copies of these genes than C3T. hassleriana, and also that the G. gynandra C4 genes might have been under positive selection pressure. Both whole-genome and single-gene duplication were found to contribute to the expansion of the aforementioned gene families in G. gynandra. Collectively, this study enhances our understanding of the polyploidy history, gene duplication and retention, as well as their impact on the evolution of C4 photosynthesis in Cleomaceae.

Stigma receptors control intraspecies and interspecies barriers in Brassicaceae.

Flowering plants have evolved numerous intraspecific and interspecific prezygotic reproductive barriers to prevent production of unfavourable offspring1. Within a species, self-incompatibility (SI) is a widely utilized mechanism that rejects self-pollen2,3 to avoid inbreeding depression. Interspecific barriers restrain breeding between species and often follow the SI × self-compatible (SC) rule, that is, interspecific pollen is unilaterally incompatible (UI) on SI pistils but unilaterally compatible (UC) on SC pistils1,4-6. The molecular mechanisms underlying SI, UI, SC and UC and their interconnections in the Brassicaceae remain unclear. Here we demonstrate that the SI pollen determinant S-locus cysteine-rich protein/S-locus protein 11 (SCR/SP11)2,3 or a signal from UI pollen binds to the SI female determinant S-locus receptor kinase (SRK)2,3, recruits FERONIA (FER)7-9 and activates FER-mediated reactive oxygen species production in SI stigmas10,11 to reject incompatible pollen. For compatible responses, diverged pollen coat protein B-class12-14 from SC and UC pollen differentially trigger nitric oxide, nitrosate FER to suppress reactive oxygen species in SC stigmas to facilitate pollen growth in an intraspecies-preferential manner, maintaining species integrity. Our results show that SRK and FER integrate mechanisms underlying intraspecific and interspecific barriers and offer paths to achieve distant breeding in Brassicaceae crops.

Sequencing of Camelina neglecta, a diploid progenitor of the hexaploid oilseed Camelina sativa.

Camelina neglecta is a diploid species from the genus Camelina, which includes the versatile oilseed Camelina sativa. These species are closely related to Arabidopsis thaliana and the economically important Brassica crop species, making this genus a useful platform to dissect traits of agronomic importance while providing a tool to study the evolution of polyploids. A highly contiguous chromosome-level genome sequence of C. neglecta with an N50 size of 29.1 Mb was generated utilizing Pacific Biosciences (PacBio, Menlo Park, CA) long-read sequencing followed by chromosome conformation phasing. Comparison of the genome with that of C. sativa shows remarkable coincidence with subgenome 1 of the hexaploid, with only one major chromosomal rearrangement separating the two. Synonymous substitution rate analysis of the predicted 34 061 genes suggested subgenome 1 of C. sativa directly descended from C. neglecta around 1.2 mya. Higher functional divergence of genes in the hexaploid as evidenced by the greater number of unique orthogroups, and differential composition of resistant gene analogs, might suggest an immediate adaptation strategy after genome merger. The absence of genome bias in gene fractionation among the subgenomes of C. sativa in comparison with C. neglecta, and the complete lack of fractionation of meiosis-specific genes attests to the neopolyploid status of C. sativa. The assembled genome will provide a tool to further study genome evolution processes in the Camelina genus and potentially allow for the identification and exploitation of novel variation for Camelina crop improvement.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae.

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

The coexpression of two desaturases provides an optimized reduction of saturates in camelina oil.

Reducing the saturate content of vegetable oils is key to increasing their utility and adoption as a feedstock for the production of biofuels. Expression of either the FAT5 16 : 0-CoA desaturase from Caenorhabditis elegans, or an engineered cyanobacterial 16 : 0/18 : 0-glycerolipid desaturase, DES9*, in seeds of Arabidopsis (Arabidopsis thaliana) substantially lowered oil saturates. However, because pathway fluxes and regulation of oil synthesis are known to differ across species, translating this transgene technology from the model plant to crop species requires additional investigation. In the work reported here, we found that high expression of FAT5 in seeds of camelina (Camelina sativa) provided only a moderate decrease in saturates, from 12.9% of total oil fatty acids in untransformed controls to 8.6%. Expression of DES9* reduced saturates to 4.6%, but compromised seed physiology and oil content. However, the coexpression of the two desaturases together cooperatively reduced saturates to only 4.0%, less than one-third of the level in the parental line, without compromising oil yield or seedling germination and establishment. Our successful lowering of oil saturates in camelina identifies strategies that can now be integrated with genetic engineering approaches that reduce polyunsaturates to provide optimized oil composition for biofuels in camelina and other oil seed crops.

The genome of Orychophragmus violaceus provides genomic insights into the evolution of Brassicaceae polyploidization and its distinct traits.

Orychophragmus violaceus, referred to as "eryuelan" (February orchid) in China, is an early-flowering ornamental plant. The high oil content and abundance of unsaturated fatty acids in O. violaceus seeds make it a potential high-quality oilseed crop. Here, we generated a whole-genome assembly for O. violaceus using Nanopore and Hi-C sequencing technologies. The assembled genome of O. violaceus was ∼1.3 Gb in size, with 12 pairs of chromosomes. Through investigation of ancestral genome evolution, we determined that the genome of O. violaceus experienced a tetraploidization event from a diploid progenitor with the translocated proto-Calepineae karyotype. Comparisons between the reconstructed subgenomes of O. violaceus identified indicators of subgenome dominance, indicating that subgenomes likely originated via allotetraploidy. O. violaceus was phylogenetically close to the Brassica genus, and tetraploidy in O. violaceus occurred approximately 8.57 million years ago, close in time to the whole-genome triplication of Brassica that likely arose via an intermediate tetraploid lineage. However, the tetraploidization in Orychophragmus was independent of the hexaploidization in Brassica, as evidenced by the results from detailed phylogenetic analyses and comparisons of the break and fusion points of ancestral genomic blocks. Moreover, identification of multi-copy genes regulating the production of high-quality oil highlighted the contributions of both tetraploidization and tandem duplication to functional innovation in O. violaceus. These findings provide novel insights into the polyploidization evolution of plant species and will promote both functional genomic studies and domestication/breeding efforts in O. violaceus.

Ablation of glucosinolate accumulation in the oil crop Camelina sativa by targeted mutagenesis of genes encoding the transporters GTR1 and GTR2 and regulators of biosynthesis MYB28 and MYB29.

Camelina sativa is an oil crop with low input costs and resistance to abiotic and biotic stresses. The presence of glucosinolates, plant metabolites with adverse health effects, restricts the use of camelina for human and animal nutrition. Cas9 endonuclease-based targeted mutagenesis of the three homeologs of each of the glucosinolate transporters CsGTR1 and CsGTR2 caused a strong decrease in glucosinolate amounts, highlighting the power of this approach for inactivating multiple genes in a hexaploid crop. Mutagenesis of the three homeologs of each of the transcription factors CsMYB28 and CsMYB29 resulted in the complete loss of glucosinolates, representing the first glucosinolate-free Brassicaceae crop. The oil and protein contents and the fatty acid composition of the csgtr1csgtr2 and csmyb28csmyb29 mutant seeds were not affected. The decrease and elimination of glucosinolates improves the quality of the oil and press cake of camelina, which thus complies with international standards regulating glucosinolate levels for human consumption and animal feeding.

Responses of transcriptome and metabolome in the roots of Pugionium cornutum (L.) Gaertn to exogenously applied phthalic acid.

BACKGROUND: The yield and quality of Pugionium cornutum (L.) Gaertn., a healthy, green vegetable with low sugar and high protein contents and high medicinal value, is severely affected by autotoxicity, which is a leading factor in the formation of plant disease. To help characterize the autotoxicity mechanism of P. cornutum (L.) Gaertn., we performed transcriptomic and metabolic analysis of the roots of P. cornutum (L.) Gaertn. response to phthalic acid, an autotoxin from P. cornutum (L.) Gaertn. RESULTS: In this study, high-throughput sequencing of nine RNA-seq libraries generated from the roots.of P. cornutum (L.) Gaertn. under different phthalic acid treatments yielded 37,737 unigenes. In total, 1085 (703 upregulated and 382 downregulated) and 5998 (4385 upregulated and 1613 downregulated) DEGs were identified under 0.1 and 10 mmol·L- 1 phthalic acid treatment, respectively, compared with the control treatment. Glutathione metabolism was among the top five important enriched pathways. In total, 457 and 435 differentially accumulated metabolites were detected under 0.1 and 10 mmol·L- 1 phthalic acid treatment compared with the control, respectively, of which 223 and 253, respectively, increased in abundance. With the increase in phthalic acid concentration, the accumulation of ten metabolites increased significantly, while that of four metabolites decreased significantly, and phthalic acid, dambonitol, 4-hydroxy-butyric acid, homocitrulline, and ethyl ß-D-glucopyranoside were 100 times more abundant under the 10 mmol·L- 1 phthalic acid treatment than under the control. Seventeen differentially expressed genes significantly associated with phthalic acid content were identified. In addition, the L-histidinol content was highest under 0.1 mmol·L- 1 phthalic acid, and a total of eleven differentially expressed genes were significantly positively correlated with the L-histidinol content, all of which were annotated to heat shock proteins, aquaporins and cysteine proteases. CONCLUSIONS: Accumulation of autotoxins altered the metabolic balance in P. cornutum (L.) Gaertn. and influenced water absorption and carbon and nitrogen metabolism. These important results provide insights into the formation mechanisms of autotoxicity and for the subsequent development of new control measures to improve the production and quality of replanted plants.

Hybridization order is not the driving factor behind biases in duplicate gene losses among the hexaploid Solanaceae.

We model the post-hexaploidy evolution of four genomes from the Solanaceae, a group of flowering plants comprising tomatoes, potatoes and their relatives. The hexaploidy that these genomes descend from occurred through two sequential allopolyploidy events and was marked by the unequal losses of duplicated genes from the different progenitor subgenomes. In contrast with the hexaploid Brassiceae (broccoli and its relatives), where the subgenome with the most surviving genes arrived last in the hexaploidy, among the Solanaceae the most preserved subgenome descends from one of the original two tetraploid progenitors. In fact, the last-arriving subgenome in these plants actually has the fewest surviving genes in the modern genomes. We explore whether the distribution of repetitive elements (REs) in these genomes can explain the biases in gene losses, but while the signals we find are broadly consistent with a role for high RE density in driving gene losses, the REs turn over so quickly that little signal of the RE condition at the time of paleopolyploidy is extant in the modern genomes.