Comparative transcriptome profile of the leaf elongation zone of wild barley (Hordeum spontaneum) eibi1 mutant and its isogenic wild typeAbstract
The naturally occurring wild barley mutant eibi1/hvabcg31 suffers from severe water loss due to the permeable leaf cuticle. Eibi1/HvABCG31 encodes a full ATP-binding cassette (ABC) transporter, HvABCG31, playing a role in cutin deposition in the elongation zone of growing barley leaves. The eibi1 allele has pleiotropic effects on the appearance of leaves, plant stature, fertility, spike and grain size, and rate of germination. Comparative transcriptome profile of the leaf elongation zone of the eibi1 mutant as well as its isogenic wild type showed that various pathogenesis-related genes were up-regulated in the eibi1 mutant. The known cuticle-related genes that we analyzed did not show significant expression difference between the mutant and wild type. These results suggest that the pleiotropic effects may be a compensatory consequence of the activation of defense genes in the eibi1 mutation. Furthermore, we were able to find the mutation of the eibi1/hvabcg31 allele by comparing transcript sequences, which indicated that the RNA-Seq is useful not only for researches on general molecular mechanism but also for the identification of possible mutant genes.
Keywords: ABC transporter; defense genes; desiccation tolerance; plant cuticle; RNA-Seq
Introduction
The cuticle covers aerial plant organs and acts as an effective barrier against pests and pathogens (Reina-Pinto and Yephremov, 2009; Metraux et al., 2014). The cuticle mainly consists of wax and cutin. Cutin is made up of polyesters whose monomer composition is mainly C16 and C18 ω-hydroxylated fatty acid, and is typically modified by hydroxy- or epoxy groups in mid-chain positions (Yeats and Rose, 2013). Moreover, cutin contains glycerol and a small amount of phenolic compounds (Nawrath, 2006; Yeats and Rose, 2013).
A spontaneous cuticle mutant (eibi1) from wild barley (Hordeum spontaneum) genotype 23-19 has been cloned via a map-based approach (Chen et al., 2009, 2011b,). The eibi1 mutant suffers from severe water loss related to a reduced cuticle thickness and a decreased amount of cutin monomers, and shows twisted leaves, dwarf plant, low fertility, reduced spike and grain size, and germination delay. The gene, named as Eibi1/HvABCG31 is mapped to a pericentromeric region on chromosome 3H, where a mutation in an ABCG31 transporter gene is associated with the eibi1 mutant phenotype. This PDR (pleiotropic drug resistance) transporter family, as well as its homologs in Arabidopsis and rice, AtABCG32 and OsABCG31, plays an important role in cutin deposition during the development of a functional cuticle (Bessire et al., 2011; Chen et al., 2011a; Garroum et al., 2016). Analysis of Eibi1 gene expression shows abundant transcripts in the elongation zone (EZ) but only traces in non-elongation zones (NEZ) and emerged blade (EmBL) of a growing leaf, and none in the mature root (Chen et al., 2011b). Cutin but not wax deposition occurs mostly in the EZ (Richardson et al., 2007). Cutin deposition is already established by the time the EZ is formed, and the defect of the eibi1 mutant cuticle is apparent in the EZ (Chen et al., 2011b). However, gene transcription in eibi1 EZ is poorly understood.
A comparative transcriptomic analysis of the second leaves of near-isogenic eibi1 and the wild type lines has been conducted using the 22-k Barley1 Affymetrix microarray and the results revealed a pleiotropic effect of HvABCG31 gene on cuticle biogenesis and stress responsive pathways (Yang et al., 2013). In the present study, we conducted the transcriptome analysis (RNA-Seq) in the EZ of the third leaf of eibi1 mutant as well as its wild type, 23-19, with four purposes: (a) to reveal the transcriptome of the wild barley leaf EZ; (b) to test if RNA-Seq can be used for eibi1 mutation identification; (c) to test if the mutation affects the expression of other cuticle-related genes; and (d) to identify the major effects of eibi1 mutation on the other genes at the expression level.
Materials and Methods
Plant material
The wild barley accession 23-19 was selected from the Wadi Qilt population maintained at the Gene Bank of the Institute of Evolution, University of Haifa, Israel. The eibi1 mutant arose from 23-19 (Chen et al., 2004). Caryopses of the eibi1 mutant and wild type 23-19 were sown in 2.5-L pots filled with commercial compost, and the plants were grown in an incubator at 65% relative humidity, 22/16 °C (12-hour day/night cycle) under fluorescent light. The EZ of leaf three was sampled from a seedling at three-leaf stage. Six independent replicate samples were mixed for RNA extraction for both mutant and wild type. The EZ in this research is the region within 25 mm from leaf three insertion point.
RNA extraction, quality determination and validation of RNA-Seq by qRT-PCR
Total RNA of mutant or wild type was isolated using RNA kit (Tiangen, Beijing, China) following manufacturer’s instructions. The concentration and quality of RNA were verified by absorption 260/280 nm ratio between 1.8 and 2.0 using Nanodrop2000 spectrophotometer (Thermo Scientific, USA). RNA quality was checked on a Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA) and RNA integrity number (PRIN) values were greater than 8.6 for mutant and wild type.
Total RNA from 23-19 wild type and eibi1 mutant were treated with DNase and reverse-transcribed using a reverse transcription kit (RR047A) (Takara Biomedical Technology, Beijing, China). Quantitative RT-PCR was performed in an Agilent Stratagene MX3000P (Agilent) using Power SYBR green chemistry (RR820A) (Takara). All qRT-PCR reactions were performed in triplicate following the reaction conditions: 95 °C for 10 min followed by 40 cycles of 95 °C for 10 s and 60 °C for 31 s, and the results were analyzed with the relative quantification system based on the 2-ΔΔCt method. Actin primers were 5’-AAGTACAGTGTCTGGATTGGAGGG-3’ (sense) and 5’-TCGCAACTTAGAAGCACTTCCG-3’ (antisense).
Illumina cDNA library preparation and sequencing
For cDNA synthesis and sequencing, 20 μg of total RNA was used, at a concentration of ≥ 400 ng/μL, one technical replicate for both mutant and wild type. The cDNA libraries were constructed using an mRNA-Seq assay with a fragment length range of 200 bp (± 25 bp). The libraries were sequenced for paired-end reads of 90 bp using Illumina HiSeq 2000 platform (Illumina, San Diego, CA, USA) at the Beijing Genomics Institute (Shenzhen, China).
De novo assembly and assessment
Reads from each library were assembled separately. Raw reads that contained adapters and unknown or low-quality bases were discarded by (1) finding the reads containing adapters and to removing the adapters; (2) removing reads with more than 5% unknown nucleotides; and (3) removing low-quality reads using in-house Perl scripts. Transcriptome de novo assembly was conducted with Trinity (version 20130225) (Grabherr et al., 2011). Trinity-assembled reads were grouped to form longer fragments without N, named as contigs. Then, in order to detect contigs in the same transcripts and the distances between them, we mapped reads back to contigs. Next, Trinity connected the contigs using N to represent unknown sequences between each two contigs, and then scaffolds were made. Finally, we got sequences without Ns that could not be extended, which were defined as Unigenes. TGICL (TIGR Gene Indices clustering tools) (Pertea et al., 2003) was then used to assemble all the unigenes from different samples to form a single set of non-redundant unigenes. Single-nucleotide polymorphism (SNP) profiling to compare 23-19 wildtype and eibi1 mutant was done using SOAPsnp (Short Oligonucleotide Alignment Program for SNP detection) (http://soap.genomics.org.cn) (Li et al., 2009) after the Unigene sequence was assembled.
Functional annotation
We identified possible protein coding regions within the assembled Unigenes using the TransDecoder program implemented in the Trinity software and then analyzed protein sequences for Pfam matches to obtain Pfam annotation and alignments (Finn et al., 2016) in Pfam 31.0 database (http://pfam.xfam.org/) (Schaeffer et al., 2017). The pfam2go mapping (http://www.geneontology.org/external2go/pfam2go) was used to map Pfam annotation to GO terms for All-Unigenes.
Analysis of differentially expressed genes
To identify differentially expressed genes (DEG), we examined counts (the number of overlapping reads for each coding region) in the eibi1 and 23-19 libraries by empirical analysis of DEG in EdgeR, which is part of the Bioconductor project (Robinson et al., 2009; McCarthy et al., 2012). Gene expression levels were measured in the RNA-Seq analyses as numbers of reads per kilobase of exon region in a given gene per million mapped reads (RPKM), a normalized measure of exonic read density that allows transcript levels to be compared both within and between samples (Mortazavi et al., 2008). The results of all statistical tests were corrected for multiple testing with the Benjamini–Hochberg false discovery rate (FDR). In our analysis, we chose those with FDR ≤ 0.05 and ratio larger than 2. A gene is defined as absent when its RPKM value is zero.
BinGO (Cytoscape 3.3.0) (Maere et al., 2005) was used for Gene Ontology (GO) functional classification. The calculated p-value was submitted to the Bonferroni correction, using the corrected p-value ≤ 0.05 as a threshold. GO terms fulfilling this condition were defined as significantly enriched GO terms in DEGs.
Identification of EST polymorphisms
The sequences of expressed sequence tags (EST) in the Morex eibi1 (155F2) genetic map were retrieved from the NCBI nucleotide database (http://www.ncbi.nlm.nih.gov/nucleotide/). Corresponding eibi1 Unigenes were identified by local blastn against Morex contigs that were found by blastn with EST sequences as queries against the Morex sequence database: assembly3_WGSMorex_renamed_blastable_POPSEQ.fasta (http://webblast.ipk-gatersleben.de/barley). EST polymorphisms between eibi1 and Morex were identified by the comparison between eibi1 Unigenes and Morex contigs.
Results
The transcriptome of the barley leaf elongation zone
About 54 million and 51 million clean reads were generated from the 23-19 and eibi1 cDNA libraries, respectively (Table 1). A total of 78,799 Unigenes with the average length of 645bp were generated. The size distribution of the assembly is shown in Table S1 available as Supplementary Data. Distribution of unique-mapped reads of the assembled unigenes of 23-19 and eibi1 are shown in Table S2. The entire Unigene sets were then annotated with Pfam database. Among the 78,799 unique sequences, 16,845 Unigenes were annotated using the public databases Pfam with GO terms. BinGO tool was used to perform the enrichment analysis. In the biological process category, 10,077 Unigenes (60%) were involved in cellular and metabolic process. In the cell component category, more than 3,721 Unigenes (22%) played roles in cell, membrane and organelle. In the molecular function category, about 14,717 Unigenes (87%) were involved in metabolic and cellular processes (Table S3, Figure 1).
The identification of eibi1 mutation via RNA-Seq
We compared the sequences between mutant and its isogenic wild type and found 5 unigenes with SNPs (Table S4). Among those, Unigene19902_All (the combination of Unigene20340_23-19 and Unigene66600_eibil) is an ABCG transporter gene which may be involved in cutin secretion (Chen et al., 2009, 2011b). Thus, Unigene19902_All is likely the candidate gene for Eibi1. Actually, Unigene19902_All was the corresponding sequence of Eibi1; it exhibited a SNP exactly as found in our previous study (Figure 2) (Chen et al., 2011b). The depth of the Unigene19902_All was 63 in eibi1 and 320 in wild type. The coverage of this Unigene was 97% and 50% in eibi1 and 23-19, respectively. The base substitution (G → A) at position 294 mapped 56 and 255 reads of the assembled Unigenes of ebi1 and 23-19, respectively (Table S4). This result suggests that the RNA-Seq data can be used for the cloning of a mutant gene.
Identification of eibi1 mutation by RNA-Seq. Unigene20340_23-19 and Unigene66600_eibi1 were assembled from RNA-Seq of wild-type 23-19 and mutant eibi1, respectively. The SNPs are indicated by arrows.
Because a mutant gene was not expressed in the tissue sampled for the RNA-Seq, we could not find the corresponding mutation. However, RNA-Seq data can be used to find polymorphisms in Unigenes between parental lines of an F2 population for genetic mapping of the mutant. As an example, in Morex eibi1 (155F2) genetic map (Figure 3) (Chen et al., 2009), six out of eight EST-derived markers were confirmed by eibi1 Unigenes and Morex contigs (Table S5), indicating that the results of the present RNA-Seq study were reliable.
An example of marker development by RNA-Seq. (A) eibi1 genetic map. Bmag828 and Bmag603 are SSR markers. The others are EST markers. A total of 155 F2 individuals derived from the cross of Morex eibi1 were used for the genetic mapping. The numbers at the left side indicate genetic distances in centimorgan (cM). (B) Corresponding marker polymorphism between parental line Morex and eibi1. The polymorphisms were identified by PCR sequencing against parental genomic DNA. (C) Corresponding eibi1 Unigenes (derived from RNA-Seq) for the EST markers. Inside the parenthesis is a nucleotide followed by its position in the corresponding Unigene sequence. (D) Corresponding Morex contigs for the EST markers were retrieved from the database: assembly3_WGSMorex_renamed_blastable_POPSEQ.fasta (http://webblast.ipk-gatersleben.de/barley). Inside the parenthesis is a nucleotide followed by its position in the corresponding contig sequence.
Differentially expressed genes found in eibi1 mutant
About 269 Unigenes were identified as DEGs between eibi1 and 23-19 samples (Figure 4). Only 110 DEGs were annotated in the Pfam with GO terms (Table S7), which covered only half of differentially expressed Unigenes. To explore DEG results more deeply, within the top 20 up- and down-regulated DEGs that failed Pfam2go mapping were analyzed (Table S8) via barley genome annotation (International Barley Genome Sequencing Consortium 2016) (Beier et al., 2017).Table S6https://minio.scielo.br/documentstore/1678-4685/QqPPK7bzxqYLk49tNxdXxSr/e8fb07bea82a974c93078c40832cabb80966ea47.xls), 148 up-regulated and 121 down-regulated in the mutant (
Expression level of DEGs of 23-19 and eibi1 libraries. All Unigenes were examined for differences in their expression across the two samples. Red (up) and green (down) dots denote that the expression was significantly different (FDR ≤ 0.05, 2-fold difference) and the blank dots indicate no change between the expression levels in the two varieties.
In order to find out if cuticle-related genes were affected by eibi1 mutation, we searched 75 known cuticle genes (Table S9) (Nawrath, 2006; Yeats and Rose, 2013) against the RNA-Seq data by BLAST and examined their expression level. We found 66 genes that had at least one homolog in the RNA-Seq data (E-value e < 0.00001). But no homolog showed significant differential expression between eibi1 and its wild parental line 23-19.
Validation of differentially expressed genes
To assess the accuracy of DEGs, a subset of 10 DEGs in eibi1 and 23-19 were analyzed by quantitative real-time PCR (qRT-PCR). Their primers were designed based on barley sequences obtained from RNA-Seq, and are listed in Table S10. We tested their expression results obtained from qRT-PCR with those generated by the RNA-Seq, which were calculated based on the median of three repeats, and showed good consistency for all transcripts in both analyses, with a correlation coefficient R2 = 0.84638 (Figure 5). Hence, results of the DEG analysis are trustworthy, and its results are suitable for further investigation.
Comparison of log2 (FC) of 10 select transcripts using RNA-Seq and qRT-PCR. FC: fold change for eibi1 vs 23-19.
Functional annotation of differentially expressed genes
As shown in Table 2, Unigene65566_All and Unigene55562_All, annotated as GDSL esterase/lipase, were significantly down-regulated in eibi1. Unigene21078_All encoding Lipase as well as Unigene57532_All and Unigene60434_All encoding glycerophosphodiester phosphodiesterase 1 was slightly down-regulated in mutant. Meanwhile, lipid-transfer protein (Unigene75732_All), lipoxygenase (Unigene72225_All) and fatty acid elongases (Unigene6358_All) were up-regulated in eibi1.
The most up-regulated DEGs were defense-related genes (Table 3), such as pathogenesis-related (PR) genes (Unigene_78697_All and Unigene_78238_All), WIR1A (wheat induced resistance 1) (Unigene66257_All) (Coram et al., 2010), horcolin (Hordeum vulgare coleoptile lectin) (Unigene77099_All) (Grunwald et al., 2007), thionin (Unigene21824_All and Unigene43623_All) (Apel et al., 1992) and jasmonate-induced genes (Unigene78255_All). Moreover, genes involved in defense responses that were also highly up-regulated in eibi1 including Lipoxygenase (LOX) (Unigene72225_All) (Bell and Mullet, 1991; Yang et al., 2012), heat-shock protein (Unigene78052_All) (Vierling, 1991), Cysteine-rich secretory protein family (Unigene68601_All) and RNase S-like protein (Unigene19875_All and Unigene66056_All) (Zheng et al., 2014). But aspartyl protease family protein in insect digestion (Unigene78045_All) (Ryan 1990) and proteinase inhibitor I in response to wounding or insect attacks (Unigene57171_All) (Lee et al., 1986) were down-regulated in eibi1 .
Unigene69289_All and Unigene78082_All encoding peroxidase to defend plants against reactive oxygen species (ROS) (Wang et al., 1999; Mittler, 2002) were clearly up-regulated in eibi1. Unigene33356_All involved in response to water deficit stress (WDS) mediated by ABA was down-regulated (Padmanabhan et al., 1997) (Table 3). In addition, about four DEGs were associated with sugar transportation, seven DEGs were related to protein tyrosine kinase, eight DEGs were involved in cytochrome P450 and another nine DEGs were annotated as cellulose synthase (Table S7).
Discussion
We carried out a transcriptome analysis of the hvabcg31/eibi1 mutant and its isogenic wild barley 23-19 in the leaf elongation zone. About 78,799 Unigenes with average length of 645 bp were assembled and analyzed for gene expression in both mutant and wild type. A SNP filtering analysis between eibi1 and 23-19 Unigenes identified a SNP between Unigene20340_23-19 and Unigene66600_eibi1, which were partial sequences of an ABCG transporter gene. Actually, this ABCG transporter gene was Eibi1 and the SNP was the eibi1 mutation (Chen et al., 2011b). In the Morex eibi1 (155F2) genetic map (Chen et al., 2009), the six marker polymorphisms of all eight EST sequences were confirmed by the comparison of eibi1 Unigenes and Morex contigs (http://webblast.ipk-gatersleben.de/barley). Two markers had polymorphisms in introns, so they could not be confirmed by RNA-Seq data. These results infer that RNA-Seq may be used for mutation identification and/or mutant genetic mapping, similar to SNP calling (Piskol et al., 2013).
It is of interest that both studies failed to identify the most known differentially expressed cuticle genes. Using the second leaves of eibi1 mutant and wild-type seedlings 10 days after germination at 22 °C in growth chamber for transcriptome analysis by a 22-k Barley1 Affymetrix microarray, similar results to the present study were found (Yang et al., 2013). Only two Unigenes and one gene closely linked to cuticle showed differential expression in the present and previous study, respectively. Two GDSL Unigenes are close homologs of cutin synthase-like2 (CUS2) in tomato (Solanum lycopersicum) (Yeats et al., 2014). CUS protein is predicted to catalyze cutin polymerization and is located extracellularly. HvABCG31/Eibi1 is assumed to function in transporting substrates for CUS. The reduced amount of cutin monomers transported out of the epidermal cells in eibi1 mutant (Chen et al., 2011b) may be responsible for the down-regulation of the two GDSL Unigenes. In the previous study, one cytochrome P450 monooxygenase gene (CYP450) homologous to At4g39490 was up regulated. At4g39490 is a paralog of At1G57750 (AtMAH1) that is involved in the cutin monomer synthesis (Greer et al., 2007). The plant cuticle is an extracellular lipid structure. The limited impact on lipid metabolic networks is also found in Arabidopsis (Arabidopsis thaliana) cuticle mutant glossyhead1 (gsd1) (Lü et al., 2011). The eibi1 mutant showed only three lipid genes up-regulated and five lipid genes down-regulated, whereas the gsd1 mutant exhibited three lipid genes up-regulated and eight lipid genes down-regulated (Table 2). The defective cuticle leads to a particularly severe water loss from eibi1 mutant leaves. Thus, eibi1 plants suffer from dehydration stress (Chen et al., 2011b). HvABCG31/Eibi1 was not differentially expressed in the eibi1 mutant in both studies, indicating that the leaf phenotype may be responsible for differentially expressed stress- and defense-related genes. Similarly, the Arabidopsis cuticle mutant gsd1 has a large effect on genes responsible for abiotic and biotic stress (Lü et al., 2011).
eibi1 mutant plants are unable to export cuticular lipids (mostly 16- and 18-carbon fatty acids) from the epidermis cells, leading to an accumulation of intracellular lipids (Chen et al., 2011b). Both 16- and 18-carbon fatty acids are involved in defense to modulate basal, effector-triggered, and systemic immunity in plants (Kachroo and Kachroo, 2009). In a previous study, certain cuticle breakdown monomers were shown to function as elicitors of defense reactions in plants (Kolattukudy 1985; Schweizer et al., 1994, 1996b). Treatment of cutin monomers can trigger alkalinization, higher ethylene (ET), and accumulation of defense-related genes in potato (Schweizer et al., 1996a), and production of H2O2 in cucumber (Fauth et al., 1998). Permeable Arabidopsis cuticle mutants (pec/atabcg32) and rice plants compromised in OsABCG31 expression also showed increased resistance to Botrytis cinerea and Magnaporthe oryzae, respectively (Reina-Pinto and Yephremov, 2009; Bessire et al., 2011; Garroum et al., 2016;). In the study on the osabcg31 knockout mutant, and hairpin RNA interference (RNAi)-down-regulated OsABCG31 plant lines, genes involved in pathogen resistance were constitutively up-regulated. A link to increased expression of defense genes and dwarfism was also revealed (Garroum et al., 2016). A growth reduction has been described as a typical trade-off for plant defenses (Bowling et al., 1994; Yan et al., 2007). In another study, cutinase-expressing plants (CUTE plants) displayed almost complete immunity toward Botrytis cinerea, although this was not found to correlate with the induction of genes coding for various pathogenesis-related (PR) proteins (Chassot et al., 2007). It was also reported that cutin deficiency significantly affected the fruit response to biotic stress and significantly enhanced fruit sensitivity to the fungal post-harvest pathogen Colletotrichum coccodes (Shi et al., 2013). However, not all mutants affected in the cuticle structure showed an enhanced resistance to certain pathogens. It was shown that only certain kinds of cutin monomers induced the expression of defense genes in tomato (Buxdorf et al., 2014; Metraux et al., 2014).
The use of an inducible defense system, named acquired resistance (AR) has been well demonstrated in plants (Heath, 2000; Conrath, 2006; Winter et al., 2014). AR or systemic acquired resistance (SAR) is a form of inducible resistance that is triggered systemically in plants. The mechanism is associated with a systemic accumulation of salicylic acid (SA). Another potential lipid-derived AR signal is the oxylipin-derived defense hormone jasmonic acid (JA), which may be an early signal establishing systemic immunity (Vlot et al., 2008). AR is also associated with the coordinated activation of genes encoding PR proteins (Walters et al., 2014). How does the whole system respond after perceiving signals like these over-produced or over-accumulated cuticle monomers? It is very likely that JA or even other lipid signals play important roles in the generation and/or transmission of mobile signals for AR (Conrath, 2006), although it is also possible that SA initiates the process (Colebrook et al., 2012).
The role played by ABA is intriguing. It has been shown that no significant change was found in the eibi1 mutant (Chen et al., 2004). Yet, ABA-related genes were differentially expressed in the mutant. It was suggested that ABA is likely to be involved in the suppression of wound-induced ROS, when plants are kept under dry conditions (Apel and Hirt, 2004). Several reports showed a link between ABA and increased susceptibility to pathogens that was mostly explained by antagonistic interactions of ABA with defense signaling controlled by SA, JA, or ET (Mauch-Mani and Mauch, 2005).
Overall, this mutant could be a model of self-induced AR (Figure 6). The mutation of eibi1 leads to an excessive amount of monomers or precursors of cutin. Then, possibly through unknown plant immune mechanisms, the plant reacts with increased JA, which further up-regulates PR gene families. This leads to a broad spectrum response, and a metabolic regulation possibly as compensatory process (Metraux et al., 2014). Further applications on molecular mechanisms and agriculture should be studied.
Acknowledgments
This work was supported by grants from the National Natural Science Foundation of China (41621001 and 31170369), the Key Project of the Chinese National Programs for Fundamental Research and Development (grant number 2013CB429904), the China-Israel cooperation program grants from the Ministry of Science and Technology in China (grant number 2013DFA32200). We thank BGI Inc. for contributing to the RNA-Seq analysis, and Mrs. Yan Ma for maintaining plant samples, as well as other helpful contributions.
References
- Apel K, Andresen I, Becker W, Schluter K, Burges J and Parthier B (1992) The identification of leaf thionin as one of the main jasmonate-induced proteins of barley (Hordeum vulgare). Plant Mol Biol 19:193-204.
- Apel K and Hirt H (2004) Reactive oxygen species: Metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol 55:373-399.
- Beier S, Himmelbach A, Colmsee C, Zhang X-Q, Barrero RA, Zhang Q, Li L, Bayer M, Bolser D, Taudien S, et al. (2017) Construction of a map-based reference genome sequence for barley, Hordeum vulgare L Sci Data 4:170044.
- Bell E and Mullet JE (1991) Lipoxygenase gene expression is modulated in plants by water deficit, wounding, and methyl jasmonate. Mol Gen Genet 230:456-462.
- Bessire M, Borel S, Fabre G, Carraça L, Efremova N, Yephremov A, Cao Y, Jetter R, Jacquat A-C and Métraux J-P (2011) A member of the pleiotropic drug resistance family of ATP binding cassette transporters is required for the formation of a functional cuticle in Arabidopsis. Plant Cell 23:1958-1970.
- Bowling SA, Guo A, Cao H, Gordon AS, Klessig DF and Dong XN (1994) A mutation in Arabidopsis that leads to constitutive expression of systemic acquired resistance. Plant Cell 6:1845-1857.
- Buxdorf K, Rubinsky G, Barda O, Burdman S, Aharoni A and Levy M (2014) The transcription factor SlSHINE3 modulates defense responses in tomato plants. Plant Mol Biol 84:37-47.
- Chassot C, Nawrath C and Métraux JP (2007) Cuticular defects lead to full immunity to a major plant pathogen. Plant J 49:972-980.
- Chen GX, Komatsuda T, Ma JF, Li C, Yamaji N and Nevo E (2011a) A functional cutin matrix is required for plant protection against water loss. Plant Signal Behav 6:1297.
- Chen GX, Komatsuda T, Ma JF, Nawrath C, Pourkheirandish M, Tagiri A, Hu Y-G, Sameri M, Li XR and Zhao X (2011b) An ATP-binding cassette subfamily G full transporter is essential for the retention of leaf water in both wild barley and rice. Proc Natl Acad Sci U S A 108:12354-12359.
- Chen GX, Komatsuda T, Pourkheirandish M, Sameri M, Sato K, Krugman T, Fahima T, Korol AB and Nevo E (2009) Mapping of the eibi1 gene responsible for the drought hypersensitive cuticle in wild barley (Hordeum spontaneum). Breed Sci 59:21-26.
- Chen GX, Sagi M, Weining S, Krugman T, Fahima T, Korol A and Nevo E (2004) Wild barley eibi1 mutation identifies a gene essential for leaf water conservation. Planta 219:684-693.
- Colebrook E, Creissen G, McGrann G, Dreos R, Lamb C and Boyd L (2012) Broad-spectrum acquired resistance in barley induced by the Pseudomonas pathosystem shares transcriptional components with Arabidopsis systemic acquired resistance. Mol Plant Microbe In 25:658-667.
- Conrath U (2006) Systemic acquired resistance. Plant Signal Behav 1:179.
- Coram TE, Huang XL, Zhan GG, Settles ML and Chen XM (2010) Meta-analysis of transcripts associated with race-specific resistance to stripe rust in wheat demonstrates common induction of blue copper-binding protein, heat-stress transcription factor, pathogen-induced WIR1A protein, and ent-kaurene synthase transcripts. Funct Integr Genomic 10:383-392.
- Fauth M, Schweizer P, Buchala A, Markstädter C, Riederer M, Kato T and Kauss H (1998) Cutin monomers and surface wax constituents elicit H2O2 in conditioned Cucumber Hypocotyl segments and enhance the activity of other H2O2 elicitors. Plant Physiol 117:1373-1380.
- Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, Potter SC, Punta M, Qureshi M, Sangrador-Vegas A, et al. (2016) The Pfam protein families database: Towards a more sustainable future. Nucleic Acids Research 44:D279-D285.
- Garroum I, Bidzinski P, Daraspe J, Mucciolo A, Humbel BM, Morel J-B and Nawrath C (2016) Cuticular defects in Oryza sativa ATP-binding cassette transporter G31 mutant plants cause Dwarfism, elevated defense responses and pathogen resistance. Plant Cell Physiol 57:1179-1188.
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R and Zeng QD (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29:644-652.
- Greer S, Wen M, Bird D, Wu X, Samuels L, Kunst L and Jetter R (2007) The Cytochrome P450 Enzyme CYP96A15 is the midchain alkane hydroxylase responsible for formation of secondary alcohols and ketones in stem cuticular wax of Arabidopsis. Plant Physiol 145:653-667.
- Grunwald I, Heinig I, Thole HH, Neumann D, Kahmann U, Kloppstech K and Gau AE (2007) Purification and characterisation of a jacalin-related, coleoptile specific lectin from Hordeum vulgare Planta 226:225-234.
- Heath MC (2000) Hypersensitive response-related death. In: Lam E, Fukuda H and Greenberg J (eds) Programmed Cell Death in Higher Plants. Springer, Dordrecht, pp 77-90.
- Kachroo A and Kachroo P (2009) Fatty Acid-derived signals in plant defense. Annu Rev Phytopathol 47:153-176.
- Kolattukudy P (1985) Enzymatic penetration of the plant cuticle by fungal pathogens. Annu Rev Phytopathol 23:223-250.
- Lü S, Zhao H, Parsons EP, Xu C, Kosma DK, Xu X, Chao D, Lohrey G, Bangarusamy DK, Wang G, et al. (2011) The glossyhead1 allele of ACC1 reveals a principal role for multidomain Acetyl-Coenzyme A carboxylase in the biosynthesis of cuticular waxes by Arabidopsis. Plant Physiol 157:1079-1092.
- Lee J, Brown W, Graham J, Pearce G, Fox E, Dreher T, Ahern K, Pearson G and Ryan C (1986) Molecular characterization and phylogenetic studies of a wound-inducible proteinase inhibitor I gene in Lycopersicon species. Proc Natl Acad Sci USA 83:7277-7281.
- Li RQ, Li YR, Fang XD, Yang H, Wang J, Kristiansen K and Wang J (2009) SNP detection for massively parallel whole-genome resequencing. Genome Res 19:1124-1132.
- Maere S, Heymans K and Kuiper M (2005) BiNGO: A Cytoscape plugin to assess overrepresentation of Gene Ontology categories in Biological Networks. Bioinformatics 21:3448-3449.
- Mauch-Mani B and Mauch F (2005) The role of abscisic acid in plant-pathogen interactions. Curr Opin Plant Biol 8:409-414.
- McCarthy DJ, Chen Y and Smyth GK (2012) Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res 40:4288-4297.
- Metraux J, Serrano M, Torres M, Coluccia F and L’Haridon F (2014) The cuticle and plant defense to pathogens. Front Plant Sci 5:274.
- Mittler R (2002) Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci 7:405-410.
- Mortazavi A, Williams BA, McCue K, Schaeffer L and Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5:621-628.
- Nawrath C (2006) Unraveling the complex network of cuticular structure and function. Curr Opin Plant Biol 9:281-287.
- Padmanabhan V, Dias DM and Newton RJ (1997) Expression analysis of a gene family in loblolly pine (Pinus taeda L) induced by water deficit stress. Plant Mol Biol 35:801-807.
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F and Parvizi B (2003) TIGR Gene Indices clustering tools (TGICL): A software system for fast clustering of large EST datasets. Bioinformatics 19:651-652.
- Piskol R, Ramaswami G and Li Jin B (2013) Reliable identification of genomic variants from RNA-seq data. Am J Hum Genet 93:641-651.
- Reina-Pinto JJ and Yephremov A (2009) Surface lipids and plant defenses. Plant Physiol Bioch 47:540-549.
- Richardson A, Wojciechowski T, Franke R, Schreiber L, Kerstiens G, Jarvis M and Fricke W (2007) Cuticular permeance in relation to wax and cutin development along the growing barley (Hordeum vulgare) leaf. Planta 225:1471-1481.
- Robinson MD, McCarthy DJ and Smyth GK (2009) edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26:139-140.
- Ryan CA (1990) Protease inhibitors in plants: Genes for improving defenses against insects and pathogens. Annu Rev Phytopathol 28 425-449.
- Schaeffer RD, Liao Y, Cheng H and Grishin NV (2017) ECOD: New developments in the evolutionary classification of domains. Nucleic Acids Res 45:D296-D302.
- Schweizer P, Felix G, Buchala A, Müller C and Métraux JP (1996a) Perception of free cutin monomers by plant cells. Plant J 10:331-341.
- Schweizer P, Jeanguenat A, Mösinger E and Métraux J-P (1994) Plant protection by free cutin monomers in two cereal pathosystems. In: Daniels, Michael J, Downie, J Allan, Osbourn and Anne E (eds) Advances in Molecular Genetics of Plant-Microbe Interactions. Springer, Netherlands, pp 371-374.
- Schweizer P, Jeanguenat A, Whitacre D, Métraux J-P and Mösinge E (1996b) Induction of resistance in barley against Erysiphe graminis f.sp. hordei by free cutin monomers. Physiol Mol Plant P 49:103-120.
- Shi JX, Adato A, Alkan N, He YH, Lashbrooke J, Matas AJ, Meir S, Malitsky S, Isaacson T and Prusky D (2013) The tomato SlSHINE3 transcription factor regulates fruit cuticle formation and epidermal patterning. New Phytol 197:468-480.
- Vierling E (1991) The roles of heat shock proteins in plants. Ann Rev Plant Biol 42:579-620.
- Vlot AC, Klessig DF and Park S-W (2008) Systemic acquired resistance: The elusive signal (s). Curr Opin Plant Biol 11:436-442.
- Walters DR, Havis ND, Paterson L, Taylor J, Walsh DJ and Sablou C (2014) Control of foliar pathogens of spring barley using a combination of resistance elicitors. Front Plant Sci 5:241.
- Wang J, Zhang H and Allen RD (1999) Overexpression of an Arabidopsis peroxisomal ascorbate peroxidase gen in tobacco increases protection against oxidative stress. Plant Cell Physiol 40:725-732.
- Winter PS, Bowman CE, Villani PJ, Dolan TE and Hauck NR (2014) Systemic acquired resistance in moss: Further evidence for conserved defense mechanisms in plants. PloS One 9:e101880.
- Yan Y, Stolz S, Chételat A, Reymond P, Pagni M, Dubugnon L and Farmer EE (2007) A downstream mediator in the growth repression limb of the jasmonate pathway. Plant Cell 19:2470-2483.
- Yang X-Y, Jiang W-J and Yu H-J (2012) The expression profiling of the lipoxygenase (LOX) family genes during fruit development, abiotic stress and hormonal treatments in cucumber (Cucumis sativus L.). Int J Mol Sci 13:2481-2500.
- Yang Z, Zhang T, Lang T, Li G, Chen G and Nevo E (2013) Transcriptome comparative profiling of barley eibi1 mutant reveals pleiotropic effects of HvABCG31 gene on cuticle biogenesis and stress responsive pathways. Int J Mol Sci 14:20478-20491.
- Yeats TH, Huang W, Chatterjee S, Viart HMF, Clausen MH, Stark RE and Rose JKC (2014) Tomato Cutin Deficient 1 (CD1) and putative orthologs comprise an ancient family of cutin synthase-like (CUS) proteins that are conserved among land plants. Plant J 77:667-675.
- Yeats TH and Rose JK (2013) The formation and function of plant cuticles. Plant Physiol 163:5-20.
- Zheng J, Wang Y, He Y, Zhou J, Li Y, Liu Q and Xie X (2014) Overexpression of an S-like ribonuclease gene, OsRNS4, confers enhanced tolerance to high salinity and hyposensitivity to phytochrome-mediated light signals in rice. Plant Sci 214:99-105.
Supplementary Materials
The following online material is available for this article:
Table S1 - Quality of Unigenes.
Table S2 - Distribution of unique-mapped reads of the assembled unigenes of 23-19 and eibi1.
Table S3 - Functional classification by Gene Ontology (GO) analyses of all Unigenes.
Table S4 - Detailed SNP Statistics.
Table S5 - Polymorphisms were confirmed by eibi1 Unigenes and Morex contigs.
Table S6 - List of DEGs between 23-19 vs eibi1. https://minio.scielo.br/documentstore/1678-4685/QqPPK7bzxqYLk49tNxdXxSr/e8fb07bea82a974c93078c40832cabb80966ea47.xls
Table S7 - DEG annotated by Pfam with GO terms.
Table S8 - DEGs within the top 20 up -and down-regulated genes that failed Pfam2go mapping.
Publication Dates
-
Publication in this collection
23 Oct 2017 -
Date of issue
Oct-Dec 2017
History
-
Received
16 Dec 2016 -
Accepted
13 Aug 2017
