Next Article in Journal
Next Generation Sequencing of Chromosome-Specific Libraries Sheds Light on Genome Evolution in Paleotetraploid Sterlet (Acipenser ruthenus)
Previous Article in Journal
A Microbiomic Analysis in African Americans with Colonic Lesions Reveals Streptococcus sp.VT162 as a Marker of Neoplastic Transformation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sequencing and De Novo Assembly of the Toxicodendron radicans (Poison Ivy) Transcriptome

1
Department of Botany and Plant Pathology, Oregon State University, Corvallis, OR 97330, USA
2
Department of Plant Pathology, Physiology, and Weed Science, Virginia Tech, Blacksburg, VA 24061, USA
*
Author to whom correspondence should be addressed.
Submission received: 31 July 2017 / Revised: 12 October 2017 / Accepted: 12 October 2017 / Published: 10 November 2017
(This article belongs to the Section Plant Genetics and Genomics)

Abstract

:
Contact with poison ivy plants is widely dreaded because they produce a natural product called urushiol that is responsible for allergenic contact delayed-dermatitis symptoms lasting for weeks. For this reason, the catchphrase most associated with poison ivy is “leaves of three, let it be”, which serves the purpose of both identification and an appeal for avoidance. Ironically, despite this notoriety, there is a dearth of specific knowledge about nearly all other aspects of poison ivy physiology and ecology. As a means of gaining a more molecular-oriented understanding of poison ivy physiology and ecology, Next Generation DNA sequencing technology was used to develop poison ivy root and leaf RNA-seq transcriptome resources. De novo assembled transcriptomes were analyzed to generate a core set of high quality expressed transcripts present in poison ivy tissue. The predicted protein sequences were evaluated for similarity to SwissProt homologs and InterProScan domains, as well as assigned both GO terms and KEGG annotations. Over 23,000 simple sequence repeats were identified in the transcriptome, and corresponding oligo nucleotide primer pairs were designed. A pan-transcriptome analysis of existing Anacardiaceae transcriptomes revealed conserved and unique transcripts among these species.

1. Introduction

Poison ivy (Toxicodendron radicans (L.) Kuntze) is widely known for its ability to induce allergenic contact delayed-dermatitis (poison ivy rash) symptoms that can last for weeks on between 25–40 million people per year in North America [1,2,3,4]. The poison ivy alkylphenol natural product responsible for the allergenic dermatitis is called urushiol [5]. Urushiol is not a single compound but rather a collection of penta- and hepta-dec(en)yl-catechol congeners that can vary in degrees of unsaturation ranging from zero to three double bonds [6,7,8,9]. Urushiol is present in all poison ivy organs [10] and accumulates within resin ducts/canals [11]. Despite the notoriety of poison ivy as a human contact allergen, there are surprisingly few peer reviewed studies in which poison ivy ecology [12,13,14] or non-urushiol physiology [15,16] are the primary topic of the study.
Poison ivy shows a high degree of natural phenotypic variation across its natural range. Variation in leaf margin lobing and pubescence of both leaves and fruit that generally align with different North American geographical regions led William Gillis to assign subspecies to T. radicans ssp.: barkleyi, divaricatum, eximium, negundo, pubens, radicans, and verrucosum [12,17]. However, these subspecies designations can be complicated by apparent hybrid intergrade zones. Interestingly, two additional subspecies of poison ivy, T. radicans ssp. orientale and ssp. hispidium, are geographically disjoined in east Asia [12]. Poison ivy shows three distinct growth habits: shrub, ground creeping stolons, and climbing lianas. This flexibility in growth habit allows poison ivy to exploit a variety of habitats, including habitats that have undergone natural disruption. Poison ivy incidence has been shown to greatly increase after either hurricane-induced [18,19] or range fire-induced [12,20] habitat disturbance.
Poison ivy ecology, especially urushiol chemical ecology, is currently enigmatic. The familiarity of human allergenic dermatitis and urushiol accumulation in all poison ivy organs, suggests that urushiol is a potent chemical defense in natural habitats. While this may well prove to be the case, the presumed ecological targets of urushiol chemical defense have yet to be conclusively demonstrated. Wild and domesticated mammals (bear, deer, mice, rabbit, squirrel, and goats) eat poison ivy leaves or drupes [21,22,23]. Likewise, a variety of wild birds eat poison ivy drupes [13,21,24]. However, none of these fauna display demonstrable adverse effects from consuming these poison ivy organs. A number of insect species are reported to interact with poison ivy [25,26]; but whether urushiol acts as a chemical defense against these herbivorous insects has not been established.
Most aspects of poison ivy physiology are also poorly characterized. Poison ivy is a dioecious plant, exhibiting at least three years of establishment prior to flowering [12]. A urushiol biosynthetic pathway is proposed [27,28], but none of the suggested biosynthetic steps have been biochemically or genetically validated. The geographically distributed T. radicans subspecies across North America suggest that poison ivy is readily adapting to local environments [12]. Consistent with this notion, three geographically-isolated poison ivy accessions showed significant accession-level differences in a range of biometric traits in a common garden experiment in which light and nutrient conditions were manipulated [29]. In 2006, poison ivy was shown to respond to elevated atmospheric CO2 levels by growing faster, producing more biomass, and producing a urushiol congener composition that is significantly enriched in more unsaturated pentadecyltriene-catechol congeners [15,16] that are associated with an increased severity of allergenic dermatitis in humans [30]. These findings suggest that expected patterns of climate change will result in poison ivy that is both more weedy and more noxious.
Next Generation DNA sequencing technologies have the potential to broadly advance poison ivy research. In particular, Illumina sequencing can provide both deep sampling and de novo assembly of plant transcriptomes derived from total plant RNA. Such transcriptomes can serve as both the basis of gene discovery and the quantification of gene expression levels. Here, we report the first de novo RNA-seq transcriptomes from poison ivy leaves and roots, as well as a comparative analysis of the currently sequenced Anacardiaceae family transcriptomes.

2. Materials and Methods

2.1. Axenic Cultured Seedlings and Plant Tissue

T. radicans subsp. radicans drupes from the VARoaCo-1 liana were germinated and grown under sterile culture conditions as previously described [31]. A typical VARoaCo-1 seedling voucher herbarium specimen (#109165) was deposited in the Virginia Tech Massey Herbarium. Drupes were plated at 20 seeds per 150 mm × 15 mm Petri plate containing sterile 0.5× MS Basal Salts media pH 5.7 [32] and 0.3% w/v Phytagel (Sigma-Aldrich Co., St. Louis, MO, USA). After four days, drupes in closed Petri plates were examined under a Leica Zoom 2000 illuminated stereo microscope (20× magnification) for an absence of microbial growth. Visually determined axenic seedlings were transferred to individual Magenta boxes or Phytatray II boxes (Sigma-Aldrich Co.) containing 0.5× MS Basal Salts media with 0.3% w/v Phytagel and then placed in a Percival CU-36l4 growth chamber (Percival Scientific, Perry, IA, USA) at 28 °C and 16 h light.

2.2. RNA Purification and Sequencing

Total RNA was extracted from two independent replications of either true leaves or roots harvested from 10–12 combined axenic poison ivy plants using a phenol/chloroform/SDS extraction protocol [33]. A 0.2 M KOAc treatment was included to precipitate and remove polysaccharides prior to the precipitation of total nucleic acids. RNA was selectively purified using a 2 M LiCl2 precipitation step, followed by subsequent purification using a Qiagen RNeasy Mini Kit (Qiagen, Venlo, The Netherlands) until an A260/280 of at least 1.8 was reached. RNA samples were submitted to the Virginia Bioinformatics Institute, Blacksburg, VA, for analysis using a BioAnalyzer (Agilent Genomics, Santa Clara, CA, USA), library preparation, and sequencing on an Illumina HiSeq 2000 (Illumina, San Diego, CA). Each of the four samples was selected for poly-A sequences and run on individual lanes of the HiSeq 2000.

2.3. De Novo Assembly of Poison Ivy Transcripts

The quality of the sequenced paired-end reads produced in each sample was observed using FastQC 0.10.1 [34] before and after trimming using Trimmomatic 0.3 [35] (Table 1). Paired-end reads were trimmed to remove Illumina sequencing adaptors, as well as portions of poor quality reads using a sliding window analysis (Trimmomatic settings: “LEADING:3 TRAILING:3 SLIDING-WINDOW:4:15 MINLEN:36”). The combined paired reads from all four samples were de novo assembled into an initial transcriptome using Trinity RNA Seq release 14 August 2013 [36]. Reads were then aligned to transcripts using Bowtie 2.1.0 [37] to assess overall assembly quality and analyze individual transcripts. Sequencing contaminants and overall assembly quality were quantified using the blobtools package [38] employing BLASTX search results (NCBI nr database, e-value 10) to assign transcripts to taxa. A subsequent reference de novo assembly and analysis of the three high quality samples (without the contaminated Root A sample) was performed using Trinity version 2.4.0. Transcriptome completeness was assessed using BUSCO 3.0.1 [39] with the embryophyta_odb9 dataset.
Assembled transcripts were annotated using Trinotate v 3.0.1 [40], as well as several other programs, including HMMER [41] to predict PFAM domains [42], SignalP [43], tmHMM [44], and RNAmmer [45]. TransDecoder v. 3.0.1 [40] was used to predict open reading frames (ORFs) and predicted translated protein sequences from assembled transcripts. Each transcript was used as a query in a BLASTX search (default parameters, e-value 10) against the SwissProt and NCBI nr databases; predicted protein sequences were also used as queries in BLASTP searches. All annotations were then combined using Trinotate. GO terms and KEGG classifications were assigned to transcripts utilizing Trinotate and the SwissProt BLASTX hits.
The Trinity pipeline suggested protocol and helper scripts were used to estimate transcript abundance and compare expression across samples. Paired-end reads from each of the three high quality samples were aligned to the assembled transcripts using Bowtie2 v.2.1.0 [37] with the default parameters, and converted to sorted bam files using SAMtools v.1.3.1 [46]. Assembled transcript read counts and TMM (trimmed mean of M)-normalized abundance estimates were calculated using RSEM [47] employing the Trinity helper script align_and_estimate_abundance.pl with the bowtie2 read alignments and the abundance_estimates_to_matrix.pl script. An exploratory differential expression analysis comparing the two leaf samples with a single root replicate was performed with edgeR [48] using the suggested protocol for no replicates from the edgeR manual. Comparisons of read counts were made using an exactTest with an assumed dispersion value of 0.1. GO terms were assigned to transcripts based on BLASTX hits to the SwissProt database, and enrichment analysis was performed using the Bioconductor package GOseq [49]. GO term enrichment was visualized using the R package metacoder [50]. GO term ids were converted to GO terms using WEGO [51].
Microsatellite simple sequence repeats (SSRs) were identified in assembled transcripts using the MIcroSAtellite identification tool (MISA) (http://pgrc.ipk-gatersleben.de/misa/) [52]. The minimum repeat limits were six for dinucleotides, and five for tri-, tetra-, penta-, and hexa-nucleotides. Primers for SSRs were identified using Primer3 [53] with default configuration settings.

2.4. Comparative Transcriptome Analysis of Anacardiaceae

The poison ivy transcriptome was checked for completeness by comparing gene content with the sequenced transcriptomes of closely related species. The assembled transcriptome for mango (Mangifera indica (L.)) was downloaded from NCBI (BioProject: PRJNA192676). Illumina sequencing data of transcriptomes from Rhus chinensis mill (SRA: SRX2011023) pistachio (Pistacia chinensis; SRA: SRX357356) and the Japanese lacquer tree (Toxicodendron vernicifluum; SRA: DRX055614) were downloaded from the NCBI SRA database and individually de novo assembled using Trinity. Predicted peptide sequences were generated from transcript nucleotide sequences using the Trinity package Transdecoder. Nucleotide transcript sequences from the five transcriptomes were then clustered using get_homologues-est v. 20170302 with options for MCL clustering “-M”, report all clusters “-t 0”, and report gene composition analysis “-c” [54]. A Venn diagram of shared gene clusters was generated using the R package systemPipeR [55]. GO term enrichment of core genes found in all analyzed taxa was determined using GOseq.
Sequenced Illumina paired reads and assembled transcripts were submitted to the NCBI SRA and TSA databases, respectively, under BioProject ID PRJNA393598.

3. Results and Discussion

3.1. De Novo Assembly and Quality Control

There are currently very few publically available sequences from T. radicans. As of 4 February 2017, there were 96 nucleotide sequences and 25 protein sequences in NCBI’s public databases for T. radicans. Most of these nucleotide sequences are gene fragments or intergenic regions used in phylogenetic population studies [56]. In order to comprehensively identify all expressed genes in T. radicans primary organs (leaves and roots), we sequenced the transcriptomes of two biological replicates comprised of roots or true leaves, each derived from >10 axenic poison ivy seedlings. Each of the four samples was submitted for RNA-seq sequencing on an individual lane on an Illumina HiSeq 2000 in order to maximize the sequencing depth of each sample. Three of the four samples produced ~100 million 100 bp paired reads (Table 1).
De novo assembled Illumina and/or Roche 454 (Roche, Basel, Switzerland) high-throughput RNA-seq data for non-model plants and other organisms has rapidly replaced previous methods for gene discovery. Transcriptomes for safflower [57], parasitic weed Cuscuta pentagona [58], dogwood tree [59], Chlorophytum borivilianum [60], chili pepper [61], and many other non-model organisms have recently been published. Most studies used one or more of Trinity [36], Oases [62], or SOAPdenovo-Trans [63] assemblers, with comparable results [59]. Trinity was chosen for this analysis due to its excellent reputation and built-in support for downstream analyses; including the management of annotations as well as differential expression analysis. A published Trinity pipeline [40] provided a good starting reference with the addition of a few modifications (Figure 1).
The read quality of each of the four samples was checked using FastQC before and after quality trimming with Trimmomatic. Leaves replicates A and B, as well as Roots replicate B, each had ~100 million paired reads, while Roots replicate A had relatively fewer reads (~30 million). Trimmed reads from all four datasets (Roots A–B, and Leaves A–B) were combined (289,325,296 paired reads in total, ~100 bp each) and de novo assembled into an initial assembly using Trinity RNA-seq employing an adapted form of a published pipeline [40] for this analysis (Figure 1).
An unexpectedly large number of transcripts had top BLAST hits to Metazoa (particularly various insects) and Bacteria. This was surprising because the poison ivy plants were grown under sterile culture conditions. Blobplots of each sequenced dataset revealed that the poor quality Root A sample was primarily composed of arthropod and bacterial annotated transcripts and was thus the primary source of contamination (Supplementary Figures S1–S4). This contamination was found at much lower or non-existent levels in the other three samples (Figure 2).
Illumina sequencing is very sensitive and can produce reads corresponding to mRNAs at very low expression levels. Therefore, minimizing the contamination of samples becomes increasingly important, as is the correct identification of assembled transcripts stemming from contamination. It is difficult to germinate untreated T. radicans seedlings without fungal contamination that often results in the blighting of poison ivy seedlings [31] by Colletotrichum fioriniae [64], which infects both plant and insect hosts [65,66]. For this reason, axenic poison ivy seedlings were used for RNA isolation. Thus, it was expected that fungal-annotated transcripts might be abundant contaminants in the T. radicans transcriptomes. However, that was not the case. Instead, the predominant contaminating sequences initially seemed to be derived from metazoans.
The metazoan contamination in this dataset was largely confined to one root sample, and could have originated from a variety of possible sources. This foreign RNA (or DNA) may have been introduced during either poison ivy RNA purification or during library preparation. No metazoan samples were loaded on the same flow cell used for all the poison ivy samples, so cross contamination during flow cell loading was not likely the source of this contamination. The plants were grown in sterile-culture conditions in closed magenta boxes or phytatrays without apparent insect contamination of the visibly axenic poison ivy seedlings. A possible source of insect nucleic acid contamination could have been from the tap water used to clean the reusable polypropylene Oakridge screw capped tubes that were used during the initial plant RNA extractions, although those tubes were pre-treated with RNaseZAP (Ambion, Thermofisher Scientific, Waltham, MA, USA) to inactivate any contaminating RNase activity and were then rinsed with DEPC-treated water. Contamination may also have occurred during library preparation, as multiple samples from a variety of organisms were being simultaneously prepared for sequencing at the fee-for-service core lab facility. Whatever the source of metazoan nucleic acid contamination, the absolute levels of any given metazoan transcript were typically much lower than those of plant-annotated transcripts, as evidenced by their consistently very low coverage values in the blobplots in each of the high quality samples (Supplementary Figures S1–S4). Therefore, the proportion of metazoan contaminating sequences in the poor quality Root A sample was aberrantly large, and necessitated its exclusion from a subsequent reference de novo assembly.
The trimmed reads from the three high quality samples (Leaves A–B, and Roots B) were then assembled together using Trinity to produce a high quality poison ivy reference de novo transcriptome. In total, Trinity assembled 74,937 components (a component is defined as all transcripts derived from a single de Bruijn graph) comprised of 197,641 transcript isoforms (Table 1). The average contig length was 1381 bp and the N50 was 2363 bp. When considering transcript expression, the 12,182 most highly-expressed transcripts represented 90% of the total expression data (Ex90) and had an N50 (Ex90N50) of 2077 bp. Linear regression analysis of non-normalized transcripts per million (TPM) values supported 7753 expressed components (Figure 3B). A large proportion of assembled transcripts were between 200 and 400 bases long, close to the minimum length (200 bp) of transcripts reported by Trinity (Figure 3A). Aligning paired reads to transcripts using Bowtie2 to assess overall assembly quality revealed that 94.61–95.94% of reads from each of the three high quality datasets mapped to assembled transcripts.
The presence or absence of known single-copy plant genes was used to assess the completeness of the assembled transcriptome. The BUSCO pipeline identified the presence of greater than 90% of the single-copy plant gene orthologs in the assembled poison ivy transcriptome (BUSCO code “C:90.7%(S:31.3%,D:59.4%),F:5.9%,M:3.4%,n:1440”). A large proportion of these complete genes were identified in duplicate (59.4%), suggesting that they were present and de novo assembled into multiple transcript isoforms. Just 49 (3.4%) of the 1440 tested single copy genes were missing from the de novo assembly, while 5.9% were fragmented or partially present.

3.2. Transcriptome Annotation

BLASTX searches of each transcript were performed to annotate individual transcripts based on similarity to known protein sequences in the SwissProt and NCBI nr database. Of the 162,873 Trinity transcripts with BLAST hits against NCBI nr, blobtools annotated 100,918 transcripts as likely Streptophyta in origin. The rest were annotated as undefined Eukaryota (2574), Chordata (24,148), Arthropoda (3868), and a variety of bacteria, fungi, diatoms, and other microorganisms. Likewise, 79% of transcripts had a BLASTX hit against the SwissProt database (Table 1). The most common top BLAST hit to SwissProt for transcripts was to the model plant Arabidopsis thaliana (80,443 transcripts) (Figure 4B). Other plants had relatively large numbers of BLAST hits as would be expected for a plant transcriptome, though other non-plant organisms such as humans and mice, as well as some fungi and bacteria, had substantial numbers of BLAST hits. These are likely from contaminants, or they may be transcripts with no close homolog in other sequenced plants. E-values of BLASTX hits were generally small, with 61% of top blast hits having e-values less than 1e-05 and 28.2% had e-values less than 1e-50 (Figure 4A).
Some transcripts annotated as coming from a non-plant source, or not annotated at all, may also be novel plant genes not previously sequenced in other members of Viridiplantae. Therefore, their closest sequenced homologs could be from fungi or other Eukarya. Sequencing the T. radicans genome would reveal if these transcripts are fungal or metazoan in origin or if they are novel plant genes in poison ivy.
BLASTX hits were also used to annotate transcripts with the KEGG pathway database. KEGG annotations were assigned to 119,797 transcripts, or 60.6% of assembled transcripts. Transcripts were also annotated with a number of other annotation software tools, including signal peptide prediction with SignalP, transmembrane domain prediction with tmHMM, and ribosomal RNA prediction with RNAMMER.

3.3. Gene Oontology Term Analysis

Gene Ontology (GO) terms were assigned to transcripts using the Trinotate package and SwissProt database BLASTX results. These terms represent common, broad functions and cellular components and provide a birds-eye view of the represented genes. GO terms were successfully assigned to 31,446 transcripts, or 15.9% of the assembled transcriptome. The most abundant GO terms in each of the three main GO hierarchies (cellular components, biological processes, and molecular functions) for the poison ivy transcriptome were quantified (Figure 5). The most common cellular components included cell, organelle, and organelle parts GO terms. The most abundant biological process GO terms included cellular process, metabolic process, and biological regulation, while the most abundant molecular functions included binding and catalytic activity.

3.4. Conserved Protein Domain Analysis with InterProScan

An InterProScan analysis was also performed in order to further annotate the assembled transcripts with protein functional domains. InterProScan annotated 128,448 transcripts, or 65% of total transcripts. The 10 most abundant InterProScan domains (Table 2) include common protein functions such as kinase domains, Leucine-rich repeat (LRR) domains, and NAD(P) binding domains. As InterProScan annotates specific protein functional domains rather than the protein as a whole, and this analysis provides additional annotation information for transcripts for which there are few or no characterized homologs.

3.5. Comparison of Leaf and Root Gene Expression

After removing the poor quality root sample and transcriptome re-assembly, an exploratory differential expression analysis was performed with the three high quality samples. While not statistically significant due to the single Root sample B, this analysis identified 8602 transcripts and 3363 components that are potentially differentially expressed (Benjamini-Hochberg adjusted p-value < 0.05, log2FC > 2) or differentially present between leaves and roots. As expected, GO term enrichment analysis of these genes revealed enrichment for chloroplast and photosynthesis associated GO terms in genes up regulated in leaves vs. roots. Oxidoreductase activity, DNA replication activity, and heme binding activity-associated GO terms, among others, were enriched in genes up regulated in roots vs. leaves (Supplementary Tables S1 and S2, Supplementary Figures S5–S7). In addition, there were 11,186 components/genes present (have at least one read mapping) in roots but not in leaves, and 16,957 components/genes present in leaves but not in roots. Taken together, these data indicated that sampling from both leaves and roots provided a broader representation of poison ivy expressed genes than either organ alone.

3.6. Pan-Anacardiaceae Transcriptome Comparison

Another method to gauge transcriptome assembly quality involves clustering transcripts with those from a closely-related species. Currently, the only members of Anacardiaceae with publicly available transcriptome data are the mango tree (Mangifera indica (L.)), Chinese pistache tree (Pistacia chinensis), Rhus chinensis, and the Japanese lacquer tree (Toxicodendron vernicifluum). The Japanese lacquer tree produces urushiol, whereas mango produces related alkylresorcinols [67]. Transcriptomes for these taxa were downloaded or de novo assembled using Trinity, and compared with the T. radicans transcriptome. Get_homologues-est output produced 308,419 total gene clusters (individual clusters may contain more than one transcript isoform per taxa). There were 1,482 genes shared by all five species, and 7062 shared by four or more, including 5,195 genes shared by T. radicans, R. chinensis, P. chinensis, and T. vernicifluum, but not M. indica. The largest proportion of transcripts were those found in only one member of the five taxa and not the others (Figure 6), which is likely due to differing sequencing depth, analyzed plant tissue samples, environmental conditions, and contamination. The tissues sampled for RNA extraction differed for most of these sampled species: with T. radicans representing leaf and root tissues, M. indica and T. vernicifluum representing only leaf tissue, R. chinensis representing insect-induced gall tissue, and P. chinensis representing combined bud/leaf/flower/seed tissue. The 1482 genes shared by all members of Anacardiaceae were enriched in GO terms for essential cell components and cell metabolism, and plant associated GO terms (plastid, chloroplast, Supplementary Table S3).

3.7. SSR Identification and Primer Design

Microsatellites are simple sequence repeats (SSRs) of one to several nucleotides commonly found in the genome of Eukaryotic organisms. These repeating regions can be used to characterize differences between populations and individuals. The MIcroSAtellite identification tool (MISA) was used to screen for SSRs in the assembled transcriptome in order to build a reference of poison ivy microsatellites. When looking for repeats of two to six bases, MISA identified 25,781 SSRs in 21,075 transcripts (Table 3). A majority (94%) of repeats were di- or tri-nucleotide repeats, particularly AG/CT and AT/TA (Figure 7). The most common tri-nucleotide repeat was AAG/CCT, which was the third most abundant SSR motif. Primer3 was used to generate primer pairs targeting these SSR regions (Supplementary Table S4). The identified SSRs and primer sequences provide a resource for poison ivy population genetics and future investigations into the genetic diversity of large populations of poison ivy.

4. Conclusions

This work provides an initial transcriptome framework for investigating the under-studied allergenic noxious weed poison ivy in North America. The poison ivy transcriptome was sequenced in two primary vegetative tissues and de novo assembled for the first time, representing the first example of poison ivy to be sequenced en mass. The transcriptome was extensively annotated and filtered for contaminants. Given the established capacity for poison ivy to respond to anticipated future elevated atmospheric CO2 levels by growing faster, accumulating more biomass, and becoming more allergenic [15,16], combined with the annual 25–40 million people who are sensitized to poison ivy allergenic dermatitis [1,2], there is strong justification for molecular investigations of poison ivy ecology and physiology during the progression of the Anthropocene. To this end, the poison ivy transcriptome is a valuable resource for hypotheses-driven gene discovery of (as yet only postulated) enzyme activities involved in urushiol biosynthesis [27], and increased plant vigor in response to elevated atmospheric CO2 levels [15,16]. The identification of more than 23,000 SSRs greatly expands the number of polymorphic loci beyond the current 42 microsatellite markers for poison ivy [68]. Detailed studies of poison ivy population genetics would do much to clarify the genetic basis for the current T. radicans subspecies that are currently defined by all too variable morphological characters and biogeography [12,69]. Genetic characterization of poison ivy populations is further justified by accession-level differential biometric responses, strongly suggesting that poison ivy is differentially adapting to different local environments in North America [29].
Advances in genomic techniques have spurred novel means of understanding how aggressive and/or invasive plants grow and interact with their environment. Manipulation of weeds at the molecular level may provide future alternatives to weed control augmenting herbicide-based methods of weed control, and prompt a move towards “molecular weed science” approaches.

Supplementary Materials

Supplementary materials can be found at www.mdpi.com/2073-4425/8/11/317/s1. Figure S1: Blobplot of Leaf A, Figure S2: Blobplot of Leaf B, Figure S3: Blobplot of Root A, Figure S4: Blobplot of Root B, Figure S5: Relative gene expression between roots and leaves in genes matching GO term “Biological Processes”, Figure S6: Relative gene expression between roots and leaves in genes matching GO term “Cellular Component”, Figure S7: Relative gene expression between roots and leaves in genes matching GO term “Molecular Function”, Table S1: GO terms enriched in leaves compared to root tissue, Table S2: GO terms enriched in root vs. leaf tissue, Table S3: GO terms enriched in genes core to Anacardiaceae, Table S4: SSR primer sequences.

Acknowledgments

This research was funded by a Virginia Bioinformatics Institute Core Lab and Fralin Life Science Institute Next-Gen DNA sequencing grant, the Virginia Agricultural Experiment Station, and Molecular/Translational Plant Science graduate student research and travel awards. AJW was supported by the Plant Pathology, Physiology, and Weed Science department, and the Genetics, Bioinformatics, and Computational Biology graduate program at Virginia Tech. Elise Benhase provided the axenic seedlings used in this study.

Author contributions

J.G.J., A.J.W., J.H.W., and G.K. conceived and designed the experiments; J.G.J. and A.J.W. performed the experiments; A.J.W. and G.K. analyzed the data; J.G.J. and J.H.W. contributed reagents/materials/analysis tools; A.J.W. and J.G.J. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest. The Virginia Bioinformatics Institute Core Lab (funding sponsor) performed the Illumina DNA libraries and performed the Illumina DNA sequencing. With the notable exception of the collection of the Illumina sequence data, the founding sponsors had no role in the design of the study; in the analyses, or interpretation of the data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Epstein, W.L. Plant-induced dermatitis. Ann. Emerg. Med. 1987, 16, 950–955. [Google Scholar] [CrossRef]
  2. Epstein, W.L. Occupational poison ivy and oak dermatitis. Derrmatol. Clin. 1994, 3, 511–516. [Google Scholar]
  3. Gayer, K.D.; Burnett, J.W. Toxicodendron dermatitis. Cutis 1988, 42, 99–100. [Google Scholar] [PubMed]
  4. Mitchell, J.; Rook, A. Botanical Dermatology: Plants and Plant Products Injurious to the Skin; Greengrass: Vancouver, BC, Canada, 1979; p. 787. [Google Scholar]
  5. Hill, G.A.; Mattacotti, V.; Graha, W.D. The toxic principle of the poison ivy. J. Am. Chem. Soc. 1934, 56, 2736–2738. [Google Scholar] [CrossRef]
  6. Majima, R. Uber den hauptbestandteil des japan-lacs, viii. Mitteilung: Stellung der doppelbindungen in der seitenkette des urushiols und beweisfuhrung, dab das urushiol eine mischung ist. Ber. Dtsch. Chem. Ges. 1922, 55, 172–191. (In German) [Google Scholar] [CrossRef]
  7. Majima, R.; Cho, S. Ueber einen hauptbestandteil des japanischen lackes. Ber. Dtsch. Chem. Ges. 1907, 40, 632. (In German) [Google Scholar] [CrossRef]
  8. Markiewitz, K.H.; Dawson, C.R. On the isolation of the allegenically active components of the toxic principle of poison ivy. J. Org. Chem. 1965, 30, 1610–1613. [Google Scholar] [CrossRef] [PubMed]
  9. Symes, W.F.; Dawson, C.R. Separation and structural determination of the olefinic components of poison ivy urushiol, cardanol and cardol. Nature 1953, 171, 841–842. [Google Scholar] [CrossRef] [PubMed]
  10. Craig, J.C.; Waller, C.W.; Billets, S.; Elsohly, M.A. New GLC analysis of urushiol congeners in different plant parts of poison ivy, Toxicodendron radicans. J. Pharmacol. Sci. 1978, 67, 483–485. [Google Scholar] [CrossRef]
  11. Aziz, M.; Sturtevant, D.; Winston, J.; Collakova, E.; Jelesko, J.G.; Chapman, K.D. MALDI-MS imaging of urushiols in poison ivy stem. Molecules 2017, 22, 711. [Google Scholar] [CrossRef] [PubMed]
  12. Gillis, W.T. The systematics and ecology of poison-ivy and the poison-oaks (Toxicodendron, Anacardiaceae). Rhodora 1971, 73, 161–237. [Google Scholar]
  13. Penner, R.; Moodie, G.E.; Staniforth, R.J. The dispersal of fruits and seeds of poison-ivy, Toxicodendron radicans, by ruffed grouse, Bonasa umbellus, and squirrels, Tamiasciurus hudsonicus and Sciurus carolinensis. Can. Field-Nat. 1999, 113, 616–620. [Google Scholar]
  14. Talley, S.M.; Lawton, R.O.; Setzer, W.N. Host preferences of Rhus radicans (Anacardiaceae) in a southern deciduous harwood forest. Ecology 1996, 77, 1271–1276. [Google Scholar] [CrossRef]
  15. Mohan, J.E.; Ziska, L.H.; Schlesinger, W.H.; Thomas, R.B.; Sicher, R.C.; George, K.; Clark, J.S. Biomass and toxicity responses of poison ivy (Toxicodendron radicans) to elevated atmospheric CO2. Proc. Natl. Acad. Sci. USA 2006, 103, 9086–9089. [Google Scholar] [CrossRef] [PubMed]
  16. Ziska, L.H.; Sicher, R.C.; George, K.; Mohan, J.E. Rising atmospheric carbon dioxide and potential impacts on the growth and toxicity of poison ivy (Toxicodendron radicans). Weed Sci. 2007, 55, 288–292. [Google Scholar] [CrossRef]
  17. Gillis, W.T. Poison ivy and its kin. Arnoldia 1975, 35, 93–123. [Google Scholar]
  18. Allen, B.P.; Pauley, E.F.; Sharitz, R.R. Hurricane impacts on liana populations in an old-growth southeastern bottomland forest. J. Torrey Bot. Soc. 1997, 124, 34–42. [Google Scholar] [CrossRef]
  19. Allen, B.P.; Sharitz, R.R.; Goebel, P.C. Twelve years post-hurricane liana dynamics in an old-growth southeastern floodplain forest. For. Ecol. Manag. 2005, 218, 259–269. [Google Scholar] [CrossRef]
  20. Catling, P.M.; Sinclair, A.; Cuddy, D. Plant community composition and relationships of disturbed and undisturbed alvar woodland. Can. Field-Nat. 2002, 116, 571–579. [Google Scholar]
  21. Martin, A.C.; Zim, H.S.; Nelson, A.L. American Wildlife & Plants; Dover Publications, Inc.: New York, NY, USA, 1951; p. 500. [Google Scholar]
  22. Popay, I.; Field, R. Grazing animals as weed control agents. Weed Technol. 1996, 10, 217–231. [Google Scholar]
  23. Senchina, D.S. Fungal and animal associates of Toxicodendron spp. (Anacardiaceae) in North America. Perspect. Plant Ecol. Evol. Syst. 2008, 10, 197–216. [Google Scholar] [CrossRef]
  24. Suthers, H.B.; Bickal, J.M.; Rodewald, P.G. Use of successional habitat and fruit resources by songbirds during autumn migration in central New Jersey. Wilson Bull. 2000, 112, 249–260. [Google Scholar] [CrossRef]
  25. Senchina, D.S. Beetle interactions with poison ivy and poison oak (Toxicodendron p. Mill. Sect. Toxicodendron, Anacardiaceae). Coleopt. Bull. 2005, 59, 328–334. [Google Scholar] [CrossRef]
  26. Senchina, D.S.; Summerville, K.S. Great diversity of insect floral associates may partially explain ecological success of poison ivy (Toxicodendron radicans subps negundo greene gillis, Anacardiaceae). Gt. Lakes Entomol. 2007, 40, 120–128. [Google Scholar]
  27. Dewick, P.M. Medicinal Natural Products: A Biosynthetic Approach; John Wiley & Sons: Chichester, UK, 1997; p. 466. [Google Scholar]
  28. Giessman, T.A. The biosynthesis of phenolic plant products. In The Biogenesis of Natural Compounds, 2nd ed.; Bernfeld, P., Ed.; Pergamon Press Ltd.: Oxford, UK, 1967; p. 1209. [Google Scholar]
  29. Jelesko, J.G.; Benhase, E.B.; Barney, J.N. Differential responses to light and nutrient availability by geographically isolated poison ivy accessions. Northeast. Nat. 2017, 24, 191–200. [Google Scholar] [CrossRef]
  30. Johnson, R.A.; Baer, H.; Kirkpatrick, C.H.; Dawson, C.R.; Khurana, R.G. Comparison of the contact allergenicity of the four pentadecylcatechols derived from poison ivy urushiol in human subjects. J. Allergy Clin. Immunol. 1972, 49, 27–35. [Google Scholar] [CrossRef]
  31. Benhase, E.B.; Jelesko, J.G. Germinating and culturing axenic poison ivy seedlings. HortScience 2013, 48, 1–5. [Google Scholar]
  32. Murashige, T.; Skoog, F. A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol. Plant. 1962, 15, 473–497. [Google Scholar] [CrossRef]
  33. Ausubel, F.A.; Brent, R.; Kingston, R.E.; Moore, D.D.; Seidman, J.G.; Smith, J.A.; Struhl, K. Current Protocols in Molecular Biology; John Wiley & Sons: New York, NY, USA, 2006. [Google Scholar]
  34. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. Available online: http://www.Bioinformatics.Babraham.Ac.Uk/projects/fastqc (accessed on 26 April 2010).
  35. Bolger, A.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  36. Grabherr, M.; Haas, B.; Yassour, M.; Levin, J.; Thompson, D.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  37. Langmead, B.; Salzberg, S. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed]
  38. Kumar, S.; Jones, M.; Koutsovoulos, G.; Clarke, M.; Blaxter, M. Blobology: Exploring raw genome data for contaminants, symbionts and parasites using taxon-annotated GC-coverage plots. Front. Genet. 2013, 4, 237. [Google Scholar] [CrossRef] [PubMed]
  39. Simao, F.; Waterhouse, R.; Ioannidis, P.; Kriventseva, E.; Zdobnov, E. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef] [PubMed]
  40. Haas, B.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.; Bowden, J.; Couger, M.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-Seq using the Trinity platform for reference generation and analysis. Nat. Protocols 2013, 8, 1494–1512. [Google Scholar] [CrossRef] [PubMed]
  41. Finn, R.; Clements, J.; Eddy, S. Hmmer web server: Interactive sequence similarity searching. Nucleic Acids Res. 2011, 39, W29–W37. [Google Scholar] [CrossRef] [PubMed]
  42. Finn, R.; Bateman, A.; Clements, J.; Coggill, P.; Eberhardt, R.; Eddy, S.; Heger, A.; Hetherington, K.; Holm, L.; Mistry, J.; et al. Pfam: The protein families database. Nucleic Acids Res. 2014, 42, D222–D230. [Google Scholar] [CrossRef] [PubMed]
  43. Petersen, T.; Brunak, S.; von, H.G.; Nielsen, H. Signalp 4.0: Discriminating signal peptides from transmembrane regions. Nat. Methods 2011, 8, 785–786. [Google Scholar] [CrossRef] [PubMed]
  44. Krogh, A.; Larsson, B.; von, H.G.; Sonnhammer, E. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J. Mol. Biol. 2001, 305, 567–580. [Google Scholar] [CrossRef] [PubMed]
  45. Lagesen, K.; Hallin, P.; Rodland, E.; Staerfeldt, H.; Rognes, T.; Ussery, D. Rnammer: Consistent and rapid annotation of ribosomalRNA genes. Nucleic Acids Res. 2007, 35, 3100–3108. [Google Scholar] [CrossRef] [PubMed]
  46. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; 1000 Genome Project Data Processing Subgroup. The sequence alignment/map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  47. Li, B.; Dewey, C. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed]
  48. Robinson, M.; McCarthy, D.; Smyth, G. edgeR: A bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed]
  49. Young, M.; Wakefield, M.; Smyth, G.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef] [PubMed]
  50. Foster, Z.; Sharpton, T.; Grunwald, N. Metacoder: An R package for visualization and manipulation of community taxonomic diversity data. PLoS Comp. Biol. 2017, 13, e1005404. [Google Scholar] [CrossRef] [PubMed]
  51. Ye, J.; Fang, L.; Zheng, H.K.; Zhang, Y.; Chen, J.; Zhang, Z.J.; Wang, J.; Li, S.T.; Li, R.Q.; Bolund, L.; et al. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34, W293–W297. [Google Scholar] [CrossRef] [PubMed]
  52. Thiel, T. Misa—Microsatellite identification tool. Available online: http://pgrc.ipk-gatersleben.de/misa/ (accessed on 31 July 2017).
  53. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S.G. Primer3--new capabilities and interfaces. Nucleic Acids Res. 2012, 40, e115. [Google Scholar] [CrossRef] [PubMed]
  54. Contreras-Moreira, B.; Cantalapiedra, C.; Garcia-Pereira, M.; Gordon, S.; Vogel, J.; Igartua, E.; Casas, A.; Vinuesa, P. Analysis of plant pan-genomes and transcriptomes with GET_HOMOLOGUES-EST, a clustering solution for sequences of the same species. Front. Plant Sci. 2017, 8, 184. [Google Scholar] [CrossRef] [PubMed]
  55. Backman, T.W.H.; Girke, T. SystemPipeR: NGS workflow and report generation environment. BMC Bioinform. 2016, 17. [Google Scholar] [CrossRef] [PubMed]
  56. Nie, Z.-L.; Sun, H.; Meng, Y.; Wen, J. Phylogenetic analysis of Toxicodendron (Anacardiaceae) and its biogeographic implications on the evolution of north temperate and tropical intercontinental disjunctions. J. Syst. Evol. 2009, 47, 416–430. [Google Scholar] [CrossRef]
  57. Lulin, H.; Xiao, Y.; Pei, S.; Wen, T.; Shangqin, H. The first Illumina-based de novo transcriptome sequencing and analysis of safflower flowers. PLoS ONE 2012, 7, e38653. [Google Scholar] [CrossRef] [PubMed]
  58. Ranjan, A.; Ichihashi, Y.; Farhi, M.; Zumstein, K.; Townsley, B.; David-Schwartz, R.; Sinha, N. De novo assembly and characterization of the transcriptome of the parasitic weed dodder identifies genes associated with plant parasitism. Plant Physiol. 2014, 166, 1186–1199. [Google Scholar] [CrossRef] [PubMed]
  59. Zhang, J.; Ruhlman, T.; Mower, J.; Jansen, R. Comparative analyses of two Geraniaceae transcriptomes using next-generation sequencing. BMC Plant Biol. 2013, 13, 228. [Google Scholar] [CrossRef] [PubMed]
  60. Kalra, S.; Puniya, B.; Kulshreshtha, D.; Kumar, S.; Kaur, J.; Ramachandran, S.; Singh, K. De novo transcriptome sequencing reveals important molecular networks and metabolic pathways of the plant, Chlorophytum borivilianum. PLoS ONE 2013, 8, e83336. [Google Scholar] [CrossRef] [PubMed]
  61. Liu, S.; Li, W.; Wu, Y.; Chen, C.; Lei, J. De novo transcriptome assembly in chili pepper (Capsicum frutescens) to identify genes involved in the biosynthesis of capsaicinoids. PLoS ONE 2013, 8, e48156. [Google Scholar] [CrossRef] [PubMed]
  62. Schulz, M.; Zerbino, D.; Vingron, M.; Birney, E. Oases: Robust de novoRNA-Seq assembly across the dynamic range of expression levels. Bioinformatics 2012, 28, 1086–1092. [Google Scholar] [CrossRef] [PubMed]
  63. Xie, Y.; Wu, G.; Tang, J.; Luo, R.; Patterson, J.; Liu, S.; Huang, W.; He, G.; Gu, S.; Li, S.; et al. SOAPdenovo-Trans: De novo transcriptome assembly with short RNA-Seq reads. Bioinformatics 2014, 30, 1660–1666. [Google Scholar] [CrossRef] [PubMed]
  64. Kasson, M.T.; Pollok, J.R.; Benhase, E.B.; Jelesko, J.G. First report of seedling blight of eastern poison ivy (Toxicodendron radicans) by Colletotrichum fioriniae in Virginia. Plant Dis. 2014, 98, 995–996. [Google Scholar] [CrossRef]
  65. Marcelino, J.; Giordano, R.; Gouli, S.; Gouli, V.; Parker, B.L.; Skinner, M.; TeBeest, D.; Cesnik, R. Colletotrichum acutatum var. Fioriniae (teleomorph: Glomerella acutata var. Fioriniae var. Nov.) infection of a scale insect. Mycologia 2008, 100, 353–374. [Google Scholar] [CrossRef] [PubMed]
  66. Marcelino, J.A.; Gouli, S.; Parker, B.L.; Skinner, M.; Giordano, R. Entomopathogenic activity of a variety of the fungus, Colletotrichum acutatum, recovered from the elongate hemlock scale, Fiorinia externa. J. Insect Sci. 2009, 9, 13. [Google Scholar] [CrossRef] [PubMed]
  67. Aguilar-Ortigoza, C.J.; Sosa, V.; Aguilar-Ortigoza, M. Toxic phenols in various Anacardiaceae species. Econ. Bot. 2003, 57, 354–364. [Google Scholar] [CrossRef]
  68. Hsu, T.W.; Shih, H.C.; Ku, C.C.; Chiang, T.Y.; Chiang, Y.C. Characterization of 42 microsatellite markers from poison ivy, Toxicodendron radicans (Anacardiaceae). Int. J. Mol. Sci. 2013, 14, 20414–20426. [Google Scholar] [CrossRef] [PubMed]
  69. Barkley, F.A. A monographic study of Rhus and its immediate allies in North and Central America, including the West Indies. Ann. MO. Bot. Gard. 1937, 37, 265–498. [Google Scholar] [CrossRef]
Figure 1. Process for de novo RNA-seq assembly and annotation. RNA from two replicates each of roots and leaves was sequenced with an Illumina HiSeq 2000. Trinity RNA seq was used to de novo assemble a reference transcriptome from the combined paired reads from all four T. radicans leaves and roots samples. Assembled transcripts were then annotated using a variety of programs.
Figure 1. Process for de novo RNA-seq assembly and annotation. RNA from two replicates each of roots and leaves was sequenced with an Illumina HiSeq 2000. Trinity RNA seq was used to de novo assemble a reference transcriptome from the combined paired reads from all four T. radicans leaves and roots samples. Assembled transcripts were then annotated using a variety of programs.
Genes 08 00317 g001
Figure 2. Transcript taxonomy annotation summary as determined by blobtools and BLASTX searches of the NCBI nr database. Transcript taxonomy annotation summary as determined by blobtools and BLASTX searches of the NCBI nr database for the four-sample de novo assembly. Annotated transcripts are displayed as percent of paired reads mapped to assembled transcripts of a given taxonomic grouping. (A) Leaves A sample summary; (B) Leaves B sample summary; (C) Roots A sample summary; (D) Roots B sample summary.
Figure 2. Transcript taxonomy annotation summary as determined by blobtools and BLASTX searches of the NCBI nr database. Transcript taxonomy annotation summary as determined by blobtools and BLASTX searches of the NCBI nr database for the four-sample de novo assembly. Annotated transcripts are displayed as percent of paired reads mapped to assembled transcripts of a given taxonomic grouping. (A) Leaves A sample summary; (B) Leaves B sample summary; (C) Roots A sample summary; (D) Roots B sample summary.
Genes 08 00317 g002
Figure 3. Assembly transcript length distribution and estimation of expressed transcripts. a. Assembly transcript length distribution. (A) Histogram shows the binned counts of the lengths of individual transcripts represented by 50 bp bins. Sequences with length greater than 5000 were binned together. The minimum transcript length reported by Trinity is 200 bp. (B) Estimation of expressed number of transcripts supported by the sequenced dataset. Plotted as circles are the number of transcripts with a minimum transcripts per million (TPM) value in any one sample at each threshold (plotted as negative values). A linear regression model fit to values between –100 and –10 neg. min. TPM is shown as a green line. The y-intercept (7753.52) indicates the estimated number of expressed transcripts supported by the sequenced dataset.
Figure 3. Assembly transcript length distribution and estimation of expressed transcripts. a. Assembly transcript length distribution. (A) Histogram shows the binned counts of the lengths of individual transcripts represented by 50 bp bins. Sequences with length greater than 5000 were binned together. The minimum transcript length reported by Trinity is 200 bp. (B) Estimation of expressed number of transcripts supported by the sequenced dataset. Plotted as circles are the number of transcripts with a minimum transcripts per million (TPM) value in any one sample at each threshold (plotted as negative values). A linear regression model fit to values between –100 and –10 neg. min. TPM is shown as a green line. The y-intercept (7753.52) indicates the estimated number of expressed transcripts supported by the sequenced dataset.
Genes 08 00317 g003
Figure 4. BLAST search results and species distribution. (A) E-value distribution for best BLASTX hit of transcripts against the SwissProt database. (B) Distribution of species for best BLASTX hit of transcripts against the SwissProt database.
Figure 4. BLAST search results and species distribution. (A) E-value distribution for best BLASTX hit of transcripts against the SwissProt database. (B) Distribution of species for best BLASTX hit of transcripts against the SwissProt database.
Genes 08 00317 g004
Figure 5. Distribution of top GO term categories. Distribution of most abundant GO terms for each of three major GO term hierarchies from GO terms assigned to T. radicans transcripts.
Figure 5. Distribution of top GO term categories. Distribution of most abundant GO terms for each of three major GO term hierarchies from GO terms assigned to T. radicans transcripts.
Genes 08 00317 g005
Figure 6. Anacardiaceae pan-transcriptome Venn diagram. Counts of orthologous transcripts shared by one or more members of Anacardiaceae are represented in each segment. The count of orthologous transcripts shared by all (core) is colored red. Regions are not to scale.
Figure 6. Anacardiaceae pan-transcriptome Venn diagram. Counts of orthologous transcripts shared by one or more members of Anacardiaceae are represented in each segment. The count of orthologous transcripts shared by all (core) is colored red. Regions are not to scale.
Genes 08 00317 g006
Figure 7. Most abundant SSR sequence motifs. Distribution of the most abundant SSR repeat motif sequences identified in assembled transcripts.
Figure 7. Most abundant SSR sequence motifs. Distribution of the most abundant SSR repeat motif sequences identified in assembled transcripts.
Genes 08 00317 g007
Table 1. Paired-end read quality per sample.
Table 1. Paired-end read quality per sample.
SampleLeaves ALeaves BRoots ARoots BTotal
Total # of paired reads97,494,16393,384,79730,270,702128,496,278349,645,940
Total # of paired reads after trimming87,835,19381,637,89514,779,168105,073,040289,325,296
% of paired reads surviving trimming90.09%87.42%48.82%81.77%
GC content46%46.5%40%47%
Total # of assembled contigs197,641
# of Trinity components74,937
Total assembled bases273,066,414
Mean contig length (bp)1381.63
Contig N50 (bp)2363
Contig GC%39.05
Longest contig (bp)16,102
# contigs > 1 kb93,444
Contigs with SwissProt BLASTX hits156,935
Contigs with KEGG annotation119,797
Contigs with GO annotation132,208
Table 2. Ten most abundant InterProScan domains.
Table 2. Ten most abundant InterProScan domains.
DomainDescriptionTranscripts with Domain
IPR027417P-loop containing nucleoside triphosphate hydrolase9503
IPR000719Protein kinase domain7405
IPR032675Leucine-rich repeat domain, L domain-like3946
IPR011009Protein kinase-like domain2937
IPR016040NAD(P)-binding domain2652
IPR020683Ankyrin repeat-containing domain2482
IPR011990Tetratricopeptide-like helical domain2273
IPR029058Alpha/Beta hydrolase fold2022
IPR000504RNA recognition motif domain1867
IPR012677Nucleotide-binding alpha-beta plait domain1866
Table 3. Single sequence repeats (SSR) marker results.
Table 3. Single sequence repeats (SSR) marker results.
MISA SSR SettingsCount
Total number of sequences examined197,641
Total size of examined sequences (bp)273,066,414
Total number of identified SSRs25,781
Number of SSR containing sequences21,075
Number of sequences containing more than 1 SSR3764
Number of SSRs present in compound formation1843
Unit SizeNumber of SSRs
di-nucleotide14,044
tri-nucleotide10,211
quad-nucleotide835
tetra-nucleotide316
hexa-nucleotide375

Share and Cite

MDPI and ACS Style

Weisberg, A.J.; Kim, G.; Westwood, J.H.; Jelesko, J.G. Sequencing and De Novo Assembly of the Toxicodendron radicans (Poison Ivy) Transcriptome. Genes 2017, 8, 317. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8110317

AMA Style

Weisberg AJ, Kim G, Westwood JH, Jelesko JG. Sequencing and De Novo Assembly of the Toxicodendron radicans (Poison Ivy) Transcriptome. Genes. 2017; 8(11):317. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8110317

Chicago/Turabian Style

Weisberg, Alexandra J., Gunjune Kim, James H. Westwood, and John G. Jelesko. 2017. "Sequencing and De Novo Assembly of the Toxicodendron radicans (Poison Ivy) Transcriptome" Genes 8, no. 11: 317. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8110317

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop