Figures
Abstract
A family of putative PECTIN ESTERASE INHIBITOR (PEI) genes, which were detected in the genomic region co-segregating with the resistance gene Rrs2 against scald caused by Rhynchosporium commune in barley, were characterized and tested for their possible involvement in mediating resistance to the pathogen by complementation and overexpression analysis. The sequences of the respective genes were derived from two BAC contigs originating from the susceptible cultivar ‘Morex’. For the genes HvPEI2, HvPEI3, HvPEI4 and HvPEI6, specific haplotypes for 18 resistant and 23 susceptible cultivars were detected after PCR-amplification and haplotype-specific CAPS-markers were developed. None of the tested candidate genes HvPEI2, HvPEI3 and HvPEI4 alone conferred a high resistance level in transgenic over-expression plants, though an improvement of the resistance level was observed especially with OE-lines for gene HvPEI4. These results do not confirm but also do not exclude an involvement of the PEI gene family in the response to the pathogen. A candidate for the resistance gene Rrs2 could not be identified yet. It is possible that Rrs2 is a PEI gene or another type of gene which has not been detected in the susceptible cultivar ‘Morex’ or the full resistance reaction requires the presence of several PEI genes.
Citation: Marzin S, Hanemann A, Sharma S, Hensel G, Kumlehn J, Schweizer G, et al. (2016) Are PECTIN ESTERASE INHIBITOR Genes Involved in Mediating Resistance to Rhynchosporium commune in Barley? PLoS ONE 11(3): e0150485. https://doi.org/10.1371/journal.pone.0150485
Editor: Dragan Perovic, Julius Kuehn-Institute (JKI), GERMANY
Received: October 6, 2015; Accepted: February 15, 2016; Published: March 3, 2016
Copyright: © 2016 Marzin et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: The research was funded by a project of the Deutsche Forschungsgemeinschaft (DFG), project number Ro 1055/10-1. www.dfg.de.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Barley (Hordeum vulgare L.), one of the most important crop species worldwide is threatened by the hemibiotrophic, haploid ascomycete Rhynchosporium commune (formerly named Rhynchosporium secalis (Oudem.) J.J. Davis) [1–3], which causes a necrotic foliar disease (scald or leaf blotch) on infected plants [4]. Recent studies revealed that the genus Rhynchosporium which consists of two species R. secalis and R. orthosporum can be further divided based on phylogenetic criteria [5]. R. secalis separates into three monophyletic groups that are also characterized by their host spectra: R. secalis infecting rye and triticale, R. commune infecting cultivated barley and other Hordeum spp. as well as Bromus diandrus (brome grass) and R. agropyri, infecting Agropyron spp. [5].
Regarding the resistance of barley against R. commune, research has led to the identification of altogether 16 major resistance genes or alleles (R genes), however only nine of them have been mapped until now [6,7], which include four major resistance loci such as the Rrs1 complex on chromosome 3H with at least 11 known alleles [8,9], the Rrs2 locus on 7HS [10,11], Rrs13 on chromosome 6H [12] and the Rrs15 locus on 2H [13].
Our knowledge is still limited regarding the identity of the mapped resistance genes and their possible functions and resistance mechanisms. For the Rrs1 mediated resistance, prevention of subcuticular hyphae growth has been observed to be a major factor in the defence reaction against R. commune [1]. Jørgensen et al. [14] found that the Rrs2 mediated resistance reaction in the cultivars ‘Digger’ and ‘Osiris’, which both carry the Rrs2 gene, is characterized by the prevention of fungal penetration of the cuticle through an accelerated and larger formation of papillae and halos in the cell walls compared to susceptible cultivars. The authors also suspect the involvement of chemical compounds in the inhibition process.
Pectins account for ca. 30% of the cell walls of dicot and non-graminaceous monocot plants and for between 5–10% of the walls of grasses. They are a family of complex polysaccharides whose backbones are highly enriched in α-1,4-linked-D-galacturonic acid residues, with various degrees of methylesterification [15]. Pectin is secreted in highly methyl-esterified form and subsequently deesterified by pectin methylesterases (PMEs; EC 3.1.1.11). PME activity is controlled by pectin methylesterase inhibitors (PMEI or PEI), which belong to the family of invertase inhibitors [15]. PMEs as well as PMEIs appear to occur as large gene families. In Arabidopsis thaliana, 66 ORFs have been annotated as putative full-length PMEs and 35 ORFs in Oryza sativa [16]. The structural details of PME-PMEI interaction have been studied by crystallographic methods in several plant species and the inhibition apparently occurs through the formation of an inactive 1:1 complex between the two proteins [17–20]. The role of PMEs in plants have been implicated in several processes, such as cell wall extension, fruit maturation, senescence, pollen tube growth and seed germination [15,21,22]. Furthermore, PMEs are involved in the systemic spread of tobacco mosaic virus [23] as well as tobamovirus [24] and seem to stimulate virus-induced gene silencing [25]. An involvement of cell walls and PME action in resistance to pathogens has been discussed [15,16,26–28].
Our previous results of fine mapping, in which we used a segregating population ‘Atlas’ (resistant) × ‘Steffi’ (susceptible) of 9,179 F2 plants, assigned the scald resistance gene Rrs2 to an area of suppressed recombination on chromosome 7HS [29]. Based on these fine mapping results we decided to carry out a sequence survey of this suppressed recombination region in order to pinpoint interesting genes depicting haplotype variation between resistant and susceptible cultivars. Analysis of several genes, identified within the region of interest, led to the discovery of a family of putative PECTIN ESTERASE INHIBITORs (PEI) whose haplotypes differed between resistant and susceptible cultivars. This feature and the fact that these genes co-segregated with the resistance gene Rrs2 made them possible candidates for the Rrs2 gene. The goal of the present research work was to elucidate the role of the putative PECTIN ESTERASE INHIBITOR (PEI) genes, in the process of mediating resistance to R. commune in barley.
Results
The PEI-gene family in barley
The Rrs2 gene which confers resistance to R. commune was fine mapped by Hanemann et al. [29]. It is located in a co-segregating genomic region on barley chromosome 7HS. In previous research we had established two BAC contigs for the co-segregating region, which consisted of BACs from the susceptible cultivar ‘Morex’ [29]. In this genomic region a number of genes had been identified including genes coding for putative pectinesterase inhibitor domain containing proteins [29]. Further analysis of BAC sequences that became available during the past two years [30] led to the discovery of additional PEI family members in the neighborhood of already known PEI candidates (Fig 1). Gene ontology description predicted enzyme inhibitor activity (GO:0004857) for pectin esterase activity (GO:0030599) for genes HvPEI1, HvPEI2, HvPEI3, HvPEI4 and HvPEI6, therefore we refer to this gene family as putative pectin esterase inhibitor genes. For HvPEI6 additionally an ABC transporter related InterPro description was identified with molecular functions of ATP binding (GO:0005524) and ATPase activity (GO:0016887), while for HvPEI5 no gene ontology function was identified (S1 File).
The horizontal bar represents the genetic map based on 9,179 F2-plants with the numbers of recombinant plants between six genetic markers indicated according to [29]. The hatched areas represent the location of recombination events and could also harbor the Rrs2-gene. The blue boxes represent BAC-clones and the green arrows indicate the approximate locations of the PEI genes.
Cluster analysis revealed a defined structure of the family of PECTIN ESTERASE INHIBITOR genes found in the Rrs2 region. By combining database information from different sources (e.g. http://www.harvest-web.org), we were able to identify corresponding expressed unigenes or cDNA clones for a total of six PEI genes (Fig 2).
Phylogenetic analysis of six detected PEI family members was performed by using database information of known expressed unigenes (HarvEST Assembly 35), BAC sequence information as well as sequence information from 41 analyzed barley cultivars. In addition data from the Gene Index Project (Barley database 12.0) were also included.
Sequencing of all identified PEI genes revealed different haplotypes harboring characteristic SNPs which can be used to distinguish between resistant and susceptible cultivars of barley (Fig 3, S1 Table, S2 File) for the genes HvPEI2, HvPEI3, HvPEI4 and HvPEI6. One PEI candidate gene, HvPEI4, also showed insertions as well as deletions in the open reading frame (ORF) in addition (S1 Fig). This could cause a frame shift using the example of cv. ‘Morex’ or an introduction of a premature stop codon in the case of cv. ‘Steffi’. From the obtained sequence information, haplotype-specific CAPS-markers were developed for the genes HvPEI2, HvPEI3, HvPEI4, HvPEI6 and two additional INDEL-markers for HvPEI4 (Table 1, Fig 4).
Full-length sequences of all detected PECTIN ESTERASE INHIBITOR members were analyzed in different barley varieties either carrying the Rrs2 resistance gene or not. Information regarding susceptibility to Rhynchosporium commune is according to [29].
As a result of haplotype analysis various SNPs were discovered within the ORF of HvPEI2, HvPEI3, HvPEI4 and HvPEI6. Polymorphisms were used for CAPS-marker development that in turn could be applied for haplotype discrimination. Examples on the gel photo illustrate results of PCR amplification as well as digestion of the amplicon, with the appropriate restriction enzyme (given in brackets), for different PEI family members. As an exception for HvPEI6 only the results of digestion are presented. Primers that were used for each PCR reaction are given in brackets as well. Additional information regarding distinguished haplotypes is shown as abbreviation on top of lanes with digested fragments (e.g. H1 = haplotype 1, H2 = haplotype 2, et cetera). M = 1kb+ DNA Ladder (Thermo Scientific™ GeneRuler, catalog number #SM1333).
As exceptional case, HvPEI5 showed no differences between resistant and susceptible cultivars and hence was not included in further analyses. Although the ViroBLAST (IPK Barley BLAST Server, http://webblast.ipk-gatersleben.de/barley/) search assigned the gene to chromosome 7H we were not able to locate it within the investigated chromosomal region (Fig 1). The gene HvPEI1 had a presumed susceptible variety-specific haplotype H2 in the Rrs2 carrying (resistant) variety ‘Atlas46’ and a presumed resistant variety-specific haplotype H1 in several susceptible varieties and was therefore also excluded from further analysis.
Analysis of transgenic plants
For over-expression of identified PEI candidate genes, the barley cv. ‘Golden Promise’ (GP) was used because of its amenability to genetic transformation and its susceptibility to Rhynchosporium infection. The full-length coding sequences of HvPEI2, HvPEI3, HvPEI4 of resistant cv. ‘Atlas’ were first cloned for each candidate separately into the binary vector p6U and afterwards stably transformed by using Agrobacterium-based transformation technology. For HvPEI2, HvPEI3 and HvPEI4, 25, 41 and 17 independent primary transgenic plants were generated, respectively. For all investigated HvPEI2, HvPEI3 and HvPEI4 over-expression (OE)-plants, a significant increase in expression of the respective genes was confirmed by quantitative RT_qPCR in transgenic progeny in comparison to non-transgenic control plants (cv.’Golden Promise’). In several plants with high copy numbers, such as PEI2_OE BG361-E15 lines, the over-expression levels were reduced, which may be an effect of gene silencing. (S2 Fig). The gene HVPEI6 was only discovered later during the analysis of the respective genomic region; therefore for this gene copy only sequence and haplotype analysis was conducted, but it had not been included in the series of production of transgenic lines which had been finished for the other genes at the time point of discovery of HvPEI6.
The copy number in the progeny was assessed by applying TaqMan™ PCR using a probe specific for the HPT gene. In a T0-plant each transformed gene is usually only present as a single copy per integration site, this means it is present in heterozygous state. Therefore the T1-generation of transformed plants is segregating with respect to copy number, because the transformed sites are not genetically fixed yet. The presence of the transformed PEI genes from ‘Atlas’ was furthermore verified by applying gene and allele-specific CAPS-marker analysis.
Several T1 OE-plants for HvPEI2, HvPEI3 and HvPEI4 were tested for resistance to R. commune isolate ‘LfL 12F’ in an infection test in the greenhouse at the ‘Bayerische Landesanstalt’ in Freising, Germany. Overall, none of the transformed lines reached the resistance level of the resistant cultivar ‘Atlas’ of 0.78 (Table 2). In a disease scoring range of 0 (most resistant) to 4 (most susceptible) the transformed OE-plants ranged from 2.79 to 3.86 on average. Usually a scoring from 0 to 2 is considered as resistant, while a scoring of 3 and 4 is considered as susceptible according to Jackson and Webster [31]. Comparing the transformed OE-lines to non-transformed ‘Golden Promise’, with an average resistance level of 3.84, several of the transformed OE-lines were significantly different from ‘Golden Promise’ (Table 2). Thereby, four out of five OE-lines for HvPEI4 reached a significance level of P<0.001. However, none of the transformed plants reached a scoring level to be considered as resistant.
From these results, it can be concluded that none of the HvPEI2, HvPEI3 and HvPEI4 genes alone seems to represent the Rrs2 gene of ‘Atlas’, since the resistance level of the OE-lines is rather still in the susceptible range. Some improvements of resistance could be observed, especially with construct HvPEI4, however this may be attributed to somaclonal variation. Since the investigated PEI genes originate from the susceptible cultivar ‘Morex’, which does not carry Rrs2, it is possible that in the resistant cultivar ‘Atlas’ additional genes are present, which were not detected and therefore not tested here. These either could be additional PEI genes or other genes which are not present in ‘Morex’.
Since the individual plants of the T1-events are still segregating the estimation of copy numbers of the respective gene copies from ‘Atlas’ are indicated: mc = multi-copy, 1c = one copy. The disease scoring ranged from 0 (most resistant) to 4 (most susceptible).
Discussion
Involvement of pectin esterase inhibitors in plant defense reactions
Several reports indicate a potential role of PEI or PMEI genes in plant defense reactions against pathogens:
The PECTIN METHYL ESTERASE INHIBITOR gene CaPMEI1 from pepper (Capsicum annuum L.) was isolated and functionally characterized [32,33]. Virus-induced gene silencing of CaPMEI1 in pepper conferred enhanced susceptibility to Xanthomonas campestris pv. vesicatoria. Ectopic expression of CaPMEI1 in Arabidopsis led to enhanced resistance to Pseudomonas syringae pv. tomato, but not to the biotrophic pathogen Hyaloperonospora parasitica [32]. Furthermore, the promoter of CaPMEI1 induced a high level of GUS reporter activity in tobacco leaf tissue, driven by pathogen infection as well as by ethylene and methyl jasmonate treatment, suggesting an involvement in host defense response [33].
The over-expression of PMEIs in Arabidopsis resulted in an increase of 16% of pectin methylesterification in cell walls of the transformed plants [34]. Transformed plants were more resistant to the necrotrophic fungus Botrytis cinerea. The reduced symptoms caused by the fungus on transgenic plants were related to its impaired ability to grow on methylesterified pectins [34].
Differences in the methyl ester distribution of homogalacturons were reported for near-isogenic wheat lines resistant and susceptible to the wheat stem rust fungus Puccinia graminis f. sp. Tritici, which may result in differences for the enzymatic generation of endogenous suppressors resulting in different outcomes of the host-pathogen interaction [35].
The fruit-specific Fragaria × ananassa PECTIN METHYL ESTERASE FaPME1 was expressed in wild strawberry Fragaria vesca, resulting in a reduced degree of esterification of the oligogalacturonides from the transgenic fruits as compared to the wild type fruits [36]. The transgenic F. vesca lines had constitutively activated pathogen defense response, resulting in higher resistance to the necrotrophic fungus Botrytis cinerea. It was speculated if oligogalacturonides released from the cell wall may play a role in eliciting defense responses [36].
Transgenic durum wheat lines expressing the PMEI from Actinidia chinensis (AcPMEI) showed a significant increase in the degree of methyl esterification and a significant reduction of disease symptoms caused by the fungal pathogens Bipolaris sorokiniana and Fusarium graminearum [37]. The expression of the TaPMEI gene in wheat was monitored by RT-PCR and high transcription was detected in salicylic acid and hydrogen peroxidase treatments [38]. An enzyme assay indicated that PME was completely inhibited by TaPMEI.
From the previous examples, two possibilities for the involvement of PMEIs in plant defense mechanism become obvious. The first is an increase of PMEI resulting in repressed PME activity, reduced de-methylesterification and therefore in highly esterified cell walls, which are less prone to the attack by pectin-degrading enzymes of potential pathogens [32,34,37]. Another possibility is the involvement of elicitor molecules which may be released by partial degradation of cell walls [36] and subsequent elicitation or enhancement of the plant defense response. In the given examples, the plant PMEIs seem to interact with the own host PMEs. To our knowledge, so far no example exists where a plant PMEI interacts with a fungal or microbial PME.
The concept of ‘effector-triggered defence’ (ETD) has been suggested for apoplastic fungal pathogens, such as R. commune [39]. ETD responses are characterized by a relatively slow response to the pathogen, not associated to a fast hypersensitive host cell death response, whereby effectors of apoplastic pathogens are recognized at the cell surface [38].
Conclusions from own experiments
In previous research, the Rrs2 gene was placed in a linkage block of suppressed recombination which was found in all investigated varieties carrying Rrs2 and indicated a common origin for this resistance gene [29]. This linkage block had previously been defined by SNPs of non-specified genomic sequences [29]. The present results of haplotype analysis confirmed that genes HvPEI2, HvPEI3, HvPEI4 and HvPEI6 are present within the linkage block, which made these genes possible candidates for Rrs2 (Table 1) and the testing of this hypothesis in over-expression plants was initiated. As result, none of the tested candidate genes HvPEI2, HvPEI3 and HvPEI4 alone conferred a high resistance level in transgenic over-expression plants, though an improvement of the resistance level was observed especially with OE-lines for gene HvPEI4 (Table 2). Therefore, we consider the following hypotheses about the identity of the Rrs2 gene:
- Rrs2 is identical to HvPEI6 which has not been tested yet in a transgenic approach.
- The resistance of Rrs2 requires the action of two or more PEI genes.
- Since all available sequence information comes from the susceptible cultivar ‘Morex’, it is possible that in resistant cultivars, such as ‘Atlas’ additional PEI genes or other genes are present in the Rrs2 region which cannot be found in ‘Morex’.
In conclusion, the question posed in the title cannot conclusively be answered by the current results. We can conclude, that neither HvPEI2, HvPEI3 or HvPEI4 alone represents the resistance gene Rrs2. However the presence of additional PEI genes in the resistant cultivar ‘Atlas’ cannot be excluded. It is also possible that several linked genes are required for resistance expression.
Materials and Methods
Plant material and fungal isolates
For haplotype analysis of detected PECTIN ESTERASE INHIBITOR genes within the genomic Rrs2 region of barley chromosome 7HS, 41 different accessions of cultivated barley were chosen based on reported Rrs2 resistance phenotypes. The majority of varieties were taken from Hanemann et al. [29], whereas the following accessions, ‘Anakin’, ‘Atlas46’ (HOR3876), ‘Hobbit’, ‘Landi’, ‘Gloriacome’, ‘CI452395’ (Rrs2+) as well as ‘Barke’, ‘CI11258’, ‘Golden Promise’ (Rrs2-) were taken in addition.
Several fungal isolates were employed for practical reasons. These included Sachs 147–1 and LfL12F (both obtained from Günther Schweizer, LfL Freising, Germany) and UK7 (obtained from Wolfgang Knogge, IPB Halle, Germany). In all experiments, the two barley varieties ‘Atlas’ (carrying Rrs2) and ‘Steffi’ (susceptible), which were the parents of the original mapping population [29], were included as controls.
Haplotype analysis and marker development
New putative PEI genes were identified using two different strategies: First, a keyword search was performed within the HarvEST barley database (http://www.harvest-web.org), version 1.83 (Assembly 35) and sequences of resulting candidates were verified for chromosome position with the ViroBLAST tool (http://webblast.ipk-gatersleben.de/barley/viroblast.php), applying default BLAST parameter settings [40]. The second approach was based on a software application, an automated annotation system for rice: RiceGAAS ([41]; http://ricegaas.dna.affrc.go.jp/), which possesses programs for gene prediction as well as BLAST search options for identification of homologous sequences and gene ontology predictions. Available sequence information of different BAC contigs was analyzed regarding gene content using the graphical user interface.
All PEI candidate genes were amplified by PCR from genomic DNA of cultivar ‘Atlas’ (Rrs2+, resistant) using specific forward and reverse primers (S2 Table). Primer development was done with the program PrimerSelect™ which is part of the software suit DNASTAR® Lasergene version 10 (DNASTAR, Inc., Madison, Wisconsin, USA).
Haplotype analysis and alignment of nucleotide sequences was performed with the program Sequencher™ version 5.2.3 (Gene Codes Corporation, Ann Arbor, MI, USA). Alignments were analyzed for presence of restriction sites at SNP positions with the Sol Genomics Network tool (http://solgenomics.net/tools/caps_designer/caps_input.pl) [42].
Restriction digestions of PCR products were carried out in a 10-μl volume using 5 μl of PCR product, 1U of the respective restriction endonuclease and an incubation time of 2h at recommended optimal temperature. The restricted fragments were separated on 1–2% 1 x TAE agarose gels.
Preparation of over-expression constructs
Sequences of candidate PEI genes were amplified by PCR from genomic DNA of cultivar ‘Atlas’ (Rrs2+, resistant) using specific forward and reverse primers (S2 Table) and cloned into TOPO vector pCR2.1 (Invitrogen, Karlsruhe, Germany) according to manufacturer’s instruction manual for sequence verification. Afterwards correct sequences were transferred into the entry vector pUBI-ABM (DNA-Cloning-Service, Hamburg, Germany) via a digestion and ligation-based cloning strategy that ensured a proper in-frame orientation of PEI candidates. Clones that were tested positive for presence of the gene of interest were used for further cloning. A fragment that includes the maize UBIQUITIN1 promoter, the gene of interest as well as a transcription termination sequence (NOPALINE SYNTHASE terminator) was transferred into destination vector p6U (DNA-Cloning-Service, Hamburg, Germany) using the rare cutter SfiI, which has a 13 nucleotide recognition site including five variable nucleotides in central position [43]. Barley ‘Golden Promise’ immature embryos were stably transformed with PEI over-expression vectors using Agrobacterium tumefaciens according to [44,45].
Transcript abundance and fungal growth (RT-qPCR)
The assessment of candidate gene expression as well as fungal growth, on stably transgenic barley plants (OE), was performed by quantification via reverse transcription real-time PCR as described by Himmelbach et al. [46]. For this purpose, total RNA and DNA were extracted from control or treated primary leaves as well as first and second leaves at different time points after inoculation with R. commune spores. RNA isolation was performed using TRIzol Reagent (Life Technologies GmbH, Darmstadt, Germany), which is a monophasic solution of phenol, guanidine isothiocyanate and other components and which allows to perform sequential precipitation of RNA, DNA and proteins from a single sample [47]. After homogenization of the samples with TRIzol Reagent, chloroform was added and the homogenate was allowed to separate into a clear upper aqueous layer containing RNA, an interphase, and a red lower organic layer containing DNA and proteins. RNA was precipitated from the aqueous layer using isopropanol. The clear stratification in diverse phases prevented the contamination with genomic DNA. In case of analysis of target gene expression, cDNA was synthesized from 1 μg total RNA by using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA) in a 20 μl reaction.
An aliquot of 1 μl cDNA or DNA was taken as template for real-time PCR on a 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) by using the GoTaq® qPCR Master Mix Kit (Promega) according to the manufacturer’s manual, except that the total reaction volume was reduced to 10 μl. The primer pairs specific for the respective genes of interest (qRT2-PEI2_F/R, qRT-PEI3_F7/R7, qRT-PEI4_Fa/R) were designed by the software program Lasergene (DNASTAR) (S3 Table) and deployed in a final concentration of 1 μM. For amplification of DNA, a standard two-step, 40-cycle program, followed by a dissociation step was used. After initial incubation at 50°C (2 min) and extended activation of the polymerase at 95°C (10 min), 40 cycles with 95°C (15 sec) and 62°C (30 sec) were applied.
Transcript as well as DNA abundance was determined relative to standard curves and normalized to a gene encoding the barley UBIQUITIN-CONJUGATING ENZYME (UBC). Data analysis was performed with the sequence detection software (Version 2.2.2; Applied Biosystems) that is provided with the Fast Real-time cycler. Statistical validation was ensured by three biological replicates.
The verification of transgenic barley plants was carried out both via TaqMan® PCR for copy number estimation and RT-qPCR for candidate gene expression. The latter method was already described above whereas the multiplex real-time PCR (TaqMan®) was performed according to Yang et al. [48] by adapting the SYBR Green assay to TaqMan® chemistry: The amplification reaction (15 μl) contained 1 μl genomic DNA, 7.5 μl Maxima Probe/ROX qPCR Master Mix (Thermo Scientific), and forward and reverse primers that are specific for the HYGROMYCIN PHOSPHOTRANSFERASE selectable marker gene (HPT) or the barley UBC gene, respectively. Probes were dual-labeled with a fluorescent reporter dye attached to the 5’-end of the oligonucleotide and the highly efficient dark quencher BHQ1 (Biosearch Technologies, Inc., Petaluma, California, USA) attached to the 3’-end (S3 Table). In case of HPT the 6-carboxy-fluoroscein (FAM™) was used, whereas for UBC 6-carboxy-4',5'-dichloro-2',7'-dimethoxyfluoresceine (JOE™) was applied. Final primer concentration was adjusted to 0.5 μM and probe concentration to 0.2 μM.
Amplification conditions were derived from a two-step, 40-cycle program described above that is complemented with a third step at 72°C (30 sec), while omitting the dissociation step. Each sample was quantified in triplicate.
To calculate transgene copy numbers, a quantitative method was used which allowed comparing the absolute quantification reaction for the candidate gene HPT with the endogenous reference gene UBC. Transgene copy numbers were estimated using a homozygous HPT-transgenic control plant, verified as described by [44], as calibrator sample and UBC as endogenous control. Each reaction comprised three technical replicates. For statistical analysis the software package SigmaPlot11.0 was used.
Supporting Information
S1 Fig. Haplotypes of HvPEI3, HvPEI4 and HvPEI6.
https://doi.org/10.1371/journal.pone.0150485.s001
(PDF)
S2 Fig. Transcript abundances in OE-lines for HvPEI2, HvPEI3 and HvPEI4 in comparison to control lines.
https://doi.org/10.1371/journal.pone.0150485.s002
(PDF)
S1 File. Gene ontology annotations obtained for genes HvPEI1 –HvPEI6.
https://doi.org/10.1371/journal.pone.0150485.s003
(PDF)
S2 File. Sequences of genes HvPEI1 –HvPEI6 and respective haplotypes.
SNPs are labeled in red.
https://doi.org/10.1371/journal.pone.0150485.s004
(PDF)
S1 Table. Haplotypes and SNPs for five PEI genes sequenced in 18 resistant and 23 susceptible cultivars.
https://doi.org/10.1371/journal.pone.0150485.s005
(XLS)
S2 Table. Primer for amplification of full-length sequences of six PEI genes in barley.
https://doi.org/10.1371/journal.pone.0150485.s006
(PDF)
S3 Table. Primer for assessment of transcript abundance transgene copy number determination.
https://doi.org/10.1371/journal.pone.0150485.s007
(PDF)
Acknowledgments
We thank Dr. Wolfgang Knogge for providing Rhynchosporium isolates and sharing expert knowledge and Dr. Patrick Schweizer and Daniela Nowara for help in setting up the VIGS system. Ellen Weiß, Carola Bollmann and Sibylle Freist provided excellent technical assistance.
Author Contributions
Conceived and designed the experiments: MRS JK. Performed the experiments: SM AH SS GH GS. Analyzed the data: SM. Wrote the paper: SM MRS.
References
- 1. Lehnackers H, Knogge W. Cytological studies on the infection of barley cultivars with known resistance genotypes by Rhynchosporium secalis. Can J Bot 1990; 68: 1953–1961.
- 2. Shipton WA, Boyd WJR, Ali SM. Scald of barley. Rev Plant Path 1974; 53: 839–861.
- 3. Zaffarano PL, McDonald BA, Zala M, Linde CC. Global hierarchical gene diversity analysis suggests the fertile crescent is not the center of origin of the barley scald pathogen Rhynchosporium secalis. Phytopathology 2006; 96: 941–950. pmid:18944049
- 4. Avrova A, Knogge W. Rhynchosporium commune: a persistent threat to barley cultivation. Mol Plant Pathol 2012; 13: 986–997. pmid:22738626
- 5. Zaffarano PL, McDonald BA, Linde CC. Two new species of Rhynchosporium. Mycologia 2011; 103: 195–202. pmid:20943529
- 6. Zhan J, Fitt BDL, Pinnschmidt HO, Oxley SJP, Newton AC. Resistance, epidemiology and sustainable management of Rhynchosporium secalis populations on barley. Plant Pathology 2008; 57: 1–14.
- 7. Wagner C, Schweizer G, Krämer M, Dehmer-Badani AG, Ordon F, Friedt W. The complex quantitative barley-Rhynchosporium secalis interaction: newly identified QTL may represent already known resistance genes. Theor Appl Genet 2008; 118: 113–122. pmid:18806993
- 8. Bjørnstad Å, Patil V, Tekauz A, Marøy AG, Skinnes H, Jensen A, et al. Resistance to scald (Rhynchosporium secalis) in barley (Hordeum vulgare) studied by near-isogenic lines: I. markers and differential isolates. Phytopathology 2002; 92: 710–720. pmid:18943266
- 9. Hofmann K, Silvar C, Casas AM, Herz M, Büttner B, Gracia MP, et al. Fine mapping of the Rrs1 resistance locus against scald in two large populations derived from Spanish barley landraces. Theor Appl Genet 2013; 126: 3091–3102. pmid:24068343
- 10. Schweizer G, Baumer M, Daniel G, Rugel H, Röder M. RFLP markers linked to scald (Rhynchosporium secalis) resistance gene Rh2 in barley. Theor Appl Genet 1995; 90: 920–924. pmid:24173045
- 11. Schmidt D, Röder MS, Dargatz H, Wolf N, Schweizer GF, Tekauz A, et al.Construction of a YAC library from barley cultivar Franka and identification of YAC-derived markers linked to the Rh2 gene conferring resistance to scald (Rhynchosporium secalis). Genome 2001; 44: 1031–1040. pmid:11768206
- 12. Genger RK, Williams KJ, Raman H, Read BJ, Wallwork H, Burdon JJ, et al. Leaf scald resistance genes in Hordeum vulgare and Hordeum vulgare ssp. spontaneum: parallels between cultivated and wild barley. Aust J Agric Res 2003; 54: 1335–1342.
- 13.
Schweizer G, Herz M, Mikolajewski S, Brenner M, Hartl L, Baumer M.Genetic mapping of a novel scald resistance gene Rrs15CI8288 in barley. Proceedings of the 9th International Barley Genetics Symposium, Brno, Czech Republic, 2004, 20–26 June. pp. 258–265.
- 14. Jørgensen HJL, Neergaard ED, Smedegaard-Petersen V. Histological examination of the interaction between Rhynchosporium secalis and susceptible and resistant cultivars of barley. Physiol Mol Plant Pathol 1993; 42: 345–358.
- 15. Juge N. Plant protein inhibitors of cell wall degrading enzymes. Trends Plant Sci 2006; 11: 359–366. pmid:16774842
- 16. Pelloux J, Rustérucci C, Mellerowicz EJ. New insights into pectin methylesterase structure and function. Trends Plant Sci 2007; 12: 267–277. pmid:17499007
- 17. Balestrieri C, Castaldo D, Giovane A, Quagliuolo L, Servillo L. A glycoprotein inhibitor of pectin methylesterase in kiwi fruit (Actinidia chinensis). Eur J Biochem 1990; 193: 183–187. pmid:2226435
- 18. Giovane A, Servillo L, Balestrieri C, Raiola A, D’Avino R, Tamburrini M, et al. Pectin methyl esterase inhibitor. Biochimica et Biophysica Acta 2004; 1696: 245–252. pmid:14871665
- 19. Hothorn M, Wolf S, Aloy P, Greiner S, Scheffzek K. Structural insights into target specifity of plant invertase and pectin methylesterase inhibitory proteins. Plant Cell 2004; 16: 3437–3447. pmid:15528298
- 20. Di Matteo A, Giovane A, Raiola A, Camardella L, Bonivento , De Lorenzo G, et al. Structural basis for the interaction between pectin methylesterase and a specific inhibitor protein. Plant Cell 2005; 17: 849–858. pmid:15722470
- 21. Bosch M, Hepler PK. Pectin methylesterases and pectin dynamics in pollen tubes. Plant Cell 2005; 17: 3219–3226. pmid:16322606
- 22. Müller K, Levesque-Tremblay G, Bartels S, Weitbrecht K, Wormit A, Usadel B, et al. Demethylesterification of cell wall pectins in Arabidopsis plays a role in seed germination. Plant Phys 2013; 161: 305–316.
- 23. Chen M-H, Citovsky V. Systemic movement of a tobamovirus requires host cell pectin methylesterase. Plant J 2003; 35: 386–392. pmid:12887589
- 24. Lionetti V, Raiola A, Cervone F, Bellincampi D. Transgenic expression of pectin methylesterase inhibitors limits tobamovirus spread in tobacco and Arabidopsis. Mol Plant Pathol 2014; 15: 265–274. pmid:24127644
- 25. Dorokhov YL, Frolova OY, Skurat EV, Ivanov PA, Gasanova TV, Sheveleva AA, et al. The novel function for a ubiquitious plant enzyme pectin methylesterase: The enhancer of RNA silencing. FEBS Letters 2006; 580: 3872–3878. pmid:16797009
- 26. Vorwerk S, Somerville S, Somerville C. The role of plant cell wall polysaccharide composition in disease resistance. Trends Plant Sci 2004; 9: 203–209. pmid:15063871
- 27. Cantu D, Vicente AR, Labavitch JM, Bennett AB, Powell ALT. Strangers in the matrix: plant cell walls and pathogen susceptibility. Trend Plant Sci 2008; 13: 610–617.
- 28. Lionetti V, Cervone F, Bellincampi D. Methyl esterification of pectin plays a role during plant-pathogen interactions and affects plant resistance to diseases. Journal of Plant Phys 2012; 169: 1623–1630.
- 29. Hanemann A, Schweizer GF, Cossu R, Wicker T, Röder MS. Fine mapping, physical mapping and development of diagnostic markers for the Rrs2 scald resistance gene in barley. Theor Appl Genet 2009; 119: 1507–1522. pmid:19789848
- 30. Ariyadasa R, Mascher M, Nussbaumer T, Schulte D, Frenkel Z, Poursarebani N, et al. A sequence-ready physical map of barley anchored genetically by two million single-nucleotide polymorphisms. Plant Phys 2014; 164: 412–423.
- 31. Jackson LF, Webster RK. Race differentiation, distribution, and frequency of Rhynchosporium secalis in California. Phytopathology 1976; 66: 719–725.
- 32. An SH, Sohn KH, Choi HW, Hwang IS, Lee SC, Hwang BK. Pepper pectin methylesterase inhibitor protein CaPMEI1 is required for antifungal activity, basal disease reistsance and abiotic stress tolerance. Planta 2008; 228: 61–78. pmid:18327607
- 33. An SH, Choi HW, Hong JK, Hwang BK. Regulation and function of the pepper pectin methylesterase inhibitor (CaPMEI1) gene promoter in defense and ethylene and methyl jasmonate signaling in plants. Planta 2009; 230: 1223–1237. pmid:19777255
- 34. Lionetti V, Raiola A, Camardella L, Giovane A, Obel N, Pauly M, et al. Overexpression of pectin methylesterase inhibitors in Arabidopsis restricts fungal infection by Botrytis cinerea. Plant Phys 2007; 143: 1871–1880.
- 35. Wiethölter N, Graeßner B, Mierau M, Mort AJ, Moerschbacher BM. Differences in the methyl ester distribution of homogalacturonans from near-isogenic wheat lines resistant and susceptible to the wheat stem rust fungus. Mol Plant-Microbe Interact 2003; 10: 945–952.
- 36. Osorio S, Castillejo C, Quesada MA, Medina-Escobar Brownsey GJ, Suau R, et al. Partial demethylation of oligalacturonides by pectin methyl esterase 1 is reqired for eliciting defence response in wild strawberry (Fragaria vesca). Plant J 2008; 54: 43–55. pmid:18088306
- 37. Volpi C, Janni M, Lionetti V, Bellincampi D, Favaron F, D’Ovidio R. The ectopic expression of a pectin methyl esterase inhibitor increases pectin methyl esterification and limits fungal diseases in wheat. Mol Plant-Microbe Interaction 2011; 24: 1012–1019.
- 38. Hong MJ, Kim DY, Lee TG, Jeon WB, Seo YW. Functional characterization of pectin methylesterase inhibitor (PMEI) in wheat. Genes Genet Syst 2010; 85: 97–106. pmid:20558896
- 39. Stotz HU, Mitrousia GK, de Wit PJGM, Fitt BDL. Effector-triggered defence against apoplastic fungal pathogens. Trends Plant Sci 2014; 19: 491–500. pmid:24856287
- 40. Deng W, Nickle DC, Learn GH, Maust B, Mullins JI. ViroBLAST: A stand-alone BLAST web server for flexible queries of multiple databases and user´s datasets. Bioinformatics 2007; 23: 2334–2336. pmid:17586542
- 41. Sakata K, Nagamura Y, Numa H, Antonio BA, Nagasaki H, Idonuma A, et al. RiceGAAS: an automated annotation system and database for rice genome sequence. Nucleic Acids Res 2002; 30: 98–102. pmid:11752265
- 42. Fernandez-Pozo N, Menda N, Edwards JD, Saha S, Tecle IY, Strickler SR, et al. The Sol Genomics Network (SGN)–from genotype to phenotype to breeding. Nucl Acids Res 2015; 43: D1036–D1041. pmid:25428362
- 43.
Kumlehn J. The IPKb Vector Set: Modular Binary Plasmids for Cereal Transformation. ISB News Report. 2008.
- 44. Hensel G, Valkov V, Middlefell-Williams J, Kumlehn J; Efficient generation of transgenic barley: the way forward to modulate plant-microbe interactions. J Plant Physiol 2008; 165: 71–82. pmid:17905476
- 45. Hensel G, Kastner C, Oleszczuk S, Riechen J, Kumlehn J. Agrobacterium-mediated gene transfer to cereal crop plants: current protocols for barley, wheat, triticale, and maize. International Journal of Plant Genomics 2009; Article ID 835608, 9 pages.
- 46. Himmelbach A, Liu L, Zierold U, Altschmied L, Maucher H, Beier F, et al. Promotors of the barley germin-Like GER4 gene cluster enable strong transgene expression in response to pathogen attack. Plant Cell 2010; 22: 937–952. pmid:20305123
- 47. Chomczynski P. A reagent for the single-step simultaneous isolation of RNA, DNA and proteins from cell and tissue samples. BioTechniques 1993; 15: 532–537. pmid:7692896
- 48. Yang L, Ding J, Zhang C, Jia J, Weng H, Liu W, et al. Estimating the copy number of transgenes in transformed rice by real-time quantitative PCR. Plant Cell Rep 2005; 23: 759–763. pmid:15459795