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.

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.

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.

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.

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.

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.

Orbitides and free polyamines have similarly limited fungicidal activity against three common pathogens of flax in vitro

Fusarium oxysporum f. sp. lini and Septoria linicola are causes of fusarium wilt and pasmo in flax ( Linum usitatissimum). Members of a third fungal genus, Alternaria spp., have also been found in fiber and linseed varieties of flax, and are a source of post-harvest spoilage and mycotoxins in a wide range of crops. We performed a microdilution assay and calculated the median effective concentration (EC50) to compare the potency of cyclolinopeptides (CLPs), two polyamines (spermidine and spermine), and the fungicide carbendazimin in the control of three fungi that have potential pathogenic activity ( F. oxysporum, S. linicola, and Alternaria spp), of which the first two are particulary significant causes of disease in flax. For carbendazim, all EC50 values were <0.6 μg/mL. The observed EC50 ranged from 111 to 340 μg/mL for a mixture of six unique CLPs, 109 to 778 μg/mL for spermine, and 21 to 272 μg/mL for spermidine. Spermidine was most effective against Alternaria sp., with an EC50 of 21 μg/mL. The results presented here showed that polyamines and CLPs possess limited antifungal activities against several fungi, with spermidines the most effective naturally occurring compound tested. Our findings do not support the hypothesis that CLPs act as potent antifungals against the three species of pathogens tested.

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.

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.

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.

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. Ca²⁺-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.

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.

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 H₂S and H₂O₂. 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.

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.

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.

Insights into the Evolution of Hydroxyproline-Rich Glycoproteins from 1000 Plant Transcriptomes

The carbohydrate-rich cell walls of land plants and algae have been the focus of much interest given the value of cell wall-based products to our current and future economies. Hydroxyproline-rich glycoproteins (HRGPs), a major group of wall glycoproteins, play important roles in plant growth and development, yet little is known about how they have evolved in parallel with the polysaccharide components of walls. We investigate the origins and evolution of the HRGP superfamily, which is commonly divided into three major multigene families: the arabinogalactan proteins (AGPs), extensins (EXTs), and proline-rich proteins. Using motif and amino acid bias, a newly developed bioinformatics pipeline, we identified HRGPs in sequences from the 1000 Plants transcriptome project (www.onekp.com). Our analyses provide new insights into the evolution of HRGPs across major evolutionary milestones, including the transition to land and the early radiation of angiosperms. Significantly, data mining reveals the origin of glycosylphosphatidylinositol (GPI)-anchored AGPs in green algae and a 3- to 4-fold increase in GPI-AGPs in liverworts and mosses. The first detection of cross-linking (CL)-EXTs is observed in bryophytes, which suggests that CL-EXTs arose though the juxtaposition of preexisting SPn EXT glycomotifs with refined Y-based motifs. We also detected the loss of CL-EXT in a few lineages, including the grass family (Poaceae), that have a cell wall composition distinct from other monocots and eudicots. A key challenge in HRGP research is tracking individual HRGPs throughout evolution. Using the 1000 Plants output, we were able to find putative orthologs of Arabidopsis pollen-specific GPI-AGPs in basal eudicots.

The Evolution of HD2 Proteins in Green Plants

In eukaryotes, protein deacetylation is carried out by two well-conserved histone deacetylase (HDAC) families: RPD3/HDA1 and SIR2. Intriguingly, model plants such as Arabidopsis express an additional plant-specific HDAC family, termed type-2 HDACs (HD2s). Transcriptomic analyses from more than 1300 green plants generated by the 1000 plants (1KP) consortium showed that HD2s appeared early in green plant evolution, the first members being detected in several streptophyte green alga. The HD2 family has expanded via several rounds of successive duplication; members are expressed in all major green plant clades. Interestingly, angiosperm species express new HD2 genes devoid of a zinc-finger domain, one of the main structural features of HD2s. These variants may have been associated with the origin and/or the biology of the ovule/seed.

Functional characterization of NAC55 transcription factor from oilseed rape (Brassica napus L.) as a novel transcriptional activator modulating reactive oxygen species accumulation and cell death

NAC transcription factors (TFs) are plant-specific and play important roles in development, responses to biotic and abiotic cues and hormone signaling. So far, only a few NAC genes have been reported to regulate cell death. In this study, we identified and characterized a NAC55 gene isolated from oilseed rape (Brassica napus L.). BnaNAC55 responds to multiple stresses, including cold, heat, abscisic acid (ABA), jasmonic acid (JA) and a necrotrophic fungal pathogen Sclerotinia sclerotiorum. BnaNAC55 has transactivation activity and is located in the nucleus. BnaNAC55 is able to form homodimers in planta. Unlike ANAC055, full-length BnaNAC55, but not either the N-terminal NAC domain or C-terminal regulatory domain, induces ROS accumulation and hypersensitive response (HR)-like cell death when expressed both in oilseed rape protoplasts and Nicotiana benthamiana. Furthermore, BnaNAC55 expression causes obvious nuclear DNA fragmentation. Moreover, quantitative reverse transcription PCR (qRT-PCR) analysis identified that the expression levels of multiple genes regulating ROS production and scavenging, defense response as well as senescence are significantly induced. Using a dual luciferase reporter assay, we further confirm that BnaNAC55 could activate the expression of a few ROS and defense-related gene expression. Taken together, our work has identified a novel NAC TF from oilseed rape that modulates ROS accumulation and cell death.

The Amaranth Genome: Genome, Transcriptome, and Physical Map Assembly

Amaranth ( L.) is an emerging pseudocereal native to the New World that has garnered increased attention in recent years because of its nutritional quality, in particular its seed protein and more specifically its high levels of the essential amino acid lysine. It belongs to the Amaranthaceae family, is an ancient paleopolyploid that shows disomic inheritance (2 = 32), and has an estimated genome size of 466 Mb. Here we present a high-quality draft genome sequence of the grain amaranth. The genome assembly consisted of 377 Mb in 3518 scaffolds with an N of 371 kb. Repetitive element analysis predicted that 48% of the genome is comprised of repeat sequences, of which -like elements were the most commonly classified retrotransposon. A de novo transcriptome consisting of 66,370 contigs was assembled from eight different amaranth tissue and abiotic stress libraries. Annotation of the genome identified 23,059 protein-coding genes. Seven grain amaranths (, , and ) and their putative progenitor () were resequenced. A single nucleotide polymorphism (SNP) phylogeny supported the classification of as the progenitor species of the grain amaranths. Lastly, we generated a de novo physical map for using the BioNano Genomics' Genome Mapping platform. The physical map spanned 340 Mb and a hybrid assembly using the BioNano physical maps nearly doubled the N of the assembly to 697 kb. Moreover, we analyzed synteny between amaranth and sugar beet ( L.) and estimated, using analysis, the age of the most recent polyploidization event in amaranth.

RNA-seq Transcriptome Response of Flax (Linum usitatissimum L.) to the Pathogenic Fungus Fusarium oxysporum f. sp. lini

Fusarium oxysporum f. sp. lini is a hemibiotrophic fungus that causes wilt in flax. Along with rust, fusarium wilt has become an important factor in flax production worldwide. Resistant flax cultivars have been used to manage the disease, but the resistance varies, depending on the interactions between specific cultivars and isolates of the pathogen. This interaction has a strong molecular basis, but no genomic information is available on how the plant responds to attempted infection, to inform breeding programs on potential candidate genes to evaluate or improve resistance across cultivars. In the current study, disease progression in two flax cultivars [Crop Development Center (CDC) Bethune and Lutea], showed earlier disease symptoms and higher susceptibility in the later cultivar. Chitinase gene expression was also divergent and demonstrated and earlier molecular response in Lutea. The most resistant cultivar (CDC Bethune) was used for a full RNA-seq transcriptome study through a time course at 2, 4, 8, and 18 days post-inoculation (DPI). While over 100 genes were significantly differentially expressed at both 4 and 8 DPI, the broadest deployment of plant defense responses was evident at 18 DPI with transcripts of more than 1,000 genes responding to the treatment. These genes evidenced a reception and transduction of pathogen signals, a large transcriptional reprogramming, induction of hormone signaling, activation of pathogenesis-related genes, and changes in secondary metabolism. Among these, several key genes that consistently appear in studies of plant-pathogen interactions, had increased transcript abundance in our study, and constitute suitable candidates for resistance breeding programs. These included: an induced RPMI-induced protein kinase; transcription factors WRKY3, WRKY70, WRKY75, MYB113, and MYB108; the ethylene response factors ERF1 and ERF14; two genes involved in auxin/glucosinolate precursor synthesis (CYP79B2 and CYP79B3); the flavonoid-related enzymes chalcone synthase, dihydroflavonol reductase and multiple anthocyanidin synthases; and a peroxidase implicated in lignin formation (PRX52). Additionally, regulation of some genes indicated potential pathogen manipulation to facilitate infection; these included four disease resistance proteins that were repressed, indole acetic acid amido/amino hydrolases which were upregulated, activated expansins and glucanases, amino acid transporters and aquaporins, and finally, repression of major latex proteins.

Wild soybean roots depend on specific transcription factors and oxidation reduction related genesin response to alkaline stress

Soil alkalinity is an important environmental problem limiting agricultural productivity. Wild soybean (Glycine soja) shows strong alkaline stress tolerance, so it is an ideal plant candidate for studying the molecular mechanisms of alkaline tolerance and identifying alkaline stress-responsive genes. However, limited information is available about G. soja responses to alkaline stress on a genomic scale. Therefore, in the present study, we used RNA sequencing to compare transcript profiles of G. soja root responses to sodium bicarbonate (NaHCO3) at six time points, and a total of 68,138,478 pairs of clean reads were obtained using the Illumina GAIIX. Expression patterns of 46,404 G. soja genes were profiled in all six samples based on RNA-seq data using Cufflinks software. Then, t12 transcription factors from MYB, WRKY, NAC, bZIP, C2H2, HB, and TIFY families and 12 oxidation reduction related genes were chosen and verified to be induced in response to alkaline stress by using quantitative real-time polymerase chain reaction (qRT-PCR). The GO functional annotation analysis showed that besides "transcriptional regulation" and "oxidation reduction," these genes were involved in a variety of processes, such as "binding" and "response to stress." This is the first comprehensive transcriptome profiling analysis of wild soybean root under alkaline stress by RNA sequencing. Our results highlight changes in the gene expression patterns and identify a set of genes induced by NaHCO3 stress. These findings provide a base for the global analyses of G. soja alkaline stress tolerance mechanisms.

Development of Microsatellite Markers Derived from Expressed Sequence Tags of Polyporales for Genetic Diversity Analysis of Endangered Polyporus umbellatus

A large scale of EST sequences of Polyporales was screened in this investigation in order to identify EST-SSR markers for various applications. The distribution of EST sequences and SSRs in five families of Polyporales was analyzed, respectively. Mononucleotide was the most abundant type, followed by trinucleotide. Among five families, Ganodermataceae occupied the most SSR markers, followed by Coriolaceae. Functional prediction of SSR marker-containing EST sequences in Ganoderma lucidum obtained three main groups, namely, cellular component, biological process, and molecular function. Thirty EST-SSR primers were designed to evaluate the genetic diversity of 13 natural Polyporus umbellatus accessions. Twenty one EST-SSRs were polymorphic with average PIC value of 0.33 and transferability rate of 71%. These 13 P. umbellatus accessions showed relatively high genetic diversity. The expected heterozygosity, Nei's gene diversity, and Shannon information index were 0.41, 0.39, and 0.57, respectively. Both UPGMA dendrogram and principal coordinate analysis (PCA) showed the same cluster result that divided the 13 accessions into three or four groups.

Elucidating steroid alkaloid biosynthesis in Veratrum californicum: production of verazine in Sf9 cells

Steroid alkaloids have been shown to elicit a wide range of pharmacological effects that include anticancer and antifungal activities. Understanding the biosynthesis of these molecules is essential to bioengineering for sustainable production. Herein, we investigate the biosynthetic pathway to cyclopamine, a steroid alkaloid that shows promising antineoplastic activities. Supply of cyclopamine is limited, as the current source is solely derived from wild collection of the plant Veratrum californicum. To elucidate the early stages of the pathway to cyclopamine, we interrogated a V. californicum RNA-seq dataset using the cyclopamine accumulation profile as a predefined model for gene expression with the pattern-matching algorithm Haystack. Refactoring candidate genes in Sf9 insect cells led to discovery of four enzymes that catalyze the first six steps in steroid alkaloid biosynthesis to produce verazine, a predicted precursor to cyclopamine. Three of the enzymes are cytochromes P450 while the fourth is a γ-aminobutyrate transaminase; together they produce verazine from cholesterol.

Optimal fertilizer application for Panax notoginseng and effect of soil water on root rot disease and saponin contents

Background

Blind and excessive application of fertilizers was found during the cultivation of Panax notoginseng in fields, as well as increase in root rot disease incidence.

Methods

Both “3414” application and orthogonal test designs were performed at Shilin county, Yunnan province, China, for NPK (nitrogen, phosphorus, and potassium) and mineral fertilizers, respectively. The data were used to construct the one-, two-, and three-factor quadratic regression models. The effect of fertilizer deficiency on root yield loss was also analyzed to confirm the result predicted by these models. A pot culture experiment was performed to observe the incidence rate of root rot disease and to obtain the best range in which the highest yield of root and saponins could be realized.

Results

The best application strategy for NPK fertilizer was 0 kg/667 m², 17.01 kg/667 m², and 56.87 kg/667 m², respectively, which can produce the highest root yield of 1,861.90 g (dried root of 100 plants). For mineral fertilizers, calcium and magnesium fertilizers had a significant and positive effect on root yield and the content of four active saponins, respectively. The severity of root rot disease increased with the increase in soil moisture. The best range of soil moisture varied from 0.56 FC (field capacity of water) to 0.59 FC, when the highest yield of root and saponins could be realized as well as the lower incidence rate of root disease.

Conclusion

These results indicate that the amount of nitrogen fertilizer used in these fields is excessive and that of potassium fertilizer is deficient. Higher soil moisture is an important factor that increases the severity of the root rot disease.

Ion Torrent sequencing as a tool for mutation discovery in the flax (Linum usitatissimum L.) genome

BACKGROUND

Detection of induced mutations is valuable for inferring gene function and for developing novel germplasm for crop improvement. Many reverse genetics approaches have been developed to identify mutations in genes of interest within a mutagenized population, including some approaches that rely on next-generation sequencing (e.g. exome capture, whole genome resequencing). As an alternative to these genome or exome-scale methods, we sought to develop a scalable and efficient method for detection of induced mutations that could be applied to a small number of target genes, using Ion Torrent technology. We developed this method in flax (Linum usitatissimum), to demonstrate its utility in a crop species.

RESULTS

We used an amplicon-based approach in which DNA samples from an ethyl methanesulfonate (EMS)-mutagenized population were pooled and used as template in PCR reactions to amplify a region of each gene of interest. Barcodes were incorporated during PCR, and the pooled amplicons were sequenced using an Ion Torrent PGM. A pilot experiment with known SNPs showed that they could be detected at a frequency > 0.3% within the pools. We then selected eight genes for which we wanted to discover novel mutations, and applied our approach to screen 768 individuals from the EMS population, using either the Ion 314 or Ion 316 chips. Out of 29 potential mutations identified after processing the NGS reads, 16 mutations were confirmed using Sanger sequencing.

CONCLUSIONS

The methodology presented here demonstrates the utility of Ion Torrent technology in detecting mutation variants in specific genome regions for large populations of a species such as flax. The methodology could be scaled-up to test >100 genes using the higher capacity chips now available from Ion Torrent.

The location and translocation of ndh genes of chloroplast origin in the Orchidaceae family

The NAD(P)H dehydrogenase complex is encoded by 11 ndh genes in plant chloroplast (cp) genomes. However, ndh genes are truncated or deleted in some autotrophic Epidendroideae orchid cp genomes. To determine the evolutionary timing of the gene deletions and the genomic locations of the various ndh genes in orchids, the cp genomes of Vanilla planifolia, Paphiopedilum armeniacum, Paphiopedilum niveum, Cypripedium formosanum, Habenaria longidenticulata, Goodyera fumata and Masdevallia picturata were sequenced; these genomes represent Vanilloideae, Cypripedioideae, Orchidoideae and Epidendroideae subfamilies. Four orchid cp genome sequences were found to contain a complete set of ndh genes. In other genomes, ndh deletions did not correlate to known taxonomic or evolutionary relationships and deletions occurred independently after the orchid family split into different subfamilies. In orchids lacking cp encoded ndh genes, non cp localized ndh sequences were identified. In Erycina pusilla, at least 10 truncated ndh gene fragments were found transferred to the mitochondrial (mt) genome. The phenomenon of orchid ndh transfer to the mt genome existed in ndh-deleted orchids and also in ndh containing species.

Calcium mobilization in salicylic acid-induced Salvia miltiorrhiza cell cultures and its effect on the accumulation of rosmarinic acid

Ca(2+) serves as a second messenger in plant responses to different signals, and salicylic acid (SA) has been recognized as a signal mediating plant responses to many stresses. We recently found that SA treatment led to the cytoplasmic acidification of Salvia miltiorrhiza cells and alkalinization of extracellular medium. Here, we demonstrate that SA can rapidly induce Ca(2+) mobilization in protoplasts, but the induction can be blocked with a channel blocker of either plasma or organellar membranes. Following SA, A 23187, or 10 mmol/L Ca(2+) treatment, rosmarinic acid (RA) accumulation reached the highest level at 16 h, whereas the peak was found at 10 h if plasma membrane channel blockers were used. By contrast, the highest accumulation of RA occurred at 16 h when organellar channels were blocked, exhibiting the same tendency with SA-induced cells. In agreement with these observations, both phenylalanine ammonia-lyase (PAL) activity and its gene expression detected by real-time PCR also showed the same patterns. These results indicate that SA treatment firstly results in calcium release from internal stores, which in turn leads to PAL activity increase, RA accumulation, and a large amount of Ca(2+) influx from apoplast after 10 h of SA induction.

A Plant-Produced Bacteriophage Tailspike Protein for the Control of Salmonella

The receptor binding domain of the tailspike protein Gp9 from the P22 bacteriophage was recently shown to reduce Salmonella colonization in the chicken gut. In this study, we transiently expressed the receptor binding domain of the Gp9 tailspike protein in Nicotiana benthamiana, and targeted it to the endoplasmic reticulum (ER) or to the chloroplasts. Gp9 was also fused to either an elastin-like polypeptide (ELP) or hydrophobin I tag, which were previously described to improve accumulation levels of recombinant proteins. The highest levels of recombinant protein accumulation occurred when unfused Gp9 was targeted to the ER. Lower levels of chloroplast-targeted Gp9 were also detected. ELP-fused Gp9 was purified and demonstrated to bind to Salmonella enterica serovar Typhimurium in vitro. Upon oral administration of lyophilized leaves expressing Gp9-ELP to newly hatched chickens, we found that this tailspike protein has the potential to be used as a therapeutic to control Salmonella contamination in chickens.

Phylotranscriptomic analysis of the origin and early diversification of land plants

Reconstructing the origin and evolution of land plants and their algal relatives is a fundamental problem in plant phylogenetics, and is essential for understanding how critical adaptations arose, including the embryo, vascular tissue, seeds, and flowers. Despite advances in molecular systematics, some hypotheses of relationships remain weakly resolved. Inferring deep phylogenies with bouts of rapid diversification can be problematic; however, genome-scale data should significantly increase the number of informative characters for analyses. Recent phylogenomic reconstructions focused on the major divergences of plants have resulted in promising but inconsistent results. One limitation is sparse taxon sampling, likely resulting from the difficulty and cost of data generation. To address this limitation, transcriptome data for 92 streptophyte taxa were generated and analyzed along with 11 published plant genome sequences. Phylogenetic reconstructions were conducted using up to 852 nuclear genes and 1,701,170 aligned sites. Sixty-nine analyses were performed to test the robustness of phylogenetic inferences to permutations of the data matrix or to phylogenetic method, including supermatrix, supertree, and coalescent-based approaches, maximum-likelihood and Bayesian methods, partitioned and unpartitioned analyses, and amino acid versus DNA alignments. Among other results, we find robust support for a sister-group relationship between land plants and one group of streptophyte green algae, the Zygnematophyceae. Strong and robust support for a clade comprising liverworts and mosses is inconsistent with a widely accepted view of early land plant evolution, and suggests that phylogenetic hypotheses used to understand the evolution of fundamental plant traits should be reevaluated.

Multiple polyploidy events in the early radiation of nodulating and nonnodulating legumes

Unresolved questions about evolution of the large and diverse legume family include the timing of polyploidy (whole-genome duplication; WGDs) relative to the origin of the major lineages within the Fabaceae and to the origin of symbiotic nitrogen fixation. Previous work has established that a WGD affects most lineages in the Papilionoideae and occurred sometime after the divergence of the papilionoid and mimosoid clades, but the exact timing has been unknown. The history of WGD has also not been established for legume lineages outside the Papilionoideae. We investigated the presence and timing of WGDs in the legumes by querying thousands of phylogenetic trees constructed from transcriptome and genome data from 20 diverse legumes and 17 outgroup species. The timing of duplications in the gene trees indicates that the papilionoid WGD occurred in the common ancestor of all papilionoids. The earliest diverging lineages of the Papilionoideae include both nodulating taxa, such as the genistoids (e.g., lupin), dalbergioids (e.g., peanut), phaseoloids (e.g., beans), and galegoids (=Hologalegina, e.g., clovers), and clades with nonnodulating taxa including Xanthocercis and Cladrastis (evaluated in this study). We also found evidence for several independent WGDs near the base of other major legume lineages, including the Mimosoideae-Cassiinae-Caesalpinieae (MCC), Detarieae, and Cercideae clades. Nodulation is found in the MCC and papilionoid clades, both of which experienced ancestral WGDs. However, there are numerous nonnodulating lineages in both clades, making it unclear whether the phylogenetic distribution of nodulation is due to independent gains or a single origin followed by multiple losses.

Pectinmethylesterases (PME) and pectinmethylesterase inhibitors (PMEI) enriched during phloem fiber development in flax (Linum usitatissimum)

Flax phloem fibers achieve their length by intrusive-diffusive growth, which requires them to penetrate the extracellular matrix of adjacent cells. Fiber elongation therefore involves extensive remodelling of cell walls and middle lamellae, including modifying the degree and pattern of methylesterification of galacturonic acid (GalA) residues of pectin. Pectin methylesterases (PME) are important enzymes for fiber elongation as they mediate the demethylesterification of GalA in muro, in either a block-wise fashion or in a random fashion. Our objective was to identify PMEs and PMEIs that mediate phloem fiber elongation in flax. For this purpose, we measured transcript abundance of candidate genes at nine different stages of stem and fiber development and found sets of genes enriched during fiber elongation and maturation as well as during xylem development. We expressed one of the flax PMEIs in E. coli and demonstrated that it was able to inhibit most of the native PME activity in the upper portion of the flax stem. These results identify key genetic components of the intrusive growth process and define targets for fiber engineering and crop improvement.

Chitinase-like (CTL) and cellulose synthase (CESA) gene expression in gelatinous-type cellulosic walls of flax (Linum usitatissimum L.) bast fibers

Plant chitinases (EC 3.2.1.14) and chitinase-like (CTL) proteins have diverse functions including cell wall biosynthesis and disease resistance. We analyzed the expression of 34 chitinase and chitinase-like genes of flax (collectively referred to as LusCTLs), belonging to glycoside hydrolase family 19 (GH19). Analysis of the transcript expression patterns of LusCTLs in the stem and other tissues identified three transcripts (LusCTL19, LusCTL20, LusCTL21) that were highly enriched in developing bast fibers, which form cellulose-rich gelatinous-type cell walls. The same three genes had low relative expression in tissues with primary cell walls and in xylem, which forms a xylan type of secondary cell wall. Phylogenetic analysis of the LusCTLs identified a flax-specific sub-group that was not represented in any of other genomes queried. To provide further context for the gene expression analysis, we also conducted phylogenetic and expression analysis of the cellulose synthase (CESA) family genes of flax, and found that expression of secondary wall-type LusCESAs (LusCESA4, LusCESA7 and LusCESA8) was correlated with the expression of two LusCTLs (LusCTL1, LusCTL2) that were the most highly enriched in xylem. The expression of LusCTL19, LusCTL20, and LusCTL21 was not correlated with that of any CESA subgroup. These results defined a distinct type of CTLs that may have novel functions specific to the development of the gelatinous (G-type) cellulosic walls.

Horizontal transfer of an adaptive chimeric photoreceptor from bryophytes to ferns

Ferns are well known for their shade-dwelling habits. Their ability to thrive under low-light conditions has been linked to the evolution of a novel chimeric photoreceptor--neochrome--that fuses red-sensing phytochrome and blue-sensing phototropin modules into a single gene, thereby optimizing phototropic responses. Despite being implicated in facilitating the diversification of modern ferns, the origin of neochrome has remained a mystery. We present evidence for neochrome in hornworts (a bryophyte lineage) and demonstrate that ferns acquired neochrome from hornworts via horizontal gene transfer (HGT). Fern neochromes are nested within hornwort neochromes in our large-scale phylogenetic reconstructions of phototropin and phytochrome gene families. Divergence date estimates further support the HGT hypothesis, with fern and hornwort neochromes diverging 179 Mya, long after the split between the two plant lineages (at least 400 Mya). By analyzing the draft genome of the hornwort Anthoceros punctatus, we also discovered a previously unidentified phototropin gene that likely represents the ancestral lineage of the neochrome phototropin module. Thus, a neochrome originating in hornworts was transferred horizontally to ferns, where it may have played a significant role in the diversification of modern ferns.

Identification and functional analysis of mitogen-activated protein kinase kinase kinase (MAPKKK) genes in canola (Brassica napus L.)

Mitogen-activated protein kinase (MAPK) signalling cascades, consisting of three types of reversibly phosphorylated kinases (MAPKKK, MAPKK, and MAPK), are involved in important processes including plant immunity and hormone responses. The MAPKKKs comprise the largest family in the MAPK cascades, yet only a few of these genes have been associated with physiological functions, even in the model plant Arabidopsis thaliana. Canola (Brassica napus L.) is one of the most important oilseed crops in China and worldwide. To explore MAPKKK functions in biotic and abiotic stress responses in canola, 66 MAPKKK genes were identified and 28 of them were cloned. Phylogenetic analysis of these canola MAPKKKs with homologous genes from representative species classified them into three groups (A-C), comprising four MAPKKKs, seven ZIKs, and 17 Raf genes. A further 15 interaction pairs between these MAPKKKs and the downstream BnaMKKs were identified through a yeast two-hybrid assay. The interactions were further validated through bimolecular fluorescence complementation (BiFC) analysis. In addition, by quantitative real-time reverse transcription-PCR, it was further observed that some of these BnaMAPKKK genes were regulated by different hormone stimuli, abiotic stresses, or fungal pathogen treatments. Interestingly, two novel BnaMAPKKK genes, BnaMAPKKK18 and BnaMAPKKK19, which could elicit hypersensitive response (HR)-like cell death when transiently expressed in Nicotiana benthamiana leaves, were successfully identified. Moreover, it was found that BnaMAPKKK19 probably mediated cell death through BnaMKK9. Overall, the present work has laid the foundation for further characterization of this important MAPKKK gene family in canola.

Phylogenetic pinpointing of a paleopolyploidy event within the flax genus (Linum) using transcriptomics

BACKGROUND AND AIMS

Cultivated flax (Linum usitatissimum) is known to have undergone a whole-genome duplication around 5-9 million years ago. The aim of this study was to investigate whether other whole-genome duplication events have occurred in the evolutionary history of cultivated flax. Knowledge of such whole-genome duplications will be important in understanding the biology and genomics of cultivated flax.

METHODS

Transcriptomes of 11 Linum species were sequenced using the Illumina platform. The short reads were assembled de novo and the DupPipe pipeline was used to look for signatures of polyploidy events from the age distribution of paralogues. In addition, phylogenies of all paralogues were assembled within an estimated age window of interest. These phylogenies were assessed for evidence of a paleopolyploidy event within the genus Linum.

KEY RESULTS

A previously unknown paleopolyploidy event that occurred 20-40 million years ago was discovered and shown to be specific to a clade within Linum containing cultivated flax (L. usitatissimum) and other mainly blue-flowered species. The finding was supported by two lines of evidence. First, a significant change of slope (peak) was shown in the age distribution of paralogues that was phylogenetically restricted to, and ubiquitous in, this clade. Second, a large number of paralogue phylogenies were retrieved that are consistent with a polyploidy event occurring within that clade.

CONCLUSIONS

The results show the utility of multi-species transcriptomics for detecting whole-genome duplication events and demonstrate that that multiple rounds of polyploidy have been important in shaping the evolutionary history of flax. Understanding and characterizing these whole-genome duplication events will be important for future Linum research.

Identification, expression and interaction analyses of calcium-dependent protein kinase (CPK) genes in canola (Brassica napus L.)

BACKGROUND

Canola (Brassica napus L.) is one of the most important oil-producing crops in China and worldwide. The yield and quality of canola is frequently threatened by environmental stresses including drought, cold and high salinity. Calcium is a well-known ubiquitous intracellular secondary messenger in plants. Calcium-dependent protein kinases (CPKs) are Ser/Thr protein kinases found only in plants and some protozoans. CPKs are Ca2+ sensors that have both Ca2+ sensing function and kinase activity within a single protein and play crucial roles in plant development and responses to various environmental stresses.

RESULTS

In this study, we mined the available expressed sequence tags (ESTs) of B. napus and identified a total of 25 CPK genes, among which cDNA sequences of 23 genes were successfully cloned from a double haploid cultivar of canola. Phylogenetic analysis demonstrated that they could be clustered into four subgroups. The subcellular localization of five selected BnaCPKs was determined using green fluorescence protein (GFP) as the reporter. Furthermore, the expression levels of 21 BnaCPK genes in response to salt, drought, cold, heat, abscisic acid (ABA), low potassium (LK) and oxidative stress were studied by quantitative RT-PCR and were found to respond to multiple stimuli, suggesting that canola CPKs may be convergence points of different signaling pathways. We also identified and cloned five and eight Clade A basic leucine zipper (bZIP) and protein phosphatase type 2C (PP2C) genes from canola and, using yeast two-hybrid and bimolecular fluorescence complementation (BiFC), determined the interaction between individual BnaCPKs and BnabZIPs or BnaPP2Cs (Clade A). We identified novel, interesting interaction partners for some of the BnaCPK proteins.

CONCLUSION

We present the sequences and characterization of CPK gene family members in canola for the first time. This work provides a foundation for further crop improvement and improved understanding of signal transduction in plants.

Differential expression of miRNAs in Brassica napus root following infection with Plasmodiophora brassicae

Canola (oilseed rape, Brassica napus L.) is susceptible to infection by the biotrophic protist Plasmodiophora brassicae, the causal agent of clubroot. To understand the roles of microRNAs (miRNAs) during the post-transcriptional regulation of disease initiation and progression, we have characterized the changes in miRNA expression profiles in canola roots during clubroot disease development and have compared these to uninfected roots. Two different stages of clubroot development were targeted in this miRNA profiling study: an early time of 10-dpi for disease initiation and a later 20-dpi, by which time the pathogen had colonized the roots (as evident by visible gall formation and histological observations). P. brassicae responsive miRNAs were identified and validated by qRT-PCR of miRNAs and the subsequent validation of the target mRNAs through starBase degradome analysis, and through 5' RLM-RACE. This study identifies putative miRNA-regulated genes with roles during clubroot disease initiation and development. Putative target genes identified in this study included: transcription factors (TFs), hormone-related genes, as well as genes associated with plant stress response regulation such as cytokinin, auxin/ethylene response elements. The results of our study may assist in elucidating the role of miRNAs in post-transcriptional regulation of target genes during disease development and may contribute to the development of strategies to engineer durable resistance to this important phytopathogen.

Identification and characterization of CBL and CIPK gene families in canola (Brassica napus L.)

BACKGROUND

Canola (Brassica napus L.) is one of the most important oil-producing crops in China and worldwide. The yield and quality of canola is frequently threatened by environmental stresses including drought, cold and high salinity. Calcium is a ubiquitous intracellular secondary messenger in plants. Calcineurin B-like proteins (CBLs) are Ca2+ sensors and regulate a group of Ser/Thr protein kinases called CBL-interacting protein kinases (CIPKs). Although the CBL-CIPK network has been demonstrated to play crucial roles in plant development and responses to various environmental stresses in Arabidopsis, little is known about their function in canola.

RESULTS

In the present study, we identified seven CBL and 23 CIPK genes from canola by database mining and cloning of cDNA sequences of six CBLs and 17 CIPKs. Phylogenetic analysis of CBL and CIPK gene families across a variety of species suggested genome duplication and diversification. The subcellular localization of three BnaCBLs and two BnaCIPKs were determined using green fluorescence protein (GFP) as the reporter. We also demonstrated interactions between six BnaCBLs and 17 BnaCIPKs using yeast two-hybrid assay, and a subset of interactions were further confirmed by bimolecular fluorescence complementation (BiFC). Furthermore, the expression levels of six selected BnaCBL and 12 BnaCIPK genes in response to salt, drought, cold, heat, ABA, methyl viologen (MV) and low potassium were examined by quantitative RT-PCR and these CBL or CIPK genes were found to respond to multiple stimuli, suggesting that the canola CBL-CIPK network may be a point of convergence for several different signaling pathways. We also performed a comparison of interaction patterns and expression profiles of CBL and CIPK in Arabidospsis, canola and rice, to examine the differences between orthologs, highlighting the importance of studying CBL-CIPK in canola as a prerequisite for improvement of this crop.

CONCLUSIONS

Our findings indicate that CBL and CIPK family members may form a dynamic complex to respond to different abiotic or hormone signaling. Our comparative analyses of the CBL-CIPK network between canola, Arabidopsis and rice highlight functional differences and the necessity to study CBL-CIPK gene functions in canola. Our data constitute a valuable resource for CBL and CPK genomics.

Data access for the 1,000 Plants (1KP) project

The 1,000 plants (1KP) project is an international multi-disciplinary consortium that has generated transcriptome data from over 1,000 plant species, with exemplars for all of the major lineages across the Viridiplantae (green plants) clade. Here, we describe how to access the data used in a phylogenomics analysis of the first 85 species, and how to visualize our gene and species trees. Users can develop computational pipelines to analyse these data, in conjunction with data of their own that they can upload. Computationally estimated protein-protein interactions and biochemical pathways can be visualized at another site. Finally, we comment on our future plans and how they fit within this scalable system for the dissemination, visualization, and analysis of large multi-species data sets.

Characterization and transcript profiling of the pectin methylesterase (PME) and pectin methylesterase inhibitor (PMEI) gene families in flax (Linum usitatissimum)

Background

Pectin methylesterases (PMEs) catalyze the demethylesterification of homogalacturonans in the cell wall; their activity is regulated in part by pectin methylesterase inhibitors (PMEIs). PME activity may result in either rigidification or loosening of the cell wall, depending on the mode of demethylesterification. The activity of PMEs in the middle lamella is expected to affect intrusive elongation of phloem fibers, and their adhesion to adjacent cells. Length and extractability of phloem fibers are qualities important for their industrial uses in textiles and composites. As only three flax PMEs had been previously described, we were motivated to characterize the PME and PMEI gene families of flax.

Results

We identified 105 putative flax PMEs (LuPMEs) and 95 putative PMEIs (LuPMEIs) within the whole-genome assembly. We found experimental evidence for the transcription of 77/105 LuPMEs and 83/95 LuPMEIs, and surveyed the transcript abundance of these in 12 different tissues and stages of development. Six major monophyletic groups of LuPMEs could be defined based on the inferred relationships of flax genes and their presumed orthologs from other species. We searched the LuPMEs and LuPMEIs for conserved residues previously reported to be important for their tertiary structure and function. In the LuPMEs, the most highly conserved residues were catalytic residues while in the LuPMEIs, cysteines forming disulfude bridges between helices α2 and α3 were most highly conserved. In general, the conservation of critical residues was higher in the genes with evidence of transcript expression than in those for which no expression was detected.

Conclusions

The LuPMEs and LuPMEIs comprise large families with complex patterns of transcript expression and a wide range of physical characteristics. We observed that multiple PMEs and PMEIs are expressed in partially overlapping domains, indicative of several genes acting redundantly during most processes. The potential for functional redundancy was highlighted also by the phylogenetic analyses. We were able to identify a subset of PME and PMEIs that appeared particularly relevant to fiber development, which may provide a basis for the improvement of key traits in industrial feedstocks and a better understanding of the physiological roles of PMEs and PMEIs in general.

Ecological implications of single and mixed nitrogen nutrition in Arabidopsis thaliana

BACKGROUND

Ecologists recognize that plants capture nitrogen in many chemical forms that include amino acids. Access to multiple nitrogen types in plant communities has been argued to enhance plant performance, access to nitrogen and alter ecological interactions in ways that may promote species coexistence. However, data supporting these arguments have been limited. While it is known that plants uptake amino acids from soil, long term studies that link amino acid uptake to measures of plant performance and potential reproductive effort are not typically performed. Here, a series of experiments that link uptake of nitrate, glutamine or asparagine with lifetime reproductive effort in Arabidopsis thaliana are reported. Nitrogen was offered either singly or in mixture and at a variety of combinations. Traits related to reproductive output were measured, as was the preference for each type of nitrogen.

RESULTS

When plants were supplied with a single nitrogen type at concentrations from 0.1-0.9 mM, the ranking of nitrogen types was nitrate > glutamine > asparagine in terms of the relative performance of plants. When plants were supplied with two types of nitrogen in mixture at ratios between 0.1:0.9-0.9:0.1 mM, again plants performed best when nitrate was present, and poorly when amino acids were mixed. Additionally, stable isotopes revealed that plants preferentially captured nitrogen types matching the hierarchy of nitrate > glutamine > asparagine. Comparing between the two experiments revealed that mixed nitrogen nutrition was a net cost to the plants.

CONCLUSIONS

Plant performance on mixed nitrogen was less than half the performance on equal amounts of any single nitrogen type. We asked: why did A. thaliana capture amino acids when doing so resulted in a net cost? We argue that available data cannot yet answer this question, but hypothesize that access to lower quality forms of nitrogen may become important when plants compete.