- Research note
The remarkable plethora of infestation-responsive Q-type C2H2 transcription factors in potato
BMC Research Notesvolume 11, Article number: 398 (2018)
Q-type C2H2 transcription factors (TF) play crucial roles in the plant response to stress, often leading to regulation of downstream genes required for tolerance to these challenges. An infestation-responsive Q-type C2H2 TF (StZFP2) is induced by wounding and infestation in potato. While mining the Solanum tuberosum Group Phureja genome for additional members of this family of proteins, five StZFP2-like genes were found on a portion of chromosome 11. The objective of this work was to differentiate these genes in tissue specificity and expression upon infestation.
Examination of different tissues showed that young roots had the highest amounts of transcripts for five of the genes. Expression of their transcripts upon excision or infestation by Manduca sexta, showed that all six genes were induced. Overall, each gene showed variations in its response to infestation and specificity for tissue expression. The six genes encode very similar proteins but most likely play unique roles in the plant response to infestation. In contrast, only two homologs have been identified in Arabidopsis and tomato. Overexpression of similar genes has led to enhanced tolerance to, for example, salinity, drought and pathogen stress. Discovery of these new StZFP2 homologs could provide additional resources for potato breeders.
Plants have an array of responses to environmental stresses such as excessive salinity, drought, pathogens and insects. Stressors are rapidly sensed by the plant, triggering hormonally mediated reactions . Generally, attack by chewing insects lead to the induction of the jasmonic acid (JA) pathway and an array of direct and indirect defense compounds that presumably aid in protecting the plant under attack. Identification of an infestation responsive Q-type C2H2 transcription factors (TF), in potato, StZFP2  led to our interest in this family of TFs. Overexpression of these TFs can result in increased tolerance to stress [3,4,5,6]. C2H2 zinc finger proteins (ZFP) contain two cysteine and two histidine residues that surround a zinc ion. While Englbrecht et al. described a total of 176 C2H2 ZFPs in Arabidopsis, 64 are found in the C1 family containing either one zinc finger domain (ZFD) or two to five dispersed ZFDs [7, 8]. The C1 group mostly contains the QALGGH motif, unique to plants and a crucial part of the 30 amino acid ZFD. Many of these proteins are induced by stress [8, 9]. Q-type is short hand for the QALGGH motif. Arabidopsis has 18 C1-2i proteins that contain two Q-type ZFDs, with the ZFDs acting as DNA-binding domains [8, 9]. Another important amino acid motif is the EAR domain, which is as an active repressor with either a DLN or LxLxL profile. The EAR domain was named from its initial discovery in the ethylene-responsive element binding factor or ERF family proteins. EAR stands for ERF-associated amphiphilic repression motif [10, 11]. The Q-type C2H2 ZFPs also contain two additional motifs, an L-box, rich with leucine residues and the B-box, containing a nuclear localization signal (NLS) required for transcriptional regulation.
After identifying StZFP2 as infestation- and wound-responsive , we used the potato genome [12, 13] to discover additional StZFP genes. Six additional StZFP genes map within 28,000 bases to each other on chromosome (chr) 11 in potato. In comparison, Arabidopsis and tomato only have three such orthologs. Tomato and potato both have 12 chromosomes and the three orthologs are also in the same region on chr 11 in tomato. It is estimated that potato and tomato formed separate species ~ 7.3 million years ago , so differences within the individual genomes would have evolved later. This work focuses on the five new StZFP2-like genes in potato, demonstrating that StZFP2-like transcripts are expressed, differentiating their expression in different tissues and upon infestation by the chewing insect pest Manduca sexta.
RNA was isolated using TRIzol per the manufacturer’s instructions . After a final ethanol wash, Qiagen DNase and RNeasy columns were used to further purify the RNA. Integrity and concentration of RNA was determined with an Experion automated electrophoresis system for RNA (Bio-Rad).
Since infestation can perturb typical housekeeping genes  an exogenous control RNA was added to the cDNA synthesis step in lieu of an endogenous control. This approach has been used in several situations when no suitable reference gene can be found [17,18,19]. Luciferase control RNA (Promega) was added to the reverse transcription master mix. The cDNA synthesis reaction contained 2 µg template RNA and 2 ng luciferase RNA in 20 µl. This was diluted 1:20 with 10 µl used in the realtime reaction (50 ng RNA, and 50 pg luciferase RNA). Custom Taqman assay primer/probe sequences are listed (Additional file 1: Table S1). The qRT-PCR was performed using the 7500 Real Time PCR System (Life Technologies).
For the phylogenetic tree, proteins were aligned using CLUSTSALW and a tree was inferred using the Neighbor joining method . The evolutionary distances were computed using the Poisson correction method  and are in the units of the number of amino acid substitutions per site 1000 bootstrap replicates were performed . Evolutionary analyses were conducted in MEGA6 . Alignment of ZFP proteins was performed using CLUSTAL OMEGA .
The plant material, Solanum tuberosum var. Kennebec (a common commercial variety) and insects were prepared as follows. At 2 weeks, tissue culture maintained rooted nodal plantlets were transferred to 3.5″ Arabipots (Lehle Seeds) containing Sunshine LC1 mix (Sun Gro Horticulture, Inc) for 19 days in a Conviron CMP6050 growth chamber with a 16:8 light/dark and 25 °C day/20 °C night temperature. The two youngest fully expanded leaves were excised and used for ‘control’ and ‘infested’ in feeding experiments. Manduca sexta larvae were reared to 4th instar on artificial diet  then transferred to Kennebec “feeder” plants in Arabipots for 18 h. Larvae were removed from plant material to empty petri dishes 2 h before assay.
Three-way ANOVA was conducted on the ΔCt values  using SAS Proc MIXED  to accurately model among-genes and between-treatments. Pairwise means comparisons were conducted using the SLICE option of the LSMEANS statement to compare among times for each gene and treatment and to compare treatments for each gene at each time; using the Sidak method to adjust p values to protect against obtaining false positive comparisons. Letters indicating significant differences among means were generated using the pdmix800 SAS macro .
Results and discussion
Using the StZFP2 protein sequence as a query, Blastp was used to search the Spud database . The search identified the StZFP2, and six additional proteins. The names in the Spud DB database , and their map positions, are detailed in (Additional file 2: Table S2). StZFP2 gene maps to Chr11. The area around StZFP2 in the potato genome contains a cluster of genes that encode similar proteins (Fig. 1a). The ZFPs previously identified were StZFP1 and StZFP2 [2, 31]. While StZFP1 maps to chr 01 , the six potato StZFP genes form a 28 kb cluster on chr 11 (Fig. 1a). Two StZFP2-like orthologs are also found on chr 11 of the tomato genome and named Solyc11gZFP1 and Solyc11gZFP2. Their map position and database names are listed in (Additional file 3: Table S3). An additional StZFP gene (StZFP8) was found within the cluster on chr 11 and although this gene is not in the tomato database it was identified here, between the ZFP2-like genes and named Solyc11gZFP3 (Fig. 1a).
A phylogenetic tree (Fig. 1b) illustrates how the potato ZFPs compare to a family of similar proteins in Arabidopsis. The additional ZFPs in this work were named StZFP3 to StZFP8. The eight potato ZFPs, three tomato ZFPs residing on chr 11 and 18 Arabidopsis two fingered Q-type ZFPs (Zats) were aligned as described in “Methods”. Due to their conserved protein sequences, the 18 Arabidopsis ZFPs are considered C1-2i types, since they contain two zinc finger domains clustering into 2iA–D with an outlier X. Six of the 7 StZFPs, Solyc11gZFP1 and two clustered with the Arabidopsis Zat11 and Zat18 (2i-A), while StZFP8 grouped with tomato Solyc11gZFP3 and 2i-C Zats. StZFP1 clustered with the 2i-D Zats (Fig. 1b).
There were only two corresponding 2i-A genes in tomato (Fig. 1a) and Arabidopsis (Fig. 1b). StZFP6 and StZFP7 are orthologs of the tomato gene Solyc11gZFP1, while StZFP2, StZFP3, StZFP4 and StZFP5 are orthologs of a gene Solyc11gZFP2, which is not currently annotated on the tomato genome. The alignment of the protein sequences (Fig. 1c) confirms the phylogenetic tree, with the tomato protein Solyc11gZFP2 aligning with StZFP2-5 and the second tomato protein Solyc11gZFP1 aligning with StZFP6-7. Comparison of the protein sequence to tomato proteins showed that StZFP7 is closest to Solyc11gZFP1 while StZFP4 is closest to Solyc11gZFP2, suggesting that StZFP7 and StZFP4 existed prior to the split between the two species. The L-box, two zinc finger domains, and the EAR motif are also delineated (Fig. 1c). The EAR motif has been identified as an active repressor of genes that are bound by the ZFP zinc finger motifs  and the most common types are LxLxL and DLNxxP. The StZFP2-5 proteins all contain LDLNLP, which is a combination of these two classic EAR motifs.
Considering that there are StZFP2 and five StZFP2-like genes (2i-A) in potato, it is curious that there are only two in tomato and Arabidopsis (2i-A). The sequence read archive (SRA)  data for potato was queried and confirmed that all five transcripts were expressed. The most robust expression of these genes was in a sample of stolon tips rewatered after drought stress and, not surprisingly, the variety Igor after 24 h infestation by a chewing insect (Additional file 4: Figure S1). Thus, according to SRA data these genes are all expressed under stress conditions.
qRT-PCR was used to determine if these genes were expressed in Kennebec. Primers identified the unique transcripts for StZFP2-7, with StMYC2 and StLOX3 marker genes for JA (Additional file 1: Table S1). StZFP2 and the five StZFP2-like genes were compared (Fig. 2). The least amount of StZFP expression was in flower and the greatest in young roots. StZFP3 was an exception with very low levels in all tissues except young stem and leaves. StZFP5 and 7 had similar profiles of tissue expression. StZFP6 had the highest expression in young root and the lowest in the leaf. In mature root the expression was similar for all StZFPs. Since basal levels in the root are higher for all but StZFP3, examination of root tissue may lead to a better understanding of the role the StZFP2-like genes play in potato.
qRT-PCR was used to examine whether the StZFP2-like genes were expressed upon excision and infestation by the chewing insect pest M. sexta. Excised leaves were subjected to infestation for 20 min, the larvae removed and the leaves harvested at the times indicated. Uninfested excised leaves were compared to infested at each time point (Fig. 3). The expression in detached leaves diminished earlier at 160 min while the infested leaves were not reduced until the final time point. StLOX3 and StMYC2 expression diminished by 160 min. StZFP2, StZFP3 and StZFP4 were significantly reduced in expression by 80 min in both detached and infested leaves. The expression of StZFP6 in infested leaves was significantly lower at 160 min, while replicates in the detached leaves were so variable that significant differences were not found. The expression of StZFP5 and StZFP7 were significantly higher in infested leaves versus detached at 40 and 80 min, with StZFP7 still higher in the infested leaves at 160 min. Therefore, while expression of StZP5 and StZFP7 were affected by detachment, infestation significantly increased their expression.
While StZFP2 and StZFP3 were similarly expressed by infestation, they had differences in tissue specificity. StZFP2 has ~ 50% of the young tissue expressed in the root and StZFP3 only ~ 18%. StZFP4 and StZFP6 are similarly expressed by infestation, but StZFP6 had twice as much expression in young root compared with StZFP4. StZFP5 and StZFP7 had similar tissue specific expression, but were unique in their response to infestation with StZFP5 falling at 160 min, and StZFP7 remaining significantly high at 160 min. Sequence similarity suggests that StZFP6 and StZFP7 were closest homologs, while StZFP2, StZFP3 and StZFP5 were derived from StZFP4. The expression of these genes does not appear to fit into these lineages. For example, the amount of expression in the roots for StZFP6 is high while StZFP7 is intermediate. The expression of StZFP6 by infestation decreases by 80 min while StZFP7 remains high until at least 160 min.
Initially, StZFP2 was identified as a gene induced by insect infestation . Clearly, it is not the only ZFP in potato responsive to infestation. It is also not the most dramatically induced of the StZFP2-like genes. Why there are so many StZFP2-like genes in potato compared to Arabidopsis and tomato remains to be determined. It is predicted that the tomato and the potato genome contain approximately the same number of genes with 34,727 and 35,004 protein coding genes, respectively , so it is unlikely that the propagation of potatoes as tubers rather than through seed has simply allowed redundant genes to be retained in potato.
The StZFP2-like genes identified here are infestation induced and differ in their tissue specificity. Now that they have been identified, further research is needed to determine how the individual StZFPs affect the plants response to infestation. Perhaps over-expression of individual StZFP2-like genes may create different phenotypes with some more effective than others in improving resistance to insect pests with less negative effects on growth and development.
ethylene-responsive element-binding factor
ERF-associated amphiphilic repression
nuclear localization signal
sequence read archive
Erb M, Meldau S, Howe GA. Role of phytohormones in insect-specific plant reactions. Trends Plant Sci. 2012;17:250–9.
Lawrence SD, Novak NG, Jones RW, Farrar RR, Blackburn MB. Herbivory responsive C2H2 zinc finger transcription factor protein StZFP2 from potato. Plant Physiol Biochem. 2014;80:226–33.
Ciftci-Yilmaz S, Morsy MR, Song L, Coutu A, Krizek BA, Lewis MW, Warren D, Cushman J, Connolly EL, Mittler R. The EAR-motif of the Cys2/His2-type zinc finger protein Zat7 plays a key role in the defense response of Arabidopsis to salinity stress. J Biol Chem. 2007;282:9260–8.
Mittler R, Kim Y, Song L, Coutu J, Coutu A, Ciftci-Yilmaz S, Lee H, Stevenson B, Zhu JK. Gain- and loss-of-function mutations in Zat10 enhance the tolerance of plants to abiotic stress. FEBS Lett. 2006;580(28–29):6537–42.
Shi H, Wang X, Ye T, Chen F, Deng J, Yang P, Zhang Y, Chan Z. The Cysteine2/Histidine2-Type Transcription Factor ZINC FINGER OF ARABIDOPSIS THALIANA6 modulates biotic and abiotic stress responses by activating salicylic acid-related genes and C-REPEAT-BINDING FACTOR genes in Arabidopsis. Plant Physiol. 2014;165(3):1367–79.
Yin M, Wang Y, Zhang L, Li J, Quan W, Yang L, Wang Q, Chan Z. The Arabidopsis Cys2/His2 zinc finger transcription factor ZAT18 is a positive regulator of plant tolerance to drought stress. J Exp Bot. 2017;68(11):2991–3005.
Englbrecht CC, Schoof H, Bohm S. Conservation, diversification and expansion of C2H2 zinc finger proteins in the Arabidopsis thaliana genome. BMC Genomics. 2004;5:39.
Ciftci-Yilmaz S, Mittler R. The zinc finger network of plants. Cell Mol Life Sci. 2008;65:1150–60.
Kielbowicz-Matuk A. Involvement of plant C2H2-type zinc finger transcription factors in stress responses. Plant Sci. 2012;185:78–85.
Ohta M, Matsui K, Hiratsu K, Shinshi H, Ohme-Takagi M. Repression domains of class II ERF transcriptional repressors share an essential motif for active repression. Plant Cell. 2001;13:1959–68.
Sakamoto H, Maruyama K, Sakuma Y, Meshi T, Iwabuchi M, Shinozaki K, Yamaguchi-Shinozaki K. Arabidopsis Cys2/His2-type zinc-finger proteins function as transcription repressors under drought, cold, and high-salinity stress conditions. Plant Physiol. 2004;136:2734–46.
Potato Genome Consortium. Genome sequence and analysis of the tuber crop potato. Nature. 2011;475(7355):189–95. https://doi.org/10.1038/nature10158.
Hirsch CD, Hamilton JP, Childs KL, Cepela J, Crisovan E, Vaillancourt B, Hirsch CN, Habermann M, Neal B, Buell CR. Spud DB: a resource for mining sequences, genotypes, and phenotypes to accelerate potato breeding. Plant Genome. 2014. https://doi.org/10.3835/plantgenome2013.12.0042.
Wu F, Tanksley SD. Chromosomal evolution of the plant family Solanaceae. BMC Genomics. 2010;11:182.
Trizol user guide. https://tools.thermofisher.com/content/sfs/manuals/trizol_reagent.pdf. Accessed 17 June 2018.
Rehrig EM, Appel HM, Schultz JC. Measuring ‘normalcy’ in plant gene expression after herbivore attack. Mol Ecol Resour. 2011;11:294–304. https://doi.org/10.1111/j.1755-0998.2010.02929.x.
Ellefsen S, Stensløkken KO, Sandvik GK, Kristensen TA, Nilsson GE. Improved normalization of real-time reverse transcriptase polymerase chain reaction data using an external RNA control. Anal Biochem. 2008;376:83–93. https://doi.org/10.1016/j.ab.2008.01.028.
Johnston S, Gallaher Z, Czaja K. Exogenous reference gene normalization for real-time reverse transcription-polymerase chain reaction analysis under dynamic endogenous transcription. Neural Regen Res. 2012;7(14):1064–72. https://doi.org/10.3969/j.issn.1673-5374.2012.14.004.
Ortiz SC, Trienens M, Rohlfs M. Induced fungal resistance to insect grazing: reciprocal fitness consequences and fungal gene expression in the Drosophila Aspergillus model system. PLoS ONE. 2013;8(8):e74951. https://doi.org/10.1371/journal.pone.0074951.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Zuckerkandl E, Pauling L. Evolutionary divergence and convergence in proteins. In: Bryson V, Vogel HJ, editors. Evolving genes and proteins. New York: Academic Press; 1965. p. 97–166.
Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39:783–91.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
CLUSTAL OMEGA program. https://www.ebi.ac.uk/Tools/msa/clustalo/. Accessed 17 June 2018.
Martin PAW, Michael B, Blackburn MB. Using combinatorics to screen Bacillus thuringiensis isolates for toxicity against Manduca sexta and Plutella xylostella. Biol Control. 2007;42:226–32.
Khan-Malek R, Wang Y. Statistical analysis of quantitative RT-PCR results. In: Gautier JC, editor. Drug safety evaluation: methods and protocols, methods in molecular biology, vol. 691. Berlin: Springer Science + Business Media, LLC 20; 2011. p. 227–41. https://doi.org/10.1007/978-1-60761-849-2_13.
SAS v9.4.2016;SAS/STAT, PROC MIXED, SAS Institute, Inc., Cary, NC.
Saxton AM. A macro for converting mean separation output to letter groupings in Proc Mixed. In: Proc. 23rd SAS Users Group Intl., SAS Institute, Cary, NC. 1998. p. 1243–1246.
Spud database. http://solanaceae.plantbiology.msu.edu/integrated_searches.shtml. Accessed 17 June 2018.
Spud DB database. http://potato.plantbiology.msu.edu/. Accessed 17 June 2018.
Tian ZD, Zhang Y, Liu J, Xie CH. Novel potato C2H2-type zinc finger protein gene, StZFP1, which responds to biotic and abiotic stress, plays a role in salt tolerance. Plant Biol. 2010;12(5):689–97. https://doi.org/10.1111/j.1438-8677.2009.00276.x.
Kagale S, Rozwadowski K. EAR motif-mediated transcriptional repression in plants. An underlying mechanism for epigenetic regulation of gene expression. Epigenetics. 2011;6(2):141–6. https://doi.org/10.4161/epi.6.2.13627.
SRA data. https://www.ncbi.nlm.nih.gov/sra. Accessed 17 June 2018.
Green JM, Appel H, Rehrig EM, Harnsomburana J, Chang JF, Balint-Kurti P, Shyu CR. PhenoPhyte: a flexible affordable method to quantify 2D phenotypes from imagery. Plant Methods. 2012;8(1):45. https://doi.org/10.1186/1746-4811-8-45.
Lawrence SD, Novak NG, Ju CJ, Cooke JE. Potato, Solanum tuberosum, defense against Colorado potato beetle, Leptinotarsa decemlineata (Say): microarray gene expression profiling of potato by Colorado potato beetle regurgitant treatment of wounded leaves. J Chem Ecol. 2008;34(8):1013–25.
Tomato Genome Consortium. The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012;485(7400):635–41.
NN: performed all aspects of tissue and infestation experiments and produced all final figures. SL: preformed all data mining, sequence analysis, phylogenetics and wrote the manuscript. Both authors read and approved the final manuscript.
The authors wish to thank Dr. Bryan Vinyard for the statistical analysis of the infestation data. The authors appreciate the feedback of Drs. Donald Weber and Ann Smigoki on an early version of the manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analysed during this study are included in this published article (and its Additional files).
Consent for publication
Ethics approval and consent to participate
This work was funding by USDA-Agricultural Research Service 8042-22000-288-00-D.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.