Chestnut resistance to the blight disease: insights from transcriptome analysis
© Barakat et al; licensee BioMed Central Ltd. 2012
Received: 17 June 2011
Accepted: 19 March 2012
Published: 19 March 2012
A century ago, Chestnut Blight Disease (CBD) devastated the American chestnut. Backcross breeding has been underway to introgress resistance from Chinese chestnut into surviving American chestnut genotypes. Development of genomic resources for the family Fagaceae, has focused in this project on Castanea mollissima Blume (Chinese chestnut) and Castanea dentata (Marsh.) Borkh (American chestnut) to aid in the backcross breeding effort and in the eventual identification of blight resistance genes through genomic sequencing and map based cloning. A previous study reported partial characterization of the transcriptomes from these two species. Here, further analyses of a larger dataset and assemblies including both 454 and capillary sequences were performed and defense related genes with differential transcript abundance (GDTA) in canker versus healthy stem tissues were identified.
Over one and a half million cDNA reads were assembled into 34,800 transcript contigs from American chestnut and 48,335 transcript contigs from Chinese chestnut. Chestnut cDNA showed higher coding sequence similarity to genes in other woody plants than in herbaceous species. The number of genes tagged, the length of coding sequences, and the numbers of tagged members within gene families showed that the cDNA dataset provides a good resource for studying the American and Chinese chestnut transcriptomes. In silico analysis of transcript abundance identified hundreds of GDTA in canker versus healthy stem tissues. A significant number of additional DTA genes involved in the defense-response not reported in a previous study were identified here. These DTA genes belong to various pathways involving cell wall biosynthesis, reactive oxygen species (ROS), salicylic acid (SA), ethylene, jasmonic acid (JA), abscissic acid (ABA), and hormone signalling. DTA genes were also identified in the hypersensitive response and programmed cell death (PCD) pathways. These DTA genes are candidates for host resistance to the chestnut blight fungus, Cryphonectria parasitica.
Our data allowed the identification of many genes and gene network candidates for host resistance to the chestnut blight fungus, Cryphonectria parasitica. The similar set of GDTAs in American chestnut and Chinese chestnut suggests that the variation in sensitivity to this pathogen between these species may be the result of different timing and amplitude of the response of the two to the pathogen infection. Resources developed in this study are useful for functional genomics, comparative genomics, resistance breeding and phylogenetics in the Fagaceae.
The family Fagaceae includes more than 900 species in nine genera  found in temperate, subtropical, and tropical regions of the world . In the Western Hemisphere, the Fagaceae range from southern Canada to Colombia  where they grow as tall trees or, occasionally, as shrubs. The Fagaceae include major species such as chestnuts, oaks, and beeches . The Fagaceae have significant ecological and economic value. They are predominant species of most hardwood forests in the Northern Hemisphere. Nuts of Fagus (beeches), Castanea (chestnuts), Quercus (oaks) and most Castanopsis species are important sources of food for forest animals and produce a high quality extractable oil . The Fagaceae are important as a timber resource for construction, telephone poles, floors, furniture, cabinets, and other applications. American chestnut played an especially important economic and ecological role by providing food to various insects, birds, and mammals, and food, fiber, and wood for rural communities, prior to the introduction of the Chestnut Blight Disease (CBD) early in the 20th century .
Genomic resources for the Fagaceae also have practical applications. CBD is a stem canker disease caused by the fungus Cryphonectria parasitica. The pathogen infects stem tissues and kills the above ground portions of trees by girdling the cambium. Today, American chestnut exists primarily as an understory shrub, repeatedly sprouting from the root collar of blight-topped trees [8, 9]. Programs to breed timber-type blight resistant chestnut were initiated after appearance of the blight , but failed to produce the desired timber-type tree with acceptable levels of resistance. In 1983, a program was initiated by the American Chestnut Foundation  to introgress genes from blight-resistant Asian chestnut species into American chestnut via backcross breeding . The American Chestnut Cooperators Foundation, has made crosses between pure American chestnut trees that have limited resistance  for future selection. Introductions of hypovirulent fungal genotypes of C. parasitica , have successfully controlled the severity of blight in Europe, but have had only limited success in North America. Understanding the molecular basis of resistance to the CBD could enable the eventual restoration of American chestnut to the eastern North American forests .
The Genomic Tool Development for the Fagaceae (GTDF) project  was organized to develop genomic resources for several North American Fagaceae species to help address major forest health problems of American chestnut and other trees within this botanical family. The over 1.5 million cDNA sequences generated in this study for American and Chinese chestnut have been used to advance knowledge about the genetics of resistance to the CBD. A preliminary comparison of transcriptome sequences of healthy and blight infected American chestnut and Chinese chestnut resulted in the identification of potential candidate disease resistance genes . However, the 2009 report was based on only 13,000 and 15,000 contigs for American chestnut and Chinese chestnut, respectively, available at the time. We now report on over 83,000 contigs in total for the two species, which includes 9000 Sanger sequences from Chinese chestnut in addition to 454 sequence reads. Moreover, the assembly of transcript contigs was performed using the SeqMan™ NGen™ v1.2 software (DNAStar, Inc) which provided better assembly of the 454 reads than the version of the Newbler program available in 2008. The orthologs used for comparison of DTA genes between American and Chinese chestnut were better defined in this study through the use of a reciprocal best hits analysis. The DEGseq software package  used for assessing transcript abundance (TA) includes statistical analyses that make it a more powerful tool than the one used the previous analysis.
This paper reports on further analysis of the transcriptome of chestnut and the identification of candidate genes involved in the defense against CBD. The coverage of the transcriptome was analyzed for number of genes identified, sequence length, gene family size, and similarity to the transcriptome of other woody and herbaceous species. The transcript contigs include full length coding sequences from American and Chinese chestnut respectively which represents a valuable resource for genomic studies in the Fagaceae. The analysis of the differential transcript abundance (DTA) of chestnut genes using DEGseq  allowed the identification of hundreds of defense related candidate genes with DTA in canker versus healthy stem tissues. A small set of candidate genes for resistance to CBD were verified for their DTA using real time quantitative RT-PCR.
Sources of tissues and 454 sequencing results per tissue for the American and Chinese chestnut EST Datasets.
# 454 plates
Ave. 1 Read Length, bp
Total # Reads
Total bp Sequence
American chestnut cankers (infected stems)
American chestnut healthy stems
American chestnut healthy stems
American chestnut whole plant tissues, pooled
American chestnut whole plant tissues, pooled
5 cDNA libraries
Chinese chestnut var Nanking, cankers
Chinese chestnut var Mahogany, healthy stems
ChinesecChestnut var Nanking, healthy stems
Chinese chestnut var Nanking, whole plant tissues
Chinese chestnut var Mahogany whole plant tissues
5 cDNA libraries
Results of mass assembly of sequence reads from all libraries for each of American chestnut and Chinese chestnut into contigs1.
# 454 plates
# 454 reads
# of bp
AL2 of reads
AL. of contigs
% of reads in contigs
AL of contigs
% reads in contigs
Analysis of the American and Chinese chestnut transcriptomes
Coverage of the transcriptomes
Defense related genes in American and Chinese chestnut
List of some defense related GDTA in canker versus healthy stem in American chestnut.
# ACC Reads2
# ACHS Reads3
cinnamoyl CoA reductase.
isoform of 4-coumarate:CoA ligase (4 cl)
basic chitinase involved
manganese superoxide dismutase (MSD1)
acyl-CoA oxidase. Involved in jasmonate biosynthesis.
abc transporter family
hydrophilic protein lacking Cys residues
member of MRP subfamily
HopW1-1-Interacting protein 1 (WIN1).
carboxylesterase/methyl jasmonate esterase
hydroxy-2-methyl-2-(E)- butenyl 4-diphosphate synthase
SCF(TIR1) mediated degradation of Aux/IAA proteins
Arabidopsis dynamin-related protein 1E, DRP1E
lipid acyl hydrolase with wide substrate specificity
similar to the antifungal chitin-binding protein hevein
non-race specific resistance to bacterial and fungal pathogens
putative phosphoenolpyruvate carboxykinase
receptor-like protein kinase.
nucleoside diphosphate kinase type 1 (NDPK1) gene
26S proteasome AAA- ATPase subunit RPT2a (RPT2a)
plasma membrane proton ATPase
ATPase alpha subunit
a member of the annexin gene family
Leucine zipper transcription factor
protein phosphatase 2 C
Member of Alpha-expansin Gene Family
alcohol:NADP + oxidoreduct ase
malate dehydrogenase, cytosolic, putative.
catalytic/cation binding /hydrolase
glycosyl hydrolase family 17
R protein complex
putative apoplastic lipid transfer protein
pathogenesis-related thaumatin family protein
myb family transcription factor
List of some defense related GDTA in canker versus healthy stem tissues in Chinese chestnut.
# CCC Reads2
# CCHS Reads3
endo-1,4-beta glucanase (CEL2)
1-aminocyclopropane-1- carboxylate oxidase
MYB3R- and R2R3- type MYB- encoding genes
allene oxide cyclase
Glyceraldehyde-3-phosphate dehydrogenase C2
ABC transporter family
cinnamoyl CoA reductase
JAZ1 is a nuclear-localized protein
subunit 6b of cytochrome c oxidase
MLP-LIKE PROTEIN 423 (MLP423)
cysteine proteinase precursor- like protein
cytosolic short-chain dehydrogenase/reductase
squalene monooxygenase activity.
Regulatory Particle non-ATPase 12a (RPN12a)
NADP +- isocitrate dehydrogenase
arginine decarboxylase (ADC)
Myo-inositol-1-phosphate synthase isoform 2
cytochrome p450 family
member of Fe(II) transporter isolog family
organellar 3-ketoacyl-CoA thiolase
plasma membrane intrinsic protein (PIP2)
PR (pathogenesis-related) peptide
voltage-dependent anion channel
cytosolic GADPH (C subunit)
similar to the antifungal chitin- binding
manganese superoxide dismutase (MSD1)
A senescence-associated gene
protein with an ankyrin motif
cysteine proteinase inhibitor activity
dienelactone hydrolase family protein
sinapic acid:UDP-glucose glucosyltransferase
flavonol 7-O-glucosyltransferase (EC 126.96.36.199)
ELICITOR-ACTIVATED GENE 3-1 (ELI3-1)
aromatic alcohol:NADP + oxidoreductase
2-deoxy-D-arabino-heptulosonate 7-phosphate synthase
Transcription factor of the B-ZIP family
beta carbonic anhydrase
cytochrome p450 CYP74 gene family
protein with ABA 8'-hydroxylase activity
cytosolic heat shock protein AtHSP90.1
2-nitropropane dioxygenase family /NPD family
In this study, over 1.5 million sequencing reads from various tissues were generated and assembled into 34,800 and 48,501 contigs, for American and Chinese chestnut respectively. The low level of contamination with organelle DNA, the fraction of complete or nearly complete full length cDNA sequences, and the depth coverage of genes involved in various biological processes indicate that pyrosequencing is an excellent tool for gene discovery, EST sequencing, and transcriptome analysis in non-model tree species. Combining 454 and Sanger sequences can improve contig construction as the assembly using both data sets resulted in more reads being integrated into contigs than by assembly of pyrosequences alone, resulting in slightly fewer but longer contig sequences. Using a tool developed recently , that estimates the amount of sequencing needed to cover the transcriptome of a given species, that takes into consideration the sequencing platform and the number of contigs generated from each species, we determined that as few as two or three additional plates of 454 sequence from different developmental stages and physiological conditions should allow for 100% coverage of the transcriptomes of Chinese and American chestnut.
The sequencing effort of the GTDF increased the number of Fagaceae cDNAs available in GenBank from a few sequences to hundreds of thousands of sequences for each of the species studied. The cDNA sequences and contig sequences, which are publicly accessible at the website  hosted at Clemson University, provides tools for Fagaceae species-specific BLAST searches. The cDNA sequences have been used to select molecular markers such as Simple Sequence Repeats (SSRs), and Single Nucleotide Polymorphisms (SNPs) for genetic mapping in Chinese and American chestnut (Kubisiak et al., in preparation), for ordering the physical map of Chinese chestnut (Fang et al., in preparation), and for characterizing microRNAs (Barakat et al., unpublished data). The sequences were also used to identify candidate genes involved in resistance to CBD , some of which are being functionally characterized (Powell et al., unpublished data). The genes with DTA were also being examined for colocalization with QTLs for resistance to CBD.
BLASTX alignments to the proteomes of two model systems (Arabidopsis or Populus) showed that ~ 60% of the transcript contig sequences from the Fagaceae species studied have strong similarity to predicted proteins. The remaining contigs did not match any sequence in the Arabidopsis or Populus proteomes. A large fraction of these are short sequences that may originate from 3' or 5' untranslated regions, which tend to be highly divergent between species (Additional file 1: Figure S1). Part of these sequences may also correspond to non-coding RNAs, or potential chestnut-specific genes. BLAST searches against the proteomes of eight plant species with complete genome sequences (Vitis vinifera, Carica papaya, Medicago truncatula, Oryza sativa, Populus trichocarpa, Selaginella moellendorffii, Physcomitrella patens, and Chlamydomonas reinhardtii) showed that a large fraction of the Fagaceae contigs have better alignments with genes in woody species such as Vitis vinifera, Populus trichocarpa, and Carica papaya versus the non-woody species angiosperm species. Similar results have been reported for other tree species such as Liriodendron tulipifera and Carica papaya [20, 21]. Many long-lived woody plants, including Fagaceae species, exhibit extended phases of juvenile development before they reach flowering age, whereas most herbaceous plants reach flowering age in a single growing season. Thus, these observations could be associated with a slower evolutionary rate of genes in woody species versus herbaceous species.
In silico analysis of gene expression identified over two fold GDTAs in American than Chinese chestnut. However, most of the difference between the two DGTA sets lies in the number of induced housekeeping genes associated with the increase in resource utilization for plant defense at the infection site. The difference in the number of house-keeping genes induced in American and Chinese is maybe due to the amplitude of response of the two species to C. parasitica infection. The number of genes belonging to the category "response to biotic and abiotic stimuli" in American and Chinese chestnut from this is over 14 and 6 times higher than the ones reported previously using a partial dataset . Also, several genes identified previously were not confirmed in this analysis. However, many new genes with DTA were identified in this study. The discrepancy between these results is linked to the larger dataset, which includes both 454 and capillary sequences, the better contig assembly, and the use of a more powerful tool for DTA analysis.
Genes with DTA identified in this study belong to well known plant pathways such as phenylpropanoid metabolism, phytohormone (JA, ABA, ethylene and SA) signaling, cell wall biosynthesis, proteolysis, and others. These genes and pathways function at different times in the plant response to pathogens. The category of genes involved in phenylpropanoid metabolism act early in plant defense, serving to inhibit or to block the penetration and the progression of the plant pathogen. This category includes genes for biosynthesis of monolignol and other phenolic compounds. Previous studies [23–27] showed that lignin biosynthesis is crucial for cell wall apposition, one of the first lines of plant defense against invading fungi. Besides lignin, the biosynthesis of other polymers such as callose seems to follow infection as suggested by the increased transcript abundance of UDP-glucose:thiohydroximate S-glucosyltransferase . Other phenolic products that are involved in plant defense against pest and pathogens seem to be produced as well, as deduced by the presence of transcripts encoding genes such as flavanone 3-hydroxylase and flavonol 7-O-glucosyltransferase known to regulate flavonoid biosynthesis .
The second most important category of genes detected in response to the blight infection includes genes from phytohormone signaling pathways including JA, SA, and ethylene. These hormones trigger the activation of induced systemic resistance and systemic acquired resistance (SAR) to nectrotrophic pathogens [30, 31]. The SAR is an effective defense mechanism against a broad range of pathogens and insects. Several genes from the JA response pathway such as methyl jasmonate esterase (MES1), acyl-CoA oxidase, a phyB pathway, and ATPase transporter were identified . Genes involved in SA response such as hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate, HopW1-1-Interacting protein 1 (WIN1) were identified . The SA pathway, which is considered one of the major pathways involved in defense against nectrotrophic pathogens, regulates the expression of defense effector genes and systemic acquired resistance through the repression of the auxin signaling pathway [16, 33–36]. Another hormone that seems to play a role in the resistance of chestnut to CBD is abscissic acid (ABA). While ABA was described as a susceptibility factor, other studies [37, 38] showed that it activates plant defense by priming for callose deposition or by restricting the progression of the fungus Cochliobolus miyabeanus in the mesophyll of rice . Other signaling genes involved in SAR that induce numerous defense genes include apoplastic lipid transfer protein, and basic chitinase, etc .
The third category of genes with DTA in canker tissues includes genes involved in early response as part of the HR. Among these are transcripts encoding proteins such as ATPase transporter, kinases, carbonic anhydrase, AMMECR1, MIPS1, voltage-dependent anion channel, 2-deoxy-D-arabino-heptulosonate 7-phosphate (DAHP) synthase, and glutathione peroxidase that were reported previously to be involved in the hypersensitivity resistance (HR) and cell death in plants under pathogenic attack [18, 31, 40–43]. Reactive oxygen species (ROS) seem to be induced following C. parasitica infection as several genes involved in oxidative stress (alpha-dioxygenase, fumarase, cytosolic GADPH (C subunit), cytosolic ascorbate peroxidase APX1) had more abundant transcripts. Furthermore, several pathogenesis related (PR) genes such as elicitor-activated gene 3-1 (ELI3-1), aromatic alcohol:NADP + oxidoreductase, thaumatin, pathogenesis-related, and antifungal chitin-binding protein had differentially abundant transcripts in canker versus healthy stem tissues. PR proteins, of which some have antimicrobial functions , are mainly induced in localized pathogen attack around HR lesions. It is unknown what roles the HR and cell death play for chestnut defense against a necrotrophic pathogen such as C. parasitica. Alternatively, some of the genes involved in the HR may activate a systemic response of the plant or the pathogen may trigger HR to facilitate its colonization of the plant as reported for other pathogens . Several other genes involved in defense such as eIF(iso)4E  were also implicated.
The candidate genes identified in this study represent a valuable resource for studying the genetic basis underlying resistance to CBD and the isolation of the fungal pathogen resistance genes. Comparative mapping of the blight resistance quantitative trait loci (QTL) of chestnut with peach disease resistance QTL is revealing that the genes for several of the differentially abundant chestnut transcripts in canker versus healthy stem tissues map to disease resistance QTL regions in both species (Fang and collaborators, in preparation). This suggests that some of the genes identified in this study may play a major role in plant defense against the CBD. Because most of the defense genes and gene networks were induced in both American and Chinese chestnut canker tissues, the question that remains is what then accounts for one of these species being susceptible and the other resistant to CBD. The timing of the response to the pathogen infection, and the amplitude of the response, could result in Chinese chestnut resistance to the blight. More information from transcript profiling during the time course of infection in these two species is required to address this question.
This project has generated over 83,000 transcript contigs from American and Chinese chestnuts and identified hundreds of genes and regulation pathways which may be involved in chestnut resistance to the pathogen C. parasitica. These resources have also been used for SNP and SSR marker development, genetic mapping, physical mapping, and genome organization comparison among chestnut species and between chestnut and other closely related species. Several of the candidate genes identified in this study are in the process of being analyzed for their function using transformation in planta. The cDNA database generated in this project is also being used to map expressed genes and to annotate the proteome in the Chinese chestnut genome that is being assembled during submission of this report (John Carlson, personal communication).
Ten cDNA libraries, representing a range of tissues from American and Chinese chestnut were prepared for EST sequencing by 454 technology  (Table 1). Tissue samples were submerged in liquid nitrogen immediately upon collection and stored at - 80°C until use. To create cankers, the stems of chestnut trees were inoculated with the hypervirulent C. parasitica strain EP155 as described by Hebard and collaborators [15, 48]. For American chestnut and Chinese chestnut, canker tissue was collected at 5 and 14 days post-inoculation as previously described [15, 48]. These times correspond to early and late stages of interaction between the plant and the pathogen .
RNA preparation, cDNA library synthesis, and 454 sequencing
Total RNA was prepared as described previously , and assessed with a 2100 Bioanalyzer (Agilent Technologies). cDNA and 454 libraries were constructed as described . The 454 libraries were sequenced using the model GS20, for the American chestnut canker and Chinese chestnut canker libraries, and an FLX model 454 DNA sequencer (Roche Diagnostics), as previously described , for the remaining libraries. The sequence data was deposited into the Short Read Archive at the National Center for Biotechnology Information (study accession SRP000395).
To analyze the transcriptomes of American and Chinese chestnut, this project generated Sanger sequences for about ~9000 cDNA clones from a subtractive library enriched in genes highly expressed in canker tissues in Chinese versus American chestnut using capillary sequencing. A total of 8,101 cDNA sequences were obtained after filtering reads for quality. RNA was prepared using the method described previously  and reverse transcribed using the SMART PCR cDNA Synthesis Kit (Clontech, Mountain View, CA). Substrative libraries were constructed using Chinese chestnut as the tester and American chestnut as the driver following the manufacturer recommendation (Clontech). Sequencing of the substractive libraries was conducted at the Clemson University Genomics Institute by an automated Sanger sequencing protocol.
Transcript assembly and contig annotation
The 454 sequence reads were assembled into contigs using 454 Newbler (Roche Diagnostics) or SeqMan™ NGen™ v1.2 software (DNAStar, Inc), optimized for 454 next generation data. The new assemblies are available on the Fagaceae website . cDNA libraries were constructed using random priming which results in low poly A/T tail contamination and therefore no filtering was performed. Also, SeqMan removes low quality ends including homopolymer runs of poly(A/T) that have lower qualities in 454 sequencing. Contamination with mitochondrial and chloroplast genes was assessed by running a BLASTX search against Arabidopsis mitochondrion and chloroplast proteomes. An assembly using 454 and 9000 Sanger sequences was performed and compared the one that used 454 sequences only. Full-length contigs were identified by running a BLASTX search against the Arabidopsis thaliana proteome and comparing the lengths of the aligned portion of each contig and the putative proteins . The annotation of contigs was performed by BLASTX  against the Arabidopsis thaliana proteome (e-value = e-5) and the Gene Ontology (GO)  system . Comparison of GO annotation distribution between species was conducted using the GOstat program  set to the following parameters: GO-DB: tair; Min Sub-GO length: 3; P-Value Cutoff: 0.01; GO-Cluster Cutoff: -1; with no correction for multiple testing because the high dependence between GO terms will cause the test to be overly conservative. To determine which model species with most best hits to Fagaceae transcript contigs, BLAST alignments were conducted by querying the Fagaceae contigs against the proteomes of algal, moss and higher plant species with fully sequenced genomes (Chlamydomonas reinhardtii, Physcomitrella patens, Selaginella moellendorffii, Oryza sativa, Vitis vinifera, Populus trichocarpa, Carica papaya, and Arabidopsis thaliana) and the e-values of the best hits from each species were compared.
Identification of DTA in canker tissues
DEGseq  was used to identify gene specific differences in transcript abundance. The DEGseq package was chosen because it integrates several statistical methods, can estimate a theoretical replicate when an experimental one is not provided, and has been used routinely to identify DTA [16, 53, 54]. The number of 454 reads per contig for each gene was compared between canker and healthy stem tissues in American and Chinese chestnut separately. Similar analyses were performed for gene orthologs from both species. Orthologs were identified using a reciprocal best hit approach. DEGseq employs a random sampling model based on the read count in canker and healthy stem tissue libraries and performs a hypothesis test based on that model. Two theoretical four-fold local standard deviation lines can be drawn on the expression MA-plot to estimate the noise level of genes with different intensities and identify gene expression differences in different libraries. Genes passing the threshold are identified as exhibiting DTA. GO enrichment analyses were performed using Blast2Go software .
Validation tests of GDTA by real-time quantitative RT-PCR
Real-time quantitative RT-PCR tests were conducted to determine the extent to which the number of EST reads per gene obtained by shotgun sequencing accurately reflected transcript levels in the source tissues. RT-PCR estimates of transcript abundance were conducted on RNA from healthy and canker stem tissues from American chestnut and Chinese chestnut. RT-PCRs were performed as described previously [56, 57]. Quantitative real time PCRs (qRT-PCRs) were prepared using the SYBR Green Master Mix kit (Applied Biosystems) and run in an Applied Biosystems 7500 Fast Real-Time PCR system with default parameters. Primers were designed using Primer Express® software (Applied Biosystems). A gene encoding 18S rRNA was used as an endogenous standard to normalize template quantity. We used only one standard because we did not observe any tissue specific differences in expression of 18S rRNA gene in our study. In addition, RT-PCR analyses were performed to confirm the expression of GDTA already identified using in silico expression analysis. For each gene, three biological replicates (three different trees) and three technical replicates were performed. Statistical analyses used Statistica 6.0 software (StatSoft Poland Inc., Tulsa, OH, USA), to estimate the significance of the differences.
Chestnut blight disease
Expressed sequence tag
Differential transcript abundance
Genomic tool development for the Fagaceae
Million years ago
Programmed cell death
Quantitative trait locus
Single nucleotide polymorphism
Simple sequence repeat
We would like to thank our colleague Kerr Wall for his help with the GO annotation analysis and the Joint Genome Institute (Department of Energy) for providing us with the Cryphonectria parasitica proteome and EST sequences in advance of publication. This work was supported by grant DBI-PGPR-TRPGR 0605135 from The National Science Foundation (to RS), The Schatz Center for Tree Molecular Genetics (JEC), the World Class University Project R31-2009-000-20025-0 funded by the Ministry of Education, Science and Technology of Korea (to JEC), and The American Chestnut Foundation (FH).
- Flora of North America. [http://www.efloras.org].
- Oh S, Manos PS: Molecular phylogenetics and cupule evolution in Fagaceae as inferred from nuclear CRABS CLAW sequences. Taxon. 2008, 57: 434-451.Google Scholar
- Anagnostakis S: The effect of multiple importations of pests and pathogens on a native tree. Biological Invasions. 2001, 3: 245-254. 10.1023/A:1015205005751.View ArticleGoogle Scholar
- Hill J: Wildlife value of Castanea dentata past and present, the historical decline of the chestnut and its future use in restoration of natural areas. Proceedings of the International Chestnut Conference: 1994; Morgantown, West. Virginia. West Virginia University Press; 1994:186-193.Google Scholar
- Herendeen PS, Magallon-Puebla S, Lupia R, Crane PR, Kobylinska J: A preliminary conspectus of the Allon Flora from the Late Cretaceous (Late Santonian) of central Georgia, U.S.A. Ann Missouri Bot Garden. 1999, 86: 407-471. 10.2307/2666182.View ArticleGoogle Scholar
- Wang H, Moore MJ, Soltis PS, Bell CD, Brockington SF, Alexandre R, Davis CC, Latvis M, Manchester SR, Soltis DE: Rosid radiation and the rapid rise of angiosperm-dominated forests. Proc Natl Acad Sci USA. 2009, 106 (10): 3853-3858. 10.1073/pnas.0813376106.PubMedPubMed CentralView ArticleGoogle Scholar
- Sederoff R, Myburg A, Kirst M: Genomics, domestication, and evolution of forest trees. Cold Spring Harb Symp Quant Biol. 2009, 74: 303-317. 10.1101/sqb.2009.74.040.PubMedView ArticleGoogle Scholar
- Griffin GJ: Blight Control and Resoration of the American Chestnut. J Forestry. 2000, 99: 22-27.Google Scholar
- Andrade GM, Merkle SA: Enhancement of American chestnut somatic seedling production. Plant Cell Rep. 2005, 24 (6): 326-334. 10.1007/s00299-005-0941-0.PubMedView ArticleGoogle Scholar
- The American Chestnut Foundation. [http://www.acf.org].
- Hebard FV: The backcross breeding program of the American chestnut foundation. Restoration of American chestnut to forest lands--Proc of a conference and workshop: 2006; The North Carolina Arboretum. National Park Service;2006:Natural Res Rep NPS/NCR/CUE/NRR-2006/001Google Scholar
- Milgroom MG, Cortesi P: Biological control of chestnut blight with hypovirulence: a critical analysis. Annu Rev Phytopathol. 2004, 42: 311-338. 10.1146/annurev.phyto.42.040803.140325.PubMedView ArticleGoogle Scholar
- Wheeler N, Sederoff R: Role of genomics in the potential restoration of the American chestnut. Tree Genetics Genomes. 2008, 5: 181-187.View ArticleGoogle Scholar
- The Fagaceae Genomic Tool Project. [http://www.fagaceae.org].
- Barakat A, DiLoreto D, Zhang Y, Smith C, Baier K, Powell W, Wheeler N, Sederoff R, Carlson J: Comparison of the transcriptomes of American chestnut (Castanea dentata) and Chinese chestnut (Castanea mollissima) in response to the chestnut blight infection. BMC Plant Biology. 2009, 9 (1): 51-10.1186/1471-2229-9-51.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang L, Feng Z, Wang X, Zhang X: DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010, 26 (1): 136-138. 10.1093/bioinformatics/btp612.PubMedView ArticleGoogle Scholar
- Beissbarth T, Speed TP: GOstat: find statistically overrepresented Gene Ontologies within a group of genes. Bioinformatics. 2004, 20 (9): 1464-1465. 10.1093/bioinformatics/bth088.PubMedView ArticleGoogle Scholar
- Bovie C, Ongena M, Thonart P, Dommes J: Cloning and expression analysis of cDNAs corresponding to genes activated in cucumber showing systemic acquired resistance after BTH treatment. BMC Plant Biol. 2004, 4: 15-10.1186/1471-2229-4-15.PubMedPubMed CentralView ArticleGoogle Scholar
- Wall PK, Leebens-Mack J, Muller KF, Field D, Altman NS, de Pamphilis CW: PlantTribes: a gene and gene family resource for comparative genomics in plants. Nucleic Acids Res. 2008, D970-D976. 36 Database.
- Liang H, Ayyampalayam S, Wickett N, Barakat A, Xu Y, Landherr L, Ralph P, Jiao Y, Xu T, Schlarbaum S, et al: Generation of a large-scale genomic resource for functional and comparative genomics in Liriodendron tulipifera L. Tree Genetics & Genomes. 2011, 7 (5): 941-954. 10.1007/s11295-011-0386-2.View ArticleGoogle Scholar
- Ming R, Hou S, Feng Y, Yu Q, Dionne-Laporte A, Saw JH, Senin P, Wang W, Ly BV, Lewis KL, et al: The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature. 2008, 452 (7190): 991-996. 10.1038/nature06856.PubMedPubMed CentralView ArticleGoogle Scholar
- Barakat A, DiLoreto DS, Zhang Y, Smith C, Baier K, Powell WA, Wheeler N, Sederoff R, Carlson JE: Comparison of the transcriptomes of American chestnut (Castanea dentata) and Chinese chestnut (Castanea mollissima) in response to the chestnut blight infection. BMC Plant Biol. 2009, 9: 51-10.1186/1471-2229-9-51.PubMedPubMed CentralView ArticleGoogle Scholar
- Sibout R, Eudes A, Mouille G, Pollet B, Lapierre C, Jouanin L, Seguin A: CINNAMYL ALCOHOL DEHYDROGENASE-C and -D are the primary genes involved in lignin biosynthesis in the floral stem of Arabidopsis. Plant Cell. 2005, 17 (7): 2059-2076. 10.1105/tpc.105.030767.PubMedPubMed CentralView ArticleGoogle Scholar
- Qi X, Bakht S, Qin B, Leggett M, Hemmings A, Mellon F, Eagles J, Werck-Reichhart D, Schaller H, Lesot A, et al: A different function for a member of an ancient and highly conserved cytochrome P450 family: from essential sterols to plant defense. Proc Natl Acad Sci USA. 2006, 103 (49): 18848-18853. 10.1073/pnas.0607849103.PubMedPubMed CentralView ArticleGoogle Scholar
- Kawasaki T, Hisako K, Nakatsubo T, Hasegawa K, Wakabayashi K, Takahashi H, Umemura K, Umezawa T, Shimamoto K: Cinnamoyl-CoA reductase, a key enzyme in lignin biosynthesis, is an effector of small GTPase Rac in defense signaling in rice. PNAS. 2006, 103: 230-235. 10.1073/pnas.0509875103.PubMedPubMed CentralView ArticleGoogle Scholar
- Kawalleck P, Plesch G, Hahlbrock K, Somssich IE: Induction by fungal elicitor of S-adenosyl-L-methionine synthetase and S-adenosyl-L-homocysteine hydrolase mRNAs in cultured cells and leaves of Petroselinum crispum. Proc Natl Acad Sci USA. 1992, 89 (10): 4713-4717. 10.1073/pnas.89.10.4713.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang G, Ding X, Yuan M, Qiu D, Li X, Xu C, Wang S: Dual function of rice OsDR8 gene in disease resistance and thiamine accumulation. Plant Mol Biol. 2006, 60 (3): 437-449. 10.1007/s11103-005-4770-x.PubMedView ArticleGoogle Scholar
- Clay NK, Adio AM, Denoux C, Jander G, Ausubel FM: Glucosinolate metabolites required for an Arabidopsis innate immune response. Science. 2009, 323 (5910): 95-101. 10.1126/science.1164627.PubMedPubMed CentralView ArticleGoogle Scholar
- Nicholson RL, Hammerschmidt R: Phenolic compounds and their role in disease resistance. Annu Rev Phytopathol. 1992, 30: 369-389. 10.1146/annurev.py.30.090192.002101.View ArticleGoogle Scholar
- Feys BJ, Parker JE: Interplay of signaling pathways in plant disease resistance. Trends Genet. 2000, 16 (10): 449-455. 10.1016/S0168-9525(00)02107-7.PubMedView ArticleGoogle Scholar
- Glazebrook J: Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu Rev Phytopathol. 2005, 43: 205-227. 10.1146/annurev.phyto.43.040204.135923.PubMedView ArticleGoogle Scholar
- Vlot AC, Liu PP, Cameron RK, Park SW, Yang Y, Kumar D, Zhou F, Padukkavidana T, Gustafsson C, Pichersky E, et al: Identification of likely orthologs of tobacco salicylic acid-binding protein 2 and their role in systemic acquired resistance in Arabidopsis thaliana. Plant J. 2008, 56 (3): 445-456. 10.1111/j.1365-313X.2008.03618.x.PubMedView ArticleGoogle Scholar
- Lee MW, Jelenska J, Greenberg JT: Arabidopsis proteins important for modulating defense responses to Pseudomonas syringae that secrete HopW1-1. Plant J. 2008, 54 (3): 452-465. 10.1111/j.1365-313X.2008.03439.x.PubMedView ArticleGoogle Scholar
- Rajjou L, Belghazi M, Huguet R, Robin C, Moreau A, Job C, Job D: Proteomic investigation of the effect of salicylic acid on Arabidopsis seed germination and establishment of early defense mechanisms. Plant Physiol. 2006, 141 (3): 910-923. 10.1104/pp.106.082057.PubMedPubMed CentralView ArticleGoogle Scholar
- Kiba A, Nishihara M, Tsukatani N, Nakatsuka T, Kato Y, Yamamura S: A peroxiredoxin Q homolog from gentians is involved in both resistance against fungal disease and oxidative stress. Plant Cell Physiol. 2005, 46 (6): 1007-1015. 10.1093/pcp/pci109.PubMedView ArticleGoogle Scholar
- Gil MJ, Coego A, Mauch-Mani B, Jorda L, Vera P: The Arabidopsis csb3 mutant reveals a regulatory link between salicylic acid-mediated disease resistance and the methyl-erythritol 4-phosphate pathway. Plant J. 2005, 44 (1): 155-166. 10.1111/j.1365-313X.2005.02517.x.PubMedView ArticleGoogle Scholar
- Wiese J, Kranz T, Schubert S: Induction of pathogen resistance in barley by abiotic stress. Plant Biol (Stuttg). 2004, 6 (5): 529-536. 10.1055/s-2004-821176.View ArticleGoogle Scholar
- De Vleesschauwer D, Yang Y, Cruz CV, Hofte M: Abscisic acid-induced resistance against the brown spot pathogen Cochliobolus miyabeanus in rice involves MAP kinase-mediated repression of ethylene signaling. Plant Physiol. 2010, 152 (4): 2036-2052. 10.1104/pp.109.152702.PubMedPubMed CentralView ArticleGoogle Scholar
- Zander M, La Camera S, Lamotte O, Metraux JP, Gatz C: Arabidopsis thaliana class-II TGA transcription factors are essential activators of jasmonic acid/ethylene-induced defense responses. Plant J. 2010, 61 (2): 200-210.PubMedView ArticleGoogle Scholar
- La Camera S, Balague C, Gobel C, Geoffroy P, Legrand M, Feussner I, Roby D, Heitz T: The Arabidopsis patatin-like protein 2 (PLP2) plays an essential role in cell death execution and differentially affects biosynthesis of oxylipins and resistance to pathogens. Mol Plant Microbe Interact. 2009, 22 (4): 469-481. 10.1094/MPMI-22-4-0469.PubMedView ArticleGoogle Scholar
- Chen K, Du L, Chen Z: Sensitization of defense responses and activation of programmed cell death by a pathogen-induced receptor-like protein kinase in Arabidopsis. Plant Mol Biol. 2003, 53 (1-2): 61-74.PubMedView ArticleGoogle Scholar
- Li Y, Zou J, Li M, Bilgin DD, Vodkin LO, Hartman GL, Clough SJ: Soybean defense responses to the soybean aphid. New Phytol. 2008, 179 (1): 185-195. 10.1111/j.1469-8137.2008.02443.x.PubMedView ArticleGoogle Scholar
- Pennell RI, Lamb C: Programmed Cell Death in Plants. Plant Cell. 1997, 9 (7): 1157-1168. 10.1105/tpc.9.7.1157.PubMedPubMed CentralView ArticleGoogle Scholar
- Sels J, Mathys J, De Coninck BM, Cammue BP, De Bolle MF: Plant pathogenesis-related (PR) proteins: a focus on PR peptides. Plant Physiol Biochem. 2008, 46 (11): 941-950. 10.1016/j.plaphy.2008.06.011.PubMedView ArticleGoogle Scholar
- Govrin EM, Levine A: The hypersensitive response facilitates plant infection by the necrotrophic pathogen Botrytis cinerea. Curr Biol. 2000, 10 (13): 751-757. 10.1016/S0960-9822(00)00560-1.PubMedView ArticleGoogle Scholar
- Miyoshi H, Okade H, Muto S, Suehiro N, Nakashima H, Tomoo K, Natsuaki T: Turnip mosaic virus VPg interacts with Arabidopsis thaliana eIF(iso)4E and inhibits in vitro translation. Biochimie. 2008, 90 (10): 1427-1434. 10.1016/j.biochi.2008.03.013.PubMedView ArticleGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, et al: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437 (7057): 376-380.PubMedPubMed CentralGoogle Scholar
- Hebard F, Griffin G, Elkins J: Developmental histopathology of cankers incited by virulent and hypovirulent Endothia parasitica on susceptible and resistant chestnut trees. Phytopathology. 1984, 74: 140-149. 10.1094/Phyto-74-140.View ArticleGoogle Scholar
- Chang S, Puryear J, Cairney J: A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol Rep. 1993, 11: 113-116. 10.1007/BF02670468.View ArticleGoogle Scholar
- Wall PK, Leebens-Mack J, Chanderbali AS, Barakat A, Wolcott E, Liang H, Landherr L, Tomsho LP, Hu Y, Carlson JE, et al: Comparison of next generation sequencing technologies for transcriptome characterization. BMC Genomics. 2009, 10: 347-10.1186/1471-2164-10-347.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- Consortium GO: The Gene Ontology project in 2008. Nucleic Acids Res. 2008, 36: D440-D444.View ArticleGoogle Scholar
- Alsford S, Turner DJ, Obado SO, Sanchez-Flores A, Glover L, Berriman M, Hertz-Fowler C, Horn D: High-throughput phenotyping using parallel sequencing of RNA interference targets in the African trypanosome. Genome Res. 2011, 21 (6): 915-924. 10.1101/gr.115089.110.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen ZF, Matsumura K, Wang H, Arellano SM, Yan X, Alam I, Archer JA, Bajic VB, Qian PY: Toward an understanding of the molecular mechanisms of barnacle larval settlement: a comparative transcriptomic approach. PLoS One. 2011, 6 (7): e22913-10.1371/journal.pone.0022913.PubMedPubMed CentralView ArticleGoogle Scholar
- Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36 (10): 3420-3435. 10.1093/nar/gkn176.PubMedPubMed CentralView ArticleGoogle Scholar
- Barakat A, Bagniewska-Zadworna A, Frost CJ, Carlson JE: Phylogeny and expression profiling of CAD and CAD-like genes in hybrid Populus (P. deltoides × P. nigra): evidence from herbivore damage for subfunctionalization and functional divergence. BMC Plant Biol. 2010, 10: 100-10.1186/1471-2229-10-100.PubMedPubMed CentralView ArticleGoogle Scholar
- Barakat A, Bagniewska-Zadworna A, Choi A, Plakkat U, DiLoreto DS, Yellanki P, Carlson JE: The cinnamyl alcohol dehydrogenase gene family in Populus: phylogeny, organization, and expression. BMC Plant Biol. 2009, 9: 26-10.1186/1471-2229-9-26.PubMedPubMed CentralView ArticleGoogle Scholar