precise regulation of gene expression in response to environmental and endogenous signals is fundamental for reproductive success .
the ability to adapt gene expression in response to the changing external conditions is particularly important for plants owing to their sessile nature . in eukaryotes ,
the organization of genomic dna into chromatin and the ordered regulation of its accessibility to the transcription machinery are central to gene regulation ( li et  al .
mechanisms for regulation of chromatin structure include atp - dependent chromatin remodeling as well as post - translational histone modifications .
atp - dependent chromatin remodeling complexes alter nucleosome composition  and positioning , and thus can regulate dna accessibility via chromatin compactness ( li et  al . , 2007 , narlikar et  al . , 2013 ) .
the swr1 complex ( swr1c ) is an evolutionarily conserved swi2/snf2-related atpase - containing chromatin remodeling complex that catalyzes the replacement of h2a by the histone variant h2a.z in nucleosomes ( kobor et  al . , 2004 ,
2004 ) . through its distinct physicochemical properties , h2a.z influences nucleosome stability , and therefore chromatin structure , to modulate gene expression ( thambirajah et  al . , 2006 ,
these properties along with its incorporation into the chromatin out of mitosis have made h2a.z central to transcriptional regulation underlying development and environmental responses ( talbert and henikoff , 2014 ) . in budding yeast ,
swr1c has been characterized as a large multi - protein structure with 13 accessory subunits in complex with the atpase , swr1 ( kobor et  al . , 2004 ,
several of the non - catalytic subunits including actin - related protein 6 ( arp6 ) and swr1 complex 6 ( swc6 ) have been shown to be essential for histone replacement .
arp6 , swc6 , and swc2 act as a sub - complex , where the proteins are mutually essential for each other 's association and function within the complex ( wu et  al . , 2005 ) .
in line with the current understanding , swr1c is expected to act as a single complex in plants to catalyze h2a.z incorporation ( figure  1a ) .
the components and function of swr1c have been shown to be conserved in arabidopsis ( noh and amasino , 2003 , deal et  al . , 2005 , martin - trillo et  al . , 2006 , choi et  al . , 2007 , march - diaz et  al . , 2007 , lzaro et  al . ,
in contrast to other organisms , swr1c and h2a.z are not essential for  viability in arabidopsis ( coleman - derr and zilberman , 2012 ) .  forward and
reverse genetic analyses have identified the  swr1 homolog  photoperiod - independent early flowering 1 ( pie1 ) as well as arp6 , swc6 , and h2a.z encoding genes ( hta8 , hta9 , hta11 ) ( noh and amasino , 2003 , deal et  al . , 2005 , martin - trillo et  al . , 2006 , lzaro et  al .
, 2008 , march - daz et  al . , 2008 , coleman - derr and zilberman , 2012 ) .
as in other eukaryotes , the non - catalytic subunits of swr1c in arabidopsis are considered to have similar functions and to be essential for h2a.z deposition by the atpase subunit .
mutants defective in swr1c are highly pleiotropic ( march - daz et  al . , 2008 , coleman - derr and zilberman , 2012 ) .
swr1c components and h2a.z have been shown to function in various aspects of plant development , most importantly in the timing of the vegetative to reproductive transition ( noh and amasino , 2003 , deal et  al . , 2005 , martin - trillo et  al . , 2006 ,
arp6 has been shown to regulate meiosis during megagametogenesis through the regulation of meiotic gene expression ( qin et  al . , 2014 ) .
in addition , the role of h2a.z in modulating chromatin structure was recently shown to be also important for crossovers during meiosis .
while the chromatin landscape at the crossover hotspots is marked with h2a.z nucleosomes , the arp6 mutant is compromised in meiotic crossovers ( choi et  al . , 2013 ) .
h2a.z has recently been proposed to be critical for genome stability and dna repair in arabidopsis .
the pie1 , arp6 and swc6 mutants show constitutive dna damage and compromised somatic homologous recombination as well as hypersensitivity to genotoxic agents ( rosa et  al . , 2013 ) .
importantly , swr1c and h2a.z have also been implicated in regulating gene expression in response to environmental signals in general ( coleman - derr and zilberman , 2012 ) .
arp6 has been shown to be important for repression of genes involved in the phosphate starvation response ( smith et  al . , 2009 ) .
temperature - dependent h2a.z nucleosome dynamics have been shown to modulate thermosensory responses in arabidopsis ( kumar and wigge , 2010 ) and in grasses ( boden et  al . , 2013 ) .
mutations in swr1c components pie1 and swc6 or in h2a.z have been reported to result in constitutive activation of defense responses , demonstrating their importance in biotic interactions as well ( march - daz et  al . , 2008 ) .
although swr1c and h2a.z have been implicated in immunity in arabidopsis , the respective role of the complex subunits is still not well understood in plant defense .
moreover the function of arp6 , a key component of the swr1c subunit in h2a.z deposition , has not been previously studied in this respect . in this study , we have carried out a comprehensive analysis to unravel the role of swr1c components in immunity and gene regulation .
while resistance to biotrophic and necrotrophic pathogens is impaired in pie1 , swc6 , and hta9 hta11 mutants , arp6 showed wild - type or increased resistance .
loss of pie1 and swc6 but not arp6 , leads to impaired effector - triggered immunity ( eti ) .
genome - wide gene expression analyses have further highlighted the potentially specialized roles of h2a.z and swr1c components in gene regulation and thereby in plant defense responses .
to study the role of h2a.z in plant defense processes , we analyzed mutants defective in swr1 complex components pie1 , arp6 , and swc6 , as well as those depleted of the histone variant h2a.z .
consistent with previous studies , pie1 - 2 ( noh and amasino , 2003 ) , arp6 - 1 ( deal et  al . , 2005 ) , and swc6 - 1 ( lzaro et  al . , 2008 )
mutants showed pleiotropic growth and developmental phenotypes including early flowering , elongated petioles , and reduced fertility when grown in long photoperiod ( ld ) ( figure  1b ) .
the pie1 mutant showed the most severe growth defects and even though it displayed accelerated reproductive transition , bolting was delayed .
the arp6 and swc6 mutants phenocopied each other with characteristic early flowering and serrated leaves .
this was consistent with the biochemical interaction of arp6 and swc6 in yeast , where their existence in the complex is mutually dependent ( mizuguchi et  al . , 2004 ) .
the hta9 - 1 hta11 - 1 double mutant ( loss of two major h2a.z encoding genes out of the possible three in arabidopsis ) phenocopied arp6 and swc6 mutants with early flowering and serrated leaves .
the severe growth defects in the pie1 mutant have been previously attributed to the de - repression of immune responses characterized by spontaneous cell death and upregulation of defense genes ( march - daz et  al . , 2008 ) .
growth in short photoperiods ( sd ) largely suppressed some of the growth defects of the pie1 mutant , even though the rosette size remained significantly smaller compared with the other swr1c mutants ( figure  1b ) .
all the studies reported here , unless otherwise specified , were performed on plants grown under sd .
swr1c components have been proposed to be negative regulators of plant immunity with pie1 , swc6 , and hta9 hta11 mutants reported to display spontaneous cell death and enhanced resistance to virulent bacterial pathogens ( march - daz et  al . , 2008 ) .
the role of arp6 in this process is still not known and has been expected to be similar to pie1 , swc6 , and h2a.z . in order to investigate this , we characterized the defense responses of these mutants . in our experimental conditions ,
trypan blue staining revealed no apparent spontaneous cell death in the mutants ( figure  1c ) .
we occasionally observed patches of cell death and leaf torsion in some leaves of the pie1 mutant ( supplemental figure  1 ) .
to investigate their immunity phenotypes further , we challenged 4-week - old mutants with the hemibiotrophic bacterial pathogen p.  syringae pv tomato dc3000 ( pst dc3000 ) by spray inoculation and monitored bacterial growth .
we found that pie1 , swc6 , and hta9 hta11 mutants showed more macroscopic disease symptoms and increased susceptibility toward  pst dc3000 compared with the wild - type col-0 .
however , the arp6 mutant showed increased resistance ( figure  1d and 1e ) .
overall , pie1 , swc6 , and hta9 hta11 mutants accommodated 10 to 15 times more bacterial titers at 3  days post inoculation ( dpi ) than the wild type .
the pie1 , swc6 , and hta9 hta11 mutants displayed severely compromised resistance similar to that of the enhanced disease susceptibility1 ( eds1 - 2 ) mutant , which was used as a susceptible control ( feys et  al . , 2001 ) .
there were no differences in the bacterial titers at day 0 ( 2  h post inoculation ) between the different mutants and the wild - type col-0 plants .
a transgenic swc6 complemented line showed wild - type resistance to pst dc3000 ( supplemental figure  2 ) , suggesting that the resistance phenotype we observe is caused by the mutation .
our observations suggest that h2a.z , pie1 , and swc6 are essential for basal resistance in arabidopsis , whereas arp6 has an opposite function .
these results were intriguing as a previous study has shown that pie1 , swc6 , and hta9 hta11 mutants displayed rather constitutive defense activation marked by spontaneous lesions and disease - resistance phenotypes ( march - daz et  al . , 2008 ) .
as mentioned earlier , the pie1 mutant showed severe growth defects when grown under ld .
these phenotypes were partially suppressed by growth under sd ( figure  1b ) . in order to investigate
whether the severe growth defects in pie1 under ld were due to constitutive defense activation and therefore to test if pie1 shows photoperiod - dependent defense phenotypes , we assessed its resistance to pst dc3000 in these growth conditions . to bypass early flowering induction by constant growth under ld photoperiods , the mutants along with the wild type were pre - grown for two weeks under sd conditions before being shifted to ld for another two weeks prior to infection .
consistent with the results from the sd experiments , pie1 remained increasingly susceptible to pst dc3000 ( supplemental figure  3a ) .
moreover , it showed no spontaneous cell death in the leaves under these conditions , similarly to the other mutants ( supplemental figure  3b ) , except for the patches of dead tissue observed as in sd .
these results suggest that the growth defects in the pie1 mutant are not necessarily a consequence of enhanced defense activation .
the defense phenotypes of arp6 , swc6 , and hta9 hta11 grown in ld were also similar to the results obtained under sd conditions ( supplemental figure  3a ) further confirming that the immunity phenotypes we observed are not due to specific growth conditions , especially the photoperiod .
these results also confirm that the contrasting defense phenotypes observed in this study and by march - daz et  al .
( 2008 ) are unlikely to be caused by growth conditions such as the photoperiod . during interaction with avirulent pathogens , recognition of effectors by the receptor proteins in the host plant leads to the activation of eti . in arabidopsis , the bacterial effectors avrrps4 and avrrpt2 trigger resistance through recognition by the toll interleukin-1 receptor ( tir)-type and the coiled coil ( cc)-type nucleotide binding - leucine - rich - repeat ( nb - lrr ) proteins rps4/rrs1 and rps2 , respectively ( kunkel et  al .
2009 ) . to examine the possible roles of swr1c and h2a.z in eti , we inoculated the corresponding mutants with avirulent pst dc3000 expressing either avrrps4 or avrrpt2 effectors .
the pie1 and swc6 mutants are more susceptible compared with wild - type in both interactions ( figure  2a and 2b ) .
they were , however , not as hypersusceptible as eds1 or ndr1 ( non - race - specific disease resistance1 ) , used as susceptible controls in response to avrrps4 and avrrpt2 ( century et  al . , 1995 ) , respectively .
consistent with the earlier results , arp6 displayed full resistance like wild - type plants , suggesting that arp6 mutation does not alter rps4- and rps2-mediated resistance .
the hta9 hta11 double mutant showed mild but not significant susceptibility compared with wild - type . in arabidopsis , h2a.z is encoded by three genes hta9 , hta11 , and hta8 .
the presence of hta8 , even though generally lowly expressed ( supplemental figure  4 ) , in the hta9 hta11 double mutant could compensate for the lack of hta9 and hta11 .
interestingly , arp6 function appears to be dispensable . to understand the basis of the failure to mount eti
, we examined the level of salicylic acid ( sa ) accumulation , an important defense hormone in response to biotrophic pathogens and particularly during eti ( robert - seilaniantz et  al . , 2011 ) .
the basal level of sa did not noticeably differ between the the wild type and the mutants except for arp6 , where it was higher ( figure  2c ) .
interestingly , 24  h after infiltration with pst dc3000 avrrps4 , all  the mutants showed increased sa accumulation in comparison with col-0 .
accumulation of sa in pie1 , swc6 , and hta9 hta11 was 2-fold higher relative to col-0 , and reached 3-fold in arp6 .
our results show that in spite of sa accumulation , resistance to virulent and avirulent p.  syringae is impaired in pie1 and swc6 mutants , suggesting that swr1c activity ( especially pie1 , swc6 ) and h2a.z are required for sa - dependent downstream signaling processes .
these results are further supportive to the resistance phenotypes of the mutants to virulent pst dc3000 ( figure  1d and 1e ) .
signaling crosstalks and antagonistic interactions between different defense pathways regulated by hormones have been extensively studied ( robert - seilaniantz et  al . , 2011 ,
sa - mediated defense is activated during interaction with biotrophic pathogens while jasmonic acid ( ja)-mediated defense is active against necrotrophs and herbivores ( glazebrook , 2005 , howe and jander , 2008 ) .
as presented earlier , pie1 , swc6 , and hta9 hta11 mutant plants showed impaired sa - mediated defense responses even though they accumulated a high level of the sa hormone in response to pseudomonas infection . in order to investigate
whether these mutants are affected in the sa / ja crosstalk , we investigated their resistance to the necrotrophic pathogen botrytis cineria .
leaves of the wild - type and mutant plants were drop inoculated with spore suspension , and the growth of lesions was monitored .
average lesion size measured 3 dpi was significantly bigger in pie1 , swc6 , and hta9 hat11 .
the proportions of outgrowing lesion type iv ( lesion over 6 mm in diameter ) in pie1 , swc6 , and hta9 hat11 were higher compared with arp6 mutant plants that showed mild resistance ( figure  3 and supplemental figure  5 ) .
these results show that perturbation of h2a.z incorporation results in a compromised ja / ethylene ( et)-mediated immunity .
it also suggests that components of swr1c do not play the same function in this defense pathway , as arp6 showed wild - type resistance to botrytis , as opposed to the rest of the mutants . in order to elucidate the molecular basis of the phenotypes described above , and to dissect the role of h2a.z and swr1c subunits in gene regulation
, we carried out a global transcriptome analysis of the wild - type and mutant lines using rna - seq .
three replicates of total rna from 2-week - old seedlings from each mutant were sequenced .
genes showing at least 2-fold change from the wild type with a p value cut - off of 0.05 were used for further analyses .
we found that mutations in swr1c components result in misregulation of a large number of genes ( figure  4a and 4b and supplemental dataset 1 ) .
the pie1 mutation resulted in the largest number of gene misregulation , with 2295 genes upregulated and 1051 genes downregulated .
the genes misregulated in our study showed statistically significant ( hypergeometric test ) overlap with the previously published data ( march - daz et  al . , 2008 , coleman - derr and zilberman , 2012 ) .
we found a total of 671 genes to be commonly misregulated among all four mutants ( 473 genes upregulated and 198 genes downregulated ) .
this accounts for 20% of the genes misregulated in pie1 , and 22% , 32% , and 27% in arp6 , swc6 , and hta9 hta11 , respectively .
moreover , the commonly misregulated genes ( up and down ) showed a statistically significant overlap with h2a.z enrichment in the gene  bodies ( coleman - derr and zilberman , 2012 ) .
we found that 1262 genes were uniquely upregulated in pie1 accounting for 55% of the total number of genes upregulated in this mutant . in arp6 , swc6 , and hta9 hta11 , however ,
the proportions of uniquely downregulated genes were 25% in arp6 , 4.5% in swc6 , and 18% in hta9 hta11 .
the large proportion of uniquely misregulated genes in pie1 suggests a potential swr1c / h2a.z independent role for pie1 in gene regulation . to explore this further and to elucidate the relationship between individual mutants
, we performed principal component analysis ( pca ) of the rna - seq data .
the first three principal components ( pc1 , pc2 , and pc3 ) together accounted for 92.1% of the variance in the data ( 44.44% , 30.55% , and 17.11% , respectively ) ( figure  4c ) .
it also revealed that the four mutant genotypes segregate into three distinct groups : swc6 and hta9 hta11  mutants co - segregated while pie1 and arp6 mutants were  distinct .
pc1 clearly separated pie1 from the rest of the group , whereas pc3 defined the contribution of arp6 ( figure  4c and supplemental figure  6 ) .
this analysis has further strengthened the possibility of non - overlapping functions for the complex components .
in line with  the pca analysis , differential gene expression analyses showed a poor correlation of gene misregulation caused by arp6 and pie1 mutations ( r  = 0.16 ) ( figure  4d ) , as well as between hta9 hta11 and pie1 ( r  =  0.11 ) ( figure  4e ) , which is surprising as pie1 mutation is expected to phenocopy h2a.z depletion to a large extent .
the arp6 and swc6 mutants in arabidopsis have similar morphological and developmental phenotypes , especially early flowering and serrated leaves .
moreover , arp6 and swc6 have been shown to physically interact in arabidopsis ( march - diaz et  al .
intriguingly , in our rna - seq data , the transcriptional changes due to these mutations did not show a strong correlation ( r  = 0.42 ) ( figure  4f ) .
this observation is consistent with the defense phenotypes described above , suggesting that arp6 and swc6 could have distinct functions in arabidopsis . in agreement with the pca analysis ,
swc6 and hta9 hta11 double mutants displayed a strong correlation of transcriptional misregulation ( r  = 0.75 ) , indicating overlapping roles for swc6 and h2a.z in gene regulation ( figure  4 g ) .
this is in contrast to arp6 , which showed only a modest correlation ( r  = 0.38 ) with hta9 hta11 ( supplemental figure  7 ) . to understand the function of swr1c and h2a.z in gene regulation , particularly in immunity
we found a significant over - representation  of sa response ( p  = 4.46   10 ) and defense ( p  = 7.75   10 ) related go terms in the 473 genes commonly misregulated in the mutants ( supplemental dataset  2 ) .
this is in accordance with what was reported earlier ( march - daz et  al . ,
the commonly downregulated genes where enriched in triterpene ( thalianol , p  = 1.77   10 and tricyclic triterpenoid metabolic process , p  = 7.04   10 ) related go terms .
in addition , pie1 showed a significant over - representation of several biological processes including auxin metabolic processes , membrane organization , and response to heat , indicating the important role of pie1 in gene expression in general . in swc6 and hta9
hta11 , the go term response to heat was significantly enriched , which is consistent with h2a.z function in response to temperature ( kumar and wigge , 2010 ) .
we also observed that swc6 and hta9 hta11 mutants display an enrichment of hydrogen peroxide - related go terms like in arp6 ( supplemental dataset 2 ) .
taken together , the results of the transcriptome analysis demonstrate that the components of the swr1c have distinct roles in the regulation of gene expression in arabidopsis .
our rna - seq analysis revealed that , in the swr1c mutants , particularly pie1 , genes required for disease resistance were misregulated as previously shown ( march - daz et  al . , 2008 ) .
genes involved in systemic acquired resistance such as pr1 , pr5 , eds5 , and nimin1 were highly upregulated in pie1 ( figure  5a5d and supplemental dataset 1 ) .
only pr5 among these genes was upregulated in swc6 and hta9 hta11 , but not to the same extent as in pie1 .
interestingly , the pie1 mutant shows severely compromised resistance in spite of the enhanced expression of defense genes ( as shown above ) . to further understand the molecular basis for this paradox
we observed that genes encoding the wrky family of transcription factors are upregulated in pie1 ( figure  5e and supplemental dataset ) .
wrky transcription factors are important regulatory components during plant response to pathogen infection and abiotic stresses .
moreover , they have been implicated in plant immunity as both positive and negative regulators ( pandey and somssich , 2009 , rushton et  al . ,
wrky38 and wrky62 , in particular , are highly expressed in the pie1 background ( figure  5f and 5 g ) .
these two wrky proteins were shown to be negative regulators of resistance to pathogens . in line with these observations ,
the uniquely upregulated genes in pie1 are also enriched in go terms related to negative regulation of the defense response ( p  = 1.02   10 ) .
sa is a key regulator of signaling networks involved in defense along with other hormones such as ja  ( robert - seilaniantz et  al . , 2011 ) .
sa- and ja - mediated defense  responses are triggered in the plant depending on the nature of the pathogen ( van der does et  al . , 2013 ) .
these two pathways act antagonistically to modulate defense responses . activation of one pathway usually correlates with attenuation of the other ( thaler et  al . , 2012 ) .
consistent with this , we found that the sa - inducible glutaredoxin gene grx480 , which has been proposed to antagonize the ja - responsive transcription of pdf1.2 is nearly 10-fold upregulated in pie1 ( figure  6a ) .
transgenic plants ectopically expressing grx480 displayed reduced ja - induced gene expression ( ndamukong et  al . , 2007 ) .
in fact , uniquely upregulated genes in pie1  are  enriched in go terms related to the ja metabolic process ( p  = 9.12   10 ) and ja - mediated signaling pathway ( p  = 5.17   10 ) .
for instance , genes involved in ja biosynthesis ( figure  6b6d ) , signaling ( figure  6e and 6f ) , and response ( figure  6g6i ) were upregulated .
we further analyzed the expression of ja - induced genes ( nemhauser et  al .
, 2006 ) in our data and found that these are upregulated in pie1 compared with the other mutants ( figure  6j ) , confirming that ja - responsive pathways are derepressed in spite of the sar genes being upregulated .
these observations , together with the upregulation of negative regulators of defense ( e.g. , wrkys ) , could explain the unexpected disease susceptibility phenotypes we observe .
they also point to a potential role that pie1 plays in coordinating distinct defense signaling pathways , especially in the sa - ja antagonism .
our analyses revealed that mutants of the swr1c components show considerable differences in physiological and molecular phenotypes , leading to the assumption that the complex components may be functionally specialized .
it is also likely that the non - enzymatic subunits might not necessarily be essential for complex function , and their depletion or inclusion in the complex could modulate its function .
the components could also have functions independent of h2a.z incorporation . to further elucidate the possible specialization and to understand their role in relation to each other , we analyzed double mutant plants . since arp6 and
swc6 mutants show similar morphological and developmental phenotypes but contrasting immunity traits , we analyzed arp6 swc6 double mutants to study disease resistance without the confounding effects of contrasting developmental phenotypes .
the double mutants did not show any noticeable additional morphological phenotypes than the single mutants suggesting overlapping functions in development . as shown above , arp6 showed resistance to dc3000 while the swc6 mutant showed increased susceptibility compared with wild - type col-0 .
interestingly , the arp6 swc6 double mutants displayed increased susceptibility to pst dc3000 phenocopying swc6  ( figure 7 ) .
these results indicate that swc6 is epistatic to arp6 in disease - resistance traits . despite the overlapping functions in development
, it seems that  loss of swc6 significantly compromises the regulation of immunity genes by swr1c more than the loss of arp6 .
our gene expression analyses have shown that genes involved in systemic acquired resistance such as pr1 and pr5 were derepressed in the swr1c mutants and therefore showed increased basal expression ( figure  5 ) .
interestingly , however , pie1 , swc6 , and hta9 hta11 mutants were severely compromised in resistance to both virulent and avirulent p.  syringae strains in spite of the enhanced basal expression of defense genes . to understand the molecular basis for these paradoxical phenotypes , we analyzed defense gene expression in the mutants in response to pst dc3000 .
gene expression was analyzed in 4-week - old wild - type col-0 and the mutants by quantitative rt  pcr 6 h post inoculation with pst dc3000 or mock treatment ( 10  mm mgcl2 ) ( figure  8) .
as expected , pr1 showed very low basal expression in col-0 under control conditions , whereas bacterial infection resulted in a strong induction of gene expression ( figure  8a ) . the arp6 mutant showed derepression of pr1  and a stronger induction in response to pst dc3000 . on the other hand ,
the pie1 mutant showed a strong constitutive expression of pr1 ; however , inoculation with pst dc3000 inoculation resulted in strong downregulation ( figure  8a ) . in swc6 and hta9
hta11 , where basal expression of pr1 was modestly altered , dc3000 infection resulted in its repression ( figure  8a ) .
similarly , analysis of pr5 expression showed that while it was strongly induced in arp6 in response to pst dc3000 , there was little or no induction in pie1 , swc6 , and hta9 hta11 mutants ( figure  8b ) .
these results show that induction of pr1 and pr5 was compromised upon pathogen inoculation in pie1 , swc6 , and hta9 hta11 mutants despite the higher basal expression .
together , these results show that while it is required to maintain defense genes in repressed state in the absence of pathogen signals , swr1c function , and therefore h2a.z are essential for strong pathogen - induced defense gene activation that is essential for disease resistance .
swr1c is an evolutionarily conserved atp - dependent chromatin remodeling complex involved in the post mitotic incorporation of the histone variant h2a.z into nucleosomes .
putative swr1c subunits have been identified and shown to be involved in plant growth and development as well as environmental adaptation .
although they have been implicated in immunity , their exact role is yet to be understood .
above all , the respective roles that each subunit has in gene regulation and their influence in major physiological outcomes remain elusive .
we undertook a comprehensive analysis of swr1c and h2a.z mutants to systematically analyze their function in plant immunity and gene regulation .
our results indicate that perturbing swr1c function or depleting h2a.z levels results in severely impaired immune responses : pie1 , swc6 , and h2a.z are positive regulators of resistance in arabidopsis against both biotrophic and necrotrophic pathogens .
transcriptome analysis revealed further complexities as the swr1c subunits showed overlapping but , more importantly , distinct functions in the regulation of gene expression . in this study , we found that mutations of the swr1c subunits , pie1 and swc6 , and depletion of h2a.z severely impaired immunity suggesting that they are essential for defense responses in arabidopsis .
pie1 , swc6 , and h2a.z have been implicated in immunity in an earlier study , where mutations in the respective genes were reported to result in constitutive activation of defense responses including spontaneous cell death  and increased resistance to virulent p.  syringae ( march - daz et  al . ,
2008 ) . while the developmental phenotypes we observed were in agreement with the previous reports ( deal et  al . , 2005 ,
2008 ) , defense phenotypes were contrasting in our experiments . while p.  syringae infection  assays were performed under the same growth conditions in our study and in that of march - daz et  al .
we used appropriate genotypes as controls ( eds1 , ndr1 , etc . ) in all our experiments  and have shown the expected phenotypes , suggesting that the phenotypic differences are not due to the experimental conditions .
we have also extended our analyses to test if the swr1c mutants show photoperiod - dependent defense phenotypes and performed disease resistance assays using pst dc3000 under long photoperiod .
strikingly pie1 , swc6 , and hta9 hta11 were all increasingly susceptible independent of the photoperiod , suggesting that the difference in the resistance phenotype is likely not due to growth conditions , particularly day length .
it is still likely that other environmental factors such as nutrient availability associated with plant growth conditions could underlie the discrepancy .
however , our reproducible results with both virulent and avirulent p.  syringae strains , as well as the expected phenotypes for both positive and negative control genotypes , support our conclusions that the swr1c mutants pie1 , swc6 , and hta9 hta11 are compromised in disease resistance . as reported previously ( noh and amasino , 2003 , march - daz et  al . , 2008 )
, we found that the pie1 mutant is highly pleiotropic , displaying stunted growth  and low fertility .
moreover , as described above , it showed increased susceptibility even under long photoperiod conditions , where its growth defects were more severe .
these observations suggest that the pleiotropic phenotypes of pie1 are not necessarily a result of constitutive defense activation .
it is , however , still likely that defective regulation of certain defense - associated pathways underlies these phenotypes despite the lack of apparent disease resistance . similar examples where defense pathways are turned on without apparent enhanced resistance have been described previously .
mutation in cpk28 genes for instance , resulted in growth defects and constitutively enhanced ja accumulation and activation of ja - responsive genes , but not enhanced resistance to necrotrophs ( matschi et  al . , 2015 ) .
similar to the impaired basal immunity in response to the virulent pathogen p.  syringae , pie1 and swc6 mutants showed compromised resistance to pst dc3000 carrying avrrps4 and avrrpt2 effectors , suggesting that pie1 and swc6 are required for rps4- and rps2-mediated immunity .
the hta9 hta11 showed only a mild susceptibility , likely due to the compensatory activity of hta8 .
arp6 , however , was dispensable for eti as the mutant showed wild - type resistance .
these discrepancies in nb - lrr - mediated resistance between the different mutants were not due to an altered expression of either rps4 or rps2 ( supplemental dataset ) and further highlight the functional specialization of the complex subunits .
activation of r - mediated defense results in sa accumulation and signaling ( vlot et  al . ,
sa acts as an important intracellular signal for defense gene activation , and mutants impaired in sa accumulation usually display compromised basal defense and eti .
although pie1 and swc6 mutants were clearly defective in mounting resistance to pst dc3000 avrrps4 , we notice a robust increase in sa accumulation in response to this effector , suggesting that rps4-mediated recognition and downstream signaling are unlikely to be altered . in response to pst dc3000 avrrps4
compromised eti in pie1 and swc6 in spite of sa accumulation suggests that sa - mediated downstream responses might be affected in these mutant backgrounds .
sa - induced gene expression and defense activation are dependent on npr1 ( cao et  al . , 1994 , delaney et  al . , 1995 ) .
similar to the results observed for pie1 and swc6 , the npr1 mutant accumulates a high level of sa but fails to mount sar ( delaney et  al . , 1995 , clarke et  al .
npr1 is also required for the negative feedback regulation of sa production ( delaney et  al . , 1995 ) .
the hyperaccumulation of sa following infection with pst dc3000 avrrps4 in spite of enhanced susceptibility suggests a compromised npr1 activity in pie1 and swc6 .
it is therefore tempting to speculate that swr1c activity and h2a.z might be required for proper functioning of npr1-mediated regulation of immunity .
our results showing that pr gene induction following pst dc3000 infection was perturbed in pie1 , swc6 , and hta9 hta11 mutants also support this hypothesis .
in addition , pie1 and swc6 could also be important for npr1- and sa - independent processes required for resistance .
we found that swr1c and h2a.z are also important for resistance to the necrotrophic pathogen botrytis cineria .
depletion of pie1 , swc6 , and h2a.z , not arp6 , in mutants led to increased susceptibility .
this was rather intriguing given the antagonistic interaction between sa- and ja - mediated defenses ( robert - seilaniantz et  al .
the contrasting phenotypes of the mutants show that pie1 , swc6 , and h2a.z , not arp6 have a positive function in ja / et - mediated defense .
this result was not expected especially in the case of pie1 , which showed an increase in the ja - responsive genes , while sa - responsive genes are also upregulated .
moreover , pie1 mutation also resulted in a higher level of expression of grx480 , which is a negative regulator of ja - mediated immunity .
these results indicate that mutations of the swr1c components alter plant resistance and gene expression in general . taken together , our data point to the importance of swr1c function , and
the pleiotropic effect of compromised swr1c function leads to misregulation of a large subset of defense genes , including both positive and negative regulators , giving rise to complex physiological outcomes as a consequence of epistatic interactions .
further molecular dissection would unravel the precise regulatory role of h2a.z in transcriptional reprogramming in immunity .
the discrepant phenotypes observed when arp6 is depleted from swr1c suggest that the complex could regulate different defense pathways depending on the components recruited to it . our double mutant analysis for resistance to p.  syringae supports an epistatic interaction of arp6 and swc6 , which have contrasting roles in immunity , as we observed in the mutants .
extending the genetic dissection would further elucidate the complex genetic interactions between the complex components in relation to resistance to necrotrophs and the coordination of potentially antagonistic signaling pathways .
our study has revealed that the swr1c subunits are functionally specialized and have non - overlapping functions in immunity . the transcriptomic analysis that we performed further established this . in our analyses
pie1 is the arabidopsis homolog of swr1 in yeast , the enzymatic subunit of the complex ( noh and amasino , 2003 ) .
surprisingly , transcriptional misregulation in pie1 and hta9 hta11 double mutant showed a poor correlation ( figure  4e ) , while hta9 hta11 still has a functional hta8 , likely explaining the poor agreement .
a recent study has also shown similar non - overlapping changes while comparing the transcriptome of pie1 with that of h2a.z triple mutant ( coleman - derr and zilberman , 2012 ) .
it is likely that pie1 and h2a.z have also independent roles of one another , and that pie1 could functions independently of h2a.z incorporation .
consistent with their contrasting effects in plant immunity , arp6 and swc6 mutants clearly maintained distinct molecular phenotypes as supported by the gene expression analysis . in budding yeast , arp6 and swc6 have been shown to be essential subunits in swr1c for h2a.z deposition ( mizuguchi et  al .
, 2004 , wu et  al . , 2005 ) . along with swc2 , another component of the swr1c in yeast
, they act as a sub - complex that requires all three proteins for their association with the complex and for histone exchange ( mizuguchi et  al . , 2004 , wu et  al . , 2005 ) .
the contrasting immunity phenotypes and the divergent transcriptome we observed in our analyses suggest that arp6 and swc6 do not necessarily always act as part of the same complex and that they could have independent functions in arabidopsis , either related or not to h2a.z .
it has been previously proposed that arp6 could have h2a.z - independent roles in gene expression ( yoshida et  al . , 2010 ) .
of all the mutant pairs studied , swc6 and hta9 hta11 showed the strongest correlation in gene expression regulation ( figure  4 g ) suggesting that swc6 functions solely as a subunit of swr1c .
it remains to be seen how much of the gene expression changes that we observe in the mutants are due to changes in h2a.z incorporation .
further investigations are required to fully elucidate the role that each subunit could have in histone exchange and therefore gene regulation .
detailed analysis of the rna - seq data revealed that the pie1 mutant showed significant upregulation of both positive and negative regulators of defense .
pathogen - induced wrky transcription factors function as positive or negative regulators of immunity ( eulgem and somssich , 2007 , pandey and somssich , 2009 ) , and epistatic interaction between the two types of regulators is required for the proper control and fine tuning of defense . of particular interest
was the strong upregulation of wrky38 and wrky62 , which have been shown to be negative regulators of immunity , and their overexpression has led to reduced pr1 gene induction in response to bacterial infection ( kim et  al . , 2008 ) .
this is consistent with the compromised pathogen - induced pr1 expression in pie1 mutant ( figure  8) .
our data have shown that the pie1 mutant showed severely compromised immunity in spite of upregulation of a remarkable number of genes involved in sar such as pr1 . while constitutive activation of sar genes is indicative of pie1 playing a repressive role in sa - mediated gene expression in the absence of pathogens ,
data shown in figure  2c show  that lack of resistance is not due to the inability of pie1 to accumulate sa .
in addition , compromised defense gene induction in response to p.  syringae ( figure  8) indicates that while pie1 has a clear role in maintaining sar genes repressed in the absence of pathogens , it is essential for proper gene induction and resistance upon infection . a complex network of signaling pathways underlies defense responses in plants .
defense signaling elicited by the plant hormones sa and ja are mutually antagonistic and determine defense against biotrophic and necrotrophic pathogens , respectively ( thaler et  al . , 2012 ) .
ja biosynthesis and ja - responsive genes are highly upregulated in pie1 , suggesting that pie1 is also a negative regulator of ja - mediated responses .
remarkably , these results also implicate pie1 activity in regulating sa / ja signaling trade - off .
our data suggest  that swr1c and h2a.z might be playing an important role in maintaining the fine balance between signaling pathways in plant immunity
. therefore , in the context of pie1 , where activating and repressive factors are similarly de - repressed , the physiological outcome will be the consequence of epistatic interactions . in support of this
, we found that pie1 is compromised in resistance to both biotrophic and necrotrophic pathogens . compromised defense gene induction in the susceptible mutant backgrounds such as pie1 , swc6 , and hta9 hta11 in response to pst
further molecular dissection is needed to elucidate the regulatory dynamics of defense gene expression in the context of chromatin remodeling .
we have shown that , although they are as part of the same complex , swr1c subunits have divergent functions in gene regulation and therefore in important biological processes such as immunity . in this comprehensive study , we investigated defense responses and gene expression and have clearly shown the involvement of swr1c and h2a.z in immunity as positive regulators .
although the swr1c component proteins are not involved in defense signaling events directly , they exert their roles through gene regulation .
the genome - wide transcriptomic analysis of the corresponding mutants is consistent with the early reports showing that swr1c and h2a.z are involved in regulating defense - related genes ( march - daz et  al . , 2008 , coleman - derr and zilberman , 2012 ) .
this is particularly in line with the observation that responsive genes are enriched in h2a.z on gene bodies ( coleman - derr and zilberman , 2012 ) .
chromatin state and architecture have been implicated in the repression of sar and stress - related gene expression in general ( bartsch , 2006 , hamon and cossart , 2008 , van den burg and takken , 2009 , ma et  al . , 2011 ,
our extensive analyses highlight that , despite the evolutionary conservation of fundamental molecular processes such as chromatin remodeling , in arabidopsis swr1c components carry out non - overlapping roles
. this can be either due to some of the non - enzymatic subunits having roles outside the swr1c or the potential variation of the degree to which depletion or inclusion of each of the subunits affects the complex function .
it is likely that some of the non - enzymatic components such as swc6 and arp6 are non - essential for complex activity and that their depletion will affect the complex differently .
however , further molecular dissection is required to comprehensively define the role that each component protein plays in the assembly , recruitment to specific target , and histone exchange activity .
importantly , our results that h2a.z and swr1c are positive regulators of plant resistance to different pathogens are particularly significant in the context of environmental modulation of defense responses . elevated temperatures have been shown to suppress plant immunity ( wang et  al . , 2009 ,
nucleosome dynamics have been shown to underlie temperature responses in arabidopsis ( kumar and wigge , 2010 ) . therefore , compromised immunity in the swr1c subunit mutants where h2a.z deposition is perturbed could be reminiscent of the influence of temperature response .
further studies will unravel the role of h2a.z - mediated chromatin dynamics in the environmentally modulated defense gene expression .
we used the previously described mutants arp6 - 1 ( deal et  al . , 2005 ) , pie1 - 2 ( noh and amasino , 2003 ) , swc6 - 1 ( lzaro et  al . , 2008 ) , hta9 - 1 , hta11 - 1 ( march - daz et  al . , 2008 ) , ndr1 - 1 ( century et  al . , 1995 ) , and eds1 - 2 ( bartsch , 2006 ) , and the transgenic complementation line swc6 35s - myc : atswc6 ( choi et  al . , 2007 ) in our experiments . for bacterial infection assays ,
the ros burst assay , and trypan blue leaf staining , arabidopsis plants were grown on soil for 5  weeks in sd conditions ( 8  h light/16  h dark ) in sanyo growth cabinets at 22c and 70% relative humidity .
otherwise 2-week - old seedlings grown in half - strength murashige and skoog ( ms ) plates after seed - surface sterilization were used .
briefly , spray inoculations were performed with bacterial suspension at 2   10 and 2   10 cfu / ml of p.  syringae pv tomato dc3000 and pst dc3000 carrying avrrps4 or avrrpt2 genes , respectively , in 10  mm mgcl2 with 0.04% silwet l-77 .
infiltration experiments were performed with bacterial suspensions at 10 cfu / ml in 10  mm mgcl2 using a needless syringe .
bacterial titers where determined at day 0 ( 2 h ) and/or at day 3 post inoculation .
significant differences to col-0 are expressed using student 's t - test at p  < 0.05 .
the botrytis cinerea infection assay was performed as previously described in schoonbeek et  al .
a botrytis spore suspension ( 2.5   10 spore / ml ) prepared in quarter - strength potato dextrose broth ( pdb , 6
g / l ) was used for drop inoculation on the leaves of 5-week - old arabidopsis plants .
five leaves were inoculated with one droplet of spores of botrytis cinerea isolate  b05.10 .
free sa was extracted twice from 10  mg of freeze - dried leaves with 10% methanol and 1% acetic acid in the presence of 125  nm 2-hydroxybenzoic-3,4,5,6-d4 acid as internal standard ( cdn isotopes ) .
the samples were run on a xevo tqs tandem mass spectrometer attached to an acquity uplc system .
separation was on a 50   2.1 mm 2.7 kinetex xb c18 column ( phenomenex ) running in a gradient of acetonitrile versus 0.1% formic acid in water . for rna - seq analysis ,
total rna was extracted from 2-week - old seedlings grown on  ms solid medium using the rneasy plant mini kit ( qiagen ) according to the manufacturer instructions .
the 50-bp - long single - end reads were mapped to the arabidopsis genome ( tair10 ) and quantified to generate rpkm ( reads per kilobase of exon per million fragments mapped ) values .
the tophat - cufflinks pipeline was used to align rna - seq reads to arabidopsis thaliana reference genome assembly .
differential gene expression and pca analysis was performed with the rpkm values generated using strand ngs 2.0 ( agilent ) .
the expression values for the three biological replicates were analyzed for statistical significance using a one - way analysis of variance test applying the benjamini - hochberg false discovery rate ( fdr ) and the  multiple testing correction ( benjamini and hochberg , 1995 ) and asymptotic p value computation with an adjusted p value cut - off of p  < 0.05 .
genes with absolute fold change values of at least two were included in our analysis as differentially expressed genes . for the comparison of our transcriptome data with the published results , we considered the 15  199 genes that passed the fdr cut - off of < 0.05 as the global number .
four - week - old wild - type and mutant plants were used for analyzing defense gene expression in response to bacterial inoculation .
plants were either spray inoculated with pst dc3000 or were mock treated with 10  mm mgcl2 .
total rna was isolated as above , and 1  g of rna was used as a template for  reverse transcription with superscript iii reverse transcriptase ( invitrogen ) and oligo dt according to the manufacturer 's instructions .
ef1 ( at5g60390 ) transcript levels were used as internal reference ( oligos : 5 tacgccccagttctcgattg 3 and 5 ggcttggttggggtcatctt 3 ) for analyzing pr1 ( oligos : 5 accaggcacgaggagcggta 3 and 5 tccccgtaaggcccaccaga 3 ) and pr5 ( oligos : 5 acccacagcacagagacacaca 3 and 5 tggccataacagcaatgccgc 3 ) expression .
quantitative rt  pcr experiments were performed in a light cycler lc480 using light cycler 480 sybr green i master ( roche ) .
this work was supported by a biotechnology and biological sciences research council ( bbsrc ) grant bb / i019022/1 and institute strategic program grants bb / j004588/1 and bb / j004553/1 .
designed experiments ; s.b . performed most of the experiments and data analyses ; s.n.g .
performed experiments and data analysis for figures  7 and 8 ; s.b . and s.v.k .
 incorporation of the histone variant h2a.z into nucleosomes by the swr1 chromatin remodeling complex is a critical step in eukaryotic gene regulation . in arabidopsis , swr1c and h2a.z have been shown to control gene expression underlying development and environmental responses . although they have been implicated in defense , the specific roles of the complex subunits and h2a.z in immunity are not well understood . in this study 
 , we analyzed the roles of the swr1c subunits , photoperiod - independent early flowering1 ( pie1 ) , actin - related protein6 ( arp6 ) , and swr1 complex 6 ( swc6 ) , as well as h2a.z , in defense and gene regulation . 
 we found that swr1c components play different roles in resistance to different pathogens . 
 loss of pie1 and swc6 function as well as depletion of h2a.z led to reduced basal resistance , while loss of arp6 fucntion resulted in enhanced resistance . 
 we found that mutations in pie1 and swc6 resulted in impaired effector - triggered immunity . 
 mutation in swr1c components and h2a.z also resulted in compromised jasmonic acid / ethylene - mediated immunity . 
 genome - wide expression analyses similarly reveal distinct roles for h2a.z and swr1c components in gene regulation , and suggest a potential role for pie1 in the regulation of the cross talk between defense signaling pathways . 
 our data show that although they are part of the same complex , arabidopsis swr1c components could have non - redundant functions in plant immunity and gene regulation . 
Introduction
Results
Discussion
Methods
Funding
Author Contributions