Genomic, transcriptomic, and metabolomic analyses provide insights into the evolution and development of a medicinal plant Saposhnikovia divaricata (Apiaceae).

Saposhnikovia divaricata, 2n = 2x = 16, as a perennial species, is widely distributed in China, Mongolia, Russia, etc. It is a traditional Chinese herb used to treat tetanus, rubella pruritus, rheumatic arthralgia, and other diseases. Here, we assembled a 2.07 Gb and N50 scaffold length of 227.67 Mb high-quality chromosome-level genome of S. divaricata based on the PacBio Sequel II sequencing platform. The total number of genes identified was 42 948, and 42 456 of them were functionally annotated. A total of 85.07% of the genome was composed of repeat sequences, comprised mainly of long terminal repeats (LTRs) which represented 73.7% of the genome sequence. The genome size may have been affected by a recent whole-genome duplication event. Transcriptional and metabolic analyses revealed bolting and non-bolting S. divaricata differed in flavonoids, plant hormones, and some pharmacologically active components. The analysis of its genome, transcriptome, and metabolome helped to provide insights into the evolution of bolting and non-bolting phenotypes in wild and cultivated S. divaricata and lays the basis for genetic improvement of the species.

Advancement of Research Progress on Synthesis Mechanism of Cannabidiol (CBD).

Cannabis sativa L. is a multipurpose crop with high value for food, textiles, and other industries. Its secondary metabolites, including cannabidiol (CBD), have potential for broad application in medicine. With the CBD market expanding, traditional production may not be sufficient. Here we review the potential for the production of CBD using biotechnology. We describe the chemical and biological synthesis of cannabinoids, the associated enzymes, and the application of metabolic engineering, synthetic biology, and heterologous expression to increasing production of CBD.

Soil and seed both influence bacterial diversity in the microbiome of the Cannabis sativa seedling endosphere.

Introduction: Phytobiomes have a significant impact on plant health. The microbiome of Cannabis sativa is particularly interesting both because of renewed interest in this crop and because it is commercially propagated in two different ways (i.e. clonally and by seed). Angiosperms obtain a founding population of seed-borne endophytes from their seed-bearing parent. This study examines the influence of both seed and soil-derived bacteria on the endospheres of cannabis seedlings of both hemp- and drug-types. Methods: A multi-factorial metagenomic study was conducted with three cannabis genotypes and two soil sources, which were tested both before and after autoclave sterilization. Seedlings were grown on soil, then rinsed and surface-sterilized, and 16S rDNA amplicons from seedling endophytes were sequenced, taxonomically classified, and used to estimate alpha- and beta-diversity in Qiime2. The statistical significance of differences in seedling microbiomes across treatments was tested, and PiCRUST2 was used to infer the functional relevance of these differences. Results: Soil was found to have a profound effect on the alpha-diversity, beta-diversity, relative abundance, and functional genes of endophytic bacteria in germinating cannabis seedlings. Additionally, there was a significant effect of cannabis genotype on beta diversity, especially when genotypes were grown in sterilized soil. Gammaproteobacteria and Bacilli were the two most abundant taxa and were found in all genotypes and soil types, including sterilized soil. Discussion: The results indicated that a component of cannabis seedling endosphere microbiomes is seed-derived and conserved across the environments tested. Functional prediction of seedling endophytes using piCRUST suggested a number of important functions of seed-borne endophytes in cannabis including nutrient and amino acid cycling, hormone regulation, and as precursors to antibiotics. This study suggested both seed and soil play a critical role in shaping the microbiome of germinating cannabis seedlings.

WRKY6 transcription factor modulates root potassium acquisition through promoting expression of AKT1 in Arabidopsis.

Potassium (K+), being an essential macronutrient in plants, plays a central role in many aspects. Root growth is highly plastic and is affected by many different abiotic stresses including nutrient deficiency. The Shaker-type K+ channel Arabidopsis (Arabidopsis thaliana) K+ Transporter 1 (AKT1) is responsible for K+ uptake under both low and high external K+ conditions. However, the upstream transcription factor of AKT1 is not clear. Here, we demonstrated that the WRKY6 transcription factor modulates root growth to low potassium (LK) stress in Arabidopsis. WRKY6 showed a quick response to LK stress and also to many other abiotic stress treatments. The two wrky6 T-DNA insertion mutants were highly sensitive to LK treatment, whose primary root lengths were much shorter, less biomass and lower K+ content in roots than those of wild-type plants, while WRKY6-overexpression lines showed opposite phenotypes. A further investigation showed that WRKY6 regulated the expression of the AKT1 gene via directly binding to the W-box elements in its promoter through EMSA and ChIP-qPCR assays. A dual luciferase reporter analysis further demonstrated that WRKY6 enhanced the transcription of AKT1. Genetic analysis further revealed that the overexpression of AKT1 greatly rescued the short root phenotype of the wrky6 mutant under LK stress, suggesting AKT1 is epistatic to WRKY6 in the control of LK response. Further transcriptome profiling suggested that WRKY6 modulates LK response through a complex regulatory network. Thus, this study unveils a transcription factor that modulates root growth under potassium deficiency conditions by affecting the potassium channel gene AKT1 expression.

Genome-Wide Analysis of Flax (Linum usitatissimum L.) Growth-Regulating Factor (GRF) Transcription Factors.

Flax is an important cash crop globally with a variety of commercial uses. It has been widely used for fiber, oil, nutrition, feed and in composite materials. Growth regulatory factor (GRF) is a transcription factor family unique to plants, and is involved in regulating many processes of growth and development. Bioinformatics analysis of the GRF family in flax predicted 17 LuGRF genes, which all contained the characteristic QLQ and WRC domains. Equally, 15 of 17 LuGRFs (88%) are predicted to be regulated by lus-miR396 miRNA. Phylogenetic analysis of GRFs from flax and several other well-characterized species defined five clades; LuGRF genes were found in four clades. Most LuGRF gene promoters contained cis-regulatory elements known to be responsive to hormones and stress. The chromosomal locations and collinearity of LuGRF genes were also analyzed. The three-dimensional structure of LuGRF proteins was predicted using homology modeling. The transcript expression data indicated that most LuGRF family members were highly expressed in flax fruit and embryos, whereas LuGRF3, LuGRF12 and LuGRF16 were enriched in response to salt stress. Real-time quantitative fluorescent PCR (qRT-PCR) showed that both LuGRF1 and LuGRF11 were up-regulated under ABA and MeJA stimuli, indicating that these genes were involved in defense. LuGRF1 was demonstrated to be localized to the nucleus as expected for a transcription factor. These results provide a basis for further exploration of the molecular mechanism of LuGRF gene function and obtaining improved flax breeding lines.

Evaluation of Differentially Expressed Genes in Leaves vs. Roots Subjected to Drought Stress in Flax (Linum usitatissimum L.).

Drought stress is a common environmental challenge that plants face, severely constraining plant growth and reducing crop yield and quality. Several studies have highlighted distinct responses between monocotyledonous and dicotyledonous plants. However, the mechanisms underlying flax tolerance to abiotic stress, such as drought, remain unclear. In this study, we investigated the morphological, physiological, and biochemical characteristics and the genome-wide gene expression of oil flax and fiber flax in response to drought stress. The results revealed that drought stress caused significant wilting of flax leaves. Within the first 24 h of stress, various physiological and biochemical characteristics exhibited rapid responses. These included fresh weight, relative water content (RWC), proline, soluble protein, soluble sugar, superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) in the leaves or roots of flax. Additionally, drought stress led to a significant rise in lignin content in fiber flax. In addition, the transcriptome analysis demonstrated genome-wide variations in gene expression induced by drought stress. Specifically, genes associated with photosynthesis, proline biosynthesis, and phytohormone metabolism exhibited significant differences in expression levels under stress conditions in flax. These findings highlight the rapid response of flax to drought stress within a short-term period. Our experiment also revealed that, although there were variations in the levels of small compound content or gene expression between Longya10 and Fany under drought stress, most stress-resistance responses were similar. Furthermore, the results provide additional evidence supporting the existence of mechanisms underlying the response to drought stress in plants.

Hybrid Genome Assembly of Berkeleyomyces rouxiae, an Emerging Cannabis Fungal Pathogen Causing Black Root Rot in an Aeroponic Facility.

The resurged interest in cultivation of Cannabis sativa has presented an array of new challenges. Among them are the difficult-to-control pests and pathogens that infect cannabis plants. The limited methods for disease control available to cannabis growers necessitates early detection of plant pathogens, something that molecular techniques such as DNA sequencing has greatly improved. This study reports for the first time the fungal plant pathogen Berkeleyomyces rouxiae causing black root rot in high THC-containing cannabis. Aeroponically grown cannabis plants at a licenced production facility in Cranbrook BC, Canada, rapidly displayed root discoloration and rot symptoms despite testing negative for all commercially available pathogen tests. Developing sequencing-based disease diagnostics requires genomic information, so this study presents the first whole genome sequence of the multihost, widespread black root rot pathogen B. rouxiae. Hybrid genome assembly using Oxford Nanopore long-reads and Illumina short-reads yielded a genome size of 28.2 Mb represented over 404 contigs with an N50 of 267 kb. Genome annotation predicted 6,960 protein-coding genes with 59,477 functional annotations. The availability of this genome will assist in sequence-based diagnostic development, comparative genomics, and taxonomic resolution of this globally important plant pathogen.

Cannabis Seedlings Inherit Seed-Borne Bioactive and Anti-Fungal Endophytic Bacilli.

Throughout the hundreds of millions of years of co-evolution, plants and microorganisms have established intricate symbiotic and pathogenic relationships. Microbial communities associated with plants are in constant flux and can ultimately determine whether a plant will successfully reproduce or be destroyed by their environment. Inheritance of beneficial microorganisms is an adaptation plants can use to protect germinating seeds against biotic and abiotic stresses as seedlings develop. The interest in Cannabis as a modern crop requires research into effective biocontrol of common fungal pathogens, an area that has seen little research. This study examines the seed-borne endophytes present across 15 accessions of Cannabis grown to seed across Western Canada. Both hemp and marijuana seedlings inherited a closely related group of bioactive endophytic Bacilli. All Cannabis accessions possessed seed-inherited Paenibacillus mobilis with the capacity to solubilize mineral phosphate. Additionally, seeds were found to carry genera of fungal isolates known to be Cannabis pathogens and post-harvest molds: Alternaria, Penicillium, Cladosporium, Chaetomium, Aspergillus, Rhizopus, and Fusarium. Thirteen seed-borne endophytes showed antibiotic activity against Alternaria, Aspergillus, Penicillium, and Fusarium. This study suggests both fungal pathogens and bacterial endophytes that antagonize them are vectored across generations in Cannabis as they compete over this shared niche.

Mammalian Melatonin Agonist Pharmaceuticals Stimulate Rhomboid Proteins in Plants.

Melatonin is a human neurotransmitter and plant signalling metabolite that perceives and directs plant metabolism. The mechanisms of melatonin action in plants remain undefined. We hypothesized that roots have a melatonin-specific receptor and/or transporter that can respond to melatonin-mediating pharmaceuticals. To test this hypothesis Arabidopsis seedlings were grown with melatonin pharmaceutical receptor agonists: ramelteon and tasimelteon, and/or antagonists: luzindole and 4-P-PDOT. Ramelteon was found both to mimic and competitively inhibit melatonin metabolism in plants. Due to the higher selectivity of ramelteon for the MT1 receptor type in humans, a sequence homology search for MT1 in Arabidopsis identified the rhomboid-like protein 7 (RBL7). In physiological studies, Arabidopsis rbl7 mutants were less responsive to ramelteon and melatonin. Quantum dot visualizations of the effects of ramelteon on melatonin binding to root cell membranes revealed a potential mechanism. We propose that RBL7 is a melatonin-interacting protein that directs root architecture and growth in a mechanism that is responsive to environmental factors.

An Overview of Databases and Bioinformatics Tools for Plant Antimicrobial Peptides.

Antimicrobial Peptides (AMPs) are small, ribosomally synthesized proteins found in nearly all forms of life. In plants, AMPs play a central role in plant defense due to their distinct physicochemical properties. Due to their broad-spectrum antimicrobial activity and rapid killing action, plant AMPs have become important candidates for the development of new drugs to control plant and animal pathogens that are resistant to multiple drugs. Further research is required to explore the potential uses of these natural compounds. Computational strategies have been increasingly used to understand key aspects of antimicrobial peptides. These strategies will help to minimize the time and cost of "wet-lab" experimentation. Researchers have developed various tools and databases to provide updated information on AMPs. However, despite the increased availability of antimicrobial peptide resources in biological databases, finding AMPs from plants can still be a difficult task. The number of plant AMP sequences in current databases is still small and yet often redundant. To facilitate further characterization of plant AMPs, we have summarized information on the location, distribution, and annotations of plant AMPs available in the most relevant databases for AMPs research. We also mapped and categorized the bioinformatics tools available in these databases. We expect that this will allow researchers to advance in the discovery and development of new plant AMPs with potent biological properties. We hope to provide insights to further expand the application of AMPs in the fields of biotechnology, pharmacy, and agriculture.

Ectopic overexpression of a membrane-tethered transcription factor gene NAC60 from oilseed rape positively modulates programmed cell death and age-triggered leaf senescence.

Senescence is an integrative final stage of plant development that is governed by internal and external cues. The NAM, ATAF1/2, CUC2 (NAC) transcription factor (TF) family is specific to plants and membrane-tethered NAC TFs (MTTFs) constitute a unique and sophisticated mechanism in stress responses and development. However, the function of MTTFs in oilseed rape (Brassica napus L.) remains unknown. Here, we report that BnaNAC60 is an MTTF associated with the endoplasmic reticulum (ER) membrane. Expression of BnaNAC60 was induced during the progression of leaf senescence. Translocation of BnaNAC60 into nuclei was induced by ER stress and oxidative stress treatments. It binds to the NTLBS motif, rather than the canonical NAC recognition site. Overexpression of BnaNAC60 devoid of the transmembrane domain, but not the full-length BnaNAC60, induces significant reactive oxygen species (ROS) accumulation and hypersensitive response-like cell death in both tobacco (Nicotiana benthamiana) and oilseed rape protoplasts. Moreover, ectopic overexpression of BnaNAC60 devoid of the transmembrane domain, but not the full-length BnaNAC60, in Arabidopsis also induces precocious leaf senescence. Furthermore, screening and expression profiling identified an array of functional genes that are significantly induced by BnaNAC60 expression. Further it was found that BnaNAC60 can activate the promoter activities of BnaNYC1, BnaRbohD, BnaBFN1, BnaZAT12, and multiple BnaVPEs in a dual-luciferase reporter assay. Electrophoretic mobility shift assay and chromatin immunoprecipitation coupled to quantitative PCR assays revealed that BnaNAC60 directly binds to the promoter regions of these downstream target genes. To summarize, our data show that BnaNAC60 is an MTTF that modulates cell death, ROS accumulation, and leaf senescence.

WRKY42 transcription factor positively regulates leaf senescence through modulating SA and ROS synthesis in Arabidopsis thaliana.

Leaf senescence represents the final stage of leaf growth and development, and its onset and progression are strictly regulated; however, the underlying regulatory mechanisms remain largely unknown. In this study we found that WRKY42 was highly induced during leaf senescence. Loss-of-function wrky42 mutants showed delayed leaf senescence whereas the overexpression of WRKY42 accelerated senescence. Transcriptome analysis revealed 2721 differentially expressed genes between wild-type and WRKY42-overexpressing plants, including genes involved in salicylic acid (SA) and reactive oxygen species (ROS) synthesis as well as several senescence-associated genes (SAGs). Moreover, WRKY42 activated the transcription of isochorismate synthase 1 (ICS1), respiratory burst oxidase homolog F (RbohF) and a few SAG genes. Consistently, the expression of these genes was reduced in wrky42 mutants but was markedly increased in transgenic Arabidopsis overexpressing WRKY42. Both in vitro electrophoretic mobility shift assays (EMSAs) and in vivo chromatin immunoprecipitation and dual luciferase assays demonstrated that WRKY42 directly bound to the promoters of ICS1 and RbohF, as well as a few SAGs, to activate their expression. Genetic analysis further showed that mutations of ICS1 and RbohF suppressed the early senescence phenotype evoked by WRKY42 overexpression. Thus, we have identified WRKY42 as a novel transcription factor positively regulating leaf senescence by directly activating the transcription of ICS1, RbohF and SAGs, without any seed yield penalty.

A Rapeseed WRKY Transcription Factor Phosphorylated by CPK Modulates Cell Death and Leaf Senescence by Regulating the Expression of ROS and SA-Synthesis-Related Genes.

Salicylic acid (SA) and reactive oxygen species (ROS) are two well-defined inducers of leaf senescence. Here, we identified a novel WRKY transcription factor gene WSR1 (WRKY regulating SA and ROS 1) in Brassica napus (rapeseed) in promoting SA and ROS production, which eventually led to leaf senescence thereafter. Its expression increased in senescing leaves. Ca2+-dependent protein kinase (CPK) 5 and -6 interacted with and phosphorylated BnaWSR1. Overexpression of phosphomimic BnaWSR1 (BnaWSR1ca) in rapeseed protoplasts elicited ROS production and cell death while its ectopic expression in Arabidopsis enhanced SA and ROS levels and, hence, accelerated leaf senescence. Furthermore, BnaWSR1ca activated the expression of Isochorismate Synthase 1 (ICS1), Respiratory Burst Oxidase Homologue (Rboh) D, and Senescence-Associated Gene 14 (SAG14). Electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR) assays demonstrated that BnaWSR1ca directly bound to promoter regions of ICS1, RbohD, and SAG14. These data have identified a CPK-WSR1 module that integrates SA and ROS to control cell death and leaf senescence.

A single nucleotide polymorphism assay sheds light on the extent and distribution of genetic diversity, population structure and functional basis of key traits in cultivated north American cannabis.

BACKGROUND: The taxonomic classification of Cannabis genus has been delineated through three main types: sativa (tall and less branched plant with long and narrow leaves), indica (short and highly branched plant with broader leaves) and ruderalis (heirloom type with short stature, less branching and small thick leaves). While still under discussion, particularly whether the genus is polytypic or monotypic, this broad classification reflects putative geographical origins of each group and putative chemotype and pharmacologic effect. METHODS: Here we describe a thorough investigation of cannabis accessions using a set of 23 highly informative and polymorphic SNP (Single Nucleotide Polymorphism) markers associated with important traits such as cannabinoid and terpenoid expression as well as fibre and resin production. The assay offers insight into cannabis population structure, phylogenetic relationship, population genetics and correlation to secondary metabolite concentrations. We demonstrate the utility of the assay for rapid, repeatable and cost-efficient genotyping of commercial and industrial cannabis accessions for use in product traceability, breeding programs, regulatory compliance and consumer education. RESULTS: We identified 5 clusters in the sample set, including industrial hemp (K5) and resin hemp, which likely underwent a bottleneck to stabilize cannabidiolic acid (CBDA) accumulation (K2, Type II & III). Tetrahydrocannabinolic acid (THCA) resin (Type I) makes up the other three clusters with terpinolene (K4 - colloquial "sativa" or "Narrow Leaflet Drug" (NLD), myrcene/pinene (K1) and myrcene/limonene/linalool (K3 - colloquial "indica", "Broad Leaflet Drug" (BLD), which also putatively harbour an active version of the cannabichrometic acid Synthase gene (CBCAS). CONCLUSION: The final chemical compositions of cannabis products have key traits related to their genetic identities. Our analyses in the context of the NCBI Cannabis sativa Annotation Release 100 allows for hypothesis testing with regards to secondary metabolite production. Genetic markers related to secondary metabolite production will be important in many sectors of the cannabis marketplace. For example, markers related to THC production will be important for adaptable and compliant large-scale seed production under the new US Domestic Hemp Production Program.

Organellomic data sets confirm a cryptic consensus on (unrooted) land-plant relationships and provide new insights into bryophyte molecular evolution.

PREMISE: Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution. METHODS: We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements. RESULTS: Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses. CONCLUSIONS: A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure.

A virus-induced gene-silencing system for functional genetics in a betalainic species, Amaranthus tricolor (Amaranthaceae).

PREMISE OF THE STUDY: Research in Amaranthaceae could be accelerated by developing methods for targeted gene silencing. Most amaranths, including Amaranthus tricolor, produce betalains. However, the physiological and ecological roles of these pigments are uncertain. We sought to establish a virus-induced gene-silencing (VIGS) method for amaranths, using silencing of betalain pigments as a proof-of-principle. METHODS: We targeted AtriCYP76AD1, a putative cytochrome P450 component of the betalain biosynthetic pathway, using VIGS, and compared two different methods of introducing the VIGS construct into plants. We measured transcript abundance and concentrations of betalains and their l-DOPA precursor in VIGS-treated plants, and compared these to controls. RESULTS: We observed that when AtriCYP76AD1 was targeted by VIGS in normally red plants, AtriCYP76AD1 and the related genes AtriCYP76AD6 and AtriCYP76AD5 had diminished transcript abundance. Furthermore, newly emergent petioles and leaves of VIGS-treated plants appeared green, betacyanin accumulation was strongly reduced, and l-DOPA accumulation was increased. No betaxanthin could be detected in this variety of A. tricolor, either before or after VIGS treatment. DISCUSSION: These results help to establish the genetic basis of betalain synthesis in amaranths. Furthermore, this is the first report of VIGS in amaranths and demonstrates the potential of this technique for basic and applied research in these species.

VvWRKY30, a grape WRKY transcription factor, plays a positive regulatory role under salinity stress.

High salinity severely inhibits the growth and productivity of grape plants. However, knowledge of salt-stress regulation remains limited in WRKY members of grapes. Here, we isolated a novel VvWRKY30 gene from Vitis vinifera L. and studied its role in salt-stress resistance. The VvWRKY30 protein fused with green fluorescent protein localized to the nucleus and the transcriptional activation activity of VvWRKY30 was confirmed in yeast. Moreover, VvWRKY30 showed key transcriptional activity domain at the N-terminal and specifically binds to the W-BOX. VvWRKY30 showed the highest expression in the shoot tip and functional leaves of grape plants. VvWRKY30 expression was induced by salt as well as stress signaling molecules H2S and H2O2. Overexpression of VvWRKY30 in Arabidopsis increased its resistance to salt stress at different stages of growth. Under salinity stress, VvWRKY30 overexpressing lines had higher antioxidant activities and lower reactive oxygen species contents. Soluble sugar and proline concentrations also increased in VvWRKY30 overexpressing lines in the presence of NaCl. In addition, the transcription of genes related to antioxidant biosynthesis, glyco-metabolism and proline biosynthesis increased in the VvWRKY30 overexpressing lines. Taken together, this study confirmed the important role of VvWRKY30 in increasing salt stress resistance by regulating reactive oxygen species-scavenging and the accumulation of osmoticum.

Chromosome-scale pseudomolecules refined by optical, physical and genetic maps in flax.

Genomes of varying sizes have been sequenced with next-generation sequencing platforms. However, most reference sequences include draft unordered scaffolds containing chimeras caused by mis-scaffolding. A BioNano genome (BNG) optical map was constructed to improve the previously sequenced flax genome (Linum usitatissimum L., 2n = 30, about 373 Mb), which consisted of 3852 scaffolds larger than 1 kb and totalling 300.6 Mb. The high-resolution BNG map of cv. CDC Bethune totalled 317 Mb and consisted of 251 BNG contigs with an N50 of 2.15 Mb. A total of 622 scaffolds (286.6 Mb, 94.9%) aligned to 211 BNG contigs (298.6 Mb, 94.2%). Of those, 99 scaffolds, diagnosed to contain assembly errors, were refined into 225 new scaffolds. Using the newly refined scaffold sequences and the validated bacterial artificial chromosome-based physical map of CDC Bethune, the 211 BNG contigs were scaffolded into 94 super-BNG contigs (N50 of 6.64 Mb) that were further assigned to the 15 flax chromosomes using the genetic map. The pseudomolecules total about 316 Mb, with individual chromosomes of 15.6 to 29.4 Mb, and cover 97% of the annotated genes. Evidence from the chromosome-scale pseudomolecules suggests that flax has undergone palaeopolyploidization and mesopolyploidization events, followed by rearrangements and deletions or fusion of chromosome arms from an ancient progenitor with a haploid chromosome number of eight.

A Novel NAC-Type Transcription Factor, NAC87, from Oilseed Rape Modulates Reactive Oxygen Species Accumulation and Cell Death.

Reactive oxygen species (ROS) are thought to play a dual role in plants by functioning as signaling molecules and toxic by-products of aerobic metabolism. The hypersensitive response (HR) is a typical feature of immune responses in plants and also a type of programmed cell death (PCD). How these two processes are regulated in oilseed rape (Brassica napus L.) at the transcriptional level remains largely unknown. In this study, we report that an oilseed rape (Brassica napus L.) NAM-ATAF-CUC (NAC)-type transcription factor NAC87 modulates ROS and cell death accompanied by typical changes at the morphological and cellular levels. The BnaNAC87 gene was induced by multiple stress and hormone treatments and was highly expressed in senescent leaves by quantitative reverse transcription-PCR (qRT-PCR). BnaNAC87 is located in nuclei and has transcriptional activation activity. Expression of BnaNAC87 promoted significant ROS production, cell death as well as death of protoplasts, as indicated by histological staining. In addition, putative downstream target genes of NAC87 were identified through both qRT-PCR and dual luciferase reporter assays. We found that genes implicated in ROS generation (RbohB), cell death (VPE1a, ZEN1), leaf senescence (WRKY6, ZAT12) and defense (PR2, PR5 and HIN1) were significantly induced. Through an electrophoretic mobility shift assay (EMSA), we confirmed that BnaNAC87 directly binds to the NACRS-containing promoter fragments of ZEN1, ZAT12, HIN1 and PR5 genes. From these results, we conclude that oilseed rape NAC87 is a novel NAC transcription factor that acts as a positive regulator of ROS metabolism and cell death.

Induced Mutagenesis in UGT74S1 Gene Leads to Stable New Flax Lines with Altered Secoisolariciresinol Diglucoside (SDG) Profiles.

Flax secoisolariciresinol (SECO) diglucoside (SDG) lignan is an emerging natural product purported to prevent chronic diseases in humans. SECO, the aglycone form of SDG, has shown higher intestinal cell absorption but it is not accumulated naturally in planta. Recently, we have identified and characterized a UDP-glucosyltransferase gene, UGT74S1, that glucosylates SECO into its monoglucoside (SMG) and SDG forms when expressed in yeast. However, whether this gene is unique in controlling SECO glucosylation into SDG in planta is unclear. Here, we report on the use of UGT74S1 in reverse and forward genetics to characterize an ethyl methane sulfonate (EMS) mutagenized flax population from cultivar CDC Bethune and consisting of 1996 M2 families. EMS mutagenesis generated 73 SNP variants causing 79 mutational events in the UGT74S1 exonic regions of 93 M2 families. The mutation frequency in the exonic regions was determined to be one per 28 Kb. Of these mutations, 13 homozygous missense mutations and two homozygous nonsense mutations were observed and all were transmitted into the M3 and M4 generations. Forward genetics screening of the population showed homozygous nonsense mutants completely lacking SDG biosynthesis while the production of SMG was observed only in a subset of the M4 lines. Heterozygous or homozygous M4 missense mutants displayed a wide range of SDG levels, some being greater than those of CDC Bethune. No additional deleterious mutations were detected in these mutant lines using a panel of 10 other genes potentially involved in the lignan biosynthesis. This study provides further evidence that UGT74S1 is unique in controlling SDG formation from SECO and this is the first report of non-transgenic flax germplasm with simultaneous knockout of SDG and presence of SMG in planta.