Original Research ARTICLE
De novo Assembly of the Pokeweed Genome Provides Insight Into Pokeweed Antiviral Protein (PAP) Gene Expression
- 1Department of Biology, York University, Toronto, ON, Canada
- 2Department of Biology, Center for Genomics and Systems Biology, New York University, New York, NY, United States
Ribosome-inactivating proteins (RIPs) are RNA glycosidases thought to function in defense against pathogens. These enzymes remove purine bases from RNAs, including rRNA; the latter activity decreases protein synthesis in vitro, which is hypothesized to limit pathogen proliferation by causing host cell death. Pokeweed antiviral protein (PAP) is a RIP synthesized by the American pokeweed plant (Phytolacca americana). PAP inhibits virus infection when expressed in crop plants, yet little is known about the function of PAP in pokeweed due to a lack of genomic tools for this non-model species. In this work, we de novo assembled the pokeweed genome and annotated protein-coding genes. Sequencing comprised paired-end reads from a short-insert library of 83X coverage, and our draft assembly (N50 = 42.5 Kb) accounted for 74% of the measured pokeweed genome size of 1.3 Gb. We obtained 29,773 genes, 73% of which contained known protein domains, and identified several PAP isoforms. Within the gene models of each PAP isoform, a long 5′ UTR intron was discovered, which was validated by RT-PCR and sequencing. Presence of the intron stimulated reporter gene expression in tobacco. To gain further understanding of PAP regulation, we complemented this genomic resource with expression profiles of pokeweed plants subjected to stress treatments [jasmonic acid (JA), salicylic acid, polyethylene glycol, and wounding]. Cluster analysis of the top differentially expressed genes indicated that some PAP isoforms shared expression patterns with genes involved in terpenoid biosynthesis, JA-mediated signaling, and metabolism of amino acids and carbohydrates. The newly sequenced promoters of all PAP isoforms contained cis-regulatory elements associated with diverse biotic and abiotic stresses. These elements mediated response to JA in tobacco, based on reporter constructs containing promoter truncations of PAP-I, the most abundant isoform. Taken together, this first genomic resource for the Phytolaccaceae plant family provides new insight into the regulation and function of PAP in pokeweed.
American pokeweed, Phytolacca americana, belongs to the Phytolaccaceae family of flowering plants, which comprises 65 species of herbs, shrubs, and trees. P. americana (pokeweed) is the most well-studied Phytolacca species due to its broad agricultural and medical applications. Pokeweed synthesizes PAP, an N-glycosidase and RIP that depurinates the conserved α-sarcin loop of large rRNAs (Endo et al., 1988). PAP exhibits antiviral activity against diverse plant and animal viruses. Specifically, depurination of rRNA inactivates ribosomes in infected cells, thereby inhibiting host and viral protein synthesis (Lodge et al., 1993; Bonness et al., 1994). PAP also depurinates the genomes of some RNA viruses, interfering with multiple stages of the viral life cycle (He et al., 2008; Karran and Hudak, 2008; Mansouri et al., 2009). Transgenic plants expressing PAP acquire novel antiviral and antifungal activities, making the gene an attractive candidate for use in agricultural engineering (Zoubenko et al., 1997, 2000; Wang et al., 1998; Dai et al., 2003). Pokeweed also shows potential in phytoremediation as a heavy metal hyperaccumulator, thriving in contaminated soil that is otherwise toxic to most plants (Peng et al., 2008; Liu et al., 2010; Zhao et al., 2011). Although this non-model plant displays resistance to diverse biotic and abiotic stresses, a wealth of information remains unknown since its genome has not been sequenced.
Several isoforms of PAP have been reported, exhibiting different temporal (PAP-I, PAP-II, PAP-III) or tissue-specific (PAP-R, PAP-S, S1, S2, PAP-α) expression patterns, or identified during cell culture (PAP-H, PAP-C) (Irvin, 1975; Irvin et al., 1980; Barbieri et al., 1982, 1989; Bolognesi et al., 1990; Kataoka et al., 1992; Rajamohan et al., 1999; Park et al., 2002). Consistent with a hypothesized role in pathogen defense, we showed previously through transcriptomic analysis that expression of several PAP isoforms is up-regulated by JA (Neller et al., 2016). JA is a plant hormone that mediates resistance to insect herbivores, which are viral vectors, and necrotrophic pathogens. Others have reported an induction of RIP expression in various plants upon treatment with phytohormones [JA, salicylic acid (SA), abscisic acid (ABA)] or associated stresses, including insect feeding, pathogen infection, cold, heat, drought, salinity, and mechanical wounding (Reinbothe et al., 1994; Song et al., 2000; Iglesias et al., 2005; Jiang et al., 2008; Qin et al., 2009; Tartarini et al., 2010). Therefore, it is well-established that RIPs are induced by various stresses; however, it is not clear how RIP isoform expression is controlled, or how RIPs are integrated within stress-response pathways.
Pokeweed is tetraploid with a chromosome count of 2n = 36 and 1C-value of 1.48 pg (∼1.5 Gb) (Bennett, 2000; Rice et al., 2015). Although there are no reference genomes available for the Phytolaccaceae family, genomes of some members from the order Caryophylalles (which includes pokeweed) have been sequenced. Genome assembly and annotation has been performed for Beta vulgaris (sugar beet), Spinacia oleracea (spinach), Chenopodium quinoa (quinoa), and Amaranthus hypochondriacus (amaranth) (Dohm et al., 2014; Clouse et al., 2016; Jarvis et al., 2017; Xu et al., 2017). This genomic information is a useful reference for pokeweed, especially since RIPs are prevalent among the Caryophyllales (Di Maro et al., 2014).
Here, we present the first de novo draft genome assembly of pokeweed and an annotation of its protein-coding genes. Using this resource, we investigated the presence and gene organization of PAP isoforms. A novel feature was discovered in the 5′ UTR: a long intron that affects gene expression. Integration of RNA-Seq data from pokeweed stress treatments enabled the identification of co-expressed genes and provided insight into defense responses in the plant. Finally, differences in CREs in the promoters of PAP isoforms, combined with their unique expression profiles, suggest that the isoforms have distinct roles in pokeweed. Our study provides a workflow that integrates genome assembly, annotation, and differential expression analysis for a non-model plant. These resources will facilitate the study of pokeweed to understand its ability to survive environmental stress.
Materials and Methods
Genomic DNA Sequencing and Genome Assembly
A summary of our study is provided in Figure 1. High-quality genomic DNA was isolated from a single pokeweed plant using a CTAB-based extraction method (Healey et al., 2014) and sent to McGill University and Genome Quebec Innovation Centre (Montreal, QB, Canada) for sequencing. Shotgun sequencing was performed on one lane of an Illumina HiSeq 2500 instrument in Rapid Run mode. Paired-end reads of 250 bp were obtained from genomic DNA fragments having an average size of 400 bp; this strategy was chosen to conform with input recommendations of the downstream assembly software. Following sequencing, adapters and low-quality bases (Q < 30, averaged over four bases) were removed with Trimmomatic (v. 0.36; Bolger et al., 2014). The pokeweed genome was assembled with Discovar De Novo using default parameters (v. 52488; Love et al., 2016) on the Sharcnet high-performance computing cluster1. Raw genomic DNA sequencing reads are available at the SRA under project # PRJNA544344. The genome completion score was measured with BUSCO (v. 3; Waterhouse et al., 2018).
Figure 1. Summary of genome assembly, annotation, and differential gene expression analysis. The assembled pokeweed genome was annotated using three sources of information provided to MAKER, indicated with dashed arrows: (i) published plant protein databases, (ii) the pokeweed de novo transcriptome assembly, (iii) the pokeweed genome-guided transcriptome assembly, derived from aligning unassembled (‘cleaned’) reads to the genome (①) and joining exonic regions into transcripts (②). Reads from (①) that aligned to exonic regions of a gene (defined by the annotated gene models) were summed. Counts were used as input for differential gene expression analysis. The top differentially expressed genes were clustered, and functional enrichment testing was performed on each cluster.
Determination of Genome Size
DNA contents of somatic nuclei (pg/2C) were measured using flow cytometry. For each sample, approximately 1 cm2 of fresh pokeweed leaf and 2 cm2 of fresh Sorghum bicolor Pioneer 8695 (1.74 pg/2C; Johnston et al., 1999) were chopped with a razor blade in 0.7 mL of ice-cold LB01 buffer (Dolezel et al., 1989) containing 100 μg/mL propidium iodide and 50 μg/mL RNAse A. Samples were stained at room temperature for 20–25 min, then tested at low speed using a BD FACSCalibur flow cytometer (BD Biosciences, San José, CA, United States). The FL2 detector (585/42 nm) was used to measure fluorescence, with integrated fluorescence as the parameter of interest. 2C DNA contents were calculated based on the relative fluorescence of the pokeweed and Sorghum G0/G1 nuclei and the known DNA content of Sorghum. Means, coefficients of variation, and nuclei numbers for the nuclei fluorescence peaks were measured using flowPloidy (v. 1.7.0; Smith et al., 2018). Nine pokeweed plants were tested in total: four as individuals and five as bulk samples of two or three plants.
Stress Treatments, mRNA Sequencing, and Transcriptome Assembly
Total RNA was extracted from four-leaf pokeweed plants subjected to the following treatments: sprayed with 5 mM JA or SA (solubilized in 0.5% ET), watered every 3 days for a 7-day period with 10% PEG, or wounded with forceps (WND). Plants sprayed with 0.5% ET or WT served as controls for JA/SA and PEG/WND, respectively. Leaf tissue was flash-frozen in liquid nitrogen 24 h following treatment for JA, SA, ET, and WND samples, and 3 days after the final treatment for PEG and WT samples. RNA-Seq libraries were constructed with the TruSeq Stranded mRNA Library Preparation Kit (RS-122-2101, Illumina). For each condition, four biological replicate libraries were prepared from equal amounts of RNA pooled from three independent plants (i.e., 24 libraries derived from 72 total plants). Strand-specific, paired-end reads of 125 bp were sequenced on two lanes of an Illumina HiSeq 2500 instrument by The Centre for Applied Genomics (The Hospital for Sick Children, Toronto, ON, Canada). Raw mRNA sequencing reads are available at the SRA under project # PRJNA309999.
In addition to the data generated in the present study, mRNA transcriptome assembly incorporated two publicly available P. americana RNA-Seq datasets from the SRA: Accessions SRX2774676 and ERX2099309, which were derived from whole plants (root, stem, leaf, and flower tissue), as well as high-coverage RNA-Seq datasets (n = 3 biological replicates) from our previous study of the pokeweed leaf mRNA transcriptome (Neller et al., 2016). All datasets (32 in total) were processed with Trimmomatic as described above. Two transcriptome assembly strategies were employed. For de novo assembly, reads were combined into a single reference transcriptome using Trinity (v. 2.5.1; Grabherr et al., 2011; Haas et al., 2013). For genome-guided assembly, RNA-Seq libraries were independently aligned to the assembled genome with HISAT (v. 2.1.0; Kim et al., 2015) and transcript reconstruction was performed with Stringtie (v. 1.3.4b; Pertea et al., 2015). The independent assemblies were combined into a single, non-redundant reference transcriptome using the Stringtie ‘merge’ function.
Only contigs of minimum length 10 Kb were annotated, as these were most likely to contain full-length protein-coding genes (Campbell et al., 2014a). Prior to annotation, a pokeweed-specific repeat library was prepared with RepeatModeler (v. 1.0.9)2. Repeat-masking and gene prediction were performed using the MAKER pipeline (v. 2.31.8; Campbell et al., 2014a,b). Both pokeweed-specific and simple repeats were masked with RepeatMasker (v. 4.0.7)3. RepeatRunner was used to mask divergent protein-coding portions of retro-elements and retro-viruses not identified by RepeatMasker (Smith et al., 2007).
Gene prediction followed the protocol outlined in Campbell et al. (2014a). The Trinity (de novo) and Stringtie (genome-guided) transcriptomes described above were provided as transcript evidence. The proteomes of published Caryophyllales species (sugar beet, spinach, quinoa, amaranth) were used as homologous protein evidence (Dohm et al., 2014; Clouse et al., 2016; Jarvis et al., 2017; Xu et al., 2017). The quinoa and amaranth proteomes were obtained from Phytozome4 (v. 1.0 and v. 2.1, respectively), spinach from SpinachBase5 (v. 1), and sugar beet from The Beta vulgaris Resource (Refbeet-1.1)6. Protein evidence was further supplemented with the plant subset of the SwissProt database (obtained October 2017; Pundir et al., 2017). Three rounds of iterative gene prediction were performed with MAKER. In Round 1, predictions were inferred directly from transcript and protein evidence (est2genome, protein2genome = 1). Top-scoring gene models from Round 1 (AED ≤ 0.25, amino acids ≥ 50) were used to train the ab initio gene predictors SNAP (Korf, 2004) and AUGUSTUS (v. 3.2.3; Stanke et al., 2004). In Round 2, MAKER was re-run with ab initio gene predictors turned on, and top-scoring models were used to re-train predictors as per above. A final Round 3 was run with the twice-trained gene predictors.
Following gene prediction, an Interproscan (v. 5.28-67.0; Zdobnov and Apweiler, 2001) search was performed to identify protein-coding (Pfam) domains and obtain associated GO terms. GO annotations were supplemented with non-redundant GO terms from the SwissProt database for the best hit as per BLAST-P analysis. Each round of genome annotation was inspected visually with JBrowse (v. 1.12.3 Buels et al., 2016). In Rounds 2 and 3, we observed spurious fusion of protein-coding gene models with putative pseudogenes; this could not be resolved despite several attempts of re-training and post-annotation processing. Therefore, the final gene set was conservatively chosen to comprise only the evidence-based annotations from Round 1. The completion score of the annotated gene set was measured with BUSCO as described above, as well as relative to CoreGFs of the ‘green plants’ subset from PLAZA v. 2.5 (Veeckman et al., 2016).
Validation of PAP Isoform Gene Models
cDNAs were generated by reverse transcribing 0.5 μg of total pokeweed RNA with SuperScript III reverse transcriptase (25 units; Thermo Fisher) and isoform-specific primers. To validate PAP isoform gene models at both the genomic and mRNA levels, PCR was performed with isoform-specific forward and reverse primers using either pokeweed gDNA or cDNA as the starting template. All PCR amplifications were conducted using Q5 High-Fidelity DNA polymerase (1 unit; New England Biolabs) and 0.5 μM of each primer. PCR products were gel-purified and used as templates for a second round of PCR, this time with forward and reverse primers that contained additional sequences for cloning. Amplicons from the second round of PCR were gel-purified and cloned into the multiple cloning site of pHSG299 vector using the one-step SLIC (sequence- and ligation-independent cloning) method (Jeong et al., 2012). All constructs were sequenced and compared with computationally derived gene and mRNA models. Primers used in this study are listed in Supplementary Table 1.
The longest isoform per gene was obtained for each Caryophyllales proteome noted above, and orthogroup assignment was performed with OrthoFinder (v. 2.2.6; Emms and Kelly, 2015). Functional enrichment analysis of genes from pokeweed-specific orthogroups was conducted with the GOseq package (Young et al., 2010) in the statistical program R (R Foundation for Statistical Computing, 2016). Pokeweed-specific genes were tested for enriched GO terms relative to all pokeweed genes.
Differential Gene Expression Analysis
HTSeq (v. 0.8.0; Anders et al., 2015) was used to sum exon-level counts per gene for each independent RNA-Seq library aligned to the pokeweed genome. Gene counts were provided as input for differential expression testing in the Bioconductor package EdgeR (v. 3.7; Robinson et al., 2010). All pairwise tests were performed with four biological replicates per treatment, and genes with FDR < 0.05 in at least one comparison were considered differentially expressed. Treatment-responsive genes were identified based on the relevant tests: ET vs. JA or SA, and WT vs. PEG or WND. Top DEGs were defined as those with FDR < 0.001 and FC > 4 in at least one pairwise comparison. These genes were clustered with DEclust (v. 1.0.1; Aoto et al., 2017) based on a multi-conditional expression profile of EdgeR results and normalized abundance in TPM. GO enrichment analysis was performed on the genes in each cluster as described above.
Quantitative RT-PCR (qRT-PCR) Validation of PAP Isoform mRNA Levels
cDNAs were generated by reverse transcribing 0.5 μg of total pokeweed RNA with SuperScript III reverse transcriptase (25 units; Thermo Fisher) and gene-specific primers (Supplementary Table 1) according to the manufacturer’s instructions. Elongation factor-1-gamma (EF1G) and the cell wall protein BIIDXI served as internal controls as these transcripts were stably expressed under our stress treatments according to our RNA-Seq differential expression analysis. Following cDNA synthesis, 5 μL of the reverse transcription (RT) reaction product was combined with 0.3 μM forward primer, 0.3 μM reverse primer, and 33 μL of 2X SYBR Green qPCR Master Mix (Bimake), to a final volume of 66 μL. Each reaction was split into three technical replicates. qRT-PCRs were conducted in a QIAGEN RotorGene Q thermocycler with the following settings: hold at 50°C for 20 s, initial denaturation and hot-start DNA polymerase activation for 10 min, followed by 40 cycles alternating between denaturation (95°C; 15 s) and combined annealing/extension (68°C, 45 s). mRNA levels were quantified using the ΔΔCt method (Livak and Schmittgen, 2001). A melting curve analysis was performed to confirm the presence of a single PCR product after each reaction. At least three biological replicates per treatment were conducted for each transcript.
Generation of Promoter-GUS Reporter Constructs
The 1262 bp region upstream of the PAP-I TSS was considered the proximal PAP-I promoter. This ∼1.3 Kb promoter, along with the 5′ UTR and 1.6 Kb intron, were PCR-amplified from pokeweed genomic DNA (500 ng) with the primers PAP-I-prom-SLIC-FOR and PAP-I-prom-SLIC-REV (0.5 μM each). Amplicons were gel-purified and cloned by one-step SLIC (Jeong et al., 2012) into the multiple cloning site of pHSG299. All PCR amplifications were conducted using Q5 High-Fidelity DNA polymerase (1 unit; New England Biolabs) using the manufacturer’s instructions. The pHSG299 plasmid containing the PAP-I promoter and intron then served as the PCR template for all downstream PCRs of PAP-I promoter fragments. PAP-I promoter fragments containing the 1.6 Kb intron (1262-int and 102-int) were generated through PCR by pairing the same reverse primer (5-UTR-SLIC-REV) with different forward primers (P1-FOR or P7-FOR). Primer sequences are listed in Supplementary Table 1.
To produce intronless versions of promoter constructs, a reverse primer (5-UTR-no-int-REV) was designed to connect the two portions of the PAP-I 5′ UTR that were originally interrupted by the intron. This reverse primer was paired with different forward primers (P1-FOR to P7-FOR) to produce 5′ promoter truncations (1262, 1124, 711, 584, 432, 296, and 102). PCR products were gel-purified and used as templates for a second round of PCR to attach sequences needed for SLIC cloning.
Amplified fragments were cloned in place of the CaMV 35S promoter in pCambia 0305.2 using one-step SLIC (Jeong et al., 2012) and transformed into Escherichia coli DH5α. For constructs used to investigate the effect of the PAP-I leader intron on gene expression, the castor bean catalase intron in the GUS reporter gene was removed to determine the influence of the PAP-I intron alone. Minus catalase intron constructs were made by excluding the intron through PCR (Primer pairs: pCambia-1-no-cat-FOR and pCambia-1-no-cat-REV; pCambia-2-no-cat-FOR and pCambia-2-no-cat-REV) and reassembling the two fragments through Gibson assembly (Cat# E2611S; New England Biolabs). Positive E. coli transformants were screened through colony PCR, and all reporter gene constructs were confirmed by sequencing.
Constructs were transformed into Rhizobium radiobacter (syn. Agrobacterium) AGL1 strain by electroporation as previously described (Wise et al., 2006). Electrocompetent cells (20 μL) mixed with plasmid DNA (50 ng) were electroporated (2.5 kV, 25 μF capacitance, and 400 Ω resistance) and allowed to recover in 2 mL of non-selective YEP medium for 2 h. After recovery, 100 μL of cells were plated on selective YEP agar (50 μg/mL carbenicillin, 50 μg/mL kanamycin) and incubated at 28°C to allow colonies to form.
Measurement of Reporter Activity (GUS Histochemical and Fluorometric Assay)
Agrobacterium cultures harboring promoter-reporter gene constructs were used to agroinfiltrate leaves of four-leaf stage Nicotiana tabacum plants as previously described (Zhao et al., 2017). Cultures were grown in selective YEP liquid medium (50 μg/mL carbenicillin, 50 μg/mL kanamycin) until late log phase (OD600 = 0.7 – 1.0). Cells were then pelleted by centrifugation, washed in agroinfiltration solution (10 mM MES-KOH, pH 5.6, 10 mM MgCl2, 200 μM acetosyringone), and resuspended in agroinfiltration solution to a final OD600 of 0.5. Agrobacterium cells were injected into the abaxial surface of leaves using a needleless syringe. For the JA experiment, leaves were treated with either 0 mM JA (0.5% ET; mock) or 5 mM JA 24 h after agroinfiltration. Leaf disks from inoculated plants were harvested 72 h post agroinfiltration and either used directly for histochemical assays or stored in liquid nitrogen until processing for the fluorometric assay.
GUS histochemical assays were performed according to Jefferson et al. (1987) with some modifications. Fresh leaf disks (0.5 cm diameter) from inoculated plants (minimum of three plants per construct) were vacuum-infiltrated with 5-bromo-4-chloro-3-indolyl-b-D-glucuronic acid, cyclohexylammonium salt (X-Gluc) solution, and incubated overnight at 37°C. After GUS staining, leaf disks were cleared of chlorophyll by washing in increasing concentrations of ET (70–100%) for 48 h.
GUS fluorometric assays were performed in black, clear-bottom 96-well plates according to Côté and Rutledge (2003), with some modifications. Frozen leaf disks (1 cm diameter; 2 disks per sample) from inoculated plants (minimum of four plants per construct) were combined with 200 mg of glass beads (1.0 mm, BioSpec) and 300 μL of GUS extraction buffer (50 mM phosphate buffer, pH 7.0, 10 mM DTT, 10 mM EDTA, 0.1% SDS, 0.1% Triton X-100, 10 mM β-mercaptoethanol), and homogenized using a 3110BX MiniBeadBeater (BioSpec) for 120 s at speed 48. Tissue debris was removed by centrifugation at 4°C and cleared plant extracts (10 μL) were mixed with 720 μL of 0.1 mM 4-methylumbelliferyl-β-D-glucuronide hydrate (4-MUG) and incubated at 37°C. Beginning at 0 min, a 10-μL aliquot was taken from each reaction every 15 min for a total of 60 min, and pipetted into a well containing 180 μL of stop buffer (0.2 M Na2CO3). Three technical replicates per time point were taken for each sample. Fluorescence values were measured at room temperature using a Synergy H4 Hybrid microplate reader (excitation: 365 nm; emission: 455 nm) and compared to a previously determined 4-methylumbelliferone (4-MU) standard curve. GUS activity was calculated from the linear slope of the fluorescence readings and normalized to the total protein concentration, which was determined using a BCA Reducing Agent Compatible Protein Assay Kit (G-Biosciences). Comparisons between mock-treated and JA-treated samples (p < 0.01) were performed for each promoter construct using two-tailed t-tests.
Identification of Putative PAP Promoter CREs and JA-Responsive Transcription Factors
As with PAP-I, the 1.3 Kb sequence upstream of the TSS was considered the proximal promoter for each PAP isoform. Promoter sequence identity analysis was performed using Clustal Omega7 with default parameters. Putative plant-specific CREs were identified using PLACE8 and PlantPan 2.09 web interfaces (Higo et al., 1999; Chow et al., 2016). To reduce false positives, only sequences with ≥90% identity to the published motifs were included. Putative TFs associated with CREs were identified in the pokeweed genome based on annotated Pfam domains, and their differential expression results were assessed for JA-responsiveness (FDR < 0.05).
Assessment of the Pokeweed Genome Assembly
Statistics of the pokeweed genome assembly are shown in Table 1. Based on 1 Kb+ contigs, the assembly is 0.93 Gb in size and has a read coverage of 83X. The contig and scaffold N50 values are 35.2 and 42.5 Kb, respectively, where a scaffold represents the single highest coverage path through each line of the genome assembly graph. The scaffold assembly was used for all downstream analysis. The assembly had a BUSCO genome completion score of 84.3%; of a possible 1440 plant BUSCOs, 1214 were identified as complete (1142 single-copy and 72 duplicated), 82 were fragmented, and 144 were missing. For an additional metric of assembly completion, we determined the expected size of the pokeweed genome using flow cytometry (Supplementary Figure 1). Based on triplicate samples of four individual plants and two non-replicate bulk samples, the pokeweed plants used in this study had a genomic content (mean ± SEM) of 2.57 pg/2C ± 0.0019. Variability was low, with all 14 samples measuring either 2.57 or 2.58 pg/2C. Pokeweed and Sorghum nuclei fluorescence peaks had coefficients of variation <3.3% and included data from at least 1,000 nuclei. Given the conversion that 1 pg DNA = 980 Mb, the haploid genome size of pokeweed is estimated to be 1.26 Gb. Therefore, the present assembly represents 74% of the expected genome size.
Annotation of the Pokeweed Genome
Scaffold contigs of minimum length 10 Kb were annotated with MAKER, as these were most likely to contain protein-coding genes (Campbell et al., 2014a). Accordingly, the mapping rate of pokeweed RNA-Seq reads aligned to these 10 Kb+ contigs averaged 93% over all samples. The annotation file and sequences of annotated transcripts and proteins are provided in Supplementary Data Sheets 1–3. A summary of the annotation is shown in Table 2. From 22,292 contigs, we obtained 29,773 genes (mean length = 5,072 bp) and 56,538 mRNA transcripts (mean length = 1,720 bp). Importantly, 73% of genes contained a Pfam domain, indicating that the majority are protein-coding, and 99% of genes had an AED score < 0.5 (Supplementary Figure 2), demonstrating excellent correspondence between evidence and gene models. The latter result is consistent with our evidence-only genome annotation; that is, gene models were based solely on the provided transcript and protein information rather than ab initio prediction.
Although we had attempted to include iterative training and ab initio prediction, doing so resulted in fused gene models that could not be validated by RT-PCR and PCR. Specifically, following two iterative training rounds we obtained an increase of only 72 genes (0.2%) but a substantial (30%) increase in mean gene length. Visual inspection of the annotations suggested that protein domain-containing pseudogenes had become fused with bona fide genes. This issue of spurious gene fusions upon ab initio incorporation persisted despite varying the gene-finding parameters, using MAKER post-processing tools, and applying the AUGUSTUS species model from the well-annotated sugar beet genome. Surprisingly, we observed that ab initio incorporation still resulted in high AED scores: after two iterative training rounds, 95% of genes had an AED < 0.5 (Supplementary Figure 2); this reinforces the importance of visibly inspecting gene models and associated evidence. Given that ab initio gene prediction resulted in only a small number of novel genes at the expense of annotation accuracy, we chose to conduct all downstream analyses with our evidence-based gene set.
Evaluation of the Annotated Gene Set
The BUSCO completion score of the annotated pokeweed gene set was 75.8%. Since the genome assembly received a completion score of 84.3% (above), 90% of the BUSCOs identified in the genome were annotated. We also measured completion of the gene set relative to CoreGFs of the ‘green plants’ subset from PLAZA. The CoreGF reference set has a lower threshold of species conservation than BUSCO and is not limited to single-copy genes (Veeckman et al., 2016). Based on CoreGF analysis, the gene set completion score was 97.6%. Therefore, in the case of a non-model species, a low BUSCO score may be more reflective of high evolutionary divergence from reference species than incomplete genome assembly or annotation.
To assess how the pokeweed gene set compared to that of previously annotated Caryophyllales species, we used OrthoFinder to identify orthologous gene families (orthogroups) among pokeweed, sugar beet, quinoa, spinach, and amaranth. The species-specific distribution of orthogroups is shown in Figure 2 and Supplementary Table 2. Pokeweed genes were distributed into 14,785 orthogroups comprised of 22,901 genes. Eighteen orthogroups (116 genes) were pokeweed-specific, and these genes were enriched in the GO terms ‘far-red light signaling pathway’ (FDR = 0.0041) and ‘negative regulation of defense response’ (FDR = 0.058). The respective genes were annotated as isoforms of the TF FAR1-Related Sequence (FRS) (PHYAM_025199, PHYAM_007976, PHYAM_022646, PHYAM_016237, PHYAM_019944, PHYAM_019825) and the F-box protein Constitutive Expresser of PR genes 30 (CPR30) (PHYAM_017382, PHYAM_000825, PHYAM_000824, PHYAM_016094). For pokeweed, 77% of genes were assigned to orthogroups and 0.4% of genes were species-specific. This is in line with other Caryophyllales species, for which orthogroup-assigned genes ranged from 72% (sugar beet) to 82% (amaranth) and species-specific genes ranged from 0.1% (spinach) to 0.4% (quinoa). Taken together, these results demonstrate that the pokeweed gene set is consistent with that of well-annotated Caryophyllales species. Additionally, we have identified a subset of genes that may contribute to distinct biological relevance in pokeweed.
Figure 2. Identification of pokeweed-specific genes by orthogroup analysis of Caryophyllales species. Orthogroup assignment was performed with OrthoFinder, using the longest representative protein per gene for each species. The species distribution of orthogroups is shown. Enriched GO terms from pokeweed-specific genes (n = 116) are indicated in the box.
Identification of RIP Genes in Pokeweed
Following genome annotation, we identified genes containing a RIP domain by performing an Interproscan search of pokeweed protein sequences against the Pfam database. This analysis revealed 10 RIP domain-containing genes, summarized in Table 3. PAP isoform annotation was made based on a BLAST-P search against the SwissProt database, retaining the best-scoring hit per protein. Three RIP domain-containing genes had 100% identity and coverage with SwissProt sequences of PAP-I, PAP-α, and PAP-S, respectively. PAP-II was also present with 99% identity and 100% coverage. One gene had 77% identity and 100% coverage with PAP-I from SwissProt, but BLAST-N against pokeweed nucleotide sequences from GenBank showed 99% identity and coverage with the partial genomic clone of PAP-S2 (Accession # AB071855.1). Furthermore, a gene annotated as PAP-S above, had 99% identity and 100% coverage with the partial genomic clone of PAP-S1 (Accession # AB071854.1). Therefore, two PAP-S isoforms exist, namely PAP-S1 and PAP-S2, which agrees with a previous finding (Honjo et al., 2002) and clarifies the single PAP-S notation in SwissProt. We also identified a gene encoding a transcript that we previously reported (c18776_g1_i1) as a potential new PAP isoform (Neller et al., 2016). The gene contains a RIP domain but only has 38% sequence identity and 96% coverage with its best hit, PAP-α, which supports this gene as a novel PAP isoform.
Our analysis also led to the identification of four RIP domain-containing genes that are likely PAP pseudogenes (Table 3). Two transcribed pseudogenes were located on the same contig as PAP-S2. Their gene models were truncated relative to the associated transcript evidence, resulting in 46% and 32% coverage with the respective hits by BLAST-P. Upon closer inspection, in-frame stop codons were observed in all reading frames for both genes; this explained the shortened models, which arose from annotating the longest open reading frame per gene. Two other probable pseudogenes were identified, but their transcripts were absent in either one or both transcriptome assemblies, suggesting that the genes are transcriptionally silent. To the best of our knowledge, this is the first report of potential pseudogenes of PAP.
Presence of a Long Leader Intron in Gene Models of PAP Isoforms
Gene models of protein-coding PAP isoforms are shown in Figure 3A. With exception of PAP-II, the coding sequence of all isoforms was annotated as a single exon of ∼900 bp. In contrast, the coding sequence of PAP-II comprised two exons, separated by an intron of 736 bp, which agrees with a previous report (Poyet, 1997). Interestingly, the gene models of all isoforms revealed a long intron within the 5′ UTR, ranging from 1.5 Kb (PAP-α) to 5.7 Kb (PAP-II). Based on the distribution of intron lengths in pokeweed, the PAP leader introns were longer than 83% (PAP-α) to 94% (PAP-II) of all introns. Several isoforms (PAP-II, PAP-S1, PAP-S2) were predicted to have multiple gene models that differed only in their 5′ UTRs, which suggested the use of alternative promoters in some cases (PAP-IIA/PAP-IIB; PAP-S2A/PAP-S2B). In addition to long leader introns, all PAP transcripts had potential upstream open reading frames in their 5′ UTRs. The majority of transcripts contained a single upstream open reading frame of four codons in length (PAP-I, PAP-α, PAP-S2A, PAP-S2B), while other transcripts contained longer ones (PAP-IIA: 21 codons, 17 codons, 11 codons; PAP-IIB: 8 codons; PAP-S1: 25 codons).
Figure 3. Identification of a novel intron in the 5′ UTR of PAP genes. (A) Gene models of PAP isoforms obtained from annotation with the MAKER pipeline. Arrows indicate primer binding sites used for gene model validation. (B) Validation of PAP gene models through PCR and RT-PCR from pokeweed genomic DNA (gDNA) and cDNA, respectively. All products were validated by sequencing. (C) Effect of the 5′ UTR intron on PAP-I gene expression. GUS reporter constructs were created with the 1262 or 102 bp PAP-I promoter, with or without the 5′ UTR intron. Constructs were agroinfiltrated into tobacco and stained for GUS. CaMV 35S = positive control; untransformed (UT) Agrobacterium = negative control. Four independent plants per construct were tested.
Where possible, gene models were validated by sequenced RT-PCR and PCR products from pokeweed total RNA and genomic DNA, respectively. As shown in Figure 3B, for each gene-specific primer pair, the PCR product size was consistent with that expected from the gene model for both cDNA and genomic DNA. The gene model of the putative novel isoform (PHYAM_012451) is not indicated in Figure 3 because its 5′ and 3′ UTRs could not be annotated; this discrepancy may be solved in the future by scaffolding with longer sequence reads. Nonetheless, we report here the validated gene models of all published PAP isoforms and the finding of a novel, conserved feature: a long intron within the 5′ UTR.
To investigate if the 5′ UTR intron affected PAP gene expression, we created GUS reporter constructs containing the PAP-I proximal promoter (1262 bp) or minimal promoter (102 bp, putative CAAT and TATA boxes only), either with or without the 5′ UTR intron. A reporter construct with the 35S CaMV promoter served as the positive control, and untransformed Agrobacterium was the negative control. Agroinfiltration of the constructs into tobacco leaves and subsequent GUS staining revealed that PAP-I promoter constructs with the intron had higher expression than those without (Figure 3C). The effect of the intron was most evident for the 102 bp promoter, where presence of the intron increased the level of GUS substantially, relative to the nearly undetectable level of staining without the intron. Based on this preliminary characterization, we hypothesize that the 5′ UTR intron, which is a conserved feature of PAP gene models, enhances PAP gene expression.
Identification of Stress-Responsive Genes in Pokeweed
To gain insight into the response of pokeweed to biotic and abiotic stresses, we identified DEGs from genome-aligned RNA-Seq reads derived from several conditions. Pokeweed plants were treated with JA, SA, PEG, or WND, with ET or WT plants serving as controls for JA/SA and PEG/WND, respectively. EdgeR differential expression results and TPM-normalized expression values of all genes are provided in Supplementary Data Sheets 4, 5. Figure 4A shows a heat map of normalized expression of the top DEGs (FDR < 0.001, FC > 4 in at least one pairwise comparison). It total, 3,548 genes were differentially expressed at this threshold. The treatment-specific distribution of DEGs (FDR < 0.05) is shown in Figure 4B. The number of DEGs per treatment was as follows: SA (10,310), JA (9,088), PEG (1,549), WND (568); therefore, pokeweed was most responsive to SA and JA, which both mediate pathogen defense, and much less responsive to the abiotic stresses. Interestingly, 58 DEGs (1.6%) were common to all four treatments, including two PAP isoforms (PAP-II and PAP-α). These common DEGs were significantly enriched in the following GO terms (FDR < 0.05): ‘proline metabolic process,’ ‘ornithine metabolic process,’ and ‘putrescine biosynthetic process from arginine, using agmatinase.’ The term ‘rRNA N-glycosylase activity’ was also highly enriched (FDR = 0.28), suggesting that PAP contributes to more widespread stress responses in pokeweed than previously known. Furthermore, the finding that not all PAP isoforms were responsive to all treatments suggests that isoform expression is differentially regulated.
Figure 4. Identification of stress-responsive genes in pokeweed. (A) Heat map of normalized expression values [log2(TPM + 1), median-centered] of the top differentially expressed genes (DEGs; FDR < 0.001, FC > 4 in at least one pairwise comparison). SA, JA, ET, WND, WT, and PEG denote salicylic acid, jasmonic acid, ethanol, wounding, water, and polyethylene glycol treatments, respectively. For each condition, four RNA-Seq libraries were prepared from three independent pokeweed plants (i.e., n = 4 pooled biological replicates). (B) Treatment-specific distribution of DEGs (FDR < 0.05). (C) Gene clusters having significant functional enrichment. Top DEGs (from A) were clustered with DEclust. For each cluster, the mean expression profile was plotted as a blue line and the biological relevance of enriched GO terms (FDR < 0.05) was summarized in red font.
To further investigate the unique expression profiles of PAP isoforms and their potential roles in pokeweed, we identified gene clusters from the top DEGs described above. Using DEclust, which extracts statistically significant gene clusters from multi-conditional transcriptome data, 36 clusters were defined (Supplementary Figure 3 and Supplementary Data Sheet 6). Thirteen clusters revealed significant GO term enrichment (FDR < 0.05). Figure 4C shows the expression profiles of these functionally enriched clusters, and Table 4 provides the associated GO terms. Overall, each cluster had a discrete and unified biological theme, indicating successful resolution of co-expressed genes. We summarized the enriched GO terms per cluster into the following biological themes: oxalate synthesis, lipid transport, photosynthesis, DNA replication, SA-mediated signaling/defense, cell wall integrity, detoxification, JA-mediated signaling/defense, amino acid/carbohydrate metabolism, cell cycle regulation, and oxidative stress response. Two PAP isoforms (PAP-S1 and PAP-α) were present in Cluster 16, which was enriched in GO terms comprising JA-associated responses including ‘regulation of JA mediated signaling pathway,’ ‘response to wounding,’ and ‘terpenoid biosynthetic process.’ While Cluster 16 contained JA-upregulated genes, Cluster 20, which included the PAP-S2 isoform, consisted of JA down-regulated genes enriched in the terms ‘fructose-bisphosphate aldolase activity’ and ‘amino acid export.’ Overall, these results provide an indication of how pokeweed responds to diverse stresses and situate PAP within key defense pathways.
Differential Regulation of PAP Isoform Expression
The individual expression profiles of PAP genes, including pseudogenes, are provided in Figure 5A. PAP gene expression changes were also validated through qRT-PCR for all four stress treatments (R2 = 0.8807; Supplementary Figure 4). In addition to showing differences in stress-induced expression change, the isoforms varied greatly in abundance. Among protein-coding isoforms, average abundance in TPM across all samples ranged from 146 (PAP-α) to 9911 (PAP-I). PAP-I was the 24th most expressed gene in the plant under normal conditions (WT) and third most expressed upon JA treatment. As expected, only two of the four PAP pseudogenes showed quantitative evidence of transcriptional expression, and their abundances (TPM = 57 and 8) were much less than those of protein-coding isoforms (average TPM = 2,394). Both transcribed pseudogenes were JA-responsive, one of which (PHYAM_010465) was assigned to Cluster 17 with PAP-II (Figure 5B). Given that pseudogenes can regulate post-transcriptional expression of functional parental genes, our finding that a PAP pseudogene is co-expressed with a PAP isoform may indicate a novel mechanism by which PAP gene expression is controlled.
Figure 5. PAP gene expression profiles. (A) Expression changes of all PAP genes, including pseudogenes, in response to pokeweed stress treatments. Green and red indicate significant up- or down-regulation, respectively. The assigned cluster and average abundance of each isoform are also shown. (B) Expression profiles of PAP-containing clusters. The profile of the relevant PAP isoform (red or black line) is indicated, as well as the mean expression profile of all genes in the cluster (blue line).
PAP-I and PAP-II, the two most abundant isoforms, were assigned to clusters that lacked significant functional enrichment (Figure 5B). PAP-I was responsive to both JA and SA. It was assigned to Cluster 6, which included only eight other genes, six of which had annotated homologs: two TFs from the HD-ZIP homeobox family (HAT5 and ATHB-7), pathogenesis-related protein STH-21, and the enzymes galactinol synthase 2, strigolactone esterase DAD2, and 1,4-dihydroxy-2-naphthoyl-CoA thioesterase. PAP-II, like PAP-α, was responsive to all stresses; however, the two isoforms were assigned to different clusters (17 and 16, respectively). The mean expression profiles of both clusters were similar, exhibiting a prominent peak with JA treatment. Accordingly, Cluster 17 was enriched in the terms ‘response to wounding’ (FDR = 0.10) and ‘jasmonic acid biosynthetic process’ (FDR = 0.23). Assignment of PAP-II and PAP-α to different clusters likely reflects differences in the intensity of their responses to JA (log2FC = 3.24 and 4.87, respectively) and SA (log2FC = 0.75 and 2.53, respectively). The putative novel isoform (PHYAM_012451) was most unlike the others. It responded only to PEG, showing a small but significant decrease in expression (log2FC = −0.72). Because this isoform was not among the top DEGs, it was not included in the clustering analysis; future work will investigate treatments that induce its expression. Taken together, these results demonstrate that PAP isoforms have distinct expression profiles, suggesting that they contribute to unique functions in pokeweed.
JA-Responsiveness of the PAP-I Promoter
Since PAP-I was highly up-regulated with JA (log2FC = 4.5), we hypothesized that the PAP-I promoter contained CREs to mediate this response. PAP-I promoter fragments (ranging from 1262 to 102 bp) were placed upstream of the GUS reporter gene and transiently expressed in tobacco leaves through agroinfiltration (Figure 6A). Apart from those expressing 102:GUS, all plants bearing the PAP-I promoter:GUS constructs showed higher GUS activity following JA treatment, relative to controls without JA, as determined through GUS histochemical and fluorometric assays (Figures 6B,C). Therefore, our results suggest that the region upstream of the TSS (−296 to −103) is sufficient for the response of PAP-I to JA. As shown in Figure 6D, bioinformatic annotation of CREs in this region revealed an element (T/GBOXPINAT2) that binds the master JA signaling regulator MYC (Boter et al., 2004). Additionally, binding sites for TFs of the bHLH, bZIP, and MYB families were present. Specifically, MYB family members bind MYB motifs, while members of bHLH and bZIP families bind the T/GBOXATPIN2 element. Given that pokeweed homologs of several of these TFs were present in JA-associated Cluster 16 (bZIP2, bZIP11, bHLH14, bHLH25, bHLH35, MYB15, MYB44, MYB62, MYB305), the annotated CREs in this region likely contribute to JA-responsiveness of the PAP-I promoter.
Figure 6. Response of the PAP-I promoter to JA. (A) Schematic illustration of PAP-I promoter:GUS constructs. The PAP-I promoter was serially truncated from the 5′ end. +1 denotes the transcription start site (TSS). (B) GUS histochemical assay in tobacco leaves agroinfiltrated with PAP-I promoter:GUS constructs and treated with either 0 mM (mock) or 5 mM JA. Untransformed (UT) Agrobacterium = negative control. Three independent plants per construct were tested. (C) GUS fluorometric assay in tobacco leaves agroinfiltrated with PAP-I promoter:GUS constructs and treated with either 0 mM (mock) or 5 mM JA. At least four independent plants per constructs were tested. Error bars represent the standard error of the mean (SEM). Comparisons between mock-treated and JA-treated samples were conducted for each promoter construct using two-tailed t-tests, p < 0.01. ∗∗p < 0.01; ∗∗∗p < 0.001; n.s., not significant. (D) Stress-related CREs in a region (–296 to +1) of the PAP-I promoter. Nucleotide position is indicated on the left, relative to the validated TSS (+1). Sequences of stress-related CREs, along with the CAAT and TATA boxes, are bolded and boxed. Green font indicates a known JA-associated element.
Stress-Associated CREs in the Promoters of PAP Isoforms
To gain further insight into the function of PAP isoforms, we extended our bioinformatic annotation of CREs to include the proximal promoters of all protein-coding PAP genes, including the putative alternate promoters of PAP-II and PAP-S2 (Table 5). All promoters contained TATA boxes −35 to −25 bp upstream of the putative TSS, consistent with other eukaryotic promoters. We also identified CREs associated with diverse biotic and abiotic stresses, such as T/GBOXATPIN2 (JA), W-boxes (SA), ABRE motifs (ABA), CCAAT boxes (heat stress), GARE motifs (gibberellic acid, GA), and MYB motifs (drought). Although some CREs were present in all PAP promoters (e.g., EBOXBNNAPA, GT1CONSENSUS, MYB1AT, and WBOXATNPR1), most elements differed in abundance and distribution. For instance, most ABRE motifs were absent in the PAP-α, PAP-S1, and PAP-S2A promoters, while T/GBOXATPIN2 was only present in PAP-I, PAP-α, and PAP-S2B. Other elements were unique to a single promoter, such as ELRECOREPCRP1 (PAP-IIB), GADOWNAT (PAP-IIA), and WBOXNTCHN48 (PAP-IIB). The two PAP-II promoters had only 50.5% sequence identity to each other, and the two PAP-S2 promoters had only 47.2% identity. Therefore, stimuli-dependent promoter selection in these two isoforms may lead to the transcription of distinct populations of mRNAs differing only in their 5′ UTRs. Overall, differences in the abundance and distribution of CREs in PAP promoters suggest that the isoforms have distinct roles and that control of expression can be fine-tuned by isoforms with more than one promoter.
Here, we have presented the first de novo assembled draft genome of pokeweed and an annotation of protein-coding genes. We also identified clusters of co-expressed genes by integrating RNA-Seq data from several pokeweed stress treatments. We found that PAP isoforms localized to multiple clusters, with some isoforms clustering together, and functional enrichment analysis suggested distinct biological relevance of isoforms. Validation of PAP gene models led to the discovery of a long intron within the 5′ UTR. The sequence of the intron varied for each isoform, but its presence was consistent. For PAP-I, the intron enhanced gene expression of promoter reporter constructs in tobacco. Finally, we confirmed JA-responsiveness of the PAP-I promoter in tobacco and identified a region that mediates this response. This region, as well as the proximal promoters of all PAP isoforms, contained CREs associated with stress.
Evaluation of Assembly and Annotation Metrics of the Pokeweed Genome
The pokeweed genome was sequenced exclusively as paired-end reads from a single short-insert library (83X coverage), and the resulting de novo assembly accounted for 74% of the expected size. We estimated the genome of pokeweed to be 1.3 Gb, in agreement with the previously reported value of 1.5 Gb obtained by Feulgen microdensitometry (Bennett, 2000). However, the assembly was highly fragmented (∼850,000 total scaffolds, with 70,834 scaffolds ≥ 1 Kb; N50 = 42.5 Kb). Available plant genomes have N50 values ranging from 103 to 108 bp (Veeckman et al., 2016). Nevertheless, the contiguity of the pokeweed genome assembly is comparable to other de novo assembled genomes derived from paired-end sequencing of short-insert libraries (Polashock et al., 2014; Van Hoeck et al., 2015). An assembly can also be highly complete but fragmented, as seen for a petunia species whose assembly accounted for 93% of the expected 1.1 Gb genome but had an N50 value of 17.9 Kb (Zhuang and Tripp, 2017). We acknowledge that more advanced sequencing would improve our genome assembly; however, the main goal of our draft genome was to identify PAP genes and their proximal promoters.
Annotation metrics indicate that our assembly sufficiently captured the protein-coding gene content of pokeweed. BUSCO completion scores for the genome assembly and gene set were 84% and 76%, respectively. Since the CoreGF score of the gene set was high (98%), and this metric is less stringent in terms of species conservation (Veeckman et al., 2016), we attribute missing BUSCOs more so to the divergence of pokeweed from model plants than incomplete assembly. Indeed, BUSCO scores are known to reflect both assembly contiguity and evolutionary history of the species under study (Simão et al., 2015). Furthermore, a similar BUSCO genome completion score (77%) was obtained for a non-model plant, a seagrass species, whose assembly metrics were similar to pokeweed: 71% assembled and N50 = 36.7 Kb (Lee et al., 2016). Pokeweed annotation metrics are also consistent with standards set by MAKER developers (Campbell et al., 2014a), who consider a genome to be well-annotated if 90% of its annotations have an AED less than 0.5 and over 50% of its proteome contains a recognizable protein domain. Furthermore, annotation of a plant genome with MAKER is expected to identify at least 20,000–40,000 genes. The pokeweed gene set meets these criteria: ∼30,000 genes, 99% of which have an AED less than 0.5, and 73% of the proteome contains a Pfam domain.
Relevance of Pokeweed-Specific Orthogroups in Plant Defense
Pokeweed-specific orthogroups were enriched in the GO terms ‘far-red light signaling pathway’ and ‘negative regulation of defense response,’ with the involved genes annotated as isoforms of FAR1-Related Sequence and CPR30, respectively. FAR1 is a TF involved in a variety of processes relating to growth and development, including light signal transduction, circadian clock and flowering time regulation, chlorophyll biosynthesis, starch synthesis, and ABA responses (Ma and Li, 2018). FAR1, together with the light-signaling factor Far-Red Elongated Hypocotyl 3 (FHY3), regulates plant defense by integrating chlorophyll biosynthesis and SA signaling in Arabidopsis (Wang et al., 2016). The F-box protein CPR30 also modulates defense in Arabidopsis (Gou et al., 2009). Plants deficient in CPR30 showed resistance to pathogen infection and induction of defense-related gene expression. Despite the dependence of both FAR1 and CPR30 on SA in Arabidopsis, we did not identify any significant changes in expression of these genes in any of the pokeweed stress treatments. This may reflect the fact that our stress treatments were short-term on wild-type plants and did not simulate a mutant condition. Specifically, both cpr30 and fhy3 far1 plants had dwarf phenotypes, indicative of disruption to wider pathways of growth and development. Our identification of pokeweed-specific orthogroups enables future comparison with the agricultural crop plants used in our analysis, to identify defense strategies unique to pokeweed.
Annotation of PAP Isoforms in Pokeweed
Through PAP isoform annotation, we confirmed the existence of the following previously identified isoforms: PAP-I (Irvin, 1975), PAP-II (Irvin et al., 1980), PAP-S1 (Honjo et al., 2002), PAP-S2 (Honjo et al., 2002), and PAP-α (Kataoka et al., 1992). We did not find evidence for PAP-H (Park et al., 2002), PAP-C (Barbieri et al., 1989), PAP-R (Bolognesi et al., 1990), or PAP-III (Rajamohan et al., 1999) in the scaffold assembly. However, a sequence having high identity (95%) and coverage (97%) with the cDNA clone of PAP-H was identified in the edge assembly through BLAST-N. Therefore, a gene for PAP-H likely exists in pokeweed, but it was excluded from the scaffold assembly since it was lower in coverage than an alternative path through that region. PAP-H, with 67% identity to PAP-I at the protein level, was purified from Rhizobium rhizogenes-transformed hairy roots of pokeweed (Park et al., 2002). It is secreted as part of the root exudates and hypothesized to contribute to the inhibition of soil-borne microbe infection. We did not find evidence for PAP-C or PAP-R, which were originally purified from pokeweed cell cultures and roots, respectively. Both PAP-C and PAP-R have N-terminal sequences that are identical to PAP-I; additionally, the three isoforms are highly similar in terms of amino acid composition, molecular weight (∼29 kDa), and pI value (∼9.5) (Barbieri et al., 1989; Bolognesi et al., 1990). Based on a lack of genomic support and high biochemical similarity, we suggest that reports of PAP-C, PAP-R, and PAP-I refer to the same isoform. Minor differences in biochemical properties could be explained by experimental variability, different post-translational modifications, or allelic diversity that cannot be resolved through de novo assembly. Finally, we did not identify genomic evidence in support of PAP-III, originally purified from late summer leaves. Through BLAST-P, we determined that the sequence of PAP-II, from early summer leaves, is 95% identical to the sequence of PAP-III at the protein level. Mismatches in the alignment resulted from ambiguous nucleotides (x) in the PAP-III sequence, owing to lysine methylation that enabled protein crystallization in the original report (Kurinov and Uckun, 2003). Therefore, PAP-II and PAP-III have the same amino acid sequence. It is curious that the two proteins have different levels of antiviral activity and separate as distinct peaks in ion-exchange chromatography (Rajamohan et al., 1999). We hypothesize that differences in PAP-II and PAP-III arise from post-translational modifications that affect enzymatic activity. Further support for this idea comes from the finding that a rare form of N-glycosylation exists in PAP seed isoforms, and this modification is thought to contribute to their high cytotoxicity (Islam et al., 1991; Hogg et al., 2015).
Integration of Abiotic and Biotic Stress Responses in Pokeweed
Clustering analysis enabled the identification of genes sharing a similar expression profile. We identified 36 gene clusters, 13 of which showed significant functional enrichment. While clusters associated with SA or JA were expected since these hormones have well-established roles in plant defense, we focus here on clusters revealing potential cross-talk of key biotic and abiotic stress responses in pokeweed.
The ‘cell wall integrity’ cluster included several enzymes involved in the synthesis of cellulose and lignin, which are critical components of the cell wall. Lignin creates a physical barrier against pathogens and makes plant cells more difficult for insect herbivores to penetrate and digest (Liu et al., 2018). Multiple laccases, enzymes required for lignin polymerization, were present in this cluster. A laccase was shown to mediate broad-spectrum pathogen resistance in cotton by integrating the phenylpropanoid pathway, JA biosynthesis, and balance of JA-SA defense responses (Hu et al., 2018). In pokeweed, an increase in laccase activity during Mn treatment is thought to contribute to heavy metal tolerance by reducing the level of toxic reactive oxygen species (Gao et al., 2012). Also present in this cluster were fasciclin-like arabinogalactan proteins (FLAs). These are cell surface adhesion proteins that enable cell expansion during salt stress (Shi, 2003). Consistent with their presence, the cluster expression profile showed down-regulation with PEG, indicating sensitivity to osmotic changes. Finally, the inclusion of Defective in Induced Resistance 1 (DIR1) in this cluster, an essential signaling component of systemic acquired resistance (Maldonado et al., 2002; Champigny et al., 2013), provides further indication of integrated biotic and abiotic stress responses in pokeweed.
In the ‘lipid transport’ cluster, we identified several genes encoding non-specific lipid transfer proteins (LTPs). These are small, basic, cysteine-rich proteins that localize to the apoplast and transport various lipids. LTPs are involved in the synthesis of lipid barrier polymers, such as cuticular wax, and their expression is induced by abiotic stress (Salminen et al., 2016). LTPs also appear to have a role in JA signaling. In barley, JA biosynthesis enzymes were found to produce a covalent adduct consisting of an LTP and reactive oxylipin (Bakan et al., 2006). Furthermore, exogenous application of an LTP-JA complex to grapevine produced a higher antifungal response than either component individually (Girault et al., 2008). This gene cluster had a relatively flat expression profile apart from a spike at wounding treatment, perhaps reflecting involvement of LTPs in specific activation of the wound-response branch of JA signaling (Wasternack and Song, 2017).
The ‘oxalate synthesis’ cluster included multiple genes encoding Petal Death Protein (PDP). PDP is an isocitrate lyase in senescent flower petals that catalyzes the conversion of oxaloacetate to acetate and oxalate (Lu et al., 2005). In Medicago truncatula, calcium oxalate plays a role in defense against chewing insects, with a clear feeding preference observed for plants defective in calcium oxalate production (Korth, 2006). Authors also found that calcium oxalate crystals were abrasive to insects during feeding and interfered with digestion. Oxalate, as oxalic acid, is used to chelate and detoxify heavy metals in plants (Yang et al., 2005). Notably, pokeweed has an intrinsically high oxalate content that is sufficient for chelating Mn at high concentrations (Dou et al., 2009). In this gene cluster, we also identified a heavy metal-associated isoprenylated plant protein. These metallochaperones have a known role in cadmium detoxification (Tehseen et al., 2010). With pathogen resistance and heavy metal tolerance being the most cited applications of pokeweed, our results provide insight into genes that may mediate cross-talk between these important activities.
Biological Relevance of PAP in Pokeweed
Our identification of genes co-expressed with different PAP isoforms provides a foundation for exploring their regulation in the plant and biological relevance. Two PAP-containing clusters (16 and 20) revealed significant GO term functional enrichment. Cluster 20, which contained PAP-S2 and showed down-regulation with JA, was enriched in terms relating to glycolysis and amino acid transport. The respective genes were annotated as fructose-bisphosphate aldolases (FBAs) and WAT1-related proteins. WAT1 is a vacuolar transporter that promotes indole metabolism and transport in Arabidopsis (Ranocha et al., 2010; Denancé et al., 2013). Indole is important in plant defense as a herbivore-induced volatile priming signal (Erb et al., 2015) and precursor for secondary metabolites (Lee et al., 2015). This cluster was also enriched in FBA genes, which function in glycolysis and gluconeogenesis. Co-enrichment of genes involved in carbohydrate and indole metabolism in this cluster suggests recruitment of the shikimate pathway, which provides a route to aromatic amino acid biosynthesis from chorismate (Parthasarathy et al., 2018). Specifically, the glycolysis intermediate phosphoenol pyruvate is used as input for the shikimate pathway, and indole is synthesized from tryptophan. FBA is also involved in the Calvin cycle of photosynthesis. Interestingly, this cluster included genes encoding pentatricopeptide repeat-containing (PPR) proteins that participate in RNA editing in the chloroplast (PCMP-H61, PCMP-H51/CRR28), an activity required for accumulation of the NADH dehydrogenase-like complex of the photosynthetic electron transport chain (Okuda et al., 2007, 2009). PPRs are a large family of sequence-specific RNA-binding proteins involved in multiple aspects of RNA metabolism. Taken together, these results demonstrate PAP-S2 co-expression with genes involved in key processes of carbohydrate and amino acid metabolism, including other proteins that interact directly with RNA.
Cluster 16, which included PAP-S1 and PAP-α, contained JA-upregulated genes and had corresponding GO enrichment in JA signaling and terpenoid biosynthesis. Like Cluster 20, this cluster also indicated PAP co-expression with genes involved in metabolism, specifically with the presence of BZIP11, a sucrose-regulated TF that controls amino acid and carbohydrate metabolism and is integrated with a wider plant growth regulatory network (Hanson et al., 2008; Ma et al., 2011). Another gene in this cluster, annotated as Enhanced Downy Mildew 2 (EDM2), is an interesting candidate that links regulation of gene expression with long introns and pathogen defense. EDM2 is a chromatin regulator that promotes normal 3′ distal polyadenylation of transcripts from genes containing intronic heterochromatin, the latter arising from methylated transposons or repeats associated with long introns (Lei et al., 2013; Duan et al., 2017). EDM2 is also required for pathogen resistance in Arabidopsis by regulating transcript accumulation of an NB-LRR disease resistance (R) gene (Eulgem et al., 2007).
Clusters 6 and 17, which contained PAP-I and PAP-II, respectively, did not have significant functional enrichment (FDR < 0.05). However, assignment of these isoforms to independent clusters reinforces their distinct regulation in the plant. The PAP-II cluster contained 92 genes and was enriched in the GO terms ‘wounding’ (FDR = 0.10) and ‘JA biosynthetic process’ (FDR = 0.23). The PAP-I cluster contained only nine genes, including two TFs from the HD-ZIP homeobox family that regulate plant growth and leaf development in response to abiotic stresses (Söderman et al., 1996; Aoyama et al., 2007). Other genes in the PAP-I cluster support an integration with plant growth responses. This includes DAD2, an esterase that mediates strigolactone signaling (Hamiaux et al., 2012). Strigolactone is a hormone involved in branching and symbiotic interactions with soil microbes (Marzec, 2016). Also present is a gene annotated as DHNAT1, encoding an enzyme (1,4-dihydroxy-2-naphthoyl-CoA thioesterase 1) involved in the synthesis of phylloquinone (Widhalm et al., 2012). Phylloquinone is required for Photosystem 1 stability (Wang et al., 2017) and is synthesized from chorismate, the final product of the shikimate pathway described above. These findings support the idea that PAP expression regulation is tied to the broader metabolic state of the plant. Specifically, our results indicate co-expression of PAP isoforms with genes involved in amino acid and carbohydrate metabolism, suggesting a link to wider nutrient-sensing networks. Therefore, the ribosome-inactivating activity of PAP may contribute to and/or be affected by large-scale reprograming in response to plant stress. This hypothesis is strengthened given the direct antiviral activity of PAP, which could potentially mediate a trade-off between plant growth and defense in pokeweed.
Mechanisms of Gene Expression Regulation by Leader Introns
Investigation of PAP gene models led to the discovery of a long intron in the 5′ UTR of all PAP isoforms. The sequence of the intron was different for each isoform, but its presence was conserved. We also provided evidence to support a functional role for the intron, as the PAP-I intron enhanced reporter gene expression in tobacco. Introns can influence gene expression in both plants and animals, particularly leader introns (Laxa, 2017; Shaul, 2017). The precise mechanism of intron-mediated enhancement is not known, but hypotheses at both the transcriptional and translational levels have been proposed. Leader introns may enhance transcription by creating a favorable zone for transcription initiation; that is, a region: (i) devoid of nucleosomes, (ii) marked by activating histone modifications, or (iii) associated with a factor that recruits transcriptional machinery (Gallegos and Rose, 2015). Additionally, splicing signals in the leader intron may affect mRNA processing, export, or decay, thereby affecting translation (Gallegos and Rose, 2015; Laxa, 2017). In contrast to intron-mediated enhancement, which requires the intron to be in its native orientation and position, a leader intron may act in other ways to enhance expression. For example, the intron may function as a classical enhancer by containing CREs (Kim et al., 2006) or as an alternate promoter (Morello et al., 2002, 2006; Qi et al., 2007).
Our preliminary characterization allowed us to conclude that the PAP-I leader intron enhanced expression. Interestingly, we observed that enhancement was greater when the intron was paired with the minimal PAP-I promoter than with the proximal promoter. There may be a limit to intron-associated enhancement of gene expression, particularly for promoters that drive high expression on their own. For example, when paired with its native weak promoter, the leader intron of Arabidopsis AtMHX increased expression by 270-fold, compared to 3-fold with the strong CaMV 35S promoter (Akua et al., 2010). Alternatively, the enhancement provided by an intron may simply be more detectable in the presence of a weak promoter. Future work will investigate the mechanism by which the leader intron enhances PAP expression.
Transcriptional Control of PAP Isoform Expression
The promoters of all PAP isoforms contained CREs associated with diverse abiotic and biotic stresses, suggesting that PAP is broadly implicated in plant defense. Since PAP genes were most responsive to JA in our study, we aimed to identify CREs that could mediate this response. Promoter truncation constructs of PAP-I revealed that a region close to the TSS (−296 to −103) was sufficient for JA-responsiveness. This result agrees with a previous finding that CREs closer to the TSS tend to have a greater influence on transcription (Zou et al., 2011). The −296 to −103 region of PAP-I contains binding motifs for TFs from bHLH, bZIP, and MYB families, which were found to be overrepresented in JA-responsive promoters (Hickman et al., 2017). Additionally, this region contains a T/GBOXATPIN2 element, which binds the master JA signaling regulator MYC (Boter et al., 2004). Mutation of this element abolished JA-responsiveness of genes in tomato, Arabidopsis, and barley (Rouster et al., 1997; Boter et al., 2004).
Analysis of JA-associated CREs in the promoters of other PAP isoforms revealed that the T/GBOXATPIN2 element was also present in PAP-α and PAP-S2B, but not in PAP-S1, which was most JA-responsive. The lack of this element in some promoters of JA-responsive isoforms may be compensated by the presence of W-boxes; this element binds WRKY TFs, which are primarily SA-responsive (Dong et al., 2003). However, a substantial fraction of WRKYs (at least 30%) are JA-responsive in Arabidopsis (Schluttenhofer et al., 2014), and we identified several that were up-regulated with JA in pokeweed (homologs of WRKY3, WRKY4, WRKY22, WRKY23, WRKY24, WRKY33, WRKY40, WRKY41, WRKY49, WRKY70, and WRKY75). Combined with the knowledge that synergistic activation of gene expression can occur in the presence of both SA and JA (Mur, 2005), we hypothesize that W-boxes in the promoters of PAP isoforms contribute to their JA-responsiveness. Although the promoter of PAP-S2 contained several W-boxes, PAP-S2 expression decreased slightly but significantly with JA and was unresponsive to other treatments. It is possible that PAP-S2 has a different temporal expression profile than the other isoforms and is more responsive outside of the 24h time-point we investigated. Consistent with this hypothesis, an RNA-Seq time-course of the JA response in Arabidopsis revealed diverse expression patterns over the first 16 h following treatment, including distinct early and late responses (Hickman et al., 2017). In addition to CREs associated with JA, we identified CREs associated with the hormones SA, ABA, and GA. JA and SA have well-established roles in plant defense against pathogens and insect herbivores, while ABA contributes to the resistance of abiotic stresses such as drought, salinity, cold, and heat stress (Verma et al., 2016). GA, through cross-talk with ABA pathways, helps mediate the balance between dormancy and plant maturation during stress (Verma et al., 2016). Importantly, we identified differences in the CREs of PAP promoters, suggesting that the isoforms mediate individual responses to hormones.
The draft genome of pokeweed has provided new information on how PAP expression is controlled. The presence of diverse stress-responsive CREs in the promoters of PAP genes, combined with their distinct expression profiles, provides a foundation to examine the role of PAP isoforms in pokeweed. One-fifth of land plants synthesize RIPs, and these genes are often up-regulated during stress. Our study contributes to the growing body of evidence illuminating RIPs as important components of plant response to environmental change.
All datasets for this study are included in the manuscript and/or the Supplementary Files.
KN and KH designed the study. KN performed the bioinformatic analyses including genome assembly, annotation, and RNA-Seq analysis, and drafted the manuscript. CD performed all cloning, gene model and RNA-seq validations, measurement of gene expression from reporter constructs, and identification of promoter elements. AP provided advice on bioinformatic analyses and data interpretation. KH edited the manuscript. All authors read and approved the final manuscript.
This work was funded by a Discovery Grant to KH from the Natural Sciences and Engineering Research Council of Canada, and PGS-D and CGS-M Scholarships to KN and CD.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors are grateful to Dr. Paul Kron, Department of Integrative Biology, University of Guelph, for genome size analysis of pokeweed.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.01002/full#supplementary-material
ABA, abscisic acid; AED, annotation edit distance; CRE, cis-regulatory elements; DEG, differentially expressed gene; ET, ethanol; FC, fold change; FDR, false discovery rate; GA, gibberellic acid; GO, gene ontology; GUS, β-glucuronidase; JA, jasmonic acid; PAP, pokeweed antiviral protein; PEG, polyethylene glycol; RIP, ribosome-inactivating protein; SA, salicylic acid; TF, transcription factor; TPM, transcripts per million; TSS, transcription start site; UTR, untranslated region; WND, wounded with forceps; WT, watered normally.
- ^ https://www.sharcnet.ca
- ^ http://www.repeatmasker.org/RepeatModeler/
- ^ http://www.repeatmasker.org/
- ^ https://phytozome.jgi.doe.gov
- ^ http://www.spinachbase.org
- ^ http://bvseq.boku.ac.at
- ^ http://www.ebi.ac.uk/Tools/msa/clustalo/
- ^ http://www.dna.affrc.go.jp/htdocs/PLACE/
- ^ http://plantpan2.itps.ncku.edu.tw/
Abe, H., Urao, T., Ito, T., Seki, M., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2003). Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell. 15, 63–78. doi: 10.1105/tpc.006130
Agarwal, M., Hao, Y., Kapoor, A., Dong, C. H., Fujii, H., Zheng, X., et al. (2006). A R2R3 type MYB transcription factor is involved in the cold regulation of CBF genes and in acquired freezing tolerance. J. Biol. Chem. 281, 37636–37645. doi: 10.1074/jbc.M605895200
Akua, T., Berezin, I., and Shaul, O. (2010). The leader intron of AtMHX can elicit, in the absence of splicing, low-level intron-mediated enhancement that depends on the internal intron sequence. BMC Plant Biol. 10:93. doi: 10.1186/1471-2229-10-93
Alison Dunn, M., White, A. J., Vural, S., and Hughes, M. A. (1998). Identification of promoter elements in a low-temperature-responsive gene (blt4.9) from barley (Hordeum vulgare L.). Plant Mol. Biol. 38, 551–564. doi: 10.1023/A:1006098132352
Aoto, Y., Hachiya, T., Okumura, K., Hase, S., Sato, K., Wakabayashi, Y., et al. (2017). DEclust: a statistical approach for obtaining differential expression profiles of multiple conditions. PLoS One 12:e0188285. doi: 10.1371/journal.pone.0188285
Aoyama, T., Dong, C.-H., Wu, Y., Carabelli, M., Sessa, G., Ruberti, I., et al. (2007). Ectopic expression of the arabidopsis transcriptional activator athb-1 alters leaf cell fate in tobacco. Plant Cell 7:1773. doi: 10.2307/3870186
Bakan, B., Hamberg, M., Perrocheau, L., Maume, D., Rogniaux, H., Tranquet, O., et al. (2006). Specific adduction of plant lipid transfer protein by an allene oxide generated by 9-lipoxygenase and allene oxide synthase. J. Biol. Chem. 281, 38981–38988. doi: 10.1074/jbc.M608580200
Barbieri, L., Aron, G. M., Irvin, J. D., and Stirpe, F. (1982). Purification and partial characterization of another form of the antiviral protein from the seeds of Phytolacca americana L. (pokeweed). Biochem. J. 203, 55–59. doi: 10.1042/bj2030055
Barbieri, L., Bolognesi, A., Cenini, P., Falasca, A. I., Minghetti, A., Garofano, L., et al. (1989). Ribosome-inactivating proteins from plant cells in culture. Biochem. J. 257, 801–807. doi: 10.1042/bj2570801
Bolognesi, A., Barbieri, L., Abbondanza, A., Falasca, A. I., Carnicelli, D., Battelli, M. G., et al. (1990). Purification and properties of new ribosome-inactivating proteins with RNA N-glycosidase activity. Biochim. Biophys. Acta 1087, 293–302. doi: 10.1016/0167-4781(90)90002-J
Bonness, M. S., Ready, M. P., Irvin, J. D., and Mabry, T. J. (1994). Pokeweed antiviral protein inactivates pokeweed ribosomes; implications for the antiviral mechanism. Plant J. 5, 173–183. doi: 10.1046/j.1365-313X.1994.05020173.x
Boter, M., Ruíz-Rivero, O., Abdeen, A., and Prat, S. (2004). Conserved MYC transcription factors play a key role in jasmonate signaling both in tomato and Arabidopsis. Genes Dev. 18, 1577–1591. doi: 10.1101/gad.297704
Buchel, A. S., Brederode, F. T., Bol, J. F., and Linthorst, H. J. M. (1999). Mutation of GT-1 binding sites in the Pr-1a promoter influences the level of inducible gene expression in vivo. Plant Mol. Biol. 40, 387–396. doi: 10.1023/A:1006144505121
Campbell, M. S., Holt, C., Moore, B., and Yandell, M. (2014a). Genome annotation and curation using MAKER and MAKER-P. Curr. Protoc. Bioinforma. 2014, 4.11.1–4.11.39. doi: 10.1002/0471250953.bi0411s48
Campbell, M. S., Law, M., Holt, C., Stein, J. C., Moghe, G. D., Hufnagel, D. E., et al. (2014b). MAKER-P: a tool kit for the rapid creation, management, and quality control of plant genome annotations. Plant Physiol. 164, 513–524. doi: 10.1104/pp.113.230144
Cercós, M., Gómez-Cadenas, A., and Ho, T. H. D. (1999). Hormonal regulation of a cysteine proteinase gene, EPB-1, in barley aleurone layers: Cis- and trans-acting elements involved in the co-ordinated gene expression regulated by gibberellins and abscisic acid. Plant J. 19, 107–118. doi: 10.1046/j.1365-313X.1999.00499.x
Chakravarthy, S., Tuori, R., D’Ascenzo, M., Fobert, P., Despres, C., and Martin, G. (2003). The tomato transcription factor Pti4 Regulates defense-related gene expression via GCC Box and Non-GCC box cis elements. Plant Cell. 15, 3033–3050. doi: 10.1105/tpc.017574
Champigny, M. J., Isaacs, M., Carella, P., Faubert, J., Fobert, P. R., and Cameron, R. K. (2013). Long distance movement of DIR1 and investigation of the role of DIR1-like during systemic acquired resistance in Arabidopsis. Front. Plant Sci. 4:230. doi: 10.3389/fpls.2013.00230
Chow, C. N., Zheng, H. Q., Wu, N. Y., Chien, C. H., Huang, H., Da Lee, T. Y., et al. (2016). PlantPAN 2.0: an update of Plant Promoter Analysis Navigator for reconstructing transcriptional regulatory networks in plants. Nucleic Acids Res. 44, D1154–D1160. doi: 10.1093/nar/gkv1035
Clouse, J. W., Adhikary, D., Page, J. T., Ramaraj, T., Deyholos, M. K., Udall, J. A., et al. (2016). The amaranth genome: genome, transcriptome, and physical map assembly. Plant Genome 9. doi: 10.3835/plantgenome2015.07.0062
Côté, C., and Rutledge, R. G. (2003). An improved MUG fluorescent assay for the determination of GUS activity within transgenic tissue of woody plants. Plant Cell Rep. 21, 619–624. doi: 10.1007/s00299-002-0543-z
Dai, W. D., Bonos, S., Guo, Z., Meyer, W. A., Day, P. R., and Belanger, F. C. (2003). Expression of pokeweed antiviral proteins in creeping bentgrass. Plant Cell Rep. 21, 497–502. doi: 10.1007/s00299-002-0534-530
Denancé, N., Ranocha, P., Oria, N., Barlet, X., Rivière, M. P., Yadeta, K. A., et al. (2013). Arabidopsis wat1 (walls are thin1)-mediated resistance to the bacterial vascular pathogen, Ralstonia solanacearum, is accompanied by cross-regulation of salicylic acid and tryptophan metabolism. Plant J. 73, 225–239. doi: 10.1111/tpj.12027
Di Maro, A., Citores, L., Russo, R., Iglesias, R., and Ferreras, J. M. (2014). Sequence comparison and phylogenetic analysis by the maximum likelihood method of ribosome-inactivating proteins from angiosperms. Plant Mol. Biol. 85, 575–588. doi: 10.1007/s11103-014-0204-y
Dohm, J. C., Minoche, A. E., Holtgräwe, D., Capella-Gutiérrez, S., Zakrzewski, F., Tafer, H., et al. (2014). The genome of the recently domesticated crop plant sugar beet (Beta vulgaris). Nature 505, 546–549. doi: 10.1038/nature12817
Dou, C. M., Fu, X. P., Chen, X. C., Shi, J. Y., and Chen, Y. X. (2009). Accumulation and detoxification of manganese in hyperaccumulator Phytolacca americana. Plant Biol. 11, 664–670. doi: 10.1111/j.1438-8677.2008.00163.x
Duan, C.-G., Wang, X., Zhang, L., Xiong, X., Zhang, Z., Tang, K., et al. (2017). A protein complex regulates RNA processing of intronic heterochromatin-containing genes in Arabidopsis. Proc. Natl. Acad. Sci. 114, E7377–E7384. doi: 10.1073/pnas.1710683114
Dubouzet, J. G., Sakuma, Y., Ito, Y., Kasuga, M., Dubouzet, E. G., Miura, S., et al. (2003). OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt- and cold-responsive gene expression. Plant J. 428, 75–79. doi: 10.1046/j.1365-313X.2003.01661.x
Emms, D. M., and Kelly, S. (2015). OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 16:157. doi: 10.1186/s13059-015-0721-722
Endo, Y., Tsurugi, K., and Lambert, J. M. (1988). The site of action of six different ribosome-inactivating proteins from plants on eukaryotic ribosomes: the RNA N-glycosidase activity of the proteins. Biochem. Biophys. Res. Commun. 150, 1032–1036. doi: 10.1016/0006-291x(88)90733-4
Erb, M., Veyrat, N., Robert, C. A. M., Xu, H., Frey, M., Ton, J., et al. (2015). Indole is an essential herbivore-induced volatile priming signal in maize. Nat. Commun. 6:6273. doi: 10.1038/ncomms7273
Eulgem, T., Tsuchiya, T., Wang, X. J., Beasley, B., Cuzick, A., Tör, M., et al. (2007). EDM2 is required for RPP7-dependent disease resistance in Arabidopsis and affects RPP7 transcript levels. Plant J. 49, 829–839. doi: 10.1111/j.1365-313X.2006.02999.x
Ezcurra, I., Ellerström, M., Wycliffe, P., Stålberg, K., and Rask, L. (1999). Interaction between composite elements in the napA promoter: both the B-box ABA-responsive complex and the RY/G complex are necessary for seed-specific expression. Plant Mol. Biol. 40, 699–709. doi: 10.1023/A:1006206124512
Gao, L., Peng, K., Chen, Y., Wang, G., and Shen, Z. (2012). Roles of apoplastic peroxidases, laccases, and lignification in the manganese tolerance of hyperaccumulator Phytolacca americana. Acta Physiol. Plant. 34, 151–159. doi: 10.1007/s11738-011-0813-x
Girault, T., François, J., Rogniaux, H., Pascal, S., Delrot, S., Coutos-Thévenot, P., et al. (2008). Exogenous application of a lipid transfer protein-jasmonic acid complex induces protection of grapevine towards infection by Botrytis cinerea. Plant Physiol. Biochem. 46, 140–149. doi: 10.1016/j.plaphy.2007.10.005
Gou, M., Su, N., Zheng, J., Huai, J., Wu, G., Zhao, J., et al. (2009). An F-box gene, CPR30, functions as a negative regulator of the defense response in Arabidopsis. Plant J. 60, 757–770. doi: 10.1111/j.1365-313X.2009.03995.x
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Gu, Y. Q., Wildermuth, M. C., Chakravarthy, S., Loh, Y. T., Yang, C., He, X., et al. (2002). Tomato transcription factors pti4, pti5, and pti6 activate defense responses when expressed in Arabidopsis. Plant Cell. 14, 817–831. doi: 10.1105/tpc.000794
Gubler, F., Kalla, R., Roberts, J. K., and Jacobsen, J. V. (1995). Gibberellin-regulated expression of a myb gene in barley aleurone cells: evidence for myb transactivation of a high-pl alpha-amylase gene promoter. Plant Cell. 7, 1879–1891. doi: 10.1105/tpc.7.11.1879
Haas, B. J., Papanicolaou, A., Yassour, M., Grabherr, M., Blood, P. D., Bowden, J., et al. (2013). De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494–1512. doi: 10.1038/nprot.2013.084
Hamiaux, C., Drummond, R. S. M., Janssen, B. J., Ledger, S. E., Cooney, J. M., Newcomb, R. D., et al. (2012). DAD2 is an α/β hydrolase likely to be involved in the perception of the plant branching hormone, strigolactone. Curr. Biol. 22, 2032–2036. doi: 10.1016/j.cub.2012.08.007
Hanson, J., Hanssen, M., Wiese, A., Hendriks, M. M. W. B., and Smeekens, S. (2008). The sucrose regulated transcription factor bZIP11 affects amino acid metabolism by regulating the expression of ASPARAGINE SYNTHETASE1 and PROLINE DEHYDROGENASE2. Plant J. 53, 935–949. doi: 10.1111/j.1365-313X.2007.03385.x
He, Y.-W., Guo, C.-X., Pan, Y.-F., Peng, C., and Weng, Z.-H. (2008). Inhibition of hepatitis B virus replication by pokeweed antiviral protein in vitro. World J. Gastroenterol. 14, 1592–1597. doi: 10.3748/wjg.14.1592
Healey, A., Furtado, A., Cooper, T., and Henry, R. J. (2014). Protocol: a simple method for extracting next-generation sequencing quality genomic DNA from recalcitrant plant species. Plant Methods 10:21. doi: 10.1186/1746-4811-10-21
Hickman, R., van Verk, M. C., Van Dijken, A. J. H., Pereira Mendes, M., Vroegop-Vos, I. A., Caarls, L., et al. (2017). Architecture and dynamics of the jasmonic acid gene regulatory network. Plant Cell 29, 2086–2105. doi: 10.1105/tpc.16.00958
Hobo, T., Asada, M., Kowyama, Y., and Hattori, T. (1999). ACGT-containing abscisic acid response element (ABRE) and coupling element 3 (CE3) are functionally equivalent. Plant J. 19, 679–689. doi: 10.1046/j.1365-313x.1999.00565.x
Hogg, T., Mendel, J. T., and Lavezo, J. L. (2015). Structural analysis of a type 1 ribosome inactivating protein reveals multiple L-asparagine-N-acetyl-D-glucosamine monosaccharide modifications: implications for cytotoxicity. Mol. Med. Rep. 12, 5737–5745. doi: 10.3892/mmr.2015.4146
Honjo, E., Dong, D., Motoshima, H., and Watanabe, K. (2002). Genomic clones encoding two isoforms of pokeweed antiviral protein in seeds (PAP-S1 and S2) and the N-Glycosidase activities their recombinant proteins on ribosomes and DNA in comparison with other isoforms. J. Biochem. 231, 225–231. doi: 10.1093/oxfordjournals.jbchem.a003092
Hu, Q., Min, L., Yang, X., Jin, S., Zhang, L., Li, Y., et al. (2018). Laccase GhLac1 modulates broad-spectrum biotic stress tolerance via manipulating phenylpropanoid pathway and jasmonic acid synthesis. Plant Physiol. 176, 1808–1823. doi: 10.1104/pp.17.01628
Iglesias, R., Pérez, Y., De Torre, C., Ferreras, J. M., Antolín, P., Jiménez, P., et al. (2005). Molecular characterization and systemic induction of single-chain ribosome-inactivating proteins (RIPs) in sugar beet (Beta vulgaris) leaves. J. Exp. Bot. 56, 1675–1684. doi: 10.1093/jxb/eri164
Irvin, J. (1975). Purification and partial characterization of the antiviral protein from Phytolacca americana which inhibits eukaryotic protein synthesis. Arch. Biochem. Biophys. 169, 522–528. doi: 10.1016/0003-9861(75)90195-2
Irvin, J. D., Kelly, T., and Robertus, J. D. (1980). Purification and properties of a second antiviral protein from Phytolacca americana which inactivates eukaryotic ribosomes. Arch. Biochem. Biophys. 200, 418–425. doi: 10.1016/0003-9861(80)90372-0
Islam, M. R., Kung, S. S., Kimura, Y., and Funatsu, G. (1991). N-acetyl-D-glucosamine-asparagine structure in ribosome-inactivating proteins from the seeds of Luffa cylindrica and Phytolacca americana. Agric. Biol. Chem. 55, 1375–1381. doi: 10.1080/00021369.1991.10870763
Jefferson, R. A., Kavanagh, T. A., and Bevan, M. W. (1987). GUS fusions: beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants. EMBO J. 6, 3901–3907. doi: 10.1002/j.1460-2075.1987.tb02730.x
Jeong, J. Y., Yim, H. S., Ryu, J. Y., Lee, H. S., Lee, J. H., Seen, D. S., et al. (2012). One-step sequence-and ligation-independent cloning as a rapid and versatile cloning method for functional genomics studies. Appl. Environ. Microbiol. 78, 5440–5443. doi: 10.1128/AEM.00844-12
Jiang, S. Y., Ramamoorthy, R., Bhalla, R., Luan, H. F., Venkatesh, P. N., Cai, M., et al. (2008). Genome-wide survey of the RIP domain family in Oryza sativa and their expression profiles under various abiotic and biotic stresses. Plant Mol. Biol. 67, 603–614. doi: 10.1007/s11103-008-9342-4
Johnston, J. S., Bennett, M. D., Rayburn, A. L., Galbraith, D. W., and Price, H. J. (1999). Reference standards for determination of DNA content of plant nuclei. Am. J. Bot. 86, 609–613. doi: 10.2307/2656569
Kaplan, B., Davydov, O., Knight, H., Galon, Y., Knight, M. R., Fluhr, R., et al. (2006). Rapid Transcriptome changes induced by Cytosolic Ca2+ transients reveal ABRE-related sequences as Ca2+-responsive cis elements in arabidopsis. Plant Cell. 18, 2733–2748. doi: 10.1105/tpc.106.042713
Karran, R. A., and Hudak, K. A. (2008). Depurination within the intergenic region of Brome mosaic virus RNA3 inhibits viral replication in vitro and in vivo. Nucleic Acids Res. 36, 7230–7239. doi: 10.1093/nar/gkn896
Kataoka, J., Habuka, N., Masuta, C., Miyano, M., and Koiwai, A. (1992). Isolation and analysis of a genomic clone encoding a pokeweed antiviral protein. Plant Mol. Biol. 20, 879–886. doi: 10.1007/bf00027159
Kim, H. J., Kim, Y. K., Park, J. Y., and Kim, J. (2002). Light signalling mediated by phytochrome plays an important role in cold-induced gene expression through the C-repeat/dehydration responsive element (C/DRE) in Arabidopsis thaliana. Plant J. 29, 693–704. doi: 10.1046/j.1365-313x.2002.01249.x
Kim, M. J., Kim, H., Shin, J. S., Chung, C. H., Ohlrogge, J. B., and Suh, M. C. (2006). Seed-specific expression of sesame microsomal oleic acid desaturase is controlled by combinatorial properties between negative cis-regulatory elements in the SeFAD2 promoter and enhancers in the 5’-UTR intron. Mol. Genet. Genomics 276, 351–368. doi: 10.1007/s00438-006-0148-2
Kim, S. Y., Chung, H. J., and Thomas, T. L. (1997). Isolation of a novel class of bZIP transcription factors that interact with ABA-responsive and embryo-specification elements in the Dc3 promoter using a modified yeast one-hybrid system. Plant J. 11, 1237–1251. doi: 10.1046/j.1365-313X.1997.11061237.x
Kizis, D., and Pagès, M. (2002). Maize DRE-binding proteins DBF1 and DBF2 are involved in rab17 regulation through the drought-responsive element in an ABA-dependent pathway. Plant J. 30, 679–689. doi: 10.1046/j.1365-313X.2002.01325.x
Korth, K. L. (2006). Medicago truncatula mutants demonstrate the role of plant calcium oxalate crystals as an effective defense against chewing insects. Plant Physiol. 141, 188–195. doi: 10.1104/pp.106.076737
Laloi, C., Mestres-Ortega, D., Marco, Y., Meyer, Y., and Reichheld, J. (2004). The arabidopsis cytosolic thioredoxin h5 gene induction by oxidative stress and its w-box-mediated response to pathogen elicitor. Plant Physiol. 134, 1006–1016. doi: 10.1104/pp.103.035782
Lee, H., Golicz, A. A., Bayer, P. E., Jiao, Y., Tang, H., Paterson, A. H., et al. (2016). The genome of a southern hemisphere seagrass species (Zostera muelleri). Plant Physiol. 172, 272–283. doi: 10.1104/pp.16.00868
Lei, M., La, H., Lu, K., Wang, P., Miki, D., Ren, Z., et al. (2013). Arabidopsis EDM2 promotes IBM1 distal polyadenylation and regulates genome DNA methylation patterns. Proc. Natl. Acad. Sci. 111, 527–532. doi: 10.1073/pnas.1320106110
Liu, X., Peng, K., Wang, A., Lian, C., and Shen, Z. (2010). Cadmium accumulation and distribution in populations of Phytolacca americana L. and the role of transpiration. Chemosphere 78, 1136–1141. doi: 10.1016/j.chemosphere.2009.12.030
Lodge, J. K., Kaniewski, W. K., and Tumer, N. E. (1993). Broad-spectrum virus resistance in transgenic plants expressing pokeweed antiviral protein. Proc. Natl. Acad. Sci. U. S. A. 90, 7089–7093. doi: 10.1073/pnas.90.15.7089
Love, R. R., Weisenfeld, N. I., Jaffe, D. B., Besansky, N. J., and Neafsey, D. E. (2016). Evaluation of DISCOVAR de novo using a mosquito sample for cost-effective short-read genome assembly. BMC Genomics 17:187. doi: 10.1186/s12864-016-2531-7
Lu, Z., Feng, X., Song, L., Han, Y., Kim, A., Herzberg, O., et al. (2005). Diversity of function in the isocitrate lyase enzyme superfamily: the dianthus caryophyllus petal death protein cleaves α-keto and α-hydroxycarboxylic acids. Biochemistry 44, 16365–16376. doi: 10.1021/bi051776l
Luo, H., Song, F., Goodman, R. M., and Zheng, Z. (2005). Up-regulation of OsBIHD1, a rice gene encoding BELL homeodomain transcriptional factor, in disease resistance responses. Plant Biol. 7, 459–468. doi: 10.1055/s-2005-865851
Ma, J., Hanssen, M., Lundgren, K., Hernández, L., Delatte, T., Ehlert, A., et al. (2011). The sucrose-regulated Arabidopsis transcription factor bZIP11 reprograms metabolism and regulates trehalose metabolism. New Phytol. 191, 733–745. doi: 10.1111/j.1469-8137.2011.03735.x
Maldonado, A. M., Doerner, P., Dixonk, R. A., Lamb, C. J., and Cameron, R. K. (2002). A putative lipid transfer protein involved in systemic resistance signalling in Arabidopsis. Nature 419, 399–403. doi: 10.1038/nature00962
Mansouri, S., Choudhary, G., Sarzala, P. M., Ratner, L., and Hudak, K. A. (2009). Suppression of human T-cell leukemia virus I gene expression by pokeweed antiviral protein. J. Biol. Chem. 284, 31453–31462. doi: 10.1074/jbc.M109.046235
Maruyama-Nakashita, A., Nakamura, Y., Watanabe-Takahashi, A., Inoue, E., Yamaya, T., and Takahashi, H. (2005). Identification of a novel cis-acting element conferring sulfur deficiency response in Arabidopsis roots. Plant J. 42, 305–314. doi: 10.1111/j.1365-313X.2005.02363.x
Mena, M., Cejudo, F., Isabel-Lamoneda, I., and Carbonero, P. (2002). A role for the dof transcription factor bpbf in the regulation of gibberellin-responsive genes in barley aleurone. Plant Physiol. 130, 111–119. doi: 10.1104/pp.005561
Mohanty, B., Krishnan, S. P. T., Swarup, S., and Bajic, V. B. (2005). Detection and preliminary analysis of motifs in promoters of anaerobically induced genes of different plant species. Ann. Bot. 96, 669–681. doi: 10.1093/aob/mci219
Morello, L., Bardini, M., Cricrì, M., Sala, F., and Breviario, D. (2006). Functional analysis of DNA sequences controlling the expression of the rice OsCDPK2 gene. Planta 223, 479–491. doi: 10.1007/s00425-005-0105-z
Morello, L., Bardini, M., Sala, F., and Breviario, D. (2002). A long leader intron of the Ostub16 rice beta-tubulin gene is required for high-level gene expression and can autonomously promote transcription both in vivo and in vitro. Plant J. 29, 33–44. doi: 10.1046/j.0960-7412.2001.01192.x
Mur, L. A. J. (2005). The outcomes of concentration-specific interactions between salicylate and jasmonate signaling include synergy, antagonism, and oxidative stress leading to cell death. Plant Physiol. 140, 249–262. doi: 10.1104/pp.105.072348
Nishiuchi, T., Shinshi, H., and Suzuki, K. (2004). Rapid and transient activation of transcription of the ERF3 gene by wounding in tobacco leaves: possible involvement of NtWRKYs and autorepression. J. Biol. Chem. 79, 55355–55361. doi: 10.1074/jbc.M409674200
Ogawa, M., Hanada, A., Yamauchi, Y., Kuwahara, A., Kamiya, Y., and Yamaguchi, S. (2003). Gibberellin biosynthesis and response during arabidopsis seed germination. Plant Cell. 15, 1591–1604. doi: 10.1105/tpc.011650
Okuda, K., Chateigner-Boutin, A.-L., Nakamura, T., Delannoy, E., Sugita, M., Myouga, F., et al. (2009). Pentatricopeptide repeat proteins with the DYW motif have distinct molecular functions in RNA editing and RNA cleavage in arabidopsis chloroplasts. Plant Cell 21, 146–156. doi: 10.1105/tpc.108.064667
Okuda, K., Myouga, F., Motohashi, R., Shinozaki, K., and Shikanai, T. (2007). Conserved domain structure of pentatricopeptide repeat proteins involved in chloroplast RNA editing. Proc. Natl. Acad. Sci. U. S. A. 104, 8178–8183. doi: 10.1073/pnas.0700865104
Park, H., Kim, M., Kang, Y., Jeon, J., Yoo, J., Kim, M., et al. (2004). Pathogen- and NaCl-induced expression of the SCaM-4 promoter is mediated in part by a GT-1 box that interacts with a GT-1-like transcription factor. Plant Physiol. 135, 2150–2161. doi: 10.1104/pp.104.041442
Park, S., Lawrence, C. B., Linden, J. C., and Vivanco, J. M. (2002). Isolation and characterization of a novel ribosome- inactivating protein from root cultures of pokeweed and its mechanism of secretion from roots 1. Society 130, 164–178. doi: 10.1104/pp.000794.of
Parthasarathy, A., Cross, P. J., Dobson, R. C. J., Adams, L. E., Savka, M. A., and Hudson, A. O. (2018). A three-ring circus: metabolism of the three proteogenic aromatic amino acids and their role in the health of plants and animals. Front. Mol. Biosci. 5:29. doi: 10.3389/fmolb.2018.00029
Peng, K., Luo, C., You, W., Lian, C., Li, X., and Shen, Z. (2008). Manganese uptake and interactions with cadmium in the hyperaccumulator-Phytolacca Americana L. J. Hazard. Mater. 154, 674–681. doi: 10.1016/j.jhazmat.2007.10.080
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi: 10.1038/nbt.3122
Polashock, J., Zelzion, E., Fajardo, D., Zalapa, J., Georgi, L., Bhattacharya, D., et al. (2014). The American cranberry: first insights into the whole genome of a species adapted to bog habitat. BMC Plant Biol. 14:165. doi: 10.1186/1471-2229-14-165
Qi, X. T., Zhang, Y. X., and Chai, T. Y. (2007). The bean PvSR2 gene produces two transcripts by alternative promoter usage. Biochem. Biophys. Res. Commun. 356, 273–278. doi: 10.1016/j.bbrc.2007.02.124
Qin, X., Zheng, X., Shao, C., Gao, J., Jiang, L., Zhu, X., et al. (2009). Stress-induced curcin-L promoter in leaves of Jatropha curcas L. and characterization in transgenic tobacco. Planta 230, 387–395. doi: 10.1007/s00425-009-0956-9
Rajamohan, F., Venkatachalam, T. K., Irvin, J. D., and Uckun, F. M. (1999). Pokeweed antiviral protein isoforms PAP-I, PAP-II, and PAP-III depurinate RNA of human immunodeficiency virus (HIV)-1. Biochem. Biophys. Res. Commun. 260, 453–458. doi: 10.1006/bbrc.1999.0922
Ranocha, P., Denancé, N., Vanholme, R., Freydier, A., Martinez, Y., Hoffmann, L., et al. (2010). Walls are thin 1 (WAT1), an Arabidopsis homolog of Medicago truncatula NODULIN21, is a tonoplast-localized protein required for secondary wall formation in fibers. Plant J. 63, 469–483. doi: 10.1111/j.1365-313X.2010.04256.x
Reinbothe, S., Reinbothe, C., Lehmann, J., Becker, W., Apel, K., and Parthier, B. (1994). JIP60, a methyl jasmonate-induced ribosome-inactivating protein involved in plant stress reactions. Proc. Natl. Acad. Sci. U. S. A. 91, 7012–7016. doi: 10.1073/pnas.91.15.7012
Rice, A., Glick, L., Abadi, S., Einhorn, M., Kopelman, N. M., Salman-Minkov, A., et al. (2015). The Chromosome Counts Database (CCDB) - a community resource of plant chromosome numbers. New Phytol. 206, 19–26. doi: 10.1111/nph.13191
Rieping, M., and Schöffl, F. (1992). Synergistic effect of upstream sequences, CCAAT box elements, and HSE sequences for enhanced expression of chimaeric heat shock genes in transgenic tobacco. Mol. Gen. Genet. 231, 226–232. doi: 10.1007/BF00279795
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Rouster, J., Leah, R., Mundy, J., and Cameron-Mills, V. (1997). Identification of a methyl jasmonate-responsive region in the promoter of a lipoxygenase 1 gene expressed in barley grain. Plant J. 11, 513–523. doi: 10.1046/j.1365-313X.1997.11030513.x
Satoh, R., Fujita, Y., Nakashima, K., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2004). A novel subgroup of bZIP proteins functions as transcriptional activators in hypoosmolarity-responsive expression of the ProDH gene in Arabidopsis. Plant Cell Physiol. 45, 309–317. doi: 10.1093/pcp/pch036
Schluttenhofer, C., Pattanaik, S., Patra, B., and Yuan, L. (2014). Analyses of catharanthus roseus and arabidopsis thaliana WRKY transcription factors reveal involvement in jasmonate signaling. BMC Genomics 15:502. doi: 10.1186/1471-2164-15-502
Simpson, S. D., Nakashima, K., Narusaka, Y., Seki, M., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2003). Two different novel cis-acting elements of erd1, a clpA homologous arabidopsis gene function in induction by dehydration stress and dark-induced senescence. Plant J. 33, 259–270. doi: 10.1046/j.1365-313X.2003.01624.x
Smith, C. D., Edgar, R. C., Yandell, M. D., Smith, D. R., Celniker, S. E., Myers, E. W., et al. (2007). Improved repeat identification and masking in Dipterans. Gene 389, 1–9. doi: 10.1016/j.gene.2006.09.011
Söderman, E., Mattsson, J., and Engström, P. (1996). The Arabidopsis homeobox gene ATHB-7 is induced by water deficit and by abscisic acid. Plant J. 10, 375–381. doi: 10.1046/j.1365-313X.1996.10020375.x
Song, S. K., Choi, Y., Moon, Y. H., Kim, S. G., Choi, Y. D., and Lee, J. S. (2000). Systemic induction of a Phytolacca insularis antiviral protein gene by mechanical wounding, jasmonic acid, and abscisic acid. Plant Mol. Biol. 43, 439–450.
Sutoh, K., and Yamauchi, D. (2003). Two cis-acting elements necessary and sufficient for gibberellin-upregulated proteinase expression in rice seeds. Plant J. 34, 635–645. doi: 10.1046/j.1365-313x.2003.01753.x
Tartarini, A., Pittaluga, E., Marcozzi, G., Testone, G., Rodrigues-Pousada, R. A., Giannino, D., et al. (2010). Differential expression of saporin genes upon wounding, ABA treatment and leaf development. Physiol. Plant. 140, 141–152. doi: 10.1111/j.1399-3054.2010.01388.x
Van Hoeck, A., Horemans, N., Monsieurs, P., Cao, H. X., Vandenhove, H., and Blust, R. (2015). The first draft genome of the aquatic model plant Lemna minor opens the route for future stress physiology research and biotechnological applications. Biotechnol. Biofuels 8:188. doi: 10.1186/s13068-015-0381-1
Wang, L., Li, Q., Zhang, A., Zhou, W., Jiang, R., Yang, Z., et al. (2017). The phytol phosphorylation pathway is essential for the biosynthesis of phylloquinone, which is required for photosystem i stability in arabidopsis. Mol. Plant 10, 183–196. doi: 10.1016/j.molp.2016.12.006
Wang, P., Zoubenko, O., and Tumer, N. E. (1998). Reduced toxicity and broad spectrum resistance to viral and fungal infection in transgenic plants expressing pokeweed antiviral protein II. Plant Mol. Biol. 38, 957–964.
Wang, W., Tang, W., Ma, T., Niu, D., Jin, J. B., Wang, H., et al. (2016). A pair of light signaling factors FHY3 and FAR1 regulates plant immunity by modulating chlorophyll biosynthesis. J. Integr. Plant Biol. 58, 91–103. doi: 10.1111/jipb.12369
Waterhouse, R. M., Seppey, M., Simao, F. A., Manni, M., Ioannidis, P., Klioutchnikov, G., et al. (2018). BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 35, 543–548. doi: 10.1093/molbev/msx319
Widhalm, J. R., Ducluzeau, A. L., Buller, N. E., Elowsky, C. G., Olsen, L. J., and Basset, G. J. C. (2012). Phylloquinone (vitamin K1) biosynthesis in plants: two peroxisomal thioesterases of lactobacillales origin hydrolyze 1,4-dihydroxy-2-naphthoyl-coa. Plant J. 71, 205–215. doi: 10.1111/j.1365-313X.2012.04972.x
Wise, A. A., Liu, Z., and Binns, A. N. (2006). “Three methods for the introduction of foreign dna into agrobacterium,” in Methods in Molecular Biology Agrobacterium Protocols, ed. K. Wang (Totowa, NJ: Humana Press), 47–48.
Yamamoto, S., Nakano, T., Suzuki, K., and Shinshi, H. (2004). Elicitor-induced activation of transcription via W box-related cis-acting elements from a basic chitinase gene by WRKY transcription factors in tobacco. Biochim. Biophys. Acta - Gene Struct. Expr. 1679, 279–287. doi: 10.1016/j.bbaexp.2004.07.005
Yang, X., Feng, Y., He, Z., and Stoffella, P. J. (2005). Molecular mechanisms of heavy metal hyperaccumulation and phytoremediation. J. Trace Elements Med. Biol. 18, 339–353. doi: 10.1016/j.jtemb.2005.02.007
Yu, D., Chen, C., Chen, Z., Jiang, C.-J., Ono, K., Toki, S., et al. (2001). Evidence for an important role of WRKY DNA binding proteins in the regulation of NPR1 gene expression. Plant Cell 13, 1527–1540. doi: 10.1105/TPC.010115
Zhao, H., Tan, Z., Wen, X., and Wang, Y. (2017). An improved syringe agroinfiltration protocol to enhance transformation efficiency by combinative use of 5-azacytidine, ascorbate acid and tween-20. Plants 6:E9. doi: 10.3390/plants6010009
Zhao, L., Sun, Y., Le Cui, S. X., Chen, M., Yang, H. M., Liu, H. M., et al. (2011). Cd-induced changes in leaf proteome of the hyperaccumulator plant Phytolacca americana. Chemosphere 85, 56–66. doi: 10.1016/j.chemosphere.2011.06.029
Zou, C., Sun, K., Mackaluso, J. D., Seddon, A. E., Jin, R., Thomashow, M. F., et al. (2011). Cis-regulatory code of stress-responsive transcription in Arabidopsis thaliana. Proc. Natl. Acad. Sci. U. S. A. 108, 14992–14997. doi: 10.1073/pnas.1103202108
Zoubenko, O., Hudak, K., and Tumer, N. E. (2000). A non-toxic pokeweed antiviral protein mutant inhibits pathogen infection via a novel salicylic acid-independent pathway. Plant Mol. Biol. 44, 219–229.
Keywords: genome assembly, cis regulatory element, intron mediated enhancement, jasmonic acid, Phytolacca americana, pokeweed, pokeweed antiviral protein, ribosome inactivating protein
Citation: Neller KCM, Diaz CA, Platts AE and Hudak KA (2019) De novo Assembly of the Pokeweed Genome Provides Insight Into Pokeweed Antiviral Protein (PAP) Gene Expression. Front. Plant Sci. 10:1002. doi: 10.3389/fpls.2019.01002
Received: 11 April 2019; Accepted: 17 July 2019;
Published: 06 August 2019.
Edited by:Marcelo R. S. Briones, Federal University of São Paulo, Brazil
Reviewed by:Sudhir P. Singh, Center of Innovative and Applied Bioprocessing (CIAB), India
Guilherme Corrêa De Oliveira, Vale Technological Institute (ITV), Brazil
Copyright © 2019 Neller, Diaz, Platts and Hudak. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Katalin A. Hudak, firstname.lastname@example.org