Transcriptome-based discovery of pathways and genes related to resistance against Fusariumhead blight in wheat landrace Wangshuibai
© Xiao et al.; licensee BioMed Central Ltd. 2013
Received: 19 December 2012
Accepted: 8 March 2013
Published: 21 March 2013
Fusarium head blight (FHB), caused mainly by Fusarium graminearum (Fg) Schwabe (teleomorph: Gibberellazeae Schwble), brings serious damage to wheat production. Chinese wheat landrace Wangshuibai is one of the most important resistance sources in the world. The knowledge of mechanism underlying its resistance to FHB is still limited.
To get an overview of transcriptome characteristics of Wangshuibai during infection by Fg, a high-throughput RNA sequencing based on next generation sequencing (NGS) technology (Illumina) were performed. Totally, 165,499 unigenes were generated and assigned to known protein databases including NCBI non-redundant protein database (nr) (82,721, 50.0%), Gene Ontology (GO) (38,184, 23.1%), Swiss-Prot (50,702, 30.6%), Clusters of orthologous groups (COG) (51,566, 31.2%) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) (30,657, 18.5%), as determined by Blastx search. With another NGS based platform, a digital gene expression (DGE) system, gene expression in Wangshuibai and its FHB susceptible mutant NAUH117 was profiled and compared at two infection stages by inoculation of Fg at 24 and 48 hour, with the aim of identifying genes involved in FHB resistance.
Pathogen-related proteins such as PR5, PR14 and ABC transporter and JA signaling pathway were crucial for FHB resistance, especially that mediated by Fhb1. ET pathway and ROS/NO pathway were not activated in Wangshuibai and may be not pivotal in defense to FHB. Consistent with the fact that in NAUH117 there presented a chromosome fragment deletion, which led to its increased FHB susceptibility, in Wangshuibai, twenty out of eighty-nine genes showed changed expression patterns upon the infection of Fg. The up-regulation of eight of them was confirmed by qRT-PCR, revealing they may be candidate genes for Fhb1 and need further functional analysis to confirm their roles in FHB resistance.
KeywordsFusarium head blight Wangshuibai NAUH117 Resistance Transcriptome Digital gene expression
Fusarium head blight (FHB) caused mainly by Fusarium graminearum (Fg) Schwabe (teleomorph: Gibberellazeae Schwble) not only reduces grain yield and quality, but is also a major safety concern when human and animal consume Fusarium-contaminated wheat products . Development and utilization of wheat varieties with FHB resistance is one of the main breeding objects, especially in the warm and humid wheat-growing regions, and has been recognized as one of the most economical, environmentally safe, and effective strategies for disease control . Wheat cultivar Sumai 3 and wheat landrace Wangshuibai are the two most important sources for FHB resistance worldwide. Characterization of the mechanism underlying FHB resistance and the efforts for breeding new varieties with both high FHB resistance and good agronomic traits have been very difficult because resistance to FHB was a complex trait  and exhibited several types of resistance such as resistance to FHB initial infection (type I) and FHB symptom spread within a spike (type II) . Molecular mapping of quantitative trait loci (QTLs) controlling FHB resistance has been extensively reported in the past ten years. However, map-based cloning of major QTLs including Fhb1, which has been mapped to the specific genome region of chromosome arm 3BS, has not been successful . Defense-related genes were transformed to susceptible or moderate susceptible wheat varieties in the aim of obtaining transgenic plants with enhanced FHB resistance. These genes included beta-1, 3-glucanase (PR2), chitinase (PR3), wheatwins (PR4), thaumatin-like protein (PR5), α-1-purothionin and so on [6–11]. However, the obtained transgenic wheat lines only showed relatively low level of improvement of FHB resistance.
Transcriptome analysis of wheat spikes during infection by Fusarium pathogen will help us to understand the mechanism underlying FHB resistance and identify genes related to FHB resistance, and finally propose effective strategies to breed resistant variety for better control of FHB. In the recent years, the Fusarium-wheat interactions have been studied by differential gene expression. Using wheat microarrays, Li and Yen suggested jasmonate and ethylene signaling pathway may mediate FHB resistance in wheat . Using a wheat cDNA microarray consisting of 5739 expressed sequence tags (ESTs), higher expression levels for three PR genes, PR-2, -4 and -5, after Fg infection has been observed in Sumai 3 compared with two FHB susceptible near-isogenic lines . Using wheat microarrays, 14 wheat gene transcripts showed accumulation differences between the resistant and susceptible alleles at Fhb1. Combined proteomic and transcriptomic approach, the FHB resistance was found to be associated with coordinated and ordered expression of diverse defense signaling pathways that involves the signaling molecules SA, JA, ET, PA, Ca2+ and ROS as well as altered secondary metabolism . Despite of these achievements of knowledge about the genes and pathways associated with FHB resistance in the above transcriptome research, there may be limitations to secure the identification of key genes related to wheat FHB resistance.
Wheat EST sequencing efforts provided genomics information and EST data, and allowed the creation of commercial wheat microarrays for high throughput gene expression profiling of different tissues or organs at a certain condition. The most currently used Affymetrix Wheat Genome Array contains 61,127 probe sets representing 55,052 transcripts from diverse wheat tissues or about 22,000 unigenes for all 42 chromosomes in the wheat genome, and key genes involved in FHB resistance may not be included. Thus, a comprehensive description of the genes expression of wheat spikes during infection is necessary for the discovery of candidate genes related to FHB resistance. Sequencing and characterizing of a cDNA library constructed from Fg infected wheat spikes did enrich the understanding of gene expression profiling. However, due to the high cost and time consuming of Sanger sequencing, the number of transcripts represented in the global gene expression assay remained limited . Over the past several years, the next generation sequencing (NGS) technology has emerged as a cutting edge approach for high-throughput sequence determination, promptly improved the efficiency and speed of gene discovery, and dramatically reduced the time, labor, and cost [17–19]. With this technology, we can easily get a large number of sequence data which could represent the gene expression profile of the tissue at a particular condition.
Genetic background interference is another important factor for gene discovery based on differential gene expression analysis. Upon the infection, the differences of gene expression patterns between FHB resistant and susceptible lines were partially due to their genetic background differences rather than the differences for specific resistance locus, even for those studies using near-isogenic lines . Pairwise expression comparison between resistant wildtypes and their susceptible mutants enables us to minimize the genetic background interference and hence is of great value for simplifying such complicated traits as FHB resistance by studying dissected individual QTL or single gene effects. Chinese wheat landrace Wangshuibai, which is originally from the severe epidemic areas of FHB disease of the Middle and Lower Region of the Yangtze River Area of China, has been recognized worldwide as an excellent FHB-resistant source. Similar to that in Sumai3, the major QTL on 3BS (namely Fhb1) harboring FHB resistance was also present in Wangshuibai . The QTL has not been cloned and the resistance mechanism remains unclear. By fast neutron irradiation, we obtained a mutant (namely NAUH117) from Wangshuibai with increased FHB susceptibility, and the mutant had a 3BS chromosome fragment deletion including the Fhb1 region . Although it is difficult to clone the Fhb1 by map-based cloning due to the deletion of a large chromosome fragment in the mutant, by comparing gene expression patterns between Wangshuibai and NAUH117 during Fg infection based on transcriptome analysis, it is feasible to uncover the mechanisms underlying FHB resistance, especially mediated by Fhb1.
By NGS of the Illumina technology, a huge number of distinct sequences (designated as unigenes) from the equally-mixed total RNA from Wangshuibai spikes at four infection stages (0, 12, 24 and 48 hai of Fg) were obtained. The unigenes were then assigned to putative functions, classifications or pathways based on sequence similarity analysis against the sequences in the public database resources. The assembled and annotated transcriptome sequences will provide a global view of gene expression of Wangshuibai spikes during infection by Fusarium and are the basis for candidate gene discovery. Based on our transcriptome data, a digital gene expression (DGE) system was used to compare the gene expression profiles of Wangshuibai and NAUH117 at two infection stages (24 and 48 hai). The comparison allowed us to reveal the molecular mechanism underlying the different resistance phenotypes of different genotypes, and may facilitate better understanding of resistance mechanisms of wheat in response to Fusarium infection, especially that conferred by the Fhb1 locus. This is the first published study on global expression profiling of FHB-related genes in common wheat using NGS technology. Our results may aid the identification of pathways and genes associated with resistance to FHB in wheat.
Transcriptome characterization of Wangshuibai spikes during infection with Fusarium graminearumby high-throughput RNA sequencing
To obtain the global gene expression profile of Wangshuibai spikes during infection by Fg RNA samples from spikes at three infection stages, 12, 24 and 48 hours after inoculation (hai), as well as the non-inoculated spikes were prepared, equally-mixed and then sequenced using the Illumina sequencing platform. After raw reads filtering and quality check, about 54 millions of 90 bp reads were obtained with 52.94% GC percentage (Additional file 1: Table S1). Transcriptome de novo assembly was carried out with short reads assembling program Trinity, resulting in 165,499 unigenes, with 62,442 unigenes have a length of more than 500 bp (see Additional file 2: Figure S1 for details).
For annotation, all the unigenes were aligned by Blastx to the NCBI non-redundant protein database (nr) and the Swiss-Prot protein database, using a cut-off E-value of 10-5. In the two databases, 82721 and 50702 unigenes (49.98% and 30.64% of all unigenes) returned the above cut-off BLAST results, respectively. The number of all unigenes with match hits in the two databases is 83027, and most of unigenes (50406 out of 50702 unigenes, 99.4%) with match hits in Swiss-Prot protein database can also be found in nr database. In both databases, the proportion of unigenes with match hits increased with the increase of the length of the assembled sequence. For sequences longer than 2,000 bp, the match efficiencies was 97.2% and 86.7%, whereas it decreased to 76.3% and 53.2% for sequences ranging from 500 to 1,000 bp, and to 38.8% and 20.3% for sequences less than 500 bp (Additional file 2: Figure S2A). Those unigenes with no match hit in the protein database might be due to their too short sequence length. As for annotation of species distribution, 67.6% and 15.5% had top matches (first hit) with sequences from rice (Oryza sativa), 9.6% and 1.8% from maize (Zea mays), 5.7% and 1.4% from wheat (Triticum aestivum), 3.2% and 1.2% from barley (Hordeum vulgare) and 10.5% and 39.2% from others, respectively (Additional file 2: Figure S2B).
To characterize the complex biological behaviors for the transcriptome, KEGG (Kyoto Encyclopedia of Genes and Genomes) database was used analyze the pathway annotation of unigenes. In total, 30,657 sequences were assigned to 121 KEGG pathways (the top ten represented networks were listed in Additional file 1: Table S2). The pathways with most representation by the unigenes were Metabolic pathways (6965, 22.72%), Biosynthesis of secondary metabolites (3584, 11.69%) and Plant-pathogen interaction (3548, 11.57%). These annotations provided a valuable clue for investigating specific processes, especially those involved in plant-pathogen interaction.
Digital gene expression (DGE) library sequencing and DGE tag annotation
Statistics of DGE sequencing for 6 libraries from Wangshuibai (W) and NAUH117 (M) before inoculation (W0 and M0) and at 24 and 48 hour after inoculation (W24, W48 and M24, M48) of Fg
Tag mapped to gene
% of tag
Tag mapped to gene
% of tag
Tags mapped to a unigene
% of tag
Tags mapped to a unigene
% of tag
% of tag
To reveal the molecular events underlying the DGE profiles, the sequenced tags of the 6 DGE libraries were mapped to our transcriptome data containing 165,499 unigenes as mentioned above. Among these clean tags, 3.77–4.09 millions (66.41–68.52%) from the 6 libraries can be mapped to unigenes (Table 1). However, a proportion of tags which matched to more than one unigenes could be assigned to specific single unigene; they only represent the resulted all-inclusive expression levels of all unigenes with identical tag sequence. Those tags which could be mapped unambiguously to single unigenes indicated that the transcription data are reliable . In this case, 20,858 to 24,446 distinct tags representing about a quarter of total tags (tags numbers from 1,439,849 to 1,598,522 or percentage from 25.36% to 27.07%) in each of 6 DGE libraries could be mapped unambiguously to one unigene in the reference database (Table 1). When considering all the libraries, 35,948 (21.72%) unigenes were detected in at least one of the libraries.
Comparison of digital gene expression values with qRT-PCR analysis results
Gene expression profile of Wangshuibai and NAUH117 during infection by Fg
In order to identify genes associated with FHB resistance, especially those related to Fhb1, FHB-responsive genes of Wangshuibai and NAUH117 upon the Fusarium infection were compared. We found that more genes were down-regulated and fewer genes were up-regulated in Wangshuibai than in NAUH117 (Figure 4). Only a small proportion of common DEGs were identified for the two genotypes, including 283 (12.9%) up-regulated and 216 (10.1%) down-regulated genes (Figure 5C). The great gene expression difference of the pattern between the two genotypes upon the infection was mainly due to their different FHB resistance level caused by the deletion at the Fhb1 locus. Thus, those specific DEGs in Wangshuibai were assumed to be associated with FHB resistance mediated by or related to Fhb1.
The difference of gene expression pattern between Wangshuibai and NAUH117 revealed by clustering analysis
Differences of the molecular events in response to Fginfection in Wangshuibai and NAUH117
Previous studies have identified genes and pathways involved in or related to FHB resistance. In our data, we examined how these genes expressed in spikes of Wangshuibai and NAUH117 after infection of Fg toward the identification of genes or pathways critical for FHB resistance in Wangshuibai, especially that mediated by Fhb1.
Two classes of genes, the UDP-glucosyltransferases (UDPG) [23, 24] and ABC transporters were reported to detoxificate the Fg-produced trichothecenes, which act as virulence factor during the infection of Fg in wheat spikes . In our research, 15 unigenes encoding UDPGs were detected with no transcript change after Fg infection except for unigene32856, which was repressed in both genotypes. Among the 8 unigenes encoding ABC transporters, unigene28380 was induced in both genotypes, while unigene9838 was induced in Wangshuibai but remained unchanged in NAUH117 (Figure 7B).
Jasmonic acid-, ethylene- and salicylic acid- pathways
Plant defense in response to microbial attack is regulated through a complex network of signaling pathways that involve three signaling molecules: salicylic acid (SA), jasmonic acid (JA) and ethylene (ET) . Conflict results have been reported for the effects of these pathways, either contributed or compromised the FHB resistance in wheat [12, 15, 28].
Twenty-five genes associated with the JA pathway were identified, including 14 genes responsible for JA biosynthesis [lipoxygenase (LOX), allene oxide synthase (AOS), allene oxide cyclase (AOC) and 12-oxo-phytodienoic acid reductase (OPR3)], and 11 genes involved in JA signaling transduction [CORONATINE INSENSITIVE1 (COI1), Jasmonate ZIM domain (JAZ) and MYC2] (Figure 7C). One LOX gene (unigene159837), one AOS gene (unigene20915) exhibited increased accumulation in both genotypes. The LOX gene (unigene140988) and one OPR3 gene (unigene38264) were induced only in Wangshuibai, and another LOX gene (unigene97104) and OPR3 gene (unigene139032), and one AOC gene (unigene68036) were induced in NAUH117 while were repressed in Wangshuibai. This suggested that JA biosynthesis pathway was activated in response to the Fg infection in both genotypes. One COI1 gene (unigene29682) and all the three JAZ genes (unigene5457, unigene51573 and unigene51574) were induced in Wangshuibai while their expressions remained unchanged or decreased in NAUH117. This suggested that the JA signal transduction pathway was normal in Wangshuibai but was blocked in NAUH117.
Twenty-six genes associated with the ET pathway were identified, including 12 genes responsible for ET biosynthesis [S-adenosylmethionine decarboxylase (SAM), ACS synthesis (ACS) and ACC oxidase (ACO)] and fourteen genes responsible for ET signaling [ETHYLENE INSENSITIVE 2 (EIN2), EIN3, ER-associated receptors (ETR) and ETHYLENE RESPONSE FACTOR1 (ERF)] (Figure 7D). Two ACO genes (unigene81651 and unigene54524) showed induced transcripts in both genotypes. One SAM gene (unigene77432) was repressed in Wangshuibai and remained unchanged in NAUH117. This suggested that ET biosynthesis pathway was activated in both Wangshuibai and NAUH117 after Fg infection. One EIN3 gene (unigene20153) and one ETR gene (unigene14311) were up-regulated in NAUH117 but were repressed in Wangshuibai. ChiB gene (unigene667), the marker gene for ET pathway was also responsive to the Fg infection in NAUH117 and remained unchanged in Wangshuibai. The results showed that ET pathway was only activated in NAUH117 while not in Wangshuibai.
Ding et at (2011) reported that the SA pathway was involved in FHB resistance at early infection stage (within about 12 hai) . In our research, none of the genes in the SA-mediated defense pathway, either for biosynthesis or signaling, was found to be differentially expressed both in Wangshuibai and NAUH117 (data not shown), indicating the SA pathway was not involved in defense response to Fg infection, especially at 24 or 48 hai.
Reactive oxygen species and Nitric oxide
In plant defense response, reactive oxygen species (ROS) and Nitric oxide (NO) are thought to regulate programmed cell death (PCD) through the establishment of the hypersensitive reaction (HR) . Sixteen genes related to the ROS production/scavenging systems were involved in the defense responses of wheat to Fg infection, including NADPH oxidase (Nox), ascorbate peroxidase (APX), peroxidase (POD), glutathione peroxidase (GPX), superoxide dismutase (SOD), catalase (CAT); Three genes related to the NO production/scavenging systems were involved in the defense responses, including NO synthase (NOS), peroxiredoxin (Prx) (Figure 7E). ROS-producing gene Nox (unigene30518) and NO-producing genes NOS (unigene57606 and unigene76235) were induced in NAUH117 while were repressed or remained unchanged in Wangshuibai. Consistently, many ROS- or NO- scavenging related genes were induced in NAUH117, including two APX genes (unigene32160 and unigene76965), two POD genes (unigene1767 and unigene34625), two GPX genes (unigene19711 and unigene31707), one SOD gene (unigene7029) and one Prx (unigene53318), while all ROS- or NO- scavenging related genes were repressed or remained stable in Wangshuibai except for one GPX gene (unigene31707) and one SOD gene (unigene13914) that were induced. The results suggested that ROS and NO were accumulated in NAUH117 while not in Wangshuibai.
Calcium signaling pathway
Identification of genes that were expressed in wildtype Wangshuibai while not in NAUH117
As mentioned above, a total of 35,948 unigenes were detected to be expressed in at least one of the 6 libraries. We found that 89 genes were expressed in at least two Wangshuibai libraries, but were not expressed or showed significantly decreased expression in at least two NAUH117’s libraries (Figure 8). We deduced that these differentially expressed genes in the two genotypes were due to the fragment deletion on the distal part of 3BS in NAUH117. Using Blastn search to the subject of wheat mapped ESTs database using a cut-off E-value of 10-5, these 89 unigenes were mapped to specific chromosomes (arms). Out of the 12 unigenes mapped to wheat mapped ESTs, 9 unigenes (75%) were mapped to chromosome arm 3BS (Additional file 1: Table S4), confirming the assumption that the undetected expression of these genes in NAUH117 was ascribed to its fragment deletion.
Transcriptome of Wangshuibai spikes during Fginfection provided comprehensive knowledge for FHB resistance candidate gene discovery
Chinese Wheat landrace Wangshuibai has been recognized as an excellent FHB-resistant source worldwide. Recent studies have focused on the characterization of the mechanism of FHB resistance at molecular level . Despite of the availability of wheat ESTs database from diverse tissues, a lack of available genetic and/or genomic information specifically for spikes of Wangshuibai during infection by Fg will be a barrier to further identify FHB resistance-related genes . In this study, by Illumina sequencing, we obtained 165,499 unigenes from Wangshuibai spikes in response to Fg infection. About half of unigenes have Blastx hits in nr database, unigenes without annotation might be due to the relatively short length of their assembled sequences. However, 948 unigenes (9.1%) with sequence length longer than 1,000 bp have no match hits with any known sequence, they probably represented novel genes and will be supplemental genomic information for wheat. It should be pointed out that although a large number of potentially interesting genes were obtained from the transcriptome data, most of them were partial sequences of the corresponding genes and some unigenes were from different regions of the same gene. Due to short size or poor alignment, some sequences were even excluded from analysis, temperately. The availability of wheat genome sequence in the near future will offer an opportunity for transcriptome sequence assembly, and getting ‘real unigenes’. Anyway, our results provided an overview of gene expression profile of Fusarium-inoculated spikes of Wangshuibai and offer a valuable set of sequence data for further FHB resistant candidate gene discovery.
Advantages and defects of DGE system combined with transcriptome in gene expression profiling study
DGE profiles of wheat spikes during Fusarium infection was the immediate direct application of our transcriptome data. The whole genome-wide gene expression analysis by DGE technique was based on computational analysis on 21-bp tags derived from the 3’ends of transcripts, an approach which has been proven successful in human and other species . We used Wangshuibai transcriptome data as reference sequence to map the obtained DGE data for samples from both Wangshuibai and NAUH117. The transcriptome data from Wangshuibai is applicable for annotation of DGEs from NAUH117 because of their same genomic background except for the distal fragment deletion in 3BS of NAUH117 which leads to its FHB susceptibility . We did not use wheat ESTs or unigenes in the public database as reference sequence for two reasons. Firstly, the public ESTs have a complicated genomics background, while our transcriptome sequences are specific for Wangshuibai or its mutant and make the DGE annotation more accurate. Secondly, many genes especially that are involved in FHB resistance may be not included in the public database due to the unsaturated genomic information. Our transcriptome sequence data offer an opportunity to improve the efficiency for DGE mapping.
The combination of transcriptome with DGE system has advantage over microarray technology, especially for the species without available genome sequences. The DGE, which generates digital rather than analog gene expression measurements, avoids many of the inherent limitations of microarray analysis. We believe microarray is still useful for accurate gene expression analysis at high throughput level , however, it is relative low sensitive for detection of rare transcripts and potentially can miss many targets that is not included on the array . NGS-based transcriptome profiling is based on sequencing of relatively short reads, with extensive sequence data for new genes discovery and depth-of-coverage to detect and quantify even the rare transcripts [33, 35, 36].
This technology based on the relatively short sequenced reads (data form RNA-seq and Tag-seq) also had its drawbacks, especially for the organisms without reference genome sequence. The unigenes derived from de novo assembly of RNA-Seq data were partial sequences of the corresponding genes and gave rise to an inevitable bias for their annotation that more than two third (67.6%) of the blastx hits in nr database belonged to rice. As for Tag-seq data, the low mapping efficiency for DGE was the main problem that only a quarter of all DGE tags were mapped to the reference genes (Table 1). More thorough story could be gained when wheat genome is available.
JA pathway played important role in FHB resistance of wheat and may be regulated by Fhb1
ET pathway was not important for FHB resistance or even compromised the FHB resistance in wheat
Ethylene is a gaseous plant hormone that not only affected myriad developmental processes, but also led to fitness responses including programmed cell death . The role of ET in host resistance appears to differ, depending upon the interaction system, improving resistance towards some pathogens, but increasing susceptibility towards others . Similar to the JA pathway, the ET pathway also includes two parts, ET biosynthesis followed by ET signaling pathway . Genes for both ET biosynthesis and signal transduction were induced in NAUH117 while not in Wangshuibai when infected by Fg and it was verified by qRT-PCR analysis for the ET pathway-related genes including SAM (unigene77432), ACO (unigene81651), EIN3 (unigene20153) and Chib (unigene667) (Figure 10). It indicated that ET pathway was not important in the FHB resistance that mediated by Fhb1. The activation of ET pathways in NAUH117 may be a passive defense response to infection due to the lack of the most effective defense pathway. It was also reported that ET signaling compromised the resistance to Fg in Arabidopsis, wheat and barley through premature senescence. So, another speculation is that ET serves as an important passive regulator in the resistant lines. The inactivation of ET signaling pathways in Wangshuibai, which could decrease the number of cell death, may contribute the resistance at necrotrophic phase of Fg. However, Ding et al.  and Li and Yen  reported contrary results. They noted that ACO, the key gene in ET biosynthesis, showed different expression patterns between the resistant and susceptible genotypes, and was induced in the resistant genotype. In our research, one of the two ACO genes (unigene54524) was induced both in the resistant and susceptible genotypes. However, most of genes involved in ET pathway were responsive to the infection in susceptible rather than resistant genotype. Our results indicated ET pathway was not related to, or might even compromise the FHB resistance in wheat.
ROS/NO and calcium signaling pathways were also not important in FHB resistance
Burst of NO and ROS was the main feature in plant defense against pathogen attack, in particular through the establishment of HR, which led to PCD to restrict the invading pathogen at the infection site [43, 44]. Comparison of ROS and NO bursts in Wangshuibai and NAUH117 showed that ROS/NO production/scavenging systems were more active in NAUH117 than exhibited in Wangshuibai. We deduced ROS/NO pathway may play the similar role as ET pathway did in response to Fg. They contributed to the necrotrophic phase during Fg infection through PCD (Figure 11). Thus, the repression of NO and ROS bursts would lead to enhanced FHB resistance in wheat.
It was reported that there was relationship between ROS and calcium signaling in plant responses to stresses, that is, Ca2+ influx was important in the oxidative burst through activation of NADPH oxidase  and conversely, ROS regulated Ca2+ channels through the activity of NADPH oxidases . Similar to the findings for ROS, calcium signaling pathway was responsive to Fg infection in NAUH117 while not activated in Wangshuibai. So, this pathway which was ready to leading to plant cell death was also repressed in Wangshuibai to keep its high FHB resistance level.
Transcriptome-based discovery of candidate genes for or related to Fhb1
Transcriptome based on NGS generates absolute rather than relative gene expression measurements. Thus, our transcriptome data enabled the identification of unexpressed genes in NAUH117 compared to Wangshuibai, which was highlighted as one of the unique merits of the method. The undetected gene expressions in NAUH117 were related with the presence of the distal deletion on 3BS in NAUH117 . These genes, especially 8 genes which were validated to be up-regulated in Wangshuibai upon Fg infection by qRT-PCR may be candidate genes for Fhb1 or related to resistance mediated by Fhb1. A notable fact is that the annotation of 5 out of the 6 genes were defense-related, including the receptor-like kinase referred to R gene, the lipid transfer protein belonged to PR14 and PR2 respectively, the pathogen-related protein gene family, the CBL-interacting protein kinase was involved in calcium signaling, and a glycine-rich protein was the physical barrier to resist the pathogen infection. The NADP-dependent oxidoreductase was responsible for ROS burst, which may lead to susceptibility to FHB, thus was not important in our research. It was also worth to investigate the functions of the two genes without any annotated known function; they may represent the novelty of the molecular mechanism of FHB resistance in wheat. The cloning and functional analysis of part of these selected genes will be our next research focus for approaching the final goal of better understanding the mechanism of FHB resistance.
Our transcriptome data provided comprehensive insight into gene expression profiles at two infection stages by Fg and facilitated the molecular mechanism study of FHB resistance in wheat. A putative network underlying resistance to FHB mediated mainly by Fhb1 was proposed (Figure 11). Wangshuibai resists the spread of Fg in the infected spikelets mainly through the activation of the JA defense pathway, which was regulated by Fhb1. The ET, ROS/NO pathway and calcium signaling pathways, which were repressed by Fhb1, either were not involved in or were passive regulators of FHB resistance. SA pathway did not contribute to FHB resistance after the time point of 24 hai of Fg. Eight genes were identified to be Fhb1 candidates and it would be of great value to further characterize their roles in FHB resistance. The cloning and functional analysis of two genes, a receptor-like kinase and a lipid transfer protein gene, was in progress in our lab.
Plant materials and fungal strain
Wangshuibai, a Chinese wheat landrace with high FHB resistance, was an important resistant source from southern China. NAUH117, a fast-neutron induced mutant from Wangshuibai, showed decreased typeII FHB resistance. Characterization of NAUH117 by genetic, molecular marker and cytogenetics analysis indicated the major FHB-resistance QTL Fhb1 was deleted, and resulted in increased FHB susceptibility . The materials were maintained by Cytogenetics Institute and grown in greenhouses at Jiangpu Experimental Station, Nanjing Agricultural University, Nanjing, Jiangsu Province, China.
The monosporic isolate of F. graminearum Fg0609 (kindly provided by Dr Xu Zhang, Jiangsu Academy of Agricultural Sciences) was used for fungal inoculations.
Preparation of biological samples
Wheat spikes from Wangshuibai and NAUH117 were point-inoculated at anthesis with freshly prepared spore suspension. Each spike of 4 central spikelets was injected with 10 μl of conidial inoculant (105 macro conidia per milliliter) and covered with a plastic bag to maintain humidity until sampling. Both non-inoculated spikes and inoculated spikes at 12, 24, 48, 72 and 96 hours after inoculation (hai) were sampled. Inoculations and sampling were conducted at 7 a.m. except for the sample at 12 hai at 7 p.m. A total of 12 samples (one treatment, two genotypes and six time points) were prepared for RNA extraction. Total RNA was isolated by using TRIzol Reagent (Invitrogen) according to the manufacturer’s protocol. RNA integrity was confirmed using the 2100 Bioanalyzer (Agilent Technologies) with a minimum RNA integrated number (RIN) value of 7.
Preparation of cDNA library and Illumina sequencing for transcriptome analysis
The samples for transcriptome analysis were the mixture of equal amount of RNA from non-inoculated spikes and spikes at 12, 24 and 48 hai of Wangshuibai. Transcriptome library with fragments between 200 to 700 bp was prepared following the Illumina’s kits provided by manufacturer and sequenced on Illumina HiSeq™ 2000 using paired-end technology in a single run.
Prior to assembly, adapters and unknown or low quality bases which negatively affect bioinformatics analysis were discarded. Transcriptome de novo assembly was carried out with short reads assembling program Trinity , resulting in distinct sequence, designated as unigenes. The initial short reads data sets were available at the NCBI Short Read Archive (SRA) with the accession number: SRX212270. The assembled sequences (200 bp and above) have been deposited in the NCBI’s TSA database (TSA BioProject: 183717). The transcriptome shotgun assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GAEF00000000. The version described in this paper was the first version, GAEF01000000.
The generated unigenes were analyzed by Blastx alignment search (e-value < 10-5) against protein databases nr, Swiss-Prot, KEGG and COG, and the best aligning results were used to determine sequence direction of unigenes. If results of different databases conflicted with each other, a priority order of nr, Swiss-Prot, KEGG and COG should be followed when deciding the sequence direction of unigenes. When a unigene happened to be unaligned to none of the above databases, software named ESTScan  was introduced to decide the sequence direction. For unigenes with sequence directions, we provided their sequences from 5’ end to 3’ end; for those without any direction we provided their sequences from assembly software.
Unigene sequences were aligned by blastx to protein databases nr, Swiss-Prot, KEGG and COG retrieved proteins with the highest sequence similarity for protein functional annotation and classification. For GO analysis of all the unigenes, we used Blast2GO  program to get GO annotation and WEGO  software for GO functional classification to understand the distribution of gene functions at the macro level.
DGE library preparation and Illumina sequencing
Each sample of 6 μg of the total RNA was purified by Oligo (dT) magnetic beads adsorption. The mRNA was then used as template to synthesize the first and second-strand cDNA by Oligo (dT) primer. The 3’ ends of tags can be generated by two types of Endonuclease: NlaIII or DpnII. Usually, the bead-bound cDNA was subsequently digested with the restriction enzyme NlaIII, which recognized and cut off the CATG sites. The fragments apart from the 3’ cDNA fragments connected to Oligo (dT) beads were washed away and the Illumina adaptor 1 is ligated to the sticky 5’ end of the digested bead-bound cDNA fragments. The junction of Illumina adaptor 1 and CATG site was the recognition site of MmeI, which is a type of Endonuclease with separated recognition sites and digestion sites. It cut at 17 bp downstream of the CATG site, and produced tags with adaptor 1. After removing the 3’ fragments with magnetic beads precipitation, Illumina adaptor 2 was ligated to the 3’ ends of tags, acquiring tags with different adaptors of both ends to form a tag library. After 15 cycles of linear PCR amplification, 105 bp fragments were purified by 6% TBE PAGE Gel electrophoresis. After denaturation, the single-chain molecules were fixed onto the Illumina Sequencing Chip (flowcell) for sequencing. Each tunnel generated millions of raw reads with sequencing length of 49 bp.
Analysis and mapping of DGE tags
Raw sequences had 3’ adaptor fragments as well as a few low-quality sequences and several types of impurities. Raw sequences were transformed into clean 21 bp (CATG + 17 bp) tags by the following steps: (1) 3’ adaptor sequence was trimmed, resulting in 21 bp tags from 49 bp, of raw sequence (2) empty reads were removed (reads with only 3’ adaptor sequences but no tags); (3) low quality tags were removed (tags with ambiguous base calls); (4) tags of unusual length were removed, leaving only tags of 21 bp; (5) non-redundant tags were removed (each tag needs to be detected at least twice to be considered reliable). All clean tags were mapped to this reference database allowing no more than 1 bp mismatch. The number of mapped clean tags were calculated for each library and then normalized to TPM (number of transcripts per million clean tags).
SYBR green real-time RT-PCR assay
All selected genes were conducted for SYBR Green Real-time RT-PCR assay. A real-time RT-PCR reaction (20 μl) included 20 ng cDNA, 0.2 μM of each prime, 1 × SYBR Premix ExTaq (TaKaRa Bio Inc., Japan). Reactions were performed on a Bio-Rad IQ single-color Real-Time PCR detection System (Bio-Rad, CA, USA) under the following conditions: 94°C for 30 s, 40 cycles of 94°C for 5 s, 60°C for 15 s and 72°C for 20s to calculate cycle threshold (Ct) values, followed by 95°C for 15 s, 60°C for 1 min, then 95°C for 15 s to obtain melt curves to ensure primer specificity. Pearson correlation coefficient was calculated between the cycle threshold (Ct) value of the qRT-PCR analysis and the log2 TPM values from the DGE analysis. The gene IDs and primer sequences were listed in Additional file 1: Table S3.
Clustering analysis of differential gene expression pattern
Genes with similar expression patterns usually meant functional correlation. We performed cluster analysis of gene expression patterns with “cluster” software  and “Java Treeview” software  to construct the Hierarchical clustering (HCL) tree (using the Pearson correlation method with average linkage).
The research was supported by grants from China Transgenic Research and Commercialization Key Special Project Grant (Grant No. 2011ZX08002-001), China postdoctoral Science Foundation (Grant No. 2012M511291), the National University Doctoral Fund of China (Grant No. 20100097120036), Jiangsu Natural Science Foundation (Grant No. BK2011636), Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD) and 111 Project of Ministry of Education of China (Grant No. B08025). We are grateful to Dr. Martin B. Dickman (Texas A&M University, USA) for useful comments on our research and his critical reading of the manuscript.
- Champeil A, Dore T, Fourbet JF: Fusarium head blight: epidemiological origin of the effects of cultural practices on head blight attacks and the production of mycotoxins by Fusarium in wheat grains. Plant Sci. 2004, 166 (6): 1389-1415. 10.1016/j.plantsci.2004.02.004.View Article
- Bai GH, Shaner G: Scab of wheat - prospects for control. Plant Dis. 1994, 78 (8): 760-766. 10.1094/PD-78-0760.View Article
- Snijders CHA: Fusarium head blight and mycotoxin contamination of wheat, a review. Neth J Plant Pathol. 1990, 96 (4): 187-198. 10.1007/BF01974256.View Article
- Mesterhazy A: Types and components of resistance to Fusarium head blight of wheat. Plant Breeding. 1995, 114 (5): 377-386. 10.1111/j.1439-0523.1995.tb00816.x.View Article
- Liu SX, Pumphrey MO, Gill BS, Trick HN, Zhang JX, Dolezel J, Chalhoub B, Anderson JA: Toward positional cloning of Fhb1, a major QTL for Fusarium head blight resistance in wheat. Cereal Res Commun. 2008, 36: 195-201. 10.1556/CRC.36.2008.Suppl.B.15.View Article
- Anand A, Zhou T, Trick HN, Gill BS, Bockus WW, Muthukrishnan S: Greenhouse and field testing of transgenic wheat plants stably expressing genes for thaumatin-like protein, chitinase and glucanase against Fusarium graminearum. J Exp Bot. 2003, 54 (384): 1101-1111. 10.1093/jxb/erg110.View ArticlePubMed
- Chen WP, Chen PD, Liu DJ, Kynast R, Friebe B, Velazhahan R, Muthukrishnan S, Gill BS: Development of wheat scab symptoms is delayed in transgenic wheat plants that constitutively express a rice thaumatin-like protein gene. Theor Appl Genet. 1999, 99 (5): 755-760. 10.1007/s001220051294.View Article
- Li WL, Faris JD, Muthukrishnan S, Liu DJ, Chen PD, Gill BS: Isolation and characterization of novel cDNA clones of acidic chitinases and β-1,3-glucanases from wheat spikes infected by Fusarium graminearum. Theor Appl Genet. 2001, 102 (2–3): 353-362.View Article
- Mackintosh CA, Lewis J, Radmer LE, Shin S, Heinen SJ, Smith LA, Wyckoff MN, Dill-Macky R, Evans CK, Kravchenko S: Overexpression of defense response genes in transgenic wheat enhances resistance to Fusarium head blight. Plant Cell Rep. 2007, 26 (4): 479-488. 10.1007/s00299-006-0265-8.PubMed CentralView ArticlePubMed
- Pritsch C, Muehlbauer GJ, Bushnell WR, Somers DA, Vance CP: Fungal development and induction of defense response genes during early infection of wheat spikes by Fusarium graminearum. Mol Plant Microbe In. 2000, 13 (2): 159-169. 10.1094/MPMI.2000.13.2.159.View Article
- Pritsch C, Vance CP, Bushnell WR, Somers DA, Hohn TM, Muehlbauer GJ: Systemic expression of defense response genes in wheat spikes as a response to Fusarium graminearum infection. Physiol Mol Plant P. 2001, 58 (1): 1-12. 10.1006/pmpp.2000.0308.View Article
- Li GL, Yen Y: Jasmonate and ethylene signaling pathway may mediate Fusarium head blight resistance in wheat. Crop Sci. 2008, 48 (5): 1888-1896. 10.2135/cropsci2008.02.0097.View Article
- Golkari S, Gilbert J, Ban T, Procunier JD: QTL-specific microarray gene expression analysis of wheat resistance to Fusarium head blight in Sumai-3 and two susceptible NILs. Genome. 2009, 52 (5): 409-418. 10.1139/G09-018.View ArticlePubMed
- Jia HY, Cho SH, Muehlbauer GJ: Transcriptome analysis of a wheat near-isogenic line pair carrying Fusarium head blight-resistant and -susceptible alleles. Mol Plant Microbe In. 2009, 22 (11): 1366-1378. 10.1094/MPMI-22-11-1366.View Article
- Ding LN, Xu HB, Yi HY, Yang LM, Kong ZX, Zhang LX, Xue SL, Jia HY, Ma ZQ: Resistance to hemi-biotrophic F. graminearum infection is associated with coordinated and ordered expression of diverse defense signaling pathways. PLoS One. 2011, 6 (4): e19008-10.1371/journal.pone.0019008.PubMed CentralView ArticlePubMed
- Hill-Ambroz K, Webb CA, Matthews AR, Li WL, Gill BS, Fellers JP: Expression analysis and physical mapping of a cDNA library of Fusarium head blight infected wheat spikes. Crop Sci. 2006, 46: S15-S26.View Article
- Shendure J, Ji H: Next-generation DNA sequencing. Nat Biotechnol. 2008, 26 (10): 1135-1145. 10.1038/nbt1486.View ArticlePubMed
- Fullwood MJ, Wei CL, Liu ET, Ruan Y: Next-generation DNA sequencing of paired-end tags (PET) for transcriptome and genome analyses. Genome Res. 2009, 19 (4): 521-532. 10.1101/gr.074906.107.PubMed CentralView ArticlePubMed
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29 (7): 644-U130. 10.1038/nbt.1883.PubMed CentralView ArticlePubMed
- Zhang X, Zhou MP, Ren LJ, Bai GH, Ma HX, Scholten OE, Guo PG, Lu WZ: Molecular characterization of Fusarium head blight resistance from wheat variety Wangshuibai. Euphytica. 2004, 139 (1): 59-64. 10.1007/s10681-004-2298-9.View Article
- Xiao J, Jia XP, Wang HY, Zhao RH, Fang YH, Gao RH, Wu ZZ, Cao AZ, Wang J, Xue ZK: A fast-neutron induced chromosome fragment deletion of 3BS in wheat landrace Wangshuibai increased its susceptibility to Fusarium head blight. Chromosome Res. 2011, 19 (2): 225-234. 10.1007/s10577-011-9192-9.View ArticlePubMed
- Wang XW, Luan JB, Li JM, Bao YY, Zhang CX, Liu SS: De novo characterization of a whitefly transcriptome and analysis of its gene expression during development. BMC Genomics. 2010, 11: 400-10.1186/1471-2164-11-400.PubMed CentralView ArticlePubMed
- Poppenberger B, Berthiller F, Lucyshyn D, Sieberer T, Schuhmacher R, Krska R, Kuchler K, Glossl J, Luschnig C, Adam G: Detoxification of the Fusarium mycotoxin deoxynivalenol by a UDP-glucosyltransferase from Arabidopsis thaliana. J Biol Chem. 2003, 278 (48): 47905-47914. 10.1074/jbc.M307552200.View ArticlePubMed
- Ma L, Shang Y, Cao A, Qi Z, Xing L, Chen P, Liu D, Wang X: Molecular cloning and characterization of an up-regulated UDP-glucosyltransferase gene induced by DON from Triticum aestivum L. cv. Wangshuibai. Mol Biol Rep. 2010, 37 (2): 785-795. 10.1007/s11033-009-9606-3.View Article
- Shang Y, Xiao J, Ma L, Wang H, Qi Z, Chen P, Liu D, Wang X: Characterization of a PDR type ABC transporter gene from wheat (Triticum aestivum L.). Chinese Sci Bull. 2009, 54 (18): 3249-3257. 10.1007/s11434-009-0553-0.View Article
- Rocha O, Ansari K, Doohan FM: Effects of trichothecene mycotoxins on eukaryotic cells: a review. Food Addit Contam. 2005, 22 (4): 369-378. 10.1080/02652030500058403.View ArticlePubMed
- Dong X: SA, JA, ethylene, and disease resistance in plants. Curr Opin Plant Biol. 1998, 1 (4): 316-323. 10.1016/1369-5266(88)80053-0.View ArticlePubMed
- Chen X, Steed A, Travella S, Keller B, Nicholson P: Fusarium graminearum exploits ethylene signalling to colonize dicotyledonous and monocotyledonous plants. New Phytol. 2009, 182 (4): 975-983. 10.1111/j.1469-8137.2009.02821.x.View ArticlePubMed
- Zaninotto F, La Camera S, Polverari A, Delledonne M: Cross talk between reactive nitrogen and oxygen species during the hypersensitive disease resistance response. Plant Physiol. 2006, 141 (2): 379-383. 10.1104/pp.106.078857.PubMed CentralView ArticlePubMed
- Diaz J, ten Have A, van Kan JAL: The role of ethylene and wound signaling in resistance of tomato to Botrytis cinerea. Plant Physiol. 2002, 129 (3): 1341-1351. 10.1104/pp.001453.PubMed CentralView ArticlePubMed
- Trewavas A: Le calcium, C’est la vie: calcium makes waves. Plant Physiol. 1999, 120 (1): 1-6. 10.1104/pp.120.1.1.PubMed CentralView ArticlePubMed
- Saha S, Sparks AB, Rago C, Akmaev V, Wang CJ, Vogelstein B, Kinzler KW, Velculescu VE: Using the transcriptome to annotate the genome. Nat Biotechnol. 2002, 20 (5): 508-512. 10.1038/nbt0502-508.View ArticlePubMed
- Malone JH, Oliver B: Microarrays, deep sequencing and the true measure of the transcriptome. BMC Biol. 2011, 9: 34-10.1186/1741-7007-9-34.PubMed CentralView ArticlePubMed
- Ye WW, Wang XL, Tao K, Lu YP, Dai TT, Dong SM, Dou DL, Gijzen M, Wang YC: Digital gene expression profiling of the phytophthora sojae transcriptome. Mol Plant Microbe In. 2011, 24 (12): 1530-1539. 10.1094/MPMI-05-11-0106.View Article
- Morrissy AS, Morin RD, Delaney A, Zeng T, McDonald H, Jones S, Zhao Y, Hirst M, Marra MA: Next-generation tag sequencing for cancer gene expression profiling. Genome Res. 2009, 19 (10): 1825-1835. 10.1101/gr.094482.109.PubMed CentralView ArticlePubMed
- Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, Boer JM, van Ommen GJ, den Dunnen JT, t Hoen PA: Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36 (21): e141-10.1093/nar/gkn705.PubMed CentralView ArticlePubMed
- Jansen C, von Wettstein D, Schafer W, Kogel KH, Felk A, Maier FJ: Infection patterns in barley and wheat spikes inoculated with wild-type and trichodiene synthase gene disrupted Fusarium graminearum. Proc Natl Acad Sci U S A. 2005, 102 (46): 16892-16897. 10.1073/pnas.0508467102.PubMed CentralView ArticlePubMed
- Loake G, Grant M: Salicylic acid in plant defence–the players and protagonists. Curr Opin Plant Biol. 2007, 10 (5): 466-472. 10.1016/j.pbi.2007.08.008.View ArticlePubMed
- Fonseca S, Chico JM, Solano R: The jasmonate pathway: the ligand, the receptor and the core signalling module. Curr Opin Plant Biol. 2009, 12 (5): 539-547. 10.1016/j.pbi.2009.07.013.View ArticlePubMed
- Maldonado AM, Doerner P, Dixon RA, Lamb CJ, Cameron RK: A putative lipid transfer protein involved in systemic resistance signalling in Arabidopsis. Nature. 2002, 419 (6905): 399-403. 10.1038/nature00962.View ArticlePubMed
- Guo H, Ecker JR: The ethylene signaling pathway: new insights. Curr Opin Plant Biol. 2004, 7 (1): 40-49. 10.1016/j.pbi.2003.11.011.View ArticlePubMed
- Wang KL, Li H, Ecker JR: Ethylene biosynthesis and signaling networks. Plant Cell. 2002, 14 (Suppl): S131-S151.PubMed CentralPubMed
- Wendehenne D, Durner J, Klessig DF: Nitric oxide: a new player in plant signalling and defence responses. Curr Opin Plant Biol. 2004, 7 (4): 449-455. 10.1016/j.pbi.2004.04.002.View ArticlePubMed
- Greenberg JT, Yao N: The role and regulation of programmed cell death in plant-pathogen interactions. Cell Microbiol. 2004, 6 (3): 201-211. 10.1111/j.1462-5822.2004.00361.x.View ArticlePubMed
- Du LQ, Ali GS, Simons KA, Hou JG, Yang TB, Reddy ASN, Poovaiah BW: Ca2+/calmodulin regulates salicylic-acid-mediated plant immunity. Nature. 2009, 457 (7233): 1154-U1116. 10.1038/nature07612.View ArticlePubMed
- Foreman J, Demidchik V, Bothwell JH, Mylona P, Miedema H, Torres MA, Linstead P, Costa S, Brownlee C, Jones JD: Reactive oxygen species produced by NADPH oxidase regulate plant cell growth. Nature. 2003, 422 (6930): 442-446. 10.1038/nature01485.View ArticlePubMed
- Iseli C, Jongeneel CV, Bucher P: ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol. 1999, 138-148.
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMed
- Ye J, Fang L, Zheng HK, Zhang Y, Chen J, Zhang ZJ, Wang J, Li ST, Li RQ, Bolund L: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34: W293-W297. 10.1093/nar/gkl031.PubMed CentralView ArticlePubMed
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.PubMed CentralView ArticlePubMed
- Saldanha AJ: Java Treeview–extensible visualization of microarray data. Bioinformatics. 2004, 20 (17): 3246-3248. 10.1093/bioinformatics/bth349.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.